Terrain analysis is an important tool for modeling environmental systems. Aiming to use the cloud-based computing capabilities of Google Earth Engine (GEE), we customized an algorithm for calculating terrain attributes, such as slope, aspect, and curvatures, for different resolution and geographical extents. The calculation method is based on geometry and elevation values estimated within a 3 × 3 spheroidal window, and it does not rely on projected elevation data. Thus, partial derivatives of terrain are calculated considering the great circle distances of reference nodes of the topographic surface. The algorithm was developed using the JavaScript programming interface of the online code editor of GEE and can be loaded as a custom package. The algorithm also provides an additional feature for making the visualization of terrain maps with a dynamic legend scale, which is useful for mapping different extents: from local to global. We compared the consistency of the proposed method with an available but limited terrain analysis tool of GEE, which resulted in a correlation of 0.89 and 0.96 for aspect and slope over a near-global scale, respectively. In addition to this, we compared the slope, aspect, horizontal, and vertical curvature of a reference site (Mount Ararat) to their equivalent attributes estimated on the System for Automated Geospatial Analysis (SAGA), which achieved a correlation between 0.96 and 0.98. The visual correspondence of TAGEE and SAGA confirms its potential for terrain analysis. The proposed algorithm can be useful for making terrain analysis scalable and adapted to customized needs, benefiting from the high-performance interface of GEE.

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.

Although many Soil Spectral Libraries (SSLs) have been created globally, these libraries still have not been operationalized for end-users. To address this limitation, this study created an online Brazilian Soil Spectral Service (BraSpecS). The system was based on the Brazilian Soil Spectral Library (BSSL) with samples collected in the Visible–Near–Short-wave infrared (vis–NIR–SWIR) and Mid-infrared (MIR) ranges. The interactive platform allows users to find spectra, act as custodians of the data, and estimate several soil properties and classification. The system was tested by 500 Brazilian and 65 international users. Users accessed the platform (besbbr.com.br), uploaded their spectra, and received soil organic carbon (SOC) and clay content prediction results via email. The BraSpecS prediction provided good results for Brazilian data, but performed variably for other countries. Prediction for countries outside of Brazil using local spectra (External Country Soil Spectral Libraries, ExCSSL) mostly showed greater performance than BraSpecS. Clay R² ranged from 0.5 (BraSpecS) to 0.8 (ExCSSL) in vis–NIR–SWIR, but BraSpecS MIR models were more accurate in most situations. The development of external models based on the fusion of local samples with BSSL formed the Global Soil Spectral Library (GSSL). The GSSL models improved soil properties prediction for different countries. Nevertheless, the proposed system needs to be continually updated with new spectra so they can be applied broadly. Accordingly, the online system is dynamic, users can contribute their data and the models will adapt to local information. Our community-driven web platform allows users to predict soil attributes without learning soil spectral modeling, which will invite end-users to utilize this powerful technique.