Engineering Papers⌕ Search

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.

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

At least 397 records · Page 22

Reconciliation of asynchronous satellite-based $\mathrm{NO_2}$ and $\mathrm{XCO_2}$ enhancements with mesoscale modeling over two urban landscapes

Fossil fuel carbon dioxide (CO 2ff ), the main driver of global warming and climate change, is often co-emitted with nitrogen oxides (NO x ) and precursors to ground-level ozone from anthropogenic sources like power plants or vehicles. In urban and suburban areas, satellite-based NO 2 can be used as a proxy to track the emissions of CO 2ff . Because of NO 2 ’s shorter lifetime, urban NO 2 plumes are more distinguishable from backgrounds and more sensitive to variations in emissions. However, the combination of these two gases is limited by the asynchrony among NO 2 and CO 2 monitoring satellites. We used CO 2ff simulated by the Weather Research and Forecasting model coupled with Chemistry (WRF-Chem) model to reconcile the tropospheric NO 2 vertical column density (VCD) from the Tropospheric Monitoring Instrument (TROPOMI) and column-averaged dry-air mole fractions of carbon dioxide enhancements (ΔXCO 2 ) from Orbiting Carbon Observatory 3 (OCO-3) Snapshot Area Maps (SAMs) over a multicity area, Washington D.C.-Baltimore (DC-Balt), and a basin city, Mexico City. NO 2 /CO 2ff ratios over DC-Balt are smaller than Mexico City, indicative of stricter emission restrictions, a more combustion-efficient vehicle fleet, and higher combustion efficiency due to lower altitude in DC-Balt. For single-track cases, the spatial correlations between NO 2 and ΔXCO 2 over Mexico City are stronger than DC-Balt because the NO 2 and CO 2 are mostly trapped in the valley of Mexico City, while DC-Balt is severely affected by distant sources (i.e., US East Coast cities). Using multi-track averaging, spatial correlation coefficients increase with the number of days used for averaging. The correlations reached a maximum when averaging >12 continuous images for DC-Balt and > 10 continuous images for Mexico City. This finding indicates that multi-track averaging using modeled CO 2ff as a proxy is helpful to filter the noise in single-track images, to cancel the interference from distant sources, and to magnify correlations between NO 2 and CO 2ff . Mexico City showed stronger spatial correlations but weaker temporal correlations than DC-Balt due to biomass burning hot spots and large transport errors caused by the trapping effects of the surrounding mountains. Tracking the 20-day moving average of CO 2ff emissions using TROPOMI NO 2 seems technically feasible, considering the relationship between correlation coefficients and the number of available satellite images.

54 ENVIRONMENTAL SCIENCES↗

Tendencies of Soil Microbial NO Emissions During HI‐SCALE as Predicted by a Nitrification/Denitrification Scheme

Many atmospheric chemical processes, including the formation of secondary organic aerosol (SOA), are strongly modulated by the reactions of NO and NO 2 (NO x ). Though NO x is controlled by anthropogenic emissions near urban areas, in rural areas soil microbes can be a significant contribution to NO emission globally. The relative rates of emissions of different nitrogen-containing species (e.g., NO, N 2 O, HONO, and N 2 ) are strong functions of soil properties such as temperature, moisture content, pH, and soil carbon and nitrogen pools. However, typical large-scale biogeochemical models either express these emissions simplistically, or not at all. Here we investigate the potential impact of soil NO emissions on atmospheric chemistry and SOA formation over regional and monthly time scales, specifically the 2016 spring and summer Intensive Observational Periods (IOPs) of the Holistic Interactions of Shallow Clouds, Aerosols and Land Ecosystems (HI-SCALE) field campaign. We implement the soil NO nitrification/denitrification parameterization of Rasool et al. (2019) into the Weather Research and Forecasting model coupled with Chemistry (WRF-Chem), supplemented by a 1-km soil moisture analysis. We then simulate both IOPs over the U.S. Great Plains, evaluating against ground stations and flight data. We show that soil NO emissions can account for a large fraction of total NO x and locally increase O 3 concentrations by up to 25%, while alleviating negative biases of gases and aerosols toward observations. Soil moisture and temperature changes between IOP1 and IOP2 lead to overall differences in emissions, but with large regional variability due to heterogeneous surface characteristics.

Atmoshpheric Chemistry↗

Studying Cloud and Radiative Impacts Through an Improved Physically Based Representation of Organic Aerosol in a Large-Scale Model (WRF-Chem) (Final Report)

The work performed under this award resulted in the development of a computationally-efficient kinetic model to simulate the chemistry, thermodynamics, and microphysics of organic aerosol. Versions of the kinetic model were used to study a suite of aerosol processes including, multigenerational gas-phase chemistry, phase-state influenced gas/particle partitioning, heterogeneous oxidation, oligomerization reactions, and experimental artifacts. The kinetic model has now been implemented into a regional chemistry-climate model (WRF-Chem) to further study the formation, evolution, properties, and spatiotemporal distribution of organic aerosol at urban, regional, and global scales.

54 ENVIRONMENTAL SCIENCES↗

The Role of Low-Level, Terrain-Induced Jets in Rainfall Variability in Tigris Euphrates Headwaters

Rainfall variability in the Tigris Euphrates headwaters is a result of interaction between topography and meteorological features at a range of spatial scales. Here, the Weather Research and Forecasting (WRF) Model, driven by the NCEP-DOE AMIP-II reanalysis (R-2), has been implemented to better understand these interactions. Simulations were performed over a domain covering most of the Middle East. The extended simulation period (1983 - 2013) enables us to study seasonality, interannual variability, spatial variability, and extreme events of rainfall. Results showed that the annual cycle of precipitation produced by WRF agrees much more closely with observations than does R-2. This was particularly evident during the transition months of April and October, which were further examined to study the underlying physical mechanisms. In both months, WRF improves representation of interannual variability relative to R-2, with a substantially larger benefit in April. This improvement results primarily from WRFs ability to resolve two low-level, terrain-induced flows in the region that are either absent or weak in R-2: one parallel to the western edge of the Zagros Mountains, and one along the east Turkish highlands. The first shows a complete reversal in its direction during wet and dry days, when flowing southeasterly it transports moisture from the Persian Gulf to the region, and when flowing northwesterly it blocks moisture and transports it away from the region. The second is more directly related to synoptic-scale systems and carries moist, warm air from the Mediterranean and Red Seas toward the region. The combined contribution of these flows explains about 50 of interannual variability in both WRF and observations for April and October precipitation.

rainfall variability↗

Cold-Season Precipitation Sensitivity to Microphysical Parameterizations: Hydrologic Evaluations Leveraging Snow Lidar Datasets

Abstract Cloud microphysical processes are an important facet of atmospheric modeling, as they can control the initiation and rates of snowfall. Thus, parameterizations of these processes have important implications for modeling seasonal snow accumulation. We conduct experiments with the Weather Research and Forecasting (WRF V4.3.3) Model using three different microphysics parameterizations, including a sophisticated new scheme (ISHMAEL). Simulations are conducted for two cold seasons (2018 and 2019) centered on the Colorado Rockies’ ∼750-km 2 East River watershed. Precipitation efficiencies are quantified using a drying-ratio mass budget approach and point evaluations are performed against three NRCS SNOTEL stations. Precipitation and meteorological outputs from each are used to force a land surface model (Noah-MP) so that peak snow accumulation can be compared against airborne snow lidar products. We find that microphysical parameterization choice alone has a modest impact on total precipitation on the order of ±3% watershed-wide, and as high as 15% for certain regions, similar to other studies comparing the same parameterizations. Precipitation biases evaluated against SNOTEL are 15% ± 13%. WRF Noah-MP configurations produced snow water equivalents with good correlations with airborne lidar products at a 1-km spatial resolution: Pearson’s r values of 0.9, RMSEs between 8 and 17 cm, and percent biases of 3%–15%. Noah-MP with precipitation from the PRISM geostatistical precipitation product leads to a peak SWE underestimation of 32% in both years examined, and a weaker spatial correlation than the WRF configurations. We fall short of identifying a clearly superior microphysical parameterization but conclude that snow lidar is a valuable nontraditional indicator of model performance.

54 ENVIRONMENTAL SCIENCES↗

Orographic enhancement of rainfall over the Congo Basin

Abstract The Congo rainforest located in central equatorial Africa is an important, yet understudied part of the globe surrounded by complex orographic features. A primitive understanding of precipitation processes such as mesoscale convective dynamics magnifies uncertainties in the future climate projections of the hydrological cycle over the Congo. Furthermore, the effects of orography, which is an important forcing for convection and precipitation, are poorly resolved by climate models, and ill‐conceptualized over the Congo. To address this knowledge gap, perturbed orographic forcing experiments are conducted using the Weather Research and Forecasting (WRF) mesoscale numerical model in a high‐resolution convection‐permitting model setup. The model simulated selected dates in November 2014. The thunderstorms and rainfall simulated in the control run for the case study analyzed in this article compared reasonably well to satellite‐derived brightness temperature and rainfall data. The results from this case study show that the dynamical impact of increasing the height of the East African Highlands is the blocking of the lower‐tropospheric tropical easterlies. This weakening of the lower‐tropospheric zonal winds increases the windshear over the Congo Basin resulting in slower propagating, more intense mesoscale convective systems with enhanced rainfall.

54 ENVIRONMENTAL SCIENCES↗

Structural Uncertainty in the Sensitivity of Urban Temperatures to Anthropogenic Heat Flux

One key source of uncertainty for weather and climate models is structural uncertainty arising from the fact that these models must simplify or approximate complex physical, chemical, and biological processes that occur in the real world. However, structural uncertainty is rarely examined in the context of simulated effects of anthropogenic heat flux in cities. Using the Weather Research and Forecasting (WRF) model coupled with a single-layer urban canopy model, it is found that the sensitivity of urban canopy air temperature to anthropogenic heat flux can differ by an order of magnitude depending on how anthropogenic heat flux is released to the urban environment. Moreover, varying model structures through changing the treatment of roof-air interaction and the parameterization of convective heat transfer between the canopy air and the atmosphere can affect the sensitivity of urban canopy air temperature by a factor of 4. Urban surface temperature and 2-m air temperature are less sensitive to the methods of anthropogenic heat flux release and the examined model structural variants than urban canopy air temperature, but their sensitivities to anthropogenic heat flux can still vary by as much as a factor of 4 for surface temperature and 2 for 2-m air temperature. Our study recommends using temperature sensitivity instead of temperature response to understand how various physical processes (and their representations in numerical models) modulate the simulated effects of anthropogenic heat flux.

54 ENVIRONMENTAL SCIENCES↗

Cloud Computing Applications in Support of Earth Science Activities at Marshall Space Flight Center

Currently, the NASA Nebula Cloud Computing Platform is available to Agency personnel in a pre-release status as the system undergoes a formal operational readiness review. Over the past year, two projects within the Earth Science Office at NASA Marshall Space Flight Center have been investigating the performance and value of Nebula s "Infrastructure as a Service", or "IaaS" concept and applying cloud computing concepts to advance their respective mission goals. The Short-term Prediction Research and Transition (SPoRT) Center focuses on the transition of unique NASA satellite observations and weather forecasting capabilities for use within the operational forecasting community through partnerships with NOAA s National Weather Service (NWS). SPoRT has evaluated the performance of the Weather Research and Forecasting (WRF) model on virtual machines deployed within Nebula and used Nebula instances to simulate local forecasts in support of regional forecast studies of interest to select NWS forecast offices. In addition to weather forecasting applications, rapidly deployable Nebula virtual machines have supported the processing of high resolution NASA satellite imagery to support disaster assessment following the historic severe weather and tornado outbreak of April 27, 2011. Other modeling and satellite analysis activities are underway in support of NASA s SERVIR program, which integrates satellite observations, ground-based data and forecast models to monitor environmental change and improve disaster response in Central America, the Caribbean, Africa, and the Himalayas. Leveraging SPoRT s experience, SERVIR is working to establish a real-time weather forecasting model for Central America. Other modeling efforts include hydrologic forecasts for Kenya, driven by NASA satellite observations and reanalysis data sets provided by the broader meteorological community. Forecast modeling efforts are supplemented by short-term forecasts of convective initiation, determined by geostationary satellite observations processed on virtual machines powered by Nebula.

Molthan, Andrew L.↗

The SPoRT-WRF: Evaluating the Impact of NASA Datasets on Convective Forecasts

The Short-term Prediction Research and Transition Center (SPoRT) is a collaborative partnership between NASA and operational forecasting entities, including a number of National Weather Service offices. SPoRT transitions real-time NASA products and capabilities to its partners to address specific operational forecast challenges. One challenge that forecasters face is applying convection-allowing numerical models to predict mesoscale convective weather. In order to address this specific forecast challenge, SPoRT produces real-time mesoscale model forecasts using the Weather Research and Forecasting (WRF) model that includes unique NASA products and capabilities. Currently, the SPoRT configuration of the WRF model (SPoRT-WRF) incorporates the 4-km Land Information System (LIS) land surface data, 1-km SPoRT sea surface temperature analysis and 1-km Moderate resolution Imaging Spectroradiometer (MODIS) greenness vegetation fraction (GVF) analysis, and retrieved thermodynamic profiles from the Atmospheric Infrared Sounder (AIRS). The LIS, SST, and GVF data are all integrated into the SPoRT-WRF through adjustments to the initial and boundary conditions, and the AIRS data are assimilated into a 9-hour SPoRT WRF forecast each day at 0900 UTC. This study dissects the overall impact of the NASA datasets and the individual surface and atmospheric component datasets on daily mesoscale forecasts. A case study covering the super tornado outbreak across the Ce ntral and Southeastern United States during 25-27 April 2011 is examined. Three different forecasts are analyzed including the SPoRT-WRF (NASA surface and atmospheric data), the SPoRT WRF without AIRS (NASA surface data only), and the operational National Severe Storms Laboratory (NSSL) WRF (control with no NASA data). The forecasts are compared qualitatively by examining simulated versus observed radar reflectivity. Differences between the simulated reflectivity are further investigated using convective parameters along with model soundings to determine the impacts of the various NASA datasets. Additionally, quantitative evaluation of select meteorological parameters is performed using the Meteorological Evaluation Tools model verification package to compare forecasts to in situ surface and upper air observations.

Zavodsky, Bradley↗

Cloud-Resolved Simulations of Lightning NOx in an Observed Hector Thunderstorm

Cloud chemistry simulations are being performed for a "Hector" storm observed on 16 November 2005 during the SCOUT-03lACTIVE campaigns based in Darwin, Australia. The primary objective of these simulations is to estimate the average production of NO per lightning flash during the storm. The 3-D WRF-AqChem model is being used for these calculations. This modeling package contains the WRF nonhydrostatic cloud-resolving model, online gas- and aqueous-phase chemistry, and a lightning algorithm (Barth et al., 2007). Early morning soundings of temperature, water vapor and winds are used to initialize the model. Surface heating of the Tiwi Islands is simulated in the model to induce convection. Observations from the Egrett, Falcon, Geophysica, and Dornier aircraft in air undisturbed by the storm are used to construct composite initial condition chemical profiles. Convective transport in the model is tested using tracer species such as CO and O3. Lightning flashes observed by the LINET network are input to the model and a lightning placement scheme is used to inject the resulting NO into the simulated cloud. Various scenarios of NO production per flash are used for cloud-to-ground and intracloud flashes in a series of simulations for the storm. Resulting NO, mixing ratios from each simulation are compared with upper tropospheric anvil observations (from the Geophysica and Egrett aircraft) to determine the best fit with the mean NOx at anvil altitudes, the profile shape, and the frequency distribution of NOx values. We will compare the results for lightning NO production from this tropical thunderstorm with similar analyses conducted for several midlatitude and subtropical convective events.

Pickering, Kenneth↗

Representations of Precipitation Diurnal Cycle in the Amazon as Simulated by Observationally Constrained Cloud‐System Resolving and Global Climate Models

Abstract The ability of an observationally‐constrained cloud‐system resolving model (Weather Research and Forecasting; WRF, 4‐km grid spacing) and a global climate model (Energy Exascale Earth System Model; E3SM, 1‐degree grid spacing) to represent the precipitation diurnal cycle over the Amazon basin during the 2014 wet season is assessed. The WRF model coupled with a 3‐D variational data assimilation scheme reproduces the spatial variability of the precipitation diurnal cycle over the basin and the lifecycle of westward propagating MCSs initiated by the coastal sea‐breeze front. In contrast, a single morning peak in rainfall is produced by E3SM for simulations despite the nudging of large‐scale winds toward global reanalysis, indicating precipitation in E3SM is largely controlled by local convection associated with diurnal heating. The role of propagating MCS on the environment are discussed by using a multivariate perturbation analysis. We also find that the advection of moisture perturbations from ocean to inland regions have a higher correlation with the occurrence of MCSs in the Amazon than the intensity of colder air intrusion associated with sea breezes along the coast. Moreover, the presence of large cold pools over the central Amazon basin are responsible for the maintenance of propagating deep convection.

54 ENVIRONMENTAL SCIENCES↗

Modeling the impact of extreme weather events and future climate on the radiologically contaminated sites of Enewetak Atoll

Enewetak Atoll underwent 43 historical nuclear tests from 1948 to 1958, including the first hydrogen bomb test, resulting in a substantial nuclear material fallout contaminating the Atoll and the lagoon waters. The radionuclide fallout material deposited in lagoon sediments and land soil will remain for decades to come. With intensifying climate and extreme weather events, the possibility of redistribution of deposited radionuclide material has become a great concern. This study uses a numerical modeling approach to estimate the potential elevated radionuclide concentrations that can be distributed during storm events under current and future climates. We simulated three historical storm scenarios that are most likely to impact Atoll’s environment and remobilize the radionuclide-bound sediments. WRF-ARW was used to reconstruct these storm scenarios under current year (2015) and future year (2090) climates. Storm-induced ocean hydrodynamics conditions were generated using FVCOM. FVCOM-ICM was externally coupled to simulate the fate and transport of radionuclides. Given that the 239 Pu is the largest inventory of the lagoon and Atoll islands, the model results show the highest average 239 Pu concentration that an island may be exposed to is 3.25E-4 Bq/m 3 (becquerel per cubic meters), which is an increase of 84 times the average baseline/existing 239 Pu concentration without the storm conditions. The overall increase in 239 Pu average over all the islands of Atoll is about 20 folds relative to the baseline concentration. Despite the high relative increase ratios, the significantly low activity concentrations may not pose an immediate exposure risk. However, due to the limitations of the study and uncertainties/biases in the historical data used, further research supported by field surveys to better characterize the current contamination level may be needed to make more accurate predictions.

54 ENVIRONMENTAL SCIENCES↗

Influence of Lake Ice Biases in Reanalysis Data on Downscaled Climate Simulations over the Great Lakes Region

This data package contains observation-based and model-simulated datasets (all provided in NetCDF format) for evaluating how wintertime lake-ice representation affects regional weather and climate over the Laurentian Great Lakes (freshwater lake ecosystem) during the high–ice-cover winter of 2009. The observational component includes: (1) Stage IV gridded precipitation at 4 km, hourly resolution for January–February 2009 over the Great Lakes region (radar–gauge multisensor precipitation analyses); (2) Great Lakes Surface Environmental Analysis (GLSEA) satellite-derived lake-ice coverage at 1.3 km, daily resolution for the 2009 winter months, providing ice coverage over Lakes Superior, Michigan, Huron, Erie, and Ontario; and (3) in situ measurements at the Standard Rock site on Lake Superior from the Great Lakes Evaporation Network (GLEN) at hourly resolution, including near-surface atmospheric variables and sensible and latent heat fluxes (air–lake exchange) at a fixed point location. The modeling component provides corresponding fields from two simulations, both archived at 4 km, hourly resolution: a standalone Weather Research Forecasting model (WRF) run driven by the European Centre for Medium-Range Weather Forecasts (ECMWF) Reanalysis 5 (ERA5), and a two-way coupled model using WRF and the Finite Volume Community Ocean Model (WRF-FVCOM, a 3-D hydrodynamic lake model). These outputs include variables relevant to air–lake interaction and lake-effect processes (e.g., near-surface temperature, humidity, wind, precipitation, and surface turbulent fluxes), enabling direct comparison with the observational datasets. Users can analyze and visualize these NetCDF files with common tools such as Python (e.g., xarray, netCDF4, numpy, pandas), NCO/CDO, Panoply, or ncview; NetCDF variables can also be converted to other formats (e.g., CSV, GeoTIFF) using these utilities.

EARTH SCIENCE > ATMOSPHERE > ATMOSPHERIC TEMPERATU↗

Slow Wake Recovery and Low Turbulence Behind Wind Farms Parameterized in Mesoscale Simulations

Numerical weather prediction (NWP) and climate models equipped with wind-farm parameterizations (WFPs) can simulate cluster wake effects affecting downstream wind farms in both onshore and offshore environments. This study evaluates wake recovery behind a wind farm represented by the NWP-WFP approach in the Weather Research and Forecasting (WRF) model using either the Fitch et al. (2012) or Ma et al. (2022a, b) WFPs. Results are benchmarked against large-eddy simulations (LES) of an idealized offshore wind farm with aligned and staggered layouts under neutral atmospheric stability. Near-farm wake recovery is underestimated in NWP-WFP simulations due to its representation on a coarse mesoscale grid. This limitation leads to slow wake recovery through two interconnected mechanisms: (i) spatial gradients in the wind velocity field are weaker compared to LES and (ii) turbulence kinetic energy (TKE) remains low not because of excessive dissipation but due to insufficient shear production caused by these weakened gradients. For the scenario considered here, a wind-speed bias develops in the near-farm wake and persists into the far wake. Differences between the NWP-WFP simulations and LES emerge within a short distance downstream of the farm exit, where the mesoscale simulations recover too slowly. This reduced recovery contributes approximately 0.15-0.50 m s-1 to the near-farm wind-speed bias. The bias established in this region is not subsequently compensated for downstream but instead propagates into the far wake, where wind-speed differences of approximately 0.4-0.6 m s-1 remain up to 50 km downstream. Higher-resolution mesoscale simulations partially reduce this bias. Increasing turbine-added TKE or including subgrid wake effects provides additional improvement, but neither fully addresses the underlying cause. The slow wake recovery is not caused by limitations of the WFPs themselves, as it also occurs outside their region of influence, and adding subgrid wake effects does not significantly impact recovery. Rather, the slow wake recovery is a consequence of mesoscale flow representation. This behavior is not limited to regions downstream of the wind farm but is less visible within the farm, where wake recovery occurs simultaneously with turbine-induced momentum extraction. These results highlight the need for improved representations of wake recovery both within and downstream of wind farms. While enhanced subgrid modeling, shear-driven TKE production, and refined WFP formulations may improve intra-farm dynamics, accurately capturing near-farm wake recovery downstream remains challenging, as WFPs do not act in this region.

17 WIND ENERGY↗

Reply to Comment by Laprise on 'the Added Value to Global Model Projections of Climate Change by Dynamical Downscaling: a Case Study over the Continental U.S. Using the GISS-ModelE2 and WRF Models'

In his comment, Laprise raises several points that we agree merit consideration. His primary critique is that our study [Racherla et al., 2012] tested the ability of the WRF regional climate model to reproduce historical temperature and precipitation change relative to the driving global climate model (GCM) using only a single simulation rather than an ensemble. He asserts that the observed changes are smaller than the internal variability in the climate system (i.e., not statistically significant) and that thus a single simulation should not necessarily be able to capture the observations. Laprise points out that the statistical signal is reduced for a multi-decadal trend such as the one we analyzed in comparison with mean climatology and cites two studies showing that for particular climate parameters it can take any years for a signal to be discerned over internal variability. He states that The results of theexperiment as designed were strongly influenced by the presence of internal variability and sampling errors,which masked the rather small climate changes that may have occurred as a consequence of changes inforcing during the period considered. While Laprise discusses statistics in general terms at some length, for the actual climate trends examined in our study, he offers no evidence that the forced signal was smallcompared with internal variability. The two studies he cites [de Ela et al., 2013; Maraun, 2013] do not provide convincing evidence as they concern climate variables averaged over different times and areas. One in fact examines extreme precipitation events, which by definition are rare and thus have a lower significance level. We accept the general point that it is important to consider internal variability, and as noted in our paper we agree that an ensemble of simulations is in principle an optimal, though computationally expensive, approach. While we did not present the statistical significance of the observations in our original paper, we have now evaluated those for the regional temperature trends used in our study to evaluate the added value of WRF and thus can analyze data as to the magnitude of the trends with respect to internal variability.

climate↗

Influences of Cloud Microphysics on the Components of Solar Irradiance in the WRF-Solar Model

An accurate forecast of Global Horizontal solar Irradiance (GHI) and Direct Normal Irradiance (DNI) in cloudy conditions remains a major challenge in the solar energy industry. This study focuses on the impact of cloud microphysics on GHI and its partition into DNI and Diffuse Horizontal Irradiance (DHI) using the Weather Research and Forecasting model specifically designed for solar radiation applications (WRF-Solar) and seven microphysical schemes. Three stratocumulus (Sc) and five shallow cumulus (Cu) cases are simulated and evaluated against measurements at the US Department of Energy’s Atmospheric Radiation Measurement (ARM) user facility, Southern Great Plains (SGP) site. Results show that different microphysical schemes lead to spreads in simulated solar irradiance components up to 75% and 350% from their ensemble means in the Cu and Sc cases, respectively. The Cu cases have smaller microphysical sensitivity due to a limited cloud fraction and smaller domain-averaged cloud water mixing ratio compared to Sc cases. Cloud properties also influence the partition of GHI into DNI and DHI, and the model simulates better GHI than DNI and DHI due to a non-physical error compensation between DNI and DHI. The microphysical schemes that produce more accurate liquid water paths and effective radii of cloud droplets have a better overall performance.

54 ENVIRONMENTAL SCIENCES↗

Assessing Raindrop Breakup Parameterizations Using Disdrometer Observations

In this work, an intercomparison of raindrop mean diameter frequency distribution (RDFD) is performed for numerical simulations of precipitating cloud systems using an array of models and microphysics schemes. This includes results from the Regional Atmospheric Modeling System (RAMS) double-moment microphysics, the Hebrew University Cloud Model bin microphysics (HUCM) interfaced to the RAMS parent model, and the Weather Research and Forecasting (WRF) Model with the Thompson, Morrison, double-moment 6-class (WDM6), and National Severe Storms Laboratory (NSSL) double-moment schemes. Simulations are examined with respect to the raindrop size distribution (DSD) volume-number mean diameter (D m ) and intercept parameter (N w ). When compared to a suite of disdrometer observations, the RDFD resulting from each microphysics scheme exhibits varying degrees of mean drop size constraints and peaks in the frequency distribution of D m . A more detailed investigation of the peaked RDFD from the RAMS simulations suggests that the parameterization of raindrop collisional breakup can impose strong limitations on the evolution of simulated drop growth. As such, a summary and comparison of the drop breakup parameterizations among the aforementioned microphysics schemes is presented. While some drop breakup parameterizations are adjusted toward the observations by modifying the threshold diameter for the onset of breakup, this study explores the use of a modified maximum breakup efficiency. This method permits the parameterization to retain its threshold breakup diameter, while limiting the strength of drop breakup and permitting a broader range of drop sizes. As a result, the simulated mean drop sizes are in better agreement with observations.

54 ENVIRONMENTAL SCIENCES↗