Ice in Convective Storm Environments: Assessing Impacts of Microphysics on Simulated Reflectivities

Julia A. Shates; Derek J Posselt; Matthew D Lebsock; Joseph A Finlon; Greg M McFarquhar; Andrew J Heymsfield

Home/Archive/NASA NTRS

NASA NTRS · 20260004970

Ice in Convective Storm Environments: Assessing Impacts of Microphysics on Simulated Reflectivities

Julia A. Shates·Derek J Posselt·Matthew D Lebsock·Joseph A Finlon·Greg M McFarquhar·Andrew J Heymsfield

View original source ↗Read source document ↗Export metadata ↓

Abstract

Satellite radar retrievals of clouds and precipitation rely on assumptions about the microphysical properties of hydrometeors such as particle size distributions (PSDs) and massdimensional relationships, which have been shown to vary in distinct environments. A deeper understanding of linkages between mass-dimensional relationships, ice crystal shape, and radar reflectivities across cloud and precipitation regimes is necessary for developing new retrievals for upcoming satellite radar missions, including the INvestigation of Convective UpdraftS (INCUS), as well as future missions that will feature space-borne radars observing clouds, convection, and precipitation. This study investigates ice microphysical characteristics through aircraft in-situ and remote sensing observations from two convective events during the Mid-latitude Continental Convective Clouds Experiment (MC3E) in the Southern Great Plains, and tests implications for ice particle habit in forward modeled radar observations. During MC3E, the Ka/Ku- band High Altitude Imaging Wind and Rain Airborne Profiler (HIWRAP) was onboard the NASA ER-2 aircraft, while in-situ measurements were taken onboard the University of North Dakota Citation. In-situ PSDs and particle shape information from optical array probe data are used to forward model radar reflectivities from particle scattering databases and are compared against HIWRAP reflectivities. Various particle combinations can be used to find agreement with HIWRAP observations, and the forward modeled reflectivities were sensitive to particle type. Varying the riming on the aggregate and graupel particles resulted in larger than 10 dB differences in reflectivities. The in-situ and radar observations imply the occurrence of aggregation, size sorting and lofting particles.

Keep this discovery

BibTeX RIS

Julia A. Shates, Derek J Posselt, Matthew D Lebsock, Joseph A Finlon, Greg M McFarquhar, Andrew J Heymsfield. 2026-05-25. Ice in Convective Storm Environments: Assessing Impacts of Microphysics on Simulated Reflectivities. https://ntrs.nasa.gov/citations/20260004970

Cite the original work for its findings. Save a collection to share your selection of sources.

KEEP EXPLORING

Related reports

DOE OSTI2026

Aerosol-deep convection interaction based on joint cell-thermal tracking in Large Eddy Simulations during the TRACER campaign

In cumulus clouds, aerosol concentrations control cloud droplet concentrations, modifying cloud radiative properties, precipitation processes, and cloud electrification. However, mechanisms of aerosol-deep convection interactions are not well understood due to complex cloud dynamics and microphysics. We investigate the interaction of aerosols with isolated deep convection using Large Eddy Simulations of two cases during the TRacking Aerosol Convection interactions ExpeRiment (TRACER) near Houston, Texas, using a joint cell-thermal tracking algorithm. Cumulus thermals are droplet generators, since supersaturation and droplet nucleation coincide with thermal centers, where the strongest updrafts occur. Primary ice crystal formation does not take place inside thermals, but at layers where previous thermals detrained moisture. As subsequent thermals containing supercooled droplets penetrate these layers, hail and graupel form at or near these thermals. Higher aerosol concentrations result in higher droplet concentrations that suppress drizzle, delay warm rain processes, and transport more moisture aloft. This increases snow and ice amount, as well as graupel and hail, leading to more lightning. Polluted thermals initiate at slightly higher altitudes, and are slightly larger and faster, suggesting a weak invigoration. We also find more thermals per cell, but fewer isolated cells, since convection is more aggregated and intense, especially near the end of the 24 h simulation. Non-linear mesoscale feedback likely triggered by temperature and moisture responses to aerosol-thermal interactions causes the aggregation. Time-lagged aerosol-reinitialization experiments show that the mesoscale response is the predominant forcing for the invigoration. These changes happen within one day, on a smaller scale than previously suggested.

54 ENVIRONMENTAL SCIENCES↗

NASA NTRS2020

An Overview of CMIP5 and CMIP6 Simulated Cloud Ice, Radiation Fields, Surface Wind Stress, Sea Surface Temperatures and Precipitation over Tropical and Subtropical Oceans

The potential links between ice water path (IWP), radiation, circulation, sea surface temperature (SST) and precipitation over the Pacific and Atlantic Oceans resulting from the falling ice radiative effects (FIREs) are examined from present day model outputs of CMIP5 and CMIP6. The latter is divided into two subsets with (SON6) and without FIREs (NOS6) as more models with FIREs are included in CMIP6 than in CMIP5. Improvement in floating cloud ice (~20 g m-2) is noticeable over convective regions in CMIP6 relative to CMIP5. The inclusion of FIREs in SON6 subset may contribute to reduce biases of overestimated outgoing longwave radiation and downward surface shortwave and overestimated reflected shortwave at the top of the atmosphere (TOA) by magnitudes of 4?8 W m-2 over convective regions against CERES, compared to NOS6 subset. The reduced biases in radiative fluxes in convective regions stabilize the atmosphere and lead to circulation, SST, cloud and precipitation changes over the trade-wind regions, as seen from improved radiative fluxes (4?15 W m-2), surface wind stress biases, SST (0.2?0.8 K) and precipitation (1 mm day-1) biases. The significant improvement from NOS6 to SON6 leads to improved multi-model means for CMIP6 relative to CMIP5 for radiation fields over the trade wind regions but the degradation over convective zones is attributed to NOS6 subset. The results suggest that other sources of uncertainty and deficiencies in climate models may play significant roles for reducing discrepancies although FIREs, via radiation-circulation coupling, may be one of the factors that help to reduce regional biases.

Jui-Lin F Li↗

DOE OSTI2026

Evolution of the Antarctic Ice Sheet from 2000–2300 and beyond: model sensitivity and uncertainty analysis using MPAS-Albany Land Ice

We present a description of the Antarctic Ice Sheet model configuration submitted to the ISMIP6-Antarctica-2300 experiment using the MPAS-Albany Land Ice model, along with three new sets of simulations: (1) a set of extended simulations to 2500 for three forced experiments and to 2775 for the control experiment; (2) a sensitivity analysis of our model configuration to parameters controlling basal sliding and sub-shelf melt, and to model structural choices including the choice of the energy and stress balances; and (3) a 72-member ensemble run on graphics processing units (GPUs) and analysis of variance to determine the primary sources of uncertainty in our ice-sheet model projections. Our extended simulations predict rapid retreat beginning after 2300 for SSP1-2.6 forcing and after 2500 for present-day (control) forcing, primarily in the Amundsen Sea Embayment. We find that varying the sub-shelf melt parameter between the 5th to 95th percentile values for a mean-Antarctic calibration target results in an up to ∼ ± 40 % change in sea-level contribution relative to our baseline simulations that used the median value. Using a linear basal sliding law reduces sea-level contribution by 51 %–73 % relative to our baseline nonlinear sliding law with an exponent of 1/5. When using basal sliding law exponents of 1/3 and 1/10, the overall difference from our baseline simulations at 2300 is on the order of 10 %. The Amundsen Sea Embayment region displays a strongly non-linear dependence of mass loss on the sliding law exponent, with no discernible relationship between the sliding law exponent and the mass loss by 2300, while the sectors feeding the Ross and Filchner-Ronne ice shelves exhibit more mass loss with a more-plastic sliding law. Our model fidelity sensitivity experiments reveal a 9 %–31 % increase in sea-level contribution when using a depth-integrated stress balance approximation relative to our three-dimensional solver, while using a fixed-in-time temperature field increases sea-level contribution by 14 %–88 % relative to two thermomechanically coupled configurations. Our 72-member ensemble and analysis of variance show that the uncertainty in long-term projections is dominated by the choice of Earth system model forcing and the presence or absence of hydrofracture forcing, rather than uncertainty in sliding and sub-shelf melt parameters.

58 GEOSCIENCES↗