No abstract

No abstract

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.

The mapping of soil attributes provides support to agricultural planning and land use monitoring, which consequently aids the improvement of soil quality and food production. Landsat 5 Thematic Mapper (TM) images are often used to estimate a given soil attribute (i.e., clay), but have the potential to model many other attributes, providing input for soil mapping applications. In this paper, we aim to evaluate a Bare Soil Composite Image (BSCI) from the state of São Paulo, Brazil, calculated from a multi-temporal dataset, and study its relationship with topsoil properties, such as soil class and geology. The method presented detects bare soil in satellite images in a time series of 16 years, based on Landsat 5 TM observations. The compilation derived a BSCI for the agricultural sites (242,000 hectare area) characterized by very complex geology. Soil properties were analyzed to calibrate prediction models using 740 soil samples (0–20 cm) collected of the area. Partial least squares regression (PLSR) based on the BSCI spectral dataset was performed to quantify soil attributes. The method identified that a single image represents 7 to 20% of bare soil while the compilation of the multi-temporal dataset increases to 53%. Clay content had the best soil attribute prediction estimates (R2 = 0.75, root mean square error (RMSE) = 89.84 g kg−1, and accuracy = 74%). Soil organic matter, cation exchange capacity and sandy soils also achieved moderate predictions. The BSCI demonstrates a strong relationship with legacy geological maps detecting variations in soils. From a single composite image, it was possible to use spectroscopy to evaluate several environmental parameters. This technique could greatly improve soil mapping and consequently aid several applications, such as land use planning, environmental monitoring, and prevention of land degradation, updating legacy surveys and digital soil mapping.

Reflectance of light across the visible, near-infrared and shortwave infrared (VIS-NIR-SWIR, 0.4–2.5 µm) spectral region is very useful for investigating mineralogical, physical and chemical properties of soils, which can reduce the need for traditional wet chemistry analyses. As many collections of multispectral satellite data are available for environmental studies, a large extent with medium resolution mapping could be benefited from the spectral measurements made from remote sensors. In this paper, we explored the use of bare soil composites generated from the large historical collections of Landsat images for mapping cropland topsoil attributes across the European extent. For this task, we used the Geospatial Soil Sensing System (GEOS3) for generating two bare soil composites of 30 m resolution (named synthetic soil images, SYSI), which were employed to represent the median topsoil reflectance of bare fields. The first (framed SYSI) was made with multitemporal images (2006–2012) framed to the survey time of the Land-Use/Land-Cover Area Frame Survey (LUCAS) soil dataset (2009), seeking to be more compatible to the soil condition upon the sampling campaign. The second (full SYSI) was generated from the full collection of Landsat images (1982–2018), which although displaced to the field survey, yields a higher proportion of bare areas for soil mapping. For evaluating the two SYSIs, we used the laboratory spectral data as a reference of topsoil reflectance to calculate the Spearman correlation coefficient. Furthermore, both SYSIs employed machine learning for calibrating prediction models of clay, sand, soil organic carbon (SOC), calcium carbonates (CaCO3), cation exchange capacity (CEC), and pH determined in water, using the gradient boosting regression algorithm. The original LUCAS laboratory spectra and a version of the data resampled to the Landsat multispectral bands were also used as reference of prediction performance using VIS-NIR-SWIR multispectral data. Our results suggest that generating a bare soil composite displaced to the survey time of soil observations did not improve the quality of topsoil reflectance, and consequently, the prediction performance of soil attributes. Despite the lower spectral resolution and the variability of soils in Europe, a SYSI calculated from the full collection of Landsat images can be employed for topsoil prediction of clay and CaCO3 contents with a moderate performance (testing R2, root mean square error (RMSE) and ratio of performance to interquartile range (RPIQ) of 0.44, 9.59, 1.77, and 0.36, 13.99, 1.54, respectively). Thus, this study shows that although there exist some constraints due to the spatial and temporal variation of soil exposures and among the Landsat sensors, it is possible to use bare soil composites for mapping key soil attributes of croplands across the European extent.