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 (R² = 0.75, root mean square error (RMSE) = 89.84 g kg⁻¹, 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.

The search for sustainable land use has increased in Brazil due to the important role that agriculture plays in the country. Soil detailed classification is related with texture attribute. How can one discriminate the same soil class with different textures using proximal soil sensing, as to reach surveys, land use planning and increase crop productivity? This study aims to evaluate soil texture using a regional spectral library and its usefulness on classification. We collected 3750 soil samples covering 3 million ha within strong soil class variations in São Paulo State. The spectral analyses of soil samples from topsoil and subsoil were measured in laboratory (400–2500 nm). The potential of a regional soil spectral library was evaluated on the discrimination of soil texture. We considered two types of soil texture systems, one related with soil classification and another with soil managements. The soil line technique was used to assess differentiation between soil textural groups. Soil spectra were summarized by principal component analysis (PCA) to select relevant information on the spectra. Partial least squares regression (PLSR) was used to predict texture. Spectral curves indicated different shapes according to soil texture and discriminated particle size classes from clayey to sandy soils. In the visible region, differences were small because of the organic matter, while the short wave infrared (SWIR) region showed more differences; thus, soil texture variation could be differentiated by quartz. Angulation differences are on a spectral curve from NIR to SWIR. The statistical models predicted clay and sand levels with R² = 0.93 and 0.96, respectively. Indeed, we achieved a difference of 1.2% between laboratory and spectroscopy measurement for clay. The spectral information was useful to classify Ferralsols with different texture classification. In addition, the spectra differentiated Lixisols from Ferralsols and Arenosols. This work can help the development of computer programs that allow soil texture classification and subsequent digital soil mapping at detailed scales. In addition, it complies with requirements for sustainable land use and soil management.

Intensive cropland expansion for an increasing population has driven soil degradation worldwide. Modeling how agroecosystems respond to variations in soil attributes, relief and crop management dynamics can guide soil conservation. This research presents a new approach to evaluate soil loss by water erosion in cropland using the RUSLE model and Synthetic Soil Image (spectroscopy technique), which uses time series remotely sensed environmental, agricultural and anthropic variables, in the southeast region of São Paulo State, Brazil. The availability of the open-access satellite images of Tropical Rainfall Measuring Mission (TRMM) and Landsat satellite images provided ten years of rainfall data and 35 years of exposed soil surface. The bare soil surface and agricultural land use were extracted, and the multi-temporal rainfall erosivity was assessed. We predict soil maps’ attributes (texture and organic matter) through innovative soil spectroscopy techniques to assess the soil erodibility and soil loss tolerance. The erosivity, erodibility, and topography obtained by the Earth observations were adopted to estimate soil erosion in four scenarios of sugarcane (Saccharum spp.) residue coverage (0%, 50%, 75%, and 100%) in five years of the sugarcane cycle: the first year of sugarcane harvest and four subsequent harvesting years from 2013 to 2017. Soil loss tolerance means 4.3 Mg ha⁻¹ exceeds the minimum rate in 40% of the region, resulting in a total soil loss of ~6 million Mg yr⁻¹ under total coverage management (7 Mg ha⁻¹). Our findings suggest that sugarcane straw production has not been sufficient to protect the soil loss against water erosion. Thus, straw removal is unfeasible unless alternative conservation practices are adopted, such as minimum soil tillage, contour lines, terracing and other techniques that favor increases in organic matter content and soil flocculating cations. This research also identifies a spatiotemporal erosion-prone area that requests an immediately sustainable land development guide to restore and rehabilitate the vulnerable ecosystem service. The high-resolution spatially distribution method provided can identify soil degradation-prone areas and the cropland expansion frequency. This information may guide farms and the policymakers for a better request of conservation practices according to site-specific management variation.

Existe a necessidade de avaliar a importância do relevo associado aos parâmetros espectrais de solos como base no mapeamento. O objetivo deste trabalho foi determinar um método de detecção de limites de solos por meio da interação de dados espectrais e formas de relevo. Foram percorridas 14 topossequências representativas de uma área de 13.000 ha próxima dos municípios de São Carlos e Araraquara, SP. As amostras foram caracterizadas pelos métodos convencionais de análise química e granulométrica. Posteriormente, foram obtidos dados espectrais de 400 a 2.500 nm. As informações do relevo foram obtidas pelo emprego de técnicas de geoprocessamento. Geraram-se o modelo digital de elevação do terreno e os mapas de declividade, de curvatura, de índice topográfico composto e de Potencial de Densidade de Drenagem. Ainda, procedeu-se à validação dos métodos pontual e espacial. Na primeira validação, os pontos classificados nas topossequências foram tomados como verdadeiros e contrastados com as informações contidas no mapa de solo pré-existente, com os dados de relevo e com os dados espectrais agrupados. Na validação em nível espacial, procuraram-se avaliar em que locais os diferentes métodos indicavam mudanças nos limites dos solos e comparar com as observações reais. Verificou-se que a análise de agrupamento com cluster evidenciou-se eficiente na discriminação das unidades de solos em topossequência, quando utilizados parâmetros espectrais do solo. Já o conjunto de parâmetros de relevo isoladamente não foi o mais adequado.