SEARCH · Engineering Papers

Results for “WRF simulation”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 181 records · Page 10

All sources NASA NTRS DOE OSTI

DOE OSTI2020

Efficacy of the Cell Perturbation Method in Large-Eddy Simulations of Boundary Layer Flow over Complex Terrain

A challenge to simulating turbulent flow in multiscale atmospheric applications is the efficient generation of resolved turbulence motions over an area of interest. One approach is to apply small perturbations to flow variables near the inflow planes of turbulence-resolving simulation domains nested within larger mesoscale domains. While this approach has been examined in numerous idealized and simple terrain cases, its efficacy in complex terrain environments has not yet been fully explored. Here, we examine the benefits of the stochastic cell perturbation method (CPM) over real complex terrain using data from the 2017 Perdigão field campaign, conducted in an approximately 2-km wide valley situated between two nearly parallel ridges. Following a typical configuration for multiscale simulation using nested domains within the Weather Research and Forecasting (WRF) model to downscale from the mesoscale to a large-eddy simulation (LES), we apply the CPM on a domain with horizontal grid spacing of 150 m. At this resolution, spurious coherent structures are often observed under unstable atmospheric conditions with moderate mean wind speeds. Results from such an intermediate resolution grid are often nested down for finer, more detailed LES, where these spurious structures adversely affect the development of turbulence on the subsequent finer grid nest. We therefore examine the impacts of the CPM on the representation of turbulence within the nested LES domain under moderate mean flow conditions in three different stability regimes: weakly convective, strongly convective, and weakly stable. In addition, two different resolutions of the underlying terrain are used to explore the role of the complex topography itself in generating turbulent structures. We demonstrate that the CPM improves the representation of turbulence within the LES domain, relative to the use of high-resolution complex terrain alone. During the convective conditions, the CPM improves the rate at which smaller-scales of turbulence form, while also accelerating the attenuation of the spurious numerically generated roll structures near the inflow boundary. During stable conditions, the coarse mesh spacing of the intermediate LES domain used herein was insufficient to maintain resolved turbulence using CPM as the flow develops downstream, highlighting the need for yet higher resolution under even weakly stable conditions, and the importance of accurate representation of flow on intermediate LES grids.

17 WIND ENERGY↗

DOE OSTI2023

Terrain-Influenced Winds and Fire-Fire Interactions in Wildland Fire Simulations [Dissertation]

Ensemble-based approaches to prescribed fire planning cannot be supported by computational fluid dynamics based models like FIRETEC and the Wildland-Urban Interface Fire Dynamics Simulator (WFDS) because they are too computationally expensive and cannot leverage large eddy simulation approaches like CAWFE and WRF-SFIRE because they have too coarse of resolution. QUIC-Fire was developed to fill this gap but it cannot currently address complex terrain, that is typical for instance in the Western United States. This dissertation describes a variety of improvements made to QUIC-Fire and its various incorporated algorithms in an effort to make it a viable tool in simulating wildland and prescribed fires on terrain. The modifications made to QUIC-Fire are described in three chapters. The first chapter describes the extension of the diagnostic wind model QUIC-URB, the wind engine of QUIC-Fire, to a terrain-following coordinate system. The terraininfluenced winds it generates are analyzed and compared. In particular, this chapter presents the mathematical derivation of the wind solver leading to a linear system of equations that are solved through the successive over-relaxation method. The model is validated against a standard test used in previous works (the Askervein Hill) and against a new dataset from measurements in the Socorro Mountains, New Mexico. The terrain-following implementation captures the correct phenomenology for the isolated Askervein Hill, with a wind speed up at the top of the hill. The model agrees well with measurements on the upwind side of the peak, but overestimates speed-up on the downwind side of the hill. This is due to the inability of the model to generate flow separation and wake-eddy dynamics. On a common laptop, the divergence-free wind field is obtained in 6 s, making the solver appealing for coupled fire-atmosphere simulations. The Socorro Mountain is highly complex, with many cliff faces, peaks, and valleys. Although the model captures the magnitude and direction of inlet and outlet areas of the domain, it performs rather poorly in the valley region and in the regions near the steep cliffs. Hence, the model shows good agreement with data in areas of open sloped terrain but lacks in areas where flow separation and thermally driven effects may be present (neither effect is addressed in this work). In the second chapter the implementation of the terrain-following version of QUIC-URB into QUIC-Fire, and the necessary changes needed to include terrain are described. No changes to the underlying fire spread algorithm are made other than what is required to correctly account for the inclusion of terrain. Previously published FIRETEC results that use five different topographies that share the same centerline profile are compared to simulation results from the modified QUIC-Fire that use the same topographies and fuels. QUIC-Fire results show overall similar behaviors in terms of how the topographies affect fire shapes and trends in spread rates. Due to the terrain-following version of QUIC-URB being unable to generate flow separations at the crest of hills, fire spread rates in these regions across all topographies are over-predicted when compared to FIRETEC. Lateral fire growth shows similar trends with FIRETEC between topographies but does not capture the increase in spread due to a diagonal interface between grassland and forested fuel region of the domain. These results suggest that there are three algorithms within QUIC-Fire that could use improvement: how flame tilt angle is accounted for, the incorporation of non-local drag effects, and the inclusion of the wake-eddy parameterizations that are used in QUIC-URB. Lastly, the third chapter describes a modification to the initial guess used for the calculation of the QUIC-URB mass-conserved wind solution during fire simulations. The modification is aimed at improving fire-fire interactions in QUIC-Fire simulations. The modification consists of using the solution from the previous timestep as the starting point for the calculation of the solution for the next timestep. Fire-fire interactions is greatly improved by the change but a new source of error is introduced. Due to how plumes are modelled in QUIC-Fire the new solution contains errors where gaps in the plume structure are present. However, these errors are mostly limited to the upper atmosphere, where they do not affect fire behavior at the surface, and their magnitude isn’t significant enough to discount the amount of new fire phenomenology now captured in QUIC-Fire with the change.

58 GEOSCIENCES↗

DOE OSTI2020

Stochastic Mixing WRF Setup

Documents describing how to set up simulations using a novel stochastic mixing scheme within the Weather Research and Forecasting model.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2023

Biological particles and aerosol-cloud interactions in the Southern Great Plains (Final Report)

The United States Southern Great Plains (SGP) is the genesis for much of the warm season precipitation in the central and eastern United States. While atmospheric thermodynamics and large-scale dynamics play an important role in formation of precipitation, precipitation mechanisms are also sensitive to aerosols. Aerosols can suppress precipitation as cloud condensation nuclei (CCN) or act as ice nucleating particles (INP) in deep convective clouds. This project focused on understanding the role of primary biological aerosol particles (PBAP) in the region and its influence on cloud formation. Specifically, we focused on biological aerosol in the form of pollen, one type of PBAP that is emitted in large yet variable quantities from vegetation in the mid latitudes. Field observations provide evidence of pollen in the planetary boundary layer and pollen components in cloud droplets and fine particulate matter. Further, pollen grains can easily rupture when wet, forming smaller, sub pollen particles with sizes less than one micron. We evaluated the potential for PBAP events at the Department of Energy (DoE) SGP Atmospheric Radiation Measurement (ARM) research facility due to the rich dataset available. We also developed model simulations that accounted for pollen emission and the generation of sub pollen particles, which have been shown to be both cloud condensation nuclei and ice nucleating particles. The proposed work was designed to address the following questions: 1. What are the physical and chemical signatures of biological aerosol such as pollen in SGP ARM measurements? 2. What is the role of pollen-derived particles on deep convection and precipitation in the Central US? Using measurements from the DoE SGP ARM site and recent airborne campaigns, we evaluated the signatures of pollen and pollen-derived aerosols on optical properties, cloud properties and precipitation (Subba et al., 2021). We identified pollen-driven events over the data record and used the Weather Research and Forecasting Model with fully coupled chemistry (WRF-Chem) to conduct chemically realistic simulations of aerosols during and summer mesoscale convective events in the Southern Great Plains (Subba et al., in review). This work improved our understanding of the role of biological aerosol on clouds and precipitation in the Central United States and placed these results in context with anthropogenically-driven processes.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2023

Biological particles and aerosol-cloud interactions in the Southern Great Plains

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2026

Assessing the Impact of a Forest Canopy on Near-Surface Wind Statistics

Representing the forest canopy in atmospheric numerical models should improve simulated winds within and above the canopy up to a few hundred meters above the ground. Here, in this study, we implement a forest canopy parameterization into the Weather Research and Forecasting (WRF) Model in a large-eddy simulation (LES) mode by applying drag forces across multiple layers within the canopy height. We use unique observations from the Lidar Experiments for Assessing Flow over Forests (LEAFF) field campaign at the Wind River Experimental Forest (WREF) in the U.S. Pacific Northwest to evaluate model performance. In a 2-day case study, the canopy parameterization improved wind predictions both within and above the canopy, particularly during the daytime and at finer grid resolution. Without it, winds were frequently overpredicted above the canopy. Similarly, derived quantities such as the wind shear index also yielded estimates closer to observations with the canopy parameterization implemented. These findings suggest that representing the canopy using drag forces alone can improve simulated mean winds up to 200 m above the surface. Furthermore, second-order statistical moments of wind were more sensitive to canopy density than first-order moments, especially during the daytime. This increased sensitivity and the improved daytime performance in wind speed—evidenced by the lowest bias from observations (3% compared to 20% over diurnal cycle)—imply that winds above the canopy layer are strongly influenced by how well turbulence above the canopy is modeled. The results of this study can serve as a foundation for parameterizing forest canopy effects in coarser weather forecast models.

Energy - Wind↗

DOE OSTI2023

The June 2012 North American Derecho: A Testbed for Evaluating Regional and Global Climate Modeling Systems at Cloud‐Resolving Scales

Abstract In this paper, we introduce a testbed for evaluating and comparing climate modeling systems at cloud resolving scales using hindcasts of the June 2012 North American derecho. To demonstrate its utility for model intercomparison, the testbed is applied to two models: the regionally‐refined Simple Cloud‐Resolving E3SM Atmosphere Model (SCREAM) at 6.5, 3.25 and 1.625 km grid spacing and the Weather Research and Forecasting (WRF) model with 3.2 and 1.6 km grid spacing. We find the simulation results to be highly sensitive to the initial conditions (ICs), initialization time, and model configurations, with ICs from the Rapid Refresh producing the best simulation. Significant improvement is identified in both models as horizontal grid spacing is refined. While a propagation delay of approximately 2 hr is found in both models, SCREAM at 1.625 km simulates the observed bow echo structure of the derecho well and predicts strong surface gusts that exceed 30 m/s. In comparison, WRF has difficulty producing surface wind over 25 m/s, with wind gusts in WRF 42%–46% lower than in SCREAM. However, WRF has a lower bias in simulating cloud top temperature and extent, but overestimates precipitation intensity. Both models reproduce the observed outgoing longwave radiation spatial patterns well (Pearson correlation >0.88), but, compared with NEXRAD observations, simulate generally larger areas of composite radar reflectivity >40 dBZ and underestimate the precipitating area by ∼47%.

2012 North American derecho↗

DOE OSTI2025

Simulating Hurricane Katrina in the Simple Cloud‐Resolving E3SM Atmosphere Model v1

Climate models are important tools for advancing understanding and prediction of tropical cyclones (TCs). Traditional global climate models, however, do not have the ability to properly simulate TC intensity due to their coarse horizontal resolution. Regional models can be run at convection‐permitting resolutions, but these models are often strongly influenced by the data used in the lateral boundary forcing, and domain choice can have a large impact on the simulation. Cloud‐resolving global climate models have demonstrated great potential for realism in TC simulations, and in this study we focus specifically on the Simple Cloud‐Resolving Energy Exascale Earth System Model (E3SM) Atmosphere Model (SCREAM) v1 configuration. We evaluate SCREAMv1 against the observational record and the Weather Research and Forecasting (WRF) model run at a convection‐permitting resolution with Hurricane Katrina as our case study. We found that both models produced realistic simulations of Hurricane Katrina. SCREAMv1 demonstrated skill in simulating TC track, size, and intensity, while the model produced an excessive amount of precipitation. In comparison, WRF more accurately simulated TC precipitation and intensity, although the TC wind extent was smaller than the observations.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2021

Aerosols from anthropogenic and biogenic sources and their interactions – modeling aerosol formation, optical properties, and impacts over the central Amazon basin

The Green Ocean Amazon experiment – GoAmazon 2014–2015 – explored the interactions between natural biogenic forest emissions from central Amazonia and urban air pollution from Manaus. Previous GoAmazon 2014–2015 studies showed that nitrogen oxide (NO x = NO + NO 2 ) and sulfur oxide (SO x ) emissions from Manaus strongly interact with biogenic volatile organic compounds (BVOCs), affecting secondary organic aerosol (SOA) formation. In previous studies, ground-based and aircraft measurements provided evidence of SOA formation and strong changes in aerosol composition and properties. Aerosol optical properties also evolve, and their impacts on the Amazonian ecosystem can be significant. As particles age, some processes, such as SOA production, black carbon (BC) deposition, particle growth and the BC lensing effect change the aerosol optical properties, affecting the solar radiation flux at the surface. This study analyzes data and models SOA formation using the Weather Research and Forecasting with Chemistry (WRF-Chem) model to assess the spatial variability in aerosol optical properties as the Manaus plumes interact with the natural atmosphere. The following aerosol optical properties are investigated: single scattering albedo (SSA), asymmetry parameter (g aer ), absorption Ångström exponent (AAE) and scattering Ångström exponent (SAE). These simulations were validated using ground-based measurements at three experimental sites, namely the Amazon Tall Tower Observatory – ATTO (T0a), downtown Manaus (T1), Tiwa Hotel (T2) and Manacapuru (T3), as well as the U.S. Department of Energy (DOE) Gulfstream 1 (G-1) aircraft flights. WRF-Chem simulations were performed over 7 d during March 2014. Results show a mean biogenic SOA (BSOA) mass enrichment of 512 % at the T1 site, 450 % in regions downwind of Manaus, such as the T3 site, and 850 % in areas north of the T3 site in simulations with anthropogenic emissions. The SOA formation is rather fast, with about 80 % of the SOA mass produced in 3–4 h. Comparing the plume from simulations with and without anthropogenic emissions, SSA shows a downwind reduction of approximately 10 %, 11 % and 6 % at the T1, T2 and T3 sites, respectively. Other regions, such as those further downwind of the T3 site, are also affected. The g aer values increased from 0.62 to 0.74 at the T1 site and from 0.67 to 0.72 at the T3 site when anthropogenic emissions are active. During the Manaus plume-aging process, a plume tracking analysis shows an increase in SSA from 0.91 close to Manaus to 0.98 160 km downwind of Manaus as a result of SOA production and BC deposition.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2025

Historical and Future Windstorms in the Northeastern United States

Large-scale windstorms represent an important atmospheric hazard in the Northeastern US (NE) and are associated with substantial socioeconomic losses. Regional simulations performed with the Weather Research and Forecasting (WRF) model using lateral boundary conditions from three Earth System Models (ESMs: Geophysical Fluid Dynamics Laboratory (GFDL), Hadley Centre Global Environment Model (HadGEM) and Max Planck Institute (MPI)) are used to quantify possible future changes in windstorm characteristics and/or changes in the parent cyclone types responsible for windstorms. WRF nested within MPI ESM best represents important aspects of historical windstorms and the cyclone types responsible for generating windstorms compared with a reference simulation performed with the ERA-Interim reanalysis for the historical climate. The spatial scale and frequency of the largest windstorms in each simulation defined using the greatest extent of exceedance of local 99.9th percentile wind speeds (U > U999) plus 50-year return period wind speeds (U50,RP) do not exhibit secular trends. Projections of extreme wind speeds and windstorm intensity/frequency/geolocation and dominant parent cyclone type associated with windstorms vary markedly across the simulations. Only the MPI nested simulations indicate statistically significant differences in windstorm spatial scale, frequency and intensity over the NE in the future and historical periods. This model chain, which also exhibits the highest fidelity in the historical climate, yields evidence of future increases in 99.9th percentile 10 m height wind speeds, the frequency of simultaneous U > U999 over a substantial fraction (5–25%) of the NE and the frequency of maximum wind speeds above 22.5 ms−1. These geophysical changes, coupled with a projected doubling of population, leads to a projected tripling of a socioeconomic loss index, and hence risk to human systems, from future windstorms.

Pryor, Sara C. (ORCID:0000000348473440)↗

DOE OSTI2025

Evaluating Sea Breezes and Associated Convective Cloud Evolution in the Model Gray Zone

We characterize convective clouds associated with sea‐breeze circulations (SBC) using multi‐agency observations and multi‐case ensemble model simulations. The focus is on assessing convective cloud lifecycle properties and their merging behavior, as well as the environmental conditions they are embedded in, particularly SBC features. In total, 46 SBC days over the Houston‐Galveston region are selected and simulated using the Weather Research and Forecasting (WRF) model at a gray zone scale with a forecast‐like parameterization setup. Advanced techniques, including change‐point detection, a Lagrangian cloud tracking method, and a newly developed cell merging and splitting detection algorithm, are applied and/or developed for this study. Our findings indicate that the WRF model at 1 km grid spacing well represents the thermodynamic conditions over the region, as well as SBC timing and intensity. However, for the associated convective cells, WRF overestimates the 30‐dBZ echo top height, cell area, and maximum radar reflectivity compared to radar observations. This overestimation is potentially due to under‐resolved entrainment processes, an overestimated merging frequency, and the overestimation of updraft intensity. Furthermore, the model exhibits a deficiency in simulating congestus clouds, showing a more rapid transition from shallow to deep convection compared to observed behavior. Moreover, observations indicate stronger, deeper, and wider clouds when merging happens. Conversely, in simulations, the merging process does not necessarily lead to higher or longer‐lived cells, as many cases experience rapid and frequent merging and splitting which may result in more variance in convective updraft velocity during the convection lifetime.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2025

A Scientist-in-the-Loop Data Analytics Framework for Intelligent Simulation Model Tuning and Validation

This project developed a scientist-in-the-loop data analytics framework for intelligent simulation model tuning and validation, targeting the Weather Research and Forecasting (WRF) model and its solar energy variant, WRF-Solar-BNL. Domain experts, such as climate scientists, depend on large-scale numerical simulations for knowledge discovery and decision-making, yet the complexity of parameter tuning and the disconnect between automated optimization and domain expertise pose significant challenges. We extended an interactive visual analytics framework that enables domain experts to observe and intervene in the computational steering process by identifying disagreements between the simulation model, surrogate model, and the expert’s domain knowledge. Using Bayesian Optimization with Gaussian Process Regression as the surrogate model, our system allows users to probe parameter relationships, analyze correlation patterns, and adjust tuning parameters in real time. We developed use cases for solar irradiance forecasting through sustained collaboration with Brookhaven National Laboratory, resolving critical model configuration challenges and achieving meaningful reductions in prediction error. The project supported one PhD student, one MS student, and eight undergraduate students across three Data Science Capstone projects, resulting in one master’s thesis.

Dasgupta, Aritra [New Jersey Institute of Technolo↗

DOE OSTI2020

Summertime Post-Cold-Frontal Marine Stratocumulus Transition Processes over the Eastern North Atlantic

Abstract Marine boundary layer (MBL) cloud morphology associated with two summertime cold fronts over the eastern North Atlantic (ENA) is investigated using high-resolution simulations from the Weather Research and Forecasting (WRF) Model and observations from the Atmospheric Radiation Measurement (ARM) ENA Climate Research Facility. Lagrangian trajectories are used to study the evolution of post-cold-frontal MBL clouds from solid stratocumulus to broken cumulus. Clouds within specified domains in the vicinity of transitions are classified according to their degree of decoupling, and cloud-base and cloud-top breakup processes are evaluated. The Lagrangian derivative of the surface latent heat flux is found to be strongly correlated with that of the cloud fraction at cloud base in the simulations. Cloud-top entrainment instability (CTEI) is shown to operate only in the decoupled MBL. A new indicator of inversion strength at cloud top that employs the vertical gradients of equivalent potential temperature and saturation equivalent potential temperature, which can be computed directly from soundings, is proposed as an alternative to CTEI. Overall, results suggest that the deepening–warming hypothesis suggested by Bretherton and Wyant explains many of the characteristics of the summertime postfrontal MBL evolution of cloud structure over the ENA, thereby widening the phase space over which the hypothesis may be applied. A subset of the deepening–warming hypothesis involving warming initially dominating over moistening is proposed. It is postulated that changes in climate change–induced modifications in cold-frontal structure over the ENA may be accompanied by coincident changes in the location and timing of MBL cloud transitions in the post-cold-frontal environment.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2022

Microphysical Piggybacking in the Weather Research and Forecasting Model

Abstract This paper presents incorporation of the microphysical piggybacking into the Weather Research and Forecasting (WRF) model. Microphysical piggybacking is to run a single simulation applying two microphysical schemes, the first scheme driving the simulation and the second piggybacking this simulated flow. “Driving the simulation” means that the simulated microphysical processes, affect the cloud buoyancy and thus force the simulated flow. In contrast, the piggybacking variables are advected by the simulated flow and undergo microphysical transformation, but they do not affect the simulated flow (like in prescribed flow—kinematic—simulations). The two sets of variables (driver and piggybacker) include temperature, water vapor mixing ratio, and all microphysical variables. We provide details of implementing piggybacking into the WRF model, illustrate its applications, and demonstrate the benefits of this methodology in two idealized three‐dimensional cases: (a) a squall line case applying two microphysics schemes, the Thompson bulk microphysics scheme and the University of Pécs/NCAR bin (UPNB) scheme. The piggybacking simulations revealed that the microphysics‐dynamics interaction plays a more important role than the pure microphysical size sorting effect in the transition zone formation. (b) A case of daytime shallow‐to‐deep convective development over land. This case uses the UPNB scheme and contrasts convection developing in environments with either pristine or polluted cloud condensation nuclei (CCN). The piggybacking results indicated that the increase of cloud cover and decrease of supersaturation are mainly associated with the microphysical effect of increasing CCN while the change of precipitation on the ground is also influenced by microphysics‐dynamics interactions.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2022

Airborne Lidar Measurements of XCO 2 in Synoptically Active Environment and Associated Comparisons With Numerical Simulations

Frontal boundaries have been shown to cause large changes in CO 2 mole-fractions, but clouds and the complex vertical structure of fronts make these gradients difficult to observe. It remains unclear how the column average CO 2 dry air mole-fraction (XCO 2 ) changes spatially across fronts, and how well airborne lidar observations, data assimilation systems, and numerical models without assimilation capture XCO 2 frontal contrasts (ΔXCO 2 , i.e., warm minus cold sector average of XCO 2 ). We demonstrated the potential of airborne Multifunctional Fiber Laser Lidar (MFLL) measurements in heterogeneous weather conditions (i.e., frontal environment) to investigate the ΔXCO 2 during four seasonal field campaigns of the Atmospheric Carbon and Transport-America (ACT-America) mission. Most frontal cases in summer (winter) reveal higher (lower) XCO 2 in the warm (cold) sector than in the cold (warm) sector. During the transitional seasons (spring and fall), no clear signal in ΔXCO 2 was observed. Intercomparison among the MFLL, assimilated fields from NASA's Global Modeling and Assimilation Office (GMAO), and simulations from the Weather Research and Forecasting-—Chemistry (WRF-Chem) showed that (a) all products had a similar sign of ΔXCO 2 though with different levels of agreement in ΔXCO 2 magnitudes among seasons; (b) ΔXCO 2 in summer decreases with altitude; and (c) significant challenges remain in observing and simulating XCO 2 frontal contrasts. A linear regression analyses between ΔXCO 2 for MFLL versus GMAO, and MFLL versus WRF-Chem for summer-2016 cases yielded a correlation coefficient of 0.95 and 0.88, respectively. The reported ΔXCO 2 variability among four seasons provide guidance to the spatial structures of XCO 2 transport errors in models and satellite measurements of XCO 2 in synoptically-active weather systems.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2024

Applying double cropping and interactive irrigation in the North China Plain using WRF4.5

Abstract. Irrigated cultivation exerts a significant influence on the local climate and the hydrological cycle. The North China Plain (NCP) is known for its intricate agricultural system, marked by expansive cropland, high productivity, compact rotation, a semi-arid climate, and intensive irrigation practices. As a result, there has been considerable attention on the potential impact of this intensive irrigated agriculture on the local climate. However, studying the irrigation impact in this region has been challenging due to the lack of an accurate simulation of crop phenology and irrigation practices within the climate model. By incorporating double cropping with interactive irrigation, our study extends the capabilities of the Weather Research Forecast (WRF) model, which has previously demonstrated commendable performance in simulating single-cropping scenarios. This allows for two-way feedback between irrigated crops and climate, further enabling the inclusion of irrigation feedback from both ground and vegetation perspectives. The improved crop modeling system shows significant enhancement in capturing vegetation and irrigation patterns, which is evidenced by its ability to identify crop stages, estimate field biomass, predict crop yield, and project monthly leaf area index. The improved simulation of large-scale irrigated crops in the NCP can further enhance our understanding of the intricate relationship between agricultural development and climate change.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2025

WRF-ELM v1.0: a regional climate model to study land–atmosphere interactions over heterogeneous land use regions

Abstract. The Energy Exascale Earth System Model (E3SM) Land Model (ELM) is a state-of-the-art land surface model that simulates the intricate interactions between the terrestrial land surface and other components of the Earth system. Originating from the Community Land Model (CLM) version 4.5, ELM has been under active development, with added new features and functionality, including plant hydraulics, radiation–topography interaction, subsurface multiphase flow, and more explicit land use and management practices. This study integrates ELM v2.1 with the Weather Research and Forecasting (WRF; WRF-ELM) model through a modified Lightweight Infrastructure for Land Atmosphere Coupling (LILAC) framework, enabling affordable high-resolution regional modeling by leveraging ELM's innovative features alongside WRF's diverse atmospheric parameterization options. This framework includes a top-level driver for variable communication between WRF and ELM and Earth System Modeling Framework (ESMF) caps for the WRF atmospheric component and ELM workflow control, encompassing initialization, execution, and finalization. Importantly, this LILAC–ESMF framework demonstrates a more modular approach compared to previous coupling efforts between WRF and land surface models. It maintains the integrity of ELM's source code structure and facilitates the transfer of future developments in ELM to WRF-ELM. To test the ability of the coupled model to capture land–atmosphere interactions over regions with a variety of land uses and land covers, we conducted high-resolution (4 km) WRF-ELM ensemble simulations over the Great Lakes region (GLR) in the summer of 2018 and systematically compared the results against observations, reanalysis data, and WRF-CTSM (WRF coupled with the Community Terrestrial Systems Model). In general, the coupled WRF-ELM model has reasonably captured the spatial distribution of surface state variables and fluxes across the GLR, particularly over the natural vegetation areas. The evaluation results provide a baseline reference for further improvements in ELM in the regional application of high-resolution weather and climate predictions. Our work serves as an example to the model development community for expanding an advanced land surface model's capability to represent fully-coupled land–atmosphere interactions at fine spatial scales. The development and release of WRF-ELM marks a significant advancement for the ELM user community, providing opportunities for fine-scale regional representation, parameter calibration in coupled mode, and examination of new schemes with atmospheric feedback.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2020

Effects of Cloud Liquid-Phase Microphysical Processes in Mixed-Phase Cumuli Over the Tibetan Plateau

Numerical simulations often overpredict precipitation over the Tibetan Plateau (TP). To examine the factors causing precipitation overprediction, different parameterizations of liquid-phase microphysical processes (accretion, autoconversion, and entrainment-mixing) are implemented into the Morrison microphysics scheme to simulate a TP precipitation event in summer with the Weather Research and Forecasting (WRF) model. The general spatial distribution and temporal trend of precipitation are captured by all simulations, but the precipitation rate is overpredicted. The results from sensitivity experiments suggest that compared to other examined liquid-phase processes, the accretion process is more important in precipitation simulation over the TP region. Further investigation with the Heidke skill scores reveals that accretion parameterization that takes into account the raindrop size produces the most accurate results in terms of the total surface precipitation. This parameterization suppresses spurious accretion and does not produce liquid-phase precipitation until cloud droplets are big enough. Here, it is also confirmed that increasing the model resolution can reduce precipitation overprediction. Results from the case study are confirmed by the use of a one-month simulation.

54 ENVIRONMENTAL SCIENCES↗