 This paper explored the use of Sentinel-2, S-2, temporal mosaics of bare soil, S-tube soil, over six years combined with soil moisture products, SMPs, derived from Sentinel-1 and 2 images, soil organic carbon, SOC measurement data and other environmental covariates derived from digital elevation models, lithology maps, and airborne gamma-ray data. The authors found that the best performing models included airborne gamma-ray thorium, slope and S-2 bands, bands 6, 7, 8, 8A, and indices, calcareous sedimentary rocks, Calcol. Additionally, the authors noted that 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 pedagogical point of view. This article was authored by Diego Urbina Salazar, Emanuel Vador, and C. Richard Aforges, and others.