Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Ice Microphysical Processes”

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 37 records · Page 2
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

NSA Site Science: Use of ARM Observations from Northern Alaska to Evaluate and Improve Prediction Capabilities

The Arctic is warming at a rate nearly double that of the rest of the planet, leading to profound changes in atmospheric, oceanic, and ice processes. The U.S. Department of Energy's (DOE) Atmospheric Radiation Measurement (ARM) user facility has played a significant role in Arctic research, operating observatories in Alaska's North Slope for over 25 years. These observatories provide a rich, and wide-reaching dataset that offers insight into atmospheric processes in northern Alaska. This report details the results of a nine-year research project (2015–2024) supported by the DOE Atmospheric Systems Research (ASR) program, that leverages data from ARM’s deployment of observing facilities at Utqiaġvik (known as the North Slope of Alaska, or NSA, site) and Oliktok Point, Alaska. The project was conducted in two phases: - Phase 1 (2015–2019): Focused on understanding key atmospheric processes at Oliktok Point, including cloud formation, high-latitude precipitation, aerosol-cloud interactions, and cloud properties. - Phase 2 (2019–2024): Extended the research to the broader North Slope region, using data from both Oliktok Point and Utqiaġvik. Topics explored included surface energy budgets, atmospheric stability, ice nucleation processes, and microphysics in Arctic clouds. The project resulted in numerous research products, including 50 peer-reviewed publications and dissertations, 169 presentations, and 10 data products. These products cover a variety of topics, including: - Cloud Macro- and Microphysical Properties: Arctic clouds play a crucial role in energy transfer, and accurate representation in models is critical. The study explored cloud transitions, ice crystal shapes, and dual-wavelength radar data to understand ice crystal habits and size distributions. - Aerosol Properties and Processes: The team examined aerosol sources in the Arctic, including industrial emissions and natural sources. Observations showed significant spatial gradients in aerosol concentrations due to human activities and wildfire smoke. The influence of aerosols on cloud formation and the surface energy budget was also assessed. - Aerosol-Cloud Interactions: Research revealed that aerosols might suppress cloud ice production, affecting cloud radiative forcing and precipitation. The impact of local industrial emissions on cloud properties was also investigated. - Contextualizing the North Slope of Alaska in the context of the broader Arctic: To understand broader trends, the project evaluated large-scale circulation patterns and the influence of weather systems on the Arctic. Studies indicated that large-scale processes play a significant role in temperature patterns and the timing of snowmelt. - Advancing ARM Observational and Modeling Capabilities: The project developed new radar data products and advanced measurement techniques, including clutter mitigation and drizzle detection. Uncrewed aerial systems (UAS) and tethered balloon systems (TBS) were deployed to gather detailed atmospheric data. Additionally, the project supported 10 early career scientists, providing training and mentorship to undergraduate interns, graduate students, postdoctoral researchers, and early career researchers. These efforts contributed to the advancement of ARM research capabilities and fostered a new generation of scientists skilled in Arctic atmospheric research. Ultimately, this ASR-supported project has provided valuable insights into Arctic atmospheric processes and their broader climate implications. Recommendations for future work include continuing support for long-term observing at Arctic locations to foster additional research, further exploration of aerosol-cloud interactions and the potential impacts of enhanced industrialization of the Arctic, and expanded use of uncrewed systems to gather data in this remote and harsh environment. Additionally, the data products developed by this work, and the data products developed through the ARM infrastructure, leave a treasure-trove of additional information that should be explored for many years to come to gain additional insight into physical processes in the Arctic atmosphere that drive the rapid changes occurring in at high latitudes and their global impact.

58 GEOSCIENCES↗
DOE OSTI2021

Tropical Cirrus in Global Storm-Resolving Models. Part II: Cirrus Life Cycle and Top-of-Atmosphere Radiative Fluxes

Cirrus clouds of various thicknesses and radiative characteristics extend over much of the tropics, especially around deep convection. They are difficult to observe due to their high altitude and sometimes small optical depths. They are also difficult to simulate in conventional global climate models, which have coarse grid spacings and simplified parameterizations of deep convection and cirrus formation. We investigate the representation of tropical cirrus in global storm-resolving models (GSRMs), which have higher spatial resolution and explicit convection and could more accurately represent cirrus cloud processes. This study uses GSRMs from the DYnamics of the Atmospheric general circulation Modeled On Non-hydrostatic Domains (DYAMOND) project. The aggregate life cycle of tropical cirrus is analyzed using joint albedo and outgoing longwave radiation (OLR) histograms to assess the fidelity of models in capturing the observed cirrus cloud populations over representative tropical ocean and land regions. The proportions of optically-thick deep convection, anvils, and cirrus vary across models and are portrayed in the vertical distribution of cloud cover and top-of-atmosphere radiative fluxes. Model differences in cirrus populations, likely driven by subgrid processes such as ice microphysics, dominate over regional differences between convectively-active tropical land and ocean locations.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

The microphysics of the warm-rain and ice crystal processes of precipitation in simulated continental convective storms

Precipitation in clouds can form by either warm-rain or ice crystal processes, referred to as warm and cold formation pathways, respectively. Here, we investigate the warm and cold pathway contributions to surface precipitation in simulated continental convective storms. We analyze three contrasting convective storms that are cold-based, slightly warm-based and very warm-based. We apply tracer-tagging techniques in our aerosol-cloud model to determine simulated microphysical pathways that lead to precipitation. We find cold components of graupel and rain mass were higher than warm components in cold- and slightly warm-based clouds. By contrast, in very warm-based clouds nearly 80% of surface precipitation was formed via warm-rain processes. Lowering of cloud base altitude to levels about 10–20 K warmer switched surface precipitation to being mostly warm, due to enhanced moisture content in the planetary boundary layer and larger cloud droplets aloft intensifying raindrop freezing. Our simulations indicate that warm and cold processes co-exist in any storm and the balance between them is determined by cloud base temperature and solute aerosol conditions.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Implications of Aerosol Physicochemical Properties Including Ice Nucleation at ARM Mega Sites for Improved Understanding of Microphysical Atmospheric Cloud Processes (Final Technical Report)

Continuous, long-term measurements of atmospheric ice-nucleating particles (INPs) that influence clouds and precipitation were conducted to investigate the abundance and variability of ground-level INPs across the world. Three field campaigns were organized by the DOE Atmospheric Radiation Measurement (ARM) program, including Examining INP from Southern Great Plains (ExINP-SGP, 2019), Examining INP from Eastern North Atlantic (ExINP-ENA, 2020 – 2021), and Examining INP from North Slope of Alaska (ExINP-NSA, 2021 – 2024). Additional small-scale supporting field experiments were performed in 2019 and 2021 to collect airborne particulate matter at SGP for complementary laboratory characterization of the particles’ physical and chemical properties [Aerosol-Ice Formation Closure Pilot Study (AEROICESTUDY), 2019; ExINP-SGP II, 2021]. In these studies, the PI’s team measured INP concentration with both real-time and laboratory measurements in a wide range of freezing temperatures ($T$ from 0 °C to about –30 °C). This project elucidated spatial variability and seasonality in the abundance of immersion mode active INPs across three ARM sites using a single instrument, a Portable Ice Nucleation Experiment (PINE) chamber version 03 (PINE-03 hereafter). Collocated aerosol and meteorological data were analyzed to assess the correlation between ambient INP abundance, air mass origin region, and meteorological variability. Our findings suggest very high freezing efficiency of INPs at the NSA site across the measured temperatures (ice nucleation active surface site density, $n_s(T)$, $\approx 2 \times 10^{8} - 10^{10}$ m -2 for from –16 to –31 °C), which is a factor of 10 – 1000 times greater efficiency as compared to that found in the previous mid-latitude INP measurements in autumn using the same instrument; surprisingly high INP abundance ($\ge 1 \text{ L}^{-1}$ at –25 °C) for the examined temperatures throughout the year that PINE-03 did not measure at other sites; and high INP concentration in spring, possibly related to arctic haze episodes.

58 GEOSCIENCES↗
DOE OSTI2022

Development of a triple-moment ice-phase cloud microphysics scheme and its application to the Single Column Atmosphere Model

Parameterization of cloud microphysics is critical for accurate simulation of weather and climate, in which the characteristics of cloud particle spectrum inevitably further affect climate simulation by changing the cloud evolution and cloud radiation effects. The popular currently used double-moment cloud microphysics schemes in numerical models can predict the intercept (N 0 ) and slope (λ) parameters of cloud particle spectrum but cannot predict the spectral shape parameter (μ), which hinders accurate description of cloud physical processes in climate models. Therefore, in the present study, we built upon the ideas of previously developed triple-moment cloud microphysics scheme, considered radar reflectivity factor as the third predictor in addition to number and mass concentration, and introduced its related prediction equation to the Single Column Atmosphere Model Version 5.3 (SCAM5.3). Moreover, the relevant microphysical process formulas in the model were revised, and a triple-moment ice-cloud microphysics scheme was constructed to predict the μ of ice particle spectrum. Based on this model, a 29-day case of the Atmospheric Radiation Measurement Program in the summer of 1997 (ARM97) was simulated, and the differences in cloud fraction and radiation simulations between the double- and triple-moment schemes were analysed. The μ of ice particle spectrum predicted by the triple-moment scheme mainly ranged from 0 to 4; the peak value was around 2, and at least 85% of the μ values were greater than the default 0. Therefore, the μ=0 setting in the double-moment scheme is unreasonable. Moreover, the developed triple-moment ice-cloud microphysics scheme yielded a narrower ice particle spectrum, which was closer to the observation results of previous studies, than the double-moment cloud microphysics scheme (e.g., height=233 hPa). Furthermore, compared with the double-moment scheme, the triple-moment scheme achieved closer simulations of the cloud fraction observations, particularly for ice clouds in the upper levels. Considering the close association between cloud fraction and radiation, the error with the triple-moment scheme was smaller than that with the double-moment scheme regardless of the shortwave (surface downward shortwave radiation, net downward shortwave radiation at the top of the atmosphere, and net shortwave radiation at the surface) or longwave (surface downward longwave radiation and net longwave radiation at the surface) flux density. Finally, the improvement mechanism of the triple-moment scheme on cloud fraction and radiation simulations was explored. The major reason is that the triple-moment scheme weakens the autoconversion of ice crystal to form snow and enhances the growth process of ice crystal deposition, thereby increasing the mass concentration of ice crystals and ice cloud fraction. Overall, the developed triple-moment ice-cloud microphysics scheme can help improve the simulation ability of models for cloud-climate feedback and other processes.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

The Energy Exascale Earth System Model Version 3: 1. Overview of the Atmospheric Component

This paper describes the atmospheric component of the US Department of Energy's Energy Exascale Earth System Model (E3SM) version 3. Significant updates have been made to the atmospheric physics compared to earlier versions. Specifically, interactive gas chemistry has been implemented, along with improved representations of aerosols and dust emissions. A new stratiform cloud microphysics scheme more physically treats ice processes and aerosol‐cloud interactions. The deep convection parameterization has been largely improved with sophisticated microphysics for convective clouds, making model convection sensitive to large‐scale dynamics, and incorporating the dynamical and physical effects of organized mesoscale convection. Improvements in aerosol wet removal processes and parameter re‐tuning of key aerosol and cloud processes have improved model aerosol radiative forcing. The model's vertical resolution has increased from 72 to 80 layers with the extra eight layers added in the lower stratosphere to better simulate the Quasi‐Biennial Oscillation. These improvements have enhanced E3SM's capability to couple aerosol, chemistry, and biogeochemistry and reduced some long‐standing biases in simulating tropical variability. Compared to its predecessors, the model shows a much stronger signal for the Madden‐Julian Oscillation, Kelvin waves, mixed Rossby‐gravity waves, and eastward inertia‐gravity waves. Aerosol radiative forcing has been considerably reduced and is now better aligned with community best estimates, leading to significantly improved skill in simulating historical temperature records. Its simulated mean‐state climate is largely comparable to E3SMv2, but with some notable degradation in shortwave cloud radiative effect, precipitable water, and surface wind stress, which will be addressed in future updates.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Toward Understanding the Simulated Phase Partitioning of Arctic Single-Layer Mixed-Phase Clouds in E3SM

21 Mixed-phase clouds are widespread in the Arctic and have substantial impacts on the 22 surface energy budget, sea ice extent, and regional climate. Significant changes are found in the 23 modeled phase partitioning of Arctic mixed-phase clouds in the newly developed U.S. 24 Department of Energy (DOE) Energy Exascale Earth System Model (E3SM) Atmosphere Model 25 version 1 (EAMv1) compared to its predecessor (CAM5/EAMv0). The mixed-phase clouds 26 simulated by EAMv1 is overly dominated by super-cool liquid which is often underestimated by 27 CAM5. Compared to CAM5, EAMv1 has made several important changes in its cloud physical 28 parameterizations, including using CLUBB to unify shallow convection, PBL turbulence, and 29 cloud macrophysics, as well as updating the MG two-moment cloud microphysics scheme and 30 using the CNT heterogeneous ice nucleation scheme to replace the Meyers scheme for mixed- 31 phase clouds. To understand which of the newly-implemented physical parameterizations is 32 responsible for the change of modeled mixed-phase cloud phase partitioning in EAMv1, 33 sensitivity tests are conducted by utilizing the short-term hindcast approach and results are 34 compared to the observations from the ARM Mixed-Phase Arctic Cloud Experiment (M-PACE) 35 field campaign. 36 It is found that although each scheme contributes to the decrease of cloud ice water 37 content (IWC) and increase of cloud liquid water content (LWC), only the combined effect of 38 CLUBB and CNT can explain the significant underestimation of IWC. By examining the budgets 39 of cloud microphysical processes, CLUBB is found to play a dominant role due to the fact that 40 CLUBB does not diagnoses cloud ice phase water in its assumed joint PDF method and thus may 41 not provide sufficient cloud ice to initialize ice phase microphysical processes such as the growth 42 of ice crystals through vapor deposition. Although a more effective heterogeneous ice nucleation 43 may alleviate the underestimation of IWC, the ice phase process in CLUBB requires urgent 44 attention in the future model development.

Zhang, Meng↗
DOE OSTI2022

High concentrations of ice: Investigations using polarimetric radar observations combined with in situ measurements and cloud modeling (Final Report)

This DOE-funded joint project had the over-arching goal of understanding the reasons for observed size distributions of ice particles in cold clouds. Its approach involves the use of cloud models and field observations by radar and aircraft of real storms. The project addresses three major research objectives: (1) Utilize a novel polarimetric radar technique to retrieve size distributions and amounts of ice from radar data collected during the previous DOE ARM field campaigns; (2) Evaluate the results of ice microphysical retrievals using available in situ aircraft measurements and other remote sensors (e.g., vertically pointing cloud radars and wind profilers); (3) Improve the treatment of microphysical processes leading to cloud glaciation and ice multiplication in numerical cloud models. The study was performed by three research organizations: University of Oklahoma, The Hebrew University of Jerusalem, Israel, and Lund University, Sweden. The Israel and Swedish partners were funded by the sub-awards from the University of Oklahoma. Cloud modeling studies have been performed by the research teams at The Hebrew University of Jerusalem (HUJ) and Lund University to identify the origins of high concentrations of cloud ice in areas of high ice water content (HIWC). The Hebrew University Cloud Model (HUCM) with full spectral bin microphysics and the Lund University aerosol-cloud (AC) model with a hybrid bin / bulk microphysics scheme complementing HUCM were utilized for simulations. Both research teams had particular focus on secondary ice production (SIP) as one of the possible sources of enhanced ice concentration. The HUJ group suggested a novel concept of ice multiplication during droplet freezing. It is assumed that splintering and droplet fragmentation during droplet freezing takes place because of dendritic growth within a supercooled drop. The resulting simulations of SIP generated small ice in concentrations exceeding hundreds per liter similar to what was observed in the HIWC regions of the tropical storms. The Lund team explored the SIP mechanisms such as breakup of ice particles due to ice-ice collisions and ice sublimation that are expected to dominate the continental storms. They also quantified the impact of homogeneous nucleation of cloud droplets on the total number concentration of ice at very low temperatures near the tops of the clouds. Additionally, the Lund AC model is able to simulate the effect of aerosols of various types (including biological) on the cloud life cycle and the corresponding ice production. The HUCM / AC model was used to simulate one of the “golden” cases of the DOE MC3E campaign on 20 May 2011. The model output was converted into the fields of polarimetric radar variables using the polarimetric radar forward operator developed at the University of Oklahoma and compared with radar observations and in situ microphysical measurements onboard research aircraft. It was demonstrated that specific differential phase KDP is the best radar parameter to identify the HIWC areas and quantify the corresponding ice parameters. For the first time, the shape of the vertical profile of KDP was realistically reproduced by the cloud model with the KDP maximum in the dendritic growth layer (DGL) centered at the -15°C isotherm. The University of Oklahoma team has developed a methodology for polarimetric radar retrievals of such microphysical parameters of ice as ice water content (IWC), mean volume diameter (Dm), and total number concentration (Nt) of ice particles. These retrievals have been validated using in situ aircraft measurements during 6 field campaigns and proved to be quite robust and reliable. This allowed to build the first climatology of the vertical profiles of polarimetric radar variables and retrieved microphysical parameters for the three types of weather systems: continental MCSs, maritime MCSs, and tropical cyclones / hurricanes (Hu and Ryzhkov 2022). The data were collected by a multitude of the WSR-88D radars in 13 continental and 10 maritime MCSs and 11 landfalling hurricanes. The HIWC areas were identified within the examined storms and the corresponding “HIWC statistics” was compared with the “background” one without HIWC. An overarching conclusion of the study is that maritime tropical storms (MCSs and hurricanes) are characterized by smaller size ice in higher concentration compared to the continental MCSs. High ice water content in the HIWC areas is primarily caused by a strong jump in a number concentration of ice particles rather than the increase of their size compared to the “background” environment. This may point to the homogeneous nucleation of excessive amounts of supercooled droplets and / or secondary ice production as the possible origins of HIWC. Such a climatology provides a good observational reference for the modelers to evaluate the performance of their models. As an example, the in-depth analysis of the 20 May 2011 MC3E case shows that the advanced cloud models developed in the course of this study still tend to underestimate the number concentration of ice in the HIWC areas although they succeed in reproducing realistically looking vertical profiles of IWC and Nt. The results of the project research are summarized in 13 journal papers.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Multiple Environmental Influences on the Lightning of Cold-Based Continental Convection. Part II: Sensitivity Tests for Its Charge Structure and Land–Ocean Contrast

In Part I, an electrification scheme was described and a simulation of an observed cold-based storm from the U.S. Great Plains was validated with electrical observations. Most charge in the storm was separated by rebounding collisions of secondary ice originating from prior graupel–snow collisions. In this Part II, sensitivity tests are performed with the control simulation (Part I) and influences from environmental factors (aerosols, temperature, moisture, and shear) on lightning are elucidated. Environmental factors [e.g., convective available potential energy (CAPE)] controlling updraft speed are salient. When ascent is reduced by 30% and 70%, flashes become 70% fewer and disappear, respectively; faster ascent promotes positive cloud-to-ground (+CGs) flashes. Since cloud base is too cold (1°C) for coalescence, cloud condensation nucleus aerosol concentrations do not influence the lightning appreciably. The electrical response to varying concentrations of active ice nuclei is limited by most ice particles being secondary and less sensitive—a natural “buffer.” Imposing a maritime sounding suggests that the land–sea contrast in lightning for such storms is due to the vertical structure of environmental temperature and humidity. Weak CAPE, and both entrainment and condensate weight from a low cloud base, suppress ascent and charging. Maritime thermodynamic conditions reduce simulated flash rates by two orders of magnitude. Reducing aerosol loadings from continental to maritime only slightly reinforces this suppression. Last, a conceptual model is provided for how any simulated storm is either normal because graupel/hail is mostly positively charged or else is inverted/anomalous because graupel/hail is mostly negatively charged, with environmental factors controlling the charging. Impacts from microphysical processes, including three processes of secondary ice production, on lightning are analyzed.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Analysis of the microphysical properties of snowfall using scanning polarimetric and vertically pointing multi-frequency Doppler radars

Radar dual-wavelength ratio (DWR) measurements from the Stony Brook Radar Observatory Ka-band scanning polarimetric radar (KASPR, 35 GHz), a W-band profiling radar (94 GHz), and a next-generation K-band (24 GHz) micro rain radar (MRRPro) were exploited for ice particle identification using triple-frequency approaches. The results indicated that two of the radar frequencies (K and Ka band) are not sufficiently separated; thus, the triple-frequency radar approaches had limited success. On the other hand, a joint analysis of DWR, mean Doppler velocity (MDV), and polarimetric radar variables indicated potential in identifying ice particle types and distinguishing among different ice growth processes and even in revealing additional microphysical details. We investigated all DWR pairs in conjunction with MDV from the KASPR profiling measurements and differential reflectivity (Z DR ) and specific differential phase (K DP ) from the KASPR quasi-vertical profiles. The DWR-versus-MDV diagrams coupled with the polarimetric observables exhibited distinct separations of particle populations attributed to different rime degrees and particle growth processes. In fallstreaks, the 35–94 GHz DWR pair increased with the magnitude of MDV corresponding to the scattering calculations for aggregates with lower degrees of riming. The DWR values further increased at lower altitudes while Z DR slightly decreased, indicating further aggregation. Particle populations with higher rime degrees had a similar increase in DWR but a 1–1.5 m s –1 larger magnitude of MDV and rapid decreases in K DP and Z DR . The analysis also depicted the early stage of riming where Z DR increased with the MDV magnitude collocated with small increases in DWR. This approach will improve quantitative estimations of snow amount and microphysical quantities such as rime mass fraction. The study suggests that triple-frequency measurements are not always necessary for in-depth ice microphysical studies and that dual-frequency polarimetric and Doppler measurements can successfully be used to gain insights into ice hydrometeor microphysics.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Ice nucleation by volcanic ash greatly alters cirrus cloud properties

The formation of ice crystals in the atmosphere strongly affects cloud properties and climate. While volcanic ash (VA) has been shown to nucleate ice crystals efficiently in laboratory settings, its importance for ice formation in the atmosphere remains elusive. Here, we show evidence of cirrus modification by volcanic eruptions through ice nucleation on VA, revealed by abrupt changes in cirrus properties following volcanic eruptions based on satellite measurements. The distinct changes captured are a phenomenal decrease in number, an increase in size of ice crystals in cirrus clouds, and an increase in cirrus occurrences after ash-rich volcanic eruptions. Conversely, no such changes were detected following the ash-poor eruption. We propose a cirrus formation mechanism where VA nucleates ice heterogeneously, suppressing homogeneous freezing and resulting in fewer but larger ice crystals. This suppression of homogeneous freezing by VA is supported by process-level cloud microphysical simulations. Our findings advance the understanding of aerosol–ice cloud interactions and illuminate cirrus geoengineering.

Lin, Lin [Texas A & M Univ., College Station, TX (↗
DOE OSTI2026

Prognostic simulations of mixed-phase clouds with model AC-1D v1.0: the impact of aerosol types and freezing parameterizations on ice crystal budgets

Mixed-phase clouds at high latitudes contribute to the uncertainty in predicting cloud feedbacks and climate sensitivity, mainly due to the complexity of microphysical processes that influence the partitioning between the supercooled liquid and ice phases, and hence, cloud radiative effects on regional scales. Particularly in Arctic mixed-phase clouds, the activation of ice-nucleating particles (INPs) from various aerosol populations remains a leading source of uncertainty. We developed an aerosol-cloud one-dimensional (AC-1D) model, which provides a novel framework to prognostically treat INP and ice crystal budgets while explicitly accounting for polydisperse and multicomponent aerosol that activate INPs following different freezing parameterizations. The AC-1D model is informed by large-eddy simulations to probe the impact of INP representation on predicted ice crystal number concentrations (N i ) and ice crystal budgets in mixed-phase Arctic stratus. We apply three immersion freezing (IMF) parameterizations, two time-independent (singular) and one time-dependent (classical nucleation theory), to predict the evolution of the INP reservoir and resulting ice crystal budget from polydisperse mineral dust, organic (humic-like substances), and sea spray aerosol particle size distributions. Our analysis focuses on how variations in aerosol number concentration and cloud system parameters such as cloud cooling rate, cloud-top entrainment rate, and ice crystal fall speed influence the INP reservoir and ice crystal budgets. Furthermore, this study investigates the competitive ice nucleation dynamics in mixed aerosol environments and provides a process-level quantification of the INP budget terms, which directly controls ice crystal budgets. For all studied case scenarios, the aerosol types and associated particle size distributions significantly impact INP and N i , and the choice between a time-dependent and a singular freezing description yields orders-of-magnitude differences in the predicted INP and N i over the 10 h simulation time, reflecting typical cloud lifetimes. Our results show that the influence of cloud cooling, INP entrainment, and sedimentation varies significantly depending on the chosen freezing parameterization. These findings underscore the critical need for robust IMF parameterizations and precise cloud system observations to enhance the accuracy of models in predicting mixed-phase cloud structure and evolution.

Arctic clouds↗
DOE OSTI2024

Study of Microphysical Signatures Based on Spectral Polarimetry during the RELAMPAGO Field Experiment in Argentina

Weather radars with dual-polarization capabilities enable the study of various characteristics of hydrometeors, including their size, shape, and orientation. Radar polarimetric measurements, coupled with Doppler information, allow for analysis in the spectral domain. This analysis can be leveraged to reveal valuable insight into the microphysics and kinematics of hydrometeors in precipitation systems. This paper uses spectral polarimetry to investigate precipitation microphysics and kinematics in storm environments observed during the Remote Sensing of Electrification, Lightning, and Mesoscale/Microscale Processes with Adaptive Ground Observations (RELAMPAGO) field experiment in Argentina. This study uses range–height indicator scan measurements from a C-band polarimetric Doppler weather radar deployed during the field campaign. Here, in this work, the impact of storm dynamics on hydrometeors is studied, including the size sorting of hydrometeors due to vertical wind shear. In addition, particle microphysical processes because of aggregation and growth of ice crystals in anvil clouds, as well as graupel formation resulting from the riming of ice crystals and dendrites, are also analyzed here. The presence of different particle size distributions because of the mixing of hydrometeors in a sheared environment and resulting size sorting has been reported using spectral differential reflectivity (sZ dr ) slope. Spectral reflectivity sZ h and sZ dr have also been used to understand the signature of ice crystal aggregation in an anvil cloud. The regions of pristine ice crystals are identified from vertical profiles of spectral polarimetric variables in anvil cloud because of sZ h < 0 dB and sZ dr values around 2 dB. It is also found that the growth process of these ice crystals causes a skewed bimodal sZ h spectrum due to the presence of both pristine ice crystals and dry snow. Next, graupel formation due to riming has been studied, and it is found that the riming process produces sZ h values of about 10 dB and corresponding sZ dr values of 1 dB. This positive sZ dr indicates the presence of needle/columnar secondary ice particles formed by ice multiplication processes in the riming zones. Last, the temporal evolution of a storm is investigated by analyzing changes in hydrometeor types with time and their influence on the spectral polarimetric variables.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Quantifying Uncertainty in Ice Particle Velocity–Dimension Relationships Using MC3E Observations

Abstract Ice particle terminal fall velocity ( V t ) is fundamental for determining microphysical processes, yet remains extremely challenging to measure. Current theoretical best estimates of V t are functions of Reynolds number. The Reynolds number is related to the Best number, which is a function of ice particle mass, area ratio ( A r ), and maximum dimension ( D max ). These estimates are not conducive for use in most models since model parameterizations often take the form , where ( α , β ) depend on habit and D max . A previously developed framework is used to determine surfaces of equally plausible ( α , β ) coefficients whereby ice particle size/shape distributions are combined with V t best estimates to determine mass- ( V M ) or reflectivity-weighted ( V Z ) velocities that closely match parameterized V M ,SD or V Z ,SD calculated using the ( α , β ) coefficients using two approaches. The first uses surfaces of equally plausible ( a , b ) coefficients describing mass ( M )–dimension relationships (i.e., ) to calculate mass- or reflectivity-weighted velocity from size/shape distributions that are then used to determine ( α , β ) coefficients. The second investigates how uncertainties in A r , D max , and size distribution N ( D ) affect V M or V Z . For seven of nine flight legs flown on 20 and 23 May 2011 during the Mesoscale Continental Convective Clouds Experiment (MC3E), uncertainty from natural parameter variability—namely, the variability in ice particle parameters in similar meteorological conditions—exceeds uncertainties arising from different A r assumptions or D max estimates. The combined uncertainty between A r , D max , and N ( D ) produced smaller variability in ( α , β ) compared to varying M ( D ), demonstrating M ( D ) must be accurately quantified for model fall velocities. Primary sources of uncertainty vary considerably depending on environmental conditions. Significance Statement Ice particle fall velocity is fundamental for numerous processes within clouds, and hence is a critical property that must be accurately represented in weather and climate models. Using aircraft observations of ice particle shapes and sizes obtained in clouds behind midlatitude thunderstorms, this work develops a new framework for estimating ice particle fall velocities and their uncertainty, including quantifying the importance of different uncertainty sources from cloud microphysics measurements. Natural parameter variability contributes the most uncertainty in ice particle fall velocity estimates, although other sources can also be important contributors to uncertainty in certain conditions. Additional work examining ice particle data is needed to further understand how dependent uncertainty in certain ice particle properties are to local environmental conditions.

Meteorology & Atmospheric Sciences↗
DOE OSTI2023

Vertical Structure and Ice Production Processes of Shallow Convective Postfrontal Clouds over the Southern Ocean in MARCUS. Part II: Modeling Study

Abstract Part I of this series presented a detailed overview of postfrontal mixed-phase clouds observed during the Measurements of Aerosols, Radiation, and Clouds over the Southern Ocean (MARCUS) field campaign. In Part II, we focus on a multiday (23–26 February 2018) case with the aim of understanding ice production as well as model sensitivity to ice process parameterizations using the Weather Research and Forecasting (WRF) Model. The control simulation with the Predicted Particle Properties (P3) microphysics scheme underestimates the ice content and overestimates the supercooled liquid water content, contrary to the bias common in global climate models. The simulations targeted at ice production processes show negligible sensitivity to cloud droplet number concentrations. Further, neither increasing ice nuclei particle (INP) concentrations to an unrealistic level nor adjusting it to MARCUS field estimations alone guarantees more ice production in the model. However, the simulated clouds are found to be highly sensitive to the implementation of immersion freezing, the thresholding of condensation/deposition freezing initiation, and the rime splintering process. By increasing immersion freezing of cloud droplets, relaxing thresholds for condensation/deposition freezing, or removing rime splintering thresholds, the model significantly improves its performance in producing ice. The relaxation of the immersion freezing temperature threshold to the observed cloud-top temperature suggests an in-cloud seeder–feeder mechanism. The results of this work call for an increase in observations of INP, especially over the remote Southern Ocean and at relatively high temperatures, and measurements of ice particle size distributions to better constrain ice nucleating processes in models.

Meteorology & Atmospheric Sciences↗
DOE OSTI2020

A New Method for Ice-Ice Aggregation in the Adaptive Habit Model

In this paper, a novel methodology for modeling ice-ice aggregation is presented. This methodology combines a modified hydrodynamic collection algorithm with bulk aggregate characteristic information from an offline simulator that collects ice particles, namely, the Ice Particle and Aggregate Simulator, and has been implemented into the Adaptive Habit Microphysics scheme in the Weather Research and Forecasting Model. Aggregates, or snow, are formed via collection of cloud ice particles, where initial ice characteristics and the resulting geometry determine aggregate characteristics. Upon implementation, idealized squall-line simulations are performed to examine the new methodology in comparison with commonly used bulk microphysics schemes. It is found that the adaptive habit aggregation parameterization develops snow and reduces ice mass and number concentrations compared to other schemes. The development of aggregates through the new methodology cascades into other interesting effects, including enhancements in ice and snow growth, as well as homogeneous freezing. Further microphysical analyses reveal varying sensitivities, where snow processes are most sensitive to the new parameterization, followed by ice, then cloud, rain, and graupel processes. Further, the new scheme results in enhancements in surface precipitation due to the persistence of snow at lower altitudes. This persistence is a result of shape-dependent melting and sublimation, increasing the residence time. Moreover, these low-level enhancements are reflected in increases in radar reflectivity at the surface and its spatial distribution. Finally, the ability to predict snow shape and density allows for the simulation of polarimetric radar quantities, resulting in signature enhancements compared to schemes that do not consider spatial and temporal variations in snow shape and density.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Enrichment of Phosphates, Lead, and Mixed Soil–Organic Particles in INPs at the Southern Great Plains Site

Ice nucleating particles (INPs) are rare particles that initiate primary ice formation, a critical step required for subsequent important cloud microphysical processes that ultimately govern cloud phase and cloud radiative properties. Laboratory studies have found that organic-rich dusts, such as those found in soils, are more efficient INPs compared to mineral dust. However, the atmospheric relevance of these organic-rich dusts are not well understood, particularly in regions with significant agricultural activity. The Agricultural Ice nuclei at the Southern Great Plains field campaign (AGINSGP) was conducted in rural Oklahoma to investigate how soil dusts contribute to INP populations in the Great Plains. We present chemical characterization of ambient and ice crystal residual particles from a single day of sampling, using single particle mass spectrometry (SPMS) and scanning microscopy. Ambient particles were primarily carbonaceous or secondary aerosol, while the fraction of dust particles was higher in the residual particles. We also observed an unusual particle type consisting of a carbonaceous core mixed with dust fragments on the surface, which was found in higher proportion in residuals. Dust particles measured during residual sampling contained greater proportions of phosphate $(^{63}\text{PO}^-_2$ and $^{79}\text{PO}^-_3)$ and lead ( 206 Pb + ). Strong sulfate signals were not seen in the residual dust particles measured by the SPMS, while nitrate was slightly depleted relative to ambient dust. This study shows that organic-rich soils may be important contributors to the ambient INP population in agricultural regions.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Dataset for "Enrichment of phosphates, lead, and mixed soil-organic particles in INPs at the Southern Great Plains site"

Ice nucleating particles (INPs) are rare particles that initiate primary ice formation, a critical step required for subsequent important cloud microphysical processes that ultimately govern cloud phase and cloud radiative properties. Laboratory studies have found that organic-rich dusts, such as those found in soils, are more efficient INPs compared to mineral dust. However, the atmospheric relevance of these organic-rich dusts are not well understood, particularly in regions with significant agricultural activity. The Agricultural Ice nuclei at the Southern Great Plains field campaign (AGINSGP) was conducted in rural Oklahoma to investigate how soil dusts contribute to INP populations in the Great Plains. We present chemical characterization of ambient and ice crystal residual particles from a single day of sampling using single particle mass spectrometry (SPMS) and scanning microscopy. Ambient particle concentrations were primarily carbonaceous or secondary aerosol, while dust particles were enhanced in the residual particles. Dust particles measured during residual sampling contained elevated signals for phosphate ($^{63}$PO$_2^-$ and $^{79}$PO$_3^-$) and lead ($^{208}$Pb$_2^+$). Nitrate was slightly depleted relative to ambient dust, and strong sulfate signals were not seen in the residual particles measured by the miniSPLAT. The residual and ambient particles measured by the miniSPLAT were approximately the same size, though this may be a reflection of the SPMS size-dependent transmission efficiency and not that of physical truth. This study shows that organic-rich soils may be important contributors to the ambient INP population in agricultural regions.

Cornwell, Gavin C↗