Keywords
Field capacity, maximum retention potential, precipitation-runoff relationship, curve number method
This article is included in the Agriculture, Food and Nutrition gateway.
Field capacity, maximum retention potential, precipitation-runoff relationship, curve number method
The most important natural resources are soil (Bautista et al., 2018; Mátyás et al., 2018; Singla et al., 2018; Sukiasyan & Kirakosyan, 2020), air (DRC, 2020), and water (Bersosa et al., 2019). Therefore, protection of natural resources and their rational use are essential to provide sustainable environmental and agricultural management practices (Mátyás et al., 2020; Sándor et al., 2020; Shomana et al., 2020). The hydrological cycle is the process of water circulation in the earth, therefore, it is important to study surface and underground water, to understand the cycle in a global way and its interrelation with the environment (Pacheco Moya et al., 2018).
Also, the management of a basin must be oriented to the rational use of its resources, i.e., to order and spatially distribute their use (Sánchez-Cohen et al., 2015). Management is important, as a basin has environmental, social, economic, and cultural implications in the distribution of water resources (Sánchez-Cohen et al., 2015). As of today, in Ecuador, most of the hydrographic basins do not have proper conservation management, causing alterations in the hydrological cycle and generating for a given area little retention of humidity, erosion, and flooding over periods of rain and water scarcity in dry periods (Lizárraga-Mendiola et al., 2019).
Surface runoff is part of the hydrological cycle; it occurs as a result of excess water, which that does not filter off, and collects in the secondary channels, and finally ends up in the main channel, until it reaches the mouth of the river (Francés & Montalvo, 2017). When this process is not adequately controlled, it becomes a critical environmental problem, as it favors water erosion of the soil, rivers overflow, increasing the sediment load, and, by dragging, it causes water contamination (Alonso-Brito, 2016). Various study techniques have been developed, both empirically and in the field (Maidment et al., 1994), which allow determination of the volume of runoff in a given area, for providing adequate water management (Amy, 1974). These techniques enable appropriate design structures for the use and control of water resources. Surface and underground runoffs depend on the hydraulic characteristics of the subsoil, such as hydraulic conductance and porosity, which are distributed according to the geological and topographic soil conditions (Lee et al., 2015).
As a result, runoff calculations can be derived from precipitation events, using correlations between the characteristics of the soil and its vegetation, and also taking into account the interaction of different components of the natural system (Pacheco Moya et al., 2018). Nevertheless, the hydro-meteorological information in the hydrographic basins is limited in terms of a rainfall and river database of several years, as there is no good hydro-meteorological network in Ecuador (Ochoa-Tocachi et al., 2018), and this data is needed for hydrological modeling of water balance, when calculating water supply and demand (Zubieta et al., 2018).
Therefore, one can apply empirical models that provide approximate data on the hydrological reality of an area, e.g., for runoff the Curve Number (CN) method is used (Uwizeyimana et al., 2019).
The goal of this work is to determine the field capacity, considering the physical characteristics of the environment and the use of the soil and vegetation cover. We aim to obtain quantitative values for CN, and to determine the maximum potential retention and the precipitation-runoff relationship. These calculations are done for the Chibunga and Guano rivers micro-basins, located in the province of Chimborazo, Ecuador. This work should be helpful in quantifying the water balance in a basin, and to determine water supply and demand in a more efficient manner.
Our study was conducted from January 7 to March 29, 2019, in the area of interest, namely, the micro-basins of the Chibunga and Guano rivers, which are part of the sub-basin of the Chambo River, in the province of Chimborazo. Economic activity, which predominates in the region is agriculture, providing living means for most families residing there. Chibunga river originates from the slopes of the Chimborazo Volcano and descends through the Arenal moorland, until it reaches agricultural areas in San Juan town, and the communities including Chimborazo, Shobol Llinllin, with the name arising from Rio Chimborazo. Next, it joins the Cajabamba river at 3,238 m above sea level, and becomes Chibunga river (Mayorga & Carbonel, 2018). This micro-basin extends over 38 km, extends from northwest to southeast, and is the main tributary of the Chambo River at the Northern extremity of the Amazon region of Ecuador (at 680 km from the Amazonas River) (CNRH, 2007). The Guano River micro-basin is born on the slopes of the Chimborazo Volcano and receives contributions to its channel from streams, such as Cascajal, Chuquipogio, Abras, Puluchaca, Patulú, Igualata, Asaco, and others (Cevallos Gaibor, 2015). Its approximate area is 37,061 ha, and most of its water is used in the Guano Canton (Mendoza, 2016) until it eventually flows into the Chambo River.
Sampling points were identified by the over position of geology and elevation layer, where the most representative categories concerning land use were identified: Pa (moorland), Hu (wetlands), Bs (forest), Na (snow and water), Z (urban area), C (crops), Se (eroded soil), and Ps (grass).
Subsequently, the following criteria were considered to determine the sampling points:
The CN parameter measured in the field (in-situ) and laboratory. First, the infiltration capacity was determined by infiltrometer (model 09.04, SDEC France, France). Second, soil characterization was done for permeability, texture, porosity, and organic matter (Tailor & Shrimali, 2016).
At each sampling point, infiltration tests were performed with by means of a double ring infiltrometer (Model 09.04, SDEC France, France). Briefly, two concentric rings were set up on a portion of soil and water was poured into the outer ring. Variation of the water level in the inner ring was measured to determine infiltration of the soil (Carreras Nampulá et al., 2015). Samples should be as little disturbed as possible, for not altering the results.
Figure 1 shows a total of 23 points corresponding to the Chibunga river micro-basin and 21 points for the Guano river micro-basin, which were used for sampling. Five soil samples were taken from each sampling point, with the excavation being every 0.40 m to a maximum depth of 2.0 m, within a circular area 0.50 m in diameter. A total of 2 kg of soil were taken by a tubular soil sampler (15”L x ¾”Dia, Accuproducts International, USA) for each point, were placed in sampling bags (125 mm × 225 mm, Geology Superstore, United Kingdom), and transported to the Environmental Services Laboratory to determine the permeability, texture, porosity, and organic matter.
The red pins correspond to the Chibunga and Guano river micro-basins. Figure 1 was adapted from satellite images to identify the categories of land use and vegetation cover. These images were downloaded from the Copernicus website, using the Sentinel 2A satellite (Immitzer et al., 2016) (data is open source once a user registers for access).The images are from 2013, 2014, and 2018.
Soil texture was determined according to the USDA triangle of textural classes (Bautista Carrascosa & Oliver Talens, 2017), applying the organoleptic method. Briefly, numerical values were assigned within the range of 1–12 to define granulometry: lower numbers (1–5), represent fine granulometry, such as clay and loam clay type; mid numbers (6–9) correspond to medium granulometry, such as loams and silts soils; while high numbers (10–12), are for coarse granulometry, such as sands.
Porosity was obtained by a literature method (Liu et al., 2017), which relies on the real and apparent density of the soil. Calculations were based on Equation (1), which can be applied to each soil class, identified in the study area.
Where:
P = Total porosity of the soil sample, [%]
Da = Apparent density of the soil, [g/cm3]
Dr = Real particle density, [g/cm3]
Permeability or hydraulic conductivity (Papierowska et al., 2018) was calculated using a permeameter (with internal diameter D = 6.4 cm, height L = 15.3 cm) (Cotecno Falling Head Permeameter, Santiago, Chile), in which h values are measured, obtained for various times elapsed from the beginning of the test. Then, the soil permeability coefficient, k, is calculated using Darcy's Law, according to Equation (2) and Equation (3), (Wang et al., 2020):
Where: Q = flow, expressed in m3/s; L = length of sample (m); k = Darcy's permeability coefficient, variable as a function of the sample material, expressed in m/s; A = cross-sectional area of the sample (m2); H3 = height, above the water level in reference tube placed at the entrance of the filter layer; H4 =height, above the reference level, that reaches the water in a tube placed at the out of the filter layer.
Organic matter was determined via calcination of the samples to later calculate the weight loss, which corresponds to the organic compounds present in the samples (Papierowska et al., 2018). The following procedure was used: crucibles with samples were measured in triplicate. Samples were subjected to 105 °C for 2 h; then, the temperature was raised to 550 °C for 2h. At the end, samples were cooled for 2 h inside the muffle. Subsequently, samples are returned to the oven at 105 °C for 30 min, to stabilize their temperature and eliminate moisture remained in the samples, then cooled in the desiccator, and finally weighed for determining the weight loss of each sample. Soil organic matter (SOM) was determined gravimetrically and calculated according to (Equation 4).
CN was determined according to the SCS method, a technique developed for infiltration proofs depending on the runoff generating properties. The procedure considered: the hydrological group of soil (HGS), the previous humidity condition, the use of the land and treatment within the hydrographic basin. A value is calculated, which estimates the soil condition with regard to its usefulness (Satheeshkumar et al., 2017). CN is a curve number, obtained from cross-tabulation in spatial digital format, processed by soil and geomorphology maps, generating the hydrological groups (GH), as described in the SCS Handbook of Hydrology (NEH-4), Section-4 (USDA 1972). CN can be expressed by Equation (5).
where:
Taking as the input values of the CN parameter, the field capacity or maximum retention potential (S) was determined as the relatively constant amount of water contained by a saturated soil, after 48 h of drainage (Wróbel & Boczoń, 2020). Hence, there is a relationship between the soil and the coverage conditions within the micro-watersheds, via CN values. Numbers are shown on a scale from 0 to 100; the greater the number, the greater the direct runoff volume that can be expected from a storm (Lian et al., 2020).
Precipitation values were obtained from the interpolation by the kriging method, by using data from 6 meteorological stations (Table 1), over the time encompassing years 2010–2015. Precipitation values can be calculated using Equation (7) (Naranjo, 2014)
Z (si) = the value measured at location
λi = an unknown weight for the value measured at location
S0 = the location of the prediction
N = the number of measured values
Similarly, for runoff, the maximum retention potential was analyzed as a function of yearly precipitation. By this, one obtained the total flow, based on the vegetation cover (Sikorska et al., 2018).
For the multitemporal study the following methodology was used, all methods described below being applied for each sampling site:
Supervised Classification Criteria: Study areas should match the number of categories that are intended to be analyzed (TELEDET, 2017). The areas must be correctly identified to cover all classes and must be homogeneous. Also, in supervised classification uniband statistical analysis verifies that values resemble normal distribution. To meet these criteria, we used the following classic algorithms: classifier by minimum distance, by parallelepipeds, and by maximum probability (Acosta, 2017).
With this classifier collected data is used to determine the mean of the classes selected at the study sites. This approach assigns each unidentified pixel to the class with the closest mean. For the closest mean values one can use the Euclidean distance between the pixel and the center of each class. Although simple from the computational standpoint, this algorithm has certain limitations, including the insensitivity to different degrees of variance in the spectral responses (TELEDET, 2017).
ML or maximum probability assumes that the reflectivity values in each class follow a multivariate normal probability distribution. The vector of means and the variance-covariance matrix are then used to estimate the probability by which one given pixel belongs to all classes. Finally, pixels were assigned to the class to which it belonged with the highest probability (Tso & Mather, 2009). A Bayesian approach was used and considered was the a priori probability, whether a pixel belonged to a certain class. The greater the surface area, the greater the probability that a pixel belongs to it (Toro et al., 2013).
It consists of a double-entry table, which compares the real values with results obtained in the classification. The diagonal of the matrix reveals the number of real pixels and the classification that coincides by category, while the remaining pixels indicate the ones that may be confused with other categories. One represents vertically the percentage of real pixels that were confused in the truth-terrain, while horizontally one shows the percentage of all pixels (Borràs et al., 2017).
The accuracy measures considered by us are simple to use, being based on either the main diagonal, or the rows and columns of the contingency table. In some cases, a totally random distribution of pixels in the classes can yield apparently correct results in the contingency table (TELEDET, 2017).
Through the multi-temporal study and the application of spectral signatures, eight categories of vegetation cover were identified for each micro-basin: Pa (moorland), Hu (wetlands), Bs (forest), Na (snow and water), Z (urban area), C (crops), Se (eroded soil), and Ps (grass). The 26 sampling points were in the Chibunga river micro-basin and 23 points in the Guano river micro-basin.
Physical characteristics of the soil from each sampling point were determined, including porosity, permeability, and infiltration capacity (Table 2). These unique values are required to determine vegetation cover and soil texture, necessary to calculate CN.
Table 3 lists the values for CN and HSG (hydrological group) depending on the state of conservation and conditions present in the sampling area.
In Table 3 one finds different combinations of vegetation cover and the HSGs of the soil in the micro-basin of the Chibunga River, where each combination corresponds to a CN that represents the runoff potential, in which the highest values belong to the wetlands, urban area, snow, and water, assigning them the value of 100. In these areas rain is converted directly to runoff, as they are impervious areas. By contrast, the Guano River micro-basin revealed that the moorland class had clay soil with CN of 78, the soil texture for eucalyptus forest was silty clay, and its CN was the lowest of the data set (46).
In Figure 2, the “S” values (Lian et al., 2020) are shown for the Chibunga River micro-basin for years 2013 and 2014. One can observe a low CN value (orange color) in the range of 1–52 mm of water retained in the soil. Highest values of the water retention potential in the study area are highlighted in yellow; they correspond to moorlands and coniferous forests. In the micro-basin of the Guano River, we uncovered water retention potential values of 0 mm (red color). Such areas are characterized by the presence of urbanization, bodies of water, and snow. On the opposite end, highest values of water retention correspond to areas with eucalyptus forests (blue color), which are located in the middle part of the micro-basin. One should emphasize that the vegetation cover in this area has not changed over the three years when the multitemporal study was conducted.
In Figure 3, the variation of precipitation is shown over the 3 years of our multitemporal study. Our findings agree with a prior study of Chidichimo and co-workers, according to which highest precipitation occurs over the months of March, April, and May, while minimum precipitation values are recorded in January, November, and December (Chidichimo et al., 2018).
The values of precipitation-runoff in the micro-basins were obtained via the curve number method (CN). Values of maximum retention potential were found to the vegetation cover for years 2013, 2014, and 2018.
In Figure 4, the precipitation-runoff relationship is much higher for an urban area in 2018 with a similar value of 44.6 mm, and a CN of 100, due to its retention capacity which is 0 . On the other hand, we obtained for vegetable crops and eucalyptus forests a similar value of 13 mm of runoff in 2014 and the retention capacity was of 60 and 68 mm of water in the soil, respectively. Therefore, it was determined that the grass category in 2013 had the lowest runoff at 4.4 mm and a CN of 79, demonstrating the saturation of porous medium in silty clay soil.
Figure 5 reveals that in 2013, vegetable crops have a runoff of 12.7 mm and an average annual rainfall of 32.5 mm, while in 2014 and 2018, the urban area, wetlands, and eroded soil was in the range of 27–44 mm of runoff, showing retention of water of 0 mm in the soil. In addition, the greatest values of precipitation were found in 2018 for moorland category with 59.2 mm and runoff of 21.7 mm in the precipitation-runoff conversion, while the retention value is 72 mm of water, indicating the saturation of the porous medium in the soil.
The maximum retention potential is directly affected by the CN values that are recorded according to the conditions and use of the vegetation cover on the land. The categories snow, and water, wetlands, and urban area showed null retention of water in the soil because their transformation is directly from precipitation to runoff. This trend was opposite for the other categories, with the highest values in the maximum retention potential, while their values of runoff decreased.
To better understand the water storage capacity, we had to adequately interpret the values of field capacity. This was done by relating the water storage capacity to permeability, porosity, and the amount of organic matter. Obtained values allowed the formulation of adequate recovery or conservation measures to achieve better retention of water in the soil.
Figshare: GPS Coordinates of Sampling Points in Guano and Chibunga basins, https://doi.org/10.6084/m9.figshare.13874696.v1 (Mendoza et al., 2021).
Figshare: Raw data Chibunga, https://doi.org/10.6084/m9.figshare.13875050.v1 (Melendez, 2021).
Figshare: Supporting data: Guano (Recuperado).xlsx, https://doi.org/10.6084/m9.figshare.13299686.v1 (Dama Research Center Limited, 2020).
Data are available under the terms of the Creative Commons Attribution 4.0 International license (CC-BY 4.0).
Views | Downloads | |
---|---|---|
F1000Research | - | - |
PubMed Central
Data from PMC are received and updated monthly.
|
- | - |
Is the work clearly and accurately presented and does it cite the current literature?
Partly
Is the study design appropriate and is the work technically sound?
Yes
Are sufficient details of methods and analysis provided to allow replication by others?
Yes
If applicable, is the statistical analysis and its interpretation appropriate?
Yes
Are all the source data underlying the results available to ensure full reproducibility?
Yes
Are the conclusions drawn adequately supported by the results?
Partly
Competing Interests: No competing interests were disclosed.
Is the work clearly and accurately presented and does it cite the current literature?
Yes
Is the study design appropriate and is the work technically sound?
Yes
Are sufficient details of methods and analysis provided to allow replication by others?
Yes
If applicable, is the statistical analysis and its interpretation appropriate?
Yes
Are all the source data underlying the results available to ensure full reproducibility?
Yes
Are the conclusions drawn adequately supported by the results?
Yes
Competing Interests: No competing interests were disclosed.
Reviewer Expertise: Hydraulic and environmental engineering
Alongside their report, reviewers assign a status to the article:
Invited Reviewers | ||
---|---|---|
1 | 2 | |
Version 1 03 Mar 21 |
read | read |
Provide sufficient details of any financial or non-financial competing interests to enable users to assess whether your comments might lead a reasonable person to question your impartiality. Consider the following examples, but note that this is not an exhaustive list:
Sign up for content alerts and receive a weekly or monthly email with all newly published articles
Already registered? Sign in
The email address should be the one you originally registered with F1000.
You registered with F1000 via Google, so we cannot reset your password.
To sign in, please click here.
If you still need help with your Google account password, please click here.
You registered with F1000 via Facebook, so we cannot reset your password.
To sign in, please click here.
If you still need help with your Facebook account password, please click here.
If your email address is registered with us, we will email you instructions to reset your password.
If you think you should have received this email but it has not arrived, please check your spam filters and/or contact for further assistance.
Comments on this article Comments (0)