FORCE—Landsat + Sentinel-2 Analysis Ready Data and Beyond

The devil is in the detail: Environmental variables frequently used for habitat suitability modeling lack information for forest-dwelling bats in Germany

In response to the pressing challenges of the ongoing biodiversity crisis, the protection of endangered species and their habitats, as well as the monitoring of invasive species are crucial. Habitat suitability modeling (HSM) is often treated as the silver bullet to address these challenges, commonly relying on generic variables sourced from widely available datasets. However, for species with high habitat requirements, or for modeling the suitability of habitats within the geographic range of a species, variables at a coarse level of detail may fall short. Consequently, there is potential value in considering the incorporation of more targeted data, which may extend beyond readily available land cover and climate datasets. In this study, we investigate the impact of incorporating targeted land cover variables (specifically tree species composition) and vertical structure information (derived from LiDAR data) on HSM outcomes for three forest specialist bat species (Barbastella barbastellus, Myotis bechsteinii, and Plecotus auritus) in Rhineland-Palatinate, Germany, compared to commonly utilized environmental variables, such as generic land-cover classifications (e.g., Corine Land Cover) and climate variables (e.g., Bioclim). The integration of targeted variables enhanced the performance of habitat suitability models for all three bat species. Furthermore, our results showed a high difference in the distribution maps that resulted from using different levels of detail in environmental variables. This underscores the importance of making the effort to generate the appropriate variables, rather than simply relying on commonly used ones, and the necessity of exercising caution when using habitat models as a tool to inform conservation strategies and spatial planning efforts.

Multispectral analysis-ready satellite data for three East African mountain ecosystems

The East African mountain ecosystems are facing increasing threats due to global change, putting their unique socio-ecological systems at risk. To monitor and understand these changes, researchers and stakeholders require accessible analysis-ready remote sensing data. Although satellite data is available for many applications, it often lacks accurate geometric orientation and has extensive cloud cover. This can generate misleading results and make it unreliable for time-series analysis. Therefore, it needs comprehensive processing before usage, which encompasses multi-step operations, requiring large computational and storage capacities, as well as expert knowledge. Here, we provide high-quality, atmospherically corrected, and cloud-free analysis-ready Sentinel-2 imagery for the Bale Mountains (Ethiopia), Mounts Kilimanjaro and Meru (Tanzania) ecosystems in East Africa. Our dataset ranges from 2017 to 2021 and is provided as monthly and annual aggregated products together with 24 spectral indices. Our dataset enables researchers and stakeholders to conduct immediate and impactful analyses. These applications can include vegetation mapping, wildlife habitat assessment, land cover change detection, ecosystem monitoring, and climate change research.

Topographic barriers drive the pronounced genetic subdivision of a range-limited fossorial rodent

Due to their limited dispersal ability, fossorial species with predominantly belowground activity usually show increased levels of population subdivision across relatively small spatial scales. This may be exacerbated in harsh mountain ecosystems, where landscape geomorphology limits species' dispersal ability and leads to small effective population sizes, making species relatively vulnerable to environmental change. To better understand the environmental drivers of species' population subdivision in remote mountain ecosystems, particularly in understudied high-elevation systems in Africa, we studied the giant root-rat (Tachyoryctes macrocephalus), a fossorial rodent confined to the afro-alpine ecosystem of the Bale Mountains in Ethiopia. Using mitochondrial and low-coverage nuclear genomes, we investigated 77 giant root-rat individuals sampled from nine localities across its entire ~1000 km2 range. Our data revealed a distinct division into a northern and southern group, with no signs of gene flow, and higher nuclear genetic diversity in the south. Landscape genetic analyses of the mitochondrial and nuclear genomes indicated that population subdivision was driven by slope and elevation differences of up to 500 m across escarpments separating the north and south, potentially reinforced by glaciation of the south during the Late Pleistocene (~42,000-16,000 years ago). Despite this landscape-scale subdivision between the north and south, weak geographic structuring of sampling localities within regions indicated gene flow across distances of at least 16 km at the local scale, suggesting high, aboveground mobility for relatively long distances. Our study highlights that despite the potential for local-scale gene flow in fossorial species, topographic barriers can result in pronounced genetic subdivision. These factors can reduce genetic variability, which should be considered when developing conservation strategies.

Unveiling patterns in human dominated landscapes through mapping the mass of US built structures

Built structures increasingly dominate the Earth's landscapes; their surging mass is currently overtaking global biomass. We here assess built structures in the conterminous US by quantifying the mass of 14 stock-building materials in eight building types and nine types of mobility infrastructures. Our high-resolution maps reveal that built structures have become 2.6 times heavier than all plant biomass across the country and that most inhabited areas are mass-dominated by buildings or infrastructure. We analyze determinants of the material intensity and show that densely built settlements have substantially lower per-capita material stocks, while highest intensities are found in sparsely populated regions due to ubiquitous infrastructures. Out-migration aggravates already high intensities in rural areas as people leave while built structures remain - highlighting that quantifying the distribution of built-up mass at high resolution is an essential contribution to understanding the biophysical basis of societies, and to inform strategies to design more resource-efficient settlements and a sustainable circular economy.

The global wildland-urban interface

The wildland-urban interface (WUI) is where buildings and wildland vegetation meet or intermingle1,2. It is where human-environmental conflicts and risks can be concentrated, including the loss of houses and lives to wildfire, habitat loss and fragmentation and the spread of zoonotic diseases3. However, a global analysis of the WUI has been lacking. Here, we present a global map of the 2020 WUI at 10 m resolution using a globally consistent and validated approach based on remote sensing-derived datasets of building area4 and wildland vegetation5. We show that the WUI is a global phenomenon, identify many previously undocumented WUI hotspots and highlight the wide range of population density, land cover types and biomass levels in different parts of the global WUI. The WUI covers only 4.7% of the land surface but is home to nearly half its population (3.5 billion). The WUI is especially widespread in Europe (15% of the land area) and the temperate broadleaf and mixed forests biome (18%). Of all people living near 2003-2020 wildfires (0.4 billion), two thirds have their home in the WUI, most of them in Africa (150 million). Given that wildfire activity is predicted to increase because of climate change in many regions6, there is a need to understand housing growth and vegetation patterns as drivers of WUI change.

Ecodatacube.eu: analysis-ready open environmental data cube for Europe

The article describes the production steps and accuracy assessment of an analysis-ready, open-access European data cube consisting of 2000-2020+ Landsat data, 2017-2021+ Sentinel-2 data and a 30 m resolution digital terrain model (DTM). The main purpose of the data cube is to make annual continental-scale spatiotemporal machine learning tasks accessible to a wider user base by providing a spatially and temporally consistent multidimensional feature space. This has required systematic spatiotemporal harmonization, efficient compression, and imputation of missing values. Sentinel-2 and Landsat reflectance values were aggregated into four quarterly averages approximating the four seasons common in Europe (winter, spring, summer and autumn), as well as the 25th and 75th percentile, in order to retain intra-seasonal variance. Remaining missing data in the Landsat time-series was imputed with a temporal moving window median (TMWM) approach. An accuracy assessment shows TMWM performs relatively better in Southern Europe and lower in mountainous regions such as the Scandinavian Mountains, the Alps, and the Pyrenees. We quantify the usability of the different component data sets for spatiotemporal machine learning tasks with a series of land cover classification experiments, which show that models utilizing the full feature space (30 m DTM, 30 m Landsat, 30 m and 10 m Sentinel-2) yield the highest land cover classification accuracy, with different data sets improving the results for different land cover classes. The data sets presented in the article are part of the EcoDataCube platform, which also hosts open vegetation, soil, and land use/land cover (LULC) maps created. All data sets are available under CC-BY license as Cloud-Optimized GeoTIFFs (ca. 12 TB in size) through SpatioTemporal Asset Catalog (STAC) and the EcoDataCube data portal.

High-resolution data and maps of material stock, population, and employment in Austria from 1985 to 2018

High-resolution maps of material stocks in buildings and infrastructures are of key importance for studies of societal resource use (social metabolism, circular economy, secondary resource potentials) as well as for transport studies and land system science. So far, such maps were only available for specific years but not in time series. Even for single years, data covering entire countries with high resolution, or using remote-sensing data are rare. Instead, they often have local extent (e.g., [1]), are lower resolution (e.g., [2]), or are based on other geospatial data (e.g., [3]). We here present data on the material stocks in three types of buildings (commercial and industrial, single- and multifamily houses) and three types of infrastructures (roads, railways, other infrastructures) for a 33-year time series for Austria at a spatial resolution of 30 m. The article also presents data on population and employment in Austria for the same time period, at the same spatial resolution. Data were derived with the same method applied in a recent study for Germany [4].

Seeing the System from Above: The Use and Potential of Remote Sensing for Studying Ecosystem Dynamics

Remote sensing techniques are increasingly used for studying ecosystem dynamics, delivering spatially explicit information on the properties of Earth over large spatial and multi-decadal temporal extents. Yet, there is still a gap between the more technology-driven development of novel remote sensing techniques and their applications for studying ecosystem dynamics. Here, I review the existing literature to explore how addressing these gaps might enable recent methods to overcome longstanding challenges in ecological research. First, I trace the emergence of remote sensing as a major tool for understanding ecosystem dynamics. Second, I examine recent developments in the field of remote sensing that are of particular importance for studying ecosystem dynamics. Third, I consider opportunities and challenges for emerging open data and software policies and suggest that remote sensing is at its most powerful when it is theoretically motivated and rigorously ground-truthed. I close with an outlook on four exciting new research frontiers that will define remote sensing ecology in the upcoming decade.

A Newly Developed Algorithm for Cloud Shadow Detection—TIP Method

The masking of cloud shadows in optical satellite imagery is an important step in automated processing chains. A new method (the TIP method) for cloud shadow detection in multi-spectral satellite images is presented and compared to current methods. The TIP method is based on the evaluation of thresholds, indices and projections. Most state-of-the-art methods solemnly rely on one of these evaluation steps or on a complex working mechanism. Instead, the new method incorporates three basic evaluation steps into one algorithm for easy and accurate cloud shadow detection. Furthermore the performance of the masking algorithms provided by the software packages ATCOR (“Atmospheric Correction”) and PACO (“Python-based Atmospheric Correction”) is compared with that of the newly implemented TIP method on a set of 20 Sentinel-2 scenes distributed over the globe, covering a wide variety of environments and climates. The algorithms incorporated in each piece of masking software use the class of cloud shadows, but they employ different rules and class-specific thresholds. Classification results are compared to the assessment of an expert human interpreter. The class assignment of the human interpreter is considered as reference or “truth”. The overall accuracies for the class cloud shadows of ATCOR and PACO (including TIP) for difference areas of the selected scenes are 70.4% and 76.6% respectively. The difference area encompasses the parts of the classification image where the classification maps disagree. User and producer accuracies for the class cloud shadow are strongly scene-dependent, typically varying between 45% and 95%. The experimental results show that the proposed TIP method based on thresholds, indices and projections can obtain improved cloud shadow detection performance.

Using machine learning to produce a cost-effective national building height map of Ireland to categorise local climate zones

Abstract. ECOCLIMAP-Second Generation (ECO-SG) is the land-cover map used in the HARMONIE-AROME configuration of the shared ALADIN-HIRLAM Numerical Weather Prediction system used for short-range operational weather forecasting for Ireland. The ECO-SG urban classification implicitly includes building heights. The work presented in this paper involved the production of the first open-access building height map for the island of Ireland which complements the Ulmas-Walsh land cover map, a map which has improved the horizontal extent of urban areas over Ireland. The resulting building height map will potentially enable upgrades to ECO-SG urban information for future implementation in HARMONIE-AROME. This study not only produced the first open-access building height map of Ireland at 10 m × 10 m resolution, but assessed various types of regression models trained using pre-existing building height information for Dublin City and selected 64 important spatio-temporal features, engineered from both the Sentinel-1A/B and Sentinel-2A/B satellites. The performance metrics revealed that a Convolutional Neural Network is superior in all aspects except the computational time required to create the map. Despite the superior accuracy of the Convolutional Neural Network, the final building height map created results from the ridge regression model which provided the best blend of realistic output and low computational complexity. The method relies solely on freely available satellite imagery, is cost-effective, can be updated regularly, and can be applied to other regions depending on the availability of representative regional building height sample data.

Land Use and Land Cover Mapping Using Sentinel-2, Landsat-8 Satellite Images, and Google Earth Engine: A Comparison of Two Composition Methods

Accurate and real-time land use/land cover (LULC) maps are important to provide precise information for dynamic monitoring, planning, and management of the Earth. With the advent of cloud computing platforms, time series feature extraction techniques, and machine learning classifiers, new opportunities are arising in more accurate and large-scale LULC mapping. In this study, we aimed at finding out how two composition methods and spectral–temporal metrics extracted from satellite time series can affect the ability of a machine learning classifier to produce accurate LULC maps. We used the Google Earth Engine (GEE) cloud computing platform to create cloud-free Sentinel-2 (S-2) and Landsat-8 (L-8) time series over the Tehran Province (Iran) as of 2020. Two composition methods, namely, seasonal composites and percentiles metrics, were used to define four datasets based on satellite time series, vegetation indices, and topographic layers. The random forest classifier was used in LULC classification and for identifying the most important variables. Accuracy assessment results showed that the S-2 outperformed the L-8 spectral–temporal metrics at the overall and class level. Moreover, the comparison of composition methods indicated that seasonal composites outperformed percentile metrics in both S-2 and L-8 time series. At the class level, the improved performance of seasonal composites was related to their ability to provide better information about the phenological variation of different LULC classes. Finally, we conclude that this methodology can produce LULC maps based on cloud computing GEE in an accurate and fast way and can be used in large-scale LULC mapping.

CubeSat constellations provide enhanced crop phenology and digital agricultural insights using daily leaf area index retrievals

Satellite remote sensing has great potential to deliver on the promise of a data-driven agricultural revolution, with emerging space-based platforms providing spatiotemporal insights into precision-level attributes such as crop water use, vegetation health and condition and crop response to management practices. Using a harmonized collection of high-resolution Planet CubeSat, Sentinel-2, Landsat-8 and additional coarser resolution imagery from MODIS and VIIRS, we exploit a multi-satellite data fusion and machine learning approach to deliver a radiometrically calibrated and gap-filled time-series of daily leaf area index (LAI) at an unprecedented spatial resolution of 3 m. The insights available from such high-resolution CubeSat-based LAI data are demonstrated through tracking the growth cycle of a maize crop and identifying observable within-field spatial and temporal variations across key phenological stages. Daily LAI retrievals peaked at the tasseling stage, demonstrating their value for fertilizer and irrigation scheduling. An evaluation of satellite-based retrievals against field-measured LAI data collected from both rain-fed and irrigated fields shows high correlation and captures the spatiotemporal development of intra- and inter-field variations. Novel agricultural insights related to individual vegetative and reproductive growth stages were obtained, showcasing the capacity for new high-resolution CubeSat platforms to deliver actionable intelligence for precision agricultural and related applications.

A Flexible Multi-Temporal and Multi-Modal Framework for Sentinel-1 and Sentinel-2 Analysis Ready Data

The rich, complementary data provided by Sentinel-1 and Sentinel-2 satellite constellations host considerable potential to transform Earth observation (EO) applications. However, a substantial amount of effort and infrastructure is still required for the generation of analysis-ready data (ARD) from the low-level products provided by the European Space Agency (ESA). Here, a flexible Python framework able to generate a range of consistent ARD aligned with the ESA-recommended processing pipeline is detailed. Sentinel-1 Synthetic Aperture Radar (SAR) data are radiometrically calibrated, speckle-filtered and terrain-corrected, and Sentinel-2 multi-spectral data resampled in order to harmonise the spatial resolution between the two streams and to allow stacking with multiple scene classification masks. The global coverage and flexibility of the framework allows users to define a specific region of interest (ROI) and time window to create geo-referenced Sentinel-1 and Sentinel-2 images, or a combination of both with closest temporal alignment. The framework can be applied to any location and is user-centric and versatile in generating multi-modal and multi-temporal ARD. Finally, the framework handles automatically the inherent challenges in processing Sentinel data, such as boundary regions with missing values within Sentinel-1 and the filtering of Sentinel-2 scenes based on ROI cloud coverage.

Extension of an Open GEOBIA Framework for Spatially Explicit Forest Stratification with Sentinel-2

Spatially explicit information about forest cover is fundamental for operational forest management and forest monitoring. Although open-satellite-based earth observation data in a spatially high resolution (i.e., Sentinel-2, ≤10 m) can cover some information needs, spatially very high-resolution imagery (i.e., aerial imagery, ≤2 m) is needed to generate maps at a scale suitable for regional and local applications. In this study, we present the development, implementation, and evaluation of a Geographic Object-Based Image Analysis (GEOBIA) framework to stratify forests (needleleaf, broadleaved, non-forest) in Luxembourg. The framework is exclusively based on open data and free and open-source geospatial software. Although aerial imagery is used to derive image objects with a 0.05 ha minimum size, Sentinel-2 scenes of 2020 are the basis for random forest classifications in different single-date and multi-temporal feature setups. These setups are compared with each other and used to evaluate the framework against classifications based on features derived from aerial imagery. The highest overall accuracies (89.3%) have been achieved with classification on a Sentinel-2-based vegetation index time series (n = 8). Similar accuracies have been achieved with classification based on two (88.9%) or three (89.1%) Sentinel-2 scenes in the greening phase of broadleaved forests. A classification based on color infrared aerial imagery and derived texture measures only achieved an accuracy of 74.5%. The integration of the texture measures into the Sentinel-2-based classification did not improve its accuracy. Our results indicate that high resolution image objects can successfully be stratified based on lower spatial resolution Sentinel-2 single-date and multi-temporal features, and that those setups outperform classifications based on aerial imagery only. The conceptual framework of spatially high-resolution image objects enriched with features from lower resolution imagery facilitates the delivery of frequent and reliable updates due to higher spectral and temporal resolution. The framework additionally holds the potential to derive additional information layers (i.e., forest disturbance) as derivatives of the features attached to the image objects, thus providing up-to-date information on the state of observed forests.

Earth Observation for Phenological Metrics (EO4PM): Temporal Discriminant to Characterize Forest Ecosystems

The study of vegetation phenology has great relevance in many fields since the importance of knowing timing and shifts in periodic plant life cycle events to face the consequences of global changes in issues such as crop production, forest management, ecosystem disturbances, and human health. The availability of high spatial resolution and dense revisit time satellite observations, such as Sentinel-2 satellites, allows high resolution phenological metrics to be estimated, able to provide key information from time series and to discriminate vegetation typologies. This paper presents an automated and transferable procedure that combines validated methodologies based on local curve fitting and local derivatives to exploit full satellite Earth observation time series to produce information about plant phenology. Multivariate statistical analysis is performed for the purpose of demonstrating the capacity of the generated smoothed vegetation curve, temporal statistics, and phenological metrics to serve as temporal discriminants to detect forest ecosystems processes responses to environmental gradients. The results show smoothed vegetation curve and temporal statistics able to highlight seasonal gradient and leaf type characteristics to discriminate forest types, with additional information about forest and leaf productivity provided by temporal statistics analysis. Furthermore, temporal, altitudinal, and latitudinal gradients are obtained from phenological metrics analysis, which also allows to associate temporal gradient with specific phenophases that support forest types distinction. This study highlights the importance of integrated data and methodologies to support the processes of vegetation recognition and monitoring activities.

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.

A First Assessment of Canopy Cover Loss in Germany’s Forests after the 2018–2020 Drought Years

Central Europe was hit by several unusually strong periods of drought and heat between 2018 and 2020. These droughts affected forest ecosystems. Cascading effects with bark beetle infestations in spruce stands were fatal to vast forest areas in Germany. We present the first assessment of canopy cover loss in Germany for the period of January 2018–April 2021. Our approach makes use of dense Sentinel-2 and Landsat-8 time-series data. We computed the disturbance index (DI) from the tasseled cap components brightness, greenness, and wetness. Using quantiles, we generated monthly DI composites and calculated anomalies in a reference period (2017). From the resulting map, we calculated the canopy cover loss statistics for administrative entities. Our results show a canopy cover loss of 501,000 ha for Germany, with large regional differences. The losses were largest in central Germany and reached up to two-thirds of coniferous forest loss in some districts. Our map has high spatial (10 m) and temporal (monthly) resolution and can be updated at any time.

Land Cover and Vegetation Coverage Changes in the Mining Area—A Case Study from Slovakia

Dealing with landscape changes in space and time is an important activity in terms of the process of future development of the selected area. In particular, it is necessary to focus on territories that are exposed to the effects of extraction activities. The main objective of the paper was the mapping of spatio-temporal changes in the landscape in connection with the extraction of minerals due to mining activities on the landscape using satellite images and data from the Corine land cover (CLC) database in the environment of geographic information systems. The selected study area is specific to the presence of four mineral deposits (three of which are under active mining). The Rohožník-Konopiská deposit was abandoned and the area was subsequently reclaimed. The study used Corine land cover (CLC) data and Landsat 5, 7, 8 satellite images for selected years in the period 1990–2021. The Normalized Difference Vegetation Index (NDVI) was calculated for vegetation cover analysis, which was further combined with the forest spatial division units (FSDU) layer. Areas in the immediate vicinity of the open-pit mine were selected for detailed analysis of vegetation changes. Using the FSDU data, an average NDVI index value was calculated using the Zonal statistics function for each plot. The results showed that over the selected period there have been changes indicating an improvement in the landscape condition by reclamation operations at two deposits, Rohožník-Konopiská (inactive) and Sološnica-Hrabník (active). The analyzed CLC data detected the change at the Rohožník-Konopiská deposit, but the active deposit Sološnica-Hrabník was not detected in these data. The loss of vegetation on the other two deposits is mainly due to pre-mining preparatory work, which causes the removal of soil and vegetation layers.

Impact of High-Cadence Earth Observation in Maize Crop Phenology Classification

For farmers, policymakers, and government agencies, it is critical to accurately define agricultural crop phenology and its spatial-temporal variability. At the moment, two approaches are utilized to report crop phenology. On one hand, land surface phenology provides information about the overall trend, whereas weekly reports from USDA-NASS provide information about the development of particular crops at the regional level. High-cadence earth observations might help to improve the accuracy of these estimations and bring more precise crop phenology classifications closer to what farmers demand. The second component of the proposed solution requires the use of robust classifiers (e.g., random forest, RF) capable of successfully managing large data sets. To evaluate this solution, this study compared the output of a RF classifier model using weather, two different satellite sources (Planet Fusion; PF and Sentinel-2; S-2), and ground truth data to improve maize (Zea mays L.) crop phenology classification using two regions of Kansas (Southwest and Central) as a testbed during the 2017 growing season. Our findings suggests that high temporal resolution (PF) data can significantly improve crop classification metrics (f1-score = 0.94) relative to S-2 (f1-score = 0.86). Additionally, a decline in the f1-score between 0.74 and 0.60 was obtained when we assessed the ability of S-2 to extend the temporal forecast for crop phenology. This research highlights the critical nature of very high temporal resolution (daily) earth observation data for crop monitoring and decision making in agriculture.

The Swiss data cube, analysis ready data archive using earth observations of Switzerland

Since the opening of Earth Observation (EO) archives (USGS/NASA Landsat and EC/ESA Sentinels), large collections of EO data are freely available, offering scientists new possibilities to better understand and quantify environmental changes. Fully exploiting these satellite EO data will require new approaches for their acquisition, management, distribution, and analysis. Given rapid environmental changes and the emergence of big data, innovative solutions are needed to support policy frameworks and related actions toward sustainable development. Here we present the Swiss Data Cube (SDC), unleashing the information power of Big Earth Data for monitoring the environment, providing Analysis Ready Data over the geographic extent of Switzerland since 1984, which is updated on a daily basis. Based on a cloud-computing platform allowing to access, visualize and analyse optical (Sentinel-2; Landsat 5, 7, 8) and radar (Sentinel-1) imagery, the SDC minimizes the time and knowledge required for environmental analyses, by offering consistent calibrated and spatially co-registered satellite observations. SDC derived analysis ready data supports generation of environmental information, allowing to inform a variety of environmental policies with unprecedented timeliness and quality.

Rewetting does not return drained fen peatlands to their old selves

Peatlands have been drained for land use for a long time and on a large scale, turning them from carbon and nutrient sinks into respective sources, diminishing water regulation capacity, causing surface height loss and destroying biodiversity. Over the last decades, drained peatlands have been rewetted for biodiversity restoration and, as it strongly decreases greenhouse gas emissions, also for climate protection. We quantify restoration success by comparing 320 rewetted fen peatland sites to 243 near-natural peatland sites of similar origin across temperate Europe, all set into perspective by 10k additional European fen vegetation plots. Results imply that rewetting of drained fen peatlands induces the establishment of tall, graminoid wetland plants (helophytisation) and long-lasting differences to pre-drainage biodiversity (vegetation), ecosystem functioning (geochemistry, hydrology), and land cover characteristics (spectral temporal metrics). The Paris Agreement entails the rewetting of 500,000 km² of drained peatlands worldwide until 2050-2070. A better understanding of the resulting locally novel ecosystems is required to improve planning and implementation of peatland rewetting and subsequent management.

Whether rewetting leads to effective restoration of drained peatlands is unclear. Here the authors analyse a large number of near-natural and rewetted fen peatland sites in Europe, finding persistent differences in plant community composition and ecosystem functioning, and higher variance in the restored sites.

Mapping Crop Types and Cropping Systems in Nigeria with Sentinel-2 Imagery

Reliable crop type maps from satellite data are an essential prerequisite for quantifying crop growth, health, and yields. However, such maps do not exist for most parts of Africa, where smallholder farming is the dominant system. Prevalent cloud cover, small farm sizes, and mixed cropping systems pose substantial challenges when creating crop type maps for sub-Saharan Africa. In this study, we provide a mapping scheme based on freely available Sentinel-2A/B (S2) time series and very high-resolution SkySat data to map the main crops—maize and potato—and intercropping systems including these two crops on the Jos Plateau, Nigeria. We analyzed the spectral-temporal behavior of mixed crop classes to improve our understanding of inter-class spectral mixing. Building on the Framework for Operational Radiometric Correction for Environmental monitoring (FORCE), we preprocessed S2 time series and derived spectral-temporal metrics from S2 spectral bands for the main temporal cropping windows. These STMs were used as input features in a hierarchical random forest classification. Our results provide the first wall-to-wall crop type map for this key agricultural region of Nigeria. Our cropland identification had an overall accuracy of 84%, while the crop type map achieved an average accuracy of 72% for the five relevant crop classes. Our crop type map shows distinctive regional variations in the distribution of crop types. Maize is the dominant crop, followed by mixed cropping systems, including maize–cereals and potato–maize cropping; potato was found to be the least prevalent class. Plot analyses based on a sample of 1166 fields revealed largely homogeneous mapping patterns, demonstrating the effectiveness of our classification system also for intercropped classes, which are temporally and spatially highly heterogeneous. Moreover, we found that small field sizes were dominant in all crop types, regardless of whether or not intercropping was used. Maize–legume and maize exhibited the largest plots, with an area of up to 3 ha and slightly more than 10 ha, respectively; potato was mainly cultivated on fields smaller than 0.5 ha and only a few plots were larger than 1 ha. Besides providing the first spatially explicit map of cropping practices in the core production area of the Jos Plateau, Nigeria, the study also offers guidance for the creation of crop type maps for smallholder-dominated systems with intercropping. Critical temporal windows for crop type differentiation will enable the creation of mapping approaches in support of future smart agricultural practices for aspects such as food security, early warning systems, policies, and extension services.

European Space Agency (ESA) Calibration/Validation Strategy for Optical Land-Imaging Satellites and Pathway towards Interoperability

Land remote sensing capabilities in the optical domain have dramatically increased in the past decade, owing to the unprecedented growth of space-borne systems providing a wealth of measurements at enhanced spatial, temporal and spectral resolutions. Yet, critical questions remain as how to unlock the potential of such massive amounts of data, which are complementary in principle but inherently diverse in terms of products specifications, algorithm definition and validation approaches. Likewise, there is a recent increase in spatiotemporal coverage of in situ reference data, although inconsistencies in the used measurement practices and in the associated quality information still hinder their integrated use for satellite products validation. In order to address the above-mentioned challenges, the European Space Agency (ESA), in collaboration with other Space Agencies and international partners, is elaborating a strategy for establishing guidelines and common protocols for the calibration and validation (Cal/Val) of optical land imaging sensors. Within this paper, this strategy will be illustrated and put into the context of current validation systems for land remote sensing. A reinforced focus on metrology is the basic principle underlying such a strategy, since metrology provides the terminology, the framework and the best practices, allowing to tie measurements acquired from a variety of sensors to internationally agreed upon standards. From this general concept, a set of requirements are derived on how the measurements should be acquired, analysed and quality reported to users using unified procedures. This includes the need for traceability, a fully characterised uncertainty budget and adherence to community-agreed measurement protocols. These requirements have led to the development of the Fiducial Reference Measurements (FRM) concept, which is promoted by the ESA as the recommended standard within the satellite validation community. The overarching goal is to enhance user confidence in satellite-based data and characterise inter-sensor inconsistencies, starting from at-sensor radiances and paving the way to achieving the interoperability of current and future land-imaging systems.

A Sentinel-2 Dataset for Uganda

Earth observation data provide useful information for the monitoring and management of vegetation- and land-related resources. The Framework for Operational Radiometric Correction for Environmental monitoring (FORCE) was used to download, process and composite Sentinel-2 data from 2018–2020 for Uganda. Over 16,500 Sentinel-2 data granules were downloaded and processed from top of the atmosphere reflectance to bottom of the atmosphere reflectance and higher-level products, totalling > 9 TB of input data. The output data include the number of clear sky observations per year, the best available pixel composite per year and vegetation indices (mean of EVI and NDVI) per quarter. The study intention was to provide analysis-ready data for all of Uganda from Sentinel-2 at 10 m spatial resolution, allowing users to bypass some basic processing and, hence, facilitate environmental monitoring.

High-Resolution Maps of Material Stocks in Buildings and Infrastructures in Austria and Germany

The dynamics of societal material stocks such as buildings and infrastructures and their spatial patterns drive surging resource use and emissions. Two main types of data are currently used to map stocks, night-time lights (NTL) from Earth-observing (EO) satellites and cadastral information. We present an alternative approach for broad-scale material stock mapping based on freely available high-resolution EO imagery and OpenStreetMap data. Maps of built-up surface area, building height, and building types were derived from optical Sentinel-2 and radar Sentinel-1 satellite data to map patterns of material stocks for Austria and Germany. Using material intensity factors, we calculated the mass of different types of buildings and infrastructures, distinguishing eight types of materials, at 10 m spatial resolution. The total mass of buildings and infrastructures in 2018 amounted to ∼5 Gt in Austria and ∼38 Gt in Germany (AT: ∼540 t/cap, DE: ∼450 t/cap). Cross-checks with independent data sources at various scales suggested that the method may yield more complete results than other data sources but could not rule out possible overestimations. The method yields thematic differentiations not possible with NTL, avoids the use of costly cadastral data, and is suitable for mapping larger areas and tracing trends over time.

Comparison of Masking Algorithms for Sentinel-2 Imagery

Masking of clouds, cloud shadow, water and snow/ice in optical satellite imagery is an important step in automated processing chains. We compare the performance of the masking provided by Fmask (“Function of mask” implemented in FORCE), ATCOR (“Atmospheric Correction”) and Sen2Cor (“Sentinel-2 Correction”) on a set of 20 Sentinel-2 scenes distributed over the globe covering a wide variety of environments and climates. All three methods use rules based on physical properties (Top of Atmosphere Reflectance, TOA) to separate clear pixels from potential cloud pixels, but they use different rules and class-specific thresholds. The methods can yield different results because of different definitions of the dilation buffer size for the classes cloud, cloud shadow and snow. Classification results are compared to the assessment of an expert human interpreter using at least 50 polygons per class randomly selected for each image. The class assignment of the human interpreter is considered as reference or “truth”. The interpreter carefully assigned a class label based on the visual assessment of the true color and infrared false color images and additionally on the bottom of atmosphere (BOA) reflectance spectra. The most important part of the comparison is done for the difference area of the three classifications considered. This is the part of the classification images where the results of Fmask, ATCOR and Sen2Cor disagree. Results on difference area have the advantage to show more clearly the strengths and weaknesses of a classification than results on the complete image. The overall accuracy of Fmask, ATCOR, and Sen2Cor for difference areas of the selected scenes is 45%, 56%, and 62%, respectively. User and producer accuracies are strongly class- and scene-dependent, typically varying between 30% and 90%. Comparison of the difference area is complemented by looking for the results in the area where all three classifications give the same result. Overall accuracy for that “same area” is 97% resulting in the complete classification in overall accuracy of 89%, 91% and 92% for Fmask, ATCOR and Sen2Cor respectively.

National-scale mapping of building height using Sentinel-1 and Sentinel-2 time series

Urban areas and their vertical characteristics have a manifold and far-reaching impact on our environment. However, openly accessible information at high spatial resolution is still missing at large for complete countries or regions. In this study, we combined Sentinel-1A/B and Sentinel-2A/B time series to map building heights for entire Germany on a 10 m grid resolving built-up structures in rural and urban contexts. We utilized information from the spectral/polarization, temporal and spatial dimensions by combining band-wise temporal aggregation statistics with morphological metrics. We trained machine learning regression models with highly accurate building height information from several 3D building models. The novelty of this method lies in the very fine resolution yet large spatial extent to which it can be applied, as well as in the use of building shadows in optical imagery. Results indicate that both radar-only and optical-only models can be used to predict building height, but the synergistic combination of both data sources leads to superior results. When testing the model against independent datasets, very consistent performance was achieved (frequency-weighted RMSE of 2.9 m to 3.5 m), which suggests that the prediction of the most frequently occurring buildings was robust. The average building height varies considerably across Germany with lower buildings in Eastern and South-Eastern Germany and taller ones along the highly urbanized areas in Western Germany. We emphasize the straightforward applicability of this approach on the national scale. It mostly relies on freely available satellite imagery and open source software, which potentially permit frequent update cycles and cost-effective mapping that may be relevant for a plethora of different applications, e.g. physical analysis of structural features or mapping society's resource usage.

Sentinel-1 and Sentinel-2 Data for Savannah Land Cover Mapping: Optimising the Combination of Sensors and Seasons

Savannahs are heterogeneous environments with an important role in supporting biodiversity and providing essential ecosystem services. Due to extensive land use/cover changes and subsequent land degradation, the provision of ecosystems services from savannahs has increasingly declined over recent years. Mapping the extent and the composition of savannah environments is challenging but essential in order to improve monitoring capabilities, prevent biodiversity loss and ensure the provision of ecosystem services. Here, we tested combinations of Sentinel-1 and Sentinel-2 data from three different seasons to optimise land cover mapping, focusing in the Ngorongoro Conservation Area (NCA) in Tanzania. The NCA has a bimodal rainfall pattern and is composed of a combination savannah and woodland landscapes. The best performing model achieved an overall accuracy of 86.3 ± 1.5% and included a combination of Sentinel-1 and 2 from the dry and short-dry seasons. Our results show that the optical models outperform their radar counterparts, the combination of multisensor data improves the overall accuracy in all scenarios and this is particularly advantageous in single-season models. Regarding the effect of season, models that included the short-dry season outperform the dry and wet season models, as this season is able to provide cloud free data and is wet enough to allow for the distinction between woody and herbaceous vegetation. Additionally, the combination of more than one season is beneficial for the classification, specifically if it includes the dry or the short-dry season. Combining several seasons is, overall, more beneficial for single-sensor data; however, the accuracies varied with land cover. In summary, the combination of several seasons and sensors provides a more accurate classification, but the target vegetation types should be taken into consideration.

Automated Mosaicking of Sentinel-2 Satellite Imagery

Repeat frequencies of optical remote sensing satellites have been increasing over the last 40 years, but there is still dependence on clear skies to acquire usable imagery. To increase the quality of data, composited mosaics of satellite imagery can be used. In this paper, we develop an automated method for clearing clouds and producing different types of composited mosaics suitable for use in cloud-affected countries, such as New Zealand. We improve the Tmask algorithm for cloud detection by using a parallax method to produce an initial cloud layer and by using an object-based cloud and shadow approach to remove false cloud detections. We develop several parametric scoring approaches for choosing best-pixel composites with minimal remaining cloud. The automated mosaicking approach produced Sentinel-2 mosaics of New Zealand for five successive summers, 2015/16 through 2019/20, with remaining cloud being less than 0.1%. Contributing satellite overpasses were typically of the order of 100. In comparison, manual methods for cloud clearing produced mosaics with 5% remaining cloud and from satellite overpasses typically of the order of 20. The improvements to cloud clearing enable the use of all possible Sentinel-2 imagery to produce automatic mosaics capable of regular land monitoring, at a reasonable cost.

Land Cover Dynamics and Mangrove Degradation in the Niger Delta Region

The Niger Delta Region is the largest river delta in Africa and features the fifth largest mangrove forest on Earth. It provides numerous ecosystem services to the local populations and holds a wealth of biodiversity. However, due to the oil and gas reserves and the explosion of human population it is under threat from overexploitation and degradation. There is a pressing need for an accurate assessment of the land cover dynamics in the region. The limited previous efforts have produced controversial results, as the area of western Africa is notorious for the gaps in the Landsat archive and the lack of cloud-free data. Even fewer studies have attempted to map the extent of the degraded mangrove forest system, reporting low accuracies. Here, we map the eight main land cover classes over the NDR using spectral-temporal metrics from all available Landsat data centred around three epochs. We also test the performance of the classification when L-band radar data are added to the Landsat-based metrics. To further our understanding of the land cover change dynamics, we carry out two additional assessments: a change intensity analysis for the entire NDR and, focusing specifically on the mangrove forest, we analyse the fragmentation of both the healthy and the degraded mangrove land cover classes. We achieve high overall classification accuracies in all epochs (~79% for 1988, and 82% for 2000 and 2013) and are able to map the degraded mangroves accurately, for the first time, with user’s accuracies between 77% and 87% and producer’s accuracies consistently above 82%. Our results show that mangrove forests, lowland rainforests, and freshwater forests are reporting net and highly intense losses (mangrove net loss: ~500 km2; woodland net loss: ~1400 km2), while built-up areas have almost doubled in size (from 1990 km2 in 1988 to 3730 km2 in 2013). The mangrove forests are also consistently more fragmented, with the opposite effect being observed for the degraded mangroves in more recent years. Our study provides a valuable assessment of land cover dynamics in the NDR and the first ever accurate estimates of the extent of the degraded mangrove forest and its fragmentation.

Sensitivity analysis of the dark spectrum fitting atmospheric correction for metre- and decametre-scale satellite imagery using autonomous hyperspectral radiometry

The performance of the dark spectrum fitting (DSF) atmospheric correction algorithm is evaluated using matchups between metre- and decametre-scale satellite imagery as processed with ACOLITE and measurements from autonomous PANTHYR hyperspectral radiometer systems deployed in the Adriatic and North Sea. Imagery from the operational land imager (OLI) on Landsat 8, the multispectral instrument (MSI) on Sentinel-2 A and B, and the PlanetScope CubeSat constellation was processed for both sites using a fixed atmospheric path reflectance in a small region of interest around the system's deployment location, using a number of processing settings, including a new sky reflectance correction. The mean absolute relative differences (MARD) between in situ and satellite measured reflectances reach <20% in the Blue and 11% in the Green bands around 490 and 560 nm for the best performing configuration for MSI and OLI. Higher relative errors are found for the shortest Blue bands around 440 nm (30-100% MARD), and in the Red-Edge and near-infrared bands (35-100% MARD), largely influenced by the lower absolute data range in the observations. Root mean squared differences (RMSD) increase from 0.005 in the NIR to about 0.015-0.020 in the Blue band, consistent with increasing atmospheric path reflectance. Validation of the Red-Edge and NIR bands on Sentinel-2 is presented, as well as for the first time, the Panchromatic band (17-26% MARD) on Landsat 8, and the derived Orange contra-band (8-33% MARD for waters in the algorithm domain, and around 40-80% MARD overall). For Sentinel-2, excluding the SWIR bands from the DSF gave better performances, likely due to calibration issues of MSI at longer wavelengths. Excluding the SWIR on Landsat 8 gave good performance as well, indicating robustness of the DSF to the available band set. The DSF performance was found to be rather insensitive to (1) the wavelength spacing in the lookup tables used for the atmospheric correction, (2) the use of default or ancillary information on gas concentration and atmospheric pressure, and (3) the size of the ROI over which the path reflectance is estimated. The performance of the PlanetScope constellation is found to be similar to previously published results, with the standard DSF giving the best results in the visible bands in terms of MARD (24-40% overall, and 18-29% for the turbid site). The new sky reflectance correction gave mixed results, although it reduced the mean biases for certain configurations and improved results for the processing excluding the SWIR bands, giving lower RMSD and MARD especially at longer wavelengths (>600 nm). The results presented in this article should serve as guidelines for general use of ACOLITE and the DSF.

Recent Applications of Landsat 8/OLI and Sentinel-2/MSI for Land Use and Land Cover Mapping: A Systematic Review

Recent applications of Landsat 8 Operational Land Imager (L8/OLI) and Sentinel-2 MultiSpectral Instrument (S2/MSI) data for acquiring information about land use and land cover (LULC) provide a new perspective in remote sensing data analysis. Jointly, these sources permit researchers to improve operational classification and change detection, guiding better reasoning about landscape and intrinsic processes, as deforestation and agricultural expansion. However, the results of their applications have not yet been synthesized in order to provide coherent guidance on the effect of their applications in different classification processes, as well as to identify promising approaches and issues which affect classification performance. In this systematic review, we present trends, potentialities, challenges, actual gaps, and future possibilities for the use of L8/OLI and S2/MSI for LULC mapping and change detection. In particular, we highlight the possibility of using medium-resolution (Landsat-like, 10–30 m) time series and multispectral optical data provided by the harmonization between these sensors and data cube architectures for analysis-ready data that are permeated by publicizations, open data policies, and open science principles. We also reinforce the potential for exploring more spectral bands combinations, especially by using the three Red-edge and the two Near Infrared and Shortwave Infrared bands of S2/MSI, to calculate vegetation indices more sensitive to phenological variations that were less frequently applied for a long time, but have turned on since the S2/MSI mission. Summarizing peer-reviewed papers can guide the scientific community to the use of L8/OLI and S2/MSI data, which enable detailed knowledge on LULC mapping and change detection in different landscapes, especially in agricultural and natural vegetation scenarios.

Multispectral high resolution sensor fusion for smoothing and gap-filling in the cloud

Remote sensing optical sensors onboard operational satellites cannot have high spectral, spatial and temporal resolutions simultaneously. In addition, clouds and aerosols can adversely affect the signal contaminating the land surface observations. We present a HIghly Scalable Temporal Adaptive Reflectance Fusion Model (HISTARFM) algorithm to combine multispectral images of different sensors to reduce noise and produce monthly gap free high resolution (30 m) observations over land. Our approach uses images from the Landsat (30 m spatial resolution and 16 day revisit cycle) and the MODIS missions, both from Terra and Aqua platforms (500 m spatial resolution and daily revisit cycle). We implement a bias-aware Kalman filter method in the Google Earth Engine (GEE) platform to obtain fused images at the Landsat spatial-resolution. The added bias correction in the Kalman filter estimates accounts for the fact that both model and observation errors are temporally auto-correlated and may have a non-zero mean. This approach also enables reliable estimation of the uncertainty associated with the final reflectance estimates, allowing for error propagation analyses in higher level remote sensing products. Quantitative and qualitative evaluations of the generated products through comparison with other state-of-the-art methods confirm the validity of the approach, and open the door to operational applications at enhanced spatio-temporal resolutions at broad continental scales.

Mapping urban-rural gradients of settlements and vegetation at national scale using Sentinel-2 spectral-temporal metrics and regression-based unmixing with synthetic training data

The increasing impact of humans on land and ongoing global population growth requires an improved understanding of land cover (LC) and land use (LU) processes related to settlements. The heterogeneity of built-up areas and infrastructures as well as the importance of not only mapping, but also characterizing anthropogenic structures suggests using a sub-pixel mapping approach for analysing related LC from space. We implement a regression-based unmixing approach for mapping built-up surfaces and infrastructure, woody vegetation and non-woody vegetation for all of Germany and Austria at 10 m resolution to demonstrate the potential of sub-pixel mapping. We map LC fractions for one point in time, using all available Sentinel-2 data from 2017 and 2018 (<70% cloud cover). We combine the concept of synthetically mixed training data with statistical aggregations from spectral-temporal metrics (STM) derived from Sentinel-2 reflectance time series. We specifically examine how STM can be used for creating synthetically mixed training data. STM are known to facilitate large area mapping by being largely independent of image acquisition dates and inherently incorporate phenological information. Vegetation is an important part of settlements and time series information supports its mapping. Synthetically mixed training data facilitates a streamlined training by using pure reference spectra to generate artificial mixtures as input to regression modelling of LC fractions in mixed pixels. We here show how combining both offers great potential for wall-to-wall LC fraction mapping. We further investigate the positive effect of STM on map results by comparing the performance of different subsets of STM combinations. Our results indicate that many STM combinations containing spectral variability and vegetation indices provide suitable input to creating synthetic training data for regression-based fraction mapping. Results for built-up surfaces and infrastructure (MAE 0.13/RMSE 0.18 at 20 m resolution), woody vegetation (0.18, 0.22) and non-woody vegetation (0.14, 0.19) are highly consistent across Germany and Austria. Only a few surface types were not accurately predicted in our nation-wide mapping. Further research is required to optimize mapping of temporally invariant bare soil and rock surfaces that show spectral similarity to built-up surfaces and infrastructure. The proposed methodology combines benefits of both regression-based modelling with synthetically mixed training data and STM, and thus facilitates mapping of LC fractions on a national scale and at high resolution. Such information will allow to better characterize settlements and identifying processes such as densification that are best represented by continuous LC mapping.

Mosaicking Landsat and Sentinel-2 Data to Enhance LandTrendr Time Series Analysis in Northern High Latitude Permafrost Regions

Permafrost is warming in the northern high latitudes, inducing highly dynamic thaw-related permafrost disturbances across the terrestrial Arctic. Monitoring and tracking of permafrost disturbances is important as they impact surrounding landscapes, ecosystems and infrastructure. Remote sensing provides the means to detect, map, and quantify these changes homogeneously across large regions and time scales. Existing Landsat-based algorithms assess different types of disturbances with similar spatiotemporal requirements. However, Landsat-based analyses are restricted in northern high latitudes due to the long repeat interval and frequent clouds, in particular at Arctic coastal sites. We therefore propose to combine Landsat and Sentinel-2 data for enhanced data coverage and present a combined annual mosaic workflow, expanding currently available algorithms, such as LandTrendr, to achieve more reliable time series analysis. We exemplary test the workflow for twelve sites across the northern high latitudes in Siberia. We assessed the number of images and cloud-free pixels, the spatial mosaic coverage and the mosaic quality with spectral comparisons. The number of available images increased steadily from 1999 to 2019 but especially from 2016 onward with the addition of Sentinel-2 images. Consequently, we have an increased number of cloud-free pixels even under challenging environmental conditions, which then serve as the input to the mosaicking process. In a comparison of annual mosaics, the Landsat+Sentinel-2 mosaics always fully covered the study areas (99.9–100 %), while Landsat-only mosaics contained data-gaps in the same years, only reaching coverage percentages of 27.2 %, 58.1 %, and 69.7 % for Sobo Sise, East Taymyr, and Kurungnakh in 2017, respectively. The spectral comparison of Landsat image, Sentinel-2 image, and Landsat+Sentinel-2 mosaic showed high correlation between the input images and mosaic bands (e.g., for Kurungnakh 0.91–0.97 between Landsat and Landsat+Sentinel-2 mosaic and 0.92–0.98 between Sentinel-2 and Landsat+Sentinel-2 mosaic) across all twelve study sites, testifying good quality mosaic results. Our results show that especially the results for northern, coastal areas was substantially improved with the Landsat+Sentinel-2 mosaics. By combining Landsat and Sentinel-2 data we accomplished to create reliably high spatial resolution input mosaics for time series analyses. Our approach allows to apply a high temporal continuous time series analysis to northern high latitude permafrost regions for the first time, overcoming substantial data gaps, and assess permafrost disturbance dynamics on an annual scale across large regions with algorithms such as LandTrendr by deriving the location, timing and progression of permafrost thaw disturbances.

Sharpening the Sentinel-2 10 and 20 m Bands to Planetscope-0 3 m Resolution

Combination of near daily 3 m red, green, blue, and near infrared (NIR) Planetscope reflectance with lower temporal resolution 10 m and 20 m red, green, blue, NIR, red-edge, and shortwave infrared (SWIR) Sentinel-2 reflectance provides potential for improved global monitoring. Sharpening the Sentinel-2 reflectance with the Planetscope reflectance may enable near-daily 3 m monitoring in the visible, red-edge, NIR, and SWIR. However, there are two major issues, namely the different and spectrally nonoverlapping bands between the two sensors and surface changes that may occur in the period between the different sensor acquisitions. They are examined in this study that considers Sentinel-2 and Planetscope imagery acquired one day apart over three sites where land surface changes due to biomass burning occurred. Two well-established sharpening methods, high pass modulation (HPM) and Model 3 (M3), were used as they are multiresolution analysis methods that preserve the spectral properties of the low spatial resolution Sentinel-2 imagery (that are better radiometrically calibrated than Planetscope) and are relatively computationally efficient so that they can be applied at large scale. The Sentinel-2 point spread function (PSF) needed for the sharpening was derived analytically from published modulation transfer function (MTF) values. Synthetic Planetscope red-edge and SWIR bands were derived by linear regression of the Planetscope visible and NIR bands with the Sentinel-2 red-edge and SWIR bands. The HPM and M3 sharpening results were evaluated visually and quantitatively using the Q2n metric that quantifies spectral and spatial distortion. The HPM and M3 sharpening methods provided visually coherent and spatially detailed visible and NIR wavelength sharpened results with low distortion (Q2n values > 0.91). The sharpened red-edge and SWIR results were also coherent but had greater distortion (Q2n values > 0.76). Detailed examination at locations where surface changes between the Sentinel-2 and the Planetscope acquisitions occurred revealed that the HPM method, unlike the M3 method, could reliably sharpen the bands affected by the change. This is because HPM sharpening uses a per-pixel reflectance ratio in the spatial detail modulation which is relatively stable to reflectance changes. The paper concludes with a discussion of the implications of this research and the recommendation that the HPM sharpening be used considering its better performance when there are surface changes.

Monitoring Vegetation Change in the Presence of High Cloud Cover with Sentinel-2 in a Lowland Tropical Forest Region in Brazil

Forests play major roles in climate regulation, ecosystem services, carbon storage, biodiversity, terrain stabilization, and water retention, as well as in the economy of numerous countries. Nevertheless, deforestation and forest degradation are rampant in many parts of the world. In particular, the Amazonian rainforest faces the constant threats posed by logging, mining, and burning for agricultural expansion. In Brazil, the “Sete de Setembro Indigenous Land”, a protected area located in a lowland tropical forest region at the border between the Mato Grosso and Rondônia states, is subject to illegal deforestation and therefore necessitates effective vegetation monitoring tools. Optical satellite imagery, while extensively used for landcover assessment and monitoring, is vulnerable to high cloud cover percentages, as these can preclude analysis and strongly limit the temporal resolution. We propose a cloud computing-based coupled detection strategy using (i) cloud and cloud shadow/vegetation detection systems with Sentinel-2 data analyzed on the Google Earth Engine with deep neural network classification models, with (ii) a classification error correction and vegetation loss and gain analysis tool that dynamically compares and updates the classification in a time series. The initial results demonstrate that such a detection system can constitute a powerful monitoring tool to assist in the prevention, early warning, and assessment of deforestation and forest degradation in cloudy tropical regions. Owing to the integrated cloud detection system, the temporal resolution is significantly improved. The limitations of the model in its present state include classification issues during the forest fire period, and a lack of distinction between natural vegetation loss and anthropogenic deforestation. Two possible solutions to the latter problem are proposed, namely, the mapping of known agricultural and bare areas and its subsequent removal from the analyzed data, or the inclusion of radar data, which would allow a large amount of finetuning of the detection processes.

DATimeS: A machine learning time series GUI toolbox for gap-filling and vegetation phenology trends detection

Optical remotely sensed data are typically discontinuous, with missing values due to cloud cover. Consequently, gap-filling solutions are needed for accurate crop phenology characterization. The here presented Decomposition and Analysis of Time Series software (DATimeS) expands established time series interpolation methods with a diversity of advanced machine learning fitting algorithms (e.g., Gaussian Process Regression: GPR) particularly effective for the reconstruction of multiple-seasons vegetation temporal patterns. DATimeS is freely available as a powerful image time series software that generates cloud-free composite maps and captures seasonal vegetation dynamics from regular or irregular satellite time series. This work describes the main features of DATimeS, and provides a demonstration case using Sentinel-2 Leaf Area Index time series data over a Spanish site. GPR resulted as an optimum fitting algorithm with most accurate gap-filling performance and associated uncertainties. DATimeS further quantified LAI fluctuations among multiple crop seasons and provided phenological indicators for specific crop types.

Long-Term Land Use/Land Cover Change Assessment of the Kilombero Catchment in Tanzania Using Random Forest Classification and Robust Change Vector Analysis

Information about land use/land cover (LULC) and their changes is useful for different stakeholders to assess future pathways of sustainable land use for food production as well as for nature conservation. In this study, we assess LULC changes in the Kilombero catchment in Tanzania, an important area of recent development in East Africa. LULC change is assessed in two ways: first, post-classification comparison (PCC) which allows us to directly assess changes from one LULC class to another, and second, spectral change detection. We perform LULC classification by applying random forests (RF) on sets of multitemporal metrics that account for seasonal within-class dynamics. For the spectral change detection, we make use of the robust change vector analysis (RCVA) and determine those changes that do not necessarily lead to another class. The combination of the two approaches enables us to distinguish areas that show (a) only PCC changes, (b) only spectral changes that do not affect the classification of a pixel, (c) both types of change, or (d) no changes at all. Our results reveal that only one-quarter of the catchment has not experienced any change. One-third shows both, spectral changes and LULC conversion. Changes detected with both methods predominantly occur in two major regions, one in the West of the catchment, one in the Kilombero floodplain. Both regions are important areas of food production and economic development in Tanzania. The Kilombero floodplain is a Ramsar protected area, half of which was converted to agricultural land in the past decades. Therefore, LULC monitoring is required to support sustainable land management. Relatively poor classification performances revealed several challenges during the classification process. The combined approach of PCC and RCVA allows us to detect spatial patterns of LULC change at distinct dimensions and intensities. With the assessment of additional classifier output, namely class-specific per-pixel classification probabilities and derived parameters, we account for classification uncertainty across space. We overlay the LULC change results and the spatial assessment of classification reliability to provide a thorough picture of the LULC changes taking place in the Kilombero catchment.

Landsat Analysis Ready Data for Global Land Cover and Land Cover Change Mapping

The multi-decadal Landsat data record is a unique tool for global land cover and land use change analysis. However, the large volume of the Landsat image archive and inconsistent coverage of clear-sky observations hamper land cover monitoring at large geographic extent. Here, we present a consistently processed and temporally aggregated Landsat Analysis Ready Data produced by the Global Land Analysis and Discovery team at the University of Maryland (GLAD ARD) suitable for national to global empirical land cover mapping and change detection. The GLAD ARD represent a 16-day time-series of tiled Landsat normalized surface reflectance from 1997 to present, updated annually, and designed for land cover monitoring at global to local scales. A set of tools for multi-temporal data processing and characterization using machine learning provided with GLAD ARD serves as an end-to-end solution for Landsat-based natural resource assessment and monitoring. The GLAD ARD data and tools have been implemented at the national, regional, and global extent for water, forest, and crop mapping. The GLAD ARD data and tools are available at the GLAD website for free access.

Seasonal Crop Water Balance Using Harmonized Landsat-8 and Sentinel-2 Time Series Data

Efficient water management in agriculture requires a precise estimate of evapotranspiration ( E T ). Although local measurements can be used to estimate surface energy balance components, these values cannot be extrapolated to large areas due to the heterogeneity and complexity of agriculture environment. This extrapolation can be done using satellite images that provide information in visible and thermal infrared region of the electromagnetic spectrum; however, most current satellite sensors do not provide this end, but they do include a set of spectral bands that allow the radiometric behavior of vegetation that is highly correlated with the E T . In this context, our working hypothesis states that it is possible to generate a strategy of integration and harmonization of the Normalized Difference Vegetation Index ( N D V I ) obtained from Landsat-8 ( L 8 ) and Sentinel-2 ( S 2 ) sensors in order to obtain an N D V I time series used to estimate E T through fit equations specific to each crop type during an agricultural season (December 2017–March 2018). Based on the obtained results it was concluded that it is possible to estimate E T using an N D V I time series by integrating data from both sensors L 8 and S 2 , which allowed to carry out an updated seasonal water balance over study site, improving the irrigation water management both at plot and water distribution system scale.

Using Landsat and Sentinel-2 Data for the Generation of Continuously Updated Forest Type Information Layers in a Cross-Border Region

From global monitoring to regional forest management there is an increasing demand for information about forest ecosystems. For border regions that are closely connected ecologically and economically, a key factor is the cross-border availability and consistency of up-to-date information such as the forest type. The combination of existing forest information with Earth observation data is a rational method and can provide valuable contribution to serve the increased information demand on a transnational level. We present an approach for the remote sensing-based generation of a transnational and temporally consistent forest type information layer for the German federal states of Rhineland-Palatinate and Saarland, and the Grand Duchy of Luxembourg. Existing forest information data from different countries were merged and combined with suitable vegetation indices derived from Landsat 8 and Sentinel-2 imagery acquired in early spring. An automated bootstrap-based approximation of the optimum threshold for the distinction of “broadleaved” and “coniferous” forest was applied. The spatially explicit forest type information layer is updated annually depending on image availability. Overall accuracies between 79 and 96 percent were obtained. Every spot in the region will be updated successively within a period of expectably three years. The presented approach can be integrated in fully automated processing chains to generate basic forest type information layers on a regular basis.