Search NASA⌕ Search

SEARCH · Search NASA

Results for “WRF”

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.

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

At least 415 records · Page 23

The Contribution of Subtropical Moisture Within an Atmospheric River on Moisture Flux, Cloud Structure, and Precipitation Over the Salmon River Mountains of Idaho Using Moisture Tracers

The impact of an atmospheric river (AR) on the flux of subtropical moisture across Idaho's Salmon River Mountains and precipitation over the mountains is evaluated using the Weather, Research, and Forecasting model with water vapor tracers (WRF-WVT). The AR impacted Idaho between 17 and 19 January 2017 during the Seeded and Natural Orographic Wintertime Clouds: The Idaho Experiment (SNOWIE) campaign. WRF-WVT is configured to isolate the subtropical moisture contribution to the AR, the moisture flux, and precipitation. Subtropical water vapor advected by the AR into Idaho is tagged and tracked in three-dimensional space throughout the run. This allows the contribution of the subtropical moisture to the vertical distribution of water vapor and the precipitation to be directly calculated. The simulated cloud structure is compared with airborne radar data collected during two SNOWIE intensive operation periods. This study found that more than 70% of the moisture flux and more than 80% of the precipitation across the Idaho Mountains during SNOWIE IOP 4 could be attributed to subtropical moisture within the AR. Nearly all of the moisture flux in the upper cloud and 50% of the moisture in the lower cloud was attributable to the subtropical moisture. The subtropical moisture contribution within the AR to precipitation ranged from 35% in northern Idaho to more than 90% in southern Idaho. Across the entire period of impact of the AR, more than 60% of precipitation in Idaho was attributable to the subtropical moisture within the AR, with this percentage increasing toward the south across the state.

54 ENVIRONMENTAL SCIENCES↗

Assessing the Atmospheric Response to Subgrid Surface Heterogeneity in the Single-Column Community Earth System Model, Version 2 (CESM2)

Land-atmosphere interactions are central to the evolution of the atmospheric boundary layer and the subsequent formation of clouds and precipitation. Existing global climate models represent these connections with bulk approximations on coarse spatial scales, but observations suggest that small-scale variations in surface characteristics and co-located turbulent and momentum fluxes can significantly impact the atmosphere. Recent model development efforts have attempted to capture this phenomenon by coupling existing representations of subgrid-scale (SGS) heterogeneity between land and atmosphere models. Such approaches are in their infancy and it is not yet clear if they can produce a realistic atmospheric response to surface heterogeneity. Here, we implement a parameterization to capture the effects of SGS heterogeneity in the Community Earth System Model (CESM2), and compare single-column simulations against high-resolution Weather Research and Forecasting (WRF) large-eddy simulations (LESs), which we use as a proxy for observations. The CESM2 experiments increase the temperature and humidity variances in the lowest atmospheric levels, but the response is weaker than in WRF-LES. In part, this is attributed to an underestimate of surface heterogeneity in the land model due to a lack of SGS meteorology, a separation between deep and shallow convection schemes in the atmosphere, and a lack of explicitly represented mesoscale secondary circulations. These results highlight the complex processes involved in capturing the effects of SGS heterogeneity and suggest the need for parameterizations that communicate their influence not only at the surface but also vertically.

54 ENVIRONMENTAL SCIENCES↗

Impact of Urban Representation on Simulation of Hurricane Rainfall

Abstract Taking the examples of Hurricane Florence (2018) over the Carolinas and Hurricane Harvey (2017) over the Texas Gulf Coast, the study attempts to understand the performance of slab, single‐layer Urban Canopy Model (UCM), and Building Environment Parameterization (BEP) in simulating hurricane rainfall using the Weather Research and Forecasting (WRF) model. The WRF model simulations showed that for an intense, large‐scale event such as a hurricane, the model quantitative precipitation forecast over the urban domain was sensitive to the model urban physics. The spatial and temporal verification using the modified Kling‐Gupta efficiency and Method for Object based Diagnostic and Evaluation in Time Domain suggests that UCM performance is superior to the BEP scheme. Additionally, using the BEP urban physics scheme over UCM for landfalling hurricane rainfall simulations has helped simulate heavy rainfall hotspots.

Geology↗

Enhancing Extreme Precipitation Predictions With Dynamical Downscaling: A Convection‐Permitting Modeling Study in Texas and Oklahoma

Precipitation in the Southern Plains of the United States is relatively well depicted by the Community Earth System Model (CESM). However, despite its ability to capture seasonal mean precipitation anomalies, CESM consistently underestimates extreme pluvial and drought events, rendering it an insufficient tool for extending simulation lead times for exceptional events, such as the abnormally dry May 2011, which helped drive Texas into its worst period of drought in more than a century, and the abnormally wet May 2015, which led to widespread flooding in that state. Ensemble‐based regional climate experiments are completed for the two extreme years using Weather Research and Forecasting model (WRF) and downscaled from CESM. WRF simulations are at convection‐permitting grid resolution for improved physical representation of simulated precipitation over the Southern Great Plains. By integrating convection‐permitting models (CPMs) into each individual member of a CESM ocean data assimilation ensemble, this study demonstrates that high‐resolution dynamical downscaling can improve model skillfulness at capturing these two events and is thus a potentially useful tool for forecasting extremely high and extremely low precipitation events at subseasonal or even seasonal lead times.

54 ENVIRONMENTAL SCIENCES↗

Effects of Pollen on Hydrometeors and Precipitation in a Convective System

Anemophilous (wind-driven) pollen is one type of primary biological aerosol particle, which can rupture under high humidity conditions and form smaller sub-pollen particles (SPPs). Both pollen and SPPs can reach the upper troposphere under convective conditions, acting as cloud condensation nuclei (CCN) and ice nucleating particles (INPs), thus influencing cloud formation and precipitation. However, the impacts of these biological aerosols on cold cloud formation and local climate remain unclear as there are large uncertainties on their emission flux and ice nucleating abilities. Here, we incorporate pollen emission and rupture processes in the Weather Research and Forecasting Model with Chemistry (WRF-Chem) simulations and update the Morrison microphysics scheme within WRF-Chem using aerosol-aware INP parameterizations to account for pollen in addition to other anthropogenic and biogenic aerosol. INP parameterizations for pollen and SPP are derived from laboratory experiments. When including pollen rupture rates as observed in a series of chamber studies, SPP concentrations increase, leading to an increase of cloud ice and water by up to 50% and potentially extending the duration of the convective system. Among all simulated hydrometeors, graupel and raindrops exhibit the largest enhancements from the inclusion of SPPs, with intensifying precipitation at the backside of the convective system and a greater spatial extent. Sensitivity simulations indicate that SPPs have a greater effect on cloud microphysical processes than whole pollen grains, and further observational evidence is needed to constrain these processes.

54 ENVIRONMENTAL SCIENCES↗

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↗

The Impact of Circumsolar Diffuse Radiation on the Numerical Computation and Forecasting of DNI

Weather Research and Forecasting model with solar extensions (WRF-Solar) uses the Fast All-sky Radiation Model for Solar applications (FARMS) to predict solar radiation from the forecasted atmospheric gases, aerosols, and clouds. The Beer-Bouguer-Lambert law utilized by FARMS as well as other radiative transfer models assumes Direct Normal Irradiance (DNI) is a narrow beam only covering the sun disk. This assumption is interpreted differently in solar observations by pyrheliometers where the circumsolar diffuse radiation is accounted. This study precisely computes circumsolar diffuse radiation in all-sky conditions using a radiative transfer model based on a parameterization of precomputed cloud transmittance and reflectance. The model is implemented in WRF-Solar to provide day-ahead DNI forecast in the north America. The computation and forecasting of DNI is compared with the Beer-Bouguer-Lambert law to understand the impact of circumsolar diffuse radiation.

aerosols↗

AutoBEM-DynamicArchetypes

Automatic Building Energy Modeling (AutoBEM, https://bit.ly/AutoBEM) has been used to create an OpenStudio and EnergyPlus building energy model of 122.9 million U.S. buildings (https://bit.ly/ModelAmerica). Simulating and analyzing a model of every building for large areas (e.g. cities) is often not feasible. This dynamic archetyping capability uses a representative building and calculates a floor-space multiplier that allows millions of buildings to be represented by less than 100 buildings. This script (WRF_Archetypes_Parallel.py) calculates these building archetypes for each of the grid cells from a Weather Research and Forecasting (WRF) model in a parallel fashion. The script works by looping through each of the grid cells in the shapefile in parallel, spatially joining the building metadata table to each grid cell, aggregating relevant archetypes and calculating necessary statistics related to area and number of buildings in each cell. The output is a table (.csv) in which each row is an archetype building with properties about that building as well as statistics that relate that building to the total cell (such as an area multiplier). The following inputs are required: WRF zone shapefile (.shp) (wrf-grids-origin_Vegas_Select_100.geojson) The projection of the shapefile ("EPSG:XXXX") Input table containing building metadata for area corresponding to shapefile (.csv) (https://zenodo.org/record/4552901#.YZQEotDMJPY - ClarkCounty2.csv) The number of cores that will be parallelized (integer) The output file name for the archetype table (.csv) Sample command line inputs: python3 ~/WRF_Archetypes_Parallel.py -i ~/wrf-grids-origin_Vegas_Select_100.geojson -c ~/ClarkCounty2.csv -o ~/OutputArchetypes.csv -j 72 -e EPSG:4326 Using Geopandas Version 0.9.0

Bass, Brett (0000000240988434)↗

Impact of Meteorological Factors on the Mesoscale Morphology of Cloud Streets during a Cold-Air Outbreak over the Western North Atlantic

We report postfrontal clouds (PFC) are ubiquitous in the marine boundary layer, and their morphology is essential to estimating the radiation budget in weather and climate models. Here we examine the roles of sea surface temperature (SST) and meteorological factors in controlling the mesoscale morphology and evolution of shallow clouds associated with a cold-air outbreak that occurred on 1 March 2020 during phase I of the Aerosol Cloud Meteorology Interactions over the Western Atlantic Experiment (ACTIVATE). Our results show that the simulated PFC structure and ambient conditions by the Weather Research and Forecasting (WRF) Model are generally consistent with observations from GOES-16 and dropsonde measurements. We also examine the thermodynamical and dynamical influences in the cloud mesoscale morphology using WRF sensitivity experiments driven by two meteorological forcing datasets with different domain-mean SST and spatial gradients, which lead to dissimilar values of hydrometeor water path and cloud core fraction. The SST from ERA5 leads to weaker stability and higher inversion height than the SST from FNL does. In addition, the use of large-scale meteorological forcings from ERA5 yields a distinctive time evolution of wind direction shear in the inner domain, which favors the formation and persistence of longer cloud rolls. Both factors contribute to a change in the time evolution of domain-mean water path and cloud core fraction of cloud streets. Our study takes advantage of the simulation driven by the differences between two large-scale forcing datasets to illustrate the importance of SST and wind direction shear in the cloud street morphology in a realistic scenario.

54 ENVIRONMENTAL SCIENCES↗

Analyzing Relationships between Tropical Cyclone Intensity and Rain Rate over the Ocean Using Numerical Simulations

In this study, the relationship between tropical cyclone (TC) intensity and rain rate over the ocean is investigated using a full-physics numerical model (WRF) and a physics-based TC rainfall model (TCR). TC intensity is found to be nearly linearly correlated with the average rain rate in the inner core [~0.97 (mm h –1 m –2 )/(m s –1 )], while the correlation is weak at outer radii. This difference is induced because TC intensity is significantly correlated with both the vertical velocity and specific humidity in the inner core but is not significantly correlated with the vertical velocity in the outer radii. Further investigation shows that the intensity–rain-rate relationship at the outer radii is influenced by the TC evolution stage. The rain rate for the outer radii is positively correlated with TC intensity for nondecaying TCs, while this correlation is reduced for decaying TCs due to systematic downdrafts in the outer radii. In the context of climate change, the sensitivity of the TC rain rate to sea surface temperature (SST) is found to be +9% per 1 K increase of SST, roughly the product of the sensitivity of TC intensity to SST (+3%) and the Clausius-Clapeyron scaling (+7%). Coupled with synthetic storms, evolution of the TC rain rate over the twenty-first century under the SSP5-8.5 scenario is projected by the TCR (calibrated with the WRF simulations). The annual increase rates of averaged TC rain rate are 0.17% and 0.20% for the inner core and outer radii, respectively, larger than the annual increase rate of TC intensity (0.046%) but comparable to that of cube of intensity (0.18%).

54 ENVIRONMENTAL SCIENCES↗

Moisture Recycling through Pumping by Mesoscale Convective Systems

Abstract Moisture recycling, the contribution of local evapotranspiration (ET) to precipitation, has been studied using bulk models assuming a well-mixed atmosphere. The latter is inconsistent with a climatologically stratified atmosphere that slants across latitudes. Reconciling the two views requires an understanding of overturning associated with different weather systems. In this study, we aim to better understand moisture recycling associated with mesoscale convective systems (MCSs). Using a convection-permitting WRF simulation equipped with water vapor tracers (WRF-WVT), we tag moisture from terrestrial ET in the U.S. Southern Great Plains during May 2015, when more than 20 MCS events occurred and produced significant precipitation and flooding. Water budget analysis reveals that approximately 76% of terrestrial ET is advected away from the region while the remaining 24% of terrestrial ET is “pumped” upward within the region, accounting for 12% of precipitation. Moisture recycling peaks during early night hours (1800–2400 LT) due to the mixing of the daytime accumulated ET by active convection. By focusing on five “diurnally driven” MCSs with less large-scale circulation influence than other MCSs during the same period, we find an upright pumping of terrestrial ET at the MCS initiation and development stages, which diverges into two branches during the MCS mature and decaying stages. One branch in the upper level advects the ET-sourced moisture downstream, while the other branch in the mid-to-upper level contributes to the trailing precipitation upstream. Overall, our analysis depicts a pumping mechanism associated with MCSs that mixes local ET vertically, highlighting its specific contributions to enhancing convective precipitation processes.

54 ENVIRONMENTAL SCIENCES↗

How Frequent Will the Rarest Daily Rainfall Records of Hurricane Ida’s Remnants Be in the Future?

Abstract Gaining continued insights into the impact of global warming on the occurrence of hurricane-associated intense record downpours is essential for building climate resilient communities. This study investigates projected future changes in extreme rainfall over the Northeast United States, as represented by extreme daily amounts during Hurricane Ida in 2021. We used historical control simulations of Weather Research and Forecasting (WRF) Model generated from 40 years of weather events (1980–2014, 12 km) forced by the fifth generation European Centre for Medium-Range Weather Forecasts atmospheric reanalysis. These simulations are thermodynamically modified (2060–2100) via an imposed warming for the high-emission scenario of shared socioeconomic pathway (SSP585) from a range of general circulation models. Ground observations from the Global Historical Climatology Network (1950–2014) and WRF simulations (historical, 1980–2014, and future, 2060–2100) are integrated into a nonstationary generalized extreme value (GEV) framework to assess the frequency of Ida’s heaviest daily rain rates under the SSP585 scenario. Results show that Ida’s daily maximum rainfall recorded at different observation locations was higher than the single highest September daily maximum observed (1950–2014) for 5 out of 17 stations (∼30% of the stations). Ida-like extreme daily rain rates are projected to be, on average, more than 2 times more likely to occur at the end of the century in the simulations (with some regions as high as 5 times). This work demonstrates that integrating a high-resolution atmospheric model’s present-day and thermodynamically modified future simulations along with ground observations, within a nonstationary statistical framework, is crucial for understanding changing characteristics of extreme weather events. Significance Statement Daily scale extreme precipitation is expected to become more frequent and severe, as evidenced by observations and model simulations. While it is important to investigate how these intensifying heavy rainfall events affect current engineering standards, fewer studies have contextualized how warming impacts the most extreme rainfall from a single storm event relative to historical heavy downpours. In this study, we focused on the daily extreme rainfall associated with the extratropical transition of Hurricane Ida (2021), particularly over the northeastern United States—some of which exceeded the commonly used hydrologic design criteria for a 100-yr storm. Using a high-resolution atmospheric model simulation, we investigated how continued warming may influence the frequency of such daily rain rates. Under a high-emission scenario, these events are projected to become up to 5 times more likely at the end of the twenty-first century.

Dollan, Ishrat J↗

On the Parameterization of Convective Downdrafts for Marine Stratocumulus Clouds

The role of nonlocal transport on the development and maintenance of marine stratocumulus (Sc) clouds in coarse-resolution models is investigated, with a special emphasis on the downdraft contribution. A new parameterization of cloud-top-triggered downdrafts is proposed and validated against large-eddy simulation (LES) for two Sc cases. The applied nonlocal mass-flux scheme is part of the stochastic multiplume eddy-diffusivity/mass-flux (EDMF) framework decomposing the turbulent transport into local and nonlocal contributions. The complementary local turbulent transport is represented with the Mellor–Yamada–Nakanishi–Niino (MYNN) scheme. This EDMF version has been implemented in the Weather Research and Forecasting (WRF) single-column model (SCM) and tested for three model versions: without mass flux, with updrafts only, and with both updrafts and downdrafts. In the LES, the downdraft and updraft contributions to the total heat and moisture transport are comparable and significant. The WRF SCM results show a good agreement between the parameterized downdraft turbulent transport and LES. While including updrafts greatly improves the modeling of Sc clouds over the simulation without mass flux, the addition of downdrafts is less significant, although it helps improve the moisture profile in the planetary boundary layer.

Wu, Elynn↗

Characterization of oceanic post-cold frontal clouds and their model representation. Final report

Low-level clouds are ubiquitous over the oceans, and their impact on the Earth’ radiative budget is important. However, general circulation models (GCMs) that are used to represent our climate and its evolution still experience problems to represent their correct amount and radiative impact. In the extratropical latitudes (~30-60 N/S), cloud properties are strongly modulated by the occurrence of anti-cyclones (high pressure systems) and cyclones (low pressure systems). While deep clouds dominate in cyclones near the center and along the cold and warm fronts, low-level clouds populate the colder regions of the cyclones, and more specifically the region in the wake of the cold fronts, the post-cold frontal region (PCF). The GCMs mentioned above were found to underestimate the amount of these PCF clouds, and this caused errors in their representation of the amount of solar radiation at the surface of the southern hemisphere oceans, causing an excess in absorption and errors in long-term predictions. To better understand which processes are ill-represented in the GCMs to cause such issues, we used ground-based observations from the Eastern North Atlantic (ENA) ARM site and a high resolution regional model (called WRF) to examine the properties of PCF clouds. With the WRF model, we simulated the passage of a cold front at the ENA site using theoretically distinct representations of convection and the physics of the boundary layer. While these physical representations had little impact on the timing, structure and circulation of the cold front, we found a high sensitivity to how the large scale information was fed into the model. Errors rapidly developed if the incoming flow information was injected too far from the site, but 1000 km was found to be an optimal distance between the outer boundary and the ENA site for a realistic cold front passage. With this model, we then conducted a series of experiment where convection and boundary layer schemes were changed, and revealed that 1) the clouds were more sensitive to the convection than boundary layer schemes and 2) the interaction between the two has a large impact. Furthermore, changing the representation of the physics in the model had greater impacts on PCF clouds than changing the initial conditions. As large-scale models increase their spatial resolution, convection will eventually be resolved but it remains that their representation of the physics of the boundary layer will matter for the correct representation of shallow (low-level) clouds. Observations were used to find which large-scale information had a clear relationship with PCF cloud and precipitation properties, in view of helping with GCM parameterizations and evaluation. Both lower troposphere stability and surface forcing were found to be important for PCF clouds that experience more vigorous dynamics than their tropical or subtropical counterparts. For cloud base and top heights, a simple measure of the change in potential temperature between the surface and the 800 hPa was found to be a good predictor. This relationship was found to hold as well when using data from campaigns recently conducted in the southern ocean south of Tasmania. A recent GCM was found to succeed in representing similar relations. For cloud depth, water content, and different characteristics of precipitation both the temperature contrast between the surface and near-air and the near-surface wind speed had a strong relationship. To summarize, through observational constraints and regional model experiments, the project has provided guidance to evaluate GCMs and tools to improve their representation of low-level clouds in the extratropics. It has also highlighted differences in behavior between these clouds and those typically found in the subtropics, thus opening new avenues of research for shallow clouds in general.

54 ENVIRONMENTAL SCIENCES↗

Developing Metrics to Evaluate the Skill and Credibility of Downscaling

In this project, our overarching goal was to evaluate the skill and credibility of statistical and dynamical downscaling techniques for reconstructing the historical climate record and projecting future climate change. This report summarizes our approach and accomplishments. In the first year of the project, the priority was to generate new data sets. For Year 2, the overall goal was to focus comparisons during the historical period and for a subset of variables, especially temperature and precipitation. In Year 3 and beyond, the focus was mainly on comparisons and metrics for future climate change simulations. To quantify the uncertainty associated with gridded dataset choice, the UCLA group compared eight gridded temperature datasets, examining climatologies, variability, and trends. The Scripps group completed statistical downscaling of NARR with the LOCA method, and UCLA and Scripps collaborated on a comparison of historical and end-of-century temperature and precipitation patterns over California from WRF, LOCA-WRF, and LOCA-Livneh. In addition to downscaling temperature and precipitation, the Scripps group has spent considerable time working on near-surface wind and humidity downscaling. The Santa Clara group assembled the Variable Infiltration Capacity (VIC) for nine snow-dominated river basins draining the western slopes of the Sierra Nevada range and assessed the role of snow albedo feedback in their VIC simulations. The results will help inform future downscaling studies.

54 ENVIRONMENTAL SCIENCES↗

Developing Metrics to Evaluate the Skill and Credibility of Downscaling (Final Report)

In this project, our overarching goal was to evaluate the skill and credibility of statistical and dynamical downscaling techniques for reconstructing the historical climate record and projecting future climate change. This report summarizes our approach and accomplishments. In the first year of the project, the priority was to generate new data sets. For Year 2, the overall goal was to focus comparisons during the historical period and for a subset of variables, especially temperature and precipitation. In Year 3 and beyond, the focus was mainly on comparisons and metrics for future climate change simulations. To quantify the uncertainty associated with gridded dataset choice, the UCLA group compared eight gridded temperature datasets, examining climatologies, variability, and trends. The Scripps group completed statistical downscaling of NARR with the LOCA method, and UCLA and Scripps collaborated on a comparison of historical and end-of-century temperature and precipitation patterns over California from WRF, LOCA-WRF, and LOCA-Livneh. In addition to downscaling temperature and precipitation, the Scripps group has spent considerable time working on near-surface wind and humidity downscaling. The Santa Clara group assembled the Variable Infiltration Capacity (VIC) for nine snow-dominated river basins draining the western slopes of the Sierra Nevada range and assessed the role of snow albedo feedback in their VIC simulations. The results will help inform future downscaling studies.

54 ENVIRONMENTAL SCIENCES↗

Developing High-Resolution Constrained Variational Analysis of Vertical Velocity and Advective Tendencies within the Range of ARM Scanning Radars at the SGP

Research progress has been made in two areas. One is about the incorporation of the ARM variationally constrained objective analysis method into the WRF GSI data assimilation system. The other is the development of high resolution ARM data and its applications. Specially, we developed a new data assimilation algorithm by adding dynamical constraints to the WRF GSI data assimilation system using hybrid ensemble variational system to derive 3-D fields of atmospheric dynamics and thermodynamics over the ARM SGP sites. We also developed 4x4 km high-resolution constrained variational analysis data over the SGP during the PECAN and made them available to the community. Details are in the attached report.

54 ENVIRONMENTAL SCIENCES↗

Developing and Testing a Novel Stochastic Ice Microphysics Parameterization for Cloud and Climate Models Using ARM Field Campaign Data (Final Progress Report)

The major goals of this project were: 1) to use field campaign measurements from DOE’s Atmospheric Radiation Measurement (ARM) program to characterize variability of important parameters describing properties of ice particles in the atmosphere; 2) based on this observational analysis, to develop a parameterization scheme for weather and climate models that stochastically varies these parameters, and implement the new scheme into a weather model called the Weather Research and Forecasting model (WRF); 3) to use WRF coupled with the new stochastic scheme to simulate ARM field campaign thunderstorm cases and analyze how accounting for this parameter variability affects the model simulations. This work was performed jointly between the National Center for Atmospheric Research, University of Oklahoma, and University of Utah. To accomplish these goals, we extended an approach previously developed to characterize the variability in the size distribution of ice particles to parameters that are explicitly represented in models (i.e., relationships between ice particle mass and size, and between particle fall velocity and size). Our project was, to our knowledge, the first to apply observationally-constrained estimates of this parameter variability describing mass-size and fall velocity-size in a modeling framework. Our results showed efficacy of the approach, evaluated using ARM observations. Similarly, to our knowledge, work in this project was the first to propose and evaluate in detail a stochastic approach for unresolved turbulent mixing in high-resolution model simulations against detailed, benchmark large eddy simulations and ARM observations. Results showed some promising behavior, particularly with increased mixing and dilution of air in thunderstorm cores with surrounding environmental air, bringing the stochastic simulations closer to the benchmark large eddy simulations; however, results were somewhat degraded using stochastic mixing compared to observations from the AMIE/DYNAMO field campaign. This project also further refined and applied a modeling methodology called “piggybacking” that can robustly separate dynamical and thermodynamic impacts of model changes, and comparison studies of different models based on cases developed from ARM observations. Finally, this project directly supported three graduate students who completed their PhDs as well as a postdoctoral research fellow.

54 ENVIRONMENTAL SCIENCES↗