Since the onset of agriculture, soils have lost their organic carbon to such an extent that the soil functions of many croplands are threatened. Hence, there is a strong demand for mapping and monitoring critical soil properties and in particular soil organic carbon (SOC). Pilot studies have demonstrated the potential for remote sensing techniques for SOC mapping in croplands. It has, however, been shown that the assessment of SOC may be hampered by the condition of the soil surface. While growing vegetation can be readily detected by means of the well-known Normalized Difference Vegetation Index (NDVI), the distinction between bare soil and crop residues is expressed in the shortwave infrared region (SWIR), which is only covered by two broad bands in Landsat or Sentinel-2 imagery. Here we tested the effect of thresholds for the Cellulose Absorption Index (CAI), on the performance of SOC prediction models for cropland soils. Airborne Prism Experiment (APEX) hyperspectral images covering an area of 240 km² in the Belgian Loam Belt were used together with a local soil dataset. We used the partial least square regression (PLSR) model to estimate the SOC content based on 104 georeferenced calibration samples (NDVI < 0.26), firstly without setting a CAI threshold, and obtained a satisfactory result (coefficient of determination (R²) = 0.49, Ratio of Performance to Deviation (RPD) = 1.4 and Root Mean Square Error (RMSE) = 2.13 g kgC⁻¹ for cross-validation). However, a cross comparison of the estimated SOC values to grid-based measurements of SOC content within three fields revealed a systematic overestimation for fields with high residue cover. We then tested different CAI thresholds in order to mask pixels with high residue cover. The best model was obtained for a CAI threshold of 0.75 (R² = 0.59, RPD = 1.5 and RMSE = 1.75 g kgC⁻¹ for cross-validation). These results reveal that the purity of the pixels needs to be assessed aforehand in order to produce reliable SOC maps. The Normalized Burn Ratio (NBR2) index based on the SWIR bands of the MSI Sentinel 2 sensor extracted from images collected nine days before the APEX flight campaign correlates well with the CAI index of the APEX imagery. However, the NBR2 index calculated from Sentinel 2 images under moist conditions is poorly correlated with residue cover. This can be explained by the sensitivity of the NBR2 index to both soil moisture and residues.

Pilot studies have demonstrated the potential of remote sensing for soil organic carbon (SOC) mapping in exposed croplands. However, the use of remote sensing for SOC prediction is often hindered by disturbing factors at the soil surface, such as photosynthetic active and non-photosynthetic active vegetation, variation in soil moisture or surface roughness. With the increasing amount of freely available satellite data, recent studies have focused on stabilizing the soil reflectance by building image composites. These composites tend to minimize the disturbing effects by applying sets of criteria. Here, we aim to develop a robust method that allows selecting Sentinel-2 (S-2) pixels with minimal influence of the following disturbing factors: crop residues, surface roughness and soil moisture. We selected all S-2 cloud-free images covering the Belgian Loam Belt from January 2019 to December 2020 (in total 36 images). We then built nine exposed soil composites based on four sets of criteria: (1) lowest Normalized Burn Ratio (NBR2), (2) Normalized Difference Vegetation Index (NDVI) < 0.25, (3–5) NDVI < 0.25 and NBR2 < threshold, (6) the ‘greening-up’ period of a crop and (7–9) the ‘greening-up’ period of a crop and NBR2 < threshold. The ‘greening-up’ period was selected based on the NDVI timeline, where ‘greening-up’ is considered as the last date of acquisition where the soil is exposed (NDVI < 0.25) before the crop develops (NDVI > 0.25). We then built a partial least square regression (PLSR) model with 10-fold cross-validation to estimate the SOC content based on 137 georeferenced calibration samples on the nine composites. We obtained non-satisfactory results (R² < 0.30, RMSE > 2.50 g C kg^–1, and RPD < 1.4, n > 68) for all composites except for the composite in the ‘greening-up’ stage with a NBR2 < 0.07 (R² = 0.54 ± 0.12, RPD = 1.68 ± 0.45 and RMSE = 2.09 ± 0.39 g C kg^–1, n = 49). Hence, the ‘greening-up’ method combined with a strict NBR2 threshold allows selecting the purest exposed soil pixels suitable for SOC prediction. The limit of this method might be its coverage of the total cropland area, which in a two-year period reached 62%, compared to 95% coverage if only the NDVI threshold is applied.

Since the onset of agriculture, soils have lost their organic carbon to such an extent that the soil functions of many croplands are threatened and there is therefore a strong demand for mapping and monitoring critical soil properties and in particular soil organic carbon (SOC). Pilot studies have demonstrated the potential for remote sensing techniques for SOC mapping in croplands, given their large spatial coverage and high temporal resolution. It has however been shown that the assessment of SOC may be hampered by crop residues. In this study we tested the effect of the threshold for the cellulose absorption index (CAI), on the performance of SOC prediction models for bare cropland soils. Airborne Prism Experiment (APEX) hyperspectral images covering an area of 230 km2 in the Belgian Loam Belt were used together with a local soil dataset. We used the partial least square regression (PLSR) model to estimate the SOC content based on 104 georeferenced calibration samples, firstly without setting a CAI threshold, and obtained a satisfactory result (R²=0.49, RPD=1.4 and RMSE=2.14 g kgC-1 for cross-validation). However, a cross comparison of the estimated SOC values to grid-based measurements of SOC content within three fields revealed a systematic overestimation for fields with high residue cover. We then tested different thresholds of CAI in order to mask pixels with high residue cover, by eliminating calibration samples used in the PLSR model based on this threshold. The best model was obtained for CAI threshold of 0.8 (R²=0.59, RPD=1.5 and RMSE=1.76 g kgC-1 for cross-validation). These results reveal that the purity of the pixels needs to be assessed aforehand in order to produce reliable SOC maps. Preliminary results indicate that an index based on the SWIR bands of the MSI Sentinel 2 sensor is also capable of detecting crop residues. However, the application under moist conditions and for different types of residues needs to be confirmed.

Global efforts to restore the world’s degraded croplands require knowledge on the degree and extent of accelerated soil organic carbon (SOC) loss induced by soil erosion. However, the methods for assessing where and to what extent erosion takes place are still inadequate for precise detection of erosion hotspots at high spatial resolution. Drawing on recent advances in multitemporal Sentinel-2 remote sensing to create a bare soil composite that reflects erosion-induced variations in soil spectral signatures, this study attempted to develop a spectra-based soil erosion mapping approach to pinpoint eroded hotspots in a typical catchment located in the black soil region of northeast China as characterized by undulating landscapes. We built a ground-truth dataset consisting of three classes of soils representing Severe, Moderate and Low erosion intensity because of their inter-class contrasts in estimated erosion rates from ¹³⁷Cs tracing. The spectral separability of different erosion classes was first tested by a combined principal component and linear discriminant analysis (PCA-LDA) against laboratory hyperspectral data and then validated against Sentinel-2-derived broadband spectra. The results show that PCA-LDA produced excellent classification accuracy (Kappa coefficient > 0.9) for both data sources and even more so for Sentinel-2 spectra, highlighting the effectiveness of the multitemporal approach to extract bare soil pixels. Further investigations into the spectral curves enabled identification of distinctive spectral features representative of shifting soil albedo and biochemical composition due to erosion-induced SOC mobilization. A classification scheme comprising the spectral features was applied to the Sentinel-2 bare soil composite for pixel-wise soil erosion mapping, in which 15.9% of the cropland area was detected as erosion hotspots, while the Moderate class occupied 65.4%. Comparing the erosion map to a NDVI map demonstrated the negative impact of soil erosion on crop growth from a spatial perspective, highlighting the potential of the proposed approach to aid targeted cropland management for food security and climate.

Intra-field heterogeneity of soil properties, such as soil organic carbon (SOC), nitrogen (N), phosphorous (P), exchangeable cations, pH, or soil texture, is a function of complex interactions between biological factors, physical factors, and historic agricultural management. Mapping the crop growth and final yield heterogeneity and quantifying their link with soil properties can contribute to an optimization of amendment/fertilizer application and crop yield in a management variable zones (MVZ) approach. To this end, we studied a field of 17 ha consisting of four former fields that were merged in early 2017 and cropped with winter wheat in 2018. Historical management practices data were collected. The topsoil characteristics were analyzed by grid-based sampling and kriged to create maps. We tested the capacity of a multispectral MicaSense^® RedEdge-M^TM camera sensor embedded on an unmanned aerial vehicle (UAV) to map in-season growth of winter wheat. Relating several vegetation indices (VIs) to the plant area index (PAI) measured in the field highlighted the red-edge NDVI (RENDVI) as the most suitable to follow the crop growth throughout the growing season. The georeferenced final grain yield of the winter wheat was measured by a combine harvester. The spatial patterns in RENDVI at three phenological stages were mapped and analyzed together with the yield map. For each of these images a conditional inference forest (CI-forest) algorithm was used to identify the soil properties significantly influencing these spatial patterns. Historical management practices of the four former fields have induced significant heterogeneity in soil properties and crop growth. The spatial patterns of RENDVI are rather constant over time and their Spearman rank correlation with yield is similar along the growing season (r ≃ 0.7). Soil properties explain between 87% (mid-March) to 78% (mid-May) of the variance in RENDVI throughout the growing season, as well as 66% of the variance in yield. The pH and exchangeable K are the most significant factors explaining from 15 to 26% of the variance in crop growth. The methodology proposed in this paper to quantify the importance of soil parameters based on the CI-forest algorithm can contribute to a better management of amendment/fertilizer inputs by stressing the most important parameters to take into consideration for site-specific management. We also showed that heterogeneity induced by the soil properties can be described by a crop map early in the season and that this crop map can be used to optimize soil sampling and thus amendment/fertilizer management.