CALIPSO Lidar Calibration at 532 nm: Version 4 Nighttime Algorithm.

Data products from the Cloud-Aerosol Lidar with Orthogonal Polarization (CALIOP) on board Cloud-Aerosol Lidar and Infrared Pathfinder Satellite Observations (CALIPSO) were recently updated following the implementation of new (version 4) calibration algorithms for all of the level 1 attenuated backscatter measurements. In this work we present the motivation for and the implementation of the version 4 nighttime 532 nm parallel channel calibration. The nighttime 532 nm calibration is the most fundamental calibration of CALIOP data, since all of CALIOP's other radiometric calibration procedures - i.e., the 532 nm daytime calibration and the 1064 nm calibrations during both nighttime and daytime - depend either directly or indirectly on the 532 nm nighttime calibration. The accuracy of the 532 nm nighttime calibration has been significantly improved by raising the molecular normalization altitude from 30-34 km to 36-39 km to substantially reduce stratospheric aerosol contamination. Due to the greatly reduced molecular number density and consequently reduced signal-to-noise ratio (SNR) at these higher altitudes, the signal is now averaged over a larger number of samples using data from multiple adjacent granules. As well, an enhanced strategy for filtering the radiation-induced noise from high energy particles was adopted. Further, the meteorological model used in the earlier versions has been replaced by the improved MERRA-2 model. An aerosol scattering ratio of 1.01 ± 0.01 is now explicitly used for the calibration altitude. These modifications lead to globally revised calibration coefficients which are, on average, 2-3% lower than in previous data releases. Further, the new calibration procedure is shown to eliminate biases at high altitudes that were present in earlier versions and consequently leads to an improved representation of stratospheric aerosols. Validation results using airborne lidar measurements are also presented. Biases relative to collocated measurements acquired by the Langley Research Center (LaRC) airborne high spectral resolution lidar (HSRL) are reduced from 3.6% ± 2.2% in the version 3 data set to 1.6% ± 2.4 % in the version 4 release.

CALIPSO (IIR-CALIOP) Retrievals of Cirrus Cloud Ice Particle Concentrations.

A new satellite remote sensing method is described whereby the sensitivity of thermal infrared wave resonance absorption to small ice crystals is exploited to estimate cirrus cloud ice particle number concentration N, effective diameter De, and ice water content IWC. This method uses co-located observations from the Infrared Imaging Radiometer (IIR) and from the CALIOP (Cloud and Aerosol Lidar with Orthogonal Polarization) lidar aboard the CALIPSO (Cloud-Aerosol Lidar and Infrared Pathfinder Satellite Observation) polar orbiting satellite, employing IIR channels at 10.6 µm and 12.05 µm. Using particle size distributions measured over several flights of the TC4 (Tropical Composition, Cloud and Climate Coupling) and the mid-latitudes SPARTICUS (Small Particles in Cirrus) field campaigns, we show for the first time that N/IWC is tightly related to ßeff; the ratio of effective absorption optical depths at 12.05 µm and 10.6 µm. Relationships developed from in situ aircraft measurements are applied to ßeff derived from IIR measurements to retrieve N. This satellite remote sensing method is constrained by measurements of ßeff from the IIR and is by essence sensitive to the smallest ice crystals. Retrieval uncertainties are discussed, including uncertainties related to in situ measurement of small ice crystals (D < 15 µm), which are studied through comparisons with IIR ßeff. The method is applied here to single-layered semi-transparent clouds having a visible optical depth between about 0.3 and 3, where cloud base temperature is ≤ 235 K. Two years of CALIPSO data have been analyzed for the years 2008 and 2013, with the dependence of cirrus cloud N and De on altitude, temperature, latitude, season (winter vs. summer) and topography (land vs. ocean) described. The results for the mid-latitudes show a considerable dependence on season. In the high latitudes, N tends to be highest and De smallest, whereas the opposite is true for the tropics. The frequency of occurrence of these relatively thick cirrus clouds exhibited a strong seasonal dependence in the high latitudes, with the occurrence frequency during Arctic winter being at least twice that of any other season. Processes that could potentially explain some of these micro-and macroscopic cloud phenomena are discussed.

CALIPSO IIR Version 2 Level 1b calibrated radiances: analysis and reduction of residual biases in the Northern Hemisphere.

Version 2 of the Level 1b calibrated radiances of the Imaging Infrared Radiometer (IIR) on board the Cloud-Aerosol Lidar and Infrared Satellite Observation (CALIPSO) satellite has been released recently. This new version incorporates corrections of small but systematic seasonal calibration biases previously revealed in Version 1 data products mostly north of 30° N. These biases of different amplitudes in the three IIR channels 8.65 µm (IIR1), 10.6 µm (IIR2), and 12.05 µm (IIR3) were made apparent by a striping effect in images of IIR inter-channel brightness temperature differences (BTDs) and through seasonal warm biases of nighttime IIR brightness temperatures in the 30-60° N latitude range. The latter were highlighted through observed and simulated comparisons with similar channels of the Moderate Resolution Imaging Spectroradiometer (MODIS) on board the Aqua spacecraft. To characterize the calibration biases affecting Version 1 data, a semi-empirical approach is developed, which is based on the in-depth analysis of the IIR internal calibration procedure in conjunction with observations such as statistical comparisons with similar MODIS/Aqua channels. Two types of calibration biases are revealed: an equalization bias affecting part of the individual IIR images and a global bias affecting the radiometric level of each image. These biases are observed only when the temperature of the instrument increases, and they are found to be functions of elapsed time since night-to-day transition, regardless of the season. Correction coefficients of Version 1 radiances could thus be defined and implemented in the Version 2 code. As a result, the striping effect seen in Version 1 is significantly attenuated in Version 2. Systematic discrepancies between nighttime and daytime IIR-MODIS BTDs in the 30-60° N latitude range in summer are reduced from 0.2 K in Version 1 to 0.1 K in Version 2 for IIR1-MODIS29. For IIR2-MODIS31 and IIR3-MODIS32, they are reduced from 0.4 K to close to zero, except for IIR3-MODIS32 in June, where the night-minus-day difference is around -0.1 K.

Atmospheric boundary layer CO2 remote sensing with a direct detection LIDAR instrument based on a widely tunable optical parametric source.

We report on the capability of a direct detection differential absorption lidar (DIAL) for range resolved and integrated path (IPDIAL) remote sensing of CO2 in the atmospheric boundary layer (ABL). The laser source is an amplified nested cavity optical parametric oscillator (NesCOPO) emitting approximately 8 mJ at the two measurement wavelengths selected near 2050 nm. Direct detection atmospheric measurements are taken from the ground using a 30 Hz frequency switching between emitted wavelengths. Results show that comparable precision measurements are achieved in DIAL and IPDIAL modes (not better than a few ppm) on high SNR targets such as near range ABL aerosol and clouds, respectively. Instrumental limitations are analyzed and degradation due to cloud scattering variability is discussed to explain observed DIAL and IPDIAL limitations.

IAOOS microlidar-on-buoy development and first atmospheric observations obtained during 2014 and 2015 arctic drifts.

The in situ tests of first ever autonomous aerosol and cloud backscatter LIDAR (light detection and ranging) systems implemented on buoys for Arctic observations has been achieved in 2015 within the French EQUIPEX IAOOS project. The environmental and operational constraints were met by adopting a concept of a fibered microjoule lidar system using a laser diode. Two systems have been developed with and without polarization analysis capability. A specific optical design was used for polarization discrimination. These systems were integrated in buoys and tested in the Arctic in 2014 and 2015 at latitudes higher than 80°N. Data were transmitted through an Iridium space link. Measurements have been obtained 90% of the time from the non-polarized system in 2014 over 8 months as the first fully equipped buoy drifted from the Barneo Russian camp close to the North Pole toward Svalbard. A polarized system was then tested over a short period in winter 2015 north of Svalbard during the Norwegian campaign N-ICE. In April and May 2014, the unattended lidar measurements showed a large occurrence of aerosols and haze. The average attenuated scattering ratio for non-cloudy profiles during this period was about 2.2. Aerosols could reach an altitude of 5km on average, whereas over the rest of the period low level clouds (below 1000 m) were prevailing with an average attenuated scattering ratio of about 103. The main features of the developed lidar instruments and first results are presented here.

Retrieval of ocean subsurface particulate backscattering coefficient from space-borne CALIOP lidar measurements.

A new approach has been proposed to determine ocean subsurface particulate backscattering coefficient bbp from CALIOP 30° off-nadir lidar measurements. The new method also provides estimates of the particle volume scattering function at the 180° scattering angle. The CALIOP based layer-integrated lidar backscatter and particulate backscattering coefficients are compared with the results obtained from MODIS ocean color measurements. The comparison analysis shows that ocean subsurface lidar backscatter and particulate backscattering coefficient bbp can be accurately obtained from CALIOP lidar measurements, thereby supporting the use of space-borne lidar measurements for ocean subsurface studies.

Direct-detection wind lidar operating with a multimode laser.

A direct-detection wind lidar that operates with a multimode laser has been developed and tested. The instrument exploits the light backscattered by particles using a Mach-Zehnder interferometer with an optical path difference matched to the free spectral range of the laser longitudinal modes. In addition to requiring no monomodal emission, the system requires no frequency locking between the interferometer and the laser. We report laboratory and atmospheric measurements that show that the lidar is capable of measuring the radial wind velocity with a systematic error lower than 1 ms(-1) and a random error lower than 2 ms(-1) for a signal-to-noise ratio of 100. The development is motivated by the possibility to probe wind with a compact system in planetary atmospheres.

CALIPSO lidar ratio retrieval over the ocean.

We are demonstrating on a few cases the capability of CALIPSO to retrieve the 532 nm lidar ratio over the ocean when CloudSat surface scattering cross section is used as a constraint. We are presenting the algorithm used and comparisons with the column lidar ratio retrieved by the NASA airborne high spectral resolution lidar. For the three cases presented here, the agreement is fairly good. The average CALIPSO 532 nm column lidar ratio bias is 13.7% relative to HSRL, and the relative standard deviation is 13.6%. Considering the natural variability of aerosol microphysical properties, this level of accuracy is significant since the lidar ratio is a good indicator of aerosol types. We are discussing dependencies of the accuracy of retrieved aerosol lidar ratio on atmospheric aerosol homogeneity, lidar signal to noise ratio, and errors in the optical depth retrievals. We are obtaining the best result (bias 7% and standard deviation around 6%) for a nighttime case with a relatively constant lidar ratio (in the vertical) indicative of homogeneous aerosol type.

Lidar equation for ocean surface and subsurface.

The lidar equation for ocean at optical wavelengths including subsurface signals is revisited using the recent work of the radiative transfer and ocean color community for passive measurements. The previous form of the specular and subsurface echo term are corrected from their heritage, which originated from passive remote sensing of whitecaps, and is improved for more accurate use in future lidar research. A corrected expression for specular and subsurface lidar return is presented. The previous formalism does not correctly address angular dependency of specular lidar return and overestimates the subsurface term by a factor ranging from 89% to 194% for a nadir pointing lidar. Suggestions for future improvements to the lidar equation are also presented.

Study of the statistics of water vapor mixing ratio determined from Raman lidar measurements.

The statistical properties of atmospheric water vapor mixing ratio (WVMR) determined as the ratio of Raman lidar signals backscattered from water vapor and nitrogen molecules are studied. It is shown that WVMR estimates can be biased by a small percentage at low signal photon-counting rates due to fluctuations in the nitrogen signal in the denominator of the ratio, the magnitude of the bias being linked to the signal-to-noise ratio of the nitrogen signal. This is particularly important when unbiased estimates are required as in the case of climate studies and global positioning system (GPS) signal calibration. Different bias corrections and a modified ratio formulation are proposed in order to correct or eliminate this bias. The method is successfully applied in processing signals obtained with an experimental Raman lidar system devoted to calibrate GPS signals for slant path delays. It is shown to reduce biases into negligible values in both WVMR and wet path delay estimates in the range interval of 0-7 km.

Complementary study of differential absorption lidar optimization in direct and heterodyne detections.

A detailed study using both analytical and numerical calculations of direct and heterodyne differential absorption lidar (DIAL) techniques is conducted to complement previous studies. The DIAL measurement errors depend on key experimental parameters, some of which can be adjusted to minimize the statistical error. Accordingly, the pertinent criteria on optical thickness, the number of photons emitted at the on and off wavelengths, are discussed to reduce the relative error on the total column content or range-resolved measurements that rely on either hard target or atmospheric backscatter returns. In direct detection, the optimal optical thickness decreases from 1.3 to 0.8 when the background increases while the on-line-to-off-line optimal energy ratio decreases from 3.6 to 2.7. In heterodyne detection, the minimum error is obtained for an optical thickness of 1.2 and an energy ratio of 4.3.

Simultaneous measurements of particle backscattering and extinction coefficients and wind velocity by lidar with a Mach-Zehnder interferometer: principle of operation and performance assessment.

The development of remote-sensing instruments that can be used to monitor several parameters at the same time is important for the study of complex processes such as those that control climate and environment. In this paper the performance of a new concept of lidar receiver that allows for the direct measurement of aerosol and cloud optical properties simultaneously with wind velocity is investigated. This receiver uses a Mach-Zehnder interferometer. Two different configurations, either with four photometric output channels or with fringe imaging on a multichannel detector, are studied. Analytical expressions of the statistical errors are given under the assumption of Gaussian signal spectra. It is shown that similar accuracies can be achieved for both configurations. Performance modeling of the retrieval of semitransparent cloud optical scattering properties and wind velocity was done at different operation wavelengths for a Nd:YAG laser source. Results for such a lidar system onboard an aircraft flying at an altitude of 12 km show that for semitransparent clouds the best results were obtained at 355 nm, with relative standard deviations of 0.5% and 5% for the backscatter and extinction coefficients, respectively, together with a velocity accuracy of 0.2 ms(-1). The accuracy of optical properties retrieved for boundary layer aerosols are comparable, whereas the velocity accuracy is decreased to 1 ms(-1). Finally, an extrapolation to a large 355-nm spaceborne lidar shows accuracies in the range from 2.5% to 5% for the backscatter coefficient and from 10% to 15% for the extinction coefficient together with a vertical wind speed accuracy of better than 0.5 ms(-1) for semitransparent clouds and boundary layer, with a vertical resolution of 500 m and a 100 shot averaging.

Variational method for the retrieval of the optical thickness and the backscatter coefficient from multiangle lidar profiles.

A variational method for retrieving the aerosol optical thickness and backscatter coefficient profiles from multiangle lidar measurements is presented and discussed. A monostatic single-wavelength low-energy lidar system was operated at different zenith angles during the Indian Ocean Experiment (INDOEX) campaign in 1999 to characterize the aerosol plumes in the Indian monsoon. The variational method was applied to lidar data to retrieve profiles of optical thickness and the backscatter coefficient for nighttime and daytime measurements. Results are obtained with an uncertainty of 10% below 3 km (nighttime) and 2.8 km (daytime) and a bias of less than 0.01. During daytime the retrieval of optical parameters is indeed limited to a lower altitude owing to the sky background signal and the atmospheric inhomogeneity. In both cases the total aerosol optical thickness is consistent (+/- 10%) with the integrated value derived from sunphotometer measurements. Backscatter-to-extinction ratios estimated in different regions by two distinct methods compared well, which proves the capability of the method to assess optical measurements and account for the altitude dependence of the phase function.