scholarly journals Error sources in the retrieval of aerosol information over bright surfaces from satellite measurements in the oxygen A band

2018 ◽  
Vol 11 (1) ◽  
pp. 161-175 ◽  
Author(s):  
Swadhin Nanda ◽  
Martin de Graaf ◽  
Maarten Sneep ◽  
Johan F. de Haan ◽  
Piet Stammes ◽  
...  

Abstract. Retrieving aerosol optical thickness and aerosol layer height over a bright surface from measured top-of-atmosphere reflectance spectrum in the oxygen A band is known to be challenging, often resulting in large errors. In certain atmospheric conditions and viewing geometries, a loss of sensitivity to aerosol optical thickness has been reported in the literature. This loss of sensitivity has been attributed to a phenomenon known as critical surface albedo regime, which is a range of surface albedos for which the top-of-atmosphere reflectance has minimal sensitivity to aerosol optical thickness. This paper extends the concept of critical surface albedo for aerosol layer height retrievals in the oxygen A band, and discusses its implications. The underlying physics are introduced by analysing the top-of-atmosphere reflectance spectrum as a sum of atmospheric path contribution and surface contribution, obtained using a radiative transfer model. Furthermore, error analysis of an aerosol layer height retrieval algorithm is conducted over dark and bright surfaces to show the dependence on surface reflectance. The analysis shows that the derivative with respect to aerosol layer height of the atmospheric path contribution to the top-of-atmosphere reflectance is opposite in sign to that of the surface contribution – an increase in surface brightness results in a decrease in information content. In the case of aerosol optical thickness, these derivatives are anti-correlated, leading to large retrieval errors in high surface albedo regimes. The consequence of this anti-correlation is demonstrated with measured spectra in the oxygen A band from the GOME-2 instrument on board the Metop-A satellite over the 2010 Russian wildfires incident.

2017 ◽  
Author(s):  
Swadhin Nanda ◽  
Martin de Graaf ◽  
Maarten Sneep ◽  
Johan F. de Haan ◽  
Piet Stammes ◽  
...  

Abstract. Retrieving aerosol optical thickness and aerosol layer height over a bright surface from measured top of atmosphere reflectance spectrum in the oxygen A-band is known to be challenging, often resulting in large errors. In certain atmospheric conditions and viewing geometries, a loss of sensitivity to aerosol optical thickness has been reported in literature. This loss of sensitivity has been attributed to a phenomenon known as critical surface albedo regime, which is a range of surface albedos for which the top of atmosphere reflectance has minimal sensitivity to aerosol optical thickness. This paper extends the concept of critical surface albedo for aerosol layer height retrievals in the oxygen A-band, and discusses its implications. The underlying physics are introduced by analysing top of atmosphere reflectance spectra obtained using a radiative transfer model. Furthermore, error analysis of the aerosol layer height retrieval algorithm are conducted over dark and bright surfaces to show the dependency on surface reflectance. The analysis shows that the information on aerosol layer height from atmospheric path contribution and the surface contribution to the top of atmosphere are opposite in sign – an increase in surface brightness results in a decrease in information content. In the case of aerosol optical thickness, these contributions are anti-correlated, leading to large retrieval errors in high surface albedo regimes. The consequence of this anti-correlation is demonstrated with measured spectra in the oxygen A-band from GOME-2A instrument on board the Metop-A satellite over the 2010 Russian wildfires incident.


2016 ◽  
Author(s):  
Julien Chimot ◽  
Joris Pepijn Veefkind ◽  
Tim Vlemmix ◽  
Johan de Haan ◽  
Vassilis Amiridis ◽  
...  

Abstract. This paper presents an exploratory study on the retrieval of aerosol layer height (ALH) from the OMI 477 nm O2–O2 spectral band. We have developed algorithms based on the Multilayer Perceptron (MLP) Neural Network (NN) approach and applied them on 3-year (2005–2007) OMI cloud-free scenes over North-East Asia, collocated with MODIS-Aqua aerosol product. In addition to the importance of aerosol altitude for climate and air quality objectives, the main motivation of this study is to evaluate the possibility of retrieving ALH for potential future improvements of trace gas retrievals (e.g. NO2, HCHO, SO2, etc..) from UV-Vis air quality satellite measurements over scenes including high aerosol concentrations. ALH retrieval relies on the analysis of the O2–O2 slant column density (SCD) and requires an accurate knowledge of the aerosol optical thickness τ. Using the MODIS-Aqua aerosol optical thickness at 550 nm as a prior information, comparison with the LIdar climatology of vertical Aerosol Structure for space-based lidar simulation (LIVAS) shows that ALH average biases over scenes with MODIS τ ≥ are in the range of 260–800 m. These results depend on the assumed aerosol single scattering albedo (sensitivity up to 600 m) and the chosen surface albedo (variation less than 200 m). Scenes with τ ≤ 0.5 are expected to show too large biases due to the little impacts of particles on the O2–O2 SCD changes. In addition, NN algorithms also enable aerosol optical thickness retrieval by exploring the OMI reflectance in the continuum. Comparisons with collocated MODIS-Aqua show agreements between −0.02 ± 0.45 and −0.18 ± 0.24 depending on the season. Improvements may be obtained from a better knowledge of the surface albedo, and higher accuracy of the aerosol model. This study shows the first encouraging aerosol layer height retrieval results over land from satellite observations of the 477 nm O2–O2 spectral band.


2019 ◽  
Vol 12 (12) ◽  
pp. 6619-6634 ◽  
Author(s):  
Swadhin Nanda ◽  
Martin de Graaf ◽  
J. Pepijn Veefkind ◽  
Mark ter Linden ◽  
Maarten Sneep ◽  
...  

Abstract. To retrieve aerosol properties from satellite measurements of the oxygen A-band in the near-infrared, a line-by-line radiative transfer model implementation requires a large number of calculations. These calculations severely restrict a retrieval algorithm's operational capability as it can take several minutes to retrieve the aerosol layer height for a single ground pixel. This paper proposes a forward modelling approach using artificial neural networks to speed up the retrieval algorithm. The forward model outputs are trained into a set of neural network models to completely replace line-by-line calculations in the operational processor. Results comparing the forward model to the neural network alternative show an encouraging outcome with good agreement between the two when they are applied to retrieval scenarios using both synthetic and real measured spectra from TROPOMI (TROPOspheric Monitoring Instrument) on board the European Space Agency (ESA) Sentinel-5 Precursor mission. With an enhancement of the computational speed by 3 orders of magnitude, TROPOMI's operational aerosol layer height processor is now able to retrieve aerosol layer heights well within operational capacity.


2015 ◽  
Vol 8 (11) ◽  
pp. 4947-4977 ◽  
Author(s):  
A. F. J. Sanders ◽  
J. F. de Haan ◽  
M. Sneep ◽  
A. Apituley ◽  
P. Stammes ◽  
...  

Abstract. An algorithm setup for the operational Aerosol Layer Height product for TROPOMI on the Sentinel-5 Precursor mission is described and discussed, applied to GOME-2A data, and evaluated with lidar measurements. The algorithm makes a spectral fit of reflectance at the O2 A band in the near-infrared and the fit window runs from 758 to 770 nm. The aerosol profile is parameterised by a scattering layer with constant aerosol volume extinction coefficient and aerosol single scattering albedo and with a fixed pressure thickness. The algorithm's target parameter is the height of this layer. In this paper, we apply the algorithm to observations from GOME-2A in a number of systematic and extensive case studies, and we compare retrieved aerosol layer heights with lidar measurements. Aerosol scenes cover various aerosol types, both elevated and boundary layer aerosols, and land and sea surfaces. The aerosol optical thicknesses for these scenes are relatively moderate. Retrieval experiments with GOME-2A spectra are used to investigate various sensitivities, in which particular attention is given to the role of the surface albedo. From retrieval simulations with the single-layer model, we learn that the surface albedo should be a fit parameter when retrieving aerosol layer height from the O2 A band. Current uncertainties in surface albedo climatologies cause biases and non-convergences when the surface albedo is fixed in the retrieval. Biases disappear and convergence improves when the surface albedo is fitted, while precision of retrieved aerosol layer pressure is still largely within requirement levels. Moreover, we show that fitting the surface albedo helps to ameliorate biases in retrieved aerosol layer height when the assumed aerosol model is inaccurate. Subsequent retrievals with GOME-2A spectra confirm that convergence is better when the surface albedo is retrieved simultaneously with aerosol parameters. However, retrieved aerosol layer pressures are systematically low (i.e., layer high in the atmosphere) to the extent that retrieved values no longer realistically represent actual extinction profiles. When the surface albedo is fixed in retrievals with GOME-2A spectra, convergence deteriorates as expected, but retrieved aerosol layer pressures become much higher (i.e., layer lower in atmosphere). The comparison with lidar measurements indicates that retrieved aerosol layer heights are indeed representative of the underlying profile in that case. Finally, subsequent retrieval simulations with two-layer aerosol profiles show that a model error in the assumed profile (two layers in the simulation but only one in the retrieval) is partly absorbed by the surface albedo when this parameter is fitted. This is expected in view of the correlations between errors in fit parameters and the effect is relatively small for elevated layers (less than 100 hPa). If one of the scattering layers is near the surface (boundary layer aerosols), the effect becomes surprisingly large, in such a way that the retrieved height of the single layer is above the two-layer profile. Furthermore, we find that the retrieval solution, once retrieval converges, hardly depends on the starting values for the fit. Sensitivity experiments with GOME-2A spectra also show that aerosol layer height is indeed relatively robust against inaccuracies in the assumed aerosol model, even when the surface albedo is not fitted. We show spectral fit residuals, which can be used for further investigations. Fit residuals may be partly explained by spectroscopic uncertainties, which is suggested by an experiment showing the improvement of convergence when the absorption cross section is scaled in agreement with Butz et al. (2013) and Crisp et al. (2012), and a temperature offset to the a priori ECMWF temperature profile is fitted. Retrieved temperature offsets are always negative and quite large (ranging between −4 and −8 K), which is not expected if temperature offsets absorb remaining inaccuracies in meteorological data. Other sensitivity experiments investigate fitting of stray light and fluorescence emissions. We find negative radiance offsets and negative fluorescence emissions, also for non-vegetated areas, but from the results it is not clear whether fitting these parameters improves the retrieval. Based on the present results, the operational baseline for the Aerosol Layer Height product currently will not fit the surface albedo. The product will be particularly suited for elevated, optically thick aerosol layers. In addition to its scientific value in climate research, anticipated applications of the product for TROPOMI are providing aerosol height information for aviation safety and improving interpretation of the Absorbing Aerosol Index.


2015 ◽  
Vol 8 (6) ◽  
pp. 6045-6118 ◽  
Author(s):  
A. F. J. Sanders ◽  
J. F. de Haan ◽  
M. Sneep ◽  
A. Apituley ◽  
P. Stammes ◽  
...  

Abstract. An algorithm setup for the operational Aerosol Layer Height product for TROPOMI on the Sentinel-5 Precursor mission is described and discussed, applied to GOME-2A data, and evaluated with lidar measurements. The algorithm makes a spectral fit of reflectance at the O2 A band in the near-infrared and the fit window runs from 758 to 770 nm. The aerosol profile is parameterized by a scattering layer with constant aerosol volume extinction coefficient and aerosol single scattering albedo and with a fixed pressure thickness. The algorithm's target parameter is the height of this layer. In this paper, we apply the algorithm to observations from GOME-2A in a number of systematic and extensive case studies and we compare retrieved aerosol layer heights with lidar measurements. Aerosol scenes cover various aerosol types, both elevated and boundary layer aerosols, and land and sea surfaces. The aerosol optical thicknesses for these scenes are relatively moderate. Retrieval experiments with GOME-2A spectra are used to investigate various sensitivities, in which particular attention is given to the role of the surface albedo. From retrieval simulations with the single-layer model, we learn that the surface albedo should be a fit parameter when retrieving aerosol layer height from the O2 A band. Current uncertainties in surface albedo climatologies cause biases and non-convergences when the surface albedo is fixed in the retrieval. Biases disappear and convergence improves when the surface albedo is fitted, while precision of retrieved aerosol layer pressure is still largely within requirement levels. Moreover, we show that fitting the surface albedo helps to ameliorate biases in retrieved aerosol layer height when the assumed aerosol model is inaccurate. Subsequent retrievals with GOME-2A spectra confirm that convergence is better when the surface albedo is retrieved simultaneously with aerosol parameters. However, retrieved aerosol layer pressures are systematically low (i.e., layer high in the atmosphere) to the extent that retrieved values are not realistically representing actual extinction profiles anymore. When the surface albedo is fixed in retrievals with GOME-2A spectra, convergence deteriorates as expected, but retrieved aerosol layer pressures become much higher (i.e., layer lower in atmosphere). The comparison with lidar measurements indicates that retrieved aerosol layer heights are indeed representative of the underlying profile in that case. Finally, subsequent retrieval simulations with two-layer aerosol profiles show that a model error in the assumed profile (two layers in the simulation but only one in the retrieval) is partly absorbed by the surface albedo when this parameter is fitted. This is expected in view of the correlations between errors in fit parameters and the effect is relatively small for elevated layers (less than 100 hPa). In case one of the scattering layers is near the surface (boundary layer aerosols), the effect becomes surprisingly large such that the retrieved height of the single layer is above the two-layer profile. Furthermore, we find that the retrieval solution, once retrieval converges, hardly depends on the starting values for the fit. Sensitivity experiments with GOME-2A spectra also show that aerosol layer height is indeed relatively robust against inaccuracies in the assumed aerosol model, even when the surface albedo is not fitted. We show spectral fit residuals, which can be used for further investigations. Fit residuals may be partly explained by spectroscopic uncertainties, which is suggested by an experiment showing the improvement of convergence when the absorption cross section is scaled in agreement with Butz et al. (2012) and Crisp et al. (2012) and a temperature offset to the a priori ECMWF temperature profile is fitted. Retrieved temperature offsets are always negative and quite large (ranging between −4 and −8 K), which is not expected if temperature offsets absorb remaining inaccuracies in meteorological data. Other sensitivity experiments investigate fitting of stray light and fluorescence emissions. We find negative radiance offsets and negative fluorescence emissions, also for non-vegetated areas, but from the results it is not clear whether fitting these parameters improves the retrieval. Based on the present results, the operational baseline for the Aerosol Layer Height product currently will not fit the surface albedo. The product will be particularly suited for elevated, optically thick aerosol layers. In addition to its scientific value in climate research, anticipated applications of the product for TROPOMI are providing aerosol height information for aviation safety and improving interpretation of the Absorbing Aerosol Index.


2019 ◽  
Author(s):  
Swadhin Nanda ◽  
Martin de Graaf ◽  
J. Pepijn Veefkind ◽  
Mark ter Linden ◽  
Maarten Sneep ◽  
...  

Abstract. To retrieve aerosol properties from satellite measurements of the oxygen A-band in the near infrared, a line-by-line radiative transfer model implementation requires a large number of calculations. These calculations severely restrict a retrieval algorithm's operational capability as it can take several minutes to retrieve aerosol layer height for a single ground pixel. This paper proposes a forward modeling approach using artificial neural networks to speed up the retrieval algorithm. The forward model outputs are trained into a set of neural network models to completely replace line-by-line calculations in the operational processor. Results of comparing the forward model to the neural network alternative show encouraging results with good agreements between the two when applied to retrieval scenarios using both synthetic and real measured spectra from TROPOMI (TROPOspheric Monitoring Instrument) on board the ESA Sentinel-5 Precursor mission. With an enhancement of the computational speed by three orders of magnitude, TROPOMI's operational aerosol layer height processor is now able to retrieve aerosol layer heights well within operational capacity.


2017 ◽  
Vol 10 (3) ◽  
pp. 783-809 ◽  
Author(s):  
Julien Chimot ◽  
J. Pepijn Veefkind ◽  
Tim Vlemmix ◽  
Johan F. de Haan ◽  
Vassilis Amiridis ◽  
...  

Abstract. This paper presents an exploratory study on the aerosol layer height (ALH) retrieval from the OMI 477 nm O2 − O2 spectral band. We have developed algorithms based on the multilayer perceptron (MLP) neural network (NN) approach and applied them to 3-year (2005–2007) OMI cloud-free scenes over north-east Asia, collocated with MODIS Aqua aerosol product. In addition to the importance of aerosol altitude for climate and air quality objectives, our long-term motivation is to evaluate the possibility of retrieving ALH for potential future improvements of trace gas retrievals (e.g. NO2, HCHO, SO2) from UV–visible air quality satellite measurements over scenes including high aerosol concentrations. This study presents a first step of this long-term objective and evaluates, from a statistic point of view, an ensemble of OMI ALH retrievals over a long time period of 3 years covering a large industrialized continental region. This ALH retrieval relies on the analysis of the O2 − O2 slant column density (SCD) and requires an accurate knowledge of the aerosol optical thickness, τ. Using MODIS Aqua τ(550 nm) as a prior information, absolute seasonal differences between the LIdar climatology of vertical Aerosol Structure for space-based lidar simulation (LIVAS) and average OMI ALH, over scenes with MODIS τ(550 nm) ≥ 1. 0, are in the range of 260–800 m (assuming single scattering albedo ω0 = 0. 95) and 180–310 m (assuming ω0 = 0. 9). OMI ALH retrievals depend on the assumed aerosol single scattering albedo (sensitivity up to 660 m) and the chosen surface albedo (variation less than 200 m between OMLER and MODIS black-sky albedo). Scenes with τ ≤ 0. 5 are expected to show too large biases due to the little impact of particles on the O2 − O2 SCD changes. In addition, NN algorithms also enable aerosol optical thickness retrieval by exploring the OMI reflectance in the continuum. Comparisons with collocated MODIS Aqua show agreements between −0. 02  ±  0. 45 and −0. 18  ±  0. 24, depending on the season. Improvements may be obtained from a better knowledge of the surface albedo and higher accuracy of the aerosol model. Following the previous work over ocean of Park et al.(2016), our study shows the first encouraging aerosol layer height retrieval results over land from satellite observations of the 477 nm O2 − O2 absorption spectral band.


2021 ◽  
Author(s):  
Marta Luffarelli ◽  
Yves Govaerts

<p>The CISAR (Combined Inversion of Surface and AeRosols) algorithm is exploited in the framework of the ESA Aerosol Climate Change Initiatiave (CCI) project, aiming at providing a set of atmospheric (cloud and aerosol) and surface reflectance products derived from S3A/SLSTR observations using the same radiative transfer physics and assumptions. CISAR is an advance algorithm developed by Rayference originally designed for the retrieval of aerosol single scattering properties and surface reflectance from both geostationary and polar orbiting satellite observations.  It is based on the inversion of a fast radiative transfer model (FASTRE). The retrieval mechanism allows a continuous variation of the aerosol and cloud single scattering properties in the solution space.</p><p> </p><p>Traditionally, different approaches are exploited to retrieve the different Earth system components, which could lead to inconsistent data sets. The simultaneous retrieval of different atmospheric and surface variables over any type of surface (including bright surfaces and water bodies) with the same forward model and inversion scheme ensures the consistency among the retrieved Earth system components. Additionally, pixels located in the transition zone between pure clouds and pure aerosols are often discarded from both cloud and aerosol algorithms. This “twilight zone” can cover up to 30% of the globe. A consistent retrieval of both cloud and aerosol single scattering properties with the same algorithm could help filling this gap.</p><p> </p><p>The CISAR algorithm aims at overcoming the need of an external cloud mask, discriminating internally between aerosol and cloud properties. This approach helps reducing the overestimation of aerosol optical thickness in cloud contaminated pixels. The surface reflectance product is delivered both for cloud-free and cloudy observations.  </p><p> </p><p>Global maps obtained from the processing of S3A/SLSTR observations will be shown. The SLSTR/CISAR products over events such as, for instance, the Australian fire in the last months of 2019, will be discussed in terms of aerosol optical thickness, aerosol-cloud discrimination and fine/coarse mode fraction.</p>


2018 ◽  
Vol 11 (1) ◽  
pp. 499-514 ◽  
Author(s):  
Travis D. Toth ◽  
James R. Campbell ◽  
Jeffrey S. Reid ◽  
Jason L. Tackett ◽  
Mark A. Vaughan ◽  
...  

Abstract. Due to instrument sensitivities and algorithm detection limits, level 2 (L2) Cloud-Aerosol Lidar with Orthogonal Polarization (CALIOP) 532 nm aerosol extinction profile retrievals are often populated with retrieval fill values (RFVs), which indicate the absence of detectable levels of aerosol within the profile. In this study, using 4 years (2007–2008 and 2010–2011) of CALIOP version 3 L2 aerosol data, the occurrence frequency of daytime CALIOP profiles containing all RFVs (all-RFV profiles) is studied. In the CALIOP data products, the aerosol optical thickness (AOT) of any all-RFV profile is reported as being zero, which may introduce a bias in CALIOP-based AOT climatologies. For this study, we derive revised estimates of AOT for all-RFV profiles using collocated Moderate Resolution Imaging Spectroradiometer (MODIS) Dark Target (DT) and, where available, AErosol RObotic NEtwork (AERONET) data. Globally, all-RFV profiles comprise roughly 71 % of all daytime CALIOP L2 aerosol profiles (i.e., including completely attenuated profiles), accounting for nearly half (45 %) of all daytime cloud-free L2 aerosol profiles. The mean collocated MODIS DT (AERONET) 550 nm AOT is found to be near 0.06 (0.08) for CALIOP all-RFV profiles. We further estimate a global mean aerosol extinction profile, a so-called “noise floor”, for CALIOP all-RFV profiles. The global mean CALIOP AOT is then recomputed by replacing RFV values with the derived noise-floor values for both all-RFV and non-all-RFV profiles. This process yields an improvement in the agreement of CALIOP and MODIS over-ocean AOT.


Sign in / Sign up

Export Citation Format

Share Document