The Earth’s surface dynamics provide essential information for guiding environmental and agricultural policies. Uncovered and unprotected surfaces experience several undesirable effects, which can affect soil ecosystem functions. We developed a technique to identify global bare surface areas and their dynamics based on multitemporal remote sensing images to aid the spatiotemporal evaluation of anthropic and natural phenomena. The bare Earth’s surface and its changes were recognized by Landsat image processing over a time range of 30 years using the Google Earth Engine platform. Two additional products were obtained with a similar technique: a) Earth’s bare surface frequency, which represents where and how many times a single pixel was detected as bare surface, based on Landsat series, and b) Earth’s bare soil tendency, which represents the tendency of bare surface to increase or decrease. This technique enabled the retrieval of bare surfaces on 32% of Earth’s total land area and on 95% of land when considering only agricultural areas. From a multitemporal perspective, the technique found a 2.8% increase in bare surfaces during the period on a global scale. However, the rate of soil exposure decreased by ~4.8% in the same period. The increase in bare surfaces shows that agricultural areas are increasing worldwide. The decreasing rate of soil exposure indicates that, unlike popular opinion, more soils have been covered due to the adoption of conservation agriculture practices, which may reduce soil degradation.

There is a need to update soil maps and monitor soil organic carbon (SOC) in the upper horizons or plough layer for enabling decision support and land management, while complying with several policies, especially those favoring soil carbon storage. This review paper is dedicated to the satellite-based spectral approaches for SOC assessment that have been achieved from several satellite sensors, study scales and geographical contexts in the past decade. Most approaches relying on pure spectral models have been carried out since 2019 and have dealt with temperate croplands in Europe, China and North America at the scale of small regions, of some hundreds of km²: dry combustion and wet oxidation were the analytical determination methods used for 50% and 35% of the satellite-derived SOC studies, for which measured topsoil SOC contents mainly referred to mineral soils, typically cambisols and luvisols and to a lesser extent, regosols, leptosols, stagnosols and chernozems, with annual cropping systems with a SOC value of ~15 g·kg⁻¹ and a range of 30 g·kg⁻¹ in median. Most satellite-derived SOC spectral prediction models used limited preprocessing and were based on bare soil pixel retrieval after Normalized Difference Vegetation Index (NDVI) thresholding. About one third of these models used partial least squares regression (PLSR), while another third used random forest (RF), and the remaining included machine learning methods such as support vector machine (SVM). We did not find any studies either on deep learning methods or on all-performance evaluations and uncertainty analysis of spatial model predictions. Nevertheless, the literature examined here identifies satellite-based spectral information, especially derived under bare soil conditions, as an interesting approach that deserves further investigations. Future research includes considering the simultaneous analysis of imagery acquired at several dates i.e., temporal mosaicking, testing the influence of possible disturbing factors and mitigating their effects fusing mixed models incorporating non-spectral ancillary information.

Satellite-based soil organic carbon content (SOC) mapping over wide regions is generally hampered by the low soil sampling density and the diversity of soil sampling periods. Some unfavorable topsoil conditions, such as high moisture, rugosity, the presence of crop residues, the limited amplitude of SOC values and the limited area of bare soil when a single image is used, are also among the influencing factors. To generate a reliable SOC map, this study addresses the use of Sentinel-2 (S2) temporal mosaics of bare soil (S2Bsoil) over 6 years jointly with soil moisture products (SMPs) derived from Sentinel 1 and 2 images, SOC measurement data and other environmental covariates derived from digital elevation models, lithology maps and airborne gamma-ray data. In this study, we explore (i) the dates and periods that are preferable to construct temporal mosaics of bare soils while accounting for soil moisture and soil management; (ii) which set of covariates is more relevant to explain the SOC variability. From four sets of covariates, the best contributing set was selected, and the median SOC content along with uncertainty at 90% prediction intervals were mapped at a 25-m resolution from quantile regression forest models. The accuracy of predictions was assessed by 10-fold cross-validation, repeated five times. The models using all the covariates had the best model performance. Airborne gamma-ray thorium, slope and S2 bands (e.g., bands 6, 7, 8, 8a) and indices (e.g., calcareous sedimentary rocks, “calcl”) from the “late winter–spring” time series were the most important covariates in this model. Our results also indicated the important role of neighboring topographic distances and oblique geographic coordinates between remote sensing data and parent material. These data contributed not only to optimizing SOC mapping performance but also provided information related to long-range gradients of SOC spatial variability, which makes sense from a pedological point of view.

Soils are at the crossroads of many existential issues that humanity is currently facing. Soils are a finite resource that is under threat, mainly due to human pressure. There is an urgent need to map and monitor them at field, regional, and global scales in order to improve their management and prevent their degradation. This remains a challenge due to the high and often complex spatial variability inherent to soils. Over the last four decades, major research efforts in the field of pedometrics have led to the development of methods allowing to capture the complex nature of soils. As a result, digital soil mapping (DSM) approaches have been developed for quantifying soils in space and time. DSM and monitoring have become operational thanks to the harmonization of soil databases, advances in spatial modeling and machine learning, and the increasing availability of spatiotemporal covariates, including the exponential increase in freely available remote sensing (RS) data. The latter boosted research in DSM, allowing the mapping of soils at high resolution and assessing the changes through time. We present a review of the main contributions and developments of French (inter)national research, which has a long history in both RS and DSM. Thanks to the French SPOT satellite constellation that started in the early 1980s, the French RS and soil research communities have pioneered DSM using remote sensing. This review describes the data, tools, and methods using RS imagery to support the spatial predictions of a wide range of soil properties and discusses their pros and cons. The review demonstrates that RS data are frequently used in soil mapping (i) by considering them as a substitute for analytical measurements, or (ii) by considering them as covariates related to the controlling factors of soil formation and evolution. It further highlights the great potential of RS imagery to improve DSM, and provides an overview of the main challenges and prospects related to digital soil mapping and future sensors. This opens up broad prospects for the use of RS for DSM and natural resource monitoring.

In agronomy, soil organic carbon (SOC) content is important for the development and growth of crops. From an environmental monitoring viewpoint, SOC sequestration is essential for mitigating the emission of greenhouse gases into the atmosphere. SOC dynamics in cropland soils should be further studied through various approaches including remote sensing. In order to predict SOC content over croplands in southwestern France (area of 22,177 km²), this study addresses (i) the influence of the dates on which Sentinel-2 (S2) images were acquired in the springs of 2017–2018 as well as the influence of the soil sampling period of a set of samples collected between 2005 and 2018, (ii) the use of soil moisture products (SMPs) derived from Sentinel-1/2 satellites to analyze the influence of surface soil moisture on model performance when included as a covariate, and (iii) whether the spatial distribution of SOC as mapped using S2 is related to terrain-derived attributes. The influences of S2 image dates and soil sampling periods were analyzed for bare topsoil. The dates of the S2 images with the best performance (RPD ≥ 1.7) were 6 April and 26 May 2017, using soil samples collected between 2016 and 2018. The soil sampling dates were also analyzed using SMP values. Soil moisture values were extracted for each sample and integrated into partial least squares regression (PLSR) models. The use of soil moisture as a covariate had no effect on the prediction performance of the models; however, SMP values were used to select the driest dates, effectively mapping topsoil organic carbon. S2 was able to predict high SOC contents in the specific soil types located on the old terraces (mesas) shaped by rivers flowing from the southwestern Pyrénées.