Atmospheric Correction Inter-comparison eXercise, ACIX-II Land: An assessment of atmospheric correction processors for Landsat 8 and Sentinel-2 over land

The correction of the atmospheric effects on optical satellite images is essential for quantitative and multitemporal remote sensing applications. In order to study the performance of the state-of-the-art methods in an integrated way, a voluntary and open-access benchmark Atmospheric Correction Inter-comparison eXercise (ACIX) was initiated in 2016 in the frame of Committee on Earth Observation Satellites (CEOS) Working Group on Calibration & Validation (WGCV). The first exercise was extended in a second edition wherein twelve atmospheric correction (AC) processors, a substantially larger testing dataset and additional validation metrics were involved. The sites for the inter-comparison analysis were defined by investigating the full catalogue of the Aerosol Robotic Network (AERONET) sites for coincident measurements with satellites' overpass. Although there were more than one hundred sites for Copernicus Sentinel-2 and Landsat 8 acquisitions, the analysis presented in this paper concerns only the common matchups amongst all processors, reducing the number to 79 and 62 sites respectively. Aerosol Optical Depth (AOD) and Water Vapour (WV) retrievals were consequently validated based on the available AERONET observations. The processors mostly succeeded in retrieving AOD for relatively light to medium aerosol loading (AOD 90% of the results falling within the suggested empirical specifications and with the Root Mean Square Error (RMSE) being mostly <0.25 g/cm2. Regarding Surface Reflectance (SR) validation two main approaches were followed. For the first one, a simulated SR reference dataset was computed over all of the test sites by using the 6SV (Second Simulation of the Satellite Signal in the Solar Spectrum vector code) full radiative transfer modelling (RTM) and AERONET measurements for the required aerosol variables and water vapour content. The performance assessment demonstrated that the retrievals were not biased for most of the bands. The uncertainties ranged from approximately 0.003 to 0.01 (excluding B01) for the best performing processors in both sensors' analyses. For the second one, measurements from the radiometric calibration network RadCalNet over La Crau (France) and Gobabeb (Namibia) were involved in the validation. The performance of the processors was in general consistent across all bands for both sensors and with low standard deviations (<0.04) between on-site and estimated surface reflectance. Overall, our study provides a good insight of AC algorithms' performance to developers and users, pointing out similarities and differences for AOD, WV and SR retrievals. Such validation though still lacks of ground-based measurements of known uncertainty to better assess and characterize the uncertainties in SR retrievals

ACIX-Aqua: A global assessment of atmospheric correction methods for Landsat-8 and Sentinel-2 over lakes, rivers, and coastal waters

Atmospheric correction over inland and coastal waters is one of the major remaining challenges in aquatic remote sensing, often hindering the quantitative retrieval of biogeochemical variables and analysis of their spatial and temporal variability within aquatic environments. The Atmospheric Correction Intercomparison Exercise (ACIX-Aqua), a joint NASA – ESA activity, was initiated to enable a thorough evaluation of eight state-of-the-art atmospheric correction (AC) processors available for Landsat-8 and Sentinel-2 data processing. Over 1000 radiometric matchups from both freshwaters (rivers, lakes, reservoirs) and coastal waters were utilized to examine the quality of derived aquatic reflectances (̂ρw). This dataset originated from two sources: Data gathered from the international scientific community (henceforth called Community Validation Database, CVD), which captured predominantly inland water observations, and the Ocean Color component of AERONET measurements (AERONET-OC), representing primarily coastal ocean environments. This volume of data permitted the evaluation of the AC processors individually (using all the matchups) and comparatively (across seven different Optical Water Types, OWTs) using common matchups. We found that the performance of the AC processors differed for CVD and AERONET-OC matchups, likely reflecting inherent variability in aquatic and atmospheric properties between the two datasets. For the former, the median errors in ̂ρw(560) and ̂ρw(664) were found to range from 20 to 30% for best-performing processors. Using the AERONET-OC matchups, our performance assessments showed that median errors within the 15–30% range in these spectral bands may be achieved. The largest uncertainties were associated with the blue bands (25 to 60%) for best-performing processors considering both CVD and AERONET-OC assessments. We further assessed uncertainty propagation to the downstream products such as near-surface concentration of chlorophyll-a (Chla) and Total Suspended Solids (TSS). Using satellite matchups from the CVD along with in situ Chla and TSS, we found that 20–30% uncertainties in ̂ρw(490 ≤ λ ≤ 743 nm) yielded 25–70% uncertainties in derived Chla and TSS products for topperforming AC processors. We summarize our results using performance matrices guiding the satellite user community through the OWT-specific relative performance of AC processors. Our analysis stresses the need for better representation of aerosols, particularly absorbing ones, and improvements in corrections for sky- (or sun-) glint and adjacency effects, in order to achieve higher quality downstream products in freshwater and coastal ecosystems

Revisiting the Past: Replicability of a Historic Long-Term Vegetation Dynamics Assessment in the Era of Big Data Analytics

Open and analysis-ready data, as well as methodological and technical advancements have resulted in an unprecedented capability for observing the Earth’s land surfaces. Over 10 years ago, Landsat time series analyses were inevitably limited to a few expensive images from carefully selected acquisition dates. Yet, such a static selection may have introduced uncertainties when spatial or inter-annual variability in seasonal vegetation growth were large. As seminal pre-open-data-era papers are still heavily cited, variations of their workflows are still widely used, too. Thus, here we quantitatively assessed the level of agreement between an approach using carefully selected images and a state-of-the-art analysis that uses all available images. We reproduced a representative case study from the year 2003 that for the first time used annual Landsat time series to assess long-term vegetation dynamics in a semi-arid Mediterranean ecosystem in Crete, Greece. We replicated this assessment using all available data paired with a time series method based on land surface phenology metrics. Results differed fundamentally because the volatile timing of statically selected images relative to the phenological cycle introduced systematic uncertainty. We further applied lessons learned to arrive at a more nuanced and information-enriched vegetation dynamics description by decomposing vegetation cover into woody and herbaceous components, followed by a syndrome-based classification of change and trend parameters. This allowed for a more reliable interpretation of vegetation changes and even permitted us to disentangle certain land-use change processes with opposite trajectories in the vegetation components that were not observable when solely analyzing total vegetation cover. The long-term budget of net cover change revealed that vegetation cover of both components has increased at large and that this process was mainly driven by gradual processes. We conclude that study designs based on static image selection strategies should be critically evaluated in the light of current data availability, analytical capabilities, and with regards to the ecosystem under investigation. We recommend using all available data and taking advantage of phenology-based approaches that remove the selection bias and hence reduce uncertainties in results.Peer Reviewe

Evaluation of SEN2COR surface reflectance products over land surface with reference measurements on ground

Sen2Cor is the atmospheric correction processor selected by ESA for operational, systematic processing of Copernicus Sentinel-2 mission data. It is used for generating the Level 2A products distributed to users by the Copernicus SciHub. Accurate atmospheric correction of Sentinel-2 data and knowledge of its uncertainties are preconditions for high quality downstream applications. In this work we present the comparison of Sentinel-2 Bottom-of-Atmosphere products with measurements of surface reflectance on ground. Source of reference measurements are both surface reflectance data from RadCalNet and from dedicated field campaigns. The analysis shows, that the uncertainty of SR-retrieval with Sen2Cor is better than about 7% for bright surfaces and about 17% for darker. In addition to this performance evaluation, the data are also applied to compare the use of reference data coming from permanent operating bright RadCalNet sites and from ad-hoc field campaigns at darker sites

Comparing Atmospheric Correction Performance for Sentinel-2 and Landsat-8 data

In terms of atmospheric impact, the volcanic eruption of Mt. Pinatubo (1991) is thebest characterized large eruption on record. We investigate here the stratosphericwarming following the Pinatubo eruption derived from SAGE II extinction data includ-ing most recent improvements in the processing algorithm and a data filling procedure in the opacity-induced “gap” regions. From these data, which cover wavelengths of1.024 micrometer and shorter, we derived aerosol size distributions which properly re-produce extinction coefficients at much longer wavelength. This provides a good basisfor calculating the absorption of terrestrial infrared radiation and the resulting strato-spheric heating. However, we also show that the use of this dataset in the global chemistry-climate model (CCM) SOCOL leads to exaggerated aerosol-induced strato-spheric heating compared to observations, even partly larger than the already too highvalues found by many models in recent general circulation model (GCM) and CCMintercomparisons. This suggests that the overestimation of the stratospheric warm-ing after the Pinatubo eruption arises from deficiencies in the model radiation codes rather than an insufficient observational data basis. Conversely, our approach reducesthe infrared absorption in the tropical tropopause region, in better agreement with thepost-volcanic temperature record at these altitudes.ISSN:1680-7375ISSN:1680-736

ACIX-Aqua: A global assessment of atmospheric correction methods for Landsat-8 and Sentinel-2 over lakes, rivers, and coastal waters

Atmospheric correction over inland and coastal waters is one of the major remaining challenges in aquatic remote sensing, often hindering the quantitative retrieval of biogeochemical variables and analysis of their spatial and temporal variability within aquatic environments. The Atmospheric Correction Intercomparison Exercise (ACIX-Aqua), a joint NASA – ESA activity, was initiated to enable a thorough evaluation of eight state-of-the-art atmospheric correction (AC) processors available for Landsat-8 and Sentinel-2 data processing. Over 1000 radiometric matchups from both freshwaters (rivers, lakes, reservoirs) and coastal waters were utilized to examine the quality of derived aquatic reflectances (̂ρw). This dataset originated from two sources: Data gathered from the international scientific community (henceforth called Community Validation Database, CVD), which captured predominantly inland water observations, and the Ocean Color component of AERONET measurements (AERONET-OC), representing primarily coastal ocean environments. This volume of data permitted the evaluation of the AC processors individually (using all the matchups) and comparatively (across seven different Optical Water Types, OWTs) using common matchups. We found that the performance of the AC processors differed for CVD and AERONET-OC matchups, likely reflecting inherent variability in aquatic and atmospheric properties between the two datasets. For the former, the median errors in ̂ρw(560) and ̂ρw(664) were found to range from 20 to 30% for best-performing processors. Using the AERONET-OC matchups, our performance assessments showed that median errors within the 15–30% range in these spectral bands may be achieved. The largest uncertainties were associated with the blue bands (25 to 60%) for best-performing processors considering both CVD and AERONET-OC assessments. We further assessed uncertainty propagation to the downstream products such as near-surface concentration of chlorophyll-a (Chla) and Total Suspended Solids (TSS). Using satellite matchups from the CVD along with in situ Chla and TSS, we found that 20–30% uncertainties in ̂ρw(490 ≤ λ ≤ 743 nm) yielded 25–70% uncertainties in derived Chla and TSS products for topperforming AC processors. We summarize our results using performance matrices guiding the satellite user community through the OWT-specific relative performance of AC processors. Our analysis stresses the need for better representation of aerosols, particularly absorbing ones, and improvements in corrections for sky- (or sun-) glint and adjacency effects, in order to achieve higher quality downstream products in freshwater and coastal ecosystems

Sentinel-2 Level-2 processing Sen2Cor status and outlook of 2022

The Sentinel-2 mission is fully operating since June 2017 with a constellation of two polar orbiting satellite units. Both Sentinel-2A and Sentinel-2B are equipped with an optical imaging sensor MSI (MultiSpectral Instrument) which acquires high spatial resolution optical data products. The Sentinel-2 mission is dedicated to land monitoring, emergency management and security. It serves for monitoring of land-cover change and biophysical variables related to agriculture and forestry, monitors coastal and inland waters and is useful for risk and disaster mapping. Accurate atmospheric correction of satellite observations is a precondition for the development and delivery of high quality applications. Therefore the atmospheric correction processor Sen2Cor was developed with the objective of delivering land surface reflectance products. Sen2Cor is designed to process single tile Level-1C products, providing Level-2A surface (Bottom-of-Atmosphere) reflectance product together with Aerosol Optical Thickness (AOT), Water Vapour (WV) estimation maps and a Scene Classification (SCL) map including cloud / cloud shadow classes for further processing. Sen2Cor processor can be downloaded from ESA website as a standalone tool for individual Level-2A processing by the users. It can be run either from command line or as a plugin of the Sentinel-2 Toolbox (SNAPS2TBX). In parallel, ESA started in June 2017 to use Sen2Cor for systematic Level-2A processing of Sentinel-2 acquisitions over Europe. Since March 2018, Level-2A products are generated by the official Sentinel-2 ground segment (PDGS) and are available on the Copernicus Open Access Hub. The objective of this presentation is to provide users with an overview of the Level-2A product contents and up-to-date information about the data quality of the Level-2A products (processing baseline PB.04.00 onwards) generated by Sentinel-2 PDGS since end of January 2022, in terms of Cloud Screening and Atmospheric Correction. In addition, the presentation will give an outlook on the recent updates of Sen2Cor, which improve L2A Data Quality: updated L2A metadata, updated scene classification, updated fall-back method using meteorological information from the Copernicus Atmosphere Monitoring Service, updated Copernicus DEM

Atmospheric Correction of Sentinel 2A Data and Influence on Winter Wheat Phenology in Kabul Region, Afghanistan

No abstract availabl

Application of machine learning techniques to derive sea water turbidity from Sentinel-2 imagery

Earth Observation (EO) from satellites has the potential to provide comprehensive, rapid and inexpensive information about water bodies, integrating in situ measurements. Traditional methods to retrieve optically active water quality parameters from satellite data are based on semiempirical models relying on few bands, which often revealed to be site and season specific. The use of machine learning (ML) for remotely sensed water quality estimation has spread in recent years thanks to the advances in algorithm development and computing power. These models allow to exploit the wealth of spectral information through more flexible relationships and are less affected by atmospheric and other background factors. The present study explores the use of Sentinel-2 MultiSpectral Instrument (MSI) Level-1C Top of Atmosphere spectral radiance to derive water turbidity, through application of machine learning techniques. A dataset of 222 combination of turbidity measurements, collected in the North Tyrrhenian Sea – Italy from 2015 to 2021, and values of the 13 spectral bands in the pixel corresponding to the sample location was used. Two regression techniques were tested and compared: a Stepwise Linear Regression (SLR) and a Polynomial Kernel Regression. The two models show accurate and similar performance (R2 = 0.736, RMSE = 2.03 NTU, MAE = 1.39 NTU for the SLR and R2 = 0.725, RMSE = 2.07 NTU, MAE = 1.40 NTU for the Kernel). A band importance analysis revealed the contribution of the different spectral bands and the main role of the red-edge range. The work shows that it is possible to reach a good accuracy in turbidity estimation from MSI TOA reflectance using ML models, fed by the whole spectrum of available bands, although the possible generation of errors related to atmospheric effect in turbidity estimates was not evaluated. Comparison between turbidity estimates obtained from the models with turbidity data from Copernicus CMEMS dataset named ‘Mediterranean Sea, Bio-Geo-Chemical, L3, daily observation’ produced consistent results. Finally, turbidity maps from satellite imagery were produced for the study area, showing the ability of the models to catch extreme events