Determining the contribution of volcanic ash and boundary layer aerosol in backscatter lidar returns: a three‐component atmosphere approach

A solution of the lidar equation is discussed, that permits combining backscatter and depolarization measurements to quantitatively distinguish two different aerosol types with different depolarization properties. The method has been successfully applied to simultaneous observations of volcanic ash and boundary layer aerosol obtained in Exeter, United Kingdom, on 16 and 18 April 2010, permitting the contribution of the two aerosols to be quantified separately. First a subset of the atmospheric profiles is used where the two aerosol types belong to clearly distinguished layers, for the purpose of characterizing the ash in terms of lidar ratio and depolarization. These quantities are then used in a three‐component atmosphere solution scheme of the lidar equation applied to the full data set, in order to compute the optical properties of both aerosol types separately. On 16 April a thin ash layer, 100–400 m deep, is observed (average and maximum estimated ash optical depth: 0.11 and 0.2); it descends from ∼2800 to ∼1400 m altitude over a 6‐hour period. On 18 April a double ash layer, ∼400 m deep, is observed just above the morning boundary layer (average and maximum estimated ash optical depth: 0.19 and 0.27). In the afternoon the ash is entrained into the boundary layer, and the latter reaches a depth of ∼1800 m (average and maximum estimated ash optical depth: 0.1 and 0.15). An additional ash layer, with a very small optical depth, was observed on 18 April at an altitude of 3500–4000 m. By converting the lidar optical measurements using estimates of volcanic ash specific extinction, derived from other works, the observations seem to suggest approximate peak ash concentrations of ∼1500 and ∼1000 mg/m3,respectively, on the two observations dates.

Lidar observation campaign of sugar cane fires and industrial emissions in the State of São Paulo, Brazil

Brazil has an important role in the biomass burning, with the detection of approximately 100,000 burning spots in a single year (2007). Most of these spots occur in the southern part of the Amazon basin during the dry season (from August to november) and these emissions reach the southeast of the country, a highly populated region and with serious urban air pollution problems. With the growing demand on biofuels, sugarcane is considerably expanding in the state of São Paulo, being a strong contributor to the bad air quality in this region. In the state of São Paulo, the main land use are pasture and sugarcane crop, that covers around 50% and 10% of the total area, respectively. Despite the aerosol from sugarcane burning having reduced atmospheric residence time, from a few days to some weeks, they might get together with those aerosol which spread over long distances (hundreds to thousands of kilometers). In the period of June through February 2010 a LIDAR observation campaign was carried in the state of São Paulo, Brazil, in order to observe and characterize optically the aerosols from two distinct sources, namely, sugar cane biomass burning and industrial emissions. For this purpose 2 LIDAR systems were available, one mobile and the other placed in a laboratory, both working in the visible (532 nm) and additionally the mobile system had a Raman channel available (607 nm). Also this campaign counted with a SODAR, a meteorological RADAR specially set up to detect aerosol echoes and gas-particle analyzers. To guarantee a good regional coverage 4 distinct sites were available to deploy the instruments, 2 in the near field of biomass burning activities (Rio Claro and Bauru), one for industrial emissions (Cubatão) and others from urban sources (São Paulo). The whole campaign provide the equivalent of 30 days of measurements which allowed us to get aerosol optical properties such as backscattering/extinction coefficients, scatter and LIDAR ratios, those were used to correlate with air quality and meteorological indicators and quantities. In this paper we should focus on the preliminary results of the Raman LIDAR system and its derived aerosol optical quantities. © 2010 Copyright SPIE - The International Society for Optical Engineering.

Initial analysis from a Lidar observation campaign of sugar cane fires in the central and western portion of the S̃ao Paulo State - Brazil

The central and western portion of the S̃ao Paulo State has large areas of sugar cane plantations, and due to the growing demand for biofuels, the production is increasing every year. During the harvest period some plantation areas are burnt a few hours before the manual cutting, causing significant quantities of biomass burning aerosol to be injected into the atmosphere. During August 2010, a field campaign has been carried out in Ourinhos, situated in the south-western region of S̃ao Paulo State. A 2-channel Raman Lidar system and two meteorological S-Band Doppler Radars are used to indentify and quantify the biomass burning plumes. In addiction, CALIPSO Satellite observations were used to compare the aerosol optical properties detected in that region with those retrieved by Raman Lidar system. Although the campaign yielded 30 days of measurements, this paper will be focusing only one case study, when aerosols released from nearby sugar cane fires were detected by the Lidar system during a CALIPSO overpass. The meteorological radar, installed in Bauru, approximately 110 km northeast from the experimental site, had recorded echoes (dense smoke comprising aerosols) from several fires occurring close to the Raman Lidar system, which also detected an intense load of aerosol in the atmosphere. HYSPLIT model forward trajectories presented a strong indication that both instruments have measured the same air masss parcels, corroborated with the Lidar Ratio values from the 532 nm elastic and 607 nm Raman N2 channel analyses and data retrieved from CALIPSO have indicated the predominance of aerosol from biomass burning sources. © 2011 SPIE.

Raman Lidar monitors emissions from sugar cane fires in the State of São Paulo: a pilot-project integrating radar, sodar, aerosol and gas observations

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

Cirrus clouds observation in Santa Maria, Rio Grande do Sul during the experiment Chuva - Sul.

Cirrus clouds are an interesting point in the research of the atmosphere due their behavior and the effect on the earth radiation budget. They can affect the atmospheric radiation budget by reflecting the incoming solar radiation and absorbing the outgoing terrestrial radiation. Also, this cloud type is involved in the dehydration of the upper troposphere and lower stratosphere. So, it is interesting to increment the measurements of this type of clouds from the ground. During November and December 2012, through the CHUVA-SUL campaign, measurements with lidar in Santa Maria, Rio Grande do Sul were conducted. The system installed in Santa Maria site (29.8 °S; 53.7 °W, 100 m asl) was a single elastic-backscatter lidar using the wavelength of 532 nm. Some days with cirrus clouds lidar measurements were detected. Four days with presence of cirrus cloud are showed in the present study. These days, 7, 8, 19 and 28 November 2012, was selected due the persistence of cirrus clouds over many hours. The raw retrieval lidar signals and inverted backscatter coefficient profiles were analyzed for the selected days. Base and top height was obtained by analysis of raw signal and backscatter coefficient. Extinction coefficient profiles were obtained by the assumption of the lidar ratio. Cirrus cloud optical depth (COD) values were calculated, from the integration of the extinction coefficient between the base and top altitudes of the cirrus clouds.

LOOKING INTO CALIPSO CLIMATOLOGICAL PRODUCTS: EVALUATION AND SUGGESTIONS FROM EARLINET

CALIPSO (Cloud-Aerosol Lidar and Pathfinder Satellite Observations) Level 3 (CL3) data were compared against EARLINET (European Aerosol Research Lidar Network) monthly averages obtained by profiles during satellite overpasses. Data from EARLINET stations of Évora, Granada, Leipzig, Naples and Potenza, equipped with advanced multi-wavelength Raman lidars were used for this study. Owing to spatial and temporal differences, we reproduced the CL3 filtering rubric onto the CALIPSO Level 2 data. The CALIPSO monthly mean profiles following this approach are called CALIPSO Level 3*, CL3*. This offers the possibility to achieve direct comparable datasets. In respect to CL3 data, the agreement typically improved, in particular above the areas directly affected by the anthropogenic activities within the planetary boundary layer. However in most of the cases a subtle CALIPSO underestimation was observed with an average bias of 0.03 km-1. We investigated the backscatter coefficient applying the same screening criteria, where the mean relative difference in respect to the extinction comparison improved from 15.2% to 11.4%. Lastly, the typing capabilities of CALIPSO were assessed outlining the importance of the correct aerosol type (and associated lidar ratio value) assessment to the CALIPSO aerosol properties retrieval.

Independent evaluation of the ability of spaceborne radar and lidar to retrieve the microphysical and radiative properties of ice clouds

The combination of radar and lidar in space offers the unique potential to retrieve vertical profiles of ice water content and particle size globally, and two algorithms developed recently claim to have overcome the principal difficulty with this approach-that of correcting the lidar signal for extinction. In this paper "blind tests" of these algorithms are carried out, using realistic 94-GHz radar and 355-nm lidar backscatter profiles simulated from aircraft-measured size spectra, and including the effects of molecular scattering, multiple scattering, and instrument noise. Radiation calculations are performed on the true and retrieved microphysical profiles to estimate the accuracy with which radiative flux profiles could be inferred remotely. It is found that the visible extinction profile can be retrieved independent of assumptions on the nature of the size distribution, the habit of the particles, the mean extinction-to-backscatter ratio, or errors in instrument calibration. Local errors in retrieved extinction can occur in proportion to local fluctuations in the extinction-to-backscatter ratio, but down to 400 m above the height of the lowest lidar return, optical depth is typically retrieved to better than 0.2. Retrieval uncertainties are greater at the far end of the profile, and errors in total optical depth can exceed 1, which changes the shortwave radiative effect of the cloud by around 20%. Longwave fluxes are much less sensitive to errors in total optical depth, and may generally be calculated to better than 2 W m(-2) throughout the profile. It is important for retrieval algorithms to account for the effects of lidar multiple scattering, because if this is neglected, then optical depth is underestimated by approximately 35%, resulting in cloud radiative effects being underestimated by around 30% in the shortwave and 15% in the longwave. Unlike the extinction coefficient, the inferred ice water content and particle size can vary by 30%, depending on the assumed mass-size relationship (a problem common to all remote retrieval algorithms). However, radiative fluxes are almost completely determined by the extinction profile, and if this is correct, then errors in these other parameters have only a small effect in the shortwave (around 6%, compared to that of clear sky) and a negligible effect in the longwave.

Estimating drizzle drop size and precipitation rate using two-colour lidar measurements

A method to estimate the size and liquid water content of drizzle drops using lidar measurements at two wavelengths is described. The method exploits the differential absorption of infrared light by liquid water at 905 nm and 1.5 μm, which leads to a different backscatter cross section for water drops larger than ≈50 μm. The ratio of backscatter measured from drizzle samples below cloud base at these two wavelengths (the colour ratio) provides a measure of the median volume drop diameter D0. This is a strong effect: for D0=200 μm, a colour ratio of ≈6 dB is predicted. Once D0 is known, the measured backscatter at 905 nm can be used to calculate the liquid water content (LWC) and other moments of the drizzle drop distribution. The method is applied to observations of drizzle falling from stratocumulus and stratus clouds. High resolution (32 s, 36 m) profiles of D0, LWC and precipitation rate R are derived. The main sources of error in the technique are the need to assume a value for the dispersion parameter μ in the drop size spectrum (leading to at most a 35% error in R) and the influence of aerosol returns on the retrieval (≈10% error in R for the cases considered here). Radar reflectivities are also computed from the lidar data, and compared to independent measurements from a colocated cloud radar, offering independent validation of the derived drop size distributions.

A method for estimating the turbulent kinetic energy dissipation rate from a vertically pointing Doppler lidar, and independent evaluation from balloon-borne in situ measurements

A method of estimating dissipation rates from a vertically pointing Doppler lidar with high temporal and spatial resolution has been evaluated by comparison with independent measurements derived from a balloon-borne sonic anemometer. This method utilizes the variance of the mean Doppler velocity from a number of sequential samples and requires an estimate of the horizontal wind speed. The noise contribution to the variance can be estimated from the observed signal-to-noise ratio and removed where appropriate. The relative size of the noise variance to the observed variance provides a measure of the confidence in the retrieval. Comparison with in situ dissipation rates derived from the balloon-borne sonic anemometer reveal that this particular Doppler lidar is capable of retrieving dissipation rates over a range of at least three orders of magnitude. This method is most suitable for retrieval of dissipation rates within the convective well-mixed boundary layer where the scales of motion that the Doppler lidar probes remain well within the inertial subrange. Caution must be applied when estimating dissipation rates in more quiescent conditions. For the particular Doppler lidar described here, the selection of suitably short integration times will permit this method to be applicable in such situations but at the expense of accuracy in the Doppler velocity estimates. The two case studies presented here suggest that, with profiles every 4 s, reliable estimates of ϵ can be derived to within at least an order of magnitude throughout almost all of the lowest 2 km and, in the convective boundary layer, to within 50%. Increasing the integration time for individual profiles to 30 s can improve the accuracy substantially but potentially confines retrievals to within the convective boundary layer. Therefore, optimization of certain instrument parameters may be required for specific implementations.

A novel satellite mission concept for upper air water vapour, aerosol and cloud observations using integrated path differential absorption LiDAR limb sounding

We propose a new satellite mission to deliver high quality measurements of upper air water vapour. The concept centres around a LiDAR in limb sounding by occultation geometry, designed to operate as a very long path system for differential absorption measurements. We present a preliminary performance analysis with a system sized to send 75 mJ pulses at 25 Hz at four wavelengths close to 935 nm, to up to 5 microsatellites in a counter-rotating orbit, carrying retroreflectors characterized by a reflected beam divergence of roughly twice the emitted laser beam divergence of 15 µrad. This provides water vapour profiles with a vertical sampling of 110 m; preliminary calculations suggest that the system could detect concentrations of less than 5 ppm. A secondary payload of a fairly conventional medium resolution multispectral radiometer allows wide-swath cloud and aerosol imaging. The total weight and power of the system are estimated at 3 tons and 2,700 W respectively. This novel concept presents significant challenges, including the performance of the lasers in space, the tracking between the main spacecraft and the retroreflectors, the refractive effects of turbulence, and the design of the telescopes to achieve a high signal-to-noise ratio for the high precision measurements. The mission concept was conceived at the Alpbach Summer School 2010.

Effective LAI and CHP of a single tree from small-footprint full-waveform LiDAR

This letter has tested the canopy height profile (CHP) methodology as a way of effective leaf area index (LAIe) and vertical vegetation profile retrieval at a single-tree level. Waveform and discrete airborne LiDAR data from six swaths, as well as from the combined data of six swaths, were used to extract the LAIe of a single live Callitris glaucophylla tree. LAIe was extracted from raw waveform as an intermediate step in the CHP methodology, with two different vegetation-ground reflectance ratios. Discrete point LAIe estimates were derived from the gap probability using the following: 1) single ground returns and 2) all ground returns. LiDAR LAIe retrievals were subsequently compared to hemispherical photography estimates, yielding mean values within ±7% of the latter, depending on the method used. The CHP of a single dead Callitris glaucophylla tree, representing the distribution of vegetation material, was verified with a field profile manually reconstructed from convergent photographs taken with a fixed-focal-length camera. A binwise comparison of the two profiles showed very high correlation between the data reaching R2 of 0.86 for the CHP from combined swaths. Using a study-area-adjusted reflectance ratio improved the correlation between the profiles, but only marginally in comparison to using an arbitrary ratio of 0.5 for the laser wavelength of 1550 nm.

Validation of Canopy Height Profile methodology for small-footprint full-waveform airborne LiDAR data in a discontinuous canopy environment

A Canopy Height Profile (CHP) procedure presented in Harding et al. (2001) for large footprint LiDAR data was tested in a closed canopy environment as a way of extracting vertical foliage profiles from LiDAR raw-waveform. In this study, an adaptation of this method to small-footprint data has been shown, tested and validated in an Australian sparse canopy forest at plot- and site-level. Further, the methodology itself has been enhanced by implementing a dataset-adjusted reflectance ratio calculation according to Armston et al. (2013) in the processing chain, and tested against a fixed ratio of 0.5 estimated for the laser wavelength of 1550nm. As a by-product of the methodology, effective leaf area index (LAIe) estimates were derived and compared to hemispherical photography-derived values. To assess the influence of LiDAR aggregation area size on the estimates in a sparse canopy environment, LiDAR CHPs and LAIes were generated by aggregating waveforms to plot- and site-level footprints (plot/site-aggregated) as well as in 5m grids (grid-processed). LiDAR profiles were then compared to leaf biomass field profiles generated based on field tree measurements. The correlation between field and LiDAR profiles was very high, with a mean R2 of 0.75 at plot-level and 0.86 at site-level for 55 plots and the corresponding 11 sites. Gridding had almost no impact on the correlation between LiDAR and field profiles (only marginally improvement), nor did the dataset-adjusted reflectance ratio. However, gridding and the dataset-adjusted reflectance ratio were found to improve the correlation between raw-waveform LiDAR and hemispherical photography LAIe estimates, yielding the highest correlations of 0.61 at plot-level and of 0.83 at site-level. This proved the validity of the approach and superiority of dataset-adjusted reflectance ratio of Armston et al. (2013) over a fixed ratio of 0.5 for LAIe estimation, as well as showed the adequacy of small-footprint LiDAR data for LAIe estimation in discontinuous canopy forests.

Preliminary leaf area index estimates from airborne small footprint full-waveform LiDAR data

This study has compared preliminary estimates of effective leaf area index (LAI) derived from fish-eye lens photographs to those estimated from airborne full-waveform small-footprint LiDAR data for a forest dataset in Australia. The full-waveform data was decomposed and optimized using a trust-region-reflective algorithm to extract denser point clouds. LAI LiDAR estimates were derived in two ways (1) from the probability of discrete pulses reaching the ground without being intercepted (point method) and (2) from raw waveform canopy height profile processing adapted to small-footprint laser altimetry (waveform method) accounting for reflectance ratio between vegetation and ground. The best results, that matched hemispherical photography estimates, were achieved for the waveform method with a study area-adjusted reflectance ratio of 0.4 (RMSE of 0.15 and 0.03 at plot and site level, respectively). The point method generally overestimated, whereas the waveform method with an arbitrary reflectance ratio of 0.5 underestimated the fish-eye lens LAI estimates.

A generalised background correction algorithm for a Halo Doppler lidar and its application to data from Finland

Current commercially available Doppler lidars provide an economical and robust solution for measuring vertical and horizontal wind velocities, together with the ability to provide co- and cross-polarised backscatter profiles. The high temporal resolution of these instruments allows turbulent properties to be obtained from studying the variation in radial velocities. However, the instrument specifications mean that certain characteristics, especially the background noise behaviour, become a limiting factor for the instrument sensitivity in regions where the aerosol load is low. Turbulent calculations require an accurate estimate of the contribution from velocity uncertainty estimates, which are directly related to the signal-to-noise ratio. Any bias in the signal-to-noise ratio will propagate through as a bias in turbulent properties. In this paper we present a method to correct for artefacts in the background noise behaviour of commercially available Doppler lidars and reduce the signal-to-noise ratio threshold used to discriminate between noise, and cloud or aerosol signals. We show that, for Doppler lidars operating continuously at a number of locations in Finland, the data availability can be increased by as much as 50 % after performing this background correction and subsequent reduction in the threshold. The reduction in bias also greatly improves subsequent calculations of turbulent properties in weak signal regimes.

Dust and smoke transport from Africa to South America: Lidar profiling over Cape Verde and the Amazon rainforest

Quasi-simultaneous vertically resolved multiwavelength aerosol Raman lidar observations were conducted in the near field (Praia, Cape Verde, 15 degrees N, 23.5 degrees W) and in the far field (Manaus, Amazon basin, Brazil, 2.5 degrees S, 60 degrees W) of the long-range transport regime between West Africa and South America. Based on a unique data set (case study) of spectrally resolved backscatter and extinction coefficients, and of the depolarization ratio a detailed characterization of aerosol properties, vertical stratification, mixing, and aging behavior during the long-distance travel in February 2008 (dry season in western Africa, wet season in the Amazon basin) is presented. While highly stratified aerosol layers of dust and smoke up to 5.5 km height were found close to Africa, the aerosol over Manaus was almost well-mixed, reached up to 3.5 km, and mainly consisted of aged biomass burning smoke. Citation: Ansmann, A., H. Baars, M. Tesche, D. Muller, D. Althausen, R. Engelmann, T. Pauliquevis, and P. Artaxo (2009), Dust and smoke transport from Africa to South America: Lidar profiling over Cape Verde and the Amazon rainforest, Geophys. Res. Lett., 36, L11802, doi: 10.1029/2009GL037923.