SEARCH · Engineering Papers
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.
During this work, we adapted the BOSS microphysics scheme, originally formulated for rain-only, to include all warm-phase microphysical processes of cloud and rain, in particular the processes that mediate between these two categories, for example autoconversion and accretion.
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.
High resolution extended-range cloud condensation nuclei (CCN) spectral comparisons with cloud microphysics and drizzle of the Physics of Stratocumulus Tops (POST) field experiment confirmed results in the Marine Stratus/Stratocumulus Experiment (MASE). Both of these stratus cloud projects demonstrated that bimodal CCN spectra typically caused by cloud processing were associated with clouds that exhibited higher concentrations of smaller droplets with narrower distributions and less drizzle than clouds associated with unimodal CCN spectra. Resulting brighter clouds and increased cloudiness could enhance both indirect aerosol effects (IAE). These stratus findings are opposite of analogous measurements in two cumulus cloud projects, which showed bimodal CCN associated with fewer larger droplets more broadly distributed and with more drizzle than clouds associated with unimodal CCN. Resulting reduced cumulus brightness and cloudiness could reduce both IAE. Physics of Stratocumulus Tops (POST) flights in air masses with higher CCN concentrations, N CCN , showed more extremes of the stratus characteristics. However, POST flights with lower N CCN showed opposite droplet characteristics similar to the cumulus clouds, yet still showed less drizzle in clouds associated with bimodal CCN, but not as much less as the flights with higher N CCN . Since all MASE clouds were in polluted air masses, while the two cumulus projects were in clean air masses we deduce from these four projects that both the dynamic stratus/cumulus differences (vertical wind) and N CCN are responsible for the microphysics and drizzle differences among these projects. This is because the clean POST characteristics are a hybrid between MASE/POST high N CCN and the two cumulus projects.
Marine cold-air outbreak (MCAO) clouds begin as shallow cloud streets near the ice edge and evolve into broken open-cell convection downstream. Their rapid transitions in cloud structure and microphysics are closely linked to vertical motions, but this interplay has not previously been observed from space. The EarthCARE satellite, carrying the first spaceborne Doppler radar (Cloud Profiling Radar, CPR), now measures hydrometeor sedimentation velocity and vertical air motion. Using CPR observations over the Norwegian and Barents Seas (December 2024–May 2025), we analyze MCAO and non-MCAO clouds. We find that strong updrafts in MCAO clouds coincide with larger supercooled liquid water (SLW) path (LWP). These conditions are accompanied by increases in sedimentation velocity and its vertical gradient, reflecting rapid riming-driven ice growth. In contrast, non-MCAO clouds show lower riming efficiency at similar LWP due to relatively weaker dynamical support. This emphasizes the importance of dynamics in shaping polar cloud structure and microphysics.
The complexity and the small characteristic scales of aerosol microphysical processes pose a big challenge for accurate and efficient Earth system simulations at regional and global scales. In this work, we construct and evaluate a surrogate model: the aerosol deep operator network (ADON), a physics-inspired dual-net architecture for emulating the aerosol microphysics parameterization suite in the version 2 of the Energy Earth System Model (E3SMv2). The current version of the surrogate model is trained on a dataset comprising 9.8 million samples obtained from a global E3SMv2 simulation with the horizontal resolution of about one degree under cloud-free conditions. Incorporating domain spatial and temporal coordinates, as well as principle components extracted from training data, the dual-net surrogate model effectively captures the intricate representations of aerosol and the relationship with atmospheric state variables, achieving an R-squared score over $$95.7\%$$ for all the lognormal aerosol modes in the extrapolated regime. The validated model provides feature importance of input variables and their impact on the predictive capacity of the surrogate model in relation to the E3SM. The computational cost of online inference time deployed on CPUs and GPUs with lower precisions highlights ADON’s efficiency and potential in robust predictive modeling for large-scale Earth system computations.
Abstract Ice formation remains one of the most poorly represented microphysical processes in climate models. While primary ice production (PIP) parameterizations are known to have a large influence on the modeled cloud properties, the representation of secondary ice production (SIP) is incomplete and its corresponding impact is therefore largely unquantified. Furthermore, ice aggregation is another important process for the total cloud ice budget, which also remains largely unconstrained. In this study, we examine the impact of PIP, SIP, and ice aggregation on Arctic clouds, using the Norwegian Earth System Model, version 2 (NorESM2). Simulations with both prognostic and diagnostic PIP show that heterogeneous freezing alone cannot reproduce the observed cloud ice content. The implementation of missing SIP mechanisms (collisional breakup, drop shattering, and sublimation breakup) in NorESM2 improves the modeled ice properties, while improvements in liquid content occur only in simulations with prognostic PIP. However, results are sensitive to the description of collisional breakup. This mechanism, which dominates SIP in the examined conditions, is very sensitive to the treatment of the sublimation correction factor, a poorly constrained parameter that is included in the utilized parameterization. Finally, variations in ice aggregation treatment can also significantly impact cloud properties, mainly through their impact on collisional breakup efficiency. Overall, enhancement in ice production through the addition of SIP mechanisms and the reduction in ice aggregation (in line with radar observations of shallow Arctic clouds) result in enhanced cloud cover and decreased TOA radiation biases, compared to satellite measurements, especially during the cold months. Significance Statement Arctic clouds remain a large source of uncertainty in projections of the future climate due to the poor representation of the microphysical processes that govern their life cycle. Ice formation is among the least understood processes. While it is widely recognized that better constraints on primary ice production (PIP) are needed to improve existing parameterizations, we show that secondary ice production (SIP) and ice aggregation can have also a significant impact on ice number concentrations. Constraining ice formation through the addition of missing SIP mechanisms and reducing ice aggregation can improve the representation of the cloud macrophysical properties and enhance total cloud cover in the Arctic region, which in turn contributes to decreased TOA radiation biases in the cold months.
This technical report describes the Continuous Baseline Microphysical Retrieval (MICROBASE) Value-Added Product (VAP) produced operationally by the U.S. Department of Energy Atmospheric Radiation Measurement (ARM) User Facility. MICROBASE provides a continuous estimate of cloud microphysical properties at ARM fixed observatories and ARM Mobile Facility (AMF) sites. It is designed to run operationally and provide data to the ARM Data Center for scientific distribution. This technical report presents an overview of the VAP as a resource for data users and ongoing records for major updates to these products.
Abstract. Mixed-phase clouds affect precipitation and radiation differently from liquid and ice clouds, posing greater challenges to their representation in numerical simulations. Recent laboratory experiments using the Pi Cloud Chamber explored cloud glaciation conditions based on increased injection of ice-nucleating particles. In this study, we use two approaches to reproduce the results of the laboratory experiments: a bulk scalar mixing model and large-eddy simulation (LES) with bin microphysics. The first approach assumes a well-mixed domain to provide an efficient assessment of the mean cloud properties for a wide range of conditions. The second approach resolves the energy-carrying turbulence, the particle size distribution, and their spatial distribution to provide more details. These modeling approaches enable a separate and detailed examination of liquid and ice properties, which is challenging in the laboratory. Both approaches demonstrate that, with an increased ice number concentration, the flow and microphysical properties exhibit the same changes in trends. Additionally, both approaches show that the ice integral radius reaches the theoretical glaciation threshold when the cloud is subsaturated with respect to liquid water. The main difference between the results of the two approaches is that the bulk model allows for the complete glaciation of the cloud. However, LES reveals that, in a dynamic system, the cloud is not completely glaciated as liquid water droplets are continuously produced near the warm lower boundary and subsequently mixed into the chamber interior. These results highlight the importance of the ice mass fraction in distinguishing the mixed-phase clouds and ice clouds.
Ice microphysical processes are inherently complex because of their sensitivity to temperature and humidity, the diversity of ice crystal habits, and their interaction with supercooled liquid water (SCL) and turbulence. Long-term surface-based radar observations have been systematically used to unravel the different processes that affect ice particle growth. In this study, we present a statistical analysis of 6.5 years of Ka-band radar observations in Arctic cloud systems, combined with thermodynamic profiles derived from radiosonde measurements. For the first time, ice particle growth and sublimation – diagnosed from vertical gradients of radar reflectivity and mean Doppler velocity – are systematically mapped across a broad range of temperature and moisture conditions. These vertical gradients correspond closely to saturation levels relative to ice and exhibit a strong temperature dependence in supersaturated regions. Notably, distinct signatures near −15 °C are indicative of dendritic growth. Turbulence, quantified via the eddy dissipation rate (EDR), is most frequently observed in regions containing SCL. The co-occurrence of SCL and elevated turbulence results in significantly enhanced ice particle growth compared to conditions in which either is present alone. This work provides new observational constraints that are critical for improving the representation of ice microphysics in atmospheric models.
Cloud droplets containing immersed ice-nucleating particles (INPs) may freeze at temperatures above the homogeneous freezing threshold temperature in a process referred to as immersion freezing. In modeling studies, immersion freezing is often described using either so-called “singular” or “time-dependent” parameterizations. Here, we compare both approaches and discuss them in the context of probabilistic particle-based (super-droplet) cloud microphysics modeling. First, using a box model, we contrast how both parameterizations respond to idealized ambient cooling rate profiles and quantify the impact of the polydispersity of the immersed surface spectrum on the frozen fraction evolution. Presented simulations highlight that the singular approach, constituting a time-integrated form of a more general time-dependent approach, is only accurate under a limited range of ambient cooling rates. The time-dependent approach is free from this limitation. Second, using a prescribed-flow two-dimensional cloud model, we illustrate the macroscopic differences in the evolution in time of ice particle concentrations in simulations with flow regimes relevant to ambient cloud conditions. The flow-coupled aerosol-budget-resolving simulations highlight the benefits and challenges of modeling cloud condensation nuclei activation and immersion freezing on insoluble ice nuclei with super-particle methods. The challenges stem, on the one hand, from heterogeneous ice nucleation being contingent on the presence of relatively sparse immersed INPs, and on the other hand, from the need to represent a vast population of particles with relatively few so-called super particles (each representing a multiplicity of real particles). We discuss the critical role of the sampling strategy for particle attributes, including the INP size, the freezing temperature (for singular scheme) and the multiplicity.
A large convection cloud chamber has been proposed for exploring aerosol–cloud–drizzle interactions under well‐controlled turbulent conditions. Recent theoretical and numerical studies suggest that a convection cloud chamber with two heated and two cooled sidewalls can significantly enhance the liquid water content and thus benefit drizzle initiation. However, a chamber with such a sidewall configuration develops stable stratification and extremely weak turbulence therein. In this study, we conduct large‐eddy simulations of a tall convection chamber with five different sidewall configurations consisting of alternating warm and cold patches. For each configuration, the total surface area of warm patches equals that of cold patches, resulting in the same expected cloud‐free supersaturation based on a flux budget model. Results show that changing the sidewall configuration, while keeping all other factors constant, can substantially enhance turbulent mixing and improve the uniformity of thermodynamic and cloud microphysical properties in the bulk region of the chamber. In addition, turbulence strength is positively correlated with liquid water content and negatively correlated with cloud droplet number concentration, consistent with theoretical predictions. Our results highlight the advantage of building a large cloud chamber using modular patches with individually controllable temperature and humidity to achieve well‐mixed conditions.
Millimeter-wavelength radars, with their ability to collect observations and high spatiotemporal resolution, have become cornerstone instruments for investigating cloud and precipitation microphysics and dynamics in the last two decades. Since the1990s the US Department of Energy Atmospheric Radiation Measurement (ARM) program has invested in mm-wavelength radar technology progressively adding scanning capability, additional frequencies and polarimetric capability. In addition, atmospheric observatories including those technologies have been operated overseas. Those technologies have enhanced our understanding of cloud and precipitation processes. Only a few universities operate those instruments and use them for educational purposes. The central role of mm- wavelength radars and their role will be further highlighted in the next decade. NASA and other international space agencies are planning to launch millimeter-wavelength (i.e., Ku, Ka, W, G- bands) Doppler radar in space to improve our understanding of cloud and precipitation formation, evolution and their impacts on global climate.
The microphysical impacts of aerosol particles on scattered isolated deep convective cells near Houston, Texas, on 19 June 2013, are examined using multiple cloud-resolving model (CRM) simulations initialized with vertical profiles of low and high concentrations of cloud droplet–nucleating aerosols. These simulations formed part of the Model Intercomparison Project (MIP) conducted by the Deep Convective Working Group of the Aerosol, Cloud, Precipitation and Climate (ACPC) initiative. Each CRM generated a field of convective cells representing those observed during the case study with varying degrees of accuracy. The Tracking and Object-Based Analysis of Clouds (tobac) cell-tracking algorithm was applied to each MIP CRM simulation to track relatively long-lived convective cells (20–60 min). Most of the CRMs produced similar aerosol loading impacts on the warm phase of tracked cell properties with reduced autoconversion and accretion growth of rain, increased cloud water, reduced rainfall, and reduced near-surface evaporation of rain. The sign of aerosol impacts on the warm-phase properties of the convective cells was also quite consistent over cell lifetimes with the greatest magnitude of influence in the first half of the life cycle in most CRMs. In contrast, the ice-phase response to aerosol loading was highly variable among CRMs and included increases or decreases in ice amounts at inconsistent stages of the cell life cycle and midlevel versus upper-level changes in ice. This intermodel variability in ice is indicative both of the complex indirect interactions between aerosols and ice-phase processes in deep convection and their associated parameterizations.
We evaluate cloud simulations using satellite simulators against multiple observational data sets. These simulators have been run within the Geophysical Fluid Dynamics Laboratory's Atmosphere Model version 4.0 (AM4.0), as well as an alternative configuration where a fully two‐moment Morrison‐Gettelman cloud microphysical parameterization with prognostic precipitation (MG2) is applied, denoted as AM4‐MG2. The modeled cloud spatial distributions, vertical profiles, phase partitioning, cloud‐to‐precipitation transitions, and radiative effects compare reasonably well with satellite observations. Model biases include the under‐prediction of total and low‐level clouds, especially optically thin/intermediate clouds with cloud optical depth of less than 23, but the over‐prediction of thick clouds, indicating “too few, too bright” biases. These biases counteract each other, and give rise to reasonable estimates of cloud radiative effects. The underestimate of low‐level clouds is associated with too early and too frequent drizzle/precipitation formation. The precipitation bias is improved in AM4‐MG2, where the autoconversion scheme initiates the precipitation more realistically. There also exist discrepancies between models and observations for midlevel and high‐level clouds. Additional biases include the underestimate of liquid cloud fraction and the overestimate of ice cloud fraction.
Skillful representation of tropical variability and diurnal cycle of precipitation has remained a challenge in global atmosphere models, and often improvements in the variability lead to degradation in the mean‐state. Here, we introduce a configuration of the E3SM Atmosphere Model with a new large‐scale microphysics scheme and several enhancements to the deep convective scheme that improves the variability. The new configuration improves various modes of convectively‐coupled equatorial waves, with increased strength of Kelvin waves and more coherent eastward propagation of the Madden‐Julian Oscillation from the Indian Ocean to the central Pacific Ocean. The same configuration also improves the phase of the diurnal cycle of precipitation, particularly over the continental United States in the boreal summer and over Tropical land regions. Previous studies have shown that, individually taken, some of the deep convective enhancements can improve certain aspects of the variability, and here we show that combining their effects can lead to robust improvements in the variability. This model configuration can form the basis for future studies to examine the response of tropical and diurnal variability under various climate states and their relationships with other modes of variability.
Drizzle initiation through the collision and coalescence of cloud droplets plays a crucial role in warm cloud precipitation. Recent theoretical studies suggest that the influence of collisional growth on the droplet size distribution can be quantified by a non-dimensional drizzle number (Dz). Here, large-eddy simulations with Lagrangian microphysics are employed to evaluate the theory by simulating a tall convection-cloud chamber under various conditions. Results show that the smaller the Dz, the larger the impact of collisions on the right tail of the droplet size distribution, consistent with the theory. The simulations confirm that the collision rate can be estimated from the droplet size distribution interacting only with cloud droplets of the same size at the mode radius. This suggests that the idealized theory can be a useful tool to design a cloud chamber for drizzle investigation, as well as to represent drizzle formation in models of real atmospheric clouds.
Drizzle, a common feature of marine boundary layer clouds formed through collision coalescence, plays a key role in cloud microphysics and evolution. Yet, simultaneously retrieving cloud and drizzle properties from remote-sensing observations remains challenging because drizzle droplets often dominate radar signals, masking cloud contributions. The goal of the proposed research is to provide constraints for the process of autoconversion and accretion using ARM cloud measurements. Specifically, we provide concurrent retrievals of cloud and drizzle that allows users to derive corresponding autoconversion and accretion rates.
Aerosols influence the Earth's radiative balance through direct interactions with radiation and by affecting cloud properties. Anthropogenic aerosols have led to cooling during the industrial era through aerosol–cloud interactions (ACI), including aerosol effects on cloud microphysical properties and the subsequent adjustments. However, large uncertainties remain in Earth system models (ESMs) regarding the magnitude of this cooling. In part, ESMs substantially disagree on cloud properties, thermodynamics, the hydrological cycle, and general circulation. Reanalysis provides a useful avenue for exploring the impact of ACI on clouds and radiation because its atmosphere is forced to match realistic conditions through the assimilation of observations. Here, we explore the impact of ACI on clouds in the GiOcean reanalysis – the first to incorporate aerosol-cloud interactions. We contrast variables important for ACI between GiOcean and satellite observations and develop 2-dimensional lookup tables of ACI for both using a source-sink budget perspective to attribute the changes in cloud droplet number (Nd) and liquid water path (LWP) to aerosol and meteorology. A compositing analysis using lookup tables shows that GiOcean captures key aspects of aerosol–cloud–precipitation interactions, including (1) activation of aerosol into cloud droplets, (2) effective precipitation scavenging of Nd, (3) suppression of precipitation by high Nd in regions with heavy aerosol emissions. In contrast, satellite observations do not exhibit clear patterns for processes (2) and (3). Random Forest analysis shows that interannual variability in Nd and LWP over the Northern Hemisphere ocean in GiOcean is primarily driven by precipitation, consistent with satellite observations.