SEARCH · Search NASA
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.
An optical ice particle counter was proposed as a companion instrument to the GSFC laser nephelometer. By counting ice particles and total cloud particles (both ice and liquid water), these two instruments may be used to study the balance between ice and water in clouds.
The solar reflectance bi-spectral (SRBS) and infrared split-window (IRSpW) methods are two of the most popular techniques for passive ice cloud property retrievals from multispectral imagers. Ice clouds are usually assumed to be vertically homogeneous in global operational algorithms based on these methods, although significant vertical variations of ice particle size are typically observed in ice clouds. In this Study we investigate uncertainties in retrieved optical thickness, effective particle size, and ice water path introduced by a homogeneous cloud assumption in both the SRBS and IRSpW methods, and focus on whether the assumption can lead to significant discrepancies between the two methods. The study simulates the upwelling spectral radiance associated with vertically structured clouds and passes the results through representative SRBS and IRSpW retrieval algorithms. Cloud optical thickness is limited to values for which IRSpW retrievals are possible (optical thickness less than about 7). When the ice cloud is optically thin and yet has a significant ice particle size vertical variation, it is found that both methods tend to underestimate the effective radius and ice water path. The reason for the underestimation is the nonlinear dependence of ice particle scattering properties (extinction and single scattering albedo) on the effective radius. Because the nonlinearity effect is stronger in the IRSpW than the SRBS method, the IRSpW-based IWP tends to be smaller than the SRBS counterpart. When the ice cloud is moderately optically thick, the IRSpW method is relatively insensitive to cloud vertical structure and effective radius retrieval is weighted toward smaller ice particle size, while the weighting function makes the SRBS method more sensitive to the ice particle size in the upper portion of the cloud. As a result, when ice particle size increases monotonically toward cloud base, the two methods are in qualitative agreement; in the event that ice particle size decreases toward cloud base, the effective radius and ice water path retrievals based on the SRBS method are substantially larger than those from the IRSpW. The main findings of this Study Suggest that the homogenous cloud assumption can affect the SRBS and IRSpW methods to different extents and, consequently, can lead to significantly different retrievals. Therefore caution should be taken when comparing and combining the ice cloud property retrievals from these two methods.
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.
The ice particles are formed when liquid and/or humid gases vent to the space. These submicroscopic ice particles are potential contamination sources of the environments during Skylab operations. The critical size of ice particle and its nucleation rate based on the theory of homogeneous and heterogeneous nucleation by sublimation are analyzed. The equations which are pertinent for studying the growth and evaporation of the ice particles are formulated. The mechanisms affecting the lifetime of ice particle are discussed. The gas dynamic techniques for experimental study of ice particle formation are proposed.
Ice clouds influence the climate system by changing the radiation budget and large-scale circulation. Therefore, climate models need to have an accurate representation of ice clouds and their radiative effects. In this paper, new broadband parameterizations for ice cloud bulk scattering properties are developed for severely roughened ice particles. The parameterizations are based on a general habit mixture that includes nine habits (droxtals, hollow/solid columns, plates, solid/hollow bullet rosettes, aggregate of solid columns, and small/large aggregates of plates). The scattering properties for these individual habits incorporate recent advances in light-scattering computations. The influence of ice particle surface roughness on the ice cloud radiative effect is determined through simulations with the Fu-Liou and the GCM version of the Rapid Radiative Transfer Model (RRTMG) codes and the National Center for Atmospheric Research Community Atmosphere Model (CAM, version 5.1). The differences in shortwave (SW) and longwave (LW) radiative effect at both the top of the atmosphere and the surface are determined for smooth and severely roughened ice particles. While the influence of particle roughening on the single-scattering properties is negligible in the LW, the results indicate that ice crystal roughness can change the SW forcing locally by more than 10 W m(exp −2) over a range of effective diameters. The global-averaged SW cloud radiative effect due to ice particle surface roughness is estimated to be roughly 1-2 W m(exp −2). The CAM results indicate that ice particle roughening can result in a large regional SW radiative effect and a small but nonnegligible increase in the global LW cloud radiative effect.
There are large uncertainties in the single-scattering (i.e., morphologies) and microphysical (i.e., concentrations) properties of ice particles whose size are less than ~100 µm. Insufficient resolutions of the most advanced cloud probes (e.g., cloud particle imager) cannot resolve the micrometer-scale morphologies of small ice particles. Further, the shattering of large ice particles on probes’ inlets or tips causes uncertainties in the measurement of the concentrations of small ice particles. These uncertainties have large impacts on the single-scattering and microphysical properties of small ice particles that are utilized to quantify the bulk radiative properties of cirrus. In this study, the impacts of uncertainties in the morphologies and concentrations of small ice particles on the bulk radiative properties of tropical cirrus were calculated using measurements acquired during the Tropical Warm Pool-International Cloud Experiment. Five different models (i.e., budding Buckyball, Chebyshev particle, droxtal, Gaussian random sphere, and sphere) that represent the shapes of small ice particles were used to calculate the single-scattering properties. The bulk radiative properties, average phase-function ($\overline{P_{11}}$), and average asymmetry parameter ($\overline{g}$) were computed by combining the measured size/habit distributions and the calculated single-scattering properties of ice particles. The impacts of the selection of varying morphologies of small particles on the bulk radiative properties were quantified. For these calculations, the possible range of the concentrations of small ice particles which depend on the degree of shattered large particles were also used. The impacts of varying the single-scattering properties of small ice particles on the bulk radiative properties were the largest in the upper parts of cirrus (T < –60 °C), while they were the smallest in the lower parts of cirrus (–45 < T < –30 °C). The impacts of uncertainties in the concentrations of small ice particles on the bulk radiative properties were largest in the lower parts of cirrus (–45 < T < –30 °C), whereas they were smallest in the upper parts of cirrus (T < –60 °C). The effect of shattering was maximum in the lower parts of cirrus, whilst it was minimum in the upper parts of cirrus. The combined impacts of uncertainties in the single-scattering (i.e., morphologies) and microphysical (i.e., concentrations) properties of small ice particles revealed variations of up to 11.2% (127.1%; 67.3%) of the integrated intensity in the forward (sideward; backward) angles in $\overline{P_{11}}$ and a corresponding change in $\overline{g}$ by up to 12.61%.
Determining the total amount of water contained in an icing cloud necessitates the measurement of both the liquid droplets and ice particles. One commonly accepted method for measuring cloud water content utilizes a hot wire sensing element, which is maintained at a constant temperature. In this approach, the cloud water content is equated with the power required to keep the sense element at a constant temperature. This method inherently assumes that impinging cloud particles remain on the sensing element surface long enough to be evaporated. In the case of ice particles, this assumption requires that the particles do not bounce off the surface after impact. Recent tests aimed at characterizing ice particle impact on a thermally heated wing section, have raised questions about the validity of this assumption. Ice particles were observed to bounce off the heated wing section a very high percentage of the time. This result could have implications for Total Water Content sensors which are designed to capture ice particles, and thus do not account for bouncing or breakup of ice particles. Based on these results, a test was conducted to investigate ice particle impact on the sensing elements of the following hot-wire cloud water content probes: (1) Nevzorov Total Water Content (TWC)/Liquid Water Content (LWC) probe, (2) Science Engineering Associates TWC probe, and (3) Particle Measuring Systems King probe. Close-up video imaging was used to study ice particle impact on the sensing element of each probe. The measured water content from each probe was also determined for each cloud condition. This paper will present results from this investigation and attempt to evaluate the significance of ice particle impact on hot-wire cloud water content measurements.
This dataset includes a unique set of cryo-SEM images of captured cirrus ice particle from ARM Project ICEBall, which was conducted at and around the Southern Great Plains site from October 21- November 1, 2021. The instrument used for particle collection, Ice Crystal Encapsulation by Ball (ICE-Ball), is described in Magee et al. 2021, 'Captured cirrus ice particles in high definition", Atmos. Chem. Phys. The ice particles were captured on balloon flights directly from cirrus cloud formations between 7-12 km on seven ICEBall flights on these dates: October 21, 23, 24, 25, 26, 31, and Nov. 1.
This work presents the results of an experimental study of ice particle impacts on a flat glass plate. The experiment was conducted at the Ballistics Impact Laboratory of NASA Glenn Research Center. The main objective of the experiment was to gain understanding about the modifications needed to the experimental configuration for a future parametric study at a larger range of values for particle diameters and other parameters. This was achieved by studying the effect of the velocity of an impacting ice particle on the post-impact fragment size and distribution for a reduced range of impacting particle diameters. Pre-impact particle diameter and velocity data were captured with a high-speed side camera. Post-impact fragment data were captured in a single frame with a 29-megapixel camera located above and normal to the target. Repeat runs were conducted for ice particles with diameters ranging from 1.7 to 2.9 millimeters, impacting at velocities between 39 and 98 meters per second. The fragment areas were measured, and the corresponding equivalent diameters and histogram distributions were calculated. Analysis of the data showed that the average equivalent diameter for the fragments in a run was an order of magnitude smaller than the diameter of the impacting ice particle. The histograms for equivalent diameter distribution were nonnormal with long tails, with most of the fragments having equivalent diameters concentrated toward the minimum value of the fragment size that could be resolved. Factors affecting the accuracy of the data during the digital imaging analysis were identified. Needed modifications to the setup to handle small size ice particles and other testing conditions were also identified.
Abstract The Ice Particle and Aggregate Simulator (IPAS) is used to theoretically represent the aggregation process of ice crystals. Aggregates have a variety of formations based on initial ice particle size, shape, and falling orientation, all of which influence water phase partitioning. Aggregate dimensional properties and density changes are calculated for monomer–monomer (MON–MON), monomer–aggregate (MON–AGG), and aggregate–aggregate (AGG–AGG) collection to be used by ice-microphysical models for improvement in aggregation parameterizations. Aggregates are chosen from a database of 9 744 000 preformed combinations to be further collected (see Part II). AGG–AGG collection results in more extreme and a smaller range of aggregate aspect ratios than MON–AGG collection. A majority of aggregates are closer to prolate than oblate spheroids, regardless of collection type, except for quasi-horizontally oriented particles that have extreme aspect ratios to begin with. MON–AGG collection frequently results in an increase in density upon collection, whereas MON–MON and AGG–AGG collection almost always result in particle density decreases, with extreme reductions near 99% for MON–MON collection. MON–MON collection results in the greatest decreases in density but then quickly becomes unaffected by the addition of more monomers due to inherent size differences between monomers and aggregates. Finally, a holistic analysis to in situ observations of cloud particle images is presented. IPAS 2D aspect ratios surround a median value of 0.6 and closely follow that of previous studies while varying by no more than ≈12% on average from observed aggregates.
The 1 November 1986 FIRE I case study was used to test an ice particle growth model which predicts bimodal size spectra in cirrus clouds. The model was developed from an analytically based model which predicts the height evolution of monomodal ice particle size spectra from the measured ice water content (IWC). Size spectra from the monomodal model are represented by a gamma distribution, N(D) = N(sub o)D(exp nu)exp(-lambda D), where D = ice particle maximum dimension. The slope parameter, lambda, and the parameter N(sub o) are predicted from the IWC through the growth processes of vapor diffusion and aggregation. The model formulation is analytical, computationally efficient, and well suited for incorporation into larger models. The monomodal model has been validated against two other cirrus cloud case studies. From the monomodal size spectra, the size distributions which determine concentrations of ice particles less than about 150 mu m are predicted.
This work presents the results of an experimental study of ice particle impacts on a moving wedge. The experiment was conducted in the Adverse Environment Rotor Test Stand (AERTS) facility located at Penn State University. The wedge was placed at the tip of a rotating blade. Ice particles shot from a pressure gun intercepted the moving wedge and impacted it at a location along its circular path. The upward velocity of the ice particles varied from 7 to 12 meters per second. Wedge velocities were varied from 0 to 120 meters per second. Wedge angles tested were 0, 30, 45, and 60. High speed imaging combined with backlighting captured the impact allowing observation of the effect of velocity and wedge angle on the impact and the post-impact fragment behavior. It was found that the pressure gun and the rotating wedge could be synchronized to consistently obtain ice particle impacts on the target wedge. It was observed that the number of fragments increase with the normal component of the impact velocity. Particle fragments ejected immediately after impact showed velocities higher than the impact velocity. The results followed the major qualitative features observed by other researchers for hailstone impacts, even though the reduced scale size of the particles used in the present experiment as compared to hailstones was 4:1.
Abstract. We provide the first dedicated laboratory study of collisions of supercooled water drops with ice particles as a secondary ice production mechanism. We experimentally investigated collisions of supercooled water drops (∼ 5 mm in diameter) with ice particles of a similar size (∼ 6 mm in diameter) placed on a glass slide at temperatures >-12 ∘C. Our results showed that secondary drops were generated during both the spreading and retraction phase of the supercooled water drop impact. The secondary drops generated during the spreading phase were emitted too fast to quantify. However, quantification of the secondary drops generated during the retraction phase with diameters >0.1 mm showed that 5–10 secondary drops formed per collision, with approximately 30 % of the secondary drops freezing over a temperature range between −4 and −12 ∘C. Our results suggest that this secondary ice production mechanism may be significant for ice formation in atmospheric clouds containing large supercooled drops and ice particles.
This study builds upon the microphysical modeling described in Part 1 by deriving formulations for the extinction and absorption coefficients in terms of the size distribution parameters predicted from the micro-physical model. The optical depth and single scatter albedo of a cirrus cloud can then be determined, which, along with the asymmetry parameter, are the input parameters needed by cloud radiation models. Through the use of anomalous diffraction theory, analytical expressions were developed describing the absorption and extinction coefficients and the single scatter albedo as functions of size distribution parameters, ice crystal shapes (or habits), wavelength, and refractive index. The extinction coefficient was formulated in terms of the projected area of the size distribution, while the absorption coefficient was formulated in terms of both the projected area and mass of the size distribution. These properties were formulated as explicit functions of ice crystal geometry and were not based on an 'effective radius.' Based on simulations of the second cirrus case study described in Part 1, absorption coefficients predicted in the near infrared for hexagonal columns and rosettes were up to 47% and 71% lower, respectively, than absorption coefficients predicted by using equivalent area spheres. This resulted in single scatter albedos in the near-infrared that were considerably greater than those predicted by the equivalent area sphere method. Reflectances in this region should therefore be underestimated using the equivalent area sphere approach. Cloud optical depth was found to depend on ice crystal habit. When the simulated cirrus cloud contained only bullet rosettes, the optical depth was 142% greater than when the cloud contained only hexagonal columns. This increase produced a doubling in cloud albedo. In the near-infrared (IR), the single scatter albedo also exhibited a significant dependence on ice crystal habit. More research is needed on the geometrical properties of ice crystals before the influence of ice crystal shape on cirrus radiative properties can be adequately understood. This study provides a way of coupling the radiative properties of absorption, extinction, and single scatter albedo to the microphysical properties of cirrus clouds. The dependence of extinction and absorption on ice crystal shape was not just due to geometrical differences between crystal types, but was also due to the effect these differences had on the evolution of ice particle size spectra. The ice particle growth model in Part 1 and the radiative properties treated here are based on analytical formulations, and thus represent a computationally efficient means of modeling the microphysical and radiative properties of cirrus clouds.
This work presents the results of an experimental study of ice particle impacts on a moving wedge. The experiment was conducted in the Adverse Environment Rotor Test Stand (AERTS) facility located at Penn State University. The wedge was placed at the tip of a rotating blade. Ice particles shot from a pressure gun intercepted the moving wedge and impacted it at a location along its circular path. The upward velocity of the ice particles varied from 7 to 12 meters per second. Wedge velocities were varied from 0 to 120 meters per second. Wedge angles tested were 0 deg, 30 deg, 45 deg, and 60 deg. High speed imaging combined with backlighting captured the impact allowing observation of the effect of velocity and wedge angle on the impact and the post-impact fragment behavior. It was found that the pressure gun and the rotating wedge could be synchronized to consistently obtain ice particle impacts on the target wedge. It was observed that the number of fragments increase with the normal component of the impact velocity. Particle fragments ejected immediately after impact showed velocities higher than the impact velocity. The results followed the major qualitative features observed by other researchers for hailstone impacts, even though the reduced scale size of the particles used in the present experiment as compared to hailstones was 4:1.
The behavior of spherical water-ice particles in interplanetary space is investigated theoretically using an equation which gives the sublimation rate as a function of the solar energy distribution at 1 AU, solar distance, the temperature of the ice-particle surface, the temperature of the surrounding background, particle absorption efficiency, particle radius, and the complex refractive index of ice. The variation of sublimation rate with particle size is examined, and it is shown that the sublimation rate depends critically on the absorptivity of ice in the near-IR spectral region. Two different spectral-absorption curves are plotted on the basis of two different choices for the complex refractive index of ice, and an attempt is made to decide experimentally which energy distribution is more appropriate. Results of direct laboratory measurements of the sublimation rate of smooth spherical ice particles are compared with the theoretical results obtained from the cited sublimation-rate equation. Agreement between theory and experiment is found to be excellent when the complex refractive indices given by Bertie et al., (1969) are employed. Sublimation rates and particle lifetimes are then computed as a function of particle radius for heliocentric distances of 0.5, 0.75, and 1.0 AU.
Abstract Bulk ice-microphysical models parameterize the dynamic evolution of ice particles from advection, collection, and sedimentation through a cloud layer to the surface. Frozen hydrometeors can grow to acquire a multitude of shapes and sizes, which influence the distribution of mass within cloud systems. Aggregates, defined herein as the collection of ice particles, have a variety of formations based on initial ice particle size, shape, falling orientation, and the number of particles that collect. This work focuses on using the Ice Particle and Aggregate Simulator (IPAS) as a statistical tool to repetitively collect ice crystals of identical properties to derive bulk aggregate characteristics. A database of 9 744 000 aggregates is generated with resulting properties analyzed. After 150 single ice crystals (monomers) collect, the most extreme aggregate aspect ratio calculations asymptote toward and ϕ ca ≈ 0.50 for aggregates composed of quasi-horizontally oriented and randomly oriented monomers, respectively. The results presented are largely consistent with both a previous theoretical study and estimates derived from ground-based observations from two different geographic locations. Particle falling orientation highly influences newly formed aggregate aspect ratios from the collection of particles with extreme aspect ratios; quasi-horizontally oriented particles can produce aggregate aspect ratios an order of magnitude more extreme than randomly oriented particles but can also produce near-spherical aggregates as the number of monomers comprising the aggregate reach approximately 100. Finally, a majority of collections result in aggregates that are closer to prolate than oblate spheroids.
A flow and ice particle trajectory analysis was performed for the booster of the Honeywell ALF502 engine. The analysis focused on two closely related conditions one of which produced an icing event and another which did not during testing of the ALF502 engine in the Propulsion Systems Lab (PSL) at NASA Glenn Research Center. The flow analysis was generated using the NASA Glenn GlennHT flow solver and the particle analysis was generated using the NASA Glenn LEWICE3D v3.63 ice accretion software. The inflow conditions for the two conditions were similar with the main differences being that the condition that produced the icing event was 6.8 K colder than the non-icing event case and the inflow ice water content (IWC) for the non-icing event case was 50% less than for the icing event case. The particle analysis, which considered sublimation, evaporation and phase change, was generated for a 5 micron ice particle with a sticky impact model and for a 24 micron median volume diameter (MVD), 7 bin ice particle distribution with a supercooled large droplet (SLD) splash model used to simulate ice particle breakup. The particle analysis did not consider the effect of the runback and re-impingement of water resulting from the heated spinner and anti-icing system. The results from the analysis showed that the amount of impingement for the components were similar for the same particle size and impact model for the icing and non-icing event conditions. This was attributed to the similar aerodynamic conditions in the booster for the two cases. The particle temperature and melt fraction were higher at the same location and particle size for the non-icing event than for the icing event case due to the higher incoming inflow temperature for the non-event case. The 5 micron ice particle case produced higher impact temperatures and higher melt fractions on the components downstream of the fan than the 24 micron MVD case because the average particle size generated by the particle breakup was larger than 5 microns which yielded less warming and melting. The analysis also showed that the melt fraction and wet bulb temperature icing criterion developed during tests in the Research Altitude Test Facility (RATFac) at the National Research Council (NRC) of Canada were useful in predicting icing events in the ALF502 engine. The development of an ice particle impact model which includes the effects of particle breakup, phase change, and surface state is necessary to further improve the prediction of ice particle transport with phase change through turbomachinery.