Contribution of Sentinel-2 spring seedbed spectra to the digital mapping of soil organic carbon concentration

•Higher model performance for granulometry than SOC% using seedbed spectra.•PLSR outperformed other algorithms for SOC% estimation, CR for granulometry.•Low contribution of seedbed spectra to the DSM of SOC% using environmental covariates. Soil organic carbon (SOC) is central to the functioning of terrestrial ecosystems, has climate mitigation potential and provides several benefits for soil health. Understanding the spatial distribution of SOC can help formulate sustainable soil management practices. Digital soil mapping (DSM) uses advanced statistical and geostatistical methods to estimate soil properties across large areas. DSM integrates climate data, topographic features, geology, legacy soil maps, land management and remote sensing data. Bare soil spectra may reflect the presence of particular soil components, making satellite derived spectra suitable predictors of SOC. Bare soil spectra derived from Sentinel-2 were used to estimate SOC concentration (SOC%) and granulometric fractions in the plough layer (0–30 cm) of agricultural parcels in northern Belgium. Thereafter, the estimation performance of SOC% was compared for three DSM models: one with bare soil spectra, one with environmental covariates (topography, granulometry and vegetation), and a combined model with bare soil spectra and environmental covariates. The estimation performance of sand, silt and clay fractions using bare soil spectra from the spring seedbed (R2: 0.53–0.74; RPD: 1.49–2.05; RPIQ: 1.52–2.39) was higher than that of SOC% (R2: 0.16; RPD: 1.08; RPIQ: 1.32). The highest estimation performance of SOC% was obtained for a DSM model including all covariates (R2: 0.28; RPD: 1.18; RPIQ: 1.44), but the contribution of spring seedbed spectra to a model containing environmental covariates was small. The results provide valuable insights for refining soil property estimation using DSM with spectral and environmental covariates.