Sentinel-2 Exposed Soil Composite for Soil Organic Carbon Prediction

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). 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 (R2 2.50 g C kg–1, and RPD 68) for all composites except for the composite in the ‘greening-up’ stage with a NBR2 < 0.07 (R2 = 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.