Thank you for visiting nature.com. You are using a browser version with limited support for CSS. To obtain the best experience, we recommend you use a more up to date browser (or turn off compatibility mode in Internet Explorer). In the meantime, to ensure continued support, we are displaying the site without styles and JavaScript.
Advertisement
Scientific Reports volume 13, Article number: 10897 (2023)
2478
3
Metrics details
The pressure for food production has expanded agriculture frontiers worldwide, posing a threat to water resources. For instance, placing crop systems over hydromorphic soils (HS), have a direct impact on groundwater and influence the recharge of riverine ecosystems. Environmental regulations improved over the past decades, but it is difficult to detect and protect these soils. To overcome this issue, we applied a temporal remote sensing strategy to generate a synthetic soil image (SYSI) associated with random forest (RF) to map HS in an 735,953.8 km2 area in Brazil. HS presented different spectral patterns from other soils, allowing the detection by satellite sensors. Slope and SYSI contributed the most for the prediction model using RF with cross validation (accuracy of 0.92). The assessments showed that 14.5% of the study area represented HS, mostly located inside agricultural areas. Soybean and pasture areas had up to 14.9% while sugar cane had just 3%. Here we present an advanced remote sensing technique that may improve the identification of HS under agriculture and assist public policies for their conservation.
Over the past decades multiple environmental challenges were addressed and worldwide initiatives have been established to promote sustainable development practices1. The world’s growing population and the increase demand for food and water, while minimizing the impact on climate raised awareness on the effort to achieve Food, Water and Energy Security, Climate Change Abatement, Biodiversity Protection and Ecosystem Service Delivery2,3. The soil has an important role on the achievement of such goals, but so far has been poorly applied in models to investigate these global challenges4. The lack of soil knowledge and the advance of its degradation caused by agriculture pose a global threat as the population is estimated to be 9 billion by the middle of the twenty-first century5,6.
Brazil has 41% of its area dedicated to agriculture (351 million ha)7, which increased 71 million ha in the last 33 years to cattle ranching and agriculture activities8. Naturally, fragile environments are pressured by this expansion, since agriculture activity releases chemicals considered toxic to flora, fauna, and human health9,10. In this regard, HS stand out as a fragile ecosystem responsible for hydrological and biogeochemical cycles with fauna and flora11,12. These soils are connected with the water table and represent a supplier for water recharge nutrients and sediments for riverine ecosystems13.
Saturated and waterlogged soils occupy around 6% of the Earth’s surface (2,1 million km2)14,15 and is the ecosystem with the highest damage rates16. Multiple studies presented ground water contamination through agricultural activities involving fertilizers, pesticides, and other chemical agents17,18, implying the fragility of such ecosystems and the impact of unregulated anthropic activities19. According to the Brazilian Forestry Code, the margin areas of a river, lake, and water source must be preserved with the natural vegetation20. However, the areas to be preserved start counting from the riverbed at the drought season, which end up excluding potential wetland soils that will only be affected during the flood season. These soils will naturally be included in agricultural sites at risk of getting contaminated, since they are difficult to map.
Identification and mapping of HS is vital to preserve these ecosystems from degradation regarding agriculture and other anthropic activity expansion. Multiple works combined digital geospatial information (i.e. satellite data) with machine learning algorithms to explore the relationship between soil and hydrology21. Some works focused on how the soil and landscape attributes affect channels morphometry22,23,24,25,26, while others explored the behavior within the soil27,28,29,30. Thompson 199731 presented an early effort to quantitatively map HS through a soil color index using field sampling and terrain attributes derived from a digital elevation model. As remote sensing tools advanced, the mapping of wetlands32,33,34,35 and HS36,37 were possible due to the soil moisture effect on the reflectance intensity38.
To this end, agriculture expansion is bringing environmental pressure over HS and wetlands, demanding more detailed soil survey and mapping approaches. To our knowledge, there is no use of a remote sensing tool to assess HS located inside cultivated and not cultivated areas, enabling a better regulation of agriculture activity. Here we present a new remote sensing technique that allows the identification of HS at 30 m of spatial resolution for large areas. The method described as GEOS339 uses a time-series of Landsat images to extract pixels with bare soil and aggregate them into a single synthetic soil image (SYSI). We used the Landsat 4 Thematic Mapper (TM) (1982–1993), Landsat 5 TM (1984–2012), Landsat 7 Enhanced Thematic Mapper Plus (ETM +) (1999–2018), and the Landsat 8 Operational Land Manager (OLI) (2013–2018). The technique combined more than 30 years of data in order to extract all the areas with bare soil, enabling a better regulation of agriculture activity and public policies towards environmental protection.
We obtained 4954 georeferenced soil field samples of hydromorphic and not HS from the Brazilian Soil Spectral Library (BSSL)40 in an 863,577.9 km2 area in Brazil (Supplementary Fig. S1). This dataset contains samples with laboratory spectral analysis (350–2500 nm) and respective spectra acquired from bare soil satellite spectra from a synthetic soil image (SYSI)39. We also obtained 1579 additional synthetic observations of HS by combining field visits and qualitative analysis of SYSI, which served as reference for the digital sampling of locations with hydromorphic conditions (up to 40 cm depth). We selected field samples with laboratory spectra and classified as hydromorphic (Gleysols, Planosols, and Fluvisols) to evaluate their spectral behavior and to compare them with the SYSI spectra (Fig. 1a,b) and, to relate the SYSI spectra with agricultural areas covered with vegetation and uncovered hydromorphic features under the vegetation (Fig. 1b–d, respectively) (“Methods”).
Expression of SYSI indicating hydromorphic soils across the study area, inside agricultural. (a) Study area with soil locations classified as hydromorphic and not hydromorphic. (b) Synthetic soil image (SYSI) for the study area. (c) Agricultural area covered with vegetation. (d) Uncovered hydromorphic features under the vegetation. Map created with ESRI ArcGIS 10.4.
The main characteristic observed in the laboratory spectra was the absence of a concave feature in the region of 900 nm, typical of the presence of iron oxyhydroxides41 (Fig. 1a). There is also an attenuation of features caused by the presence of organic matter in the soil between 350 and 1300 nm, with some of the curves showing a concave-rectilinear pattern between these bands (Fig. 2a). In addition, some samples presented a small convex feature between 350 and 450 nm. Therefore, the main spectral signature characteristics of HS were observed in the range from 350 to 1350 nm (Fig. 2a).
Laboratory and satellite topsoil spectra from hydromorphic and not hydromorphic soils. (a) Visible, near infrared, and short-wave infrared (350–2500 nm) laboratory topsoil spectra from A horizon of hydromorphic soils. (b) The spectral curves were convolved into the Landsat/TM bands to analyze the satellite spectral behavior. (c) Laboratory and satellite topsoil spectra from hydromorphic and not hydromorphic locations found in SYSI. (d) SYSI spectra collected at the same field locations from 1c. (e) Second derivative of the spectra to highlight the lower features of iron oxides at the hydromorphic locations.
The spectral signatures presented in Fig. 1a were convolved for the Landsat bands in order to compare the spectra from the hydromorphic field samples with the SYSI spectral data (“Methods”). There was a concave-rectilinear pattern between Landsat band 1 to band 4 (Fig. 1b), indicating that part of the characteristics presented in the spectral signatures of laboratory data are reflected in the SYSI signatures, as, for instance, the attenuation caused by organic matter and absence of iron oxides, which are normal conditions for HS.
Four soil field samples were selected to illustrate the spectral differentiation between hydromorphic and not hydromorphic soils (Fig. 2c). Two samples were located in a toposequence with clayey soils and the other two sites were located in a toposequence with sandy soils (“Methods”). The HS had greater reflectance than not hydromorphic ones and did not present the typical concave shape around 900 nm regarding presence of Fe oxides (Fig. 2c). We applied the second derivative of the Kubelka–Munk function to the soil spectra42, which highlighted the absence of the typical hematite amplitude located between 520 and 580 nm (Fig. 2e). The result showed a slight peak around 560 nm, but not an absorption at 525 nm for HS samples and a higher reflectance between 350 and 450 nm for sandy soils (Fig. 2e).
When we observed the SYSI spectral signatures at the same four locations, we verified a higher reflectance for HS, and a slightly more concave shape between band 1 and band 4 (Fig. 2d). This was caused by a low presence of Fe oxides and higher contents of quartz (sand and silt fractions) in the soil surface, a result from permanent or periodic saturation of the soil by water43. The reflectance detected by the satellite only retrieves information about the surface, although it can be related to subsurface characteristics44. The surface of Gleysols (Hydromorphic) are not necessarily wet, contrary to the subsurface where water saturation promotes anoxic conditions. Nonetheless, the surface mineral and textural properties can be affected by the water saturation from below (groundwater) or from above (rain or irrigation water) and removal or reduction of Fe3+ to Fe2+45.
Based on the absence of Fe oxides features in the laboratory and SYSI spectra at the hydromorphic sites, we identified specific areas as hydromorphic in the synthetic image, which enabled its use for the classification of hydromorphic areas. The hydromorphic sites were often near drainage network systems, which are zones normally affected by groundwater level fluctuations and that concentrate the drained water from upper positions of the watershed (Supplementary Fig. S3). This observation suggested the influence of soil moisture in the reflectance response, causing differences in color intensity similar to the features found by46.
We combined the soil dataset with SYSI and a set of terrain attributes calculated with the Terrain Analysis in Google Earth Engine package47 (Supplementary Table S1). The resulting dataset was used to fit a random forest (RF) model in order to classify the areas as hydromorphic across the study area (Methods). The optimal model used RF and cross-validation as the resampling method, reaching an overall accuracy of 0.92 and Kappa coefficient of 0.77 for the binary classification of pixels as hydromorphic and not hydromorphic. We also evaluated the model with a confusion matrix, reaching a producer accuracy of 87% for the hydromorphic and 94% for the not hydromorphic class (Supplementary Table S2).
The model was able to correctly classify 6086 out of 6533 soil observations displaced across the study area. The boxplot analysis indicated a significant difference in slope, elevation, and horizontal curvature regarding hydromorphic and not hydromorphic classes (Supplementary Fig. S2a), which also manifested in the prediction performance with 98, 38, and 29% of importance for the model prediction (Supplementary Fig. S2b). SySI had a similar prediction importance for bands 2, 3, 4, 6, and 7, while the boxplot analysis showed significant differences for all satellite bands between hydromorphic and not hydromorphic soils (Supplementary Fig. S2a). SYSI band 1 (Blue) had the second highest contribution for the model with 69%, highlighting the influence of soil moisture in the reflectance response.
We selected four locations from the predicted map (Fig. 3a) and displayed them in true color composition, RGB 543 and 321 to observe and analyze the areas classified as hydromorphic (Fig. 2b). It was possible to identify some bright and darker features in the images, showing differences in the reflectance factor and the occurrence of HS at different landscape positions (Fig. 2b).
Results of the modelling and the relief pattern of areas classified as hydromorphic. (a) Predicted map of hydromorphic soils for the study area. (b) Indication of SYSI as a tool to identify hidden hydromorphic soils. (c) Toposequence extracted from area (i) indicating the relief positions of hydromorphic soils, the geology and soil classes. Map created with ESRI ArcGIS 10.4.
The differences in the reflectance factor usually indicates textural variation, Fe oxides presence, higher organic carbon content, soil moisture, and others48. Figure 2bi presented three possible HS, one surrounding a channel’s water source enhancing the possible area of Gleysols. In the same location there is an intermittent channel of second order, which is an overland flow path. Finally, a closed depression with concave landform that functions as an accumulation zone (Fig. 3bi). Figure 3bii represents a footslope, a flat landform next to a thalweg and normally influenced by groundwater level fluctuations. The Landsat 321 RGB composite showed an area with lower reflectance intensity at the footslope, indicating hydromorphism (Fig. 3bii). A flat area at the summit and an overland flow path with the same features of lower reflectance intensity were observed in Fig. 3biii, iv. These areas are displaced across the landscape, retaining water within the soil due to landforms or a soil characteristic that hampers water infiltration.
We plotted the relief profile in a toposequence to evaluate how the relief contributes for the formation of a HS (Fig. 3c). The HS occurred over two soil types, a Ferralsol and a Gleysol (Fig. 3c). The Ferralsols are weathered soils normally located at flat surfaces, which favor water infiltration and prevent the formation of drainage channels45. The Ferralsol was also located at the hillslope with convex surface, favoring surface runoff and drainage channel formation (Fig. 3c). The Gleysols are normally located at lower relief positions that constantly receives and accumulates sediment and water, favoring redoximorphic activity in the soils49.
Figure 4a shows an agricultural site with a Permanent Protected Area (PPA) surrounding the channels and the water source. However, a prolonged area with HS is under agriculture exploration (Fig. 4b). The HS is hidden under an agricultural area in contact with chemical fertilizers and pesticides, common for management practices. Agricultural practices over these soils can accelerate nutrient loss, affect particle aggregation, distribution and mineralogy of Fe oxides between particle-size fractions, and the interaction with organic matter stabilization50,51. Also, these soils are more fragile to receive pollution from herbicides and pesticides, reaching ground water. Two Landsat images from the same area in 2014 and 2017 indicate a different location for the water source and an intermittent channel connecting with the actual PPA (Fig. 4c,d). The technique identified these soils and showed that multiple hydromorphic areas are located at agricultural fields, promoting degradation (Table 1).
Example of hydromorphic soils within an agricultural field. (a) The supposed water source is protected with vegetation. (b) SYSI identifies abrupt difference in reflectance indicating hydromorphic conditions. (c,d) Landsat true color in 2014 and 2017 with exposed soil showing accumulation of water in the actual water source. Map created with ESRI ArcGIS 10.4.
We found differences in the hydromorphic areas between federal states and land uses, requiring further investigations for the correct application Brazilian Forestry Code. Only 6% of the studied area belonging to the state of São Paulo were classified as hydromorphic, while Mato Grosso do Sul, Minas Gerais, and Goiás had up to 17.8% (Table 1). Mato Grosso had only 5.8% of its area included in this work, but had 28% classified as hydromorphic. In the world about 26% of the gleysols are agricultural lands52. The state of São Paulo is considered the Brazilian state with the greatest legal framework for environmental protection53,54. In 1994, it was the first state to regulate the use of floodplain areas. In addition, it has a law that regulates the use and conservation of agricultural soil (“Soil law”) since 1988. In addition to São Paulo, only the states of Paraná (1984), Espírito Santo (2001) and Rio Grande do Sul (2015), have similar regulation55.
The sugar cane areas had 3.1% of hydromorphic areas, an indication of adequate ap-plication of PPA for a major agricultural activity in the country56. On the other hand, pastures across the study area had 14.9% of hydromorphic areas, enhancing the risks of contamination since this land use is normally degraded due to due to low grass productivity and inadequate grazing management56. The soybean areas also had around 14% of areas classified as hydromorphic, followed by forest plantation and temporary crops with 10 and 12%, respectively (Table 1). Soybean is the predominant crop system in the states of Mato Grosso, Mato Grosso do Sul, and Goiás57, which were the states with higher percentage of HS (Table 1). Sugar cane is the predominant agricultural system in São Paulo state, which explains the lower occurrence of HS in this state.
The observed differences in HS regarding federal states and land use suggest a further investigation on the application of the Brazilian Forest Code. The modification of hydromorphic environments increases the emission of gases and changes the soil dynamics58, being important areas for springs and water bodies59, suggesting stronger preservation of these areas. Multiple initiatives discussed the real benefits and limitations of the defined areas of preservation, discussing whether they should be larger in order to preserve natural resources. However, this technique was able to identify multiple channel networks (intermittent and perennial) and water sources inside agricultural sites (Supplementary Fig. S3). In the current Forest Code, the delimitation of PPAs occurs from the regular riverbank, unlike the previous code (Forest Code of 1965), in which the largest riverbank was considered60. This change may contribute to an increase in the agricultural use of HS areas, especially in floodplain areas. Finally, more remote sensing tools must be applied to monitor the control of natural resources, which should be preserved and well managed in order to avoid future degradation.
We obtained 4715 georeferenced soil field samples of hydromorphic and not hydromorphic soils from the Brazilian Soil Spectral Library (BSSL) in an 735,953.8 km2 area across the Southeast and Midwest regions of Brazil (Supplementary Fig. S1). The study area comprises tropical and subtropical climates classified as Savanna (Aw), Subtropical highlands (Cwb), and Humid subtropical (Cfa, and Cwa) according to the Köppen climate classification. The rainfall varies between 1000 to 2200 mm per year and the mean annual temperature between 18 and 24 °C61.
Most of these soil samples were acquired from traditional soil surveys, which consist of a soil specialist using conventional soil surveying methods to select locations based on pedogeomorphological relationships for profile description and sampling62. We excluded the Gleysols and Planosols observations from the dataset to keep only not hydromorphic soils. We also obtained 1341 additional synthetic observations of HS by combining field visits and visual analysis of a synthetic soil image (SYSI)39, which served as reference for the digital sampling of locations with hydromorphic conditions (up to 40 cm depth). These synthetic samples were based on geographical location and SYSI’s reflectance intensity (indicating soil moisture). Finally, we reached a dataset with two categories, as follows: (i) not hydromorphic soils, composed by Acrisols, Cambisols, Chernozems, Podzols, Ferralsols, Luvisols, Leptosols, Arenosols, Regosols, and Nitisols; (ii) HS, composed by synthetic samples (field and digital).
We used an external soil dataset with proximal spectral information (350–2500 nm) of Gleysols and Planosols (Hydromorphic) from three works40,63,64. We used these data to identify the features related to hydromorphism, such as organic matter presence, absence of Fe oxide features, and the spectra intensity. We analyzed 10 samples collected at 0–20 cm depth to compare the laboratory and satellite spectra. The method for the spectroscopic analysis is described by Demattê et al.65. Thus, we convolved the soil spectra for the Landsat bands42, in order to compare the spectra from the hydromorphic field samples with the SySI spectral signature.
After establishing the spectral signatures of HS, we selected two samples at locations where SySI had the hydromorphic conditions. These samples contained SYSI and laboratory spectral analysis (350–2500 nm) and were used to define if the changes in spectral intensity was in fact related to hydromorphism. We also applied the second derivative of the Kubelka–Munk function, as a way to highlight the changes in the spectra42. After confirming the spectral pattern (350–2500 nm) for HS, we selected 1341 locations with the same spectral features. SYSI showed a change in hue and intensity of colors at the locations with possible hydromorphic conditions. This pattern was normally observed next to drainage channels, at flat surfaces, and close to the water sources.
Soil formation is a result from the interaction of multiple factors such as climate, organisms (biota), relief (landscape processes), parent material (geology), and time66. Thus, we combined the soil dataset with a set of environmental variables related to soil forming factors (Supplementary Table S1) to fit a RF model.
We used two categories of environmental variables, the first is related to bare soil reflectance from satellite images and the second is related to the terrain. We implemented the Geospatial Soil Sensing System (GEOS3)39 to a time-series of Landsat images using the Google Earth Engine (GEE) platform67 in order to generate a 30 m spatial resolution Synthetic Soil Image (SYSI). SYSI has six spectral bands from blue to short-wave infrared regions at 30 m resolution (Supplementary Table S1).
The terrain attributes were calculated using the Terrain Analysis in GEE (TAGEE) package47. The package calculates multiple topographic variables from a Shuttle Radar Topography Mission (SRTM) digital elevation model (DEM) with a spatial resolution of 30 × 30 m. Finally, we overlapped the soil dataset with the TAGEE and SYSI attributes and sampled the matching locations to perform a boxplot analysis. The sampled data were plotted by class (hydromorphic and not hydromorphic) in order to analyze the distribution of values. The graphic was performed using the “ggplot” package in R software.
The RF algorithm was selected to perform the hydromorphic areas mapping, since its relevance for digital soil mapping68,69. RF estimates a user-specified number of decision trees by randomly sampling an existing dataset70. However, at each node construction, a random sample of the dependent variables is used. The resulting decision tree is used to estimate the out-of-bag error rate by predicting the value of the remaining unsampled data and comparing with the known results.
We performed a grid search to select the optimal hyperparameters, which were the maximum depth (150), maximum features (3), minimum samples leaf (1), minimum samples split (10), and number of trees (300). These parameters regulate the number of variables that can be randomly sampled in each split of the trees, the tree depth by setting the minimal number of samples for the terminal nodes, and the number of trees.
In order to calibrate the RF model, we tested three resampling methods coupled with the RF model using the caret package in R software. The first test used k-fold cross-validation (CV) method to fit the prediction models. CV is a resampling method used to fix optimistic results of the predictive effectiveness of regression equations. The method randomly divided the data in k groups, using k − 1 groups to fit a model, and one for validation71,72. The procedure is repeated k times, always leaving one group out of the calibration dataset. Afterwards, the results are summarized with the mean of the model scores. The second method was the bootstrapping, which is a data resampling technique for estimating the statistical parameters of an unknown distribution and a robust method for optimal model selection73. We also tested the out-of-bag resampling method which is a method of measuring the prediction error of random forests74.
The prediction performance of the data was accessed using the three default parameters of caret for classification models, being the number of randomly selected predictors (mtry), accuracy, and kappa coefficient. The mtry regulates the number of variables that can be randomly sampled in each split of the trees, which resulted in 2, 20, and 39. We used 300 trees for stable variable estimates.
As many environmental information were used as covariables to fit the RF model, we analyzed the variables’ importance through the mean decrease Gini index. The analysis helped indicating how was the contribution of the terrain, climate, and remote sensing variables.
After testing the resampling methods and fitting the RF model, we selected the optimal model and used its parameters to predict the hydromorphic and not hydromorphic classes across the study area with the caret package in R software 3.4.0 (https://cran.r-project.org/bin/windows/base/)75,76. The resulting binary map was generated using the caret and raster packages in R software, and exported to ArcGIS 10.4. For further analysis. We calculated a confusion matrix for the optimal model and analyzed the errors of inclusion or commission errors (CE), errors of exclusion or omission errors (OE), user accuracy (UA), producer accuracy (PA), and global accuracy (GA)77.
We also evaluated the areas mapped as hydromorphic according to their relief position (infiltration or surface runoff environment), proximity to the channel network or water source, and according to the pattern registered by SYSI. We selected an area to analyze the distribution of HS across a toposequence, based on the soil-landscape relationship rules and the channel network patterns.
With the digital map of HS, we were able to analyze the current land use situation of these soils. First, we masked the pixels classified as hydromorphic and exported them to a new raster in ESRI ArcGIS 10.4. The new raster was plotted over a land use and landcover map from the MapBiomas project. MapBiomas is a governmental initiative aimed to reconstruct annual land use and land cover information between 1985 and 2017 for Brazil based on random forest applied to Landsat archive using Google Earth Engine8. The dataset is available for download at their repository in GEE and at their website (https://mapbiomas.org/)8.
Finally, we quantified the areas of HS and identified the land uses at the areas. We also computed the total areas of HS for each federal state included in the study area. The result was presented in a table with the total area of HS for each land use class and federal state.
The hydromorphic soils map generated during the current study is available in raster format in the link: https://esalqgeocis.wixsite.com/english/felipe-mello-hydromorphic-soils.
The code to generate the final product of the hydromorphic soils map is available from the corresponding author upon reasonable request for the following reasons: (a) the script/code creates a proxy from which our research group is working on the development of new products that will be published in the near future.
Wu, H. & Clark, H. The sustainable development goals: 17 goals to transform our world. in Furthering the Work of the United Nations 36–54 (United Nations, 2016). https://doi.org/10.18356/69725e5a-en.
Godfray, H. C. J. et al. Food security: The challenge of feeding 9 billion people. Science (1979) 327, 812 LP–818 (2010).
Bouma, J. & McBratney, A. Framing soils as an actor when dealing with wicked environmental problems. Geoderma 200–201, 130–139 (2013).
Article ADS Google Scholar
McBratney, A., Field, D. J. & Koch, A. The dimensions of soil security. Geoderma 213, 203–213 (2014).
Article ADS Google Scholar
Baker, R. M. Soil resilience and sustainable land use. Experimental agriculture vol. 31 (Cambridge University Press, 1994).
Lal, R. & Stewart, B. A. Food security and soil quality (CRC Press, 2010). https://doi.org/10.1201/EBK1439800577.
Teixeira, G. O Censo Agropecuário 2017. Revista NECAT-Revista do Núcleo de Estudos de Economia Catarinense 8, 8–39 (2019).
Google Scholar
Souza, C. M. et al. Reconstructing three decades of land use and land cover changes in Brazilian biomes with Landsat archive and earth engine. Remote Sens. 12. https://doi.org/10.3390/rs12172735 (2020).
Rosolen, V., De-Campos, A. B., Govone, J. S. & Rocha, C. Contamination of wetland soils and floodplain sediments from agricultural activities in the Cerrado Biome (State of Minas Gerais, Brazil). Catena (Amst) 128, 203–210 (2015).
Article CAS Google Scholar
Amendola, D., Mutema, M., Rosolen, V. & Chaplot, V. Soil hydromorphy and soil carbon: A global data analysis. Geoderma 324, 9–17 (2018).
Article ADS CAS Google Scholar
Santana, S. E. & Barroso, G. F. Integrated ecosystem management of river basins and the coastal Zone in Brazil. Water Resour. Manage 28, 4927–4942 (2014).
Article Google Scholar
Lehrback, B. D., Neto, R. R., Barroso, G. F. & Bernardes, M. C. Sources and distribution of sedimentary organic matter in the northwestern portion of Victoria Bay, ES. Braz. J. Aquat. Sci. Technol. 20, 79–92 (2016).
CAS Google Scholar
Pennock, D., Bedard-Haughn, A., Kiss, J. & van der Kamp, G. Application of hydropedology to predictive mapping of wetland soils in the Canadian Prairie Pothole Region. Geoderma 235–236, 199–211 (2014).
Article ADS Google Scholar
Zhang, Z., Zimmermann, N. E., Kaplan, J. O. & Poulter, B. Modeling spatiotemporal dynamics of global wetlands: Comprehensive evaluation of a new sub-grid TOPMODEL parameterization and uncertainties. Biogeosciences 13, 1387–1408 (2016).
Article ADS Google Scholar
Secretariat, R. The Ramsar handbooks for the wise use of wetlands. Preprint at (2010).
Davidson, N. C. & Finlayson, C. M. Earth observation for wetland inventory, assessment and monitoring. Aquat. Conserv. 17, 219–228 (2007).
Article Google Scholar
Zou, Y. et al. Impacts of agricultural and reclamation practices on wetlands in the Amur River Basin Northeastern China. Wetlands 38, 383–389 (2018).
Article Google Scholar
Gebresllassie, H., Gashaw, T. & Mehari, A. Wetland degradation in Ethiopia: Causes, consequences and remedies. J. Environ. Earth Sci. 4, 40–48 (2014).
Google Scholar
Ritter, W. F. Pesticide contamination of ground water in the United States—A review. J. Environ. Sci. Health B 25, 1–29 (1990).
Article CAS PubMed Google Scholar
Piedade, M. T. F. et al. As áreas úmidas no âmbito do Código Florestal brasileiro. Código Florestal ea ciência: o que nossos legisladores ainda precisam saber. Sumários executivos de estudos científicos sobre impactos do projeto de Código Florestal 9–17 (2012).
Lin, H. Hydropedology: Synergistic integration of soil science and hydrology. (Academic Press, 2012).
Easton, Z. M. et al. Re-conceptualizing the soil and water assessment tool (SWAT) model to predict runoff from variable source areas. J. Hydrol. (Amst) 348, 279–291 (2008).
Article ADS Google Scholar
Jung, K., Niemann, J. D. & Huang, X. Under what conditions do parallel river networks occur?. Geomorphology 132, 260–271 (2011).
Article ADS Google Scholar
Filipović, V., Gerke, H. H., Filipović, L. & Sommer, M. Quantifying subsurface lateral flow along sloping horizon boundaries in soil profiles of a hummocky ground moraine. Vadose Zone J. 17, (2018).
Mello, F. A. O. et al. Expert-based maps and highly detailed surface drainage models to support digital soil mapping. Geoderma 384, 114779 (2021).
Article ADS Google Scholar
Mello, F. A. O. et al. Complex hydrological knowledge to support digital soil mapping. Geoderma 409, 115638 (2022).
Article ADS Google Scholar
Galzki, J. C., Birr, A. S. & Mulla, D. J. Identifying critical agricultural areas with three-meter LiDAR elevation data for precision conservation. J. Soil Water Conserv. 66, 423 LP–430 (2011).
Seelan, S. K., Laguette, S., Casady, G. M. & Seielstad, G. A. Remote sensing applications for precision agriculture: A learning community approach. Remote Sens. Environ. 88, 157–169 (2003).
Article ADS Google Scholar
Vepraskas, M. J. & Caldwell, P. V. Interpreting morphological features in wetland soils with a hydrologic model. Catena (Amst) 73, 153–165 (2008).
Article Google Scholar
Vepraskas, M. J., Heitman, J. L. & Austin, R. E. Future directions for hydropedology: Quantifying impacts of global change on land use. Hydrol. Earth Syst. Sci. 13, 1427–1438 (2009).
Article ADS Google Scholar
Thompson, J. A., Bell, J. C. & Butler, C. A. Quantitative soil-landscape modeling for estimating the areal extent of hydromorphic soils. Soil Sci. Soc. Am. J. 61, 971–980 (1997).
Article ADS CAS Google Scholar
Anaya-Acevedo, J. A. et al. Identification of wetland areas in the context of agricultural development using Remote Sensing and GIS. Dyna (Medellin) 84, 186–194 (2017).
Article Google Scholar
Salinas, J. B. G. et al. Wetland mapping with multitemporal sentinel radar remote sensing in the southeast region of Brazil. in 2020 IEEE Latin American GRSS & ISPRS Remote Sensing Conference (LAGIRS) 669–674. https://doi.org/10.1109/LAGIRS48042.2020.9165593 (2020).
Valencia, E. et al. Wetland monitoring technification for the Ecuadorian Andean region based on a multi-agent framework. Heliyon 8, e09054 (2022).
Article PubMed PubMed Central Google Scholar
Furlan, L. M., Rosolen, V., Moreira, C. A., Bueno, G. T. & Ferreira, M. E. The interactive pedological-hydrological processes and environmental sensitivity of a tropical isolated wetland in the Brazilian Cerrado. SN Appl. Sci. 3, 144 (2021).
Article CAS Google Scholar
Chaplot, V., Walter, C. & Curmi, P. Improving soil hydromorphy prediction according to DEM resolution and available pedological data. Geoderma 97, 405–422 (2000).
Article ADS Google Scholar
Chaplot, V., Walter, C. & Curmi, P. Testing quantitative soil-landscape models for predicting the soil hydromorphic index at a regional scale. Soil Sci. 168 (2003).
Lobell, D. B. & Asner, G. P. Moisture effects on soil reflectance. Soil Sci. Soc. Am. J. 66, 722–727 (2002).
Article ADS CAS Google Scholar
Demattê, J. A. M., Fongaro, C. T., Rizzo, R. & Safanelli, J. L. Geospatial Soil Sensing System (GEOS3): A powerful data mining procedure to retrieve soil spectral reflectance from satellite images. Remote Sens. Environ. 212, 161–175 (2018).
Article ADS Google Scholar
Demattê, J. A. M. et al. The Brazilian Soil Spectral Library (BSSL): A general view, application and challenges. Geoderma 354, 113793 (2019).
Article ADS Google Scholar
Sherman, D. M. & Waite, T. D. Electronic spectra of Fe3+ oxides and oxide hydroxides in the near IR to near UV. Am. Miner. 70, 1262–1269 (1985).
ADS CAS Google Scholar
Scheinost, A. C., Chavernas, A., Barrón, V. & Torrent, J. Use and limitations of second-derivative diffuse reflectance spectroscopy in the visible to near-infrared range to identify and quantify Fe oxide minerals in soils. Clays Clay Miner. 46, 528–536 (1998).
Article ADS CAS Google Scholar
Duchaufour, P. Hydromorphic soils BT—Pedology: Pedogenesis and classification. in (ed. Duchaufour, P.) 335–372 (Springer Netherlands, 1982). https://doi.org/10.1007/978-94-011-6003-2_12.
Gallo, B. C. et al. Multi-temporal satellite images on topsoil attribute quantification and the relationship with soil classes and geology. Remote Sens (Basel) 10, (2018).
van Breemen, N. & Buurman, P. Soil Formation—Second Edition. (2002).
Haubrock, S.-N., Chabrillat, S., Lemmnitz, C. & Kaufmann, H. Surface soil moisture quantification models from reflectance data under field conditions. Int. J. Rem. Sens. 29, 3–29 (2008).
Article Google Scholar
Safanelli, J. L. et al. Terrain analysis in google earth engine: A method adapted for high-performance global-scale analysis. ISPRS Int. J. Geo-Inf. 9. https://doi.org/10.3390/ijgi9060400 (2020).
Naimi, S. et al. Spatial prediction of soil surface properties in an arid region using synthetic soil image and machine learning. Geocarto Int. 1–22. https://doi.org/10.1080/10106049.2021.1996639 (2021).
Schaetzl, R. & Anderson, S. Soils. Genesis and Geomorphology. (Cambridge University Press, 2005).
Shaheen, S. M. et al. Stepwise redox changes alter the speciation and mobilization of phosphorus in hydromorphic soils. Chemosphere 288, 132652 (2022).
Article ADS CAS PubMed Google Scholar
Giannetta, B., Oliveira de Souza, D., Aquilanti, G., Celi, L. & Said-Pullicino, D. Redox-driven changes in organic C stabilization and Fe mineral transformations in temperate hydromorphic soils. Geoderma 406, 115532 (2022).
Amundson, R. Soil biogeochemistry and the global agricultural footprint. Soil Secur. 6, 100022 (2022).
Article Google Scholar
Kaiser, I. M., Bezerra, B. S. & Castro, L. I. S. Is the environmental policies procedures a barrier to development of inland navigation and port management? A case of study in Brazil. Transp. Res. Part A Policy Pract. 47, 78–86 (2013).
Article Google Scholar
Setzer, J. Environmental paradiplomacy: The engagement of the Brazilian state of São Paulo in international environmental relations. Preprint at (2013).
Tornquist, C. G. & Broetto, T. Protection of the Soil Resource in the Brazilian Environmental Legislation BT – Global Soil Security. in (eds. Field, D. J., Morgan, C. L. S. & McBratney, A. B.) 397–401 (Springer International Publishing, 2017). https://doi.org/10.1007/978-3-319-43394-3_36.
Bordonal, R. de O. et al. Sustainability of sugarcane production in Brazil. A review. Agron Sustain. Dev. 38, 13 (2018).
Bonato, E. R. & Dall’Agnol, A. Soybean in Brazil-production and research. in World Soybean Research Conference III: Proceedings 1248–1256 (CRC Press, 2022).
Kögel-Knabner, I. et al. Biogeochemistry of paddy soils. Geoderma 157, 1–14 (2010).
Article ADS Google Scholar
Pott, A., Pott, V. J., Catian, G. & Scremin-Dias, E. Floristic elements as basis for conservation of wetlands and public policies in Brazil: The case of veredas of the Prata River. Oecol. Austral. 23, (2019).
Valera, C. A. et al. The buffer capacity of riparian vegetation to control water quality in anthropogenic catchments from a legally protected area: A critical view over the Brazilian new forest code. Water 11. https://doi.org/10.3390/w11030549 (2019).
Alvares, C. A., Stape, J. L., Sentelhas, P. C., De Moraes Gonçalves, J. L. & Sparovek, G. Köppen’s climate classification map for Brazil. Meteorologische Zeitschrift 22, 711–728 (2013).
Santos, H. G. et al. Procedimentos normativos de levantamentos pedológicos. (EMBRAPA, 1995).
Demattê, J. A. M. et al. Genesis and properties of wetland soils by VIS-NIR-SWIR as a technique for environmental monitoring. J. Environ. Manage 197, 50–62 (2017).
Article PubMed Google Scholar
Marques, K. P. et al. How qualitative spectral information can improve soil profile classification?. J. Near. Infrared. Spectrosc. 27, 156–174 (2019).
Article ADS CAS Google Scholar
Demattê, J. A. M., Bellinaso, H., Romero, D. J. & Fongaro, C. T. Morphological Interpretation of Reflectance Spectrum (MIRS) using libraries looking towards soil classification. Sci. Agric. 71, 509–520 (2014).
Article Google Scholar
Jenny, H. Factors of soil formation: a system of quantitative pedology. (Courier Corporation, 1994).
Gorelick, N. et al. Google earth engine: Planetary-scale geospatial analysis for everyone. Remote Sens Environ 202, 18–27 (2017).
Article ADS Google Scholar
Teng, H., Viscarra Rossel, R. A., Shi, Z. & Behrens, T. Updating a national soil classification with spectroscopic predictions and digital soil mapping. Catena (Amst) 164, 125–134 (2018).
Zeraatpisheh, M., Ayoubi, S., Jafari, A., Tajik, S. & Finke, P. Digital mapping of soil properties using multiple machine learning in a semi-arid region, central Iran. Geoderma 338, 445–452 (2019).
Article ADS CAS Google Scholar
Breiman, L. Random forests. Mach. Learn. 45, 5–32 (2001).
Article MATH Google Scholar
Mosier, C. I. I. Problems and designs of cross-validation 1. Educ. Psychol. Meas. 11, 5–11 (1951).
Article Google Scholar
Browne, M. W. Cross-validation methods. J. Math. Psychol. 44, 108–132 (2000).
Article MathSciNet CAS PubMed MATH Google Scholar
Efron, B. Bootstrap methods: Another look at the jackknife. in Breakthroughs in statistics 569–593 (Springer, 1992).
Bischl, B., Mersmann, O., Trautmann, H. & Weihs, C. Resampling methods for meta-model validation with recommendations for evolutionary computation. Evol. Comput. 20, 249–275 (2012).
Article CAS PubMed Google Scholar
Kuhn, M. Building predictive models in R using the caret package. J. Stat. Softw. 1(5), 1. https://doi.org/10.18637/jss.v028.i05 (2008).
Article Google Scholar
R Core Team. R: A language and environment for statistical computing (2013).
Congalton, R. G. & Green, K. Assessing the accuracy of remotely sensed data principles and practices (Taylor & Francis Group, 2019). https://doi.org/10.1017/CBO9781107415324.004.
Download references
This research was funded by Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES)- Programa de Excelência Acadêmica (PROEX)—Brazil. This research was also funded by São Paulo Research Foundation, grants number 2014/2262-2. The authors are grateful to the Soil Science Department in Luiz de Queiroz College of Agriculture, University of São Paulo for funding and providing the facilities for the development of this research. The authors are also grateful to the Geotechnologies in Soil Science group for providing the soil database. The third author thanks the Coordination of Integrate Technical Assistance (CATI/SAA).
Department of Soil Science, Luiz de Queiroz College of Agriculture, University of São Paulo, Pádua Dias Av., 11, Postal Box 09, Piracicaba, São Paulo, 13416-900, Brazil
Fellipe A. O. Mello, José A. M. Demattê, Henrique Bellinaso, Raul R. Poppiel, Rodnei Rizzo, Nícolas Augusto Rosin, Jorge T. F. Rosas, Nélida E. Q. Silvero & Heidy S. Rodríguez-Albarracín
Coordination of Integrate Technical Assistance of Secretariat of Agriculture and Supply-CATI/SAA, Piracicaba Regional, Campos Salles Street, 507, Piracicaba, São Paulo State, 13400-200, Brazil
Henrique Bellinaso
Department of Soil Science, Federal University of Viçosa, Peter Henry Rolfs Av.-University Campus, Viçosa, Minas Gerais, 36570-900, Brazil
Danilo C. de Mello
You can also search for this author in PubMed Google Scholar
You can also search for this author in PubMed Google Scholar
You can also search for this author in PubMed Google Scholar
You can also search for this author in PubMed Google Scholar
You can also search for this author in PubMed Google Scholar
You can also search for this author in PubMed Google Scholar
You can also search for this author in PubMed Google Scholar
You can also search for this author in PubMed Google Scholar
You can also search for this author in PubMed Google Scholar
You can also search for this author in PubMed Google Scholar
Conceptualization: F.A.O.M. and J.A.M.D. Project development and writing: F.A.O.M., J.A.M.D. and H.B. Methodology: R.R.P., R.R., D.C.M., N.A.R., J.T.F.R. and H.S.R.A. data curation: R.R.P., J.T.F.R., R.R. and H.S.R.A. Supervision: J.A.M.D. Review and discussion: F.A.O.M., J.A.M.D. Editing: F.A.O.M., R.R.P., R.R. Funding acquisition: J.A.M.D. All authors have read and agreed to the published version of the manuscript.
Correspondence to José A. M. Demattê.
The authors declare no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.
Reprints and Permissions
Mello, F.A.O., Demattê, J.A.M., Bellinaso, H. et al. Remote sensing imagery detects hydromorphic soils hidden under agriculture system. Sci Rep 13, 10897 (2023). https://doi.org/10.1038/s41598-023-36219-9
Download citation
Received:
Accepted:
Published:
DOI: https://doi.org/10.1038/s41598-023-36219-9
Anyone you share the following link with will be able to read this content:
Sorry, a shareable link is not currently available for this article.
Provided by the Springer Nature SharedIt content-sharing initiative
By submitting a comment you agree to abide by our Terms and Community Guidelines. If you find something abusive or that does not comply with our terms or guidelines please flag it as inappropriate.
Advertisement
© 2023 Springer Nature Limited
Sign up for the Nature Briefing newsletter — what matters in science, free to your inbox daily.