skip to main content


Title: Oriented Ice Crystals: A Single-Scattering Property Database for Applications to Lidar and Optical Phenomenon Simulations
Abstract

A database (TAMUoic2019) of the scattering, absorption, and polarization properties of horizontally oriented hexagonal plates (HOPs) and horizontally oriented hexagonal columns (HOCs) at three wavelengths (355, 532, and 1064 nm) is developed for applications to radiative transfer simulations and remote sensing implementations involving oriented ice crystals. The maximum dimension of oriented ice crystals ranges from 50 to 10 000 μm in 165 discrete size bins. The database accounts for 94 incident directions. The single-scattering properties of oriented ice crystals are computed with the physical-geometric optics method (PGOM), which is consistent with the invariant-imbedding T-matrix method for particles with size parameters larger than approximately 100–150. Note that the accuracy of PGOM increases as the size parameter increases. PGOM computes the two-dimensional phase matrix as a function of scattering polar and azimuth angles, and the phase matrix significantly varies with the incident direction. To derive the bulk optical properties of ice clouds for practical radiative transfer applications, the optical properties of individual HOPs and HOCs are averaged over the probability distribution of the tilting angle of oriented ice crystals based on the use of the TAMUoic2019 database. Simulations of lidar signals associated with ice clouds based on the bulk optical properties indicate the importance of the fraction of oriented ice crystals and the probability distribution of the tilting angle. Simulations of optical phenomena caused by oriented ice crystals demonstrate that the computed single-scattering properties of oriented ice crystals are physically rational.

 
more » « less
Award ID(s):
1826936
NSF-PAR ID:
10109581
Author(s) / Creator(s):
 ;  
Publisher / Repository:
American Meteorological Society
Date Published:
Journal Name:
Journal of the Atmospheric Sciences
Volume:
76
Issue:
9
ISSN:
0022-4928
Page Range / eLocation ID:
p. 2635-2652
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract

    Ice crystals commonly adopt a horizontal orientation under certain aerodynamic and electrodynamic conditions that occur in the atmosphere. While the radiative impact of horizontally oriented ice crystals (HOIC) has been theoretically studied with respect to their impact on shortwave cloud albedo, the longwave impact remains unexplored. This work analyzes the occurrence of HOIC at Summit, Greenland, from July 2015 to June 2017. Using polarization lidar and ancillary atmospheric sensors, ice crystal orientations are identified and used to interpret cloud radiative impact on the surface radiation budget. We find HOIC occur in at least 25.6% of all ice‐only column observations. We find that the shortwave impact of HOIC is to increase cloud radiative effect by approximately 22% for a given solar zenith angle. We also find that the longwave impact of HOIC compared to randomly oriented ice crystals are statistically different at the p < 0.01 significance level, increasing the surface radiative effect by approximately 8% for clouds with infrared optical depths < ~1. We suggest that the observed difference between the surface radiative effect for clouds containing randomly oriented ice crystals and HOIC may be due to enhanced scattering, but this hypothesis needs to be further explored with more detailed observations and modeling.

     
    more » « less
  2. null (Ed.)
    Abstract. Ice crystal submicron structures have a large impact on the opticalproperties of cirrus clouds and consequently on their radiative effect.Although there is growing evidence that atmospheric ice crystals are rarelypristine, direct in situ observations of the degree of ice crystal complexityare largely missing. Here we show a comprehensive in situ data set of icecrystal complexity coupled with measurements of the cloud angular scatteringfunctions collected during a number of observational airborne campaigns atdiverse geographical locations. Our results demonstrate that an overwhelmingfraction (between 61 % and 81 %) of atmospheric ice crystals sampledin the different regions contain mesoscopic deformations and, as aconsequence, a similar flat and featureless angular scattering function isobserved. A comparison between the measurements and a database of opticalparticle properties showed that severely roughened hexagonal aggregatesoptimally represent the measurements in the observed angular range. Based onthis optical model, a new parameterization of the cloud bulk asymmetry factorwas introduced and its effects were tested in a global climate model. Themodelling results suggest that, due to ice crystal complexity, ice-containingclouds can induce an additional short-wave cooling effect of−1.12 W m2 on the top-of-the-atmosphere radiative budget that hasnot yet been considered. 
    more » « less
  3. The microstructure of the uppermost portions of a melting Arctic sea ice cover has a disproportionately large influence on how incident sunlight is reflected and absorbed in the ice/ocean system. The surface scattering layer (SSL) effectively backscatters solar radiation and keeps the surface albedo of melting ice relatively high compared to ice with the SSL manually removed. Measurements of albedo provide information on how incoming shortwave radiation is partitioned by the SSL and have been pivotal to improving climate model parameterizations. However, the relationship between the physical and optical properties of the SSL is still poorly constrained. Until now, radiative transfer models have been the only way to infer the microstructure of the SSL. During the MOSAiC expedition of 2019–2020, we took samples and, for the first time, directly measured the microstructure of the SSL on bare sea ice using X-ray micro-computed tomography. We show that the SSL has a highly anisotropic, coarse, and porous structure, with a small optical diameter and density at the surface, increasing with depth. As the melting surface ablates, the SSL regenerates, maintaining some aspects of its microstructure throughout the melt season. We used the microstructure measurements with a radiative transfer model to improve our understanding of the relationship between physical properties and optical properties at 850 nm wavelength. When the microstructure is used as model input, we see a 10–15% overestimation of the reflectance at 850 nm. This comparison suggests that either a) spatial variability at the meter scale is important for the two in situ optical measurements and therefore a larger sample size is needed to represent the microstructure or b) future work should investigate either i) using a ray-tracing approach instead of explicitly solving the radiative transfer equation or ii) using a more appropriate radiative transfer model.

     
    more » « less
  4. Abstract. Regions with high ice water content (HIWC), composed of mainly small ice crystals, frequently occur over convective clouds in the tropics. Such regions can have median mass diameters (MMDs) <300 µm and equivalent radar reflectivities <20 dBZ. To explore formation mechanisms for these HIWCs, high-resolution simulations of tropical convective clouds observed on 26 May 2015 during the High Altitude Ice Crystals – High Ice Water Content (HAIC-HIWC) international field campaign based out of Cayenne, French Guiana, are conducted using the Weather Research and Forecasting (WRF) model with four different bulk microphysics schemes: the WRF single‐moment 6‐class microphysics scheme (WSM6), the Morrison scheme, and the Predicted Particle Properties (P3) scheme with one- and two-ice options. The simulations are evaluated against data from airborne radar and multiple cloud microphysics probes installed on the French Falcon 20 and Canadian National Research Council (NRC) Convair 580 sampling clouds at different heights. WRF simulations with different microphysics schemes generally reproduce the vertical profiles of temperature, dew-point temperature, and winds during this event compared with radiosonde data, and the coverage and evolution of this tropical convective system compared to satellite retrievals. All of the simulations overestimate the intensity and spatial extent of radar reflectivity by over 30 % above the melting layer compared to the airborne X-band radar reflectivity data. They also miss the peak of the observed ice number distribution function for 0.1 more » « less
  5. Abstract. Scattering codes are used to study the optical properties of polar stratospheric clouds (PSCs). Particle backscattering and depolarization coefficients can be computed with available scattering codes once the particle size distribution (PSD) is known and a suitable refractive index is assumed. However, PSCs often appear as external mixtures of supercooled ternary solution (STS) droplets, solid nitric acid trihydrate (NAT) and possibly ice particles, making the assumption of a single refractive index and a single morphology to model the scatterers questionable.Here we consider a set of 15 coincident measurements of PSCs above McMurdo Station, Antarctica, using ground-based lidar, a balloon-borne optical particle counter (OPC) and in situ observations taken by a laser backscattersonde and OPC during four balloon stratospheric flights from Kiruna, Sweden. This unique dataset of microphysical and optical observations allows us to test the performances of optical scattering models when both spherical and aspherical scatterers of different composition and, possibly, shapes are present. We consider particles as STS if their radius is below a certain threshold value Rth and NAT or possibly ice if it is above it. The refractive indices are assumed known from the literature. Mie scattering is used for the STS, assumed spherical. Scattering from NAT particles, considered spheroids of different aspect ratio (AR), is treated with T-matrix results where applicable. The geometric-optics–integral-equation approach is used whenever the particle size parameter is too large to allow for a convergence of the T-matrix method.The parameters Rth and AR of our model have been varied between 0.1 and 2 µm and between 0.3 and 3, respectively, and the calculated backscattering coefficient and depolarization were compared with the observed ones. The best agreement was found for Rth between 0.5 and 0.8 µm and for AR less than 0.55 and greater than 1.5.To further constrain the variability of AR within the identified intervals, we have sought an agreement with the experimental data by varying AR on a case-by-case basis and further optimizing the agreement by a proper choice of AR smaller than 0.55 and greater than 1.5 and Rth within the interval 0.5 and 0.8 µm. The ARs identified in this way cluster around the values 0.5 and 2.5.The comparison of the calculations with the measurements is presented and discussed. The results of this work help to set limits to the variability of the dimensions and asphericity of PSC solid particles, within the limits of applicability of our model based on the T-matrix theory of scattering and on assumptions on a common particle shape in a PSD and a common threshold radius for all the PSDs. 
    more » « less