ALBERT

Description: The angular distribution of the light reflected by the Earth's surface influences top-of-atmosphere (TOA) reflectance values. This surface reflectance anisotropy has implications for UV/Vis satellite retrievals of albedo, clouds, and trace gases such as nitrogen dioxide (NO2). These retrievals routinely assume the surface to reflect light isotropically. Here we show that cloud fractions retrieved from GOME-2A and OMI with the FRESCO and OMCLDO2 algorithms have an east–west bias of 10 % to 50 %, which are highest over vegetation and forested areas, and that this bias originates from the assumption of isotropic surface reflection. To interpret the across-track bias with the DAK radiative transfer model, we implement the bidirectional reflectance distribution function (BRDF) from the Ross–Li semi-empirical model. Testing our implementation against state-of-the-art RTMs LIDORT and SCIATRAN, we find that simulated TOA reflectance generally agrees to within 1 %. We replace the assumption of isotropic surface reflection in the equations used to retrieve cloud fractions over forested scenes with scattering kernels and corresponding BRDF parameters from a daily high-resolution database derived from 16 years' worth of MODIS measurements. By doing this, the east–west bias in the simulated cloud fractions largely vanishes. We conclude that across-track biases in cloud fractions can be explained by cloud algorithms that do not adequately account for the effects of surface reflectance anisotropy. The implications for NO2 air mass factor (AMF) calculations are substantial. Under moderately polluted NO2 and backward-scattering conditions, clear-sky AMFs are up to 20 % higher and cloud radiance fractions up to 40 % lower if surface anisotropic reflection is accounted for. The combined effect of these changes is that NO2 total AMFs increase by up to 30 % for backward-scattering geometries (and decrease by up to 35 % for forward-scattering geometries), which is stronger than the effect of either contribution alone. In an unpolluted troposphere, surface BRDF effects on cloud fraction counteract (and largely cancel) the effect on the clear-sky AMF. Our results emphasise that surface reflectance anisotropy needs to be taken into account in a coherent manner for more realistic and accurate retrievals of clouds and NO2 from UV/Vis satellite sensors. These improvements will be beneficial for current sensors, in particular for the recently launched TROPOMI instrument with a high spatial resolution.

Description: The FRESCO (Fast Retrieval Scheme for Clouds from the Oxygen A band) algorithm is a simple, fast and robust algorithm used to retrieve cloud information in operational satellite data processing. It has been applied to GOME-1 (Global Ozone Monitoring Experiment), SCIAMACHY (Scanning Imaging Absorption Spectrometer for Atmospheric Chartography), GOME-2 and more recently to TROPOMI (Tropospheric Monitoring Instrument). FRESCO retrieves effective cloud fraction and cloud pressure from measurements in the oxygen A band around 761 nm. In this paper, we propose a new version of the algorithm, called FRESCO-B, which is based on measurements in the oxygen B band around 687 nm. Such a method is interesting for vegetated surfaces where the surface albedo is much lower in the B band than in the A band, which limits the ground contribution to the top-of-atmosphere reflectances. In this study we first perform retrieval simulations. These show that the retrieved cloud pressures from FRESCO-B and FRESCO differ only between −10 and +10 hPa, except for high, thin clouds over vegetation where the difference is larger (about +15 to +30 hPa), with FRESCO-B yielding higher pressure. Next, inter-comparison between FRESCO-B and FRESCO retrievals over 1 month of GOME-2B data reveals that the effective cloud fractions retrieved in the O2 A and B bands are very similar (mean difference of 0.003), while the cloud pressures show a mean difference of 11.5 hPa, with FRESCO-B retrieving higher pressures than FRESCO. This agrees with the simulations and is partly due to deeper photon penetrations of the O2 B band in clouds compared to the O2 A-band photons and partly due to the surface albedo bias in FRESCO. Finally, validation with ground-based measurements shows that the FRESCO-B cloud pressure represents an altitude within the cloud boundaries for clouds that are not too far from the Lambertian reflector model, which occurs in about 50 % of the cases.

Description: In this work, a ∼21-year global dataset from four different satellite sensors with a mid-morning overpass (GOME/ERS-2, SCIAMACHY/ENVISAT, GOME-2/Metop-A, and GOME-2/Metop-B) is compiled to study the long-term tropospheric NO2 patterns and trends. The Global Ozone Monitoring Experiment (GOME) and GOME-2 data are “corrected” relative to the SCanning Imaging Absorption spectroMeter for Atmospheric CHartographY (SCIAMACHY) data to produce a self-consistent dataset that covers the period April 1996–September 2017. The highest tropospheric NO2 concentrations are seen over urban, industrialized, and highly populated areas and over ship tracks in the oceans. Tropospheric NO2 has generally decreased during the last 2 decades over the industrialized and highly populated regions of the western world (a total decrease of the order of ∼49 % over the US, the Netherlands, and the UK; ∼36 % over Italy and Japan; and ∼32 % over Germany and France) and increased over developing regions (a total increase of ∼160 % over China and ∼33 % over India). It is suggested here that linear trends cannot be used efficiently worldwide for such long periods. Tropospheric NO2 is very sensitive to socioeconomic changes (e.g., environmental protection policies, economic recession, warfare, etc.) which may cause either short-term changes or even a reversal of the trends. The application of a method capable of detecting the year when a reversal of trends happened shows that tropospheric NO2 concentrations switched from positive to negative trends and vice versa over several regions around the globe. A country-level analysis revealed clusters of countries that exhibit similar positive-to-negative or negative-to-positive trend reversals, while 29 out of a total of 64 examined megacities and large urban agglomerations experienced a trend reversal at some point within the last 2 decades.

Description: Ozone profile retrieval from nadir-viewing satellite instruments operating in the ultraviolet–visible range requires accurate calibration of Level-1 (L1) radiance data. Here we study the effects of calibration on the derived Level-2 (L2) ozone profiles for three versions of SCanning Imaging Absorption spectroMeter for Atmospheric ChartograpHY (SCIAMACHY) L1 data: version 7 (v7), version 7 with m-factors (v7mfac) and version 8 (v8). We retrieve nadir ozone profiles from the SCIAMACHY instrument that flew on board Envisat using the Ozone ProfilE Retrieval Algorithm (OPERA) developed at KNMI with a focus on stratospheric ozone. We study and assess the quality of these profiles and compare retrieved L2 products from L1 SCIAMACHY data versions from the years 2003 to 2011 without further radiometric correction. From validation of the profiles against ozone sonde measurements, we find that the v8 performs better than v7 and v7mfac due to correction for the scan-angle dependency of the instrument's optical degradation. Validation for the years 2003 and 2009 with ozone sondes shows deviations of SCIAMACHY ozone profiles of 0.8–15 % in the stratosphere (corresponding to pressure range ∼ 100–10 hPa) and 2.5–100 % in the troposphere (corresponding to pressure range ∼ 1000–100 hPa), depending on the latitude and the L1 version used. Using L1 v8 for the years 2003–2011 leads to deviations of ∼ 1–11 % in stratospheric ozone and ∼ 1–45 % in tropospheric ozone. The SCIAMACHY L1 v8 data can still be improved upon in the 265–330 nm range used for ozone profile retrieval. The slit function can be improved with a spectral shift and squeeze, which leads to a few percent residue reduction compared to reference solar irradiance spectra. Furthermore, studies of the ratio of measured to simulated reflectance spectra show that a bias correction in the reflectance for wavelengths below 300 nm appears to be necessary.

Description: Air mass factor (AMF) calculation is the largest source of uncertainty in NO2 and HCHO satellite retrievals in situations with enhanced trace gas concentrations in the lower troposphere. Structural uncertainty arises when different retrieval methodologies are applied within the scientific community to the same satellite observations. Here, we address the issue of AMF structural uncertainty via a detailed comparison of AMF calculation methods that are structurally different between seven retrieval groups for measurements from the Ozone Monitoring Instrument (OMI). We estimate the escalation of structural uncertainty in every sub-step of the AMF calculation process. This goes beyond the algorithm uncertainty estimates provided in state-of-the-art retrievals, which address the theoretical propagation of uncertainties for one particular retrieval algorithm only. We find that top-of-atmosphere reflectances simulated by four radiative transfer models (RTMs) (DAK, McArtim, SCIATRAN and VLIDORT) agree within 1.5 %. We find that different retrieval groups agree well in the calculations of altitude resolved AMFs from different RTMs (to within 3 %), and in the tropospheric AMFs (to within 6 %) as long as identical ancillary data (surface albedo, terrain height, cloud parameters and trace gas profile) and cloud and aerosol correction procedures are being used. Structural uncertainty increases sharply when retrieval groups use their preference for ancillary data, cloud and aerosol correction. On average, we estimate the AMF structural uncertainty to be 42 % over polluted regions and 31 % over unpolluted regions, mostly driven by substantial differences in the a priori trace gas profiles, surface albedo and cloud parameters. Sensitivity studies for one particular algorithm indicate that different cloud correction approaches result in substantial AMF differences in polluted conditions (5 to 40 % depending on cloud fraction and cloud pressure, and 11 % on average) even for low cloud fractions (

Description: Retrievals of cloud properties from geostationary satellite sensors offer extensive spatial and temporal coverage and resolution. The high temporal resolution allows the observation of diurnally resolved cloud properties. However, retrievals are sensitive to varying illumination and viewing geometries, including cloud glory and cloud bow conditions, which can lead to irregularities in the diurnal data record. In this study, these conditions and their effects on liquid cloud optical thickness and effective radius retrievals are analyzed using the Cloud Physical Properties (CPP) algorithm. This analysis is based on the use of Spinning Enhanced Visible and Infrared Imager (SEVIRI) reflectances and products from Meteosat-8 and Meteosat-10, which are located over the Indian and Atlantic Ocean, respectively, and cover an extensive common area under different viewing angles. Comparisons of the retrievals from two full days, over ocean and land, and using different spectral combinations of visible and shortwave-infrared channels, are performed, to assess the importance of these factors in the retrieval process. The sensitivity of the cloud-bow- and cloud-glory-related irregularities to the width of the assumed droplet size distribution is analyzed by using different values of the effective variance of the size distribution. The results suggest for marine stratocumulus clouds an effective variance of around 0.05, which implies a narrower size distribution than typically assumed in satellite-based retrievals. For the case with continental clouds a broader size distribution (effective variance around 0.15) is obtained. This highlights the importance of appropriate size distribution assumptions and provides a way to improve the quality of cloud products in future climate data record releases.

Description: 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.

Description: This paper presents a weighted least squares approach to retrieve aerosol layer height from top-of-atmosphere reflectance measurements in the oxygen A band (758–770 nm) over bright surfaces. A property of the measurement error covariance matrix is discussed, due to which photons travelling from the surface are given a higher preference over photons that scatter back from the aerosol layer. This is a potential source of biases in the estimation of aerosol properties over land, which can be mitigated by revisiting the design of the measurement error covariance matrix. The alternative proposed in this paper, which we call the dynamic scaling method, introduces a scene-dependent and wavelength-dependent modification in the measurement signal-to-noise ratio in order to influence this matrix. This method is generally applicable to other retrieval algorithms using weighted least squares. To test this method, synthetic experiments are done in addition to application to GOME-2A and GOME-2B measurements of the oxygen A band over the August 2010 Russian wildfires and the October 2017 Portugal wildfire plume over western Europe.

Description: The Ozone Monitoring Instrument (OMI) is a push-broom imaging spectrometer, observing solar radiation backscattered by the Earth's atmosphere and surface. The shape of an OMI pixel is not quadrangular, which is common for scanning instruments, but rather Gaussian-shaped as light from neighbouring pixels enters the Field of View (FoV). This has consequences for pixel-area dependent applications, like cloud fraction products, and visualisation. The shape and sizes of OMI pixels were determined pre-flight by theoretical and experimental tests, but never verified after launch. In this paper the OMI point spread function (PSF) is characterised using collocated MODerate resolution Imaging Spectroradiometer (MODIS) reflectance measurements. MODIS measurements have a much higher spatial resolution than OMI measurements and spectrally overlap at 469 nm. The optimal OMI PSF was determined by finding the highest correlation between MODIS and OMI reflectances for both cloud-free and partially clouded scenes. Our results show that the semi-official OMPIXCOR product 75FOV corner coordinates accurately fix the Full Width at Half Maximum (FWHM) of a super-Gaussian PSF, when this pixel shape is assumed. The exponent of the super-Gaussian PSF is dependent on OMI pixel row number, from about n = 2 at nadir to 3.5 at the swath edges, due to the increase in pixel size. The optimal Gaussian exponent depends on scene changes between overpasses and reduces to about n = 1 for partially clouded scenes before 2008. Then, the time difference between Aqua and Aura was 15 minutes instead of 8 minutes after 2008. Between overpasses, clouds change the scene, reducing the correlation and changing the shape of the optimal overlap.