The Midwest region in Brazil has the largest and most recent agricultural frontier in the country where there is no currently detailed soil information to support the agricultural intensification. Producing large-extent digital soil maps demands a huge volume of data and high computing capacity. This paper proposed mapping surface and subsurface key soil attributes with 30 m-resolution in a large area of Midwest Brazil. These soil maps at multiple depth increments will provide adequate information to guide land use throughout the region. The study area comprises about 851,000 km2 in the Cerrado biome (savannah) in the Brazilian Midwest. We used soil data from 7908 sites of the Brazilian Soil Spectral Library and 231 of the Free Brazilian Repository for Open Soil Data. We selected nine key soil attributes for mapping and aggregated them into three depth intervals: 0–20, 20–60 and 60–100 cm. A total of 33 soil predictors were prepared using Google Earth Engine (GEE), such as climate and geologic features with 1 km-resolution, terrain and two new covariates with 30 m-resolution, based on satellite measurements of the topsoil reflectance and the seasonal variability in vegetation spectra. The scorpan model was adopted for mapping of soil variables using random forest regression (RF). We used the model-based optimization by tuning RF hyperparameters and calculated the scaled permutation importance of covariates in R software. Our results were promising, with a satisfactory model performance for physical and chemical attributes at all depth intervals. Elevation, climate and topsoil reflectance were the most important covariates in predicting sand, clay and silt. In general, for predicting soil chemical attributes, climatic variables, elevation and vegetation reflectance provided to be the most important of predictive components, while for organic matter it was a combination of climatic dynamics and reflectance bands from vegetation and topsoil. The multiple depth maps showed that soil attributes largely varied across the study area, from clayey to sandy, suggesting that less than 44% of the studied soils had good natural fertility. We concluded that key soil attributes from multiple depth increments can be mapped using Earth observations data and machine learning methods with good performance.

No abstract

Geospatial soil information is critical for agricultural policy formulation and decision making, land-use suitability analysis, sustainable soil management, environmental assessment, and other research topics that are of vital importance to agriculture and economy. Proximal and Remote sensing technologies enables us to collect, process, and analyze spectral data and to retrieve, synthesize, visualize valuable geospatial information for multidisciplinary uses. We obtained the soil class map provided in this article by processing and analyzing proximal and remote sensed data from soil samples collected in toposequences based on pedomorphogeological relashionships. The soils were classified up to the second categorical level (suborder) of the Brazilian Soil Classification System (SiBCS), as well as in the World Reference Base (WRB) and United States Soil Taxonomy (ST) systems. The raster map has 30 m resolution and its accuracy is 73% (Kappa coefficient of 0.73). The soil legend represents a soil class followed by its topsoil color.

Geotechnologies allow natural resources to be surveyed more quickly and cheaply than traditional methods. This paper aimed to produce a digital soil map (DSM) based on Landsat time series data. The study area, located in the eastern part of the Brazilian Federal District (Rio Preto hydrographic basin), comprises a representative basin of the Central Brazil plateau in terms of pedodiversity. A spectral library was produced based on the soil spectroscopy (from the visible to shortwave infrared spectral range) of 42 soil samples from 0–15 cm depth using the Fieldspec Pro equipment in a laboratory. Pearson’s correlation and principal component analysis of the soil attributes revealed that the dataset could be grouped based on the texture content. Hierarchical clustering analysis allowed for the extraction of 13 reference spectra. We interpreted the spectra morphologically and resampled them to the Landsat 5 Thematic Mapper satellite bands. Afterward, we elaborated a synthetic soil/rock image (SySI) and a soil frequency image (number of times the bare soil was captured) from the Landsat time series (1984–2020) in the Google Earth Engine platform. Multiple Endmember Spectral Mixture Analysis (MESMA) was used to model the SySI, using the endmembers as the input and generating a DSM, which was validated by the Kappa index and the confusion matrix. MESMA successfully modeled 9 of the 13 endmembers: Dystric Rhodic Ferralsol (clayic); Dystric Rhodic Ferralsol (very clayic); Dystric Haplic Ferralsol (loam-clayic); Dystric Haplic Ferralsol (clayic); Dystric Petric Plinthosol (clayic); Dystric Petric Plinthosol (very clayic); Dystric Regosol (clayic); Dystric Regosol (very clayic); and Dystric, Haplic Cambisol (clayic). The root mean squared error (RMSE) varied from 0 to 1.3%. The accuracy of DSM achieved a Kappa index of 0.74, describing the methodology’s effectiveness to differentiate the studied soils.