Search NASA⌕ Search

SEARCH · Search NASA

Results for “Ice crystal microphysics”

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 181 records · Page 10
All sourcesNASA NTRSDOE OSTI
NASA NTRS2007

Multi-Layer Arctic Mixed-Phase Clouds Simulated by a Cloud-Resolving Model: Comparison with ARM Observations and Sensitivity Experiments

A cloud-resolving model (CRM) is used to simulate the multiple-layer mixed-phase stratiform (MPS) clouds that occurred during a three-and-a-half day subperiod of the Department of Energy-Atmospheric Radiation Measurement Program s Mixed-Phase Arctic Cloud Experiment (M-PACE). The CRM is implemented with an advanced two-moment microphysics scheme, a state-of-the-art radiative transfer scheme, and a complicated third-order turbulence closure. Concurrent meteorological, aerosol, and ice nucleus measurements are used to initialize the CRM. The CRM is prescribed by time-varying large-scale advective tendencies of temperature and moisture and surface turbulent fluxes of sensible and latent heat. The CRM reproduces the occurrences of the single- and double-layer MPS clouds as revealed by the M-PACE observations. However, the simulated first cloud layer is lower and the second cloud layer thicker compared to observations. The magnitude of the simulated liquid water path agrees with that observed, but its temporal variation is more pronounced than that observed. As in an earlier study of single-layer cloud, the CRM also captures the major characteristics in the vertical distributions and temporal variations of liquid water content (LWC), total ice water content (IWC), droplet number concentration and ice crystal number concentration (nis) as suggested by the aircraft observations. However, the simulated mean values differ significantly from the observed. The magnitude of nis is especially underestimated by one order of magnitude. Sensitivity experiments suggest that the lower cloud layer is closely related to the surface fluxes of sensible and latent heat; the upper cloud layer is probably initialized by the large-scale advective cooling/moistening and maintained through the strong longwave (LW) radiative cooling near the cloud top which enhances the dynamical circulation; artificially turning off all ice-phase microphysical processes results in an increase in LWP by a factor of 3 due to interactions between the excessive LW radiative cooling and extra cloud water; heating caused by phase change of hydrometeors could affect the LWC and cloud top height by partially canceling out the LW radiative cooling. It is further shown that the resolved dynamical circulation appears to contribute more greatly to the evolution of the MPS cloud layers than the parameterized subgrid-scale circulation.

Luo, Yali↗
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 NTRS2002

Microphysics and Radiative Properties of Cirrus: Instrumentation and Analysis

Work under this grant has involved further development of a new aircraft instrument (the cloudscope) for real time characterization of atmospheric particulates together with field observations of such particulates, both in the form of ice and also as nuclei responsible for nucleation of both ice and water cloud particles. Part of the work involving assessment of the frequency of ice crystal shapes has been carried out in collaboration with the Meteorological Service of Canada; part of the work in a field program with the NCAR C-130. Part of the work has been interpreted in terms of laboratory simulation of ice crystal growth under a wide variety of conditions carried out under a grant from Physical Meteorology Program, National Science Foundation.

Hallett, John↗
DOE OSTI2020

The Role of Ice Splintering on Microphysics of Deep Convective Clouds Forming Under Different Aerosol Conditions: Simulations Using the Model With Spectral Bin Microphysics

Observations during the Ice in Clouds Experiment-Tropical (ICE-T) field experiment show that the ice particles concentration in a developing deep convective clouds at the level of T = -15 °C reached about 500 L -1 , that is, many orders higher than that of ice-nucleating particle. To simulate microphysics of these clouds, the 2-D Hebrew University Cloud model (HUCM) with spectral bin microphysics was used in which two main types of ice multiplication mechanisms were included in addition to the Hallet-Mossop mechanism. In the first ice multiplication mechanism ice splinters form by drop freezing and drop-ice collisions. Ice multiplication of this type dominates during developing stage of cloud evolution, when liquid water content is significant. At later stage when clouds become nearly glaciated, ice crystals are produced largely by ice splintering during ice-ice collisions (the second ice multiplication mechanism). Simulations show that droplet size distributions, as well as size distributions of ice particles, agree well with the measurements during ICE-T. Simulations with different cloud condensation nuclei concentrations show the existence of the “optimum” cloud condensation nuclei concentration (or droplet concentration), at which concentration of ice splinters reaches maximum. In these simulations ice nucleation caused by the direct formation of ice crystals upon ice-nucleating particles, as well as the Hallett-Mossop process, has a negligible contribution to the ice crystal concentration.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS1993

Cirrus properties deduced from CO2 lidar observations of zenith-enhanced backscatter from oriented crystals

Many lidar researchers have occasionally observed zenith-enhanced backscatter (ZEB) from middle and high clouds. The ZEB signature consists of strong backscatter when the lidar is pointed directly at zenith and a dramatic decline in backscatter as the zenith angle dips slightly off zenith. Mirror-like reflection from horizontal facets of oriented crystals (especially plates) is generally accepted as the cause. It was found during a 3-year observation program that approximately 50 percent of ice clouds had ZEB, regardless of cloud height. The orientation of crystals and the ZEB they cause are important to study and understand for several reasons. First, radiative transfer in clouds with oriented crystals is different than if the same particles were randomly oriented. Second, crystal growth depends partly on the orientation of the particles. Third, ZEB measurements may provide useful information about cirrus microphysical and radiative properties. Finally, the remarkable effect of ZEB on lidar signals should be understood in order to properly interpret lidar data.

Eberhard, Wynn L.↗
NASA NTRS1993

Inference of cirrus cloud properties using satellite-observed visible and infrared radiances. I - Parameterization of radiance fields

The impact of using phase functions for spherical droplets and hexagonal ice crystals to analyze radiances from cirrus is examined. Adding-doubling radiative transfer calculations are employed to compute radiances for different cloud thicknesses and heights over various backgrounds. These radiances are used to develop parameterizations of top-of-the-atmosphere visible reflectance and IR emittance using tables of reflectances as a function of cloud optical depth, viewing and illumination angles, and microphysics. This parameterization, which includes Rayleigh scattering, ozone absorption, variable cloud height, and an anisotropic surface reflectance, reproduces the computed top-of-the-atmosphere reflectances with an accruacy of +/- 6 percent for four microphysical models: 10-micron water droplet, small symmetric crystal, cirrostratus, and cirrus uncinus. The accuracy is twice that of previous models.

Minnis, Patrick↗
NASA NTRS2001

The Effects of Aircraft Wake Dynamics on Contrail Development

Results of large-eddy simulations of the development of young persistent ice contrails are presented, concentrating on the interactions between the aircraft wake dynamics and the ice cloud evolution over ages front a few seconds to approx. 30 min. The 3D unsteady evolution of the dispersing engine exhausts, trailing vortex pair interaction and breakup, and subsequent Brunt-Vaisala oscillations of the older wake plume are modeled in detail in high-resolution simulations, coupled with it bulk microphysics model for the contrail ice development. The simulations confirm that the early wake dynamics can have a strong influence on the properties of persistent contrails even at late times. The vortex dynamics are the primary determinant of the vertical extent of the contrail (until precipitate ton becomes significant): and this together with the local wind shear largely determines the horizontal extent. The ice density, ice crystal number density, and a conserved exhaust tracer all develop and disperse in different fashions from each other. The total ice crystal number can be significantly reduced due to adiabatic compression resulting from the downward motion of the vortex system, even for ambient conditions that are substantially supersaturated with respect to ice. The fraction of the initial ice crystals surviving, their spatial distribution and the ice mass distribution are all sensitive to the aircraft type, ambient humidity, assumed initial ice crystal number, and ambient turbulence conditions. There is a significant range of conditions for which a smaller transport such as a B737 produces as significant a persistent contrail as a larger transport such as a B747, even though the latter consumes almost five times as much fuel. The difficulties involved in trying to minimize persistent contrail production are discussed.

Lewellen, D. C.↗
NASA NTRS2008

Arctic Mixed-phase Clouds Simulated by a Cloud-Resolving Model: Comparison with ARM Observations and Sensitivity to Microphysics Parameterizations

Single-layer mixed-phase stratiform (MPS) Arctic clouds, which formed under conditions of large surface heat flux combined with general subsidence during a subperiod of the Atmospheric Radiation Measurement (ARM) Program Mixed-Phase Arctic Cloud Experiment (M-PACE), are simulated with a cloud resolving model (CRM). The CRM is implemented with either an advanced two-moment (M05) or a commonly used one-moment (L83) bulk microphysics scheme and a state-of-the-art radiative transfer scheme. The CONTROL simulation, that uses the M05 scheme and observed aerosol size distribution and ice nulei (IN) number concentration, reproduces the magnitudes and vertical structures of cloud liquid water content (LWC), total ice water content (IWC), number concentration and effective radius of cloud droplets as suggested by the M-PACE observations. It underestimates ice crystal number concentrations by an order of magnitude and overestimates effective radius of ice crystals by a factor of 2-3. The OneM experiment, that uses the L83 scheme, produces values of liquid water path (LWP) and ice plus snow water path (ISWP) that were about 30% and 4 times, respectively, of those produced by the CONTROL. Its vertical profile of IWC exhibits a bimodal distribution in contrast to the constant distribution of IWC produced in the CONTROL and observations.

Xu, Kuan-Man↗
NASA NTRS2009

Diagnosing the Ice Crystal Enhancement Factor in the Tropics

Recent modeling studies have revealed that ice crystal number concentration is one of the dominant factors in the effect of clouds on radiation. Since the ice crystal enhancement factor and ice nuclei concentration determine the concentration, they are both important in quantifying the contribution of increased ice nuclei to global warming. In this study, long-term cloud-resolving model (CRM) simulations are compared with field observations to estimate the ice crystal enhancement factor in tropical and midlatitudinal clouds, respectively. It is found that the factor in tropical clouds is ~10 3-104 times larger than that of mid-latitudinal ones, which makes physical sense because entrainment and detrainment in the Tropics are much stronger than in middle latitudes. The effect of entrainment/detrainment on the enhancement factor, especially in tropical clouds, suggests that cloud microphysical parameterizations should be coupled with subgrid turbulence parameterizations within CRMs to obtain a more accurate depiction of cloud-radiative forcing.

Zeng, Xiping↗
NASA NTRS2015

The Dependence of Cirrus Gamma Size Distributions Expressed as Volumes in N(sub 0)-Lambda-Mu Phase Space and Bulk Cloud Properties on Environmental Conditions: Results from the Small Ice Particles in Cirrus Experiment (SPARTICUS)

The variability of cirrus ice microphysical properties is investigated using observations obtained during the Small Particles in Cirrus (SPARTICUS) campaign. An existing approach that represents a size distribution (SD) as a single gamma function using an ellipsoid of equally realizable solutions in (N(sub 0), lambda, mu) phase space is modified to automatically identify multiple modes in SDs and characterize each mode by such an ellipsoid. The modified approach is applied to ice crystals with maximum dimension D greater than15 micrometers collected by the 2-D stereo and 2-D precipitation probes on the Stratton Park Engineering Company Learjet. The dependencies of N(sub 0), mu, and lambda from each mode, total number concentration, bulk extinction, ice water content (IWC), and mass median maximum dimension D(sub mm) as a function of temperature T and cirrus type are then analyzed. The changes in the observed codependencies between N(sub 0), mu, and lambda, bulk extinction, IWC, and D(sub mm) with environmental conditions indicate that particles were larger at higher T during SPARTICUS. At most two modes were observed in any SD during SPARTICUS, with the average boundary between them at 115 micrometers, similar to past studies not using probes with shatter mitigating tips and artifact removal algorithms. The bimodality of the SDs increased with T. This and the differences in N(sub 0), mu, and lambda between the modes suggest that particles with smaller D nucleated more recently than particles with larger D, which grew via vapor deposition and aggregation. Because smaller crystals, whose concentrations are uncertain, make marginal contributions to higher order moments, the use of higher moments for evaluating model fields is suggested.

Ellipsoids↗
NASA NTRS1992

Cirrus microphysics and radiative transfer - Cloud field study on 28 October 1986

An analysis of remote-sensing measurements which had detected an inhomogeneous cloud structure, at 6-11 km altitude, in the data acquired during the FIRE 86 experiment in a 75 x 50 km cirrus cloud field, has derived fluxes whose comparison with modeled fluxes imply a modeling-based underestimation of both solar reflectivity/attenuation and downward IR fluxes. Reconciliation of model results with measurements can be achieved either by adding large concentrations of ice crystals or by altering the backscattering properties of ice crystals.

Kinne, Stefan↗
NASA NTRS2013

Saharan Dust Event Impacts on Cloud Formation and Radiation over Western Europe

We investigated the impact of mineral dust particles on clouds, radiation and atmospheric state during a strong Saharan dust event over Europe in May 2008, applying a comprehensive online-coupled regional model framework that explicitly treats particle-microphysics and chemical composition. Sophisticated parameterizations for aerosol activation and ice nucleation, together with two-moment cloud microphysics are used to calculate the interaction of the different particles with clouds depending on their physical and chemical properties. The impact of dust on cloud droplet number concentration was found to be low, with just a slight increase in cloud droplet number concentration for both uncoated and coated dust. For temperatures lower than the level of homogeneous freezing, no significant impact of dust on the number and mass concentration of ice crystals was found, though the concentration of frozen dust particles reached up to 100 l-1 during the ice nucleation events. Mineral dust particles were found to have the largest impact on clouds in a temperature range between freezing level and the level of homogeneous freezing, where they determined the number concentration of ice crystals due to efficient heterogeneous freezing of the dust particles and modified the glaciation of mixed phase clouds. Our simulations show that during the dust events, ice crystals concentrations were increased twofold in this temperature range (compared to if dust interactions are neglected). This had a significant impact on the cloud optical properties, causing a reduction in the incoming short-wave radiation at the surface up to -75Wm-2. Including the direct interaction of dust with radiation caused an additional reduction in the incoming short-wave radiation by 40 to 80Wm-2, and the incoming long-wave radiation at the surface was increased significantly in the order of +10Wm-2. The strong radiative forcings associated with dust caused a reduction in surface temperature in the order of -0.2 to -0.5K for most parts of France, Germany, and Italy during the dust event. The maximum difference in surface temperature was found in the East of France, the Benelux, and Western Germany with up to -1 K. This magnitude of temperature change was sufficient to explain a systematic bias in numerical weather forecasts during the period of the dust event.

Bangert, M.↗
NASA NTRS2008

On the Sensitivity of Atmospheric Ensembles to Cloud Microphysics in Long-Term Cloud-Resolving Model Simulations

Month-long large-scale forcing data from two field campaigns are used to drive a cloud-resolving model (CRM) and produce ensemble simulations of clouds and precipitation. Observational data are then used to evaluate the model results. To improve the model results, a new parameterization of the Bergeron process is proposed that incorporates the number concentration of ice nuclei (IN). Numerical simulations reveal that atmospheric ensembles are sensitive to IN concentration and ice crystal multiplication. Two- (2D) and three-dimensional (3D) simulations are carried out to address the sensitivity of atmospheric ensembles to model dimensionality. It is found that the ensembles with high IN concentration are more sensitive to dimensionality than those with low IN concentration. Both the analytic solutions of linear dry models and the CRM output show that there are more convective cores with stronger updrafts in 3D simulations than in 2D, which explains the differing sensitivity of the ensembles to dimensionality at different IN concentrations.

Zeng, Xiping↗
NASA NTRS2011

Dynamical States of Low Temperature Cirrus

Low ice crystal concentration and sustained in-cloud supersaturation, commonly found in cloud observations at low temperature, challenge our understanding of cirrus formation. Heterogeneous freezing from effloresced ammonium sulfate, glassy aerosol, dust and black carbon are proposed to cause these phenomena; this requires low updrafts for cirrus characteristics to agree with observations and is at odds with the gravity wave spectrum in the upper troposphere. Background temperature fluctuations however can establish a dynamical equilibrium between ice production and sedimentation loss (as opposed to ice crystal formation during the first stages of cloud evolution and subsequent slow cloud decay) that explains low temperature cirrus properties. This newly-discovered state is favored at low temperatures and does not require heterogeneous nucleation to occur (the presence of ice nuclei can however facilitate its onset). Our understanding of cirrus clouds and their role in anthropogenic climate change is reshaped, as the type of dynamical forcing will set these clouds in one of two preferred microphysical regimes with very different susceptibility to aerosol.

Barahona, D.↗
NASA NTRS2010

Comparing Aircraft Observations of Snowfall to Forecasts Using Single or Two Moment Bulk Water Microphysics Schemes

High resolution weather forecast models with explicit prediction of hydrometeor type, size distribution, and fall speed may be useful in the development of precipitation retrievals, by providing representative characteristics of frozen hydrometeors. Several single or double-moment microphysics schemes are currently available within the Weather Research and Forecasting (WRF) model, allowing for the prediction of up to three ice species. Each scheme incorporates different assumptions regarding the characteristics of their ice classes, particularly in terms of size distribution, density, and fall speed. In addition to the prediction of hydrometeor content, these schemes must accurately represent the vertical profile of water vapor to account for possible attenuation, along with the size distribution, density, and shape characteristics of ice crystals that are relevant to microwave scattering. An evaluation of a particular scheme requires the availability of field campaign measurements. The Canadian CloudSat/CALIPSO Validation Project (C3VP) obtained measurements of ice crystal shapes, size distributions, fall speeds, and precipitation during several intensive observation periods. In this study, C3VP observations obtained during the 22 January 2007 synoptic-scale snowfall event are compared against WRF model output, based upon forecasts using four single-moment and two double-moment schemes available as of version 3.1. Schemes are compared against aircraft observations by examining differences in size distribution, density, and content. In addition to direct measurements from aircraft probes, simulated precipitation can also be converted to equivalent, remotely sensed characteristics through the use of the NASA Goddard Satellite Data Simulator Unit. Outputs from high resolution forecasts are compared against radar and satellite observations emphasizing differences in assumed crystal shape and size distribution characteristics.

Molthan, Andrew L.↗
NASA NTRS2021

Comparisons of Cloud In-Situ Microphysical Properties of Deep Convective Clouds to Appendix D/P using Data from the HAIC-HIWC and HIWC-RADAR I Flight Campaigns

In-situ cloud data from three international flight campaigns are compared to the Federal Aviation Administration Title 14 Code of Federal Regulations Part 33 Appendix D mixed-phase/glaciated environmental envelope, and the corresponding identical European Aviation Safety Agency CS-25 Appendix P envelope. The appendices consist of a temperature-altitude envelope, a 99th percentile total water content envelope at the 17.4 Nm distance scale, a distance factor for estimation at other distance scales, ice crystal median mass diameter, and recommended liquid water content levels in mixed-phase icing conditions. The data were collected during 54 flights out of one subtropical and two tropical locations, with 472 runs from about 17,000’ to 39,000’ in approximately 115 clouds. The campaigns provide about 29,600 Nm of in situ data in deep convection over four targeted temperature intervals: -10, -30, -40, and -50, all ± 5 C. The dataset is a modern and unique documentation of the ice crystal icing environment, and results described in this article will contribute to regulatory and industry assessment of Appendices D and P.

Ice Crystal Icing↗
NASA NTRS2012

Sensitivity of Cirrus and Mixed-phase Clouds to the Ice Nuclei Spectra in McRAS-AC: Single Column Model Simulations

The salient features of mixed-phase and ice clouds in a GCM cloud scheme are examined using the ice formation parameterizations of Liu and Penner (LP) and Barahona and Nenes (BN). The performance of LP and BN ice nucleation parameterizations were assessed in the GEOS-5 AGCM using the McRAS-AC cloud microphysics framework in single column mode. Four dimensional assimilated data from the intensive observation period of ARM TWP-ICE campaign was used to drive the fluxes and lateral forcing. Simulation experiments where established to test the impact of each parameterization in the resulting cloud fields. Three commonly used IN spectra were utilized in the BN parameterization to described the availability of IN for heterogeneous ice nucleation. The results show large similarities in the cirrus cloud regime between all the schemes tested, in which ice crystal concentrations were within a factor of 10 regardless of the parameterization used. In mixed-phase clouds there are some persistent differences in cloud particle number concentration and size, as well as in cloud fraction, ice water mixing ratio, and ice water path. Contact freezing in the simulated mixed-phase clouds contributed to transfer liquid to ice efficiently, so that on average, the clouds were fully glaciated at T approximately 260K, irrespective of the ice nucleation parameterization used. Comparison of simulated ice water path to available satellite derived observations were also performed, finding that all the schemes tested with the BN parameterization predicted 20 average values of IWP within plus or minus 15% of the observations.

Betancourt, R. Morales↗
NASA NTRS2013

Evaluation of Model Microphysics Within Precipitation Bands of Extratropical Cyclones

Recent studies evaluating the bulk microphysical schemes (BMPs) within cloud resolving models (CRMs) have indicated large uncertainties and errors in the amount and size distributions of snow and cloud ice aloft. The snow prediction is sensitive to the snow densities, habits, and degree of riming within the BMPs. Improving these BMPs is a crucial step toward improving both weather forecasting and climate predictions. Several microphysical schemes in the Weather Research and Forecasting (WRF) model down to 1.33‐km grid spacing are evaluated using aircraft, radar, and ground in situ data from the Global Precipitation Mission Coldseason Precipitation Experiment (GCPEx) experiment, as well as a few years (15 winter storms) of surface measurements of riming, crystal habit, snow density, and radar measurements at Stony Brook, NY (SBNY on north shore of Long Island) during the 2009-2012 winter seasons. Surface microphysical measurements at SBNY were taken every 15 to 30 minutes using a stereo microscope and camera, and snow depth and snow density were also recorded. During these storms, a vertically-pointing Ku-band radar was used to observe the vertical evolution of reflectivity and Doppler vertical velocities. A Particle Size and Velocity (PARSIVEL) disdrometer was also used to measure the surface size distribution and fall speeds of snow at SBNY. For the 15 cases at SBNY, the WSM6, Morrison (MORR), Thompson (THOM2), and Stony Brook (SBU-YLIN) BMPs were validated. A non-spherical snow assumption (THOM2 and SBU-YLIN) simulated a more realistic distribution of reflectivity than spherical snow assumptions in the WSM6 and MORR schemes. The MORR, WSM6, and SBU-YLIN schemes are comparable to the observed velocity distribution in light and moderate riming periods. The THOM2 is ~0.25 meters per second too slow with its velocity distribution in these periods. In heavier riming, the vertical Doppler velocities in the WSM6, THOM2, and MORR schemes were ~0.25 meters per second too slow, while the SBU-YLIN was 0.25 to 0.5 meters per second too fast. Overall, the BMPs simulate a size distribution close to the observed for D < 4 mm in the dendritic, plates, and mixed habit periods. The model BMPs underestimate the size distribution when large aggregates were observed. For D > 6 mm in the dendrites, side planes, and mixed habit periods, the BMPs are likely not simulating enough aggregation to create a larger size distribution, although the MORR (double moment) scheme seemed to perform best. These SBNY results will be compared with some results from GCPEx for a warm frontal snow band observed at 18 February 2012.

Colle, Brian A.↗