scholarly journals Retrieval of temperature from a multiple channel pure rotational Raman backscatter lidar using an optimal estimation method

2019 ◽  
Vol 12 (11) ◽  
pp. 5801-5816 ◽  
Author(s):  
Shayamila Mahagammulla Gamage ◽  
Robert J. Sica ◽  
Giovanni Martucci ◽  
Alexander Haefele

Abstract. We present a new method for retrieving temperature from pure rotational Raman (PRR) lidar measurements. Our optimal estimation method (OEM) used in this study uses the full physics of PRR scattering and does not require any assumption of the form for a calibration function nor does it require fitting of calibration factors over a large range of temperatures. The only calibration required is the estimation of the ratio of the lidar constants of the two PRR channels (coupling constant) that can be evaluated at a single or multiple height bins using a simple analytic expression. The uncertainty budget of our OEM retrieval includes both statistical and systematic uncertainties, including the uncertainty in the determination of the coupling constant on the temperature. We show that the error due to calibration can be reduced significantly using our method, in particular in the upper troposphere when calibration is only possible over a limited temperature range. Some other advantages of our OEM over the traditional Raman lidar temperature retrieval algorithm include not requiring correction or gluing to the raw lidar measurements, providing a cutoff height for the temperature retrievals that specifies the height to which the retrieved profile is independent of the a priori temperature profile, and the retrieval's vertical resolution as a function of height. The new method is tested on PRR temperature measurements from the MeteoSwiss RAman Lidar for Meteorological Observations system in clear and cloudy sky conditions, compared to temperature calculated using the traditional PRR calibration formulas, and validated with coincident radiosonde temperature measurements in clear and cloudy conditions during both daytime and nighttime.

2019 ◽  
Author(s):  
Shayamila Mahagammulla Gamage ◽  
Robert J. Sica ◽  
Giovanni Martucci ◽  
Alexander Haefele

Abstract. We present a new method for retrieving temperature from Pure Rotational Raman (PRR) lidar measurements. Our Optimal Estimation Method (OEM) used in this study uses the full physics of PRR scattering and does not require any assumption of the form for a calibration function nor does it require fitting of calibration factors over a large range of temperatures. The only calibration required is the estimation of the ratio of the lidar constants of the two PRR channels (coupling constant) that can be evaluated at a single or multiple height bins using a simple analytic expression. The uncertainty budget of our OEM retrieval includes both statistical and systematic uncertainties, including the uncertainty in the determination of the coupling constant on the temperature. We show that the error due to calibration can be reduced significantly using our method, in particular in the upper troposphere when calibration is only possible over a limited temperature range. Some other advantages of our OEM over the traditional Raman lidar temperature retrieval algorithm include not requiring correction or gluing to the raw lidar measurements, providing a cutoff height for the temperature retrievals that specifies the height to which the retrieved profile is independent of the a priori temperature profile, and the retrieval's vertical resolution as a function of height. The new method is tested on PRR temperature measurements from the MeteoSwiss Raman Lidar for Meteorological Observations system in different sky conditions, compared to temperature calculated using the traditional PRR calibration formulas, and validated with coincident radiosonde temperature measurements in clear and cloudy conditions during both day and night time.


2020 ◽  
Vol 37 (11) ◽  
pp. 2051-2064
Author(s):  
S. Mahagammulla Gamage ◽  
R. J. Sica ◽  
G. Martucci ◽  
A. Haefele

AbstractWe present a one-dimensional variational (1D Var) retrieval of fifth-generation European Centre for Medium-Range Weather Forecast reanalysis (ERA5) temperature and relative humidity profiles above Payerne, Switzerland, assimilating raw backscatter measurements from the MeteoSwiss Raman Lidar for Meteorological Observations (RALMO). Our reanalysis is called ERA5-reRH. We use an optimal estimation method to perform the 1D Var data retrieval. The forward model combines the Raman lidar equation with the Hyland and Wexler expression for water vapor saturation pressure. The error covariance matrix of ERA5 was derived from the differences between ERA5 and a set of 50 special radiosoundings that have not been assimilated into ERA5. We validate ERA5-reRH, ERA5, and RALMO temperature and relative humidity profiles against the same set of special radiosoundings and found the best agreement was with our reanalysis, with a bias of less than 2% relative humidity with respect to water (%RHw) and a spread of less than 8%RHw below 8 km in terms of relative humidity. Improvements for temperature in our reanalysis are only found in the boundary layer, as ERA5 assimilates a large number of upper-air temperature observations. Our retrieval also provides a full uncertainty budget of the reanalyzed temperature and relative humidity including both random and systematic uncertainties.


2018 ◽  
Vol 176 ◽  
pp. 01025
Author(s):  
R. J. Sica ◽  
A. Haefele ◽  
A. Jalali ◽  
S. Gamage ◽  
G. Farhani

The optimal estimation method (OEM) has a long history of use in passive remote sensing, but has only recently been applied to active instruments like lidar. The OEM’s advantage over traditional techniques includes obtaining a full systematic and random uncertainty budget plus the ability to work with the raw measurements without first applying instrument corrections. In our meeting presentation we will show you how to use the OEM for temperature and composition retrievals for Rayleigh-scatter, Ramanscatter and DIAL lidars.


2018 ◽  
Vol 176 ◽  
pp. 01011
Author(s):  
S. Mahagammulla Gamage ◽  
A. Haefele ◽  
R.J. Sica

We present the application of the Optimal Estimation Method (OEM) to retrieve atmospheric temperatures from pure rotational Raman (PRR) backscatter lidar measurements. A forward model (FM) is developed to retrieve temperature and tested using synthetic measurements. The OEM offers many advantages for this analysis, including eliminating the need to determine temperature calibration coefficients.


2021 ◽  
Author(s):  
Arno Keppens ◽  
Jean-Christopher Lambert ◽  
Daan Hubert ◽  
Steven Compernolle ◽  
Tijl Verhoelst ◽  
...  

<p>Part of the space segment of EU’s Copernicus Earth Observation programme, the Sentinel-5 Precursor (S5P) mission is dedicated to global and European atmospheric composition measurements of air quality, climate and the stratospheric ozone layer. On board of the S5P early afternoon polar satellite, the imaging spectrometer TROPOMI (TROPOspheric Monitoring Instrument) performs nadir measurements of the Earth radiance within the UV-visible and near-infrared spectral ranges, from which atmospheric ozone profile data are retrieved. Developed at the Royal Netherlands Meteorological Institute (KNMI) and based on the optimal estimation method, TROPOMI’s operational ozone profile retrieval algorithm has recently been upgraded. With respect to early retrieval attempts, accuracy is expected to have improved significantly, also thanks to recent updates of the TROPOMI Level-1b data product. This work reports on the initial validation of the improved TROPOMI height-resolved ozone data in the troposphere and stratosphere, as collected both from the operational S5P Mission Performance Centre/Validation Data Analysis Facility (MPC/VDAF) and from the S5PVT scientific project CHEOPS-5p. Based on the same validation best practices as developed for and applied to heritage sensors like GOME-2, OMI and IASI (Keppens et al., 2015, 2018), the validation methodology relies on the analysis of data retrieval diagnostics – like the averaging kernels’ information content – and on comparisons of TROPOMI data with reference ozone profile measurements. The latter are acquired by ozonesonde, stratospheric lidar, and tropospheric lidar stations performing network operation in the context of WMO's Global Atmosphere Watch and its contributing networks NDACC and SHADOZ. The dependence of TROPOMI’s ozone profile uncertainty on several influence quantities like cloud fraction and measurement parameters like sun and scan angles is examined and discussed. This work concludes with a set of quality indicators, enabling users to verify the fitness-for-purpose of the S5P data.</p>


2020 ◽  
Author(s):  
Arno Keppens ◽  
Daan Hubert ◽  
Jean-Christopher Lambert ◽  
Steven Compernolle ◽  
Tijl Verhoelst ◽  
...  

<p>Part of the space segment of EU’s Copernicus Earth Observation programme, the Sentinel-5 Precursor (S5P) mission is dedicated to global and European atmospheric composition measurements of air quality, climate and the stratospheric ozone layer. On board of the S5P early afternoon polar satellite, the imaging spectrometer TROPOMI (TROPOspheric Monitoring Instrument) performs nadir measurements of the Earth radiance within the UV-visible and near-infrared spectral ranges, from which atmospheric ozone profile data are retrieved. Developed at the Royal Netherlands Meteorological Institute (KNMI) and based on the optimal estimation method, TROPOMI’s operational ozone profile retrieval algorithm has recently been upgraded. With respect to early retrieval attempts, accuracy is expected to have improved significantly, also thanks to recent updates of the TROPOMI Level-1b data product. This work reports on the initial validation of the improved TROPOMI height-resolved ozone data in the troposphere and stratosphere, as collected both from the operational S5P Mission Performance Centre/Validation Data Analysis Facility (MPC/VDAF) and from the S5PVT scientific project CHEOPS-5p. Based on the same validation best practices as developed for and applied to heritage sensors like GOME-2, OMI and IASI (Keppens et al., 2015, 2018), the validation methodology relies on the analysis of data retrieval diagnostics – like the averaging kernels’ information content – and on comparisons of TROPOMI data with reference ozone profile measurements. The latter are acquired by ozonesonde, stratospheric lidar, and tropospheric lidar stations performing network operation in the context of WMO's Global Atmosphere Watch and its contributing networks NDACC and SHADOZ. The dependence of TROPOMI’s ozone profile uncertainty on several influence quantities like cloud fraction and measurement parameters like sun and scan angles is examined and discussed. This work concludes with a set of quality indicators enabling users to verify the fitness-for-purpose of the S5P data.</p>


2018 ◽  
Author(s):  
Ghazal Farhani ◽  
Robert J. Sica ◽  
Sophie Godin-Beekmann ◽  
Alexander Haefele

Abstract. This paper provides a detailed description of the first principle Optimal Estimation Method (OEM) which is applied to ozone retrieval analysis using Differential Absorption Lidar (DIAL) measurements. The air density, detector dead times, background coefficients, and lidar constants are simultaneously retrieved along with ozone density profiles. Using an averaging kernel, the OEM provides the vertical resolution of the retrieval as a function of altitude. A maximum acceptable height at which the a priori has a small contribution to the retrieval is calculated for each profile as well. Moreover, a complete uncertainty budget including both systematic and statistical uncertainties is given for each individual retrieved profile. Long term stratospheric DIAL ozone measurements have been carried out at the Observatoire de Haute-Provence (OHP) since 1985. The OEM is applied to 3 nights of measurements at OHP during an intensive ozone campaign in July 2017 where coincident lidar-ozonesonde measurements are available. The retrieved ozone density profiles are in good agreement with both traditional analysis and the ozonesonde measurements. For the three nights of measurements, below 15 km the difference between the OEM and the sonde profiles is less than 25 %, at altitudes between 15 km to 25 km the difference is less than 10 %, and the OEM can successfully catch many variations of ozone which are detected in the sonde profiles due to its ability to adjust its vertical resolution as the signal varies. Above 25 km the difference between the OEM and the sonde profiles does not exceed 20 %.


2015 ◽  
Vol 15 (12) ◽  
pp. 16615-16654 ◽  
Author(s):  
U. Jeong ◽  
J. Kim ◽  
C. Ahn ◽  
O. Torres ◽  
X. Liu ◽  
...  

Abstract. An online version of the OMI (Ozone Monitoring Instrument) near-ultraviolet (UV) aerosol retrieval algorithm was developed to retrieve aerosol optical thickness (AOT) and single scattering albedo (SSA) based on the optimal estimation (OE) method. Instead of using the traditional look-up tables for radiative transfer calculations, it performs online radiative transfer calculations with the Vector Linearized Discrete Ordinate Radiative Transfer (VLIDORT) model to eliminate interpolation errors and improve stability. The OE-based algorithm has the merit of providing useful estimates of uncertainties simultaneously with the inversion products. The measurements and inversion products of the Distributed Regional Aerosol Gridded Observation Network campaign in Northeast Asia (DRAGON NE-Asia 2012) were used to validate the retrieved AOT and SSA. The retrieved AOT and SSA at 388 nm have a correlation with the Aerosol Robotic Network (AERONET) products that is comparable to or better than the correlation with the operational product during the campaign. The estimated retrieval noise and smoothing error perform well in representing the envelope curve of actual biases of AOT at 388 nm between the retrieved AOT and AERONET measurements. The forward model parameter errors were analyzed separately for both AOT and SSA retrievals. The surface albedo at 388 nm, the imaginary part of the refractive index at 354 nm, and the number fine mode fraction (FMF) were found to be the most important parameters affecting the retrieval accuracy of AOT, while FMF was the most important parameter for the SSA retrieval. The additional information provided with the retrievals, including the estimated error and degrees of freedom, is expected to be valuable for future studies.


2014 ◽  
Vol 7 (6) ◽  
pp. 1547-1570 ◽  
Author(s):  
C. Viatte ◽  
K. Strong ◽  
K. A. Walker ◽  
J. R. Drummond

Abstract. We present a five-year time series of seven tropospheric species measured using a ground-based Fourier transform infrared (FTIR) spectrometer at the Polar Environment Atmospheric Research Laboratory (PEARL; Eureka, Nunavut, Canada; 80°05' N, 86°42' W) from 2007 to 2011. Total columns and temporal variabilities of carbon monoxide (CO), hydrogen cyanide (HCN) and ethane (C2H6) as well as the first derived total columns at Eureka of acetylene (C2H2), methanol (CH3OH), formic acid (HCOOH) and formaldehyde (H2CO) are investigated, providing a new data set in the sparsely sampled high latitudes. Total columns are obtained using the SFIT2 retrieval algorithm based on the optimal estimation method. The microwindows as well as the a priori profiles and variabilities are selected to optimize the information content of the retrievals, and error analyses are performed for all seven species. Our retrievals show good sensitivities in the troposphere. The seasonal amplitudes of the time series, ranging from 34 to 104%, are captured while using a single a priori profile for each species. The time series of the CO, C2H6 and C2H2 total columns at PEARL exhibit strong seasonal cycles with maxima in winter and minima in summer, in opposite phase to the HCN, CH3OH, HCOOH and H2CO time series. These cycles result from the relative contributions of the photochemistry, oxidation and transport as well as biogenic and biomass burning emissions. Comparisons of the FTIR partial columns with coincident satellite measurements by the Atmospheric Chemistry Experiment Fourier Transform Spectrometer (ACE-FTS) show good agreement. The correlation coefficients and the slopes range from 0.56 to 0.97 and 0.50 to 3.35, respectively, for the seven target species. Our new data set is compared to previous measurements found in the literature to assess atmospheric budgets of these tropospheric species in the high Arctic. The CO and C2H6concentrations are consistent with negative trends observed over the Northern Hemisphere, attributed to fossil fuel emission decrease. The importance of poleward transport for the atmospheric budgets of HCN and C2H2 is highlighted. Columns and variabilities of CH3OH and HCOOH at PEARL are comparable to previous measurements performed at other remote sites. However, the small columns of H2CO in early May might reflect its large atmospheric variability and/or the effect of the updated spectroscopic parameters used in our retrievals. Overall, emissions from biomass burning contribute to the day-to-day variabilities of the seven tropospheric species observed at Eureka.


2019 ◽  
Vol 12 (7) ◽  
pp. 3699-3716 ◽  
Author(s):  
Shannon Hicks-Jalali ◽  
Robert J. Sica ◽  
Alexander Haefele ◽  
Giovanni Martucci

Abstract. Raman lidars have been designated as potential candidates for trend studies by the Network for the Detection of Atmospheric Composition Change (NDACC) and GCOS (Global Climate Observing System) Reference Upper Air Network (GRUAN); however, for such studies improved calibration techniques are needed as well as careful consideration of the calibration uncertainties. Trend determinations require frequent, accurate, and well-characterized measurements. However, water vapour Raman lidars produce a relative measurement and require calibration in order to transform the measurement into a mixing ratio, a conserved quantity when no sources or sinks for water vapour are present. Typically, the calibration is done using a reference instrument such as a radiosonde. We present an improved trajectory technique to calibrate water vapour Raman lidars based on the previous work of Whiteman et al. (2006), Leblanc and Mcdermid (2008), Adam et al. (2010), and Herold et al. (2011), who used radiosondes as an external calibration source and matched the lidar measurements to the corresponding radiosonde measurement. However, they did not consider the movement of the radiosonde relative to the air mass and fronts. Our trajectory method is a general technique which may be used for any lidar and only requires that the radiosonde report wind speed and direction. As calibrations can be affected by a lack of co-location with the reference instrument, we have attempted to improve their technique by tracking the air parcels measured by the radiosonde relative to the field of view of the lidar. This study uses GRUAN Vaisala RS92 radiosonde measurements and lidar measurements taken by the MeteoSwiss RAman Lidar for Meteorological Observation (RALMO), located in Payerne, Switzerland, from 2011 to 2016 to demonstrate this improved calibration technique. We compare this technique to the traditional radiosonde–lidar calibration technique which does not involve tracking the radiosonde and uses the same integration time for all altitudes. Both traditional and our trajectory methods produce similar profiles when the water vapour field is homogeneous over the 30 min calibration period. We show that the trajectory method reduces differences between the radiosonde and lidar by an average of 10 % when the water vapour field is not homogeneous over a 30 min calibration period. We also calculate a calibration uncertainty budget that can be performed on a nightly basis. The calibration uncertainty budget includes the uncertainties due to phototube paralysis, aerosol extinctions, the assumption of the Ångström exponent, and the radiosonde. The study showed that the radiosonde was the major source of uncertainty in the calibration at 4 % of the calibration value. This trajectory method showed small improvements for RALMO's calibration but would be more useful for stations in different climatological regions or when non-co-located radiosondes are the only available calibration source.


Sign in / Sign up

Export Citation Format

Share Document