Recent years have seen the increasing inclusion of per-retrieval prognostic (predictive) uncertainty estimates within satellite aerosol optical depth (AOD) data sets, providing users with quantitative tools to assist in the optimal use of these data. Prognostic estimates contrast with diagnostic (i.e. relative to some external truth) ones, which are typically obtained using sensitivity and/or validation analyses. Up to now, however, the quality of these uncertainty estimates has not been routinely assessed. This study presents a review of existing prognostic and diagnostic approaches for quantifying uncertainty in satellite AOD retrievals, and it presents a general framework to evaluate them based on the expected statistical properties of ensembles of estimated uncertainties and actual retrieval errors. It is hoped that this framework will be adopted as a complement to existing AOD validation exercises; it is not restricted to AOD and can in principle be applied to other quantities for which a reference validation data set is available. This framework is then applied to assess the uncertainties provided by several satellite data sets (seven over land, five over water), which draw on methods from the empirical to sensitivity analyses to formal error propagation, at 12 Aerosol Robotic Network (AERONET) sites. The AERONET sites are divided into those for which it is expected that the techniques will perform well and those for which some complexity about the site may provide a more severe test. Overall, all techniques show some skill in that larger estimated uncertainties are generally associated with larger observed errors, although they are sometimes poorly calibrated (i.e. too small or too large in magnitude). No technique uniformly performs best. For powerful formal uncertainty propagation approaches such as optimal estimation, the results illustrate some of the difficulties in appropriate population of the covariance matrices required by the technique. When the data sets are confronted by a situation strongly counter to the retrieval forward model (e.g. potentially mixed land–water surfaces or aerosol optical properties outside the family of assumptions), some algorithms fail to provide a retrieval, while others do but with a quantitatively unreliable uncertainty estimate. The discussion suggests paths forward for the refinement of these techniques.

Characterizing changes in landscape fire activity at better than hourly temporal resolution is achievable using thermal observations of actively burning fires made from geostationary Earth Observation (EO) satellites. Over the last decade or more, a series of research and/or operational "active fire" products have been developed from geostationary EO data, often with the aim of supporting biomass burning fuel consumption and trace gas and aerosol emission calculations. Such Fire Radiative Power (FRP) products are generated operationally from Meteosat by the Land Surface Analysis Satellite Applications Facility (LSA SAF) and are available freely every 15 min in both near-real-time and archived form. These products map the location of actively burning fires and characterize their rates of thermal radiative energy release (FRP), which is believed proportional to rates of biomass consumption and smoke emission. The FRP-PIXEL product contains the full spatio-temporal resolution FRP data set derivable from the SEVIRI (Spinning Enhanced Visible and Infrared Imager) imager onboard Meteosat at a 3 km spatial sampling distance (decreasing away from the west African sub-satellite point), whilst the FRP-GRID product is an hourly summary at 5° grid resolution that includes simple bias adjustments for meteorological cloud cover and regional underestimation of FRP caused primarily by underdetection of low FRP fires. Here we describe the enhanced geostationary Fire Thermal Anomaly (FTA) detection algorithm used to deliver these products and detail the methods used to generate the atmospherically corrected FRP and per-pixel uncertainty metrics. Using SEVIRI scene simulations and real SEVIRI data, including from a period of Meteosat-8 "special operations", we describe certain sensor and data pre-processing characteristics that influence SEVIRI's active fire detection and FRP measurement capability, and use these to specify parameters in the FTA algorithm and to make recommendations for the forthcoming Meteosat Third Generation operations in relation to active fire measures. We show that the current SEVIRI FTA algorithm is able to discriminate actively burning fires covering down to 10⁻⁴ of a pixel and that it appears more sensitive to fire than other algorithms used to generate many widely exploited active fire products. Finally, we briefly illustrate the information contained within the current Meteosat FRP-PIXEL and FRP-GRID products, providing example analyses for both individual fires and multi-year regional-scale fire activity; the companion paper (Roberts et al., 2015) provides a full product performance evaluation and a demonstration of product use within components of the Copernicus Atmosphere Monitoring Service (CAMS).

This paper presents the simultaneous retrieval of aerosol optical thickness and surface properties from the CISAR algorithm applied both to geostationary and polar-orbiting satellite observations. The theoretical concepts of the CISAR algorithm have been described in Govaerts and Luffarelli (2018). CISAR has been applied to SEVIRI and PROBA-V observations acquired over 20 AERONET stations during the year 2015. The CISAR retrieval from the two sets of observations is evaluated against independent data sets such as the MODIS land product and AERONET data. The performance differences resulting from the two types of orbit are discussed, and the information content of SEVIRI and PROBA-V observations is analysed and compared.

The Libya-4 desert area, located in the Great Sand Sea, is one of the most important bright desert CEOS pseudo-invariant calibration sites by its size and radiometric stability. This site is intensively used for radiometer drift monitoring, sensor intercalibration and as an absolute calibration reference based on simulated radiances traceable to the SI standard. The Libya-4 morphology is composed of oriented sand dunes shaped by dominant winds. The effects of sand dune spatial organization on the surface bidirectional reflectance factor is analyzed in this paper using Raytran, a 3D radiative transfer model. The topography is characterized with the 30 m resolution ASTER digital elevation model. Four different regions-of-interest sizes, ranging from 10 km up to 100 km, are analyzed. Results show that sand dunes generate more backscattering than forward scattering at the surface. The mean surface reflectance averaged over different viewing and illumination angles is pretty much independent of the size of the selected area, though the standard deviation differs. Sun azimuth position has an effect on the surface reflectance field, which is more pronounced for high Sun zenith angles. Such 3D azimuthal effects should be taken into account to decrease the simulated radiance uncertainty over Libya-4 below 3% for wavelengths larger than 600 nm.

This Letter suggests the substitution of the spectral response of the Meteosat-5 and -6 radiometer solar channel by the one of the Meteosat-7 radiometer. This replacement is motivated by an inaccurate spectral characterisation of the former two radiometers. The impact of this correction on the vicarious calibration of the solar channel is assessed. When the original sensor response is used, a difference of about 20% is found between calibration coefficients derived over sea and bright desert. This difference drops to about 1% when the sensor spectral response of the Meteosat-7 radiometer is used for the calibration of the Meteosat-5 and -6 radiometer solar channel.

This paper presents a new algorithm for the joint retrieval of surface reflectance and aerosol properties with continuous variations of the state variables in the solution space. This algorithm, named CISAR (Combined Inversion of Surface and AeRosol), relies on a simple atmospheric vertical structure composed of two layers and an underlying surface. Surface anisotropic reflectance effects are taken into account and radiatively coupled with atmospheric scattering. For this purpose, a fast radiative transfer model has been explicitly developed, which includes acceleration techniques to solve the radiative transfer equation and to calculate the Jacobians. The inversion is performed within an optimal estimation framework including prior information on the state variable magnitude and regularisation constraints on their spectral and temporal variability. In each processed wavelength, the algorithm retrieves the parameters of the surface reflectance model, the aerosol total column optical thickness and single-scattering properties. The CISAR algorithm functioning is illustrated with a series of simple experiments.