1. Introduction
The primary objective of the National Aeronautics and Space Administration (NASA)’s Micropulse Lidar Network (MPLNET; Welton et al. 2001), a member of the Global Atmospheric Watch (GAW) Aerosol Lidar Observation Network (GALION; Hoff and Boesenberg 2008), is the depiction of cloud and aerosol particle distributions, utilizing a global and federated network of ground-based lidar instruments. Enhancing the project, efforts are presently ongoing toward the integration of ultraviolet (UV; 355 nm), eye-safe, and autonomous lidar produced by Leosphere (Lolli et al. 2011) to complement the historical array of visible (VIS; 523, 527, and 532 nm) and eye-safe micropulse lidar instruments (Spinhirne 1993; Spinhirne et al. 1995; Campbell et al. 2002). The integration of these instruments would facilitate the emergence of a denser global network. When possible, collocated UV and VIS lidar profiling would further provide multispectral information on aerosol particle and cloud scattering.
Covering a large fraction of the planet surface at any given time, stratocumulus and stratus clouds play an important role in planetary radiation budget (e.g., Bennartz et al. 2013). These clouds are important because they redistribute the column water liquid content due to the transformation of water falling from clouds into vapor when precipitating (Harrison et al. 1990; O’Connor et al. 2005). Numerical models require an accurate representation of the microphysical properties of these boundary layer clouds in order to explicitly resolve their impact on climate (Slingo and Slingo 1991; O’Connor et al. 2005). The purpose of this paper is to describe and evaluate a method for measuring light rain and drizzle drop median equivolumetric size diameter D0 from collocated MPLNET UV and VIS lidar profiles by taking advantage of the differential backscatter during precipitation at the two wavelengths, respectively. Raindrop size is of particular importance for controlling the rate at which they deposit through the atmosphere and how quickly they evaporate below cloud base (Westbrook et al. 2010).
Though atmospheric lidar profiling is attenuation limited under many common aerosol- and cloud-observing conditions (Sassen and Cho 1992), a scenario exacerbated by liquid water droplet presence given the relatively high scattering cross sections and droplet concentrations (e.g., Sassen 1978), such data have been effectively used for investigating melting-layer thermodynamics and the scattering properties of precipitating low-level liquid water clouds (e.g., Sassen et al. 2005, 2007; Campbell and Shiobara 2008; Di Girolamo et al. 2012). Our goal is to transpose a two-color lidar-based droplet measurement technique designed for infrared bands (IR; Westbrook et al. 2010) to the UV–VIS region and to evaluate performance so as to extend, when possible, MPLNET (http://mplnet.gsfc.nasa.gov) capabilities for studying stratocumulus and stratus cloud precipitation size distributions.
2. Methods






Lidar Qbkg as defined in Eq. (3), (top) derived from Mie theory for spherical water drops in 1-μm droplet size intervals and (bottom) Qbkg computed over a size distribution defined in Eq. (1) with μ = 2.
Citation: Journal of Atmospheric and Oceanic Technology 30, 12; 10.1175/JTECH-D-13-00062.1
Corresponding with Fig. 1 Qext as defined in Eq. (4) (top) derived from Mie theory for spherical water drops in 1-μm droplet size intervals and (bottom) Qext computed over a size distribution defined in Eq. (1) with μ = 2.
Citation: Journal of Atmospheric and Oceanic Technology 30, 12; 10.1175/JTECH-D-13-00062.1
The 355-/527-nm lidar color ratio, defined in Eq. (5), solved as a function of D0.
Citation: Journal of Atmospheric and Oceanic Technology 30, 12; 10.1175/JTECH-D-13-00062.1




CR retrievals and relative measurement uncertainty
To retrieve D0 from Eqs. (5) and (11), the transmission due to aerosol particles and rain and molecular and aerosol backscattering at two different wavelengths must be evaluated. In Eq. (11), for D0 ≥ 50 μm, attenuation due to precipitation or





The uncertainty of the measured range-corrected signal
The error in aerosol transmission
The effects of the aerosol contamination to backscattering are discussed further below.
3. Performance evaluation
Light rain showers were observed late in the afternoon (starting at 2100 UTC) on 9 May 2012 at the Goddard Space Flight Center in Greenbelt, Maryland [38.99°N, 76.83°W; 84 m, all heights above mean sea level (MSL)]. Shown in Fig. 4 are the logarithms of the normalized attenuated VIS and UV MPLNET lidar backscatter [Eq. (13)] composites for 1340–2359 UTC, respectively. Both instruments recorded data at 60-s temporal and 75-m vertical resolutions. The cloud base is shown lowering in both profiles from 5.0 to 3.0 km, with light rain first appearing near 2100 UTC (Fig. 4). The melting layer remains constant in height through the period; this is also confirmed by Dulles International Airport (38.9°N, 77.45°W; 95 m) radiosonde isotherm at 0°C, at about 3.1 km. A lidar dark band (e.g., Sassen et al. 2005; Di Girolamo et al. 2012) appears in both profile composites near 2.9 km, indicating a cold rain process as the precipitation originates from melting ice. Signal normalization was conducted above 3.0 km from 1900 to 2000 UTC (60 1-min-averaged profiles).
On 9 May 2012 at NASA GSFC from 0 to 5 km MSL, MPLNET (top) 527- and (bottom) 355-nm logarithmic normalized attenuated lidar backscatter coefficient [log(Mm−1 sr−1)] from 1340 to 2359 UTC.
Citation: Journal of Atmospheric and Oceanic Technology 30, 12; 10.1175/JTECH-D-13-00062.1
After 2100 UTC, precipitating fall streaks are distinct in both signal composites. From calculations described in section 2, backscatter signal profiled from precipitation at 355 nm is greater than at 527 nm. Composite depictions of CR and D0 from 2024 to 2248 UTC are shown in Fig. 5. In the cloud layer, CR is roughly 0 dB, meaning that for cloud droplets (<30 μm in diameter) attenuation due to scattering is effectively the same at both wavelengths. The derived CR from clear sky (molecular scattering only), being the attenuated backscattering profiles normalized to the attenuated molecular backscattering, is about 0 dB.
On 9 May 2012 at NASA GSFC from 0 to 5 km MSL, (top) MPLNET 355-/527-nm backscatter color ratio (dB) and (bottom) retrieved median droplet diameter (μm) from 2024 to 2248 UTC.
Citation: Journal of Atmospheric and Oceanic Technology 30, 12; 10.1175/JTECH-D-13-00062.1
The value of D0 is retrieved through bivariate analysis. The range of retrieved D0 values shown (Fig. 5, bottom) is consistent with light rain and drizzle, sized between 100 and 500 μm (e.g., Frisch et al. 1995).
Figure 6a shows how the uncertainty relative to the color ratio (profile taken at 2143 UTC), calculated using Eqs. (14)–(18), translates into uncertainty in D0 retrieval (Fig. 6b). The uncertainty in the CR is mostly dominated by error in lidar measurement and degrades as the signal-to-noise ratio. The error in D0 is instead related to the unknown drop size spectrum (parameter μ) and directly to the shape of the theoretical CR [Eq. (5); Fig. 3] that flattens as the raindrop diameter increases: the uncertainty is then bigger for larger drops (near 50% for 600 μm, Fig. 6b).
(a) Measured CR profile at 2143 UTC with relative uncertainty. (b) Retrieved D0 from CR profile and relative uncertainty. As the theoretical CR flattens, increasing the raindrop diameter, the uncertainty is bigger for larger drops, near 50% for 600 μm.
Citation: Journal of Atmospheric and Oceanic Technology 30, 12; 10.1175/JTECH-D-13-00062.1
Errors because of multiple scattering and aerosol returns
Though the 527-nm MPL instrument features a very narrow field of view, the UV lidar design includes one order of magnitude larger, and thus multiple scattering effects must be taken into account. O’Connor et al. (2004) evaluate the multiple scattering effect for a UV lidar with similar optical characteristics in a liquid water cloud consisting of relatively small droplets. Considering some slight differences between the two UV lidars, multiple scattering was found to reach a maximum 15% at roughly 1 km under heavy drizzle conditions. This is an extreme case, and being that the attenuated backscattering coefficients in rain are lower than 30 Mm−1 sr−1, a value of 10% for the multiple scattering is more likely. This static value has been applied for solving Eq. (12).
At 355 nm, backscatter from aerosols is at least one order of magnitude smaller than that from precipitation. At 527 nm, the smaller backscattered power from precipitation may be comparable to the returns from the aerosol particles, especially in the boundary layer. This will lead to D0 being underestimated. We estimate the aerosol return by looking at precipitation-free areas of the boundary layer in Fig. 4 (top). The ratio of precipitation to aerosol backscattered power is 5:1. For a measured color ratio of 4 dB, the retrieved median equivolumetric drop size is D0 = 140 μm. In this situation, the aerosol backscattered power accounts for one-fifth of the signal and the true ratio for the rain signal alone would be 4.6 dB, corresponding to a true value for the raindrop size of D0 = 162 μm. This value is taken into account to retrieve the correct D0 in the aerosol-contaminated pixels mostly in the boundary layer.
When considering each of the aforementioned sources of algorithm uncertainty, in tandem with that due to signal processing, solutions for D0 correspond with a maximum total uncertainty of roughly 50%. Though this is a considerable amount of relative error, it still implies that there is some reasonable measure of skill exhibited by the technique. Recall that, on most rainy days, continuous lidar profiling measurements are almost immediately discarded from the negative impact of source pulse attenuation. Thus, the potential use of MPLNET measurements for estimating D0, and with the thought that some uncertainty could be reduced with continuing instrument improvements and more stable calibrations, is encouraging.
4. Conclusions
In this paper, collocated ultraviolet and visible (355 and 527 nm, respectively) lidar measurements are used to estimate the precipitation drop median equivolumetric diameter D0 over a 3-h period of light rain showers profiled at the NASA Micropulse Lidar Network (MPLNET)-observing site at the Goddard Space Flight Center in Greenbelt, Maryland, on 9 May 2012. Scattering calculations are performed for spherical liquid drops from light rain samples using Mie theory to link D0 to the ratio of backscatter signals at the two wavelengths. These calculations are then applied to retrieve D0 within precipitating fall streaks, such that signal attenuation effects are not limiting, including the parent cloud layers to the extent of their profiling depth.
We estimate that total error in solutions for D0 is at about 50% for larger raindrops. The primary sources of uncertainty in the retrieval are attributable to the lidar signal measurement and the unknown shape of raindrop size distributions. Molecular and aerosol scattering contributions and relative transmission losses due to the various atmospheric constituents were evaluated to derive D0 from the observed color ratio profiles. Multiple scattering is important for the UV lidar, but it exhibits lower overall uncertainty with respect to other error sources identified.
Westbrook et al. (2010) describe the basis for the methodology investigated here. We find that transposing their technique from the near infrared to the UV–VIS region increases retrieval uncertainties; that is, the two wavelengths (355 and 527 nm) are relatively close, and the color ratio curve shapes become flat, especially for large drops (>500 μm), which increase relative error.
Despite this, the demonstration of this method shows that the technique exhibits reasonable skill. With knowledge of D0, higher moments of the drizzle droplet distribution, including liquid water content and rainfall rate, can be further estimated, which will add value to the retrieval overall. At that stage, a value-added MPLNET product reporting these parameters for light rain events can be considered, with the goal of an improved understanding of the microphysics and physical scattering properties of precipitation at near-visible wavelengths from global ground-based measurements.
Acknowledgments
The Micropulse Lidar Network is coordinated through the support of the NASA Radiation Sciences Program. The authors gratefully acknowledge the University of Wyoming’s Department of Atmospheric Sciences for access to its radiosonde visualization tools. Author JRC acknowledges the support of NASA Interagency Agreement NNG12HG05I on behalf of the Micropulse Lidar Network.
REFERENCES
Beard, K. V., 1976: Terminal velocity and shape of cloud and precipitation drops aloft. J. Atmos. Sci., 33, 851–864.
Bennartz, R., and Coauthors, 2013: July 2012 Greenland melt extent enhanced by low-level liquid clouds. Nature, 496, 83–86, doi:10.1038/nature12002.
Bohren, C. F., and Huffman D. R. , 1983: Absorption and Scattering of Light by Small Particles. John Wiley and Sons, 530 pp.
Campbell, J. R., and Shiobara M. , 2008: Glaciation of a mixed-phase boundary layer cloud at a coastal arctic site as depicted in continuous lidar measurements. Polar Sci., 2, 121–127.
Campbell, J. R., Hlavka D. L. , Welton E. J. , Flynn C. J. , Turner D. D. , Spinhirne J. D. , Scott V. S. , and Hwang I. H. , 2002: Full-time, eye-safe cloud and aerosol lidar observation at Atmospheric Radiation Measurement Program sites: Instrument and data processing. J. Atmos. Oceanic Technol., 19, 431–442.
Campbell, J. R., Welton E. J. , Spinhirne J. D. , Ji Q. , Tsay S.-C. , Piketh S. J. , Barenbrug M. , and Holben B. N. , 2003: Micropulse lidar observations of tropospheric aerosols over northeastern South Africa during the ARREX and SAFARI 2000 dry season experiments. J. Geophys. Res., 108, 8497, doi:10.1029/2002JD002563.
Campbell, J. R., Sassen K. , and Welton E. J. , 2008: Elevated cloud and aerosol layer retrievals from micropulse lidar signal profiles. J. Atmos. Oceanic Technol., 25, 685–700.
COESA, 1976: U.S. Standard Atmosphere, 1976. NOAA, 227 pp.
Di Girolamo, P., Summa D. , Cacciani M. , Norton E. G. , Peters G. , and Dufournet Y. , 2012: Lidar and radar measurements of the melting layer: Observations of dark and bright band phenomena. Atmos. Chem. Phys., 12, 4143–4157.
Fernald, F. G., 1984: Analysis of atmospheric lidar observations: Some comments. Appl. Opt., 23, 652–653.
Frisch, A. S., Fairall C. W. , and Snider J. B. , 1995: Measurement of stratus cloud and drizzle parameters in ASTEX with a Kα-band Doppler radar and a microwave radiometer. J. Atmos. Sci., 52, 2788–2799.
Harrison, E. F., Minnis P. , Barkstrom B. R. , Ramanathan V. , Cess R. D. , and Gibson G. G. , 1990: Seasonal variation of cloud radiative forcing derived from the Earth Radiation Budget Experiment. J. Geophys. Res., 95 (D11), 18 687–18 704.
Hoff, R., and Boesenberg J. , 2008: GAW Aerosol Lidar Observations Network (GALION). Preprints, Symposium on Recent Developments in Atmospheric Applications of Radar and Lidar, New Orleans, LA, Amer. Meteor. Soc., 5.2. [Available online at https://ams.confex.com/ams/88Annual/techprogram/paper_131629.htm.]
Lolli, S., Sauvage L. , Stachlewska I. , and Coulter R. , 2008: Assessment of the EZ LIDAR and micro pulse lidar (MPL) performances at ARM Southern Great Plains (SGP) Central Facility for the measurement of clouds and aerosols. Geophysical Research Abstracts, Vol. 10, Abstract EGU2008–A-11091. [Available online at http://www.cosis.net/abstracts/EGU2008/11091/EGU2008-A-11091.pdf.]
Lolli, S., Sauvage L. , Loaec S. , and Lardier M. , 2011: EZ lidar: A new compact autonomous eye-safe scanning aerosol lidar for extinction measurements and PBL height detection. Validation of the performances against other instruments and intercomparison campaigns. Opt. Pura Apl., 44, 33–41.
Lynch, D. K., and Livingston W. , 2001: Color and Light in Nature. 2nd ed. Cambridge University Press, 277 pp.
O’Connor, E. J., Illingworth A. J. , and Hogan R. J. , 2004: A technique for autocalibration of cloud lidar. J. Atmos. Oceanic Technol., 21, 777–786.
O’Connor, E. J., Hogan R. J. , and Illingsworth A. J. , 2005: Retrieving stratocumulus drizzle parameters using Doppler radar and lidar. J. Appl. Meteor., 44, 14–27.
Russell, P., Swissler T. , and McCormick M. , 1979: Methodology for error analysis and simulation of lidar aerosol measurements. Appl. Opt., 18, 3783–3797.
Sassen, K., 1978: Backscattering cross sections for hydrometeors: Measurements at 6328 Å. Appl. Opt., 17, 804–806.
Sassen, K., and Cho B. S. , 1992: Subvisual–thin cirrus lidar dataset for satellite verification and climatological research. J. Appl. Meteor., 31, 1275–1285.
Sassen, K., Campbell J. R. , Zhu J. , Kollias P. , Shupe M. , and Williams C. , 2005: Lidar and triple-wavelength Doppler radar measurements of the melting layer: A revised model for dark- and brightband phenomena. J. Appl. Meteor., 44, 301–312.
Sassen, K., Matrosov S. , and Campbell J. , 2007: CloudSat spaceborne 94 GHz radar bright bands in the melting layer: An attenuation-driven upside-down lidar analog. Geophys. Res. Lett., 34, L16818, doi:10.1029/2007GL030291.
Slingo, A., and Slingo J. M. , 1991: Response of the National Center for Atmospheric Research Community Climate Model to improvements in the representation of clouds. J. Geophys. Res., 96 (D9), 15 341–15 357.
Spinhirne, J. D., 1993: Micro pulse lidar. IEEE Trans. Geosci. Remote Sens., 31, 48–55.
Spinhirne, J. D., Rall J. A. R. , and Scott V. S. , 1995: Compact eye safe lidar systems. Rev. Laser Eng., 23, 112–118.
Welton, E. J., and Campbell J. R. , 2002: Micropulse lidar signal uncertainties. J. Atmos. Oceanic Technol., 19, 2089–2094.
Welton, E. J., Campbell J. R. , Spinhrine J. D. , and Scott III V. S. , 2001: Global monitoring of clouds and aerosols using a network of micropulse lidar systems. Lidar Remote Sensing for Industry and Environment Monitoring, U. N. Singh, T. Itabe, N. Sugimoto, Eds., International Society for Optical Engineering (SPIE Proceedings, Vol. 4153), 151–158.
Welton, E. J., and Coauthors, 2002: Measurements of aerosol vertical profiles and optical properties during INDOEX 1999 using micropulse lidars. J. Geophys. Res., 107, 8019, doi:10.1029/2000JD000038.
Westbrook, C. D., Hogan R. J. , O’Connor E. J. , and Illingworth A. J. , 2010: Estimating drizzle drop size and precipitation rate using two-colour lidar measurements. Atmos. Meas. Tech., 3, 671–681, doi:10.5194/amt-3-671-2010.
Willis, P. T., 1984: Functional fits to some observed drop size distributions and parameterization of rain. J. Atmos. Sci., 41, 1648–1661.
Wolf, S., and Voshchinnikov N. V. , 2004: Mie scattering by ensembles of particles with very large size parameters. Comput. Phys. Commun., 162, 113–123.