Medium Resolution Sensors Research Articles

Spatial structure of in situ reflectance in coastal and inland waters: implications for satellite validation

Validation of satellite-derived aquatic reflectance involves relating meter-scale in situ observations to satellite pixels with typical spatial resolution ∼ 10–100 m within a temporal “match-up window” of an overpass. Due to sub-pixel variation these discrepancies in measurement scale are a source of uncertainty in the validation result. Additionally, validation protocols and statistics do not normally account for spatial autocorrelation when pairing in situ data from moving platforms with satellite pixels. Here, using high-frequency autonomous mobile radiometers deployed on ships, we characterize the spatial structure of in situRrs in inland and coastal waters (Lake Balaton, Western English Channel, Tagus Estuary). Using variogram analysis, we partition Rrs variability into spatial and intrinsic (non-spatial) components. We then demonstrate the capacity of mobile radiometers to spatially sample in situRrs within a temporal window broadly representative of satellite validation and provide spatial statistics to aid satellite validation practice. At a length scale typical of a medium resolution sensor (300 m) between 5% and 35% (median values across spectral bands and deployments) of the variation in in situRrs was due to spatial separation. This result illustrates the extent to which mobile radiometers can reduce validation uncertainty due to spatial discrepancy via sub-pixel sampling. The length scale at which in situRrs became spatially decorrelated ranged from ∼ 100–1,000 m. This information serves as a guideline for selection of spatially independent in situRrs when matching with a satellite image, emphasizing the need for either downsampling or using modified statistics when selecting data to validate high resolution sensors (sub 100 m pixel size).

Satellite Remote Sensing for Wildlife Research in the Polar Regions

Abstract Wildlife research in the polar regions has historically been limited by the logistical constraints of site access, but recent developments in the use of satellite imagery for animal detection has unlocked new possibilities for pan-Arctic and pan-Antarctic monitoring of animal populations. A range of different sensor systems have been used for wildlife research, but most have focused on optical sensors that collect data in the visible spectrum and can be directly interpreted similar to a photograph. These include medium-resolution sensors like Landsat (30 m) and Sentinel-2 (10 m) and very high-resolution sensors such as Maxar's Worldview-2 (51 cm) and Worldview-3 (31 cm). These long-established satellite systems have been joined more recently by constellations of smaller satellites (so-called “Small Sats”) that offer imagery of comparable spatial and spectral resolution to those operated by Maxar. This rapidly expanding portfolio of earth observation satellites offers the potential for a radical transformation of wildlife research in polar regions, but the sheer volume of data being collected now eclipses our capacity for manual imagery interpretation. To meet this challenge, researchers are now harnessing advances in computer vision that, coupled with improvements in computing capacity, promise to deliver a new era in our ability to monitor polar wildlife.

Global and Regional Snow Cover Decline: 2000–2022

Snow cover affects the global surface energy balance and, with its high albedo, exerts a cooling effect on the Earth’s climate. Decreases in snow cover alter the flow of solar energy from being reflected away from Earth to being absorbed, increasing the Earth’s surface temperature. To gain a global understanding of snow cover change, in situ measurements are too few and far between, so remotely sensed data are needed. This research used the medium-resolution sensor MODIS on the Terra satellite, which has been observing global snow cover almost daily since the year 2000. Here, the MOD10C2 eight-day maximum value composite time series data from February 2000 to March 2023 were analyzed to detect global and regional trends in snow cover extent for the first 23 years of the 21st century. Trends in snow cover change during different time periods (seasons and snow-year) were examined using the Mann—Kendall test and the univariate differencing analysis. Both methods produced similar results. Globally, snow cover declined two to ten times as much as it increased, depending on the season of analysis, and annually, global snow cover decreased 5.12% (not including Antarctica or Greenland) based on the Mann—Kendall test at the 95th percentile (p < 0.05). Regionally, Asia had the greatest net area decline in snow cover, followed by Europe. Although North America has the second-largest extent of snow cover, it had the least amount of net decreasing snow cover relative to its size. South America had the greatest local decline in snow cover, decreasing 20.60% of its annual (snow-year) snow cover area. The Australia–New Zealand region, with just 0.34% of the global snow cover, was the only region to have a net increase in snow cover, increasing 3.61% of its annual snow cover area.

Development of Global Snow Cover—Trends from 23 Years of Global SnowPack

Globally, the seasonal snow cover is the areal largest, the most short-lived and the most variable part of the cryosphere. Remote sensing proved to be a reliable tool to investigate their short-term variations worldwide. The medium-resolution sensor MODIS sensor has been delivering daily snow products since the year 2000. Remaining data gaps due to cloud coverage or polar night are interpolated using the DLR’s Global SnowPack (GSP) processor which produces daily global cloud-free snow cover. With the conclusion of the hydrological year 2022 in the northern hemisphere, the snow cover dynamics of the last 23 hydrological years can now be examined. Trends in snow cover development over different time periods (months, seasons, snow seasons) were examined using the Mann–Kendall test and the Theil–Sen slope. This took place as both pixel based and being averaged over selected hydrological catchment areas. The 23-year time series proved to be sufficient to identify significant developments for large areas. Globally, an average decrease in snow cover duration of −0.44 days/year was recorded for the full hydrological year, even if slight increases in individual months such as November were also found. Likewise, a large proportion of significant trends could also be determined globally at the catchment area level for individual periods. Most drastic developments occurred in March, with an average decrease in snow cover duration by −0.16 days/year. In the catchment area of the river Neman, which drains into the Baltic Sea, there is even a decrease of −0.82 days/year.

Near real-time mapping of burned area by synergizing multiple satellites remote-sensing data

ABSTRACT Wildfires have significant impacts on human lives, critical infrastructures, and Earth’s ecosystems. Accurate and timely information on burned area (BA) affected by wildfires is vital to better understand the drivers of wildfire events, as well as its relevance for biogeochemical cycles, climate, and air quality, and to aid wildfire management. Single satellite data have been used to detect the characteristics of wildfires, retrospectively mapping BAs at a variety of spatial resolutions in previous studies. However, due to the trade-off between spatial and temporal resolutions, single-source satellite data are not sufficient to characterize the explicit dynamics of BAs at high resolutions in both space and time. Thus, a two-stage near real-time BA mapping method was developed in this study to take advantage of the high temporal frequency of coarse resolution sensors and the fine spatial resolution of medium resolution sensors in BA mapping by synergizing freely available coarse and medium spatial resolution (MSR) sensors. First, high temporal frequency sensors such as MODIS and VIIRS were used to identify wildfires and potential BAs. Then, multiple MSR sensors such as Sentinel-2A/2B, Landsat OLI, and Resourcesat AWiFS were synthesized for extracting the BAs with more spatial details in near real-time. We applied the method in California, USA, where wildfires occurred in northern and southern parts in 2017. The results showed that the proposed method is promising for BA mapping with an overall accuracy of 0.84 and 0.85 for wildfires in northern and southern California, respectively. Additionally, the proposed method greatly improved the frequency and reduced the latency, with an average interval of 3.5 days (3 days) and latency of 4 days (sub-daily) for wildfires in southern (northern) California. The extracted BAs illustrated accurate spatial details with MSR sensors. Our method can significantly take advantage of multi-source remote-sensing observations to accurately map the BAs of active wildfires in near real-time. More importantly, the method can be applied to other geographic regions where wildfires risk humans and ecosystems.

Sentinel-2 high-resolution data for river discharge monitoring

River monitoring is an open issue due to many intrinsic problems of the ground monitoring network. Over the last few decades, the role of satellite sensors in river discharge estimation is significantly increased thanks to the strong growth in technologies and applications. Focusing on daily river discharge measurements, a non-linear regression model has been used to link the near-infrared (NIR) reflectance ratio between a dry and a wet pixel around the section of a river to the ground measurements of river discharge. The use of medium-resolution satellite data, such as those from MODIS sensors, enables to monitor high and low flows in medium-sized catchments (<100,000 km2), thanks to satellite frequent revisit time and wide spatial coverage. However, such sensors are not suitable to provide information for medium-narrow rivers (< 250 m wide), nor to study river features and patterns that are averaged within a single pixel. Here, we investigated the use of Sentinel-2 NIR reflectances to support the hypothesis that a higher spatial resolution, i.e. 10 m, is able to better identify the wet pixels, more related to the river dynamics, with obvious advantages for river discharge estimation compared to the medium resolution sensors (e.g., MODIS at 250 m). Moreover, it also allows both a finer distinction between vegetation, soil and water and the characterization of water turbidity in the river area. A new formulation enriched by the sediment component is proposed together with a first step toward an uncalibrated procedure to select the wet pixels. Google Earth Engine (GEE) platform has been employed for the data analysis, allowing to avoid the download of big amounts of data, fostering the reproducibility of the analysis in different locations. The accuracy of the river discharges derived from Sentinel-2 reflectances is evaluated against the in-situ observations from selected gauging stations along two Italian rivers, Po and Tiber. The results confirm the good performances obtained with high-resolution images over the Po River, with average Nash-Sutcliffe efficiency ranging between 0.39 and 0.56 for the different configurations adopted. Relatively worse results were obtained over the Tiber River where the Nash-Sutcliffe efficiency ranged between 0.2 and 0.61, due to an issue on the registration of Sentinel-2 images.

JOINT SUPER-RESOLUTION AND IMAGE RESTORATION FOR PLÉIADES NEO IMAGERY

Abstract. Modern Earth Observation optical satellite systems, such as Airbus’s Pleiades Neo (PNeo) push the boundaries of high spatial resolution by providing commercial imagery products with up to 30cm ground sampling distance (GSD). To further enhance the quality of the images, the in-space imaging system is usually complemented by on-ground image restoration processing, such as deconvolution and denoising (Latry et al., 2012). Recent advances leverage Convolutional Neural Networks (CNNs) to improve the image restoration quality (K. Zhang et al., 2021a).Single Image Super-Resolution (SISR), or Zoom, the process of obtaining a higher resolution (HR) image from a lower resolved (LR) source, has recently gained traction for both medium resolution sensors such as Sentinel 2 (Lanaras et al., 2018) and high resolution such as Pléiades and GeoEye-1 (Zhu et al., 2020). This process further enhances the resolution of the image to improve downstream applications such as mapping (L. Zhang et al., 2021) and small objects recognition (Shermeyer and Van Etten, 2019). While SISR for remote sensing has been successfully tackled using CNNs (Rohith and Kumar, 2021) the main challenge for reaching acceptable image quality performance lies in the generation of realistic LR/HR training pairs (K. Zhang et al., 2021b). In this paper, we propose:a dedicated simulation chain leveraging extremely-high-resolution (EHR) aerial imagery to generate realistic 30cm Pléiades Neo images and their corresponding fully restored HR equivalent at 15cm GSDA residual-based CNN architecture which we train to jointly restore and zoom the images All contributions are assessed on real PNEO images.We deployed the trained models in a production context, to enhance the full Pléiades Neo products – with a swath of 47k pixels – in an efficient and scalable manner.

Improving the spatiotemporal resolution of remotely sensed ET information for water management through Landsat, Sentinel-2, ECOSTRESS and VIIRS data fusion

Robust information on consumptive water use (evapotranspiration, ET) derived from remote sensing can significantly benefit water decision-making in agriculture, informing irrigation schedules and water management plans over extended regions. To be of optimal utility for operational usage, these remote sensing ET data should be generated at the sub-field spatial resolution and daily-to-weekly timesteps commensurate with the scales of water management activities. However, current methods for field-scale ET retrieval based on thermal infrared (TIR) imaging, a valuable diagnostic of canopy stress and surface moisture status, are limited by the temporal revisit of available medium-resolution (100 m or finer) thermal satellite sensors. This study investigates the efficacy of a data fusion method for combining information from multiple medium-resolution sensors toward generating high spatiotemporal resolution ET products for water management. TIR data from Landsat and ECOSTRESS (both at ~ 100-m native resolution), and VIIRS (375-m native) are sharpened to a common 30-m grid using surface reflectance data from the Harmonized Landsat-Sentinel dataset. Periodic 30-m ET retrievals from these combined thermal data sources are fused with daily retrievals from unsharpened VIIRS to generate daily, 30-m ET image timeseries. The accuracy of this mapping method is tested over several irrigated cropping systems in the Central Valley of California in comparison with flux tower observations, including measurements over irrigated vineyards collected in the GRAPEX campaign. Results demonstrate the operational value added by the augmented TIR sensor suite compared to Landsat alone, in terms of capturing daily ET variability and reduced latency for real-time applications. The method also provides means for incorporating new sources of imaging from future planned thermal missions, further improving our ability to map rapid changes in crop water use at field scales.

Evaluation of Satellite-Derived Bathymetry from High and Medium-Resolution Sensors Using Empirical Methods

This study evaluates the accuracy of bathymetric maps generated from multispectral satellite datasets acquired from different multispectral sensors, namely the Worldview 2, PlanetScope, and the Sentinel 2, in the bay of Elounda in Crete. Image pre-processing steps were implemented before the use of the three empirical methods for estimating bathymetry. A dedicated correction and median filter have been applied to minimize noise from the sun glint and the sea waves. Due to the spectral complexity of the selected study area, statistical correlation with different numbers of bands was applied. The analysis indicated that blue and green bands obtained the best results with higher accuracy. Then, three empirical models, namely the Single Band Linear Algorithm, the Multiband Linear Algorithm, and the Ratio Transform Algorithm, were applied to the three multispectral images. Bathymetric and error distribution maps were created and used for the error assessment of results. The accuracy of the bathymetric maps estimated from different empirical models is compared with on-site Single beam Echo Sounder measurements. The most accurate bathymetric maps were obtained using the WorldView 2 and the empirical model of the Ratio Transform algorithm, with the RMSE reaching 1.01 m.

Land use/land cover mapping from airborne hyperspectral images with machine learning algorithms and contextual information

Land use and Land cover (LULC) mapping is one of the most important application areas of remote sensing which requires both spectral and spatial resolutions in order to decrease the spectral ambiguity of different land cover types. Airborne hyperspectral images are among those data which perfectly suits to that kind of applications because of their high number of spectral bands and the ability to see small details on the field. As this technology has newly developed, most of the image processing methods are for the medium resolution sensors and they are not capable of dealing with high resolution images. Therefore, in this study a new framework is proposed to improve the classification accuracy of land use/cover mapping applications and to achieve a greater reliability in the process of mapping land use map using high resolution hyperspectral image data. In order to achieve it, spatial information is incorporated together with spectral information by exploiting feature extraction methods like Grey Level Co-occurrence Matrix (GLCM), Gabor and Morphological Attribute Profile (MAP) on dimensionally reduced image with highest accuracy. Then, machine learning algorithms like Random Forest (RF) and Support Vector Machine (SVM) are used to investigate the contribution of texture information in the classification of high resolution hyperspectral images. In addition to that, further analysis is conducted with object based RF classification to investigate the contribution of contextual information. Finally, overall accuracy, producer’s/user’s accuracy, the quantity and allocation based disagreements and location and quantity based kappa agreements are calculated together with McNemar tests for the accuracy assessment. According to our results, proposed framework which incorporates Gabor texture information and exploits Discrete Wavelet Transform based dimensionality reduction method increase the overall classification accuracy up to 9%. Amongst individual classes, Gabor features boosted classification accuracies of all the classes (soil, road, vegetation, building and shadow) to 7%, 6%, 6%, 8%, 9%, and 24% respectively with producer’s accuracy. Besides, 17% and 10% increase obtained in user’s accuracy with MAP (area) feature in classifying road and shadow classes respectively. Moreover, when the object based classification is conducted, it is seen that the OA of pixel based classification is increased further by 1.07%. An increase between 2% and 4% is achieved with producer’s accuracy in soil, vegetation and building classes and an increase between 1% and 3% is achieved by user’s accuracy in soil, road, vegetation and shadow classes. In the end, accurate LULC map is produced with object based RF classification of gabor features added airborne hyperspectral image which is dimensionally reduced with DWT method.

A new method for monitoring start of season (SOS) of forest based on multisource remote sensing

As a sensitive indicator of climate change, forest phenology (e.g., the start of season (SOS)) has profound impacts on the global carbon cycle. Traditional phenological observations are based on surface observation networks. Generally, the outcomes of this manner are less representative and hard to implement in wide forested regions. Remote sensing based observation of forest SOS has currently been a popular way. Those sensors with coarse spatial resolution have been widely used to estimate forest SOS, but they create serious estimation errors in areas of high heterogeneity. Medium-resolution sensors, such as Landsat, face significant challenges in SOS monitoring due to the long revisit period. In this study, we aimed to develop a new method to estimate forest SOS from 2013 to 2019. First, we collected all available Landsat and Sentinel-2 images, and then redefined the linear regression coefficients for the bandpass adjustment to weaken the surface reflectance (SR) differences in different sensors. Subsequently, we improved and developed the modified continuous change detection and classification (MCCDC) model to generate daily vegetation index curves. Finally, we adopted the logistic regression model to test the potential of the enhanced vegetation index (EVI), normalized difference vegetation index (NDVI) and land surface water index (LSWI) in evaluating the annual SOS. The reduced root mean square error (RMSE) for all bands after the integration indicated that the adjustment was successful. We visually compared Landsat’s synthetic images with the actual acquired images and found that their respective false-colour composites were highly similar. Assessing the SOS from the EVI, NDVI and LSWI showed different estimated results. By comparing the annual SOS derived from the three indices with the field observations, it was found that the SOS based on the EVI maintained a low consistency with the field observations. The SOS accuracy from LSWI was the highest and most forest SOS from LSWI in the study area were mainly concentrated in 80–150 days. These three indices all showed that the SOS always fluctuated by ± 4 days from 2013 to 2019. Facing the lack of clear remote sensing images with medium spatial resolution in cloudy and rainy areas, this study proposes an improved method to generate clear daily vegetation index images at a spatial resolution of 30 m, making annual SOS monitoring promising and feasible.

Potential and Challenges of Harmonizing 40 Years of AVHRR Data: The TIMELINE Experience

Earth Observation satellite data allows for the monitoring of the surface of our planet at predefined intervals covering large areas. However, there is only one medium resolution sensor family in orbit that enables an observation time span of 40 and more years at a daily repeat interval. This is the AVHRR sensor family. If we want to investigate the long-term impacts of climate change on our environment, we can only do so based on data that remains available for several decades. If we then want to investigate processes with respect to climate change, we need very high temporal resolution enabling the generation of long-term time series and the derivation of related statistical parameters such as mean, variability, anomalies, and trends. The challenges to generating a well calibrated and harmonized 40-year-long time series based on AVHRR sensor data flown on 14 different platforms are enormous. However, only extremely thorough pre-processing and harmonization ensures that trends found in the data are real trends and not sensor-related (or other) artefacts. The generation of European-wide time series as a basis for the derivation of a multitude of parameters is therefore an extremely challenging task, the details of which are presented in this paper.

Land use/land cover mapping from airborne hyperspectral images with machine learning algorithms and contextual information

Land use and Land cover (LULC) mapping is one of the most important application areas of remote sensing which requires both spectral and spatial resolutions in order to decrease the spectral ambiguity of different land cover types. Airborne hyperspectral images are among those data which perfectly suits to that kind of applications because of their high number of spectral bands and the ability to see small details on the field. As this technology has newly developed, most of the image processing methods are for the medium resolution sensors and they are not capable of dealing with high resolution images. Therefore, in this study a new framework is proposed to improve the classification accuracy of land use/cover mapping applications and to achieve a greater reliability in the process of mapping land use map using high resolution hyperspectral image data. In order to achieve it, spatial information is incorporated together with spectral information by exploiting feature extraction methods like Grey Level Co-occurrence Matrix (GLCM), Gabor and Morphological Attribute Profile (MAP) on dimensionally reduced image with highest accuracy. Then, machine learning algorithms like Random Forest (RF) and Support Vector Machine (SVM) are used to investigate the contribution of texture information in the classification of high resolution hyperspectral images. In addition to that, further analysis is conducted with object based RF classification to investigate the contribution of contextual information. Finally, overall accuracy, producer’s/user’s accuracy, the quantity and allocation based disagreements and location and quantity based kappa agreements are calculated together with McNemar tests for the accuracy assessment. According to our results, proposed framework which incorporates Gabor texture information and exploits Discrete Wavelet Transform based dimensionality reduction method increase the overall classification accuracy up to 9%. Amongst individual classes, Gabor features boosted classification accuracies of all the classes (soil, road, vegetation, building and shadow) to 7%, 6%, 6%, 8%, 9%, and 24% respectively with producer’s accuracy. Besides, 17% and 10% increase obtained in user’s accuracy with MAP(area) feature in classifying road and shadow classes respectively. Moreover, when the object based classification is conducted, it is seen that the OA of pixel based classification is increased further by 1.07%. An increase between 2% and 4% is achieved with producer’s accuracy in soil, vegetation and building classes and an increase between 1% and 3% is achieved by user’s accuracy in soil, road, vegetation and shadow classes. In the end, accurate LULC map is produced with object based RF classification of gabor features added airborne hyperspectral image which is dimensionally reduced with DWT method.

EVALUATING THE IMPACT OF LASRC AND SEN2COR ATMOSPHERIC CORRECTION ALGORITHMS ON LANDSAT-8/OLI AND SENTINEL-2/MSI DATA OVER AERONET STATIONS IN BRAZILIAN TERRITORY

Abstract. Accurate and consistent Surface Reflectance estimation from optical remote sensor observations is directly dependant on the used atmospheric correction processor and the differences caused by it may have implications on further processes, e.g. classification. Brazil is a continental scale country with different biomes. Recently, new initiatives, as the Brazil Data Cube Project, are emerging and using free and open data policy data, more specifically medium spatial resolution sensor images, to create image data cubes and classify the Brazilian territory crops. For this reason, the purpose of this study is to verify, on Landsat-8 and Sentinel-2 images for the Brazilian territory, the suitability of the atmospheric correction processors maintained by their image providers, LaSRC from USGS and Sen2cor from ESA, respectively. To achieve this, we tested the surface reflectance products from Landsat-8 processed through LaSRC and Sentinel-2 processed through LaSRC and Sen2cor comparing to a reference dataset computed by ARCSI and AERONET. The obtained results point that Landsat-8/OLI images atmospherically corrected using the LaSRC corrector are consistent to the surface reflectance reference and other atmospheric correction processors studies, while for Sentinel-2/MSI images, Sen2cor performed best. Although corrections over Sentinel-2/MSI data weren’t as consistent as in Landsat-8/OLI corrections, in comparison to the surface reflectance references, most of the spectral bands achieved acceptable APU results.

Atmospheric and sunglint correction for retrieving chlorophyll-a in a productive tropical estuarine-lagoon system using Sentinel-2 MSI imagery

Remote monitoring of chlorophyll-a (chla) has been widely used to evaluate the trophic state of inland and coastal waters, however, there is still much uncertainty in the algorithms applied in different optical water types. The influence of different atmospheric correction (AC) processors, which can also provide correction for sunglint and adjacency effects, on the retrieved chla is poorly understood. In this study, state-of-the-art atmospheric correction and chla algorithms are evaluated using Sentinel-2 MSI imagery in the Mundaú-Manguaba Estuarine-Lagoon System (MMELS), a productive tropical system that consists of two turbid lagoons of different optical water types (OWT). We compared the performance of six AC processors, with the addition of sunglint correction for two of them, with field measured water reflectance. There was difficulty in correcting for the atmospheric effects, especially for bands 2, 3 and 8A. Overall, C2X showed the best performance over MMELS, but with sunglint correction, ACOLITE and GRS provided the most consistent water reflectance (ρw). Sunglint correction might be essential for retrieving accurate ρw in most low-latitude water bodies. We also found that in Mundaú, the dense urban area surrounding it likely caused heavy adjacency effects in the satellite-retrieved reflectance, and thus correction for it is necessary. We also compared the performance of six chla algorithms recommended for the OWTs present in MMELS in addition to a widely applied and a global chla algorithm in retrieving this variable using both field and satellite reflectance, in this case corrected with the three best performing processors. For the in situ data, most algorithms performed well in Manguaba lagoon, while in Mundaú lagoon the semi-analytical NIR-red ratio (2SAR) algorithm was the most consistent model, and in both cases the locally calibrated algorithms outperformed the global algorithm. When retrieving chla with the satellite-derived ρw, considerably poorer results were produced, especially in Mundaú lagoon. The global algorithm was found to be especially sensitive to the atmospheric effects. We also found that the quality of AC provided by the algorithms is not a general predictor of the performance of the chla models, even when analysing individual bands separately, while the relationship between chla concentration and the ratio of bands used by most algorithms can be. Despite containing distinct water characteristics, chla can be modelled using a single algorithm, 2SAR, calibrated for MMELS as a whole, with r2 of 0.77 and nRMSE of 38.7%, and we consider that 2SAR has the potential to be a global algorithm for these OWTs, provided that it is recalibrated for a large dataset of satellite-derived BOA reflectance. We recommend that further studies explore the impacts of AC, sunglint and adjacency effects on the performance of chla algorithms, in order to delineate the most suitable combinations of AC + chla algorithms for the variable OWTs, in an effort to provide the basis for global-scale retrievals of this pigment using medium-resolution sensors such as MSI and OLI.

Decameter Cropland LAI/FPAR Estimation From Sentinel-2 Imagery Using Google Earth Engine

Leaf area index (LAI) and fraction of photosynthetically active radiation (FPAR) products at regional and global scales have already been extensively and routinely generated from medium-resolution sensors. However, there is a lack of high-resolution LAI/FPAR product, which is especially essential for crop growth and drought monitoring of cropland in patches. This article proposes a processing framework for the derivation of decameter cropland LAI and FPAR in the Northern China plain from Sentinel-2 surface reflectance data with a random forest (RF) algorithm by exploiting the capabilities of the Google Earth Engine (GEE) cloud platform. The training database is generated from the spatially aggregated Sentinel-2 surface reflectance and the corresponding Moderate Resolution Imaging Spectroradiometer (MODIS) LAI/FPAR product over homogeneous cropland, and the training samples are strictly filtered for the best quality. RF is then trained over the processed Sentinel-2 surface reflectance and the filtered MODIS LAI/FPAR under two input groups—one group is for Sentinel-2 spectral bands of 10-m resolution only, and the other group supplements the Sentinel-2 red-edge (RE) and shortwave infrared (SWIR) bands of 20-m resolution. Extensive comparisons and validation are carried out, and they demonstrate that the new method can generate spatial and temporal consistent LAI/FPAR with MODIS at high spatial resolution. The retrieval accuracy is slightly better for 20-m input groups than that for 10-m input groups, confirming the value of RE and/or SWIR in cropland LAI/FPAR estimate. This article also demonstrates that GEE is a suitable high-performance processing tool for high-resolution biophysical variables estimation.

Surface Albedo Retrieval from 40-Years of Earth Observations through the EUMETSAT/LSA SAF and EU/C3S Programmes: The Versatile Algorithm of PYALUS

Land surface albedo quantifies the fraction of the sunlight reflected by the surface of the Earth. This article presents the algorithm concepts for the remote sensing of this variable based on the heritage of several developments which were performed at Méteo France over the last decade and described in several papers by Carrer et al. The scientific algorithm comprises four steps: an atmospheric correction, a sensor harmonisation (optional), a BRDF (Bidirectional Reflectance Distribution Function) inversion, and the albedo calculation. At the time being, the method has been applied to 11 sensors in the framework of two European initiatives (Satellite Application Facility on Land Surface Analysis—LSA SAF, and Copernicus Climate Change Service—C3S): NOAA-7-9-11-14-16-17/AVHRR2-3, SPOT/VGT1-2, Metop/AVHRR-3, PROBA-V, and MSG/SEVIRI. This work leads to a consistent archive of almost 40 years of satellite-derived albedo data (available in 2020). From a single sensor, up to three different albedo products with different characteristics have been developed to address the requirements of both, near real-time (NRT) (weather prediction with a demand of timeliness of 1 h) and climate communities. The evaluation of the algorithm applied to different platforms was recently made by Lellouch et al. and Sánchez Zapero et al. in 2020 which can be considered as companion papers. After a summary of the method for the retrieval of these surface albedos, this article describes the specificities of each retrieval, lists the differences, and discusses the limitations. The plan of continuity with the next European satellite missions and perspectives of improvements are introduced. For example, Metop/AVHRR-3 albedo will soon become the medium resolution sensor product with the longest NRT data record, since MODIS is approaching the end of its life-cycle. Additionally, Metop-SG/METimage will ensure its continuity thanks to consistent production of data sets guaranteed till 2050 by the member states of the European Organisation for the Exploitation of Meteorological Satellites (EUMETSAT). In the end, the common strategy which we proposed through the different programmes may offer an unprecedented opportunity to study the temporal trends affecting surface properties and to analyse human-induced climate change. Finally, the access to the source code (called PYALUS) is provided through an open access platform in order to share with the community the expertise on the satellite retrieval of this variable.

BURNT AREA DETECTION USING MEDIUM RESOLUTION SENTINEL 2 AND LANDSAT 8 SATELLITES

Abstract. Forest fire is one of the most serious environmental problems in Kenya that influences human activities, climate change and biodiversity. The main goal of this study is to apply medium resolution sensors (Landsat 8 OLI and Sentinel 2 MSI) to produce burnt area severity maps that will include small fires (&lt; 100 ha) in order to improve burnt area detection and mapping in Kenya. Normalized burnt area indices were generated for specified pre- and post-fire periods. The difference between pre- and post-fire Normalized Burnt Ration (NBR) was used to compute δNBR index depicting forest disturbance by fire events. Thresholded classes were derived from the computed δNBR indices to obtain burnt severity maps. The spatial and temporal agreements of the Burnt area detection dates were validated by comparing against the MODIS MCD641 500 m products and MODIS Fire Information for Resource Management System (FIRMS) 1 km daily product hot-spot acquisition dates. This approach was implemented on Google Earth Engine (GEE) platform with a simple user interface that allows users to auto-generate burnt area maps and statistics. The operational GEE application developed can be used to obtain burnt area severity maps and statistics that allow for initial accurate approximation of fire damage.

Analysis of Spatial and Temporal Variability in Libya-4 with Landsat 8 and Sentinel-2 Data for Optimized Ground Target Location

Pseudo-Invariant Calibration Sites (PICS) have been widely used by the remote sensing community in recent decades for post-launch absolute calibration, cross-calibration, and the monitoring of radiometric stability. The Committee on Earth Observation Satellites (CEOS) has established several official PICS for these purposes. Of these, Libya-4 is the most commonly used, due to its high uniformity and stability. The site was chosen as a large-area site for medium resolution sensors, and with high-resolution sensors now common, smaller sites are being identified. This work has identified an improved area of interest (AOI) within Libya-4 by using combined Landsat 8 and Sentinel 2 data. The Optimized Ground Target (OGT) was determined by calculating the coefficient of variation along with the use of a quasi-Newton optimization algorithm combined with the Basin–Hopping global optimization technique to constrain a search area small enough to perform a final brute-force refinement. The Coefficient of Variation CV of the proposed OGT is significantly lower than that in the original CEOS area, with differences between the CV of both zones in the order of 1% in the visible near-infrared (VNIR) bands. This new AOI has the potential to improve the cross-calibration between high-resolution sensors using the PICS methodology through an OGT with more homogeneous and stable characteristics.

Retrieving Total and Inorganic Suspended Sediments in Amazon Floodplain Lakes: A Multisensor Approach

Remote sensing imagery are fundamental to increasing the knowledge about sediment dynamics in the middle-lower Amazon floodplains. Moreover, they can help to understand both how climate change and how land use and land cover changes impact the sediment exchange between the Amazon River and floodplain lakes in this important and complex ecosystem. This study investigates the suitability of Landsat-8 and Sentinel-2 spectral characteristics in retrieving total (TSS) and inorganic (TSI) suspended sediments on a set of Amazon floodplain lakes in the middle-lower Amazon basin using in situ Remote Sensing Reflectance (Rrs) measurements to simulate Landsat 8/OLI (Operational Land Imager) and Sentinel 2/MSI (Multispectral Instrument) bands and to calibrate/validate several TSS and TSI empirical algorithms. The calibration was based on the Monte Carlo Simulation carried out for the following datasets: (1) All-Dataset, consisting of all the data acquired during four field campaigns at five lakes spread over the lower Amazon floodplain (n = 94); (2) Campaign-Dataset including samples acquired in a specific hydrograph phase (season) in all lakes. As sample size varied from one season to the other, n varied from 18 to 31; (3) Lake-Dataset including samples acquired in all seasons at a given lake with n also varying from 17 to 67 for each lake. The calibrated models were, then, applied to OLI and MSI scenes acquired in August 2017. The performance of three atmospheric correction algorithms was also assessed for both OLI (6S, ACOLITE, and L8SR) and MSI (6S, ACOLITE, and Sen2Cor) images. The impact of glint correction on atmosphere-corrected image performance was assessed against in situ glint-corrected Rrs measurements. After glint correction, the L8SR and 6S atmospheric correction performed better with the OLI and MSI sensors, respectively (Mean Absolute Percentage Error (MAPE) = 16.68% and 14.38%) considering the entire set of bands. However, for a given single band, different methods have different performances. The validated TSI and TSS satellite estimates showed that both in situ TSI and TSS algorithms provided reliable estimates, having the best results for the green OLI band (561 nm) and MSI red-edge band (705 nm) (MAPE &lt; 21%). Moreover, the findings indicate that the OLI and MSI models provided similar errors, which support the use of both sensors as a virtual constellation for the TSS and TSI estimate over an Amazon floodplain. These results demonstrate the applicability of the calibration/validation techniques developed for the empirical modeling of suspended sediments in lower Amazon floodplain lakes using medium-resolution sensors.