Linear particle depolarization ratio is presented for three case studies from the NASA Langley airborne High Spectral Resolution Lidar-2 (HSRL-2). Particle depolarization ratio from lidar is an indicator of non-spherical particles and is sensitive to the fraction of non-spherical particles and their size. The HSRL-2 instrument measures depolarization at three wavelengths: 355, 532, and 1064 nm. The three measurement cases presented here include two cases of dust-dominated aerosol and one case of smoke aerosol. These cases have partial analogs in earlier HSRL-1 depolarization measurements at 532 and 1064 nm and in literature, but the availability of three wavelengths gives additional insight into different scenarios for non-spherical particles in the atmosphere. A case of transported Saharan dust has a spectral dependence with a peak of 0.30 at 532 nm with smaller particle depolarization ratios of 0.27 and 0.25 at 1064 and 355 nm, respectively. A case of aerosol containing locally generated wind-blown North American dust has a maximum of 0.38 at 1064 nm, decreasing to 0.37 and 0.24 at 532 and 355 nm, respectively. The cause of the maximum at 1064 nm is inferred to be very large particles that have not settled out of the dust layer. The smoke layer has the opposite spectral dependence, with the peak of 0.24 at 355 nm, decreasing to 0.09 and 0.02 at 532 and 1064 nm, respectively. The depolarization in the smoke case may be explained by the presence of coated soot aggregates. We note that in these specific case studies, the linear particle depolarization ratio for smoke and dust-dominated aerosol are more similar at 355 nm than at 532 nm, having possible implications for using the particle depolarization ratio at a single wavelength for aerosol typing.
The impact of aerosols on climate depends on their horizontal and vertical distribution and microphysical properties. Lidar is an important tool for remote sensing of aerosol, because it provides vertically resolved information on aerosol abundance and aerosol type. One extremely useful lidar aerosol measurement is the linear particle depolarization ratio, an indicator of non-spherical particles. Polarization lidar is a large and active field, with recent contributions from ground-based networks, such as the European Aerosol Research Lidar Network (EARLINET; Pappalardo et al., 2014; Mamouri and Ansmann, 2014; Nisantzi et al., 2014) and the National Institute of Environmental Studies (NIES) East Asian network of lidars (Sugimoto et al., 2005; Nishizawa et al., 2011), and directed field campaigns, such as the Saharan Mineral Dust Experiment (SAMUM; Freudenthaler et al., 2009; Tesche et al., 2011) and the Saharan Aerosol Long-range Transport and Aerosol-Cloud Experiment (SALTRACE; Groß et al., 2015; Haarig et al., 2015); and others.
There is also considerable interest in global lidar observations from satellites. Global lidar observations of aerosol have been provided by the Cloud–Aerosol Lidar and Infrared Pathfinder Satellite Observations (CALIPSO) satellite since 2006 (Winker et al., 2007). Another satellite lidar, the experimental Cloud–Aerosol Transport System (CATS) instrument on the International Space Station (ISS) (McGill et al., 2012) was recently launched in January 2015, and the Earth Clouds Aerosols and Radiation Explorer (EarthCARE) satellite (Illingworth et al., 2015) is due to launch in 2018. CALIPSO linear particle depolarization ratio data have been used, for example, to assess the global distribution and transport of dust (e.g., Johnson et al., 2012; Liu et al., 2013; Yang et al., 2013). This measurement will also be part of the suite of measurements made by the ATLID (atmospheric lidar) on EarthCARE; however, CALIPSO measures depolarization at 532 nm and ATLID will measure it at 355 nm (Groß et al., 2014; Illingworth et al., 2015).
NASA Langley airborne High Spectral Resolution Lidars, HSRL-1 and HSRL-2, have participated in many process-oriented field campaigns, have provided validation and calibration data for CALIPSO since 2006 (Rogers et al., 2011, 2014), and will also be useful for validating the EarthCARE lidar measurements. Since the airborne HSRL-2 measures the particle depolarization ratio at both the CALIPSO and EarthCARE wavelengths and also at 1064 nm, observations from this instrument are useful for assessing how the measurements from the two satellite instruments will correspond. NASA's airborne HSRL-2 is the first HSRL system making depolarization measurements at three wavelengths. A ground-based Raman system operated by the Leibniz Institute of Tropospheric Research has also been recently upgraded to make three-wavelength depolarization measurements (Haarig et al., 2015).
Aerosol classification is one specific application of aerosol polarization measurements (Burton et al., 2012, 2013; Groß et al., 2013, 2014). The aerosol particle depolarization ratio from lidar is of key importance for the detection and assessment of dust and volcanic ash since it is a clear indicator of non-spherical particles. The particle depolarization ratio is also used to infer the amount of dust or ash in a mixture (Sugimoto and Lee, 2006; Tesche et al., 2009a, 2011; Ansmann et al., 2011, 2012; David et al., 2013; Burton et al., 2014; Mamouri and Ansmann, 2014). It is also sensitive to the size of the non-spherical particles (Ansmann et al., 2009; Sakai et al., 2010; Gasteiger et al., 2011; Gasteiger and Freudenthaler, 2014).
While a significant amount of study has been made of depolarization by dust and ash, smoke has also been observed to produce significant depolarization of lidar light in some cases (e.g., Fiebig et al., 2002; Sassen and Khvorostyanov, 2008; Sugimoto et al., 2010; Dahlkötter et al., 2014), but not in others (e.g., Müller et al., 2005). Even for cases with significant depolarization, the depolarization signature for smoke is generally smaller than for dust, at the wavelengths of 532 and 1064 nm where most lidar depolarization measurements of smoke have been made.
We will describe two dust-dominated cases and a smoke-dominated case where depolarizing aerosol was observed simultaneously at three wavelengths by the NASA Langley airborne HSRL-2 instrument. We show consistency between the three HSRL-2 cases and three previously published cases from the predecessor HSRL-1 instrument in which similar measurements were made at 532 and 1064 nm, and we also discuss similarities and differences with published lidar measurements globally. We find that the three cases each have a different spectral dependence of the particle depolarization ratio. Accordingly, we discuss possible explanations for these differences with reference to published studies. We also point out implications for future space-based observations of aerosol depolarization. We begin in Sect. 2 with a description of the NASA Langley airborne HSRL instruments and the methodology for polarization measurements, including an assessment of systematic uncertainty. In Sect. 3 we describe and discuss the dust cases and in Sect. 4 we describe and discuss the smoke case. We summarize the discussion and conclude in Sect. 5. In the Appendix we give more details about the estimation of systematic uncertainty.
The NASA Langley second-generation airborne HSRL-2 uses the HSRL technique (Shipley et al., 1983) to independently
measure aerosol extinction and backscatter at 355 and 532 nm and the
standard backscatter technique (Fernald, 1984) to measure aerosol
backscatter at 1064 nm. It also measures linear depolarization ratio at all
three wavelengths. It is a follow-on to the successful airborne HSRL-1
instrument (Hair et al., 2008), which has made measurements at 532 and
1064 nm since 2006 (Rogers et al., 2009). For measurements at 532 and
1064 nm, HSRL-2 is essentially identical to HSRL-1. HSRL measurements of
extinction and backscatter at 355 nm are made using an interferometer rather
than an iodine filter. For 355 nm measurements of depolarization discussed
here, the setup is very similar to the other channels; the small differences
are explained in Sect. 2.1. Data are sampled at 0.5 s temporal and 30 m
vertical resolutions. Aerosol backscatter and depolarization products are
averaged 10 s horizontally (
In this paper, we will focus on the measurements of the linear particle
depolarization ratio. Figure 1 shows a simplified diagram of the optics of
the transmission system that are relevant to the measurement of
depolarization. The primary optical components for the polarization of the
transmitted beams are Glan laser polarizers, which have a very high
polarization transmittance ratio of
A simplified block diagram of the transmitter optics relevant for linear depolarization measurements by HSRL-2. The two Glan laser polarizers are the primary components ensuring the transmitted laser light is polarized. The motorized calibration wave plates are used to align the output polarization to the receiver polarization analyzers. The 532 and 1064 nm laser beams are maintained as a single beam. Since they do not exit the laser at the same polarization, a co-alignment wave plate is used to align the polarization from the two wavelengths so that the Glan laser polarizer does not significantly reduce the amount of light transmitted at one of the wavelengths. The attenuator wave plate is used to attenuate the 532 nm beam for eye safety considerations when flying at low altitudes, and for maximizing the power output otherwise.
A simplified block diagram of the receiver optics relevant for
linear depolarization measurements by HSRL-2. The abbreviations are as follows:
PMT is photomultiplier tube; APD is silicon avalanche photodetector; DS is dichroic beam splitter; PBS is polarizing beam splitter;
co-pol is co-polarized channel (with respect to the transmitted light); X-pol is cross-polarized channel. The collimated light arrives from the telescope and
is split into the three wavelengths using dichroic beam splitters. The first
optical component filters solar background using either an interference
filter (indicated “Filter”) or an interference filter and etalon in
combination (indicated “Filter
The receiver optics relevant to depolarization measurements are shown in Fig. 2. The collimated light arrives from the telescope and is split into
the three wavelengths using dichroic beam splitters. Each beam is then
passed through an interference filter (1064 nm) or a combination of
interference filter and etalon (355 and 532 nm) to remove background
scattering. The effective full-width half-max (FWHM) bandwidths for the
three channels are 0.4 nm (3.5 cm
The volume (or total) linear depolarization ratio is the ratio of the signal
in the cross-polarized channel to that in the co-polarized channel,
normalized by the measured polarization gain ratio.
The particle depolarization ratio is calculated from the volume
depolarization ratio using the following (Cairo et al., 1999):
Systematic error can be a concern for polarization measurements. Potential sources of systematic error in volume depolarization ratio arise in the polarization optics and calibration. The retrieval of the particle depolarization ratio can potentially introduce additional systematic error related to the total aerosol scattering ratio or uncertainty in the molecular depolarization ratio value. We will provide an overview of the potential systematic errors here, including systematic uncertainty for volume depolarization ratio and propagated systematic uncertainty for the particle depolarization ratio. More details about these potential errors and the means of estimating the systematic uncertainty are given in the Appendix.
The linear volume depolarization ratio, given by Eq. (1), is the more basic measurement. Systematic errors in the volume depolarization ratio can arise from various sources, including calibration errors either in the polarization angle calibration or the polarization gain ratio calibration. A major concern for the measurement of depolarization is the potential for cross talk, which can arise from a number of sources, including imperfect polarization angle alignment, signal impurities due to imperfections in the polarization beam splitter (particularly the reflected channel), or other optics, including the aircraft window. Considering these sources, we estimate the systematic uncertainty in the volume depolarization ratio measurement to be the larger of 4.7 % (relative) or 0.001 (absolute) in the 355 nm channel, the larger of 5 % (relative) or 0.007 (absolute) in the 532 nm channel, and the larger of 2.6 % (relative) or 0.007 (absolute) in the 1064 nm channel. Further discussion of these estimates is given in the Appendix.
As can be seen in Eq. (2), the particle depolarization ratio,
The estimates given above are intended to be conservative estimates of the
systematic uncertainty confidence limit, such that we expect a high
probability that the systematic error is less than this value. The
systematic uncertainties on the three quantities,
Curtains of aerosol backscatter and extinction coefficients from HSRL-2 for observations on 13 July 2014 for a flight segment in Missouri and Kansas in the midwestern United States.
Measured properties for specific dust and smoke samples. To obtain
these values, samples were taken at specific times and altitudes comprising
400–4500 distinct measurement points. For the dust cases, values were chosen
near the peak value of the 532 nm particle depolarization ratio, where it
can be inferred that the aerosol is nearly pure dust. The values are
reported as mean
In this section we discuss two case studies in which HSRL-2 made three-wavelength measurements of the depolarization of dust.
On 13 July 2014, HSRL-2 aboard the B200 made measurements at three
wavelengths on a transit flight from Virginia to Colorado for the
DISCOVER-AQ field mission (
The linear particle depolarization ratio at three wavelengths measured by the HSRL-2 for the same flight segment shown in Fig. 3.
Figure 5 shows both the particle depolarization ratio and volume depolarization ratio measurements for all altitudes at 17:12 UT (17.2 UT). The particle depolarization ratio random and systematic uncertainties are also shown.
The predecessor of the HSRL-2, the NASA Langley HSRL-1 instrument, observed
several cases of transported Saharan dust in the Caribbean in August 2010,
for example the case on 18 August 2010 that is shown by Burton et al. (2012). For that case, the particle depolarization ratios at 532
and 1064 nm are 0.33
Line plots illustrating the volume and linear particle depolarization ratio profile for the HSRL-2 measurements at 17:12 UT (17.2 UT) on 13 July 2014. The volume depolarization ratio is shown as a thin black line. The error bars on the volume depolarization ratio represent random error (most are small and mostly obscured except 1064 nm). The particle depolarization ratio is shown as a thick colored line. Colored error bars indicate random error (most are small enough to be obscured by the line) while gray error bars indicate systematic uncertainty, estimated as described in the text. Systematic uncertainty is not shown for the volume depolarization ratio; see text for estimate. The vertical resolution of these measurements is 30 m and the horizontal resolution is 10 s for all wavelengths.
Measurement curtain of aerosol backscatter coefficient at 532 nm from the HSRL-2 instrument for a 280 km flight segment over the southwestern US on 8 February 2013, showing locally generated dust in approximately the first kilometer above the surface, as well as very tenuous smoke plumes at higher altitude. This flight segment was part of a transit flight from a field mission in California back to the B200 home base in Virginia. The selected flight segment is approximately 280 km and begins (at the left margin) on the slopes of the Dos Cabezas Mountains of Arizona and ends (at the right margin) at the Franklin Mountains in New Mexico. The ground surface is marked with a white line.
The spectral dependence of the particle depolarization ratio can also be compared to measurements of the Saharan dust particle depolarization ratio reported by Freudenthaler et al. (2009) for the SAMUM I campaign. For the four dates presented in their Fig. 7, the minimum and maximum values for 355 nm from the portable lidar system (POLIS) were 0.21–0.31 (compare 0.25 for NASA HSRL-2 case); at 532 nm they were 0.29–0.33 from the German Aerospace Center (DLR) HSRL (compare 0.30 and 0.33 from NASA HSRL-2 and HSRL-1) and at 1064 nm they were 0.22–0.29 from the DLR HSRL (compare 0.27 and 0.28 from NASA HSRL-2 and HSRL-1). Again, the reported values at 532 nm exceed those at the other wavelengths. Indeed all three of these case studies, Saharan dust close to the source (Freudenthaler et al., 2009), transported Saharan dust observed in the Caribbean by HSRL-1, and transported Saharan dust observed by HSRL-2 in the midwestern US, have similar wavelength dependence of the particle depolarization ratio.
A less typical observation of dust-dominated aerosol was made by the HSRL-2 instrument on 8 February 2013. On a transit flight to Virginia at the conclusion of the DISCOVER-AQ California field campaign, HSRL-2 aboard the B200 made three-wavelength measurements of a locally produced dust layer very close to the source in the US Southwest. Figure 6 shows the lidar curtain of the aerosol backscatter coefficient at 532 nm for a segment of approximately 280 km in Arizona and New Mexico. The highest backscatter values are near the surface and are associated with the dust layer. More tenuous layers are also visible between 3 and 5 km, which are probably smoke. The discussion will focus primarily on the dust layer for this example. Figure 7 shows the particle depolarization ratio at three wavelengths for the same flight segment.
The linear particle depolarization ratio at three wavelengths observed by HSRL-2 on 8 February 2013 in the southwestern US for the flight segment shown in Fig. 6.
The maximum backscatter values occur within 400 m of the ground at about
17:08 UT (17.14 UT), near the Potrillo volcanic fields in New Mexico in the
Chihuahuan Desert. The layer is shallower than the previous case, and the
layer AOD is only about 0.02 at 532 nm, but it is very strongly scattering,
with a 532 nm total aerosol scattering ratio of 3.1. The peak particle
depolarization ratio is 0.24
Figure 8 shows line plots of the volume depolarization ratio and the
particle depolarization ratio profiles, plus error bars. The systematic uncertainties
are generally larger at 355 nm. This error magnification at 355 nm occurs in
both dust-dominated cases because of the spectral dependence of the
scattering and consequent small total aerosol scattering ratio at 355 nm
(
Line plots illustrating the volume and linear particle depolarization ratio profile for the HSRL-2 measurements at 17:08 UT (17.14 UT) on 8 Feb 2013. Error bars and resolutions as described for Fig. 5.
Figure 9 shows the linear particle depolarization ratio at all three wavelengths for the four HSRL-1 and HSRL-2 cases discussed so far. The HSRL-2 observations of transported Saharan aerosol have spectral dependence consistent with the elevated Saharan dust-dominated aerosol reported by Freudenthaler et al. (2009) for the DLR airborne HSRL and ground-based lidar. However, the NASA HSRL-1 and HSRL-2 observations of North American dust at low altitude close to the source appear to fall into a different category. Although all of the observations discussed here from the NASA HSRL-2 and those of Saharan desert aerosol in Africa (Freudenthaler et al., 2009) have particle depolarization ratios at 355 nm that are less than those at 532 nm, there is a large difference at the longest wavelength, with larger 1064 nm particle depolarization ratios for the local dust-dominated cases.
Furthermore, the backscatter Ångström exponents in the Chihuahuan Desert observation on 8 February 2013 and on Pico de Orizaba on 12 March 2006 are much smaller compared to 0.45–0.68 for the cases of transported Saharan dust. These smaller values are an indication of larger particle sizes (Sasano and Browell, 1989) (although it must be noted that the backscatter Ångström exponent is also sensitive to other factors besides particle size, such as relative humidity (Su et al., 2008)). Maring et al. (2003) shows measured size distributions for dust layers over the Canary Islands and Puerto Rico at different stages of transport, and concluded by modeling of these distributions that a combination of Stokes gravitational settling and an offset upward velocity would explain these observations. According to that model, the volume mean diameter decreases only 20 % after 10 days of atmospheric transport, but 80 % of that change occurs within the first 2 days. In other words, the size distributions for transported dust-dominated aerosol are similar whether transported long distances or short distances, but even layers transported short distances probably have already lost the largest particles to settling. This model applies to Saharan dust transport, but it raises the possibility that dust-dominated aerosol size distributions immediately over the source, such as the North American dust cases presented here, will have some proportion of particles significantly larger than those found in the transported layers.
The linear particle depolarization ratio measured by HSRL-2 and HSRL-1 for the four dust cases discussed in the text. Note the spectral dependence (and in particular the 1064 nm channel) is different for the two local dust-dominated aerosol cases compared to the transported Saharan dust-dominated aerosol cases.
The spectral dependence of the particle depolarization ratio is known to be related to the size of the non-spherical particles (Mishchenko and Sassen, 1998). We infer that the difference in depolarization spectral dependence, and in particular the 1064 nm values, are due to larger particles in the observations of windblown dust close to the surface. Ground-based lidar observations by Ansmann et al. (2009) of convective plumes of dust and sand being lifted from the surface in Morocco included extremely large particle depolarization ratios of 0.50–1.0 at 710 nm. This supports the connection between large particle depolarization ratios at the longer wavelengths and large particles sizes. However, the long-wavelength values in the current case study are not nearly as extreme, suggesting perhaps that the particle sizes are not as large.
Theoretical calculations to date have shown that it is difficult to quantitatively predict the spectral dependence of the particle depolarization ratio for dust (Gasteiger et al., 2011; Wiegner et al., 2009; Gasteiger and Freudenthaler, 2014), due in part to the need for parameterizing the shape of the dust aerosols as spheroids or other simplified shapes.
In a theoretical treatment of a particular measurement case, Gasteiger et al. (2011) modeled the particle depolarization ratio at multiple wavelengths using size distributions and refractive indices appropriate for SAMUM measurements, using mixtures of various irregular shapes to represent the dust particles. For their reference distribution, the modeled particle depolarization ratio reflects a spectral dependence with a peak in the middle of the wavelength range. Calculated values at 355, 532, 710, and 1064 nm were 0.275, 0.306, 0.311, 0.298, respectively, consistent with the measurements we report for the Saharan dust-dominated cases from the NASA HSRL-1 and HSRL-2. However, for the dust-dominated cases in the immediate vicinity of North American sources, the measured maximum shifts to longer wavelengths, and there is no longer agreement with the modeled values at 1064 nm.
Gasteiger et al. (2011) do not show results for size distributions with different size particles, but Gasteiger and Freudenthaler (2014) perform theoretical calculations using spheroids for various size parameters (single particles). These calculations show that the first peak in the spectral depolarization ratio shifts to larger wavelengths as particle size increases. This result, based on highly simplified modeling of dust aerosol, should be used only cautiously, but in general supports the notion that the spectral particle depolarization ratio is sensitive to particle size.
Our third three-wavelength case study is an observation of a smoke plume
with large particle depolarization ratios measured during the Colorado
deployment of the DISCOVER-AQ field mission on 17 July 2014 at about 8 km
altitude. At this time, wildfire smoke from fires in the Pacific Northwest
of the United States blanketed much of the region, visible in a composited
Moderate Resolution Imaging Spectroradiometer (MODIS) true color image in Fig. 10. (The smoke situation on that day is also discussed in University
of Maryland Baltimore County's US Air Quality Smog Blog; see
MODIS Aqua true color images of much of North America on 17 July 2014, composited from four granules at 19:45, 19:50, 21:25, and 21:30 UT. The approximate location of the HSRL-2 observations discussed in the text (Denver, Colorado) is marked with a yellow dot. The bright white is clouds and snow cover and the gray is smoke. Several distinct smoke plumes indicate sources in the US Pacific Northwest and in western Canada within the cloud-free area on the western part of the continent. Significant smoke layers from these fires blanket the mid-continent cloud-free areas in the northern portion of the image. The HSRL-2 measurements are close to the southern edge of the extensive smoke field.
Figures 12 and 13 show HSRL-2 measurements of 532 nm backscatter and
three-wavelength linear particle depolarization ratios as time-height
curtains and Fig. 14 illustrates a profile at 19:18 UT (19.3 UT) as a line
plot with random and systematic uncertainty error bars. The pictured flight
segment began near the Boulder Atmospheric Observatory tall tower north of
downtown Denver and proceeded south for about 70 km to Chatfield State Park, then
turned north again on a parallel track for 135 km to Fort Collins. The layer
optical thickness is about 0.05 at 532 nm and the total aerosol scattering
ratio is 2.9. This layer was at high altitude near the aircraft, in the
overlap region, where there is a range dependence of the detected
backscattered light (Hair et al., 2008). While backscattered light from a
distant target is fully imaged in the detector, light from a near-field
target is focused beyond the field stop, resulting in overfilling of the
field stop at small range from the lidar. This loss of signal is range
dependent and prevents the retrieval of aerosol extinction. For this reason,
the layer optical depth given above is an estimate using the backscatter
measurements and an assumed lidar ratio of 70 sr, which is typical for
smoke. Volume depolarization ratio measurements and total aerosol scattering
ratio measurements are ratios of two channels that are equally affected and
therefore have no range-dependent overlap function. For this layer, the
particle depolarization ratio is greatest at 355 nm, about 0.24
View of the smoke plume aloft on 17 July 2013 taken from the B200. Photo credit: Tim Berkoff.
532 nm aerosol backscatter coefficient measurement curtain from HSRL-2 for a portion of a flight on 17 July 2014 in and around Denver, Colorado. Approximately the first third of the pictured curtain is a southbound track between the Boulder Atmospheric Observatory Tall Tower and Chatfield State Park, CO. The remainder of the flight is a northbound leg between Chatfield State Park and Fort Collins. The blank region indicates a tight turn at Chatfield State Park where the lasers were shuttered. Scattered clouds are visible at the top of the boundary layer. Some of these have off-scale backscatter values (tan color) and some are thick enough to cause significant attenuation of the beam; beneath these, data are blanked out due to low signal. The white line indicates underlying terrain. The smoke layer at approximately 8 km is discussed in the text.
The linear particle depolarization ratio measurement curtains for the flight segment shown in Fig. 12.
The pattern of larger particle depolarization ratios at 532 nm compared to
1064 nm has regularly been observed for smoke with the HSRL-1 instrument;
indeed the HSRL-1 aerosol classification methodology (Burton et al., 2012)
takes advantage of this spectral dependence. One such example is the aged
southwest Canadian smoke plume observed on the eastern seaboard of the US
on 2 August 2007 that was shown by Burton et al. (2012). For that prior
case, the particle depolarization ratios were 0.07
The observations by HSRL-2 in Figs. 12 and 13 are, to our knowledge, the first reported three-wavelength measurements of linear particle depolarization ratio from pure smoke. Note that while the 532 nm particle depolarization ratio for the smoke case is only about 25–30 % of the value for pure dust, the large particle depolarization ratio at 355 nm for the smoke layer is quite comparable to the 355 nm value for dust.
Observed linear particle depolarization ratios for smoke are quite variable. Frequently, the 532 nm particle depolarization ratio is observed to be only a few percent at most and often discounted as negligibly close to zero (e.g., Mattis et al., 2003; Müller et al., 2005). For example, small values of about 2–3 % at 532 nm were observed by HSRL-1 in smoke plumes in Mexico City during the Megacity Initiative: Local and Global Research Observations (MILAGRO) field mission described by de Foy et al. (2011). In the published images of HSRL-1 measurements, the smoke plumes are obvious as local minima in the particle depolarization ratio compared to the relatively high ambient values, which are due to regional dusty background conditions. However, higher values of the particle depolarization ratio at 532 nm of 0.05–0.11 have sometimes been observed in aged smoke layers, e.g., 0.07 observed by HSRL-1 on 2 August 2007 noted above (Burton et al., 2012), 0.06–0.11 for transported Canadian smoke reported in Lindenberg, Germany, during the Lindenberg Aerosol Characterization Experiment (LACE) 1998 by Fiebig et al. (2002), 0.06 to 0.08 for 3–4 day old smoke from North America observed over Germany in 2011 by Dahlkötter et al. (2014), 0.06 for transported Siberian smoke observed in Tokyo in 2003 by Murayama et al. (2004), and 0.05 for Alaskan forest fire smoke observed by Sassen and Khvorostyanov (2008) in 2004. Sugimoto et al. (2010) discuss a case in which much higher 532 nm particle depolarization ratios were observed for smoke from a Mongolian forest fire transported to Japan in 2007. The particle depolarization ratios measured for this smoke were 0.12, 0.14, and 0.15 for layers at two different altitudes observed at Nagasaki and Tsukuba. Nisantzi et al. (2014) observe values of 0.09 to 0.18 at 532 nm for aerosol from Turkish fires observed in Cyprus after 1 to 4 days of transport.
Line plots illustrating the volume and linear particle depolarization ratio profile for the HSRL-2 measurements at 19:18 UT (19.3 UT) on 17 July 2014. Error bars and resolutions as described for Fig. 5.
The causes of depolarization by smoke are not well understood. Two possible explanations are frequently cited in literature: lifting and entrainment of surface soil into the smoke plume and asymmetry of smoke particles themselves.
For example, the smoke observed by Sugimoto et al. (2010) was associated with pyrocumulonimbus and therefore it is inferred that strong convection lifted soil particles from the surface into the smoke plume, explaining the unusually large particle depolarization ratios. Lifting of soil particles is also cited as a possible explanation of the more moderate but still non-negligible particle depolarization ratios reported by Fiebig et al. (2002), since chemical composition analysis of this plume reveals the presence of aluminosilicates and iron oxides/hydroxides. Nisantzi et al. (2014) assume that the depolarization is due to fine-mode dust and infer the mass fraction of dust mixed in the smoke plumes using lab measurements of fine and coarse dust by Sakai et al. (2010). However, this explanation is not sufficient in every case. Murayama et al. (2004) discount soil lifting for their observations of depolarizing smoke, since no signature of mineral dust is found in a chemical analysis of this plume. Instead, they cite non-sphericity of smoke particle aggregates as the probable cause. Martins et al. (1998) discuss the non-sphericity of smoke particles observed by scanning electron microscope images and an electro-optical aerosol asymmetry analyzer for a variety of smoke types during the Smoke, Clouds, and Radiation – Brazil (SCAR-B) project in 1995. They concluded that most of the non-spherical particles in the observed smoke were chain aggregates of small black carbon particles, and that the non-sphericity tends to increase with the black carbon ratio. Young smoke (< 1 h) is composed of open clusters of high non-sphericity while aged smoke is composed of tighter clusters with lesser non-sphericity. They also point out that flaming fires (high combustion efficiency) tend to produce more non-spherical particles than smoldering fires.
Referring back to the theoretical calculations of spectral depolarization for spheroids discussed in Sect. 3.3., the larger particle depolarization ratio at 355 nm compared to longer wavelengths may indicate a smaller size for the non-spherical particles than the dust cases, although of course these results may be only qualitatively applicable to more general particle shapes.
Theoretical calculations of the linear particle depolarization ratio by aggregates of soot are given by, e.g., Sorensen (2001), Bescond et al. (2013), and Kahnert et al. (2012). Calculations for bare carbon aggregates (Sorensen, 2001; Bescond et al., 2013) tend to produce small values of the particle depolarization ratio, much smaller than what was measured in the HSRL-2 case study. However, Kahnert et al. (2012) modeled the scattering properties of a more realistic particle morphology, light absorbing carbon (LAC) aggregates embedded in a sulfate shell, and obtain larger values. They use the discrete dipole method to calculate the depolarization ratio of the aggregate particle in the backscatter direction at 304.0, 533.2, and 1010.1 nm. They show that the particle depolarization ratio generally increases with aggregate particle radius (defined as volume-equivalent radius) and with the volume fraction of LAC in the aggregate. The values also increase with decreasing wavelength for aggregate volume-equivalent radii of 400 nm and smaller; conversely, for 500 nm particles, the particle depolarization ratio peaks at the middle wavelength, 533.2 nm. The maximum calculated particle depolarization ratios for 7 % LAC fraction by volume is 0.08–0.11 for 500 nm particles at 533.2 nm. This is comparable to the 532 nm measurement on 17 July 2014; however, the calculated value at 304.0 nm for the same size and LAC volume fraction is 0.05–0.07, much smaller than the measured value (at 355 nm) of 0.24. The largest values calculated for the 304.0 nm wavelength are about 0.12–0.21, occurring for the case of 400 nm particle volume-equivalent radius and 20 % LAC volume fraction. The full set of theoretical calculations of the particle depolarization ratio for 20 % LAC volume fraction are replotted in Fig. 15 for all three wavelengths to highlight the wavelength dependence. Figure 15 also indicates the HSRL-2 observed particle depolarization ratio in the 17 July smoke plume (at 355, 532, and 1064 nm). The calculated particle depolarization ratios are roughly comparable in magnitude to the HSRL-2 measurements for volume-equivalent radii in the 400–500 nm range, but the wavelength dependence matches better for smaller particle sizes. A LAC volume fraction of 20 % is quite high and may be unrealistic for this smoke layer and the modeled single scattering albedos for a 20 % LAC volume fraction, shown by Kahnert et al. (2012), are quite low (below 0.7 at 533.2 nm), indicating exceptionally absorbing particles; therefore, this model is probably not an exact match for the observation in this case. Yet, it is encouraging that an estimate of the particle depolarization ratio of the right magnitude can be made by modeling coated soot aggregates. The model results were for a constant fractal dimension of 2.6, structural prefactor of 1.2, and a monomer radius of 25 nm, values chosen to be consistent with the findings for soot aerosol in Mexico City (Adachi and Buseck, 2008). In the HSRL-2 case study, there could be a different fractal dimension, different size monomer component, different coating, or a different fraction of soot per aggregate. In addition, the spectral dependence of the refractive index is not well-known, and this will have a significant effect on the spectral dependence of the particle depolarization ratio. While the current state of knowledge is not sufficient to perform a retrieval of particle size using the depolarization measurements alone, it is certainly worth noting that the particle depolarization ratio at three wavelengths is sensitive to and contains some information about the particle size of smoke particles, information that may play a role in future microphysical retrievals.
The linear particle depolarization ratio at three wavelengths for soot aggregates embedded in a sulfate shell reproduced from Kahnert et al. (2012), for 20 % LAC volume fraction. Dots indicate five realizations with randomly generated geometries, per aggregate volume-equivalent particle radius. The colored lines connect the averages of the five for each wavelength. The legend shows the aggregate volume-equivalent particle radii at which the calculation was performed. The thick black line indicates the particle depolarization ratios measured by airborne HSRL-2 within a smoke plume observed on 17 July 2014 at 355, 532, and 1064 nm.
Illustrates the systematic uncertainty in the linear particle
depolarization ratio propagated from the systematic uncertainties in total
aerosol scattering ratio, linear volume depolarization ratio, and linear
molecular depolarization ratio. Benchmark values of
We have presented three case studies of depolarizing aerosol observed at three wavelengths (355, 532 and 1064 nm) by the NASA airborne HSRL-2 instrument. These three aerosol layers, two dust-dominated layers and a smoke layer, each have a different spectral dependence of the linear particle depolarization ratio, but in each case, the 532 and 1064 nm values agree well with prior analogs in the long record of observations by the predecessor instrument, HSRL-1, and with comparable measurements in literature. The first case, transported Saharan desert aerosol, has a peak in the spectral dependence of the particle depolarization ratio at 532 nm. This is in accordance with prior measurements of Saharan desert aerosol aloft both close to the source and transported to the Caribbean Sea. The second case, also a dust-dominated measurement, but near the surface and very close to the source, has a spectral dependence increasing monotonically with wavelength, differing from the previous case primarily at the longest wavelength, 1064 nm. We infer the cause of this difference to be a greater fraction of very large particles due to its proximity to the source region; we believe that the largest particles have settled out of the observed Saharan layers but not the locally produced North American dust plumes in this case and a prior HSRL-1 case. Our third case study is of an elevated, transported smoke layer and has a spectral depolarization ratio decreasing monotonically with wavelength. Again we infer that the difference in spectral dependence is due to the size of the non-spherical particles, and, specifically, that the depolarization is probably due to smoke aerosols and may be explained by soot aggregates.
Microphysical retrievals (Müller et al., 2014) were not available for these HSRL-2 measurement cases, because the current state of these retrievals is limited to spherical particles. However, as suggested by Gasteiger and Freudenthaler (2014) for dust and ash, these observations suggest the possibility that the particle depolarization ratio measurements may aid in retrievals of particle size of non-spherical dust and smoke particles in the future.
More immediately, since the upcoming EarthCARE satellite mission will include a lidar instrument that measures the particle depolarization ratio and lidar ratio at 355 nm only, it is valuable to have measurements of the spectral dependence of depolarization ratio for depolarizing aerosol types. These data will help to build the basis for comparing observations from EarthCARE to existing measurements at 532 nm from the CALIPSO satellite. Studying such correspondence is particularly motivated by the desire to identify different aerosol types observed by the EarthCARE satellite. The particle depolarization ratio is hoped to be particularly useful for distinguishing dust and ash from smoke and other aerosol types (Groß et al., 2014; Illingworth et al., 2015), as it already is for CALIPSO (Omar et al., 2009).
However, as illustrated by the case studies presented here, there is not a single consistent spectral dependence of the particle depolarization ratio. On the positive side (from the perspective of corresponding CALIPSO and EarthCARE measurements), for aerosols dominated by dust the 355 and 532 nm particle depolarization ratios appear to be fairly consistent even for different particle sizes and may be relatively easily converted. Variation in the 532 and 355 nm particle depolarization ratios for dusty aerosols has been primarily linked to the fraction of dust particles in a sample (Sugimoto et al., 2003); therefore, there is no reason to think that inferences of dust mixing ratios (e.g., Sugimoto and Lee, 2006; Tesche et al., 2009a; Nishizawa et al., 2011; Burton et al., 2014) may not be done with 355 nm measurements. However, in the case of dust-dominated aerosol, the 355 nm signal consistently is significantly both smaller and more difficult to measure accurately than the 532 nm signal, and therefore the signature of dust may be harder to detect from space at 355 nm than at 532 nm for dilute dust mixtures.
On the other hand, the third case study presented here showed that the smoke particle depolarization ratio can be significantly larger at 355 nm than at 532 nm, and in fact the particle depolarization ratio at 355 nm for this smoke case was quite comparable to the dust-dominated cases. If this is not an isolated case, and this signature proves typical for some subsets of smoke aerosol in particular conditions, the EarthCARE satellite may observe significant particle depolarization in some types of smoke as well as in dust-dominated aerosol. If this is the case, global observations of smoke depolarization will present an exciting opportunity for improving our understanding of the optical properties of smoke and how they change with age and processing; however, it will also present a challenge. That is, a significant particle depolarization ratio signature at the single wavelength of 532 nm has been sufficient for distinguishing dust-dominated aerosol from smoke aerosol, but at 355 nm this signature by itself is more ambiguous, if the smoke case presented here is not an isolated case. EarthCARE will also measure the lidar ratio at 355 nm; this is related to absorption but has significant variability for smoke (Groß et al., 2014). EarthCARE will not have backscatter or extinction measurements at a second wavelength to give an indicator of particle size. Therefore, for any cases where the particle depolarization ratio is ambiguous, smoke and dust may not be easily separable.
In Sect. 2.2 we provided systematic uncertainties in the linear volume
depolarization ratio of the larger of 4.7 % (relative) or 0.001 (absolute)
in the 355 nm channel, the larger of 5 % (relative) or 0.007 (absolute) in
the 532 nm channel, and the larger of 2.6 % (relative) or 0.007 (absolute)
in the 1064 nm channel. For
For the linear volume depolarization ratio, potential sources of systematic error include an error in the polarization gain ratio calibration or cross talk between the co-polarized and cross-polarized signals. The polarization gain ratio calibration generally occurs once or twice per flight as described above in Sect. 2.1. Since gain ratios can potentially change during a flight, due to temperature changes for example, our best estimate of uncertainty in the gain ratio during a flight is obtained by examining the change in the gain ratio between successive calibrations in the same flight. Conservatively choosing the mean difference plus 2 standard deviations (calculated for all flights with at least two calibrations per flight in the most recent field campaign) as a realistic limit on the probable polarization gain ratio systematic error yields 4.7 % uncertainty for the 355 nm channel, 5.0 % for the 532 nm channel, and 2.6 % for the 1064 nm channel. The relative systematic uncertainty from the polarization gain ratio propagates directly to the volume depolarization ratio, since the volume depolarization ratio is linearly related to the polarization gain ratio.
Residual cross talk is known to occur in polarization lidars, and must be
carefully characterized and eliminated as much as possible. A well-known
potential source of cross talk occurs in the reflected channel from a
polarization beam splitter. Therefore, this system has been designed with
extra polarization beam splitters to eliminate that potential concern, as
described in Sect. 2.1 and illustrated in Fig. 2. Clear-air studies have
found a small residual cross talk, which appears as a value of the clear-air volume depolarization ratio that exceeds the theoretical (molecular-only) value. As described in Sect. 2.1, the narrow bandwidths in the
system completely eliminate the rotational Raman scattering sidebands, and
so the molecular depolarization ratio is temperature independent and is
calculated to be 0.0036 using N
We start with the polarization Stokes vector (Born and Wolf, 1999)
We believe that the polarization angle error is much smaller than the
inferred angle of 5.8
Note that not only the volume depolarization ratio measurement itself but
also the polarization gain ratio calibration depend on the correct alignment
of the calibration wave plates in Fig. 1. The polarization gain ratio
assessment depends on a polarization alignment of 45
The calculated particle depolarization ratio,
The 355 and 1064 nm channels are somewhat more complicated, because it is
not possible to calibrate them directly in the same way as 532 nm. The
iodine filter for the 532 nm HSRL channel allows for essentially complete
separation of the aerosol signal from the total (aerosol plus molecular)
signal, but this is not the case for the interferometer used at 355 nm, and
the 1064 nm channel has only one total channel with no separation. So for
these channels, the calibration is transferred from 532 nm in a cloud-free
region in the free troposphere, as described by Hair et al. (2008). In the
calibration transfer region, we do not assume that there is no aerosol, but
do look for regions where the aerosol backscatter ratio is small and can be
inferred from the value at 532 nm assuming a constant backscatter
Ångström exponent. By using a range of reasonable backscatter
Ångström exponents, we conservatively estimate an uncertainty of
3 % in total aerosol scattering ratio for the 355 nm channel. The 1064 nm
aerosol backscatter ratio is also affected by the assumption of the lidar
ratio to use for separating the aerosol and molecular part; this sensitivity
is relatively small for backscatter at 1064 nm, compared to shorter
wavelengths or compared to the sensitivity of extinction. Taking these
sources into account, we conservatively use 20 % as the uncertainty in
total aerosol scattering ratio,
For the 355 nm channel, the system implements an interferometer to spectrally separate the aerosol and molecular scattering components. The ratio of the aerosol signal in the aerosol-dominated channel to the aerosol signal in the molecular-dominated channel is referred to as the contrast ratio, which needs to be determined to accurately derive the total aerosol scattering ratio. For the HSRL-2 system, fairly high contrast ratios of 15–20 are routinely achieved. Our estimate of the error in the contrast ratio definition is usually a few percent but can be up to 20 %. A 20 % error in the contrast ratio for the smoke case presented here would produce an error in the total aerosol scattering ratio of less than 4 %. Adding the contrast ratio uncertainty, 4 %, and the calibration transfer uncertainty, 3 %, in quadrature yields an uncertainty of 5 % for the 355 nm total aerosol scattering ratio.
The uncertainties given above are intended to be an upper bound on the
probable systematic errors. The systematic errors on the three quantities,
The partial derivatives are calculated easily from Eq. (2), which relates the
particle depolarization ratio to the factors
These factors vary with total aerosol scattering ratio and volume
depolarization ratio but do not depend on the systematic uncertainties. To
illustrate the behavior of the particle depolarization ratio systematic
uncertainty, Table 2 gives the value of the particle depolarization ratio and
its propagated systematic uncertainty (as a percent error) for benchmark
values of the total aerosol scattering ratio and the volume and molecular
depolarization ratios, plus their estimated systematic uncertainties. It
also gives the propagation factors,
Funding for this research came from the NASA HQ Science Mission Directorate
Radiation Sciences Program and the DISCOVER-AQ project. M. Kahnert
acknowledges funding by the Swedish Research Council (