Fragmentation in Collisions of Snow with Graupel/Hail: New Formulation from Field Observations
Explore the source record for details and available documents.
Engineering topics
Publications and source records attributed to Bansemer, Aaron.
Explore the source record for details and available documents.
The NASA Investigation of Microphysics and Precipitation for Atlantic Coast Threatening Snowstorms (IMPACTS) field campaign provides high-quality, high-altitude aircraft lidar (532 nm), radar (W-band) and in-cloud microphysical aircraft data taken during wintertime storm events impacting the United States. This study evaluates two mass-dimensional relationships (Brown and Francis (1995, BF95); Heymsfield (2014, H14) and two lidar-radar microphysical retrieval algorithms (Cloudsat and CALIPSO Ice Cloud Property Product (2C-ICE); VarPy (a variational method derived from the satellite lidar-radar data community)) to estimate aircraft-retrieved volume extinction coefficient (σ), ice water content (IWC), and effective radius (r e ) during the 2020 IMPACTS deployment. BF95 and H14 have a close 1:1 correlation (R 2 = 0.98) with in-situ observations of σ. However, only BF95 displays a linear, consistent, and almost temperature-independent low bias for IWC and r e , which likely arises from the environmental conditions used to determine each. Unlike the field-campaign-derived BF95 and H14 relationships, VarPy and 2C-ICE directly ingest the aircraft-based lidar and radar data to simulate σ, IWC, and r e . For all three microphysical parameters, VarPy and 2C-ICE retrieval errors became notably more pronounced around the dendritic growth zone (-15°C to -10°C) and near freezing (≥-5°C), which suggests that both algorithms experience difficulty addressing riming and aggregation processes and with larger particles (dendrites and plates) due in part to their simplified ice particle assumptions. However, the mean-melt diameter ice-particle assumption did yield more accurate IWC estimates, which led to slightly better overall results for VarPy.
This 3-year DoE-funded joint project had the over-arching aim of understanding how ice is initiated in clouds of various types. Focus was given to processes of fragmentation of pre-existing ice, which can occur in positive feedback loops (‘ice multiplication’). A basic question to address was which fragmentation processes prevail in which basic cloud-types. The approach was to use cloud models and field observations, while pioneering our own lab observations of ice initiation to break the deadlock from the past lack of lab observations. Historically, the tendency of the cloud physics community to avoid doing lab observations has allowed a vast gap in knowledge about ice initiation to persist for decades. During the first part of the project, new formulations were created to treat two overlooked types of fragmentation of ice. First, sublimational breakup of ice was treated based on a theoretical formula that we fitted to a pooled dataset of lab observations published previously in the literature. Second, a new mode of fragmentation of freezing raindrops was treated, which involves a supercooled drop being hit by a more massive ice particle. Some of the secondary droplets from the impact freeze. This work was done at Manchester University by Co-I Connolly. Then during the second part, both formulations were implemented in our ‘aerosol-cloud model’ (AC). AC has a hybrid bin/bulk microphysics scheme, and now represents four processes of SIP. The accuracy of AC was evaluated for four cases typifying four basic cloud-types: slightly cold-based stratiform cloud and cold-, warm- and very warm-based convective clouds. We discovered that the warmth of cloud-base, especially in the tropics, promotes SIP processes of raindrop-freezing fragmentation and rime-splintering, and surprisingly, sublimational breakup too. It was found that breakup in ice-ice collisions is ubiquitous. Finally, a portable laboratory chamber was constructed at Lund and deployed in northern Sweden to observe breakup in graupel-snow collisions outdoors. This was seen to be even more prolific than treated in our 2018 formulation. Papers describing results are either published or soon to be published.
NASA and the FAA conducted two flight campaigns to quantify onboard weather radar measurements with in-situ measurements of high concentrations of ice crystals found in deep convective storms. The ultimate goal of this research was to improve the understanding and develop onboard weather radar processing to detect regions of high ice water content ahead of an aircraft and enable tactical avoidance of the potentially hazardous conditions. Both High Ice Water Content (HIWC) RADAR campaigns utilized the NASA DC-8 Airborne Science Laboratory which was equipped with a Honeywell RDR-4000 weather radar and icing instruments to characterize the ice crystal clouds. The purpose of this paper is to summarize how these campaigns were conducted and highlight key results.
NASA and the FAA (Federal Aviation Administration) conducted two flight campaigns to quantify onboard weather radar measurements with in-situ measurements of high concentrations of ice crystals found in deep convective storms. The ultimate goal of this research was to improve the understanding and develop onboard weather radar processing to detect regions of high ice water content ahead of an aircraft and enable tactical avoidance of the potentially hazardous conditions. Both High Ice Water Content (HIWC) RADAR campaigns utilized the NASA DC-8 Airborne Science Laboratory which was equipped with a Honeywell RDR-4000 weather radar and icing instruments to characterize the ice crystal clouds. The purpose of this paper is to summarize how these campaigns were conducted and highlight key results.
NASA and the FAA conducted two flight campaigns to quantify onboard weather radar measurements with in-situ measurements of high concentrations of ice crystals found in deep convective storms. The ultimate goal of this research was to improve the understanding and develop onboard weather radar processing to detect regions of high ice water content ahead of an aircraft and enable tactical avoidance of the potentially hazardous conditions. Both High Ice Water Content (HIWC) RADAR campaigns utilized the NASA DC-8 Airborne Science Laboratory which was equipped with a Honeywell RDR-4000 weather radar and icing instruments to characterize the ice crystals clouds. The purpose of this paper is to summarize how these campaigns were conducted and highlight key results.
Occurrence frequency and dynamical conditions of ice supersaturation (ISS, where relative humidity with respect to ice (RHi) greater than 100%) are examined in the upper troposphere around convective activity. Comparisons are conducted between in situ airborne observations and the Weather Research and Forecasting model simulations using four double-moment microphysical schemes at temperatures less than or or equal to -40degdegC. All four schemes capture both clear-sky and in-cloud ISS conditions. However, the clear-sky (in-cloud) ISS conditions are completely (significantly) limited to the RHi thresholds of the Cooper parameterization. In all of the simulations, ISS occurrence frequencies are higher by approximately 3-4 orders of magnitude at higher updraft speeds (greater than 1 m s(exp -1) than those at the lower updraft speeds when ice water content (IWC) greater than 0.01 gm(exp -3), while observations show smaller differences up to approximately 1-2 orders of magnitude. The simulated ISS also occurs less frequently at weaker updrafts and downdrafts than observed. These results indicate that the simulations have a greater dependence on stronger updrafts to maintain/generate ISS at higher IWC. At lower IWC (less than or equal or 0.01 gm(exp -3), simulations unexpectedly show lower ISS frequencies at stronger updrafts. Overall, the Thompson aerosol-aware scheme has the closest magnitudes and frequencies of ISS greater than 20% to the observations, and the modified Morrison has the closest correlations between ISS frequencies and vertical velocity at higher IWC and number density. The Cooper parameterization often generates excessive ice crystals and therefore suppresses the frequency and magnitude of ISS, indicating that it should be initiated at higher ISS (e.g.,lees than or equal to 25%).
We use aircraft observations combined with the reanalysis data to investigate the radiative effects of ice supersaturation (ISS). Our results show that although the excess water vapor over ice saturation itself has relatively small radiative effects, mistaking it as ice crystals in climate models would lead to considerable impacts: on average, +2.49 W/m(exp 2) change in the top of the atmosphere (TOA) radiation, -2.7 W/m(exp 2) change in surface radiation, and 1.47 K/d change in heating rates. The radiative effects of ISS generally increase with the magnitudes of supersaturation. However, there is a strong dependence on the preexisting ice water path, which can even change the sign of the TOA radiative effect. It is therefore important to consider coexistence between ISS and ice clouds and to validate their relationship in the parameterizations of ISS in climate models.
In this study, two different particle models describing the structure and electromagnetic properties of snow are developed and evaluated for potential use in satellite combined radar-radiometer precipitation estimation algorithms. In the first model, snow particles are assumed to be homogeneous ice-air spheres with single-scattering properties derived from Mie theory. In the second model, snow particles are created by simulating the self-collection of pristine ice crystals into aggregate particles of different sizes, using different numbers and habits of the collected component crystals. Single-scattering properties of the resulting nonspherical snow particles are determined using the discrete dipole approximation. The size-distribution-integrated scattering properties of the spherical and nonspherical snow particles are incorporated into a dual-wavelength radar profiling algorithm that is applied to 14- and 34-GHz observations of stratiform precipitation from the ER-2 aircraft-borne High-Altitude Imaging Wind and Rain Airborne Profiler (HIWRAP) radar. The retrieved ice precipitation profiles are then input to a forward radiative transfer calculation in an attempt to simulate coincident radiance observations from the Conical Scanning Millimeter-Wave Imaging Radiometer (CoSMIR). Much greater consistency between the simulated and observed CoSMIR radiances is obtained using estimated profiles that are based upon the nonspherical crystal/aggregate snow particle model. Despite this greater consistency, there remain some discrepancies between the higher moments of the HIWRAP-retrieved precipitation size distributions and in situ distributions derived from microphysics probe observations obtained from Citation aircraft underflights of the ER-2. These discrepancies can only be eliminated if a subset of lower-density crystal/aggregate snow particles is assumed in the radar algorithm and in the interpretation of the in situ data.
This study explores relationships between lightning, cloud microphysics, and tropical cyclone (TC) storm structure in Hurricane Karl (16 September 2010) using data collected by the NASA DC-8 and Global Hawk (GH) aircraft during NASA's Genesis and Rapid Intensification Processes (GRIP) experiment. The research capitalizes on the unique opportunity provided by GRIP to synthesize multiple datasets from two aircraft and analyze the microphysical and kinematic properties of an electrified TC. Five coordinated flight legs through Karl by the DC-8 and GH are investigated, focusing on the inner-core region (within 50km of the storm center) where the lightning was concentrated and the aircraft were well coordinated. GRIP datasets are used to compare properties of electrified and nonelectrified inner-core regions that are related to the noninductive charging mechanism, which is widely accepted to explain the observed electric fields within thunderstorms. Three common characteristics of Karl's electrified regions are identified: 1) strong updrafts of 10-20ms21, 2) deep mixed-phase layers indicated by reflectivities.30 dBZ extending several kilometers above the freezing level, and 3) microphysical environments consisting of graupel, very small ice particles, and the inferred presence of supercooled water. These characteristics describe an environment favorable for in situ noninductive charging and, hence, TC electrification. The electrified regions in Karl's inner core are attributable to a microphysical environment that was conducive to electrification because of occasional, strong convective updrafts in the eyewall.
One of the main goals of the Tropical Composition, Cloud and Climate Coupling Experiment (TC(sup 4)) during July and August 2007 was to gain a better understanding of the formation and life cycle of cirrus clouds in the upper troposphere and lower stratosphere and how their presence affects the exchange of water vapor between these layers. Additionally, it is important to compare in situ measurements taken by aircraft instruments with products derived from satellite observations and find a meaningful way to interpret the results. In this study, cloud properties derived using radiance measurements from the Geostationary Operational Environmental Satellite (GOES) imagers are compared to similar quantities from aircraft in situ observations and are examined for meaningful relationships. A new method using dual \angle satellite measurements is used to derive the ice water content (IWC) for the top portion of deep convective clouds and anvils. The results show the in situ and remotely sensed mean microphysical properties agree to within approx.10 microns in the top few kilometers of thick anvils despite the vastly different temporal and spatial resolutions of the aircraft and satellite instruments. Mean particle size and IWC are shown to increase with decreasing altitude in the top few kilometers of the cloud. Given these relationships, it may be possible to derive parameterizations for effective particle size and IWC as a function of altitude from satellite observations
An analysis of two days of in situ observations of ice particle size spectra, in convectively generated cirrus, obtained during NASA s Tropical Composition, Cloud, and Climate Coupling (TC4) mission is presented. The observed spectra are examined for their fit to the exponential, gamma, and lognormal function distributions. Characteristic particle size and concentration density scales are determined using two (for the exponential) or three (for the gamma and lognormal functions) moments of the spectra. It is shown that transformed exponential, gamma, and lognormal distributions should collapse onto standard curves. An examination of the transformed spectra, and of deviations of the transformed spectra from the standard curves, shows that the lognormal function provides a better fit to the observed spectra.
The effective radius(r(sub e)) is a crucial variable in representing the radiative properties of cloud layers in general circulation models. This parameter is proportional to the condensed water content (CWC) divided by the extinction (sigma). For ice cloud layers, parameterizations for r(sub e), have been developed from aircraft in-situ measurements 1) indirectly, using data obtained from particle spectrometer probes and assumptions or observations about particle shape and mass to get the ice water content (IWC) and area to get sigma, and recently 2) from probes that measure IWC and sigma directly. This study compares [IWC/sigma] derived from the two methods using data sets acquired from comparable instruments on two aircraft, one sampling clouds at mid-levels and the other at upper-levels during the CRYSTAL-FACE field program in Florida in 2002. The sigma and IWC derived by each method are compared and evaluated in different ways for each aircraft data set. Direct measurements of sigma exceed those derived indirectly by a factor of two to two and a half. The IWC probes, relying on ice sublimation, appear to measure accurately except when the IWC is high or the particles too large to sublimate completely during the short transit time through the probe. The IWC estimated from the particle probes are accurate when direct measurements are available to provide constraints and useful information in high IWC/large particle situations. Because of the discrepancy in sigma estimates between the direct and indirect approaches, there is a factor of 2 to 3 difference in [IWC/sigma] between them. Although there are significant uncertainties involved in its use, comparisons with several independent data sources suggest that the indirect method is the more accurate of the two approaches. However, experiments are needed to resolve the source of the discrepancy in sigma.
Simulations of a cumulonimbus cloud observed in the Cirrus regional Study of Tropical Anvils and Cirrus Layers-Florida Area Cirrus Experiment (CRYSTAL-FACE) with an advanced version of the Explicit Microphysics Model (EMM) are presented. The EMM has size-resolved aerosols and predicts the time evolution of sizes, bulk densities and axial ratios of ice particles. Observations by multiple aircraft in the troposphere provide inputs to the model, including observations of the ice nuclei and of the entire size distribution of condensation nuclei. Homogeneous droplet freezing is found to be the source of almost all of the ice crystals in the anvil updraught of this particular model cloud. Most of the simulated droplets that freeze to form anvil crystals appear to be nucleated by activation of aerosols far above cloud base in the interior of the cloud ("secondary" or "in cloud" droplet nucleation). This is partly because primary droplets formed at cloud base are invariably depleted by accretion before they can reach the anvil base in the updraught, which promotes an increase with height of the average supersaturation in the updraught aloft. More than half of these aerosols, activated far above cloud base, are entrained into the updraught of this model cloud from the lateral environment above about 5 km above mean sea level. This confirms the importance of remote sources of atmospheric aerosol for anvil glaciation. Other nucleation processes impinge indirectly upon the anvil glaciation by modifying the concentration of supercooled droplets in the upper levels of the mixed-phase region. For instance, the warm-rain process produces a massive indirect impact on the anvil crystal concentration, because it determines the mass of precipitation forming in the updraught. It competes with homogeneous freezing as a sink for cloud droplets. The effects from turbulent enhancement of the warm-rain process and from the nucleation processes on the anvil ice properties are assessed.
A customized product for analysis of microphysics data collected from aircraft during field campaigns in support of the TRMM program is described. These Common Microphysics Products (CMP's) are designed to aid in evaluation of TRMM spaceborne precipitation retrieval algorithms. Information needed for this purpose (e.g., particle size spectra and habit, liquid and ice water content) was derived using a common processing strategy on the wide variety of microphysical instruments and raw native data formats employed in the field campaigns. The CMP's are organized into an ASCII structure to allow easy access to the data for those less familiar with and without the tools to accomplish microphysical data processing. Detailed examples of the CMP show its potential and some of its limitations. This approach may be a first step toward developing a generalized microphysics format and an associated community-oriented, non-proprietary software package for microphysics data processing, initiatives that would likely broaden community access to and use of microphysics datasets.
In this study, we report on the evolution of particle size distributions (PSDs) and habits as measured during slow, Lagrangian-type spiral descents through deep subtropical and tropical cloud layers in Florida, Brazil, and Kwajalein, Marshall Islands, most of which were precipitating. The objective of the flight patterns was to learn more about how the PSDs evolved in the vertical and to obtain information of the vertical structure of microphysical properties. New instrumentation yielding better information on the concentrations of particles in the size (D) range between 0.2 and 2 cm, as well as improved particle imagery, produced more comprehensive observations for tropical stratiform precipitation regions and anvils than have been available previously. Collocated radar observations provided additional information on the vertical structure of the cloud layers sampled. Most of the spirals began at cloud top, with temperatures (T) as low as -50 C, and ended at cloud base or below the melting layer (ML). The PSDs broadened from cloud top towards cloud base, with the largest particles increasing in size from several millimeters at cloud top to one centimeter or larger towards cloud base. Some continued growth was noted in the upper part of the ML. Concentrations of particles less than 1 mm in size decreased with decreasing height. The result was a consistent change in the PSDs in the vertical. Similarly, systematic changes in the size dependence of the particle cross-sectional area was noted with decreasing height. Aggregation-as ascertained from both the changes in the PSDs and evolution of particle habits as observed in high detail with the cloud particle imager (CPI) probe-was responsible for these trends. The PSDs were generally well-represented by gamma distributions of the form N = N0 gamma D microns e- lambda gamma D that were fitted to the PSDs over 1-km horizontal intervals throughout the spirals. The intercept (N0 gamma), slope (lambda gamma), and dispersion (microns) values were derived for each PSD. Exponential curves (N = N0e- lambdaD; micron = 0) were also fitted to the distributions. The lambda gamma values for given spirals varied systematically with temperature as did the values of lambda (exponential), and the data generally conformed to values found in previous studies involving exponential fits to size distributions in mid-latitude frontal and cirrus layers. Considerable variability often noted in the PSD properties during the loops of individual spirals was manifested primarily in large changes in N0 gamma and N0, but micron, lambda gamma and lambda remained fairly stable. Temperature is not found to be the sole factor controlling lambda gamma or lambda but is a primary one. Direct relationships were found between lambda gamma and N0 gamma or lambda gamma and micron for the gamma distributions and lambda and N0 for the exponential. The latter relationship was not found as distinctly in earlier studies; observed PSDs in this study had better fidelity with less scatter. The micron values changed monotonically with T over the range of temperatures and were directly related to N0 gamma or lambda gamma, thereby reducing the number of variables in the PSD functional equation to two. In the upper part of the ML, N0 and lambda continued to decrease, and in the lower part these values began to increase as the largest particles melted. We developed general expressions relating various bulk microphysical, radar, and radiative transfer-related variables to N0 gamma and lambda gamma, useful for both tropical and mid-latitude clouds. These relationships facilitate the specification of a number of bulk properties in cloud and climate models. The results presented in this paper apply best to temperatures between 0 and -40 C, for which the measured radar reflectivities fall in the range of 0 to 25 dBZe.
This study derives effective ice particle densities from data collected from the NASA WB-57F aircraft near the tops of anvils during the Cirrus Regional Study of Tropical Anvils and Cirrus Layers (CRYSTAL) Florida Area Cirrus Experiment (FACE) in southern Florida in July 2002. The effective density, defined as the ice particle mass divided by the volume of an equivalent diameter liquid sphere, is obtained for particle populations and single sizes containing mixed particle habits using measurements of condensed water content and particle size distributions. The mean effective densities for populations decrease with increasing slopes of the gamma size distributions fitted to the size distributions. The population-mean densities range from near 0.91 g/cu m to 0.15 g/cu m. Effective densities for single sizes obey a power-law with an exponent of about -0.55, somewhat less steep than found from earlier studies. Our interpretations apply to samples where particle sizes are generally below 200-300 microns in maximum dimension because of probe limitations.