Sea And Land Surface Temperature Radiometer Research Articles

The Operational and Climate Land Surface Temperature Products from the Sea and Land Surface Temperature Radiometers on Sentinel-3A and 3B

Land Surface Temperature (LST) is integral to our understanding of the radiative energy budget of the Earth’s surface since it provides the best approximation to the thermodynamic temperature that drives the outgoing longwave flux from surface to atmosphere. Since 5 July 2017, an operational LST product has been available from the Sentinel-3A mission, with the corresponding product being available from Sentinel-3B since 17 November 2018. Here, we present the first paper describing formal products, including algorithms, for the Sea and Land Surface Temperature Radiometer (SLSTR) instruments onboard Sentinel-3A and 3B (SLSTR-A and SLSTR-B, respectively). We evaluate the quality of both the Land Surface Temperature Climate Change Initiative (LST_cci) product and the Copernicus operational LST product (SL_2_LST) for the years 2018 to 2021. The evaluation takes the form of a validation against ground-based observations of LST across eleven well-established in situ stations. For the validation, the mean absolute daytime and night-time difference against the in situ measurements for the LST_cci product is 0.77 K and 0.50 K, respectively, for SLSTR-A, and 0.91 K and 0.54 K, respectively, for SLSTR-B. These are an improvement on the corresponding statistics for the SL_2_LST product, which are 1.45 K (daytime) and 0.76 (night-time) for SLSTR-A, and 1.29 K (daytime) and 0.77 (night-time) for SLSTR-B. The key influencing factors in this improvement include an upgraded database of reference states for the generation of retrieval coefficients, higher stratification of the auxiliary data for the biome and fractional vegetation, and enhanced cloud masking.

Estimates of Crop Yield Anomalies for 2022 in Ukraine Based on Copernicus Sentinel-1, Sentinel-3 Satellite Data, and ERA-5 Agrometeorological Indicators.

The study explores the feasibility of adapting the EOStat crop monitoring system, originally designed for monitoring crop growth conditions in Poland, to fulfill the requirements of a similar system in Ukraine. The system utilizes satellite data and agrometeorological information provided by the Copernicus program, which offers these resources free of charge. To predict crop yields, the system uses several factors, such as vegetation condition indices obtained from Sentinel-3 Ocean and Land Color Instrument (OLCI) optical and Sea and Land Surface Temperature Radiometer (SLSTR). It also incorporates climate information, including air temperature, total precipitation, surface radiation, and soil moisture. To identify the best predictors for each administrative unit, the study utilizes a recursive feature elimination method and employs the Extreme Gradient Boosting regressor, a machine learning algorithm, to forecast crop yields. The analysis indicates a noticeable decrease in crop losses in 2022 in certain regions of Ukraine, compared to the previous year (2021) and the 5-year average (2017-2021), specifically for winter crops and maize. Considering the reduction in yield, it is estimated that the decline in production of winter crops in 2022 was up to 20%, while for maize, it was up to 50% compared to the decline in production.

High-resolution (1 km) all-sky net radiation over Europe enabled by the merging of land surface temperature retrievals from geostationary and polar-orbiting satellites

Abstract. Surface net radiation (SNR) is a vital input for many land surface and hydrological models. However, most of the current remote sensing datasets of SNR come mostly at coarse resolutions or have large gaps due to cloud cover that hinder their use as input in models. Here, we present a downscaled and continuous daily SNR product across Europe for 2018–2019. Long-wave outgoing radiation is computed from a merged land surface temperature (LST) product in combination with Meteosat Second Generation emissivity data. The merged LST product is based on all-sky LST retrievals from the Spinning Enhanced Visible and InfraRed Imager (SEVIRI) onboard the geostationary Meteosat Second Generation (MSG) satellite and clear-sky LST retrievals from the Sea and Land Surface Temperature Radiometer (SLSTR) onboard the polar-orbiting Sentinel-3A satellite. This approach makes use of the medium spatial (approx. 5–7 km) but high temporal (30 min) resolution, gap-free data from MSG along with the low temporal (2–3 d) but high spatial (1 km) resolution of the Sentinel-3 LST retrievals. The resulting 1 km and daily LST dataset is based on an hourly merging of both datasets through bias correction and Kalman filter assimilation. Short-wave outgoing radiation is computed from the incoming short-wave radiation from MSG and the downscaled albedo using 1 km PROBA-V data. MSG incoming short-wave and long-wave radiation and the outgoing radiation components at 1 km spatial resolution are used together to compute the final daily SNR dataset in a consistent manner. Validation results indicate an improvement of the mean squared error by ca. 7 % with an increase in spatial detail compared to the original MSG product. The resulting pan-European SNR dataset, as well as the merged LST product, can be used for hydrological modelling and as input to models dedicated to estimating evaporation and surface turbulent heat fluxes and will be regularly updated in the future. The datasets can be downloaded from https://doi.org/10.5281/zenodo.8332222 (Rains, 2023a) and https://doi.org/10.5281/zenodo.8332128 (Rains, 2023b).

Retrieving land surface reflectance anisotropy with Sentinel-3 observations and prior BRDF model constraints

Sentinel-3 supports governmental policies and European programmes such as the Copernicus Global Land Service with daily monitoring of the land surface by measuring the Earth's surface with its multi-band wide-swath visible and near-infrared radiometers. To fully exploit Sentinel-3 data sets, performing a correction of directional effects caused by the land surface reflectance anisotropy properties is mandatory. This correction requires first an estimation of the bidirectional reflectance distribution function (BRDF) to be properly achieved. The quality and robustness of global estimates of key land surface variables derived from space-borne sensors depend clearly on an accurate assessment of the spectral BRDF.Here, we present an algorithm aiming at harnessing the Sentinel-3's wide imaging swath to retrieve the land surface BRDF, approximating it via a kernel-driven semi-empirical model which can then be used to normalise the observations to a common Sun-sensor geometry. Our algorithm, named ReBeLS (Regularised BRDF inversion for Land Surface), uses a temporally regularised BRDF inversion approach and, together with prior knowledge of the land surface BRDF obtained from the Moderate Resolution Imaging Spectroradiometer (MODIS) sensors, assimilates Sentinel-3 surface reflectance observations. We focus on the results derived from observations provided by the Sentinel-3 Ocean and Land Colour Instrument (OLCI) sensors; however, the algorithm can also be applied to reflectance observations acquired by the Sea and Land Surface Temperature Radiometer (SLSTR). The retrieved BRDF model parameters reproduce within uncertainties the directional signature as seen by the Visible Infrared Imaging Radiometer Suite (VIIRS) sensor. The generated Sentinel-3 normalised surface reflectance products are a significant advance compared with time series of directional surface reflectance and maximum value composite images, showing a statistically significant (p≪0.05) correlation (r>0.9) with VIIRS nadir BRDF-adjusted surface reflectance. Furthermore, the ReBeLS Sentinel-3 BRDF products offer enough flexibility to allow users to compute spectral surface albedos or to normalise surface reflectance to any Sun-sensor configuration.

Sentinel-3 SLSTR active fire (AF) detection and FRP daytime product - Algorithm description and global intercomparison to MODIS, VIIRS and landsat AF data

The Sea and Land Surface Temperature Radiometer (SLSTR) senses the Earth from onboard two concurrently operating European Copernicus Sentinel-3 (S3) satellites. As the Terra platform carrying the Moderate Resolution Imaging Spectroradiometer (MODIS) is reaching its end of life, the S3 Active Fire Detection and FRP products generated from data captured by S3 SLSTR are expected to soon become the main global active fire (AF) product for the mid-morning and evening low Earth orbit timeslots. The S3 night-time AF product issued by the European Space Agency (ESA) has been operational since March 2020, and here we report on the significant adjustments made to enable the generation of a complimentary daytime product. Similar to MODIS, SLSTR possesses two middle infrared channels, both a ‘standard’ (normal gain; S7) channel and a ‘fire’ (low-gain; F1) channel - but in contrast to MODIS by day even the ambient background land surface is often saturated in the SLSTR standard gain MIR (S7) channel. This saturation necessitates far greater use of the F1 channel data by day for active fire detection than by night, even though F1 has characteristics which make its data more challenging to combine with that from the other SLSTR thermal infrared channels. Here we report on the approaches used to combine S7 and F1 data for optimized daytime AF detection, and also detail the other algorithm adjustments found necessary to include in the daytime AF product algorithm. We compare the resulting daytime SLSTR AF product data to that generated from near-simultaneous views provided by MODIS onboard Terra. When both sensors detect the same active fire cluster at similar time, there is minimal bias shown between the two FRP retrievals (the ordinary least squares linear best fit between matched SLSTR and MODIS per-fire FRP matchups has a slope of 0.97). At the regional scale, the S3 product detects 70% of the AF pixels that the matching MODIS product reports, but also provides a further (16%) set of unique AF pixel detections. Regional FRP totals derived from SLSTR appear slightly lower than those from MODIS, and the OLS linear best fit between these regional FRP matchup datasets has a slope of 0.91. This is largely due to SLSTR performing less well in detecting the lowest FRP fires by day, whereas by night the S3 product performs a little better than MODIS due to the increased night-time use of S7 in the earlier AF pixel detection stages. Global fire mapping at a 0.25° grid cell resolution shows very similar daytime fire patterns and FRP totals from S3 and Terra MODIS, with SLSTR detecting around twice the number of AF pixels due to the algorithm being more effective at identifying low FRP pixels at the edges of fire clusters. Regional time series case studies also show very similar temporal patterns between S3 and Terra MODIS. Longer-term intercomparisons such as these will provide the knowledge necessary to use MODIS and SLSTR AF products together to analyse long-term AF trends. Comparing near simultaneous observations of fires by SLSTR and from the 30 m spatial resolution Landsat Operational Land Image (OLI) data, we find that once there are around 150 OLI active fire pixels detected within the area of an SLSTR pixel, the chances of that SLSTR pixel being classed as an active fire by the daytime algorithm rises to almost 100%. The daytime SLSTR AF Detection and FRP product based on the algorithm described herein has been fully operational since March 2022 and is available from the Sentinel-3 Science Hub (https://scihub.copernicus.eu/).

Systematic Evaluation of Four Satellite AOD Datasets for Estimating PM2.5 Using a Random Forest Approach

The latest epidemiological studies have revealed that the adverse health effects of PM2.5 have impacts beyond respiratory and cardio-vascular diseases and also affect the development of the brain and metabolic diseases. The need for accurate and spatio-temporally resolved PM2.5 data has thus been substantiated. While the selective information provided by station measurements is mostly insufficient for area-wide monitoring, satellite data have been increasingly applied to comprehensively monitor PM2.5 distributions. Although the accuracy and reliability of satellite-based PM2.5 estimations have increased, most studies still rely on a single sensor. However, several datasets have become available in the meantime, which raises the need for a systematic analysis. This study presents the first systematic evaluation of four satellite-based AOD datasets obtained from different sensors and retrieval methodologies to derive ground-level PM2.5 concentrations. We apply a random forest approach and analyze the effect of the resolution and coverage of the satellite data and the impact of proxy data on the performance. We examine AOD data from the Moderate resolution Imaging spectroradiometer (MODIS) onboard Terra and Aqua satellites, including Dark Target (DT) algorithm products and the Multi-Angle Implementation of Atmospheric Correction (MAIAC) product. Additionally, we explore more recent datasets from the Sea and Land Surface Temperature Radiometer (SLSTR) onboard Sentinel-3a and from the Tropospheric Monitoring Instrument (TROPOMI) operating on the Sentinel-5 precursor (S5p). The method is demonstrated for Germany and the year 2018, where a dense in situ measurement network and relevant proxy data are available. Overall, the model performance is satisfactory for all four datasets with cross-validated R2 values ranging from 0.68 to 0.77 and excellent for MODIS AOD reaching correlations of almost 0.9. We find a strong dependency of the model performance on the coverage and resolution of the AOD training data. Feature importance rankings show that AOD has less weight compared to proxy data for SLSTR and TROPOMI.

Improving the combined use of reflectance and thermal channels for ocean and coastal cloud detection for the Sea and Land Surface Temperature Radiometer (SLSTR)

Reflectance imagery is used to aid daytime cloud detection in thermal remote sensing. This paper presents new approaches to utilization of reflectance for sea surface temperature (SST) remote sensing for the Sea and Land Surface Temperature Radiometer (SLSTR). SLSTR is an along-track scanning sensor with a complex instrumental viewing geometry, and no co-registration of the fields-of-view of reflectance and thermal pixels. Reflectance channels have twice the spatial resolution of thermal channels, and observations are placed on compatible “image grids” for the convenience of users of Level-1 data. We highlight limitations of simple methods, based on these image grids, of using reflectance imagery to inform cloud detection at the thermal resolution. We present improvements from averaging the N-nearest reflectance observations directly to the infrared instrument geometry, where N = 10 is chosen in this study when using the A and B stripes together. We show that the standard deviation of the N-nearest reflectance observations is another calculable quantity of use to improve the discrimination of clouds in the infrared image, beneficially reducing the weight placed on coarser-scale thermal spatial variability. The developments are illustrated by case studies in coastal zones over optically bright waters and around strong ocean fronts, and the benefit for SST products is quantified by the impacts on coverage and validation statistics. In a case study over optically bright waters, the clear-sky fraction increases from 46.2% to 93.1%. Coastal zone validation shows a 23.8-31.6% reduction in the false alarm rate and a corresponding 27.1-33.3% increase in the cloud detection true skill score. Globally the robust standard deviation for clear-sky matches between SLSTR-A and drifting buoys reduces from 0.3 to 0.29 with a 6% reduction in data due to improved screening of scattered cloud.

The FastVRP automatic platform for the thermal monitoring of volcanic activity using VIIRS and SLSTR sensors: FastFRP to monitor volcanic radiative power

Satellite thermal remote sensing is widely used to detect and quantify the high-temperature vol- canic features produced during an eruption, e.g. released radiative power. Some space agencies provide Fire Radiative Power (FRP) Products to characterize any thermal anomaly around the world. In particular, Level-2 FRP Products of the Visible Infrared Imaging Radiometer Suite (VIIRS) and the Sea and Land Surface Temperature Radiometer (SLSTR) are freely available online and they allow to monitor high-temperature volcanic features related to the dynamics of volcanic activity. Here, we propose the FastVRP platform developed in Google Colab to process automatically the FRP Products provided by the National Aeronautics and Space Administration (NASA) and the European Space Agency (ESA) space agencies. FastVRP was designed to monitor the volcanic radiative power (VRP) related to eruptive activity of Mt. Etna (Sicily, Italy). We compared the quality of these FRP Products during a number of recent paroxysmal lava fountains occurred at Etna volcano between February and March 2021. We highlighted the advantages and the limits of each sensor in monitor- ing intense volcanic eruptions lasting a few hours. Furthermore, we combine the mid-high spatial/ low temporal resolution VIIRS and SLSTR with the low spatial-high temporal resolution SEVIRI (Spinning Enhanced Visible and Infrared Radiometer Imager) to improve estimates of the energies released from each paroxysmal episode. In particular, we propose a fitting approach to enhance the accuracy of SEVIRI low spatial-high temporal resolution measurements exploiting the few acqui- sitions from VIIRS and SLSTR high spatial-low temporal resolution during lava fountain cooling phase. We validated the radiative power values forecasted from VIIRS and SLSTR with the radiative power values retrieved using MODIS (Moderate Resolution Imaging Spectroradiometer) sensor.

Global estimates of surface albedo from Sentinel-3 OLCI and SLSTR data for Copernicus Climate Change Service: Algorithm and preliminary validation

The aim of Copernicus Climate Change Service (C3S) is to supply reliable climate data in support of strategies to adaptation and mitigation to climate change. The C3S provides access to high-quality climate data through its Climate Data Records (CDRs) of atmospheric, marine and land Essential Climate Variables (ECVs). Global Earth Surface Albedo (SA) satellite-based products are included in the land (biosphere) portfolio. SA is a magnitude which quantifies the fraction of solar energy reflected by the surface of the Earth. This paper details the retrieval methodology and preliminary validation results for global estimates of surface albedo based on Sentinel-3 observations for the C3S ECVs data (C3S SA v3.0). The retrieval algorithm exploits the synergistic use of the Ocean and Land Colour Instrument (OLCI) and the Sea and Land Surface Temperature Radiometer (SLSTR) on-board Sentinel-3 A and B satellites. Firstly, the atmospherically corrected reflectances are generated in the Copernicus Global Land Service framework. After that, the Bidirectional Reflectance Distribution Function (BRDF) inversion module concludes the BRDF model parameters, which are transferred to the angular integration module in order to generate spectral albedo quantities for the selected OLCI (Oa03, Oa04, Oa07, Oa17 and Oa21) and SLSTR (S1, S2, S5 and S6) bands. At the end, the spectral integration module generates broadband albedo quantities in three different standard broadband spectral regions (visible [0.4μm − 0.7μm ], near infrared [0.7μm − 4μm] and total shortwave [0.3μm − 4μm]). Preliminary validation results over 10-months demonstration period (July 2018–April 2019) show, in terms of spatial and temporal consistency, that C3S Sentinel-3 SA global estimates reached in general good agreement as compared to other satellite operational references derived from MODIS (MCD43A3 C6) and PROBA-V (C3S PROBA-V SA v1.0) acquisitions. The comparison with ground data shows similar results to the MCD43A3 C6 comparisons but opposite sign in differences (marginally positive in case of Sentinel-3), with accuracy of 0.005 (3.7%), precision of 0.016 (11.3%) and uncertainty of 0.032 (22.7%). Our results have demonstrated the feasibility to estimate global fields of SA from Sentinel-3 observations, with similar quality of existing operational products. These Sentinel-3 based SA datasets will give the continuity to the existing C3S SA CDR, introducing improvements in terms of spatial resolution (300 m) and spectral information (9 spectral albedos) in contrast to previous datasets based on Advanced Very High Resolution Radiometer (AVHRR; 4 km, 4 channels) and Vegetation instruments (VGT; 1 km, 4 channels).

Land Surface Temperature Retrieval From Sentinel-3A SLSTR Data: Comparison Among Split-Window, Dual-Window, Three-Channel, and Dual-Angle Algorithms

Land surface temperature (LST) is a vital parameter for studying global ecological, climatic, and environmental changes. Although various LST retrieval algorithms have been proposed, including split-window (SW), dual-window (DW), three-channel (TC), and dual-angle (DA) algorithms, few studies have compared these algorithms using the same satellite observations. The Sea and Land Surface Temperature Radiometer (SLSTR) onboard Sentinel-3A provides a unique opportunity to conduct this comparison owing to its dual-angle viewing capability and multiple thermal infrared (TIR) and mid-infrared (MIR) channels. Here, we implemented two SW algorithms, one DW algorithm, two TC algorithms and one DA algorithm for the SLSTR data. The LST retrievals from these six algorithms were validated, along with the SLSTR operational LST product based on an emissivity-implicit SW algorithm. Temperature-based and radiance-based validation methods were used to evaluate different LST retrievals across different land cover types. The results indicated that the proposed SW algorithm had the highest accuracy, followed by the Pérez-Planells SW and the official algorithms. The overall root-mean-square errors (RMSEs) of these three SW algorithms were 1.42 K, 1.79 K and 2.05 K, respectively. The three algorithms involving the MIR channel (one DW and two TC algorithms) were more suitable for nighttime LST retrieval and had similar performances to the three SW algorithms, with a nighttime RMSE of approximately 1.36 K. The LST retrieval accuracy of the DA algorithm had the highest uncertainty and was closely related to the angular variation in surface emissivity and brightness temperature. The findings of this study contribute to a better understanding of the different LST retrieval algorithms and facilitate potential improvements in the official LST retrieval algorithm for SLSTR.

Evaluation of the Ability of SLSTR (Sentinel-3B) and MODIS (Terra) Images to Detect Burned Areas Using Spatial-Temporal Attributes and SVM Classification

Forest fires are considered one of the major dangers and environmental issues across the world. In the Cerrado biome (Brazilian savannas), forest fires have several consequences, including increased temperature, decreased rainfall, genetic depletion of natural species, and increased risk of respiratory diseases. This study presents a methodology that uses data from the Sea and Land Surface Temperature Radiometer (SLSTR) sensor of the Sentinel-3B satellite and the Moderate Resolution Imaging Spectroradiometer (MODIS) of the Terra satellite to analyze the thematic accuracy of burned area maps and their sensitivity under different spectral resolutions in a large area of 32,000 km2 in the Cerrado biome from 2019 to 2021. The methodology used training and the Support Vector Machine (SVM) classifier. To analyze the spectral peculiarities of each orbital platform, the Transformed Divergence (TD) index separability statistic was used. The results showed that for both sensors, the near-infrared (NIR) band has an essential role in the detection of the burned areas, presenting high separability. Overall, it was possible to observe that the spectral mixing problems, registration date, and the spatial resolution of 500 m were the main factors that led to commission errors ranging between 15% and 72% and omission errors between 51% and 86% for both sensors. This study showed the importance of multispectral sensors for monitoring forest fires. It was found, however, that the spectral resolution and burning date may gradually interfere with the detection process.

Volcanic cloud detection using Sentinel-3 satellite data by means of neural networks: the Raikoke 2019 eruption test case

Abstract. Accurate automatic volcanic cloud detection by means of satellite data is a challenging task and is of great concern for both the scientific community and aviation stakeholders due to well-known issues generated by strong eruption events in relation to aviation safety and health impacts. In this context, machine learning techniques applied to satellite data acquired from recent spaceborne sensors have shown promising results in the last few years. This work focuses on the application of a neural-network-based model to Sentinel-3 SLSTR (Sea and Land Surface Temperature Radiometer) daytime products in order to detect volcanic ash plumes generated by the 2019 Raikoke eruption. A classification of meteorological clouds and of other surfaces comprising the scene is also carried out. The neural network has been trained with MODIS (Moderate Resolution Imaging Spectroradiometer) daytime imagery collected during the 2010 Eyjafjallajökull eruption. The similar acquisition channels of SLSTR and MODIS sensors and the comparable latitudes of the eruptions permit an extension of the approach to SLSTR, thereby overcoming the lack in Sentinel-3 products collected in previous mid- to high-latitude eruptions. The results show that the neural network model is able to detect volcanic ash with good accuracy if compared to RGB visual inspection and BTD (brightness temperature difference) procedures. Moreover, the comparison between the ash cloud obtained by the neural network (NN) and a plume mask manually generated for the specific SLSTR images considered shows significant agreement, with an F-measure of around 0.7. Thus, the proposed approach allows for an automatic image classification during eruption events, and it is also considerably faster than time-consuming manual algorithms. Furthermore, the whole image classification indicates the overall reliability of the algorithm, particularly for recognition and discrimination between volcanic clouds and other objects.

Monitoring multiple satellite aerosol optical depth (AOD) products within the Copernicus Atmosphere Monitoring Service (CAMS) data assimilation system

Abstract. The Copernicus Atmosphere Monitoring Service (CAMS) provides near-real-time forecast and reanalysis of aerosols using the ECMWF Integrated Forecasting System with atmospheric composition extension, constrained by the assimilation of MODIS and the Polar Multi-Sensor Aerosol Optical Properties (PMAp) aerosol optical depth (AOD). The objective of this work is to evaluate two new near-real-time AOD products to prepare for their assimilation into CAMS, namely the Copernicus AOD (collection 1) from the Sea and Land Surface Temperature Radiometer (SLSTR) on board Sentinel 3-A/B over ocean and the NOAA EPS AOD (v2.r1) from VIIRS on board S-NPP and NOAA20 over both land and ocean. The differences between MODIS (C6.1), PMAp (v2.1), VIIRS (v2.r1), and SLSTR (C1) AOD as well as their departure from the modeled AOD were assessed at the model grid resolution (i.e., level-3) using the 3-month AOD average (December 2019–February 2020 and March–May 2020). VIIRS and MODIS show the best consistency across the products, which is explained by instrument and retrieval algorithm similarities. VIIRS AOD is frequently lower over the ocean background and higher over biomass burning and dust source land regions compared to MODIS. VIIRS shows larger spatial coverage over land and resolves finer spatial structures such as the transport of Australian biomass burning smoke over the Pacific, which can be explained by the use of a heavy aerosol detection test in the retrieval algorithm. Our results confirm the positive offset over ocean (i) between Terra/MODIS and Aqua/MODIS due to the non-corrected radiometric calibration degradation of Terra/MODIS in the Dark Target algorithm and (ii) between SNPP/VIIRS and NOAA20/VIIRS due to the positive bias in the solar reflective bands of SNPP/VIIRS. SLSTR AOD shows much smaller level-3 values than the rest of the products, which is mainly related to differences in spatial representativity at the IFS grid spatial resolution due to the stringent cloud filtering applied to the SLSTR radiances. Finally, the geometry characteristics of the instrument, which drive the range of scattering angles sampled by the instrument, can explain a large part of the differences between retrievals such as the positive offset between PMAp datasets from MetOp-B and MetOp-A.

Extended validation and evaluation of the OLCI–SLSTR SYNERGY aerosol product (SY_2_AOD) on Sentinel-3

Abstract. We present the first extended validation of a new SYNERGY global aerosol product (SY_2_AOD), which is based on synergistic use of data from the Ocean and Land Color Instrument (OLCI) and the Sea and Land Surface Temperature Radiometer (SLSTR) sensors aboard the Copernicus Sentinel-3A (S3A) and Sentinel-3B (S3B) satellites. Validation covers period from 14 January 2020 to 30 September 2021. Several approaches, including statistical analysis, time series analysis, and comparison with similar aerosol products from the other spaceborne sensor, the Moderate Resolution Imaging Spectroradiometer (MODIS), were applied for validation and evaluation of S3A and S3B SY_2 aerosol products, including aerosol optical depth (AOD) provided at different wavelengths, AOD pixel-level uncertainties, fine-mode AOD, and Angström exponent. Over ocean, the performance of SY_2 AOD (syAOD) retrieved at 550 nm is good: for S3A and S3B, Pearson correlation coefficients with the Maritime Aerosol Network (MAN) component of the AErosol RObotic NETwork (AERONET) are 0.88 and 0.85, respectively; 88.6 % and 89.5 % of pixels fit into the MODIS error envelope (EE) of ±0.05 ± 0.2 × AOD. Over land, correlation coefficients with AERONET AOD (aAOD) are 0.60 and 0.63 for S3A and S3B, respectively; 51.4 % and 57.9 % of pixels fit into MODIS EE. Reduced performance over land is expected since the surface reflectance and angular distribution of scattering are higher and more difficult to predict over land than over ocean. The results are affected by a large number of outliers. Evaluation of the per-retrieval uncertainty with the χ2 test indicates that syAOD prognostic uncertainties (PU) are slightly underestimated (χ2 = 3.1); if outliers are removed, PU describes the syAOD error well (χ2 = 1.6). The regional analysis of the Angström exponent, which relates to the aerosol size distribution, shows spatial correlation with expected sources. For 40 % of the matchups with AERONET in the Northern Hemisphere (NH) and for 60 % of the matchups in the Southern Hemisphere (SH), which fit into the AE size range of [1, 1.8], an offset between SY_2 AE (syAE) and AERONET AE (aAE) is within ±0.25. General overestimation of low (< 0.5) syAE and underestimation of high (> 1.8) syAE results in high (0.94, globally) overall bias. Good agreement (bias < 0.03) was observed between Sy_2 fine-mode AOD (syFMAOD) and AERONET fine-mode AOD (aFMAOD) for aFMAOD < 1. At aFMAOD > 1, syFMAOD is considerably underestimated (by 0.3–0.5 in different aFMAOD ranges) in the NH. In the SH, only a few aFMAOD values above 1 are measured. The fine-mode fraction (FMF) in the SY_2 AOD product (syFMF) in the range of [0, 0.7] is overestimated; the positive offset of 0.3–0.5 for low (< 0.25) FMF gradually decreases. Differences between the annual and seasonal AOD values from SY_2 and MODIS (mod) Dark Target and Deep Blue products are within 0.02 for the study area (30∘ S–60∘ N, 80∘ W–45∘ E). The agreement is better over ocean; however, a difference up to 0.6 exists between syFMF and modFMF. Over bright land surface (Saharan desert) the difference in AOD between the two products is highest (up to 0.11); the sign of the difference varies over time and space. For both S3A and S3B AOD products, validation statistics are often slightly better in the Southern Hemisphere. In general, the performance of S3B is slightly better.

Synergistic Use of the SRAL/MWR and SLSTR Sensors on Board Sentinel-3 for the Wet Tropospheric Correction Retrieval

The Sentinel-3 satellites are equipped with dual-band Microwave Radiometers (MWR) to derive the wet tropospheric correction (WTC) for satellite altimetry. The deployed MWR lack the 18 GHz channel, which mainly provides information on the surface emissivity. Currently, this information is considered using additional parameters, one of which is the sea surface temperature (SST) extracted from static seasonal tables. Recent studies show that the use of a dynamic SST extracted from Numerical Weather Models (ERA5) improves the WTC retrieval. Given that Sentinel-3 carries on board the Sea and Land Surface Temperature Radiometer (SLSTR), from which SST observations are derived simultaneously with those of the Synthetic Aperture Radar Altimeter and MWR sensors, this study aims to develop a synergistic approach between these sensors for the WTC retrieval over open ocean. Firstly, the SLSTR-derived SSTs are evaluated against the ERA5 model; secondly, their impact on the WTC retrieval is assessed. The results show that using the SST input from SLSTR, instead of ERA5, has no impact on the WTC retrieval, both globally and regionally. Thus, for the WTC retrieval, there seems to be no advantage in having collocated SST and radiometer observations. Additionally, this study reinforces the fact that the use of dynamic SST leads to a significant improvement over the current Sentinel-3 WTC operational algorithms.

Bayesian Cloud Detection over Land for Climate Data Records

Cloud detection is a necessary step in the generation of land surface temperature (LST) climate data records (CDRs) and affects data quality and uncertainty. We present here a sensor-independent Bayesian cloud detection algorithm and show that it is suitable for use in the production of LST CDRs. We evaluate the performance of the cloud detection with reference to two manually masked datasets for the Advanced Along-Track Scanning Radiometer (AATSR) and find a 7.9% increase in the hit rate and 4.9% decrease in the false alarm rate when compared to the operational cloud mask. We then apply the algorithm to four instruments aboard polar-orbiting satellites, which together can produce a global, 25-year LST CDR: the second Along-Track Scanning Radiometer (ATSR-2), AATSR, the Moderate Resolution Spectroradiometer (MODIS Terra) and the Sea and Land Surface Temperature Radiometer (SLSTR-A). The Bayesian cloud detection hit rate is assessed with respect to in situ ceilometer measurements for periods of overlap between sensors. The consistency of the hit rate is assessed between sensors, with mean differences in the cloud hit rate of 4.5% for ATSR-2 vs. AATSR, 4.9% for AATSR vs. MODIS, and 2.5% for MODIS vs. SLSTR-A. This is important because consistent cloud detection performance is needed for the observational stability of a CDR. The application of a sensor-independent cloud detection scheme in the production of CDRs is thus shown to be a viable approach to achieving LST observational stability over time.

Estimating PM2.5 surface concentrations from AOD: A combination of SLSTR and MODIS

Compared to surface in-situ observations, satellite data on aerosol optical depth (AOD) enables area-wide monitoring of tropospheric aerosols. However, coverage and reliability of satellite data products depend on atmospheric conditions and surface concentrations have to be retrieved from AOD. This study investigates the potential to produce reliable maps of PM2.5 surface concentrations for Germany and parts of the surrounding countries using AOD based on observations by three different satellite sensors. For the first time, AOD retrievals from the Sea and Land Surface Temperature Radiometer (SLSTR) onboard Sentinel-3A are used together with those from the Moderate Resolution Imaging Spectroradiometer (MODIS) onboard the two NASA platforms Terra and Aqua. We investigate the differences and similarities of the three different satellite products in terms of coverage, resolution and algorithmic performances. Based on this analysis we examine the suitability and advantage of a combination of these data sets. We can substantiate an increase in mean daily coverage from a maximum of 10.2% for the individual products to 16.7% for the ensemble product. Using a semi-empirical linear regression model, we derive surface-level PM2.5 concentrations and attain an overall correlation of 0.76 between satellite-derived and in-situ measured PM2.5 concentrations. By considering surface measurements, the systematic error (bias) and the root mean square error (RMSE) can be significantly reduced. The general model performance is evaluated by a 5-fold cross validation and the relative prediction error (RPE).

Towards a novel approach for Sentinel-3 synergistic OLCI/SLSTR cloud and cloud shadow detection based on stereo cloud-top height estimation

Sentinel-3 is an Earth observation satellite constellation launched by the European Space Agency. Each satellite carries two optical multispectral instruments: the Ocean and Land Colour Instrument (OLCI) and the Sea and Land Surface Temperature Radiometer (SLSTR). OLCI and SLSTR sensors produce images covering the visible and infrared spectrum that can be collocated in order to generate synergistic products. In Earth observation, a particular weakness of optical sensors is their high sensitivity to clouds and their shadows. An incorrect cloud and cloud shadow detection leads to mistakes in both land and ocean retrievals of biophysical parameters. In order to exploit both OLCI and SLSTR capabilities, image co-registration at ground level is needed. However, applying such collocation of the images results in cloud location mismatches due to the different viewing angles of OLCI and SLSTR, which complicates the synergistic cloud detection. This study seeks to provide a solution to correctly obtain the projected clouds based on the estimation of cloud top heights in order to better collocate clouds between sensors and detect their shadows. The study presents a forward and backward method to estimate the real nadir position of a cloud on the satellite image starting from an existing cloud mask, as well as the corresponding cloud projections on the surface depending on the solar and sensor viewing angles. The estimation of cloud top heights is based on differences in the cloud projections from SLSTR nadir and oblique views. Experimental results show that the stereo cloud matching based on maximum cross-correlation between SLSTR nadir and oblique spectra was the most robust method to match SLSTR clouds for both nadir and oblique views as compared to spectral distance and spectral angle minimization. We test the method over several images around the world, leading to higher overall accuracy (OA) as compared to Sentinel-3 official products, both in detecting SLSTR clouds and OLCI cloud shadows (SLSTR nadir OA = 93.6%, SLSTR oblique OA = 88.7%, OLCI cloud shadow OA = 93.9% for the stereo matcher, against 82.2%, 81.3% and 90.5%, respectively, for the official Sentinel-3 products). This study also provides a starting point in the development of a cloud screening approach for the upcoming Fluorescence Explorer (FLEX) satellite mission, expected to fly in tandem with Sentinel-3.

Drought Index Mapping in Java Island Using Sentinel-3 SLSTR

As the most populous island globally, Java Island has various vulnerabilities to disasters ranging from geological to hydro-meteorological. One of the most common hydro-meteorological disasters is the drought that occurs every year in the dry season. This disaster causes crop failure, land and forest fires, and clean water shortages. In this study, the Sea and Land Surface Temperature Radiometer (SLSTR) instrument onboard the Sentinel-3 platform was used to map drought using the Vegetation Temperature Condition Index (VTCI) algorithm based on the scattering space technique of Land Surface Temperature (LST) and Normalized Different Vegetation Index (NDVI). During 2018, the highest and lowest LST occurred on September 28 (309 oK), and February 1 (278 oK); the highest and lowest NDVI occurred on July 1 (0.466) and November 13 (0.221). In comparison, the driest conditions indicated the lowest VTCI (0.162) on July 17 and the wettest on August 28 (0.508). The driest and wettest situation co-occurred: at the end of the dry session and the rainy session, it shows that the rain greatly contributes to high VTCI. A regular drought mapping needs to be performed as an effort to disaster risk reduction. Drought maps are then used as a spatial recommendation in reforestation intervention to reduce drought in the future.

The retrieval of snow properties from SLSTR Sentinel-3 – Part 2: Results and validation

Abstract. To evaluate the performance of the eXtensible Bremen Aerosol/cloud and surfacE parameters Retrieval (XBAER) algorithm, presented in the Part 1 companion paper to this paper, we apply the XBAER algorithm to the Sea and Land Surface Temperature Radiometer (SLSTR) instrument on board Sentinel-3. Snow properties – snow grain size (SGS), snow particle shape (SPS) and specific surface area (SSA) – are derived under cloud-free conditions. XBAER-derived snow properties are compared to other existing satellite products and validated by ground-based and aircraft measurements. The atmospheric correction is performed on SLSTR for cloud-free scenarios using Modern-Era Retrospective Analysis for Research and Applications (MERRA) aerosol optical thickness (AOT) and the aerosol typing strategy according to the standard XBAER algorithm. The optimal SGS and SPS are estimated iteratively utilizing a look-up-table (LUT) approach, minimizing the difference between SLSTR-observed and SCIATRAN-simulated surface directional reflectances at 0.55 and 1.6 µm. The SSA is derived for a retrieved SGS and SPS pair. XBAER-derived SGS, SPS and SSA have been validated using in situ measurements from the recent campaign SnowEx17 during February 2017. The comparison shows a relative difference between the XBAER-derived SGS and SnowEx17-measured SGS of less than 4 %. The difference between the XBAER-derived SSA and SnowEx17-measured SSA is 2.7 m2/kg. XBAER-derived SPS can be reasonably explained by the SnowEx17-observed snow particle shapes. Intensive validation shows that (1) for SGS and SSA, XBAER-derived results show high correlation with field-based measurements, with correlation coefficients higher than 0.85. The root mean square errors (RMSEs) of SGS and SSA are around 12 µm and 6 m2/kg. (2) For SPS, aggregate SPS retrieved by XBAER algorithm is likely to be matched with rounded grains while single SPS in XBAER is possibly linked to faceted crystals. The comparison with aircraft measurements, during the Polar Airborne Measurements and Arctic Regional Climate Model Simulation Project (PAMARCMiP) campaign held in March 2018, also shows good agreement (with R=0.82 and R=0.81 for SGS and SSA, respectively). XBAER-derived SGS and SSA reveal the variability in the aircraft track of the PAMARCMiP campaign. The comparison between XBAER-derived SGS results and the Moderate Resolution Imaging Spectroradiometer (MODIS) Snow-Covered Area and Grain size (MODSCAG) product over Greenland shows similar spatial distributions. The geographic distribution of XBAER-derived SPS over Greenland and the whole Arctic can be reasonably explained by campaign-based and laboratory investigations, indicating a reasonable retrieval accuracy of the retrieved SPS. The geographic variabilities in XBAER-derived SGS and SSA both over Greenland and Arctic-wide agree with the snow metamorphism process.