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.

Soil organic carbon (SOC) stocks are a remarkable property for soil and environmental monitoring. The understanding of their dynamics in crop soils must go forward. The objective of this study was to determine the impact of temporal environmental controlling factors obtained by satellite images over the SOC stocks along soil depth, using machine learning algorithms. The work was carried out in São Paulo state (Brazil) in an area of 2577 km². We obtained a dataset of boreholes with soil analyses from topsoil to subsoil (0–100 cm). Additionally, remote sensing covariates (30 years of land use history, vegetation indexes), soil properties (i.e., clay, sand, mineralogy), soil types (classification), geology, climate and relief information were used. All covariates were confronted with SOC stocks contents, to identify their impact. Afterwards, the abilities of the predictive models were tested by splitting soil samples into two random groups (70 for training and 30% for model testing). We observed that the mean values of SOC stocks decreased by increasing the depth in all land use and land cover (LULC) historical classes. The results indicated that the random forest with recursive features elimination (RFE) was an accurate technique for predicting SOC stocks and finding controlling factors. We also found that the soil properties (especially clay and CEC), terrain attributes, geology, bioclimatic parameters and land use history were the most critical factors in controlling the SOC stocks in all LULC history and soil depths. We concluded that random forest coupled with RFE could be a functional approach to detect, map and monitor SOC stocks using environmental and remote sensing data.

Soil color and mineralogy are used as diagnostic criteria to distinguish different soil types. In the literature, 350–2500 nm spectra were successfully used to predict soil color and mineralogy, but these attributes currently are not mapped for most Brazilian soils. In this paper, we provided the first large-extent maps with 30 m resolution of soil color and mineralogy at three depth intervals for 850,000 km² of Midwest Brazil. We obtained soil 350–2500 nm spectra from 1397 sites of the Brazilian Soil Spectral Library at 0–20 cm, 20–60, and 60–100 cm depths. Spectra was used to derive Munsell hue, value, and chroma, and also second derivative spectra of the Kubelka–Munk function, where key spectral bands were identified and their amplitude measured for mineral quantification. Landsat composites of topsoil and vegetation reflectance, together with relief and climate data, were used as covariates to predict Munsell color and Fe–Al oxides, and 1:1 and 2:1 clay minerals of topsoil and subsoil. We used random forest for soil modeling and 10-fold cross-validation. Soil spectra and remote sensing data accurately mapped color and mineralogy at topsoil and subsoil in Midwest Brazil. Hematite showed high prediction accuracy (R² > 0.71), followed by Munsell value and hue. Satellite topsoil reflectance at blue spectral region was the most relevant predictor (25% global importance) for soil color and mineralogy. Our maps were consistent with pedological expert knowledge, legacy soil observations, and legacy soil class map of the study region.