Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “microphysics parameterization”

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 91 records · Page 5
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Warm-phase microphysical evolution in large-eddy simulations of tropical cumulus congestus: evaluating drop size distribution evolution using polarimetry retrievals, in situ measurements, and a thermal-based framework

Owing to uncertainties in convective microphysics processes, improving parameterizations in Earth system models (ESMs) can benefit from observationally constrained cases suitable for scaling between cloud-resolving models and ESMs. We propose a benchmark large-eddy simulation (LES) cumulus congestus case study from the NASA Cloud, Aerosol, and Monsoon Processes Philippines Experiment (CAMP 2 Ex) for evaluating and improving ESMs in single-column model (SCM) mode. We seek observational constraints using novel polarimetric retrievals and in situ cloud microphysics measurements. Simulations using bulk and bin microphysics initialized with observed aerosol profiles are compared to cloud-top retrievals of cloud droplet effective radius (R eff ), effective variance (ν eff ), and number concentration (N d ) from the airborne Research Scanning Polarimeter (RSP). Both schemes reproduce characteristics of cloud-top N d and R eff that increase and decrease with altitude, respectively. Cloud-top N d is low-biased relative to RSP retrievals in both schemes, potentially due to limitations in both simulations and retrieval assumptions. Cloud-top R eff is low-biased in the bulk scheme but reasonably reproduced by the bin scheme. Profiles of N d and R eff are sensitive to the collision–coalescence process and the vertical variation in aerosol size distribution. Comparison of simulated and in situ droplet size distributions (DSDs) shows that, to first order, integrated moments are always sensitive to sizes < ~ 30 µm and can also be sensitive to larger sizes if the DSDs are sufficiently broad, with implications for the assumed maximum observed size retrieved by the RSP. The bin scheme captures the observed extended tail of the DSD, while the bulk scheme is unable to due to parametric constraints. Differences in expected relationships between in situ measurements of cloud cores and cloud-top retrievals by RSP demonstrate difficulty in constraining well the case presented herein. Finally, a thermal-tracking framework demonstrates that the dilution of N d throughout a thermal's lifetime is heavily determined by collision–coalescence and the height-varying aerosol distribution and that, in the absence of these, the impact of entrainment on diluting N d is largely offset by secondary activation. Implications for evaluating warm-phase convective microphysics schemes in ESMs and translating results for use on global, space-based polarimetry platforms are discussed.

Stanford, McKenna Wallace [Columbia Univ., New Yor↗
DOE OSTI2025

Using ARM Observations to Evaluate Process-Interactions in MCS Simulations Across Scales (Final Progress Report)

This project, funded by DOE Atmospheric System Research (DE-SC0020050), focused on improving the representation of mesoscale convective systems (MCSs) in numerical weather and climate models by leveraging high-resolution observations from the DOE Atmospheric Radiation Measurement (ARM) program. The research aimed to evaluate model sensitivities to grid spacing, microphysics, and planetary boundary layer (PBL) schemes, with a particular emphasis on improving convection parameterization for high-resolution modeling. Findings from this work highlight several key advancements. Model validation against ARM radar wind profiler data from the Southern Great Plains (SGP) and Manaus (MAO) sites revealed systematic biases in simulated convective mass flux profiles, leading to the development of an observationally constrained evaluation framework for diagnosing and improving model performance. Sensitivity analyses demonstrated that the representation of Amazonian MCSs was highly dependent on PBL scheme selection, while mid-latitude MCSs were more strongly influenced by microphysics parameterizations. A series of high-resolution WRF simulations, ranging from 4 km to 125 m grid spacing, provided insight into the behavior of convective drafts across scales. While updraft properties converged at sub-kilometer resolutions, biases in downdraft intensity persisted even at the finest resolution tested, emphasizing the need for further refinements in model physics. Additionally, comparisons of MCS vertical structures between mid-latitude and tropical environments revealed stronger updrafts and larger mass flux in mid-latitude MCSs, providing critical insights for improving climate model representations of storm-scale dynamics. The project’s findings have already contributed to advancing numerical modeling capabilities, particularly in WRF, MPAS, ICON, and DOE’s SCREAM model, by refining how convective processes are represented in high-resolution climate simulations. Results were disseminated through peer-reviewed publications, conference presentations, and ARM/ASR Research Highlights, engaging the broader scientific community. The project also provided valuable training opportunities for two postdoctoral researchers, who played central roles in model development, analysis, and dissemination of results. Their work contributed to several publications and conference presentations, helping prepare them for careers in atmospheric modeling. By improving the simulation of MCSs, this research directly supports the development of next-generation climate models capable of more accurately representing extreme precipitation and convective processes. The insights gained will inform future improvements in convective parameterization and guide the design of high-resolution weather and climate simulations, ultimately enhancing the reliability of climate projections and weather forecasts.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Tracking Aerosol Convection Interactions Experiment (TRACER) Field Campaign Report

Convective clouds serve a critical role in the Earth’s energy and water cycles through their transport of heat, moisture, momentum, and chemical species through the troposphere driving the global circulation (e.g., Hartmann et al. 1984, Del Genio et al. 2012, Su et al. 2014). On more local scales, convective clouds impact the atmospheric heating profile through diabatic heating effects, removal of water from the atmospheric column through precipitation, and conditioning of the local environment impacting further development of clouds (e.g., Sullivan and Voigt 2021). These critical roles underscore the importance of realistic representation of convective processes across scales of models from large-eddy simulation (LES), to convection-permitting models (CPM; e.g., Kendon et al. 2020, Marinescu et al. 2021), to numerical weather prediction (NWP) models used for operational weather forecasting, to Earth system models used to predict climate sensitivity (Sanderson et al. 2011, Sherwood et al. 2014, Tomassini et al. 2014, Zhao et al. 2016, Cronin et al. 2017). A key component of improving model representation of convective clouds is better quantification and parameterization of updraft microphysics and dynamics, including their interactions with the surrounding environment and storm organization (Bony et al. 2015, Hagos and Houze 2016, Donner et al. 2016, Morrison et al. 2020). Aerosol is an important environmental factor that could affect convective clouds and precipitation since cloud droplet and ice formation processes are initiated by it. Andrae et al. (2004) hypothesized that aerosols associated with increased biomass burning particles acting as cloud condensation nuclei (CCN) result in smaller and more monodisperse cloud droplets leading to suppression of warm rain formation, ultimately leading to more cloud water being lofted above the freezing level based on observations in the Amazon region. The subsequent increase in latent heat release increases the buoyancy of rising convective parcels invigorating the deep convection. This work was followed by a description of the theoretical basis for this “cold-phase invigoration” by Rosenfeld et al. (2008), who argued that it could have a significant effect for deep convective clouds with warm cloud-bases. Several modeling studies (e.g., Khain et al. 2005, 2009, van den Heever et al. 2006, Fan et al. 2007, 2009, 2012, Lee et al. 2008, Storer et al. 2010, Lebo et al. 2012, Storer and van den Heever 2013, Chen et al. 2020, Dagan et al. 2022) have investigated these aerosol-convection interactions and the environmental factors that influence their relative importance and magnitude. More recently, several studies have indicated that “warm-phase invigoration”, the enhancement of convection through condensational heating, also appears to play a role in enhancing both shallow cumuli (Seiki and Nakajima 2014, Saleeby et al 2015) and deeper tropical convection (Lebo and Seinfeld 2011, Khain et al. 2012, Sheffield et al 2015, Fan et al. 2018, Igel and van den Heever 2021), as well as Houston thunderstorms (Fan et al. 2007, 2020). However, still other studies have provided additional evidence of systematic biases in simulated convective outflow ice size distribution properties, which are consistent with a lack of poorly understood secondary ice production within convective updrafts (e.g., Fridlind et al. 2017). To help address these critical gaps in our understanding of cloud processes, aerosol processes and aerosol-cloud interactions, the Tracking Aerosol Convection Interactions Experiment was designed building upon efforts by the Aerosol, Cloud, Precipitation and Climate (ACPC) Initiative (http://acpcintiative.org/), a joint effort of the International Geosphere-Biosphere Programme (IGBP) and the World Climate Research Program (WCRP) that focused on resolving uncertainties in the interactions between aerosol and clouds towards better understanding the role that these interactions play in the climate system. The TRACER campaign was motivated by recommendations from a number of pilot studies undertaken by ACPC (van den Heever et al. 2017, Fridlind et al. 2019, Hu et al. 2019, Fan et al. 2020, Marinescu et al. 2021, Hernandez-Deckers et al. 2022) that pointed towards the southeastern Texas region as a locale where aerosol-convection interactions could be studied owing to the copious occurrence of isolated convection during the summer months accompanied by diverse and significant sources of aerosols from both anthropogenic and natural sources. The TRACER campaign began on 01 October 2021 and extended through 30 September 2022 with an intensive operational period (IOP) during June-September 2022. Three main sites (Table 1) were managed by the U.S. Department of Energy (DOE) Atmospheric Radiation Measurement (ARM) user facility.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

Examining the Impacts of Microphysical-Dynamical Feedbacks on Convective Clouds in Different Aerosol Environments Using Enhanced Observational and Modeling Strategies

This project investigated how microphysical–dynamical feedbacks influence deep convective clouds across a range of aerosol conditions, storm lifecycles, and meteorological regimes, and how these sensitivities depend on the modeling framework used to represent convection and microphysics. Using the Aerosol, Cloud, Precipitation, and Climate (ACPC) Working Group model intercomparison simulations of isolated deep convection in the Houston region, we found a robust warm-phase aerosol response in most cloud-resolving models. Increased aerosol loading tended to suppress warm-rain production, increase cloud water, and reduce surface rainfall and near-surface evaporation. These impacts were typically strongest early in the convective lifecycle, with the largest aerosol-driven differences occurring during the first half of storm evolution. The ice-phase response, on the other hand, varied widely across modeling frameworks, and differed in sign, timing, and vertical structure. Ice microphysical pathways and parameterizations are therefore leading sources of uncertainty in aerosol–deep convection interactions. Theoretical analyses further suggested that aerosol-driven invigoration through cold-phase processes was much weaker than previously hypothesized for cold-based storms, and in warm-based storms could even reduce updraft strength. These results imply that any invigoration signal is more likely linked to warm-phase processes and peaks early in storm development.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

Differences in Lightning—Storm Microphysical Relationships in Subtropical South America Versus the Continental U.S. and Their Impact on Model Lightning Parameterizations

Satellite studies have shown that some of the deepest convection with high lightning flash rates occurs downstream of the Andes in subtropical South America (SSA). Using high resolution ground‐based observations in SSA, we explore the relationship between lightning and radar‐inferred storm microphysics. Results show these storms have significantly larger 35‐dBZ volumes and more graupel per lightning flash compared to storms in the U.S. We then test microphysics‐based lightning parameterizations for storms in SSA on a 1‐km convection‐permitting simulation of a mesoscale convective system. Results from the simulation match well with observed flash rates. If parameterizations based on U.S. storms had been applied, flash rates would have been overestimated by a factor of eight. These results suggest that storms in SSA are fundamentally different compared to storms in the U.S., and caution should be exercised when applying microphysics‐based lightning parameterizations in different regions of the world.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Numerically Relevant Timescales in the MG2 Microphysics Model

Climate models rely on parameterizations of a variety of processes in the atmospheric physics, but a common concern is that the temporal resolution is too coarse to consistently resolve the behavior that individual parameterizations are designed to capture. This study examines timescales numerically derived from the Morrison-Gettelman (MG2) microphysics as implemented within the Energy Exascale Earth System Model, version 1 (E3SMv1). Numerically relevant timescales in MG2 are derived by computing the eigenspectrum of its Jacobian. These timescales are found to often be smaller than the default 5 min time step used for MG2. The fast timescales are then heuristically connected to individual microphysics processes. By substepping a few particular rain processes within MG2, the time discretization error for those processes was considerably reduced with minimal additional expense to the overall microphysics. While this improvement has a substantial effect on the target processes and on the vertical distribution of stratiform-derived rain within E3SMv1, the overall model climate is found to not be sensitive to MG2 time step. Here, we hypothesize that this is because the surface climate does not depend strongly on certain process rates, especially MG2's rain evaporation rate.

58 GEOSCIENCES↗
DOE OSTI2022

Toward Better Understanding of Microphysical Processes and Resulting Precipitation Physics: A Merger Of Observations and Cloud Models (Final Report)

A large database of global disdrometer observations and a diverse set of simulations from the Regional Atmospheric Modeling System (RAMS) were compared using principal component analysis (PCA) in order to better understand warm and ice-based precipitation processes. The analysis demonstrated that six distinct precipitation groups (PGs) with common characteristics were revealed in both the observations and model simulations. These PGs were defined by similar co-variability of rain parameters. However, the model showed large concentrations of drops with a given size compared to the observations. Detailed investigations determined the parameterization of drop breakup was forcing the drops to an equilibrium size with a stronger constraint than was indicated by the observations. A series of sensitivity studies permutating the rain shape parameter in the RAMS 2-moment bulk scheme demonstrated how the assumed shape of the rain distribution influences the microphysics as well as precipitation characteristics. Results were compared to the same simulations performed with a bin microphysics scheme where the rain shape parameter is freely evolving. In the bin scheme, a wide range of shape parameters were found, with horizontal and vertical variability, in addition to being a function of storm lifetime. The results of this study highlight the limitations of a fixed assumed rain distribution in 2-moment microphysics schemes. Model simulations were used to probe the microphysical origins of the six PGs. Rain budgets from a supercell case showed that several of the groups hypothesized to be associated with strong ice microphysical processes had significant contributions from melting hail, but the complexity of ice and warm-rain processes in the supercell precluded definitive determination of microphysical origins for many of the PGs. A case study from the Mid-latitude Continental Convection and Clouds Experiment (MC3E) compared the model PGs to disdrometer observations, and demonstrated the spatial and temporal variability of the PGs which is not possible from disdrometer point measurements alone. The analyses performed during this project demonstrated how statistical analysis can provide a robust method for comparing model simulations and observations, which ultimately resulted identification of limitations in some current microphysics parameterizations as well as insights into precipitation variability which ultimately benefits both observations and modeling efforts.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Quantifying physical parameterization uncertainties associated with land-atmosphere interactions in the WRF model over Amazon

The Weather Research and Forecasting (WRF) model can be used to diagnose regional land-atmosphere (L-A) coupling strength in the absence of sufficient observations but subjected to uncertainties associated with model physical parameterizations. In this study, we propose a framework to quantify and reduce model physical parameterization uncertainties associated with surface fluxes and L-A coupling. An ensemble of WRF simulations with different physical schemes is used to simulate surface fluxes and land-atmosphere coupling strength over the Amazon region. The physical parameterizations investigated include cloud microphysics (MP), land surface processes (LSM), planetary boundary layer (PBL), surface layer (SL), and cumulus (CU). We perform 120 ensemble simulations using the WRF model and different combinations of six MPs, three LSMs, six PBLs and SLs and three CUs. The measurements from the GoAMAZON field campaign and satellite data are used to evaluate model performance. A Multi-way analysis of variance (ANOVA) approach is applied to quantify the relative importance of different physics processes on L-A coupling. The Tukey’s test is used to sort schemes that have no significant differences into one group. The suite of physics that result in the best simulations of the corresponding variables are selected based on the Taylor skill score. Results show that the relative importance of processes and their interaction vary with the variables of interest. For example, CU was the most important process in modulating soil moisture, 2 m-humidity, latent heat, and net radiation. LSM showed dominant effects on 2 m-temperature and also has the largest impact on sensible heat and the lifting condensation level. The best physical parameterization ensembles show much narrower ranges of the variables of interest than the priori ensemble. Results of this study show the roles of different physical processes in modulating L-A interactions, quantify model uncertainties from physical processes, and provide insights for improving the model physics parameterizations.

Wang, Chen↗
DOE OSTI2024

HSW-V v1.0: localized injections of interactive volcanic aerosols and their climate impacts in a simple general circulation model

Abstract. A new set of standalone parameterizations is presented for simulating the injection, evolution, and radiative forcing by stratospheric volcanic aerosols against an idealized Held–Suarez–Williamson (HSW) atmospheric background in the Energy Exascale Earth System Model version 2 (E3SMv2). In this model configuration (HSW with enabled volcanism, HSW-V), sulfur dioxide (SO2) and ash are injected into the atmosphere with a specified profile in the vertical, and they proceed to follow a simple exponential decay. The SO2 decay is modeled as a perfect conversion to a long-living sulfate aerosol which persists in the stratosphere. All three species are implemented as tracers in the model framework and are transported by the dynamical core's advection algorithm. The aerosols contribute simultaneously to local heating of the stratosphere and cooling of the surface by a simple plane-parallel Beer–Lambert law applied on two zonally symmetric radiation broadbands in the longwave and shortwave ranges. It is shown that the implementation parameters can be tuned to produce realistic temperature anomaly signatures of large volcanic events. In particular, results are shown for an ensemble of runs that mimic the volcanic eruption of Mt. Pinatubo in 1991. The design requires no coupling to microphysical subgrid-scale parameterizations and thus approaches the computational affordability of prescribed aerosol forcing strategies. The idealized simulations contain a single isolated volcanic event against a statistically uniform climate, where no background aerosols or other sources of externally forced variability are present. HSW-V represents a simpler-to-understand tool for the development of climate source-to-impact attribution methods.

Hollowed, Joseph P. (ORCID:0000000286581672)↗
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↗
DOE OSTI2021

How Snow Aggregate Ellipsoid Shape and Orientation Variability Affects Fall Speed and Self-Aggregation Rates

Snow aggregate shapes and orientations have long been known to exhibit substantial variability. Despite this observed variability, most weather and climate prediction models use fixed power-law functions that deterministically map particle size to mass and fall speed. As such, integrated quantities like precipitation and self-aggregation rates currently ignore nonlinear effects resulting from variation in shape and orientation for aggregates of the same size. This study therefore develops an analytic framework that couples an empirically based bivariate distribution of ellipsoid shapes to classical hydrodynamic theory so as to capture an appropriate dispersion of masses, projected areas, and fall speeds for an assumed size distribution. For a fixed aggregate size, shape variations produce approximately ±0.13 m s -1 standard deviation of fall speed which increases the mass flux fall speed dispersion by more than 100% over traditional microphysics models. This increased fall speed dispersion results predominantly from shape-induced mass dispersion whereas orientation and drag dispersion play a lesser role. Shape variations can increase mass- and reflectivity-weighted fall speeds by up to 60% of traditional models whereas self-aggregation rates can increase by a factor of 100 for very small slope parameters. This implies that aggregate shape variations effectively forestall the theorized onset of fall speed distribution narrowing and subsequent quenching of the aggregation process. As a result, it is likely that secondary ice formation is necessary to prevent an ever decreasing slope parameter. The mathematical theory presented in this study is used to develop simple correction factors for snow forecast and climate models.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Radar-Based Bayesian Estimation of Ice Crystal Growth Parameters within a Microphysical Model

The potential for polarimetric Doppler radar measurements to improve predictions of ice microphysical processes within an idealized model–observational framework is examined. In an effort to more rigorously constrain ice growth processes (e.g., vapor deposition) with observations of natural clouds, a novel framework is developed to compare simulated and observed radar measurements, coupling a bulk adaptive-habit model of vapor growth to a polarimetric radar forward model. Bayesian inference on key microphysical model parameters is then used, via a Markov chain Monte Carlo sampler, to estimate the probability distribution of the model parameters. The statistical formalism of this method allows for robust estimates of the optimal parameter values, along with (non-Gaussian) estimates of their uncertainty. To demonstrate this framework, observations from Department of Energy radars in the Arctic during a case of pristine ice precipitation are used to constrain vapor deposition parameters in the adaptive habit model. The resulting parameter probability distributions provide physically plausible changes in ice particle density and aspect ratio during growth. A lack of direct constraint on the number concentration produces a range of possible mean particle sizes, with the mean size inversely correlated to number concentration. Consistency is found between the estimated inherent growth ratio and independent laboratory measurements, increasing confidence in the parameter PDFs and demonstrating the effectiveness of the radar measurements in constraining the parameters. Furthermore, the combined Doppler and polarimetric observations produce the highest-confidence estimates of the parameter PDFs, with the Doppler measurements providing a stronger constraint for this case.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Fragmentation in Collisions of Snow with Graupel/Hail: New Formulation from Field Observations

Abstract Secondary ice production (SIP) has been attributed to the generation of most ice particles observed in precipitating clouds with cloud tops warmer than −36°C, from various aircraft- and ground-based field observations across the globe. One of the known SIP mechanisms is fragmentation during collisions among ice particles. It has been studied with our theoretical formulation, which has been applied in microphysical schemes of atmospheric models in a few studies. These have predicted an extensive impact on cloud glaciation and radiative properties. However, there has been a lack of experimental field studies, especially involving naturally falling snowflakes, to better understand this particular mechanism of SIP. This study reports the first field measurements with modern technology for fragmentation during collisions between naturally falling snowflakes and graupel/hail particles. This was observed with an innovatively designed portable chamber that was deployed outdoors in northern Sweden. Applying the observations from this field-based study, we optimized the existing formulation for predicting numbers of fragments from collisions of snow with graupel/hail. The observations show the average numbers of fragments per collision for dendritic (3–12 mm) and nondendritic (1–3 mm) snow were about 12 and 1, respectively. This represents a boost of predicted fragment numbers relative to our original formulation published in 2017. The updated formulation for breakup in ice–ice collisions can be implemented in the microphysical schemes of atmospheric models.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Estimating Drizzle Parameters by Aircraft Sampling

In-cloud supersaturation is calculated from aircraft derived measurements of the in-cloud vertical gradient of liquid water content, vertical velocity, the pt and 3 rd cloud droplet radius moments and a parameter related to the heat and mass transport during condensation and evaporation of cloud droplets. The aircraft data where taken during the VOCALS 2008 study in marine stratocumulus clouds off the coast of Chile. Approximately 16000 discrete samples, 2.5 meter in length, were used in the present analysis of in-cloud supersaturation and its variability, both of which are indicative of turbulent processes within the cloud and their significant impact on drizzle formation. Two additional drizzle parameters, the eddy persistence time and the variance of fluctuations in in-cloud droplet growth rate were also estimated from this data set. These results will help to constrain parameterizations of underlying turbulent microphysical processes of clouds and thereby advance the future development of warm cloud and precipitation models for climate simulation.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Improving Aerosol Radiative Forcing and Climate in E3SM: Impacts of New Cloud Microphysics and Improved Wet Removal Treatments

Abstract Numerous Earth system models exhibit excessive aerosol effective forcing at the top of the atmosphere (TOA), including the Department of Energy's Energy Exascale Earth System Model (E3SM). Here, in the context of the E3SM version 3 effort, the predicted particle property (P3) stratiform cloud microphysics scheme and an enhanced deep convection parameterization suite (ZM_plus) are implemented into E3SM. The ZM_plus includes a convective cloud microphysics scheme, a multi‐scale coherent structure parameterization for mesoscale convective systems, and a revised cloud base mass flux formulation considering impacts of the large‐scale environment. The P3 scheme improved cloud and radiation particularly over the Northern Hemisphere and the frequency of heavy precipitation over the tropics, and the ZM_plus improved clouds in the tropics. P3 decreases aerosol effective forcing by 0.15 W m −2 , while the ZM_plus increases it by 0.27 W m −2 , resulting from excessive direct (0.31 W m −2 ) and indirect forcing (−1.79 W m −2 ). The excessive aerosol forcings are due to aerosol overestimation associated with insufficient aerosol wet removal. By improving the physical treatments in the aerosol wet removal, we effectively mitigate anthropogenic aerosol overestimation and thus attenuate direct (0.09 W m −2 ) and indirect aerosol forcing (−1.52 W m −2 ). Adjustment to primary organic matter hygroscopicity reduces direct and indirect forcing to more reasonable values: −0.13 W m −2 and −1.31 W m −2 , respectively. On climatology, improved aerosol treatments mitigate overestimation of aerosol optical depth.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Using Machine Learning to Predict Cloud Turbulent Entrainment–Mixing Processes

Different turbulent entrainment–mixing mechanisms between clouds and environment are essential to cloud–related processes; however, accurate representation of entrainment–mixing in weather/climate models still poses a challenge. This study exploits the use of machine learning (ML) to address this challenge. Four ML (Light Gradient Boosting Machine [LGB], eXtreme Gradient Boosting, Random Forest, and Support Vector Regression) are examined and compared. It is found that LGB performs best, and thus is selected to understand the impact of entrainment–mixing on microphysics using simulation data from Explicit Mixing Parcel Model. Compared with traditional parameterizations, the trained LGB provides more accurate microphysical properties (number concentration and cloud droplet spectral dispersion). The partial dependences of predicted microphysics on features exhibit a strong alignment with physical mechanisms and expectations, as determined by the interpreting method, thus overcoming the limitations of the “black box” scheme. The underlying mechanisms are that the smaller number concentration and larger spectral dispersion correspond to more inhomogeneous entrainment–mixing. Specifically, number concentration after entrainment–mixing is positively correlated with adiabatic number concentration and liquid water content affected by entrainment–mixing, and inversely correlated with adiabatic volume mean radius. Spectral dispersion after entrainment–mixing is negatively correlated with liquid water content affected by entrainment–mixing, turbulent dissipation rate and relative humidity of entrained air. Sensitivity analysis further suggests that number concentration is mainly determined by cloud microphysical properties whereas spectral dispersion is influenced by both cloud microphysical properties and environmental variables. The results indicate that the LGB scheme has the potential to enhance the representation of entrainment–mixing in weather/climate models.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Understanding the Microphysical Control and Spatial‐Temporal Variability of Warm Rain Probability Using CloudSat and MODIS Observations

Abstract By combining measurements from MODIS and the CloudSat radar, we develop a parameterization scheme to quantify the combined microphysical controls by liquid water path (LWP) and cloud droplet number concentration (CDNC) of the probability of precipitation (PoP) in marine low cloud over tropical oceans. We demonstrate that the spatial‐temporal variation of grid‐mean in‐cloud can be largely explained by the variation of the joint probability density function of LWP and CDNC in the phase space specified by the bivariate PoP (LWP and CDNC) function. Through a series of sensitivity tests guided by this understanding, we find that in the Southeastern Pacific and Atlantic the stratocumulus to cumulus transition of the is mainly due to the variation of CDNC while the annual cycle is mainly due to the variation of LWP. The results of this study provide a viable way to diagnose the root cause of warm rain problems in global climate models.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Formulation of autoconversion and drop spectra shape in shallow cumulus clouds

Two-moment autoconversion parameterizations as compared to accretion parameterizations exhibit significant errors suggesting that additional moments are needed to increase their accuracy. We develop a three-moment autoconversion parameterization using output from an LES model with size-resolved microphysics. Adding the third moment decreases the errors of parameterization and improves precipitation prediction. However, the errors are still significantly larger than errors of accretion rate. An analysis of the cloud drop size distributions (DSDs) in the simulated tropical convective cloud system reveals that most of DSDs have a significant fraction of cloud liquid water content (qc) in the mid-size droplet range (radii from 20 to 40 microns). Our data indicates that more than 30% of DSDs have over half of qc contained in the mid-size range and about 60% of spectra have, at least, one third of qc in this range. Even when the rain/drizzle mode is small (radar reflectivity Z < -10 dBZ), there is a significant number of spectra in which fraction of qc in the mid-size range is as large as 60%. These DSDs are more complex than the frequently used single-mode Gamma or Log-normal distributions, which usually do not have a considerable tail extending to the mid-size range and can be defined by three microphysical moments. The need to define DSDs by more than three moments explains the large errors in the three-moment autoconversion parameterization. The limitation of three-parameter Gamma or Log-normal distributions should be kept in mind when applying them in precipitating shallow Cu clouds.

Kogan, Yefim↗