Next Article in Journal
Remote Sensing Monitoring of Rice and Wheat Canopy Nitrogen: A Review
Next Article in Special Issue
Drought Propagation in Brazilian Biomes Revealed by Remote Sensing
Previous Article in Journal
Oceanic Mesoscale Eddies Identification Using B-Spline Surface Fitting Model Based on Along-Track SLA Data
Previous Article in Special Issue
Mining Is a Growing Threat within Indigenous Lands of the Brazilian Amazon
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Use of Airborne Radar Images and Machine Learning Algorithms to Map Soil Clay, Silt, and Sand Contents in Remote Areas under the Amazon Rainforest

by
Ana Carolina de S. Ferreira
1,
Marcos B. Ceddia
2,3,*,
Elias M. Costa
3,
Érika F. M. Pinheiro
2,
Mariana Melo do Nascimento
4 and
Gustavo M. Vasques
5
1
Instituto de Agronomia, Universidade Federal Rural do Rio de Janeiro, BR 465, km 7, Seropédica 23890-000, Brazil
2
Department of AgroTechnologies and Sustainability, Institute of Agronomy, Federal Rural University of Rio de Janeiro, BR 465, km 7, Seropédica 23890-000, Brazil
3
Laboratory of Water and Soils in Agroecosystem, Universidade Federal Rural do Rio de Janeiro, BR 465, km 7, Seropédica 23890-000, Brazil
4
Agronomic Engineering, Universidade Federal Rural do Rio de Janeiro, BR 465, km 7, Seropédica 23890-000, Brazil
5
Embrapa Soils, Rua Jardim Botânico 1024, Rio de Janeiro 22460-000, Brazil
*
Author to whom correspondence should be addressed.
Remote Sens. 2022, 14(22), 5711; https://0-doi-org.brum.beds.ac.uk/10.3390/rs14225711
Submission received: 18 September 2022 / Revised: 24 October 2022 / Accepted: 4 November 2022 / Published: 11 November 2022
(This article belongs to the Special Issue Remote Sensing of the Amazon Region)

Abstract

:
Soil texture has a great influence on the physical–hydric and chemical behavior of soils. In the Amazon regions, due to the presence of dense forest cover and limited access to roads, carrying out surveys and mapping of soils is challenging. When data exist, they are relatively sparse and the distribution is quite uneven. In this context, machine learning algorithms (ML) associated with remote sensor covariates offer a framework to derive digital maps of soil attributes. The objective of this study was to produce maps of surface and subsurface soil clay, silt, and sand contents in a 13.440 km2 area in the Amazon. The specific objectives were to (a) evaluate the gain in prediction accuracy when using the P-band of airborne radar as a covariate; (b) evaluate two sampling approaches (Reference Area—RA and Total Area—TA); and (c) evaluate the transferability and performance of three ML algorithms: regression tree (RT), random forest (RF), and support vector machine (SVM). The study site was divided into three blocks, called Urucu, Araracanga, and Juruá, respectively. The soil dataset consisted of 151 surface and subsurface sand, silt, and clay observations and 21 covariates (20 relief variables and the backscattering coefficient from the P-band). Both the RA and TA sampling approach used 114 observations for training the prediction models (75%) and 37 for validation (25%). The RA approach was better for the development of sand and silt models. Overall, RF derived the most accurate predictions for all variables. The effect of introducing the P-band backscattering coefficient improved the sand prediction accuracy at the surface and subsurface in Araracanga, which had the highest sand content, with relative improvements (RI) of the R2, root mean square error (RMSE), and mean absolute error (MAE) of 46%, 3%, and 4% at the surface, respectively, and 66.7%, 4.4%, and 5.2% at the subsurface, respectively. For silt, the P-band improved the predictions at the surface in Araracanga, which had the lowest silt contents among the blocks. For clay, adding the P-band improved the RF predictions at the subsurface, with RI of the R2, RMSE, and MAE of 29%, 5%, and 5%, respectively. Despite the low observation density, inherently hindered by the low accessibility of the area and high costs of sampling thereof, the results showed the potential of ML algorithms boosted by airborne radar P-band to map soil clay, silt, and sand contents in the Amazon.

1. Introduction

Soil texture is a fundamental physical property that strongly influences many other soil properties. The soil particle size fractions, namely clay, silt, and sand, influence soil fertility, water infiltration and retention capacity, soil organic matter dynamics, and, thus, the ability of soils to support plants, animals and life, and secure biodiversity [1,2,3]. Soil sand, silt, and clay contents are input data needed for most hydrological, climatic, and environmental models. They are also used to estimate hard-to-measure soil properties such as bulk density, hydraulic conductivity, and water-holding capacity [4,5].
The Brazilian Amazon rainforest represents a major challenge for the development of systematic soil mapping studies. The region covers an immense area (59% of the Brazilian territory) and has a large portion covered by dense evergreen forest [6,7]. Additionally noteworthy is the low density of roads, with most of the territory accessed only by boat and air transport. In this region, the constant presence of clouds makes it difficult to use satellite images and aerial photos obtained by passive (optical, infrared) remote sensors [7,8]. This condition makes active sensors, such as radar, potential alternatives for observing/surveying the land, serving as support for mapping environmental patterns and resources, including soils, hydrology, geology, and geomorphology. In fact, the climatic characteristics of the Amazon region and the intense land cover by native vegetation motivated the first project of systematic mapping of the Amazon region using radar images, the RADAM Project [9], which was a pioneering effort by the Brazilian government in the 1970s to survey natural resources using airborne radar imagery. At the time, the use of side-looking airborne radar (SLAR) represented a technological advance, because the radar images could be obtained both during the day and at night and in cloudy conditions, as radar microwaves penetrate most clouds. In the RADAM project [9], the X band was used (wavelengths close to 3 cm and frequency between 8 and 12.5 GHz) and image mosaics were generated at a scale of 1:250,000. Despite the advancement in the RADAM project as a source of important maps for the Brazilian Amazon region (geological, geomorphological, soil and vegetation maps), there is still a growing demand for more detailed maps of soil attributes to support projects for different purposes, including research in soil water and carbon [7].
Among the available radar bands, for soil studies in the Amazon region under native forest, the P band is ideal because the waves can pass the clouds and the tree canopies. Most of the radar research found in the literature concentrates on forestry studies [10,11,12,13,14,15,16]; however, recently there has been an increase in the application of radar remote sensing for soil assessment, mainly focusing on soil moisture [17,18,19,20,21,22]. As the dielectric behavior of the soil is affected by the particle size distribution, by assessing the soil dielectric properties, radar remote sensing indirectly assesses soil particle size distribution [23]. In the Brazilian Amazon region [8], the addition of relief and vegetation covariates derived from multispectral images with distinct spatial and spectral resolutions (Landsat 8 and RapidEye) and L-band radar images (ALOS PALSAR) were evaluated for the prediction of soil organic carbon stock (CS) and particle size fractions. Overall, the results showed that, even under forest coverage, the ALOS PALSAR L-band backscattering coefficient improved the accuracy of subsurface clay content predictions (8.2% higher) from regression kriging (RK) [8].
In addition to the limited availability of P-band radar images, especially in the Amazon, the execution of soil surveys in this region faces challenges inherent to its remoteness (low accessibility, little infrastructure, high transportation costs) [7,8]. Therefore, using existing data and knowledge from soil databases and previous surveys is essential to build predictive models of attributes such as soil particle size fractions. In this sense, the Reference Area (RA) approach in association with machine learning (ML) techniques becomes strategic. The RA approach assumes that a small area, if strategically chosen, can be surveyed to build a detailed soil map or soil prediction models with the potential to be extended or applied to other (ideally larger) areas with similar soil and landscape characteristics [24,25]. In this case, the RA approach would significantly reduce mapping costs, requiring only new field studies to assess the accuracy of the predictions in the new area.
On the other hand, as soil databases are limited in remote areas, the available data do not always present density and spatial distribution of soil observations that allow the use of techniques commonly used in digital soil mapping, such as models based on multivariate statistics and geostatistics. As an alternative, machine learning (ML) algorithms have been shown to be promising for mapping soil types and their attributes in large areas [1,4,26,27,28,29,30,31]. They refer to a large class of data-driven algorithms, some of which not following any statistical assumptions. As such, ML algorithms have the capacity of handling a large number of cross-correlated covariates (collinearity) as predictors [32].
The objectives of this study were to combine machine learning with remote sensor data to map soil surface and subsurface clay, silt, and sand fractions in the Brazilian Amazon, aiming specifically to (a) evaluate the gain in prediction accuracy from adding the P-band of airborne radar as covariate; (b) evaluate two sampling approaches (Reference Area—RA and Total Area—TA); and (c) evaluate the transferability and performance of regression tree (RT), random forest (RF), and support vector machine (SVM) models.

2. Materials and Methods

2.1. Study Area

The study area is located in the central region of the Amazonas state (at about 640 km from Manaus), covering an area of about 13.440 km2 between the municipalities of Coari and Tefé (Figure 1). The area is remote and practically all covered by equatorial Amazon rainforest. The elevation ranges from 23 to 112 m above mean sea level and the climate is equatorial (Af), according to Köppen classification, with the temperature of the coldest month higher than 20 °C, mean annual precipitation of 2500 mm, and no pronounced dry period.
According to ref. [8], most soils in the region have low base content, high aluminum content, and medium-to-high sand content. Some soils in the region have hydromorphic characteristics, especially those close to the floodplain of water courses and flat tops. The study area was divided into three blocks, which represent the petroleum exploration blocks by Petrobras (Brazilian Oil Company), namely Urucu (~4514 km2), Araracanga (~3751 km2), and Juruá (~4703 km2), respectively (Figure 1A). The project database comprises data from 151 soil profiles surveyed in two field campaigns (year 2008 and 2018, respectively).

2.2. Soil Sampling Designs

The development of soil prediction models and maps involves financial and logistical investments to support field soil surveys and laboratory and office work. Field sampling in the Amazon is restricted by the low accessibility due to the absence of roads and limited or no infrastructure to provide essential goods and services (e.g., lodging, food, and medical services). This characteristic of the region makes the execution of soil surveys complex, especially the more detailed ones.
The Reference Area (RA) for the study was the Geólogo Pedro de Moura Support Base (BOGPM), which belongs to Petrobras (Petróleo Brasileiro S.A.) and spans across circa 80 km2. The area is only accessed by air or river transport. In 2008, a detailed soil survey was carried out at BOGPM. In this area, in addition to the soil map, a database was organized containing 114 observations that included soil taxonomic class (Table 1), chemical and physical, as well as co-located relief covariates. From these data, prediction models of soil types and attributes have been developed for other areas, considering the BOGPM as an RA.
As an RA, the BOGPM serves as a base for soil sampling, for understanding the soil–landscape relationships of the region, and for training the prediction models aiming to transfer this knowledge and derived models to a larger region expanding the soil maps and its attribute maps to remote areas at a lower cost. However, the use of the RA approach assumes that the soil and landscape data observed in the RA represent the new areas where the prediction models are intended to be applied for deriving digital maps of soils and their attributes.
In 2018, a field campaign was carried out to visit 37 new soil sites as model and map validation sites for the RA approach. In this campaign, 16 remote clearings that allowed the landing and take-off of helicopters were identified. At each clearing, soil sites located within a 2000 m buffer were visited and sampled, expanding the original soil database from 114 to 151 soil profiles (Figure 1; Table 1).
With this data set, two sampling approaches were tested to develop soil clay, silt, and sand content prediction models for the whole area (13.440 km2), which encompasses three exploration blocks (Figure 1). It is important to note that, for purposes of organizing the cartographic bases, the area was divided into exploration blocks by Petrobras. In this study, the same logic was followed for prediction and map generation. Thus, throughout this study, the names adopted for each block will be used (Urucu, Araracanga, and Juruá, as presented in Figure 1). In the first approach—Reference Area—all 114 soil profiles occurring in the RA (Figure 1C) were used for model training, while the other 37 samples outside the RA were used for external validation of the models and maps. In the second approach—Total Area—the existence of an RA was ignored inasmuch as all 151 samples were pooled together, and the 114 training (75% of the samples) and 37 validation samples (25%) were randomly drawn from the pooled database of 151 samples.
The methodological strategy to predict sand, silt and clay for each soil depth (surf and sub) is presented in the flowchart (Figure 2).

2.3. Soil Particle Size Fractions

During the soil survey, the soil profiles were described morphologically with the separation of horizons/layers (A, AB, BA, B, C, AC, and CB, for example). For each of the horizons/layers, samples were collected for chemical and physical analyses. The sand, silt, and clay contents were determined from these samples using the Pipette method [35]. The dataset with values of sand, silt, and clay of what is called the surface layer (surf) is the weighted average of these fractions at horizons A, AB, AC, and AE (0–30 cm), while the dataset of sand, silt, and clay of the subsurface layer refers to the weighted average of these fractions in the BA, BE, and B horizons (0–100 cm) (Equation (1)). The values of the sand, silt, and clay fractions in the BC, CB, and C horizons/layers were not considered in the calculation, whereas CA and C were included when there was no B horizon, that is, for soils such as Quartzipsamments and Fluvents.
PSF surf / sub = i = 1 n PSF i T i / i = 1 n T i
where: PSFSurf/Sub is the particle size fraction (clay, silt, or sand content) in the desired layer (surface or subsurface), in g kg−1; PSFi is the PSF at horizon i, in g kg−1; Ti is the thickness, in m, of the portion of the horizon i that lies within the de-sired layer; and n is the number of horizons that have a portion within the desired layer.

2.4. Radar-Derived P-Band and Relief Covariates

The use of a radar sensor is important in the Amazon region mainly due to atmospheric conditions that include long rainy periods and the presence of clouds that often limit the use of passive remote sensors. The exclusive use of P-band (72 cm wavelength) microwave radar images in large regions covered by dense vegetation, such as the Amazon rainforest, is essential to generate thematic and relief maps. The longer wavelengths (P-band) can penetrate treetops and generate sufficiently strong reflections from the terrain below them to be more sensitive to biomass variations than other bands such as X, C, and L, and can be used to generate Digital Elevation Models (DEM).
A mosaic and a DEM of the study area were obtained from 84 Synthetic Aperture Interferometric Radar OrbiSAR-1 images, developed by Orbisat. All appropriate treatments were carried out, aiming to derive a mosaic and a DEM without interpolation failures, resulting in a hydrologically consistent DEM with 20 m spatial resolution. Primary and secondary relief derivatives were derived from the DEM using SAGA GIS version 7.7.0 [36], including Convergence Index, Topographic Wetness Index, Relative Slope Position, Channel Network Distance, Channel Network Base Level, LS-factor, Multiresolution Index of Valley Bottom Flatness, Multiresolution Index of the Ridge Top Flatness, Convexity Index, Aspect, Landforms, Profile Curvature, Plan curvature, Valley Depth, Slope Height, Mid Slope Position, Slope Gradient, Melton Ruggedness Number, and Flow Accumulation. All the data layers were brought to the same projection in ArcGIS (ESRI, Redlands, CA, USA).
The backscatter coefficient (σ°) of the HH polarization of the P-band was derived from the radar image mosaic. Reflector points in the ground were used for radiometric calibration. All calibration and radiometric corrections were performed using ENVI (L3Harris Geospatial, Broomfield, CO, USA).

2.5. Covariate Selection

The development of prediction models is a complex process that involves several steps. In the specific case of developing prediction models based on ML algorithms, as highlighted by ref. [32], conventionally, the choice of covariates is based on minimizing errors in input and output values. That is, a priori, no conceptual model of soil processes is contextualized. Only the processes that are transmitted by the input data are represented on the map.
In this study, two ways of covariate selection to develop ML models were tested: the “wrapper method” (WM) and “previous covariate selection” (PCS). In the first case (WM), all the covariates were made available for the training of the ML algorithms. In the second case (PCS), two steps were followed: (a) evaluation of Pearson’s correlation between particle size fractions and relief covariates, preferentially keeping the covariates with highest correlations; and (b) expert pedological knowledge was used to choose which covariates to keep on, in a case-by-case basis, aiming to better explain the soil–relief–vegetation relationships (SRV) in the region, as proposed by ref. [7].
The existence of multicollinearity was also considered both in the WM ant PCS method to make the covariates available to ML algorithms. The assessment of multicollinearity, which assesses the increase in variance due to the presence of multicollinearity, was performed based on the Variance Inflation Factor (VIF) [37], preferably keeping the covariates with VIF < 10 Equation (2).
VIF = 1 / 1 R 2

2.6. Dissimilarities in Covariates between the Reference Area and Total Area

The similarity of the landscape between the areas is important for the adequate transferability of the models. To examine the constraining effect of the relief characteristics on the transfer of the models between the reference area and the Urucu Araracanga and Juruá blocks, the descriptive statistics of the covariates were compared and the Gower similarity index (GSI; Equation (3)) [38,39] was calculated between the RA and each block, respectively.
Sij = 1 p k = 1 p 1 x ik x jk range   k
where Sij is the GSI between sites i and j; k represents the relief variables; p is the number of variables; range k is the range of variable k.
The GSI ranges between 0 and 1. A value of 1 means maximum similarity between the sites, that is, that the sites differ in no variable, whereas 0 means that the sites differ maximally in all their variables. In the literature, the GSI is generally used in its inverted form (1—GSI), or the Gower Dissimilarity Index (GDI). In this case, the interpretation is the opposite, that is, GDI values close to 0 mean that the two sites are similar, whereas values close to 1 mean that they are dissimilar in their variables. The GDI (1—GSI) was calculated from the relief covariates plus the backscatter coefficient derived from the radar images.

2.7. Model Training

The soil surface and subsurface sand, clay, and silt contents were modeled by regression tree (RT) [40], random forest (RF) [41], and support vector machine (SVM) [42]. The regression tree represents a set of rules over a hierarchical sequence for the purpose of partitioning the data. Its most important feature is the ability to convert complex decision processes into a series of simple decisions [40]. The purpose of RT is to separate observations into smaller and homogeneous groups in relation to the result of interest, such as soil class or attributes [40].
Random forest consists of a large number of individual RT models trained from bootstrap samples of the data [41]. The results of all individual trees are aggregated to make a single prediction. This method can also rank the predictor variable’s relative importance based on the regression prediction error of out-of-bag (OOB) predictions [41].
Support vector machine aims to determine decision limits among categories or continuous values by fitting optimal hyperplanes in the feature space that separates the samples minimizing prediction errors [42]. It can be used for classification and regression tasks. Table 2 summarizes the hyperparameters of each ML algorithms used in this study, R software environment [43].

2.8. Evaluation of the Accuracy of Interpolation Methods

The coefficient of determination (R2; Equation (4)) was used to evaluate the goodness-of-fit of the RT, RF, and SVM models for soil sand, clay, and silt content, and the mean absolute error (MAE; Equation (5)), and the root mean square error (RMSE; Equation (6)) were used to assess their prediction accuracy.
R 2 = 1 i = 1 n O i P i 2 i = 1 n ( O i O ¯ ) 2
MAE = 1 n   i = 1 n O i P i
RMSE = 1 n i = 1 n ( O i P i ) 2
where n is the number of observations, Oi and Pi are the observed and predicted values, respectively, and O ¯ is the mean of observed values.

2.9. Evaluation of the Importance of P-Band to Model’s Performance

To evaluate the importance of adding the backscattering coefficient of the P-band in the model, the Relative Improvements (RI) of the R2, RMSE, and MAEwere calculated, respectively (Equation (7)).
RI = Accuracy In Accuracy Out Accuracy Out × 100
where: RI is the relative improvement, in %, accuracy is the R2, MAE, or RMSE, respectively, in is the error value using the P-band, and Out is the error value without using the P-band.
The evaluation of the importance of the P-band was made for the ML models with the best performance and the covariate selection method with the best result. It was also evaluated according to the best approach (RA or TA) for each soil attribute.

3. Results

3.1. Summary Statistics

The soil sand, silt, and clay particle size fractions at the surface and subsurface layers present a frequency distribution similar to the standard normal (both skewness and excess kurtosis approximately 0), except surface clay (Table 3, whole dataset). The training and validation datasets follow the same pattern (close to normal distribution), differing in terms of minimum and maximum values, which is expected due to data partitioning. Based on the mean and median values of the particle size fractions, taken together, the textural classes vary from loan at the surface to clay loam at the subsurface. The mean and median values of sand, silt, and clay in the validation data dataset of the RA approach (V(RA)) indicate that the soils visited in remote areas outside the reference area (accessed from the 16 clearings) present the same textural classes as those observed in the reference area.
The large coefficients of variation (CV) values (>28%) characterize the heterogeneity of sample sets in both training and validation datasets. The range of sand, silt, and clay values was high. Sand contents ranged from 80 to 918 g kg−1 and from 44 to 855 g kg−1 at the surface and subsurface layers, respectively.
Clay contents had similar amplitude in the two layers (4.67 to 500 on the surface and 13 to 573 on the subsurface); however, in average terms, the clay contents in the subsurface practically doubled in relation to the surface (from 152 to 308 g kg−1). In the opposite direction, both the average levels of sand and silt tended to decrease with increasing depth (from 458 to 353 g kg−1 for sand and from 389 to 339 g kg−1 for silt).
The feasibility of prediction models that are based on the RA approach depends on the transferability of these models to other target areas. Thus, the statistics of the validation data of sand, silt, and clay in the three blocks (Urucu—VU, Araracanga—VA, and Juruá—VJ) separately allow a view of the similarity of the soils. The RA is located in the Urucu block, and the ideal is that the training data used there captures the great diversity of values found in all blocks. Comparing the minimum, maximum, and average sand values in the Araracanga (VA) block, both on the surface and in the subsurface, it is noted that in this region the soils had higher sand values than in the Urucu and Juruá blocks. In the first case (surface), the average sand (507 g kg−1) was 19% higher than in the Urucu block (425 g kg−1), while in the subsurface this difference was even greater (34%, 459 g kg−1 in Araracanga and 342 g kg−1 in Urucu). The statistics of silt data for the Juruá block (VJ) highlight the significant superiority of this fraction, both on the surface and in the subsurface, in relation to the other blocks. Specifically, in relation to the Urucu block (VU), the average value of silt in Juruá was 40% higher (668 g kg−1 against 476 g kg−1) and 34% higher (480 g kg−1 against 359 g kg−1), considering the surface and subsurface layers, respectively.
Additionally, in the Juruá block, the average clay content was 35% lower on the surface in relation to the data from the Urucu block (64 g kg−1 against 98 g kg−1). However, in the subsurface this relationship was reversed. The average clay content was 10% higher (327 g kg−1) than that found in the Urucu block (298 g kg−1). This inversion explains another distinction in the clay data of the Juruá block in relation to the other blocks. In Juruá, the average value of clay in the subsurface layer was 5 times higher than on the surface (64 g kg−1 against 327 g kg−1). In the other blocks, the increase in clay content with increasing depth was also marked but reached lower rates (3 and 2 times higher in the Urucu and Araracanga blocks, respectively).
Analyzing the statistics of sand, silt, and clay content of the validation dataset (V) using the TA approach (dataset 2), it is noted that differences of the average values in relation to training dataset (T) were lower. Only the average values of sand and silt, both at subsurface layer, presented values 10% higher than in the training dataset. In the first case (sand at subsurface) the average value was 19% higher (337 g kg−1 against 402 g kg−1). In the second case (silt at subsurface), the average values of the validation dataset were 11% lower than the training dataset (309 g kg−1 against 349 g kg−1).
Considering the evaluation of the statistics of the different granulometric fractions, in the different depths and approaches (RA and TA), it can be considered that the data present a frequency distribution close to the standard normal and that the textural classes of the soils of the reference area and the other regions visited present the same textural class (Loan and Clay loan). However, there were greater differences between the mean values of the sand, silt, and clay fractions of the validation dataset in relation to the training dataset when using the RA approach. The effect of these differences on the development and validation of prediction models is presented below, as well as the relationship between the granulometric fractions and the relief and radar covariates.

3.2. Similarity among the Reference Area and Exploration Blocks

Table 4 presents the statistics of the prediction covariates. Comparing the data between the blocks, it is noted that the region of the Juruá block was the one with the greatest discrepancy in relation to the reference area. Some relief covariates in the Juruá block had very different minimum, maximum, average, and median values compared with the RA, which reinforce the dissimilarity between these landscapes (Table 4). The covariates CNBL, CND, MRRTF, and MRVBF stand out as those with the most different relief statistics in the Juruá region in relation to RA (Table 4).
In Figure 3A–C, graphs are presented with the general GDI (red bars) and the same index for each covariate (gray bars). According to the GDI values, the RA and the Urucu, Araracanga, and Juruá blocks were similar in their relief variables, with GDI values of 0.155, 0.164, and 0.171, respectively (Figure 3). The dissimilarity increased by about 10% departing from the Urucu block towards Juruá (farthest from the reference area). The areas with the highest GDI were those associated with lowland areas (hydromorphic lowlands—black arrows on maps) and higher regions located at watershed upper boundaries (pixels with more discrete values highlighted with blue arrows on maps). The relief covariates that contributed most to differentiate the blocks in relation to RA were MRVBF, MRRTF, and RSP. These covariates were also the ones that had the highest correlations with the soil particle size fractions under study (Figure 4). From the results seen in Figure 3, the GDI can be used to both support the choice or to change a previously selected RA. In this study, the RA was imposed because it is the only accessible area in the region. However, it is possible to conjecture that if we were to change the RA, this change should be in the sense of including regions that expand the expression of the covariates that most differentiated the exploration blocks in relation to the RA (in this case, MRVBF, MRRTF, and RSP). It is important to highlight that these areas are the most difficult to access and cause the most undersampling in these environments.
From the data obtained (Table 4 and Figure 3), it appears that although there were differences in the statistics of the covariates of the blocks in relation to the reference area, the Gower index of similarity showed that the blocks had a very low dissimilarity value, indicating that the models developed in the reference area have the potential to be transferred to other areas.

3.3. Remote Sensing Covariates and Soil Particle Size Fractions Relationships

In both the RA and TA training datasets, all covariates had correlations lower than 0.50 against soil particle size fractions (Figure 4). In the RA dataset (Figure 4A), the highest correlation values for each particle size fraction were found between the topographic wetness index (TWI) and clay at the surface (−0.47) and subsurface (−0.45), surface silt (0.33) and multiresolution index of ridge top flatness (MRRTF), subsurface silt and TWI (0.30), and relative slope position (RSP) and sand at the surface (−0.26) and subsurface (−0.35). In the TA dataset (Figure 4B), the highest correlations were surface clay against slope (0.39) or TWI (−0.39), subsurface clay against TWI (−0.32), surface silt against channel network base level (CNBL) (0.49) or MRRTF (0.49), subsurface silt against TWI (0.44), surface sand against CNBL (−0.34), and subsurface sand against CNBL (−0.28). Overall, sand content had the lowest correlations against remote sensing covariates.
The results of the general Gower index (Figure 3) showed that there was little dissimilarity between the RA and the Urucu, Araracanga, and Juruá blocks, with GDI (values of 0.155, 0.164, and 0.171, respectively). However, even though these dissimilarity values are low, most of the covariates that had higher correlations (Figure 4) also had greater contributions of dissimilarity index values in relation to the general Gower index (RSP, CI, MRVBF, MRRTF, LF) (Figure 3).
Importance of predictor covariates for the attributes evaluated in the RF model is seen in (Figure 5).
The source material, relief, vegetation, and climate act in tandem to explain the spatial distribution patterns of soil types in the region. These same covariates were contextualized in the soil–relief–vegetation model (SRV) (Figure 6).
In general, the covariates convergence index (CI), landforms (LF), radar P-band backscattering coefficient (P-band), profile curvature (ProfC), RSP, and MRRTF were the most important for sand prediction by RF (Figure 5A,B). The CI represents the behavior of the surface runoff, which was influenced by the shapes of the terrain, represented by LF. The sand contents were higher close to river channels, where CI values were negative, meaning converging terrains towards lowland channels. Positive CI values indicate divergent areas, where well-drained tops and flatter slopes predominate, from which surface runoff occurs in all directions. In these areas the sand contents were lower. The RSP was applied to identify topographical features and its values ranged from 0 to 1. The values closer to 0 were characterized by lowland regions, that is, the V- and U-shaped valleys, which have high levels of sand. Values closer to 1 represent upper slopes and ridge tops with low sand contents. The profile curvature (ProfC) expressed the difference between convex curvatures of the concave ones, influencing the surface flow velocity from the higher to the lower parts (Figure 6). It also allowed greater distinction between well-drained soils on ridge tops (convex surfaces) and imperfectly drained soils on concave to flat surfaces, for instance in V- and U-shaped valleys.
The covariates MRRTF, TWI, multiresolution index of valley bottom flatness (MRVBF), and ProfC had positive correlations with silt. The flat tops on the uplands were represented by high MRRTF values, whereas the valley bottoms had the highest MRVBF values. These covariates, associated with TWI, characterize the spatial distribution of soil saturation zones, adding important information to locate hydromorphic soils. In turn, these zones had the higher silt contents (Figure 6) and are the zones where the lowlands (MU2—Aquents, Aquepts) and uplands with flat tops (MU4—Aquults, Aquents) occur. Again, the ProfC helped to separate the areas of well-drained soils (convex surfaces) from those with imperfectly drained ones (concave to flat surfaces), mainly at the subsurface.
The combination of the slope and TWI covariates allowed identifying the regions with the highest clay contents, where the MU1 (Ultisols, Inceptisols) and MU3 (Ultisols, Inceptisols) units are found (Figure 6). The MU1 regions were represented by steeper slopes generally closer to large drainage networks where the slope influences the speed of surface and subsurface flows. The slope has great potential to help in the identification of Ultisols areas where the highest clay contents predominate. The MU3 unit occurs on well-drained tops with smoother slopes and relatively flat to smoothly wavy relief with good drainage, also with high clay contents (Figure 6).

3.4. Model Prediction Performance

Random forest derived the best predictions, with the least errors, for all soil particle size fractions at both layers, followed by SVM (Table 5, Table 6, Table 7 and Table 8). Regression tree is the simplest among the three methods tested. It creates a series of decision rules based on the covariates to make a prediction at a terminal leaf. As such, it was uncapable of outperforming RF, which is a combination of RTs, and SVM. On the other hand, RF outperformed SVM, meaning that decision rules derived from a series of RTs are better than a single set of hyperplanes. In fact, in general the prediction errors were more similar between RT and SVM than between SVM and RF. Favoring RF is the fact that it uses random selections of covariates and training and validation (OOB) sets for building each tree, which control overfitting minimizing validation errors.
The fitted model R2 varied from 0.34 to 0.62 for RT models, from 0.91 to 0.95 for RF, and from 0.39 to 0.81 for SVM models. The validation RMSE, considering all 37 validation samples, varied across all sampling approach and methods of covariate selection in the ranges of 144 to 198 for the surface sand and 162 to 202 (g kg−1) for the sand subsurface layer. For silt, the range was from 141 to 182 at the surface layer and from 89 to 102 (g kg−1) at the subsurface layer. The RMSE range of clay was from 65 to 111 at the surface and 107 to 141 (g kg−1) at the subsurface layer.
The RA sampling approach outperformed the TA approach for the surface and subsurface sand and silt contents, whereas surface and subsurface clay contents were best predicted using TA approach. The PCS covariate selection method was the best option to predict surface sand, and surface and subsurface silt and clay contents, whereas WM was the preferred choice only for subsurface sand prediction.

3.5. Relative Improvement (RI%) from Adding the Radar P-Band

Considering the combination of best results (the algorithms—RF, RT and SVM, the approach—RA or TA, and the covariate selection method—WM or PCS), the gain in accuracy of the models, with and without the P-band, was evaluated applying the RI index (%) on the R2, RMSE and MAE metrics at surface and sub surface layers (Figure 7 and Figure 8 respectively). Considering the surface layer (Figure 7), in the prediction of sand and silt, the RA approach had better results and so the metrics were separated by blocks (Urucu, Araracanga, and Jurua), how much the P-band influences the accuracy when the model generated in the RA is transferred to other blocks was evaluated. For clay, as the TA approach performed better, the metrics do not distinguish between blocks. Note that the introduction of the P-band had a greater effect on the R2 results. For the sand fraction, the introduction of the P-band allowed the R2 (the proportion of the variation of a response variable is explained by the variation of other explanatory variables) to increase by 41%, 46%, and 24% for the Urucu, Araracanga, and Juruá blocks, respectively. However, when analyzing the RMSE and MAE metrics, the gain was low (<5%). In the case of silt, the introduction of P-band also increased R2, but to a lesser extent (7.4%, 12%, and 10.6% for Urucu, Araracanga, and Juruá, respectively). As in the case of sand, the change in the RMSE and MAE metrics for silt prediction was low (between 0% and 1.8%). In the case of the clay attribute, the introduction of the P band did not change the metric values (RI% = 0).
Analyzing the subsurface layer (Figure 8), the pattern observed on the surface was maintained. In other words, the use of radar images is important to generate maps of covariates (in this case, the relief and hydrographic attributes) under native forest cover; however, the effect of the backscatter coefficient with polarization HH, by itself, did not bring a significant gain (≥10%) in the accuracy of the models (RMSE and MAE). For example, adding the P-band improved the RF predictions of clay content at the subsurface layer, with RI of the R2, RMSE, and MAE of 29%, 5%, and 5%, respectively.

3.6. Soil Particle Size Fraction Maps

In the study area, sand contents ranged from 303 to 721 g kg−1 at the surface (Figure 9), and from 212 to 635 g kg−1 at the subsurface (Figure 10), decreasing slightly with depth. The lowest sand values were predicted in hydromorphic flat tops and areas with steeper slopes (Figure 9 and Figure 10). The highest levels of sand were present in the floodplain regions, close to the channels of the large rivers and streams, and on terraces around the main watercourse (U-shaped valleys). Large sand contents were also found in the more embedded valleys (V-shaped valleys) of slope regions. These environments are characterized by the accumulation of sandy sediments from natural erosive processes, making the lowlands clogged. In these areas, the predominant soils were classified as Aquents or Aquepts (MU2 unit).
Predicted silt contents varied from 209 to 577 g kg−1 at the surface (Figure 11), and from 215 to 517 g kg−1 at the subsurface (Figure 12). The largest silt contents were found in the areas of hydromorphic flat tops (Figure 11 and Figure 12). These areas usually occur at the highest elevations of the study area, at the upland watershed boundaries. Flat relief and insufficient drainage characterize these areas, where there is a predominance of Hapludults, Aquults, and Aquents (MU4 unit) (Figure 6). Relevant silt values were also found in lowland regions, where Aquents (MU2) occur.
Predicted clay contents ranged from 47 to 303 g kg−1 at the surface (Figure 13), and increased at the subsurface, ranging from 154 to 458 g kg−1 at the subsurface (Figure 14). The increase of clay with depth is consistent with the occurrence of Ultisols, which present a diagnostic argillic B horizon at the subsurface. The highest clay contents occur in areas with steep slopes and well-drained tops (Figure 13 and Figure 14). These regions were represented by the mapping units MU1 and MU3 where there is a predominance of Ultisols.

4. Discussion

The challenge of mapping soil fraction in the Amazon rainforest comes from the difficulties in obtaining soil data that are related with the fact that a major portion of the area is covered by a dense evergreen forest, a low density of roads, with most of the territory accessed only by boat and air transport. Additionally, the difficulty of obtaining data on representative environmental covariates, because of the constant presence of clouds in the region, makes it difficult to use satellite images and aerial photos obtained by passive (optical, infrared) remote sensors. Despite all these limitations, the results of this study illustrate the potential advantages of using ML algorithms associated with remote sensor covariates (terrain attributes and P-band of airborne radar) and RA approach to map particle size fractions in this region.
The comparison of these approaches highlighted that the non-linear model introduced significant improvements in the prediction of soil texture fractions and consequently ML are potentially superior to linear methods of spatial prediction of soil texture [44]. Additionally [45], ML algorithms, in this case Support Vector Regression (SVR), produced the best prediction accuracy compared with the geostatistical interpolation techniques. The results of this study, with best the prediction for the RF model, corroborate those of ref. [46], which also used radar data to estimate soil texture and obtained better results with RF than SVM. As already highlighted by refs. [8,47], the maximum silt values are relatively high when compared with the average contents found in Brazilian soils. According to ref. [8], in the Amazon region, silt greater than 400 g kg−1 are manly found in hydromorphic soils in the region of hydromorphic soils, which are not only found on lowlands but also in broad plateaus located in higher-altitude regions [7]. These regions have specific environmental characteristics (Figure 6) that allowed a good capture of patterns by the environmental covariates, which resulted in good prediction results for this fraction.
In general, both the correlation coefficients (Figure 4) and the most important covariates used to predict and map soil particle size fractions by RF (Figure 5) coincide with the hypotheses raised in ref. [9], as well as with previous studies in the region [7,8,48].
Some of these covariates also appear as important predictors of soil particle size fractions in ref. [49], where slope and TWI predictors had 80% of the importance for predicting surface clay (0 to 30 cm), and TWI and MRVBF were important covariates for silt prediction. In Iran, ref. [2] found TWI as one of the most important covariates for clay prediction, and similarly TWI and MRVBF were important ones for silt prediction.
The spatial patterns of the soil particle size fractions found in this study corroborate the results of ref. [8] carried out in the same study region.
A few studies have investigated the potential of P-band in mapping soil properties, most of them focus on the soil moisture and soil dielectric variations [20,22]. It is even rarer to study the P-band in the soil mapping or vegetation in the Brazilian Amazon [50]; for the authors, P-band data can make a substantial contribution to the development of models in tropical rainforest regions, especially in those areas where it is difficult to obtain data from optical sensors. Although it is not possible to compare the results with other studies, as there has been no work conducted on the use of P-band to predict soil texture, our results showed that it has great potential to improve the predictions of clay, silt, and sand fractions at the surface and subsurface, and new studies with more soil data are required to formulate better conclusions. Besides, if the VV polarization of the P-band image were available, perhaps it would be possible to extract greater knowledge of the interaction of the ratios and crosses of polarizations with granulometric fractions. For example, ref. [51], working with the X-Band, found that the sensitivity of soil texture is better observed at higher-incidence angles than lower-incidence angles in both polarizations, i.e., HH- and VV-pol. Besides, changes in soil texture are also sensitive to polarization and it was observed that VV-pol is more sensitive than HH-pol for different soil texture field. On the other hand, ref. [52], also working with the X-Band, found that a strong change in specular scattering coefficient is observed by changing the sand percentage in soil for HH polarization, while in the case of VV-polarization a lesser change is observed. It is difficult to observe the change in specular scattering coefficient with change in soil texture when the surface is considered as rough. Finally, the authors highlighted that it is important to minimize the roughness effect while observing the texture with specular scattering and that for higher-incidence angles (P50°), the distinction in soil texture fields are clearly observable on the basis of the copolarization ratio.
The Amazon region has peculiar characteristics that demand an enormous logistical, financial, and personal effort to survey soils. It is not by chance that the major soil surveys date from the 70′s and 80′s [9] and they are exploratory or reconnaissance types. Despite all the limitations imposed by the condition of the region, this study showed that the RA approach can reduce logistical, financial, and personnel costs. In addition, the use of covariables such as P-band, which is able to surpass the tree canopy and suffers little or no interference from clouds, combined with covariate selection methods and the training of robust ML algorithms can greatly increase the prediction results, producing more detailed and very useful maps.

5. Conclusions

This work investigated the use of remote sensing covariates derived from airborne synthetic aperture interferometric radar images to predict soil surface and subsurface sand, silt, and clay contents in the Brazilian Central Amazon. A Reference Area sampling design was proposed to reduce costs and expedite soil survey was contrasted against a random sampling design (that is, Total Area sampling), and combined with three machine learning methods (RT, RF, and SVM) and two covariate selection approaches (WM and PCS).
The RA approach was the best sampling option, deriving the least errors, for surface and subsurface silt and sand content prediction. Total Area random sampling was preferred for surface and subsurface clay content prediction, though the errors were similar to those from the RA approach. The RA was 80 km2, whereas the whole area to be mapped was 13.440 km2. This means that a tiny fraction of 0.6% of the total area served to collect soil and remotely sensed relief and P-band data to train soil particle size prediction methods, and transfer them to the whole area, composed by three relatively huge exploration blocks. Thus, the RA approach combined with remote sensing is recommended for expediting soil mapping and saving costs, especially in large areas.
From the relief attributes derived from the DEM, it was possible to establish relationships between the soil particle size fractions and the landscape. The selection of covariates (PCS) obtained, in general, better results than the all-in WM option that is commonly employed in digital soil mapping studies. The most important covariates to predict the soil particle size fractions in the Central Amazon region were CI, LF, MRRTF, MRVBF, TWI, slope, and ProfC for all fractions, in addition to the radar P-band backscatter coefficient for surface sand and clay contents.
Random forest outperformed RT and SVM for all soil particle size fractions and both layers. It is recommended for its robustness and ease to implement in free and open-source software. The P-band backscatter coefficient was considered an important covariate for the prediction of surface sand and clay contents by RF, showing its potential use for mapping these attributes.

Author Contributions

Conceptualization: M.B.C. and A.C.d.S.F.; methodology: M.B.C. and A.C.d.S.F.; software: A.C.d.S.F. and E.M.C.; validation: M.B.C.; formal analysis: M.B.C. and A.C.d.S.F.; investigation: M.B.C., A.C.d.S.F., É.F.M.P. and G.M.V.; resources: M.B.C.; data curation: M.B.C., A.C.d.S.F. and E.M.C.; writing—original draft preparation: A.C.d.S.F.; writing—review and editing: M.B.C., A.C.d.S.F., E.M.C., É.F.M.P., M.M.d.N. and G.M.V.; visualization: M.B.C., A.C.d.S.F. and E.M.C.; supervision: M.B.C.; project administration: M.B.C.; funding acquisition: M.B.C. All authors reviewed the manuscript. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Acknowledgments

The authors acknowledge the National Petroleum Agency (ANP) for funding the project “Digital Mapping of Soils in Oil and Gas Exploration and Production Areas—Case Studies of the North and Northeast Brazilian Fields” under the agreement number 5850.0105881.17.9 (PETROBRAS/FAPUR/UFRRJ) and Evaluation of Graduate Education (CAPES, finance code 001).

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Akpa, S.I.C.; Odeh, I.O.A.; Bishop, T.F.A.; Hartemink, A.E. Digital Mapping of Soil Particle-Size Fractions for Nigeria. Soil Sci. Soc. Am. J. 2014, 78, 1953–1966. [Google Scholar] [CrossRef] [Green Version]
  2. Mehrabi-Gohari, E.; Matinfar, H.R.; Jafari, A.; Taghizadeh-Mehrjardi, R.; Triantafilis, J. The Spatial Prediction of Soil Texture Fractions in Arid Regions of Iran. Soil Syst. 2019, 3, 65. [Google Scholar] [CrossRef] [Green Version]
  3. Dos Santos, H.G.; Jacomine, P.K.T.; Dos Anjos, L.H.; De Oliveira, V.A.; Lumbreras, J.F.; Coelho, M.R.; De Almeida, J.A.; De Araujo Filho, J.C.; De Oliveira, J.B.; Cunha, T.J.F. Sistema Brasileiro de Classificação de Solos; Embrapa: Brasília, Brazil, 2018; ISBN 978-85-7035-817-2. [Google Scholar]
  4. Ließ, M.; Glaser, B.; Huwe, B. Uncertainty in the Spatial Prediction of Soil Texture: Comparison of Regression Tree and Random Forest Models. Geoderma 2012, 170, 70–79. [Google Scholar] [CrossRef]
  5. Minasny, B.; Hartemink, A.E. Predicting Soil Properties in the Tropics. Earth-Sci. Rev. 2011, 106, 52–62. [Google Scholar] [CrossRef]
  6. Amazônia Legal | IBGE. Available online: https://www.ibge.gov.br/geociencias/organizacao-do-territorio/estrutura-territorial/15819-amazonia-legal.html?=&t=o-que-e (accessed on 21 October 2022).
  7. Ceddia, M.B.; Villela, A.L.O.; Pinheiro, É.F.M.; Wendroth, O. Spatial Variability of Soil Carbon Stock in the Urucu River Basin, Central Amazon-Brazil. Sci. Total Environ. 2015, 526, 58–69. [Google Scholar] [CrossRef]
  8. Ceddia, M.B.; Gomes, A.S.; Vasques, G.M.; Pinheiro, É.F.M. Soil Carbon Stock and Particle Size Fractions in the Central Amazon Predicted from Remotely Sensed Relief, Multispectral and Radar Data. Remote Sens. 2017, 9, 124. [Google Scholar] [CrossRef] [Green Version]
  9. BRASIL Departamento Nacional de Produção Mineral. Projeto Radambrasil, 1973–1987. (Levantamento de Recursos Naturais, 38 volumes).
  10. Santos, J.; Spinelli-Araujo, L.; Kuplich, T.; Da Costa Freitas, D.; Dutra, L.; Sant’Anna, S.; Gama, F. Tropical Forest Biomass and Its Relationship with P-Band SAR Data. Rev. Bras. De Cartogr. 2006, 1, 58. [Google Scholar]
  11. Gama, F.; Santos, J.; Mura, J.; Rennó, C. Estimativa de Parâmetros Biofísicos de Povoamentos de Eucalyptus Através de Dados SAR Estimation of Biophysical Parameters in the Eucalyptus Stands by SAR Data. Ambiência 2009, 2, 29–42. [Google Scholar]
  12. Sambatti, J.B.M.; Leduc, R.; Lubeck, D.; Moreira, J.R.; dos Santos, J.R. Assessing Forest Biomass and Exploration in the Brazilian Amazon with Airborne InSAR: An Alternative for REDD. Open Remote Sens. J. 2012, 5, 21–36. [Google Scholar] [CrossRef]
  13. Saatchi, S.; Marlier, M.; Chazdon, R.L.; Clark, D.B.; Russell, A.E. Impact of Spatial Variability of Tropical Forest Structure on Radar Estimation of Aboveground Biomass. Remote Sens. Environ. 2011, 115, 2836–2849. [Google Scholar] [CrossRef]
  14. Freitas, C.D.C.; Soler, L.D.S.; Sant’Anna, S.J.S.; Dutra, L.V.; dos Santos, J.R.; Mura, J.C.; Correia, A.H. Land Use and Land Cover Mapping in the Brazilian Amazon Using Polarimetric Airborne P-Band SAR Data. IEEE Trans. Geosci. Remote Sens. 2008, 46, 2956–2970. [Google Scholar] [CrossRef]
  15. Neeff, T.; Dutra, L.; Santos, J.; Da Costa Freitas, D.; Spinelli-Araujo, L. Tropical Forest Measurement by Interferometric Height Modeling and P-Band Radar Backscatter. For. Sci. 2005, 51, 585–594. [Google Scholar]
  16. Alemohammad, S.H.; Jagdhuber, T.; Moghaddam, M.; Entekhabi, D. Soil and Vegetation Scattering Contributions in L-Band and P-Band Polarimetric SAR Observations. IEEE Trans. Geosci. Remote Sens. 2019, 57, 8417–8429. [Google Scholar] [CrossRef]
  17. Zribi, M.; Sahnoun, M.; Baghdadi, N.; Le Toan, T.; Ben Hamida, A. Analysis of the Relationship between Backscattered P-Band Radar Signals and Soil Roughness. Remote Sens. Environ. 2016, 186, 13–21. [Google Scholar] [CrossRef] [Green Version]
  18. Zribi, M.; Sahnoun, M.; Dusséaux, R.; Afifi, S.; Baghdadi, N.; Ben Hamida, A. Analysis of P Band Radar Signal Potential to Retrieve Soil Moisture Profile. In Proceedings of the 2016 2nd International Conference on Advanced Technologies for Signal and Image Processing (ATSIP), Monastir, Tunisia, 21 March 2016; pp. 591–595. [Google Scholar]
  19. Blumberg, D.G.; Freilikher, V.; Kaganovskii, Y.; Maradudin, A.A. Subsurface Microwave Remote Sensing of Soil-Water Content: Field Studies in the Negev Desert and Optical Modelling. Int. J. Remote Sens. 2002, 23, 4039–4054. [Google Scholar] [CrossRef]
  20. Du, J.; Kimball, J.S.; Moghaddam, M. Theoretical Modeling and Analysis of L- and P-Band Radar Backscatter Sensitivity to Soil Active Layer Dielectric Variations. Remote Sens. 2015, 7, 9450–9472. [Google Scholar] [CrossRef] [Green Version]
  21. Etminan, A.; Tabatabaeenejad, A.; Moghaddam, M. Retrieving Root-Zone Soil Moisture Profile From P-Band Radar via Hybrid Global and Local Optimization. IEEE Trans. Geosci. Remote Sens. 2020, 58, 5400–5408. [Google Scholar] [CrossRef]
  22. Tabatabaeenejad, A.; Burgin, M.; Duan, X.; Moghaddam, M. P-Band Radar Retrieval of Subsurface Soil Moisture Profile as a Second-Order Polynomial: First AirMOSS Results. IEEE Trans. Geosci. Remote Sens. 2015, 53, 645–658. [Google Scholar] [CrossRef]
  23. Srivastava, H.S.; Patel, P.; Navalgund, R.R. Incorporating Soil Texture in Soil Moisture Estimation from Extended Low-1 Beam Mode RADARSAT-1 SAR Data. Int. J. Remote Sens. 2006, 27, 2587–2598. [Google Scholar] [CrossRef]
  24. Lagacherie, P.; Robbez-Masson, J.M.; Nguyen-The, N.; Barthès, J.P. Mapping of Reference Area Representativity Using a Mathematical Soilscape Distance. Geoderma 2001, 101, 105–118. [Google Scholar] [CrossRef]
  25. Taghizadeh-Mehrjardi, R.; Sheikhpour, R.; Zeraatpisheh, M.; Amirian-Chakan, A.; Toomanian, N.; Kerry, R.; Scholten, T. Semi-Supervised Learning for the Spatial Extrapolation of Soil Information. Geoderma 2022, 426, 116094. [Google Scholar] [CrossRef]
  26. Arruda, G.P.D.; Demattê, J.A.; Chagas, C.D.S.; Fiorio, P.R.; Fongaro, C.T. Digital Soil Mapping Using Reference Area and Artificial Neural Networks. Sci. Agric. 2016, 73, 266–273. [Google Scholar] [CrossRef] [Green Version]
  27. Grinand, C.; Arrouays, D.; Laroche, B.; Martin, M.P. Extrapolating Regional Soil Landscapes from an Existing Soil Map: Sampling Intensity, Validation Procedures, and Integration of Spatial Context. Geoderma 2008, 143, 180–190. [Google Scholar] [CrossRef]
  28. Chagas, C.D.S.; Junior, W.D.C.; Bhering, S.B.; Filho, B.C. Spatial Prediction of Soil Surface Texture in a Semiarid Region Using Random Forest and Multiple Linear Regressions. Catena 2016, 139, 232–240. [Google Scholar] [CrossRef]
  29. Bhering, S.B.; Chagas, C.D.S.; Junior, W.D.C.; Pereira, N.R.; Filho, B.C.; Pinheiro, H.S.K. Mapeamento digital de areia, argila e carbono orgânico por modelos Random Forest sob diferentes resoluções espaciais. Pesq. Agropec. Bras. 2016, 51, 1359–1370. [Google Scholar] [CrossRef] [Green Version]
  30. Wolski, M.S.; Dalmolin, R.S.D.; Flores, C.A.; Moura-Bueno, J.M.; Caten, A.T.; Kaiser, D.R. Digital Soil Mapping and Its Implications in the Extrapolation of Soil-Landscape Relationships in Detailed Scale. Pesqui. Agropecuária Bras. 2017, 52, 633–642. [Google Scholar] [CrossRef] [Green Version]
  31. Silva, S.H.G.; de Menezes, M.D.; Owens, P.R.; Curi, N. Retrieving Pedologist’s Mental Model from Existing Soil Map and Comparing Data Mining Tools for Refining a Larger Area Map under Similar Environmental Conditions in Southeastern Brazil. Geoderma 2016, 267, 65–77. [Google Scholar] [CrossRef]
  32. Wadoux, A.; Minasny, B.; Mcbratney, A. Machine Learning for Digital Soil Mapping: Applications, Challenges and Suggested Solutions. Earth-Sci. Rev. 2020, 210, 103359. [Google Scholar] [CrossRef]
  33. World Reference Base for Soil Resources 2014: International Soil Classification System for Naming Soils and Creating Legends for Soil Maps; FAO: Rome, Italy, 2014; ISBN 978-92-5-108369-7.
  34. Keys to Soil Taxonomy | NRCS Soils. Available online: https://www.nrcs.usda.gov/wps/portal/nrcs/detail/soils/survey/class/taxonomy/?cid=nrcs142p2_053580 (accessed on 4 September 2022).
  35. Infoteca-e: Manual de Métodos de Análise de Solo. Available online: https://www.infoteca.cnptia.embrapa.br/handle/doc/1085209 (accessed on 5 September 2022).
  36. SAGA—System for Automated Geoscientific Analyses. Available online: https://saga-gis.sourceforge.io/en/index.html (accessed on 4 September 2022).
  37. Garson, G.D. Testing statistical assumptions; Blue Book Series; Statistical Associates Publishing: Hillsborough, CA, USA, 2012. [Google Scholar]
  38. Gower, J.C. A General Coefficient of Similarity and Some of Its Properties. Biometrics 1971, 27, 857. [Google Scholar] [CrossRef]
  39. Mallavan, B.P.; Minasny, B.; McBratney, A.B. Homosoil, a Methodology for Quantitative Extrapolation of Soil Information Across the Globe. In Digital Soil Mapping: Bridging Research, Environmental Application, and Operation; Boettinger, J.L., Howell, D.W., Moore, A.C., Hartemink, A.E., Kienast-Brown, S., Eds.; Progress in Soil Science; Springer: Dordrecht, The Netherlands, 2010; pp. 137–150. ISBN 978-90-481-8863-5. [Google Scholar]
  40. Breiman, L.; Friedman, J.H.; Olshen, R.A.; Stone, R.A. Classification and Regression Trees (CART); Wadsworth International: Belmont, CA, USA, 1984; Available online: https://www.routledge.com/Classification-and-Regression-Trees/Breiman-Friedman-Stone-Olshen/p/book/9780412048418 (accessed on 27 May 2021).
  41. Breiman, L. Random Forests. Mach. Learn. 2001, 45, 5–32. [Google Scholar] [CrossRef] [Green Version]
  42. Nascimento, R.F.F.; Alcântara, E.H.; Kampel, M.; Stech, J.L.; Moraes, E.M.L.; Fonseca, L.M.G. O Algoritmo Support Vector Machines (SVM): Avaliação Da Separação Ótima de Classes Em Imagens CCD-CBERS-2. Available online: http://marte.sid.inpe.br/col/dpi.inpe.br/sbsr%4080/2008/10.20.10.59/doc/2079-2086.pdf (accessed on 7 January 2021).
  43. R Development Core Team. R: A Language and Environment for Statistical Computing. Version 3.1.1; R Foundation for Statistical Computing: Vienna, Austria, 2018; ISBN 3-900051-07-0. Available online: http://www.R-Project.Org (accessed on 1 June 2022).
  44. Maino, A.; Alberi, M.; Anceschi, E.; Chiarelli, E.; Cicala, L.; Colonna, T.; De Cesare, M.; Guastaldi, E.; Lopane, N.; Mantovani, F.; et al. Airborne Radiometric Surveys and Machine Learning Algorithms for Revealing Soil Texture. Remote Sens. 2022, 14, 3814. [Google Scholar] [CrossRef]
  45. Niang, M.A.; Nolin, M.C.; Jégo, G.; Perron, I. Digital Mapping of Soil Texture Using RADARSAT-2 Polarimetric Synthetic Aperture Radar Data. Soil Sci. Soc. Am. J. 2014, 78, 673–684. [Google Scholar] [CrossRef]
  46. Bousbih, S.; Zribi, M.; Pelletier, C.; Gorrab, A.; Lili-Chabaane, Z.; Baghdadi, N.; Ben Aissa, N.; Mougenot, B. Soil Texture Estimation Using Radar and Optical Data from Sentinel-1 and Sentinel-2. Remote Sens. 2019, 11, 1520. [Google Scholar] [CrossRef]
  47. Tomasella, J.; Hodnett, M.; Rossato, L. Pedotransfer Functions for the Estimation of Soil Water Retention in Brazilian Soils. Soil Sci. Soc. Am. J. 2000, 64, 327–338. [Google Scholar] [CrossRef]
  48. Villela, A.L.O. Mapeamento Digital de Solos da Formação Solimões Sob Floresta Tropical Amazônica. Ph.D. Thesis, Agronomia-Ciência do Solo, Universidade Federal Rural do Rio de Janeiro, Rio de Janeiro, Brazil, 2013. [Google Scholar]
  49. Adhikari, K.; Kheir, R.B.; Greve, M.B.; Bøcher, P.K.; Malone, B.P.; Minasny, B.; McBratney, A.B.; Greve, M.H. High-Resolution 3-D Mapping of Soil Texture in Denmark. Soil Sci. Soc. Am. J. 2013, 77, 860–876. [Google Scholar] [CrossRef]
  50. Santos, J.R.; Freitas, C.C.; Araujo, L.S.; Dutra, L.V.; Mura, J.C.; Gama, F.F.; Soler, L.S.; Sant’Anna, S.J.S. Airborne P-Band SAR Applied to the Aboveground Biomass Studies in the Brazilian Tropical Rainforest. Remote Sens. Environ. 2003, 87, 482–493. [Google Scholar] [CrossRef]
  51. Tiwari, R.; Singh, R.K.; Chauhan, D.S.; Singh, O.P.; Prakash, R.; Singh, D. Microwave Scattering for Soil Texture at X-Band and Its Retrieval Using Genetic Algorithm. Adv. Remote Sens. 2014, 3, 120–127. [Google Scholar] [CrossRef] [Green Version]
  52. Prakash, R.; Singh, D.; Pathak, N.P. Microwave Specular Scattering Response of Soil Texture at X-Band. Adv. Space Res. 2009, 44, 801–814. [Google Scholar] [CrossRef]
Figure 1. (A) Location of the study area in Central Amazon, Brazil; (B) Total area (TA) sampling, showing the 75% training and 25% validation random samples; and (C) Reference area (RA) sampling, with the 75% training samples concentrated at the Geólogo Pedro de Moura Support Base, and the 25% validation samples lying outside the RA.
Figure 1. (A) Location of the study area in Central Amazon, Brazil; (B) Total area (TA) sampling, showing the 75% training and 25% validation random samples; and (C) Reference area (RA) sampling, with the 75% training samples concentrated at the Geólogo Pedro de Moura Support Base, and the 25% validation samples lying outside the RA.
Remotesensing 14 05711 g001
Figure 2. Flowchart of the methodology used for mapping soil surface (Surf) and subsurface (Sub) clay, silt, and sand contents. T—training; V—validation; RT—regression tree; RF—random forest; SVM—support vector machine; R2—coefficient of determination.
Figure 2. Flowchart of the methodology used for mapping soil surface (Surf) and subsurface (Sub) clay, silt, and sand contents. T—training; V—validation; RT—regression tree; RF—random forest; SVM—support vector machine; R2—coefficient of determination.
Remotesensing 14 05711 g002
Figure 3. Gower index by covariate (blue bars) and general Gower index (red bar and corresponding value) among the reference area and Juruá, Araracanga and Urucu (AC, respectively). The covariates that contributed with the greatest dissimilarity (GDI > 0.25) are highlighted in yellow. DEM—Digital Elevation Model; CI—Convergence Index; TWI—Topographic Wetness Index; RSP—Relative Slope Position; CND—Channel Network Distance; CNBL—Channel Network Base Level; LFf—LS-factor; MRVBF—Multiresolution Index of Valley Bot-tom Flatness; MRRFT—Multiresolution Index of the Ridge Top Flatness; CXI—Convexity Index; ASP—Aspect; LF—Landforms; ProfC—Profile Curvature; PlanC—Plan curvature; VD—Valley Depth; SH—Slope Height; MSP—Mid Slope Position; S—Slope Gradient; MR—Melton Ruggedness; FC—Flow Accumulation; P-band.
Figure 3. Gower index by covariate (blue bars) and general Gower index (red bar and corresponding value) among the reference area and Juruá, Araracanga and Urucu (AC, respectively). The covariates that contributed with the greatest dissimilarity (GDI > 0.25) are highlighted in yellow. DEM—Digital Elevation Model; CI—Convergence Index; TWI—Topographic Wetness Index; RSP—Relative Slope Position; CND—Channel Network Distance; CNBL—Channel Network Base Level; LFf—LS-factor; MRVBF—Multiresolution Index of Valley Bot-tom Flatness; MRRFT—Multiresolution Index of the Ridge Top Flatness; CXI—Convexity Index; ASP—Aspect; LF—Landforms; ProfC—Profile Curvature; PlanC—Plan curvature; VD—Valley Depth; SH—Slope Height; MSP—Mid Slope Position; S—Slope Gradient; MR—Melton Ruggedness; FC—Flow Accumulation; P-band.
Remotesensing 14 05711 g003
Figure 4. Correlation matrices of remote sensing covariates against soil particle size fractions in the training datasets from the reference area (A) and total area (B). ClayA, SiltA, SandA: surface soil particle sizes fractions; ClayB, SiltB, SandB: subsurface soil particle sizes fractions. A strong blue circle has a maximum positive correlation, a strong red circle has a maximum negative correlation. Between these two, the colour tone decreases as the correlation decreases. CI—Convergence Index; TWI—Topographic Wetness Index; RSP—Relative Slope Position; CND—Channel Network Distance; CNBL—Channel Network Base Level; MRVBF—Multiresolution Index of Valley Bottom Flatness; MRRFT—Multiresolution Index of the Ridge Top Flatness; CXI—Convexity Index; ASP—Aspect; LF—Landforms; ProfC—Profile Curvature; PlanC—Plan curvature; SH—Slope Height; MSP—Mid Slope Position; S—Slope Gradient; MR—Melton Ruggedness; FC—Flow Accumulation; P-band.
Figure 4. Correlation matrices of remote sensing covariates against soil particle size fractions in the training datasets from the reference area (A) and total area (B). ClayA, SiltA, SandA: surface soil particle sizes fractions; ClayB, SiltB, SandB: subsurface soil particle sizes fractions. A strong blue circle has a maximum positive correlation, a strong red circle has a maximum negative correlation. Between these two, the colour tone decreases as the correlation decreases. CI—Convergence Index; TWI—Topographic Wetness Index; RSP—Relative Slope Position; CND—Channel Network Distance; CNBL—Channel Network Base Level; MRVBF—Multiresolution Index of Valley Bottom Flatness; MRRFT—Multiresolution Index of the Ridge Top Flatness; CXI—Convexity Index; ASP—Aspect; LF—Landforms; ProfC—Profile Curvature; PlanC—Plan curvature; SH—Slope Height; MSP—Mid Slope Position; S—Slope Gradient; MR—Melton Ruggedness; FC—Flow Accumulation; P-band.
Remotesensing 14 05711 g004
Figure 5. Importance of predictor covariates for the attributes evaluated in the RF model. (A) Sand Surf (B); Sand Sub; (C) Silt Surf; (D) Silt Sub; (E) Clay Surf; (F) Clay Sub. Surf—surface; Sub—subsurface.
Figure 5. Importance of predictor covariates for the attributes evaluated in the RF model. (A) Sand Surf (B); Sand Sub; (C) Silt Surf; (D) Silt Sub; (E) Clay Surf; (F) Clay Sub. Surf—surface; Sub—subsurface.
Remotesensing 14 05711 g005
Figure 6. Soil–relief–vegetation relationships for soil sand, silt, and clay contents in the study area. Green arrow—positive correlation with the covariate; red arrow—negative correlation with the covariate. MU—mapping unit; Fac—Flooded Plain Open Tropical Forest; FDA—Dense Highland Tropical Forest; Fdb—Planalto Open Tropical Forest; APf—River plains; C11—Well-drained flat top areas; T21—Tabular Interfluves; EP2—Biplain-plain surfaces; H.S.—Holocene Sediments; P.S.—Pleistocene Sediments. (Source: modified from ref. [7]).
Figure 6. Soil–relief–vegetation relationships for soil sand, silt, and clay contents in the study area. Green arrow—positive correlation with the covariate; red arrow—negative correlation with the covariate. MU—mapping unit; Fac—Flooded Plain Open Tropical Forest; FDA—Dense Highland Tropical Forest; Fdb—Planalto Open Tropical Forest; APf—River plains; C11—Well-drained flat top areas; T21—Tabular Interfluves; EP2—Biplain-plain surfaces; H.S.—Holocene Sediments; P.S.—Pleistocene Sediments. (Source: modified from ref. [7]).
Remotesensing 14 05711 g006
Figure 7. Accuracy with and without radar P-band for surface sand, silt, and clay prediction with the best model and approach.
Figure 7. Accuracy with and without radar P-band for surface sand, silt, and clay prediction with the best model and approach.
Remotesensing 14 05711 g007
Figure 8. Accuracy with and without radar P-band for subsurface sand, silt, and clay prediction with the best model and approach.
Figure 8. Accuracy with and without radar P-band for subsurface sand, silt, and clay prediction with the best model and approach.
Remotesensing 14 05711 g008
Figure 9. Map of the sand content at the surface layer. (Map generated using Reference Area approach, Random Forest, and Previous Covariate Selection), (A) Urucu block, (B) Araracanga block, (C) Jurua block.
Figure 9. Map of the sand content at the surface layer. (Map generated using Reference Area approach, Random Forest, and Previous Covariate Selection), (A) Urucu block, (B) Araracanga block, (C) Jurua block.
Remotesensing 14 05711 g009
Figure 10. Map of the sand content at the subsurface layer. (Map generated using the Reference Area sampling design, Random Forest, and Wrapper Method), (A) Urucu block, (B) Araracanga block, (C) Jurua block.
Figure 10. Map of the sand content at the subsurface layer. (Map generated using the Reference Area sampling design, Random Forest, and Wrapper Method), (A) Urucu block, (B) Araracanga block, (C) Jurua block.
Remotesensing 14 05711 g010
Figure 11. Map of the silt content at the surface layer. (Map generated using the Reference Area sampling design, Random Forest, and Previous Covariate Selection), (A) Urucu block, (B) Araracanga block, (C) Jurua block.
Figure 11. Map of the silt content at the surface layer. (Map generated using the Reference Area sampling design, Random Forest, and Previous Covariate Selection), (A) Urucu block, (B) Araracanga block, (C) Jurua block.
Remotesensing 14 05711 g011
Figure 12. Map of the silt content at the subsurface layer. (Map generated using the Reference Area sampling design, Random Forest, and Previous Covariate Selection), (A) Urucu block, (B) Araracanga block, (C) Jurua block.
Figure 12. Map of the silt content at the subsurface layer. (Map generated using the Reference Area sampling design, Random Forest, and Previous Covariate Selection), (A) Urucu block, (B) Araracanga block, (C) Jurua block.
Remotesensing 14 05711 g012
Figure 13. Map of the clay content at the surface layer. (Map generated using the Total Area sampling design, Random Forest, and Previous Covariate Selection).
Figure 13. Map of the clay content at the surface layer. (Map generated using the Total Area sampling design, Random Forest, and Previous Covariate Selection).
Remotesensing 14 05711 g013
Figure 14. Map of the clay content at the subsurface layer. (Map generated using the Total Area sampling design, Random Forest, and Previous Covariate Selection).
Figure 14. Map of the clay content at the subsurface layer. (Map generated using the Total Area sampling design, Random Forest, and Previous Covariate Selection).
Remotesensing 14 05711 g014
Table 1. Number (n) and percent of soil taxonomic classes in the 151 field observations.
Table 1. Number (n) and percent of soil taxonomic classes in the 151 field observations.
SiBCS aSoil Taxonomy bWRB bnPercent (%)
Argissolo AmareloUltisolsAcrisols; Lixisols4127.15
Argissolo VermelhoUtisols (Typic Rhodustults)Acrisols; Lixisols21.32
Argissolo Vermelho AmareloUltisolsAcrisols; Lixisols2919.20
Argissolo AcizentadoUltisol (Hapludult)Haplic Lixisol31.98
Cambissolo HáplicoInceptisolsCambisols4932.45
Cambissolo FlúvicoEntisols (Fluvents)Fluvisols21.32
Espodossolos HumilúvicosSpodosols (Alorthods)Podzols10.66
Espodossolos Ferri-HumilúvicosSpodosols (Orthods)Podzols42.65
Neossolo QuartzarênicoEntisols (Quartzipsamments)Arenosols10.66
Neossolos FlúvicosEntisols (Fluvents)Fluvisols21.32
Planossolo HáplicoUltisols (Albaquults)Planosols21.32
Gleissolos HáplicosEntisols (Aquents)Gleysols; Stagnosols149.27
Gleissolos MelânicosEntisols (Fluvaquentic Humaquepts)Umbric Gleysols10.66
Total 151100
a Brazilian Soil Classification System [3]. b Partial equivalence of the soil classes to WRB [33] and Soil Taxonomy [34].
Table 2. Hyperparameters of machine learning algorithms used in this study.
Table 2. Hyperparameters of machine learning algorithms used in this study.
AlgorithmsHyperparametersDefinitionTuning
RTcpA non-negative number for complexity parameter.0.001–0.01
methodANOVAanova
RFmtrynumber of variables used to produce each tree1–10
ntreethe number of trees (default: 500)100–1000
nodesizethe minimum number of data points in each terminal node5
SVMKernel typethe kernel functionpolynomial
typesvm can be used as a classification machine, as a regression machine, or for novelty detection. Depending on whether y is a factor or not, the default setting for type is C-classification or eps-regression, respectively, but may be overwritten by setting an explicit value.‘nu-regression’ or ‘eps-regression’
degreeparameter needed for kernel of type polynomial (default: 3)2–3
costThe cost of predicting a sample within or on the wrong side of the margin.0–10
gammaparameter needed for all kernels except linear (default: 1/(data dimension))1
coef0parameter needed for kernels of type polynomial and sigmoid (default: 0)0
tolerancetolerance of termination criterion (default: 0.001)0.001
RT: regression tree; RF: random forest; SVM: support vector machine.
Table 3. Descriptive statistics of soil texture.
Table 3. Descriptive statistics of soil texture.
VariablesDatasetnMinMaxMeanMedianSDSkKCV (%)
Sand Surf
(g kg−1)
W151809184584371560.36−0.1134
T(RA)1141829184684501540.48−0.0732
V(RA)37807934284091620.11−0.6337
VU212257214254011440.46−0.99-
VA1180793507549176−0.890.77-
VJ515136026727375−0.36−1.38-
T(TA)114808834514351500.21−0.3633
V(TA)372089184814601730.59−0.1935
Sand Sub
(g kg−1)
W151448553533141600.50−0.1645
T(RA)114818553513071550.650.2444
V(RA)37446953573381780.16−1.0949
VU21866743423141690.41−1.00-
VA1144695460493172−0.970.51-
VJ59927919220164−0.12−1.45-
T(TA)114446953373081450.24−0.7543
V(TA)371028554023811930.54−0.6548
Silt Surf
(g kg−1)
W151267923893751450.16−0.2737
T(RA)114266873643511310.03−0.1236
V(RA)37155792466481160−0.11−0.9434
VU21155688476481142−0.42−0.59-
VA112025343543211220.19−1.70-
VJ5597792668643780.56−1.59-
T(TA)114587923983781390.21−0.4035
V(TA)37266963643501600.17−0.3244
Silt Sub
(g kg−1)
W151846003393401050.05−0.2131
T(RA)11484600332328101−0.040.0130
V(RA)371685703613491150.14−1.0532
VU211915703593431130.39−0.87-
VA111684863093031040.19−1.42-
VJ538855148047961−0.30−1.61-
T(TA)114846003493491000.07−0.2129
V(TA)371125823093061160.23−0.4337
Clay Surf
(g kg−1)
W1514500152140860.871.1256
T(RA)11434500169155820.791.0848
V(RA)3744239978771.995.8278
VU2162039886510.23−0.65-
VA114423118731211.340.83-
VJ5271306457400.66−1.37-
T(TA)1144500152139900.871.1059
V(TA)3739351154142740.810.3348
Clay Sub
(g kg−1)
W15113573308326111−0.28−0.2736
T(RA)11413530314330108−0.600.0034
V(RA)37705732882671200.52−0.4342
VU21705732982881310.36−0.76-
VA111505322502001171.120.20-
VJ5259410327340600.13−1.86-
T(TA)11470573314327105−0.09−0.5733
V(TA)3713530289317127−0.49−0.4644
Surf: surface; Sub: subsurface; W: whole dataset; T: training dataset; V: validation dataset; RA: reference area approach; TA: total area approach; VU: Urucu block data set; VA: Araracanga block data set; VJ: Jurua block data set n: number of observations; Min: minimum; Max: maximum; SD: standard deviation; Sk: skewness; K: kurtosis; CV: coefficient of variation.
Table 4. Descriptive statistics of the covariates in the study area by blocks.
Table 4. Descriptive statistics of the covariates in the study area by blocks.
Reference Area (199,167 Pixels)Urucu (11,209,198 Pixels)
Covariates (Unity)MeanMedianSDMinMaxMeanMedianSDMinMax
CI (d)0.030.5916.80−94.5196.07−0.00020.5416.41−98.0898.91
TWI (d)7.667.561.064.6112.308.077.981.234.3312.54
RSP (0–1)0.480.510.30010.440.450.3001
CND (m)6.406.154.01025.395.414.883.95029.64
CNBL (m)61.7261.165.9546.5679.5963.4764.077.1623.0383.16
MRVBF (d)5.739.384.5209.986.699.824.3309.98
MRRFT (d)2.841.972.6707.934.024.763.0907.99
CXI (d)51.3452.417.630.1569.1950.2951.858.89073.19
ASP (°)177.10175.22106.810360173.78171.04107.030360
LF (d)5.325.002.411.0010.005.185.002.111.0010.00
ProfC (m−1)−0−00−0.0090.01−000−0.0130.011
PlanC (m−1)0.03.400.0−0.0070.01000−0.0100.013
SH (m)4.083.551.851.4718.943.843.361.791.1325.51
MSP (%)0.270.250.170.000.820.250.230.160.000.85
S (%)6.235.154.870.0048.865.163.704.770.0067.20
MR (d)0.250.160.290.002.490.210.100.270.002.95
FC (d)2451299630904008120723471449295640014170
P-band (σ°)0.430.430.0700.990.440.440.0600.90
Araracanga (9,364,993 Pixels)Juruá (11,730,902 Pixels)
Covariates (Unity)MeanMedianSDMinMaxMeanMedianSDMinMax
CI (d)00.4916.45−98.7899.010.000.7818.10−99.2199.40
TWI (d)7.927.721.414.3612.377.587.381.283.8612.01
RSP (0–1)0.410.410.31010.350.320.2901
CND (m)6.015.324.83033.924.453.424.08040.50
CNBL (m)63.9365.288.8534.1685.9776.0377.628.4049.8895.63
MRVBF (d)4.964.774.1309.963.703.892.8209.65
MRRFT (d)3.372.673.1509.736.539.364.1909.98
CXI (d)48.3250.9211.07073.4039.5841.138.27063.48
ASP (°)171.04168.26109.060360168.08166.38109.740360
LF (d)5.265.002.321.0010.005.325.002.031.0010.00
ProfC (m−1)−0.00.00.0−0.0110.012−0.0−0.00−0.0140.016
PlanC (m−1)0.00.00−0.0120.0110.00.00−0.0130.018
SH (m)4.183.592.111.1627.333.623.121.731.1432
MSP (%)0.310.290.2000.880.220.180.1600.89
S (%)5.814.255.34050.215.394.025.24076.92
MR (d)0.240.110.3203.010.180.000.2704.23
FC (d)233214212993400133041609105917354006948
P-band (σ°)0.450.450.1100.930.430.430.1000.94
RA: Referencea area; U: block Urucu; A: block Araracanga; J: block Jurua; n: number of observations; Referencea area: (n = 199,167); Urucu: (n = 11,209,198); n Araracanga: (n = 9,364,993); n Jurua: (n = 11,730,902); Min: minimum; Max: maximum; SD: standard deviation. d: dimensionless. CI—Convergence Index; TWI—Topographic Wetness Index; RSP—Relative Slope Position; CND—Channel Network Distance; CNBL—Channel Network Base Level;; MRVBF—Multiresolution Index of Valley Bottom Flatness; MRRFT—Multiresolution Index of the Ridge Top Flatness; CXI—Convexity Index; ASP—Aspect; LF—Landforms; ProfC—Profile Curvature; PlanC—Plan curvature; SH—Slope Height; MSP—Mid Slope Position; S—Slope Gradient; MR—Melton Ruggedness; FC—Flow Accumulation; P-band.
Table 5. Accuracy assessment soil surface and subsurface sand content predictions using the Reference Area (RA) sampling design.
Table 5. Accuracy assessment soil surface and subsurface sand content predictions using the Reference Area (RA) sampling design.
RTRFSVM
AtributtesDataR2RMSEMAER2RMSEMAER2RMSEMAE
Sand Surf
PCS
T0.34124960.9367530.4711391
VU0.091441170.241311130.07141125
VUA0.031651290.191401160.01162140
V0.011761350.241441200.08173149
VUJ0.031661280.311381190.12162141
Sand Surf
WM
T0.36122960.9467530.5710686
VU0.211321030.201321140.04144129
VUA0.061641280.181411160.20140118
V0.031751320.191481230.19145123
VUJ0.091571130.221431250.08151134
Sand Sub
PCS
T0.45114900.9263500.4711395
VU0.091611370.241471260.20148126
VUA0.011811520.151631400.18166141
V0.001901550.111651430.24168139
VUJ0.021791450.171541320.21155127
Sand Sub
WM
T0.48111870.9264510.5710586
VU0.141591330.361371130.13154128
VUA0.051811550.251581340.15173146
V0.001941630.161621380.17167141
VUJ0.031831490.251471230.17148124
PCS—Previous Covariate Selection; WM—Wrapper Method; T: Training dataset; V: validation data set; VU: Urucu block validation dataset; VUA: Urucu/Araracanga block validation dataset; V: Urucu/Araracanga/Jurua block validation dataset; VUJ: Urucu/Jurua block validation dataset.
Table 6. Accuracy assessments of soil surface and subsurface silt content predictions using the Reference Area (RA).
Table 6. Accuracy assessments of soil surface and subsurface silt content predictions using the Reference Area (RA).
RTRFSVM
AtributtesDataR2RMSEMAER2RMSEMAER2RMSEMAE
Silt Surf
PCS
T0.4993720.9156430.509371
VU0.191631380.581301120.33130107
VUA0.071541290.37120990.17175123
V0.071751440.361411140.28185133
VUJ0.181891580.521551310.38156124
Silt Surf
WM
T0.4695730.925542.0.588765
VU0.261631440.461391200.24143122
VUA0.061571360.261281060.13149119
V0.081741490.261491220.22159128
VUJ0.261861610.421641400.26158134
Silt Sub
PCS
T0.4773580.9143320.397961
VU0.3690720.5189710.389177
VUA0.3886720.4188730.3311191
V0.2699800.4689740.39131101
VUJ0.22106830.5690730.3912693
Silt Sub
WM
T0.4972570.9243320.537256
VU0.3589720.4293740.428467
VUA0.3389730.3193760.399176
V0.22102810.3794780.3911589
VUJ0.21106830.5095780.3712088
PCS—Previous Covariate Selection; WM—Wrapper Method; T: Training dataset; V: validation data set; VU: Urucu block validation dataset; VUA: Urucu/Araracanga block validation dataset; V: Urucu/Araracanga/Jurua block validation dataset; VUJ: Urucu/Jurua block validation dataset.
Table 7. Accuracy assessment of soil surface and subsurface clay content predictions using the Reference Area (RA) sampling design.
Table 7. Accuracy assessment of soil surface and subsurface clay content predictions using the Reference Area (RA) sampling design.
RTRFSVM
AtributtesDATAR2RMSEMAER2RMSEMAER2RMSEMAE
Clay Surf
PCS
T0.5355410.9131230.476145
VU0.0973590.2471590.216753
VUA0.0490700.0292720.0811573
V0.0390730.0292730.0411173
VUJ0.0678650.1976640.176957
Clay Surf
WM
T0.5454400.9231230.565641
VU0.0874590.1871590.276550
VUA0.0489700.0291710.178261
V0.0390730.0291720.109672
VUJ0.0578660.1575630.159168
Clay Sub
PCS
T0.6167530.9139300.587052
VU0.16119900.20114860.1412093
VUA0.021361010.08122950.1711795
V0.02130930.07116890.1311392
VUJ0.15113810.18107800.1211490
Clay Sub
WM
T0.6265520.9238290.656549
VU0.021381030.18115870.0712899
VUA0.001461110.08120930.1411893
V0.001411040.07114880.03152116
VUJ0.03131950.17108810.02170131
PCS—Previous Covariate Selection; WM—Wrapper Method; T: Training dataset; V: validation data set; VU: Urucu block validation dataset; VUA: Urucu/Araracanga block validation dataset; V: Urucu/Araracanga/Jurua block validation dataset; VUJ: Urucu/Jurua block validation dataset.
Table 8. Accuracy assessment of soil surface and subsurface sand, silt, and clay content predictions using the Total Area (TA) sampling design.
Table 8. Accuracy assessment of soil surface and subsurface sand, silt, and clay content predictions using the Total Area (TA) sampling design.
RTRFSVM
AtributtesDataR2RMSEMAER2RMSEMAER2RMSEMAE
Sand Surf PCST1140.51104790.9362490.5210584
V370.001981520.111611240.15163127
Sand Surf WMT1140.51104790.9464500.777344
V370.001981520.131591240.03209158
Sand Sub PCST1140.5497800.9358470.4011394
V370.032021480.231741370.21180138
Sand Sub WMT1140.5597800.9559480.816441
V370.032021480.221771450.19207162
Silt Surf PCST1140.5889690.9153410.509878
V370.041821400.141471130.20142108
Silt Surf WMT1140.9254420.9254420.608972
V370.171441110.171441110.14147112
Silt Sub PCST1140.4971570.9139310.427661
V370.06123970.03120980.2910279
Silt Sub WMT1140.5169550.9238300.556954
V370.04126990.06116940.2110784
Clay Surf PCST1140.5658440.9134250.596046
V370.2371580.2365500.157052
Clay Surf
WM
T1140.5857430.9233250.655642
V370.2074620.2165480.128062
Clay Sub PCST1140.5470550.9338300.577056
V370.19117940.31107810.2911492
Clay Sub WMT1140.5173580.9339300.616853
V370.21116930.30107820.2612294
PCS—Previous Covariate Selection; WM—Wrapper Method; T: Training dataset; V: validation data set; RT: regression tree; RF: random forest; SVM: support vector machine.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Ferreira, A.C.d.S.; Ceddia, M.B.; Costa, E.M.; Pinheiro, É.F.M.; Nascimento, M.M.d.; Vasques, G.M. Use of Airborne Radar Images and Machine Learning Algorithms to Map Soil Clay, Silt, and Sand Contents in Remote Areas under the Amazon Rainforest. Remote Sens. 2022, 14, 5711. https://0-doi-org.brum.beds.ac.uk/10.3390/rs14225711

AMA Style

Ferreira ACdS, Ceddia MB, Costa EM, Pinheiro ÉFM, Nascimento MMd, Vasques GM. Use of Airborne Radar Images and Machine Learning Algorithms to Map Soil Clay, Silt, and Sand Contents in Remote Areas under the Amazon Rainforest. Remote Sensing. 2022; 14(22):5711. https://0-doi-org.brum.beds.ac.uk/10.3390/rs14225711

Chicago/Turabian Style

Ferreira, Ana Carolina de S., Marcos B. Ceddia, Elias M. Costa, Érika F. M. Pinheiro, Mariana Melo do Nascimento, and Gustavo M. Vasques. 2022. "Use of Airborne Radar Images and Machine Learning Algorithms to Map Soil Clay, Silt, and Sand Contents in Remote Areas under the Amazon Rainforest" Remote Sensing 14, no. 22: 5711. https://0-doi-org.brum.beds.ac.uk/10.3390/rs14225711

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop