SciELO - Scientific Electronic Library Online

 
vol.88 issue218Alkaline delignification of lignocellulosic biomass for the production of fermentable sugar syrupsDesign and validation of a prototype to calibrate the pressure variable in tonometry author indexsubject indexarticles search
Home Pagealphabetic serial listing  

Services on Demand

Journal

Article

Indicators

Related links

  • On index processCited by Google
  • Have no similar articlesSimilars in SciELO
  • On index processSimilars in Google

Share


DYNA

Print version ISSN 0012-7353On-line version ISSN 2346-2183

Dyna rev.fac.nac.minas vol.88 no.218 Medellín July/Sept. 2021  Epub Mar 28, 2022

https://doi.org/10.15446/dyna.v88n218.91077 

Articles

Characterization of alluvial soil hydrodynamics in the upper Ipanema river basin using the Beerkan method

Caracterización de la hidrodinámica del suelo aluvial en la cuenca superior del río Ipanema mediante el método Beerkan

Antonio Gebson Pinheiro a  
http://orcid.org/0000-0002-5912-6906

Alexandre Maniçoba da Rosa Ferraz Jardim a  
http://orcid.org/0000-0001-7094-3635

Abelardo Antônio de Assunção Montenegro a  
http://orcid.org/0000-0002-5746-8574

Thieres George Freire da Silva a   b  
http://orcid.org/0000-0002-8355-4935

José Raliuson Inácio Silva a  
http://orcid.org/0000-0003-04883-0514

a Department of Agricultural Engineering, Federal Rural University of Pernambuco, Recife, Brazil. gebson10@hotmail.com, alexandremrfj@gmail.com, abelardo.montenegro@ufrpe.br, raliuson.agro@gmail.com

b Department of Agronomy, Federal Rural University of Pernambuco, Serra Talhada, Brazil. thigeoprofissional@hotmail.com


Abstract

The objective was to model the spatial distribution of the physical properties of the soil through geostatistics in the Brazilian semi-arid region. The study was carried out in the Experimental Basin of Ipanema River, Pernambuco, in an alluvial area with regular grid and samples of 40 points. Infiltration tests were carried out to assess hydraulic conductivity and sampling was performed to determine soil texture and soil organic carbon (SOC). The variables were subjected to descriptive and correlation statistical analyses, in addition to the geostatistical report, which indicated high variability for the hydraulic conductivity of the soil. The best fits were obtained with the Spherical model for the hydraulic conductivity, sand, silt, and clay variables, and the Gaussian model for SOC. Thus, with the geostatistics technique, it was possible to map and verify the main physical characteristics of the soil in the area, providing greater optimization in agricultural planning.

Keywords: geostatistics; variogram; spatial dependence; hydraulic conductivity

Resumen

El objetivo fue modelar la distribución espacial de las propiedades físicas del suelo a través de geoestadísticas en la región semiárida brasileña. El estudio se realizó en la Cuenca Experimental del Río Ipanema, Pernambuco, en una zona aluvial con malla regular y muestras de 40 puntos. Se realizaron pruebas de infiltración para evaluar la conductividad hidráulica y muestreos para determinar la textura del suelo y el carbono orgánico del suelo (COS). Las variables fueron sometidas a análisis estadístico descriptivo y de correlación, además del informe geoestadístico, que indicó alta variabilidad para la conductividad hidráulica del suelo. Los mejores ajustes se realizaron con el modelo esférico para las variables de conductividad hidráulica, arena, limo y arcilla, y el modelo gaussiano para el COS. Así, con la técnica de geoestadística, fue posible mapear y verificar las principales características físicas del suelo en la zona, proporcionando una mayor optimización en la planificación agrícola.

Palabras clave: geoestadística; variograma; dependencia espacial; conductividad hidráulica

1. Introduction

Modern cultivation techniques with high land use rates and intense mechanization cause changes in soil attributes, influencing production, natural resource balance, and soil water dynamics [8,14,32]. Inadequate cultivation causes soil degradation due to compaction, degradation of organic matter, and modification of the structure, negatively altering the physical, hydraulic, chemical, and biological properties of the soil, causing damage to crop yield [2,27].

The evaluation of hydrodynamic properties in the soils of an area assists in the study of soil behavior in areas under agricultural management, degraded and recovering [19], highlighting the texture, pore size distribution, compaction, surface sealing, runoff, soil infiltration and soil hydraulic conductivity [24]. Soil physical quality can be designated by the interaction of factors that directly affect crop development through soil water infiltration [16]. Thus, knowledge of soil hydrodynamic characteristics, such as saturated hydraulic conductivity, is essential for the management of watersheds [2,26].

The characterization of the hydrodynamic properties of the soil by the Beerkan methodology has been widely used because it has field application, low cost, and simple execution, just the soil water infiltration curve and soil texture [3,27,31,32]. Currently, in the literature, there are studies with the Beerkan methodology to evaluate the spatial variability of soil hydraulic conductivity [2,10,23,25,29,32]. However, few studies are carried out in alluvial soils, which are very fertile; due to this characteristic, intensive and highly heterogeneous agriculture is practiced.

The spatial distribution of physical, chemical and biological soil properties is directly related to intrinsic variations in soil characteristics, vegetation and climate. But they can also be affected by exogenous anthropogenic factors associated with agricultural production, such as soil preparation and fertilization [5,10,17,29].

Traditional statistical methods are often used to describe the variability of soil properties in the field and find that the variation occurs at random in the samples. The application of geostatistics, based on the theory of regionalized variables, is more efficient in characterizing the spatial variability of the variables, allowing the identification of spatial dependence and the creation of maps, estimating values in unsampled locations, based on semivariogram models [4], being particularly relevant for highly heterogeneous fields.

The spatial distribution of soil attributes can be evaluated by the use of geostatistics, which is applied mainly in geological and hydrological fields and may represent regions of different permeability in soils [11,21]. Thus, the objective was to evaluate the spatial dependence of soil particle size (sand, silt, and clay) and saturated hydraulic conductivity in compacted alluvial soil, applying the Beerkan methodology.

2. Material and methods

2.1 Study area characterization and Sampling

The study was carried out in an alluvial area of Vale do Mimoso, located in the Upper Ipanema River basin, where small-scale irrigated agriculture is conducted, located in the municipality of Pesqueira, State of Pernambuco, Brazil (Fig. 1). The location is found at the following geographic coordinates: longitude −36º87', latitude -8º39' and altitude of 645 m.

Source: The Authors.

Figure 1 Geographic location of the study area in the municipality of Pesqueira-PE, Brazil, and sampling points in the experimental area. 

Upper Ipanema River Basin, next to the Hydrology Network of the Semi-Arid Region (REHISA), is a sub-basin of the Ipanema River system in its westernmost portion. This basin covers part of the municipalities of Arcoverde and Pesqueira, in the state of Pernambuco, with an area of 194.82 km2. The basin has an average yearly rainfall of 672 mm, an average annual air temperature of 27.4 ºC, and average potential evapotranspiration of 1600 mm year−1 [20]. According to the Brazilian Soil Classification System (SBCS), the basin soils are distinguished in Yellow Argisol, Red-Yellow Argisol, Fluvic Neosol, Regolithic Neosol, and Litholic Neosol. However, in the present study, the soil type of the area is Fluvic Neosol. The climate of the region is BSsh type (i.e., extremely hot, semi-arid), according to Köppen’s classification.

The experiment was conducted in an area with a regular grid spaced by 8.0 x 8.0 m, forming an 80 x 50 grid, totaling 40 sampling points.

2.2 Physical and hydraulic properties

Infiltration tests and soil sampling were performed in the sampled grid to obtain physical and hydraulic properties using the Beerkan method [18]. The infiltration tests were performed with a single-ring infiltrometer, providing three-dimensional infiltration. The tests were performed according to the following operating procedures:

  • a) A PVC (Polyvinyl chloride) cylinder with a diameter of 200 mm was positioned on the ground surface and inserted to a depth of approximately 1.0 cm. The method takes into account only the hydrodynamic characteristics of the topsoil, and the deepening of the cylinder is only to prevent lateral water loss during the infiltration process.

  • b) Fifteen volumes with 200 mL of water were prepared and poured into the cylinder. The first volume is poured into the cylinder, and immediately after the cylinder is emptied, then another volume of water is poured into it, always taking note of the time each volume takes to infiltrate the soil. This procedure is repeated successively for the 15 volumes until the last volume of water is poured into the cylinder, and the test is terminated if the infiltration rate is constant.

For soil texture (sand, silt, and clay) and organic carbon analyses, soil collections were performed at each point studied in the 0 to 0.20 m layer. The samples were placed in plastic bags and taken to the Water and Soil Laboratory of the Federal Rural University of Pernambuco.

The soil organic carbon (SOC) of the soil was determined by the oxidation of organic matter, using 0.0067 M potassium dichromate, in sulfuric medium, and titration by 0.1 mol L−1 ammoniacal ferrous sulfate [34].

The texture analyses were performed from samples that were air-dried for 72 hours, pounded to break up clods and sieved through a 2.0-mm-mesh sieve. The sand, silt, and clay granulometric fractions were determined by the Bouyoucos hydrometer method using a mechanical stirrer as a physical dispersant for 15 minutes for clayey soils and 5.0 minutes for sandy soils. 25 mL of sodium hexametaphosphate (Calgon) was used as a chemical dispersant [34].

From the infiltration data and physical properties of the soil, the Beerkan method was applied to determine the hydrodynamic parameters at the local scale in the field and, therefore, estimate the hydraulic conductivity of the soil (Ks) [10].

This method was devised by [15] and it has an advantage over other experimental methods due to its application, which is simpler, faster and requires a smaller number of people and equipment, making it more economical.

The method is used to mathematically represent the infiltration phenomenon by fitting the infiltration curves to determine the parameters of soil saturated hydraulic conductivity (Kfs) according to Eq. 1 [3].

Where, b1 is a parameter of the linearized equation of the line; r is the radius of the perforated disk (m). Regarding the values of α*, four are suggested* (0.036; 0.012; 0.004 and 0.001 mm) for the practical use of application of the Guelph permeameters and disc infiltrometers in soils ranging from coarse sands to compacted clays. An α* value of 0.012 mm was considered to be the first approximation for most field soils [22].

2.3 Data analysis

Descriptive statistics was used to determine the mean, median, standard deviation, maximum and minimum values, skewness and kurtosis coefficients and the coefficient of variation. For information on dispersion and distribution of variables, data normality was tested by the Kolmogorov-Smirnov test at the 5% probability level. Additionally, linear correlations were performed between the physical attributes of the soil by Pearson's correlation test.

The variability of physical properties was classified based on the coefficient of variation (CV) values, according to [36], as low (CV < 12%), medium (12 < CV < 62%) and high (CV > 62%).

To evaluate the presence of trend in the data, two trend surfaces were fitted, which had the following determination coefficients (R2): first-degree trend surface - 18.70%, second-degree trend surface - 25.45% [35].

2.4 Geostatistical analysis

Spatial dependence analysis was performed using geostatistics and semivariogram fitting. Based on the stationary assumption of the intrinsic hypothesis, the spatial autocorrelation between neighboring sites is calculated through the semivariance, which is estimated by Eq. 2a and 2b.

where, γ(h) is the semivariance; γZW(h) is the cross semivariance, h is the separation distance (m); N(h) is the number of data pairs separated by distance h; Z(xi), W(xi) are measured values of the variables at location xi, and Z(xi + h), W(xi + h) are measured values of the variables at a distance h from the location xi.

For construction and fitting of semivariograms, the GS+ program version 10.0 was used (Gamma Design Software, LLC, Plainwell, Michigan, USA) [27]. The criteria used and procedures performed for fitting the semivariogram model were those described in [36].

After obtaining the experimental semivariograms using the GEO-EAS software (Geostatistical Environmental Assessment Software), the Gaussian, spherical, and exponential models were tested. By fitting a mathematical model to the calculated values, the coefficients of the theoretical model for the semivariogram are estimated (i.e., the nugget effect, C0, sill, C0 + C1, and range, A).

The model that showed an adequate fit to experimental values and standardized errors with a mean close to zero and standard deviation close to one was chosen, according to the Jack-Knifing cross-validation technique [35].

To determine the degree of spatial dependence (DSD) of the variables under study (Eq. 3), the classification by [6] was performed, which is based on the analysis of the relationship between the nugget effect and the sill of its fitted semivariogram. Values below 25% characterize strong spatial dependence, between 25% and 75% moderate dependence, and above 75% weak dependence.

Where, C0 is the nugget variance and C1 is the structural variance. Effective range is the distance at which the sill (C0 + C1) is within 5% of the asymptote.

The isoline maps were constructed using the kriging algorithm, in the computer program Surfer 9.0 [12].

3. Results and discussion

The results of the descriptive statistical analysis for saturated hydraulic conductivity and texture are shown in Table 1. The Kolmogorov-Smirnov test indicated that the hydraulic conductivity and clay data did not show normality at 5% probability level.

Table 1 Descriptive statistics for data of log-normal saturated hydraulic conductivity (Kfs) and sand, silt and log-normal clay content. 

CV: coefficient of variation

Source: The Authors.

Through analysis of the coefficient of variation (Table 1), it is observed that the log of hydraulic conductivity, the log of clay and silt had coefficients of variation (CV) classified as medium (12 < and < 62%), while the sand obtained low CV (< 12%), according to the criteria of [37]. [25], analyzing the spatial variability of hydrodynamic properties of a Regolithic Neosol under pasture and Caatinga, observed the same behavior. The coefficient of variation can be an essential tool for assessing and comparing variables, as low CV variables indicate a small variation of data from the mean, showing a homogeneous distribution [19].

From the Pearson correlation matrix, it is possible to observe that there was no significant correlation for most of the physical indicators of the soil in the present study, highlighting the low correlation of clay with saturated hydraulic conductivity and organic carbon in the soil (Fig. 2). Behavior similar to that was found by [2], in a small basin in Sicily, where the saturated hydraulic conductivity data showed no significant relationship with the textural indicators of the soil. However, controversial results were reported by [6], who found significant correlations for most of the physical attributes of the soil in the Atlantic Forest biome, northeastern Brazil.

Source: The Authors.

Figure 2 Pearson's general correlation matrix for soil physical attributes. logKfs - log-normal saturated hydraulic conductivity; Sand; Silt, log-normal Clay; SOC - soil organic carbon. The largest values represent an increase in the significance of Pearson's correlation coefficient (p-value < 0.05). 

All variables evaluated showed a high degree of spatial dependence (Table 2), directly influenced by intrinsic properties of the soil such as density, porosity, soil moisture and soil texture. Attributes with low spatial dependence are more influenced by extrinsic factors, such as relief of the site, soil management, trampling of animals and compaction due to agricultural machinery traffic [7]. Thus, a moderate degree of spatial dependence (as found in the present study) shows that the soil suffers much more from the influences caused by the lack of management and traffic of agricultural machines than from those caused by the other soil characteristics (e.g. density, porosity, and texture of soil).

The mathematical model that best fitted the distribution for the organic carbon results of the soil was the Gaussian one, whereas the other physical variables of the soil were better described by the spherical model (Table 2). The spherical model has already been found in other works as expressing the best fit for some soil attributes [13,16].

Table 2 Parameters of semivariogram models and degree of spatial dependence (DSD) for data of log-normal saturated hydraulic conductivity (Kfs) and sand, silt and log-normal clay content. 

C0: Nugget Effect; C0+C1: Sill; A: Range; DSD: Degree of spatial dependence; R²: Coefficient of determination; SOC: Soil organic carbon; Sph: Spherical model; and Gau: Gaussian model.

Source: The Authors.

The experimental semivariance of the logarithms of saturated hydraulic conductivity was estimated (Table 2). The semivariogram model that best fitted the hydraulic conductivity semivariances was the Spherical, with a nugget effect of 0.0028 and a range of 12.5 m. For the particle-size fractions, the best fitted mathematical model was also the Spherical, with ranges of 26.06 m, 21.45 m and 13.00 m for sand, silt and clay, respectively (Table 2). For soil organic carbon, the mathematical model that best fitted was the Gaussian, with a nugget effect equivalent to 12.50 and a range of 22.39 m.

Although classical correlation was not observed, spatial cross-correlation occurs between log (Kfs) and texture, as verified based on cross-semivariograms. Fig. 3 shows cross-semivariograms between log Kfs and sand content and between log Kfs and silt.

Source: The Authors.

Figure 3 Cross-semivariograms between log Kfs and sand content and between log Kfs and silt. 

[25], working with data using the log of saturated hydraulic conductivity, obtained the best results of the semivariance fits with the Spherical mathematical model, whose nugget effect was 0.004 and range was 35.38 m, in a 20 x 20 m grid. Higher range values were found by [21], evaluating the spatial variability of hydraulic conductivity of plain alluvial soils, with a range of 320 m with an average spacing of 30 m between collection points.

From the soil texture and the analysis of saturated hydraulic conductivity (Ks), by the Beerkan method, [1,29] also observed the main correlation between Ks and the sand content, even with low results. However, these studies did not assess the correlation between these variables through cross-validation, so only a one-time analysis between the variables was possible and not an analysis of spatial correlation. This fact highlights the importance of the cross-validation analysis of the present study.

According to the method of [38], there was a trend in the saturated hydraulic conductivity data, requiring their fit. Fig. 4 shows the isoline maps for the saturated hydraulic conductivity of the soil. When analyzing the spatial distribution of the area, it is noted that the lowest values of soil Kfs are located near the left edge (Fig. 4). Such distribution can be explained by the high intensity of circulation of vehicles and agricultural machinery at this edge. It can be seen that the highest hydraulic conductivity values (Fig. 4A) were located in places with more significant sand fractions (see Fig. 4B). There is a positive association between sand content and soil macroporosity [32], which in turn can influence the soil to conduct more water.

Source: The Authors.

Figure 4 Maps of the spatial distribution of properties of the soil: (A) log-normal saturated hydraulic conductivity (mm s−1); (B) sand (g kg−1); (C) silt (g kg−1); (D) log-normal clay (g kg−1); and (E) soil organic carbon (g kg−1). 

High variability in hydraulic conductivity is observed, which can be explained by the influence of local soil heterogeneity. Similar responses were found by [10], who associate high variability more with soil structural factors than with textural factors.

4. Conclusions

The best fitted mathematical model for soil hydraulic conductivity was spherical, with a 12.50 meter range.

Soil physical indicators showed high spatial dependence.

Although there was no point to point relationship between sand and silt contents and hydraulic conductivity, crossed spatial variation occurs for those variables.

The highest values of hydraulic conductivity are found in regions with higher sand content in the soil.

Geostatistics made it possible to identify spatial cross-correlation between hydraulic conductivity and texture and to map regions of lower hydraulic conductivity, where specific management techniques should be adopted so that the conductivity can be increased.

Acknowledgments

The authors acknowledge the support from the Research Support Foundation of the Pernambuco State (FACEPE- APQ 0300-5.03/17), the National Council for Scientific and Technological Development (CNPq), which also funds T.G.F. Silva and A.A.A. Montenegro through the fellowship of the research productivity, and the Coordination for the Improvement of Higher Education Personnel (CAPES - Finance Code 001) for the research and study grants.

References

[1] Araújo, D.C.S., Montenegro, S.M.G.L., Montenegro, A.A.A., Santos, D.P. and Rodrigues, R.A.S., Temporal stability of soil moisture in banana cropping area in the Brazilian semiarid región, Revista Brasileira de Engenharia Agrícola e Ambiental, 23(11), pp. 852-859, 2019. DOI: 10.1590/1807-1929/agriambi.v23n11p852-859 [ Links ]

[2] Bagarello, V., Baiamonte, G., Caia, C. Variability of near-surface saturated hydraulic conductivity for the clay soils of a small Sicilian basin. Geoderma, 340, pp. 133-145, 2019. https://doi.org/10.1016/J.GEODERMA.2019.01.008Links ]

[3] Bagarello, V., Di Prima, S., Iovino, M., Provenzano, G. Estimating field-saturated soil hydraulic conductivity by a simplified Beerkan infiltration experiment. Hydrological Processes, 28 (3), pp. 1095-1103, 2014. https://doi.org/10.1002/hyp.9649Links ]

[4] Bhunia, G. S., Shit, P. K., Chattopadhyay, R. Assessment of spatial variability of soil properties using geostatistical approach of lateritic soil (West Bengal, India). Annals of Agrarian Science, 16 (4), pp. 436-443, 2018. https://doi.org/10.1016/j.aasci.2018.06.003Links ]

[5] Castellini, M., Stellaci, A. M., Sisto, D., Iovino, M. The Mechanical Impact of Water Affected the Soil Physical Quality of a Loam Soil under Minimum Tillage and No-Tillage: An Assessment Using Beerkan Multi-Height Runs and BEST-Procedure. Land, 10 (195), pp. 1-16, 2021. https://doi.org/10.3390/land10020195Links ]

[6] Cavalcanti, R. Q., Rolim, M. M., Lima, R. P., Tavares, U. E., Pedrosa, E. M. R., Cherubin, M. R. Soil physical changes induced by sugarcane cultivation in the Atlantic Forest biome, northeastern Brazil. Geoderma. 370, pp. 1-11, 2020. https://doi.org/10.1016/j.geoderma.2020.114353Links ]

[7] Cambardella, C.A., Moorman, T. B., Parkin, T. B., Karlen, D. L., Novak, J. M., Turco, R. F., Konopka, A. E. Field-Scale Variability of Soil Properties in Central Iowa Soils. Soil Science Society of America Journal, 58, pp. 1501-1511, 1994. https://doi.org/10.2136/sssaj1994.03615995005800050033xLinks ]

[8] Ferraz Jardim, A. M. R., Silva, J. R. I., Vieira Leite, M. L. M., Teixeira, V. I., Morato, M. P., Araújo Júnior, G. N., Silva, T. G. F. Symbiotic interaction in forage crop cultivations: A review. Amazonian Journal of Plant Research, 2 (1), pp. 149-160, 2018. 10.26545/ajpr.2018.b00019x [ Links ]

[9] Ferraz Jardim, A. M. R., Silva, T. G. F., Souza, L. S. B., Araújo Júnior, G. N., Alves, H. K. M. N., Souza, M. S., Araújo, G. G. L., Moura, M. S. B. Intercropping forage cactus and sorghum in a semi-arid environment improves biological efficiency and competitive ability through interspecific complementarity. Journal of Arid Environments, 188 (104464), pp 1-11, 2021. https://doi.org/10.1016/j.jaridenv.2021.104464Links ]

[10] Furtunato, O. M., Montenegro, S. M. G. L., Antonino, A. C. D., Oliveira, L. M. M., Souza, E. S., Moura, A. E. S. S. Variabilidade Espacial de Atributos Físico-hídricos de Solos em uma Bacia Experimental no Estado de Pernambuco. Revista Brasileira de Recursos Hídricos, 18 (2), pp. 135-147, 2013. https://doi.org/10.21168/rbrh.v18n2.p135-147Links ]

[11] Godoy, V. A., Zuquette, L. V., Gómez-Hernández, J. J. Spatial variability of hydraulic conductivity and solute transport parameters and their spatial correlations to soil properties. Geoderma 339, pp. 59-69, 2019. https://doi.org/10.1016/J.GEODERMA.2018.12.015Links ]

[12] Golden Software, 2010. Surfer for windows version 9.0. [ Links ]

[13] Guedes Filho, O., Vieira, S. R., Chiba, M. K., Nagumo, C. H., Dechen, S. C. F. Spatial and temporal variability of crop yield and some Rhodic Hapludox properties under no-tillage. Revista Brasileira Ciência do Solo, 34 (1), pp. 1-14, 2010. https://doi.org/10.1590/S0100-06832010000100001Links ]

[14] Haas, J., Fenner, P. T., Schack-Kirchner, H., Lang, F. Quantifying soil movement by forest vehicles with corpuscular metal tracers. Soil Tillage Research. 181, pp. 19-28, 2018. https://doi.org/10.1016/J.STILL.2018.03.012Links ]

[15] Haverkamp, R., Bouraoui, F., Angulo-Jaramillo, R., Zammit, C., Delleur, J. W. Soil properties and moisture movement in the unsaturated zone, in: The Handbook of Groundwater Engineering. CRC Press, Taylor & Francis Group, 6000 Broken Sound Parkway NW, Suite 300, Boca Raton, FL 33487-2742, pp. 413-440, 1998. https://doi.org/10.1201/9781315371801-15Links ]

[16] José, J. V., Rezende, R., Marques, P. A. A., Gonçalves, A. C. A., Souza, R. S. Variabilidade espacial de variáveis físico-hídricas de dois Latossolos da região noroeste do estado do Paraná. Irriga, 17 (2), pp. 208-219, 2012. https://doi.org/10.15809/irriga.2012v17n2p208Links ]

[17] Lan, G., Li, Y., Lesueurd, D., Wu, Z., Xie, G. Seasonal changes impact soil bacterial communities in a rubber plantation on Hainan Island, China. Science of the Total Environment. 626, pp. 826-834, 2018. https://doi.org/10.1007/s00374-010-0478-4Links ]

[18] Lassabatere, L., Angulo-Jaramillo, R., Soria Ugalde, J. M., Cuenca, R., Braud, I., Haverkamp, R. Beerkan estimation of soil transfer parameters through infiltration experiments. Soil Science Society of America Journal, 70 (2), pp. 521-532, 2006. https://doi.org/10.2136/sssaj2005.0026Links ]

[19] Lima, F. V., Silvino, G. S., Melo, R. S. S. M., Lira, E. C., Ribeiro, T. S. Variabilidade espacial de atributos físicos do solo em área de encosta sob processo de degradação. Revista Caatinga, 28 (4), pp. 53-63, 2015. https://doi.org/10.1590/1983-21252015v28n406rcLinks ]

[20] Lopes, I., Montenegro, A. A. A., Lima, J. L. M. P. Performance of Conservation Techniques for Semiarid Environments: Field Observations with Caatinga, Mulch, and Cactus Forage Palma. Water, 11 (4), pp. 792, 2019. https://doi.org/10.3390/w11040792Links ]

[21] Montenegro, A. A. A., Montenegro, S. M. G. L. Variabilidade espacial de classes de textura, salinidade e condutividade hidráulica de solos em planície aluvial. Revista Brasileira de Engenharia Agrícola e Ambiental. 10 (1), pp. 30-37, 2006. https://doi.org/10.1590/S1415-43662006000100005Links ]

[22] Moore, D.M., Reynolds, R.C., 1997. X‐Ray Diffraction and the Identification and Analysis of Clay Minerals, 2nd ed. Oxford University Press, New York. [ Links ]

[23] Mubarak, I., Angulo-Jaramillo, R., Mailhol, J. C., Ruelle, P., Khaledian, M., Vauclin, M. Spatial analysis of soil surface hydraulic properties: Is infiltration method dependent? Agricultural Water Management, 97 (10), pp. 1517-1526, 2010. https://doi.org/10.1016/J.AGWAT.2010.05.005Links ]

[24] Mügler, C., Ribolzi, O., Janeau, J.-L., Rochelle-Newall, E., Latsachack, K., Thammahacksa, C., Viguier, M., Jardé, E., Henri-Des-Tureaux, T., Sengtaheuanghoung, O., Valentin, C. Experimental and modelling evidence of short-term effect of raindrop impact on hydraulic conductivity and overland flow intensity. Journal of Hydrololy. 570, pp. 401-410, 2019. https://doi.org/10.1016/J.JHYDROL.2018.12.046Links ]

[25] Oliveira Júnior, J. A. S., Souza, E. S., Correa, M. M., Lima, J. R. de S., Souza, R. M. S., Silva Filho, L. A. Variabilidade espacial de propriedades hidrodinâmicas de um Neossolo Regolítico sob pastagem e caatinga. Revista Brasileira de Engenharia Agrícola e Ambiental, 18 (6), pp. 631-639, 2014. https://doi.org/10.1590/S1415-43662014000600010Links ]

[26] Peixoto, D. S., Silva, L. C. M., Melo, L. B. B., Azevedo, R. P., Araújo, B. C. L., Carvalho, T. S., Moreira, S. G., Curi, N., Silva, B. M. Occasional tillage in no-tillage systems: A global meta-analysis. Science of the Total Environment, 745 (140887), pp. 1-14, 2020. https://doi.org/10.1016/j.scitotenv.2020.140887Links ]

[27] Prima, S. D., Stewart, R. D., Castellini, M., Bagarello, V., Najm, M. R. A., Pirastru, M., Giadrossich, F., Lovino, M., Angulo-Jamarillo, R., Lassabatere. Estimating the macroscopic capillary length from Beerkan infiltration experiments and its impact on saturated soil hydraulic conductivity predictions. Journal of Hydrology, 589 (125159), pp. 1-11, 2020. https://doi.org/10.1016/j.jhydrol.2020.125159Links ]

[28] Robertson, G.P. GS+: Geostatistics for the Environmental Sciences; Gamma Design Software: Plainwell, MI, USA, 2008. [ Links ]

[29] Santos, C. A. G., Silva, J. F. C. B. C., Silva, R. M. Caracterização Hidrodinâmica dos Solos da Bacia Experimental do Riacho Guaraíra Utilizando o Método Beerkan. Revista Brasileira de Recursos Hídricos, 17 (4), pp. 149-160, 2012. https://doi.org/10.21168/rbrh.v17n4.p149-160Links ]

[30] Santos, D. P., Montenegro, A. A. A., Rodrigues, R. A. S., Araújo, D. C. S., Santos, C. S, Cruz Neto, J. F. Variabilidade espacial de atributos físicos do solo em Vale Aluvial na região de Pernambuco, 12 (1), pp. 2271-2282, 2018. https://doi.org/10.7127/rbai.v12n100647Links ]

[31] Souza, E. S., Antonino, A. C. D., Angulo-Jaramillo, R., Netto, A. M. Caracterização hidrodinâmica de solos: Aplicação do método Beerkan, Revista Brasileira de Engenharia Agrícola e Ambiental, 12 (2), pp. 128-135, 2008. https://doi.org/10.1590/S1415-43662008000200004Links ]

[32] Souza, R., Souza, E., Netto, A. M., Almeida, A. Q. de, Júnior, G. B., Silva, J. R. I., Lima, J. R. de S., Antonino, A. C. D. Assessment of the physical quality of a Fluvisol in the Brazilian semiarid region. Geoderma Regional 10, pp. 175-182, 2017. https://doi.org/10.1016/J.GEODRS.2017.07.008Links ]

[33] Souza, R., Hartzell, S., Ferraz, A. P. F., Almeida, A. Q., Lima, J. R. S., Antonino, A. C. D., Souza, E. S. Dynamics of soil penetration resistance in water-controlled environments, Soil & Tillage Research 205 (104768), pp. 1-8, 2021. https://doi.org/10.1016/j.still.2020.104768Links ]

[34] Teixeira, P. C., Donagemma, G. K., Fontana, A., Teixeira, W. G., 2017. Manual de Métodos de Análise de Solo, 3rd ed. Embrapa Informação Tecnológica, Brasília. [ Links ]

[35] Vauclin, M., Vieira, S. R., Vachaud, G., Nielsen, D. R. The Use of Cokriging with Limited Field Soil Observations1. Soil Science Society of America Journal 47 (2), pp. 175-184, 1983. https://doi.org/10.2136/sssaj1983.03615995004700020001xLinks ]

[36] Vieira, S. R., Hatfield, J. L., Nielsen, D. R., Biggar, J. W. Geostatistical theory and application to variability of some agronomical properties. Hilgardia, 51, pp. 1-75, 1983. https://doi.org/10.3733/hilg.v51n03p075Links ]

[37] Warrick, A. W., Nielsen, D. R., 1980. Spatial variability of soil physical properties in the field, in: Hillel, D. (Ed.), Applications of Soil Physics. Elsevier, New York, pp. 319-344. https://doi.org/10.1016/B978-0-12-348580-9.50018-3Links ]

[38] Yamamoto, J. K., Landim, P. M. B., 2013. Geoestatística conceitos e aplicações. Editora Oficina de Textos, São Paulo. [ Links ]

A.G. Pinheiro, is a BSc. in Technologist in Irrigation and Drainage of the Federal Institute of Ceará (IFCE), Brazil. MSc. and PhD. student in Agricultural Engineering, at the Federal Rural University of Pernambuco (UFRPE), Brazil. Research line: irrigation and drainage; soil and water management and conservation; hydrology; agricultural resilience practices in the semi-arid.ORCID: 0000-0002-5912-6906

A.M.R.F. Jardim, received the BSc. in Agronomy and MSc. in Crop Production from the Federal Rural University of Pernambuco, Brazil. Currently a doctoral student in Agricultural Engineering at the Federal Rural University of Pernambuco, Brazil. ORCID: 0000-0001-7094-3635.

A.A.A. Montenegro, is a full-time professor at the Federal Rural University of Pernambuco, Brazil. He is a BSc. Eng. in Civil Engineer from the Federal University of Pernambuco, a MSc. in Hydraulics and Sanitation from the University of São Paulo, Brazil, and a PhD. in Water Resources from the University of Newcastle Upon Tyne. He has a postdoctoral degree in hydrological modeling from the Center for Ecology and Hydrology of Wallingford, England. He has experience in the area of Water Resources and Agricultural and Environmental Engineering.ORCID: 0000-0002-5746-8574

T.G.F. Silva, is a full-time professor at the Federal Rural University of Pernambuco, Brazil. MSc and Ph.D. in Agricultural Meteorology. Current President of the Brazilian Society of Agrometeorology-SBAgro. ORCID: 0000-0002-8355-4935.

J.R.I. Silva, received the BSc. in Agronomy and MSc. in Crop Production from the Federal Rural University of Pernambuco, Brazil. Currently a doctoral student in Agricultural Engineering from the Federal Rural University of Pernambuco, Brazil. ORCID: 0000-0003-04883-0514.

How to cite: Gebson Pinheiro, A., Maniçoba da Rosa Ferraz Jardim, A., Antônio de Assunção Montenegro, A., Freire da Silva, T.G. and Inácio Silva, J.R., Characterization of alluvial soil hydrodynamics in the Upper Ipanema River Basin using the Beerkan Method,. DYNA, 88(218), pp. 178-184, July - September, 2021.

Received: October 20, 2020; Revised: June 28, 2021; Accepted: July 15, 2021

Creative Commons License The author; licensee Universidad Nacional de Colombia.