Next Article in Journal
Humidification–Dehumidification (HDH) Desalination and Other Volume Reduction Techniques for Produced Water Treatment
Previous Article in Journal
Analysis of Small and Medium–Scale River Flood Risk in Case of Exceeding Control Standard Floods Using Hydraulic Model
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Influence of Multi-Layered Structure of Vadose Zone on Ecological Effect of Groundwater in Arid Area: A Case Study of Shiyang River Basin, Northwest China

1
Institute of Hydrogeology and Environmental Geology, Chinese Academy of Geological Sciences, Shijiazhuang 050061, China
2
China University of Geosciences (Beijing), Beijing 100083, China
3
Key Laboratory of Groundwater Sciences and Engineering, Ministry of Natural Resources, Shijiazhuang 050061, China
4
School of Water Resources and Environment, Hebei University of Geosciences, Shijiazhuang 050031, China
*
Author to whom correspondence should be addressed.
Water 2022, 14(1), 59; https://doi.org/10.3390/w14010059
Submission received: 5 November 2021 / Revised: 24 December 2021 / Accepted: 24 December 2021 / Published: 28 December 2021
(This article belongs to the Section Hydrogeology)

Abstract

:
The natural vegetation in arid areas of northwest China is strongly dependent on the availability of groundwater. Significantly, capillary water plays an essential role in regulating the ecological groundwater level in the multilayered structure of the vadose zone. The soil-column test and field survey in the lower reaches of the Shiyang River Basin were conducted to investigate the influence of the multi-layered structure of the vadose zone on maintaining the ecological effect of groundwater. Based on the field survey, the results show that the depth of groundwater is 3.0 m, and the rising height of capillary water is 140 cm. In the soil-column test, the height of the wetting front of the column was 125 cm. During the water releasing test, the water held by the vadose zone was 182.54 mm, which would have maintained Haloxylon’s survival in a growing season. Therefore, the multi-layered structure of the vadose zone extends the ecological groundwater depth and consequently enhances the ecological function of groundwater. Importantly, with a lower groundwater level, the clay soil layer within the rising height range of the original capillary water would hold more water and maintain a higher water content for a certain period to supply surface vegetation.

1. Introduction

Inland basins of northwest China have a typical continental climate, which is characterized by low precipitation and high evaporation, as well as scarce water resources. Therefore, the natural vegetation of inland basins is highly dependent on groundwater [1,2]. Groundwater depth controls water availability to vegetation and is essential for conserving groundwater-dependent ecosystems (GDEs). Once the depth of groundwater table is lower than ecological groundwater level for a long time, the natural vegetation ecosystem will lose the water supplied by groundwater through the supporting capillary zone [3]. What is worse, it will lead to natural vegetation degradation or even death and desertification [4]. Since the 1950s, because of the rapid development of the economy and society, groundwater has been over-exploited, which led to the continuous decrease of groundwater depth, exceeding the ecological groundwater level of natural vegetation in most inland basins. Hence, how to exploit and utilize groundwater reasonably and protect ecology effectively is an urgent problem to be solved in northwest China.
The ecological groundwater level has attracted worldwide attention because of the close relationship between groundwater and groundwater-dependent ecosystems (GDEs) [5,6,7,8,9,10]. Many concepts have been proposed to define the reasonable ranges of groundwater depth for sustaining the GDEs [11], such as critical groundwater depth [12], ecological groundwater table [13], depth-to-groundwater threshold [14], and critical groundwater level [15]. Although the literal expressions are different, these statements have the same essence concerning ecological groundwater level [16]. The influence of the ecological groundwater level on different vegetation types is different. When the groundwater depth is less than 5 m, the diversity and richness of herbage species are significantly higher than that of shrub layer, while the diversity and richness of herbage species are lower than that of shrub layer when the groundwater depth is more than 5 m. The ecological groundwater level is not a single value but a series of values. The threshold of ecological groundwater level suitable for different vegetation growth has been put forward by various methods [5]. For example, the interactions among groundwater, vadose zone and vegetation are discussed based on correlation analysis between the soil water content, groundwater chemical property, plant growth and the groundwater depth [17,18], or by establishing distributed eco-hydrological model and HYDRUS solute transport model [19,20]. No matter what method is used to calculate the ecological groundwater level, the extreme ecological groundwater level should be less than the sum of the height of supporting capillary water and the depth of vegetation root system (Figure 1). As long as the root system of vegetation reaches the depth of groundwater rising through capillarity, it can continue to obtain groundwater recharge.
Different vegetation types have different depths of root development. Zhang et al. selected typical vegetation in an arid area of northwest China and statistically analyzed the average root development depth of tree, shrub and herb; they are 2.7 m, 2.2 m and 1.8 m, respectively [21]. Sheng et al. [22] studied the root system of Haloxylon ammodendron at different ages, and the average root length is 1.4 m (young seedling, 1 year), 2.5 m (young, 4 years), 8.0 m (middle age, 20 years) and 13.5 m (40 years) separately. In addition, the root system of vegetation grows with the increase of groundwater depth to absorb the deep groundwater and keep growing. The vadose zone is a key factor connecting groundwater and surface vegetation [23]. A different structure of the vadose zone supports a different rising height of capillary water, which affects the ecological groundwater level and leads to a different ecological response of vegetation on land surface [24,25,26,27]. However, when the groundwater level drops, the release of water from pore channels with different sizes in the soil is not synchronous, which results in the delay of gravity drainage. Different vadose zones have different water-holding capacities and delayed water-releasing capacities, which lead to different effects on the growth of natural vegetation [28,29].
Scholars have made outstanding contributions to the research, as mentioned above. However, when the groundwater level drops, can the water held by the vadose zone continue to be absorbed by the roots of the vegetation? And if so, for how long? The issue is rarely discussed, but it is essential for the conservation of natural vegetation. Moreover, the results can provide more accurate data to optimize the period and volume of ecological water conveyance, which positively affects elevating the groundwater table in arid areas [30,31]. Given the above problems, this paper aims to study the influence of the multi-layered structure of the vadose zone on the ecological effect of groundwater: (1) to explore the lithological structure of the vadose zone in the typical natural vegetation area; (2) based on the exploration results, to survey pit profile, including the field test of groundwater capillary rise height; (3) to discuss the availability of water-holding capacity in the vadose zone on vegetation root water uptake by a soil-column test.

2. Materials and Methods

2.1. Study Area

The Shiyang River Basin is located in the east of Hexi Corridor in northwest China, Gansu Province, which is one of the driest regions in China. In Shiyang River Basin, the upper reaches are the Qilian mountains, the middle reaches are the corridor plains, and the lower reaches are the oasis of Minqin (Figure 2). The ecological environment of Shiyang River Basin is extremely fragile, and desertification is becoming more and more serious [32]. Minqin basin, which lies in the lower reaches of the Shiyang River Basin, is in danger of becoming the second “Lop Nur” surrounded by the Tengri and Badain Jaran deserts. It is a typical desert island, known as a “Desert Oasis”. The basin geographical boundary lies between 103°02′ E to 104°02′ E longitude and 38°05′ N to 39°06′ N latitudes, and total area of the basin is 1.59 × 104 km2, in which desert and desertification area accounted for 94%; oasis area only accounted for 6%. The studied catchment experiences an arid-continental climate with an annual average temperature of 8.8 °C and 3142.2 h of sunshine. The mean annual evaporation in the watershed is 2675.6 mm, far exceeding the mean annual precipitation (113.2 mm) [33].
Note: The image is based on Landsat 8 data of 2019 and has been processed by radiation correction, atmospheric correction, geometric correction, image fusion and color synthesis. In the image, the green color indicates the vegetation area, including natural vegetation and farmland; the brown color indicates desert areas; the dark color indicates wetlands; the red box area is the study area; and the black oval area indicates the different positions of the upper, middle and lower reaches of the basin.

2.2. Materials and Methods

The field survey and soil-pit sampling (field survey) are more in line with the actual situation, but the experimental conditions cannot be controlled. While the soil-column test can control the test conditions according to the experimental purpose, the results are too ideal. Therefore, in this paper, field survey and soil-column tests were designed according to the purpose of the research, and the results of the two were compared and verified with each other. Based on the typical stratum-structure investigation, a specific area was selected to conduct the field survey, study the groundwater capillary rise height and retrieve soil samples from different stratums for soil-column tests. The result obtained by the field survey was the closest to the actual value and can be used as an accurate measure of the soil-column tests. The soil-column tests included the capillary-water rising test and the water-releasing test of the multi-layered structure of the vadose zone. The capillary-water rising test, which is the basis of the water-releasing test, aimed to study the height of groundwater capillary rise of the typical vadose zone structure, the process of water migration, and the water content of different lithological soils. The purpose of the water-releasing test was to study: (1) whether the amount of water held by different lithological soils in the vadose zone can continue to maintain the survival of the vegetation when the depth of the local groundwater drops significantly below the ecological groundwater level, and (2) how long it will last.

2.2.1. Field Survey

In order to study the capillary rising height of the typical stratum structure in the natural vegetation area, the phreatic water-depth and stratum-structure drilling survey were performed in the study area in the early stage. Two survey routes were established in the natural vegetated area, and each was approximately 10 km. At each survey line, drilling was carried out at intervals of approximately 1 to 2 km depending on the growth of vegetation and the change in vegetation type (Figure 3). Soil samples were taken every 20 cm to record soil lithology. All boreholes were manually drilled with an earth-boring auger until the groundwater table and each soil lithology were analyzed for particle composition in the laboratory.
Based on the analysis and comparison of the boreholes, the pit site was chosen to be located in the natural vegetation area about 2 km southeast of the fifth team of Fengqing village in Minqin Basin, which is between farmland and the Tengger Desert. The pit was rectangular, digging from the ground surface to the phreatic water, of which the depth was 3.0 m. The surrounding vegetation was mainly composed of Tamarix ramosissima, Salt claw and Black fruit Chinese wolfberry, while the pit was located near Tamarix chinensis.
The field survey was used to measure the rising height of capillary water in groundwater directly. In the excavation process, the soils of different depths and lithology were arranged in order, and about 100 kg of silt, clay, silty fine sand and medium-fine sand were taken back to prepare the soil-column experiment. After excavating the test pit, the field lithology was described first. The 20–30 g soil samples were taken from the section at intervals of 10 cm depths and placed in the aluminum box, and then 3 repeated samples were taken from each layer. The soil moisture content was measured by drying method in the laboratory. Then soil samples were taken from the section with a 100 cm3 soil sample ring every 20 cm depth to determine its dry bulk density. After collecting soil samples, the observation profile was covered with a shading net to avoid sun exposure and to ensure soil moisture drying under natural conditions. Repeated samples were taken again from the same location in the early morning 3 days later. After 3 days of natural evaporation, the water content in the upper part of the formation will be significantly reduced when measured again. In contrast, the water content in the lower part of the formation is rarely changed due to the continuous recharge by groundwater capillary. The intersection of the two curves is where the water content is significantly reduced. The water content below the intersection is basically unchanged due to the continuous recharge of groundwater capillary action. Therefore, the intersection of the two water content curves indicates the position of the height of groundwater capillary rise [24].

2.2.2. Soil-Column Tests

In order to study the effect of the layered structure of the vadose zone on the ecological groundwater level, the laboratory soil-column tests were designed, which included the capillary water rising test and the water releasing test of the multi-layered structure of the vadose zone.
(1)
Instrument of Soil-Column Tests
The instrument, made of plexiglass, was divided into a water supply device and a test device. The water supply device, Mariotte’s bottle, can keep the constant water head during the water supply, and has a diameter of 20 cm and a height of 80 cm. The diameter of the test device is also 20 cm, and the height is 200 cm. A series of sampling holes with a diameter of 2 cm were arranged from top to bottom at intervals of 10 cm, and were used to measure the water content of soil samples during the experiment (Figure 4). Before the test, water was added to ensure the sealing effect was good. Then, soil samples were filled in layers in the test equipment. The section structure of the soil column was optimized according to the strata of the test pit. The main body of the soil column was 150 cm. It was composed of silty soil (40 cm), clay soil (40 cm), silty fine sand (20 cm) and medium fine sand (50 cm) from top to bottom. The bottom of the soil column was coarse sand (30 cm) as the aquifer (Figure 4).
After drying, crushing and sieving, the soil samples were filled in layers according to the original dry bulk density. Different lithological soil was separated using 100 mesh nylon screens to distinguish the soils after the experiments. The soil moisture test tube was installed in the center of the soil column to monitor the change of soil moisture in the soil-column section during the test. The bottom of the column was backfilled with small gravel of 20 cm thickness as a filter layer, and 30 cm-thick coarse sand was backfilled at the top of the filter as an aquifer. The Mariotte’s bottle acted as a water supply device in the capillary rising process and a drainage device in the water release process.
(2)
Capillary Water Rising Test of the Multi-layered Structure of Vadose Zone
After finishing the soil backfill, the water content of the initial soil-column profile was measured by the soil-moisture meter and corrected by the fitting formula that is listed in the following section. A hose was used to connect the water tap at the bottom of the Mariotte’s bottle to the bottom of the soil column, from which the Mariotte’s bottle supplies water to the soil column. The water level in the soil column was determined by the water head that set in the Mariotte’s bottle. When the water supply started, the coarse sand aquifer backfilled by the soil column was gradually filled, and the time started when the water level reached the set depth. During the experiment, it was necessary to record the rising position of the wetting front of the soil column and the amount of water supplied, which was expressed by the scale of Mariotte’s bottle. After the wetting front of the soil column stopped rising, the moisture content of the soil-column section was measured again by the soil-moisture meter, and the soil-moisture content was measured by the drying method with a soil sample of about 20 g packed in an aluminum box from the sampling point.
(3)
Water-Releasing Test of the Multi-layered Structure of Vadose Zone
The water-release test was the inverse process of the capillary water-rising test. The water in the Mariotte’s bottle was emptied to act as a drainage device to record the water release. The entire soil began to release water freely by gravity after opening the water tap at the bottom of the column, and at the same time started to record the time and the amount of water entering the Mariotte’s bottle. The moisture content of the whole soil-column section was measured twice a day by a soil-moisture meter, once in the morning and once in the evening. When the water in the bottle was no longer increased, the water-release test was finished.

2.2.3. Methods

(1)
Drying Method [34]
The drying method is the most commonly used standard method for measuring soil moisture. The soil samples are dried at 105 °C to a constant weight, at which time the soil organic matter does not decompose and the free water and hygroscopic water in the soil are entirely removed. Calculate the ratio of the water quality of the soil to the quality of the drying soil, that is, the mass water content, expressed as a percentage.
The weight of the sample container (aluminum box) is weighed by balance first, and about 20 g of soil sample is put into the aluminum box and weighed again. Then the aluminum box containing the soil sample is placed in an oven and dried at 105–110 °C to the constant weight (about 8 h). Finally, it is taken out and cooled in the dryer for about 20 min and weighed immediately.
ω% = (w2 − w3)/(w3 − w1) × 100
ω—mass moisture content (%); w1—the weight of aluminum box (g); w2—the weight of aluminum box and fresh soil (g); w3—the weight of aluminum box and dry soil (g).
(2)
Dry Bulk Density [34]
The soil-sample ring method is commonly used to measure soil-bulk density, and the soil-sample ring with a volume of 100 cm3 is usually used for sampling. Before sampling, the soil-sample ring should be weighed. When sampling, the edge of the soil-sample ring should be pressed down into the soil vertically until it is filled with samples. Cut the soil around the edge of the ring with a repair knife, carefully cut off the excess soil at both ends of the soil-sample ring, and wipe off the outside soil. At the same time, take another sample in the same layer with an aluminum box to measure the natural mass moisture content by drying method. Cover both ends of the soil-sample ring immediately to prevent water evaporation and then weigh and record. The drying method is used to measure mass moisture content,
γ = (x2 − x1) × (1 − ω%)/v
where γ is dry bulk density (g/cm3); ω is mass moisture content (%); x1 is the weight of soil-sample ring (g); x2 is the weight of soil-sample ring and fresh soil (g); v is volume of soil-sample ring (cm3); ω is mass moisture content (%).
The volumetric water content (θ) of the soil mentioned in this paper can be calculated from the mass moisture content and the dry bulk density of the soil.
θ = ω × γ
(3)
Data Calibration
The soil-moisture meter used to measure the moisture content of the soil column in this test is Trime-Pico-IPH TDR (Trime), which is made in IMKO Micromodultechnik GmbH (Ettlingen, Germany). Based on the time domain reflectometry with intelligent microelements, Trime can directly measure the dielectric constant of soil or other media, which is closely related to the soil moisture content. The soil moisture content can be calculated and displayed by the reading system through the output of the analog voltage. Because the calibration of Trime is carried out by the ideal medium, there are some errors for the actual testing of rock and soil, so it is necessary to use the drying method to calibrate the Trime. Based on the collected soil samples, after gradually adding water to the target soil sample and mixing different water content from low to high, the correlation between the data measured by Trime and the data measured by drying method is established. The data of the Trime test cannot be analyzed until they are calibrated by fitting formula.

3. Results

3.1. Results of the Field Survey Soil Sample Analysis and Groundwater Depths

The survey showed that the vadose zone was a multi-layered structure composed of silty soil, clay soil, silty fine sand and medium fine sand from surface to phreatic water. The composition and dry-bulk density of different lithological particles are shown in Table 1. The overall trend of groundwater depth decreases from northwest to southeast. As shown in Figure 5, the groundwater depth near farmland is about 3.5 m, and that near Tengger Desert is about 2.0 m.

3.2. Calibration Results

It can be seen from Figure 6 that the fittings of all the soil are good, among which the silty soil calibration formula is y = 1.9577x + 0.0149, R2 = 0.9698, the clay soil calibration formula is y = 1.6073x + 0.0243, R2 = 0.9536, the silty fine sand calibration formula is y = 1.984x − 2.0331, R2 = 0.9811, and the calibration formula of medium fine sand is y = 1.0885x + 0.6578, R2 = 0.9948.
The data tested by the Trime were recalibrated by the above fitting formulas and then used for the calculation in the paper.

3.3. Capillary Water Rising Height of Typical Stratum Structure in Natural Vegetation Area

By comparing the soil-moisture curves drawn from the profile, it can be concluded that the 2 curves have an obvious intersection point at a depth of 160 cm from the surface (Figure 7). The soil moisture content of the stratum above the intersection decreased obviously after 3 days of drying, while the continuous recharge moisture content of the stratum below the intersection was little changed. It can be concluded that the depth of the point was 140 cm, which indicated the capillary rising height of the groundwater. In the field investigation, a large number of calcareous nodules were found in the silty soil at a depth of 60–170 cm, the existence of which may lead to larger pores in the soil, thus blocking the capillary rise of groundwater.
As shown from Figure 6, the overall trend of soil-moisture content decreases with the increase of capillary water height in the whole range of capillary water height, which shows that the capillary-water holding capacity decreases above the water level. Only about 20 cm above the water level can maintain a higher soil moisture content. We found that 180–220 cm of the test pit profile had the lowest soil-water content, about 10%, due to the silty fine sand with loose and uneven soil particles.

3.4. Capillary Water Rising Process in Layered Vadose Zone

The soil-column test began on 4 October 2020, and ended on 17 October 2020, with a cumulative running time of 317.45 h. During the test, the cumulative water supply of the Mariotte’s bottle was 297 mm, and the cumulative height of the wetting front of the soil column was 112.6 cm (Figure 8). The total amount of water supplied to the soil column was 291.28 mm (Table 2), calculated from the soil moisture content measured by the drying method, which was basically consistent with the water supplement record of the Mariotte’s bottle.
In the range of 0–80 cm, the wetting front increased rapidly, at about 1.6 cm/h, and then the ascending speed decreased to 0.12 cm/h during the following time. According to the scatter plot of cumulative height of wetting front with the time, the relationship was in accordance with logarithmic function, y = 14.798ln(x) + 22.762, R2 = 0.9689 (Figure 8).
As shown from Figure 9, the distribution of water content is different obviously from the groundwater level to max height of the capillary rise, and the height of the saturated-water content zone maintained by groundwater level is only about 15 cm. The water content of silt and clay soil is higher than that of sand soil, and the composition and structure of soil particles directly impact the water content. From the water-content distribution, it can be seen that the capillary-rise height of groundwater is 125 cm, which is slightly higher than the height of the wetting front (112.6 cm), which is observed visually, and the observation error is inevitable.
It can be concluded that the finer-grained formation affects the capillary height and rate from Figure 10. Region ‘A’ means the wetting front moved from medium fine sand to silty fine sand. It can be seen from the red points in region ‘A’ that when the wetting front passes through the interface of medium fine sand and silty fine sand, there is an obvious stagnation and then a rapid rise. Region ‘B’ means the wetting front moved from the silty fine sand to clay soil. It can be seen from region ‘B’ that after the wetting front passes through the interface between silty sand and clay, it continues to rise for a while (yellow points) before it comes to an obvious stagnation and then rapidly rises again (red points). The two regions have a common rule in the curve, that is, the wetting front has a short delay in the transition from coarse-grained to fine-grained lithology, and the wetting-front rate increases after entering into fine-grained lithology, in which the variation of region ‘A’ is more obvious. This phenomenon shows that the different lithology interface hinders the water migration, and the more significant the difference is between the upper and lower soil lithology, the greater is the influence.

3.5. Water Release Process in Layered Vadose Zone

In order to study the influence of the vadose zone multi-layer structure on maintaining the ecological effect under the condition of groundwater drawdown in a natural vegetation area, the water-release experiment was carried out from 17 to 24 October 2020, based on the capillary-water rising test. The water in the Mariotte’s bottle was drained and used to measure the drainage during the release experiment.
Timing and recording the drainage should be started once the soil column drained water to the Mariotte’s bottle. At 21:30 on the 17th, Trime was used to test the water content of the initial soil profile, and drainage began at 23:25. At the same time, the water content of the soil-column section was measured by Trime, and the drainage process ended at 08:45 on the 24th, with a cumulative time of 9200 min. According to the curve of cumulative drainage and time, the drainage rate decreased with the increase of time and the cumulative drainage recorded by the Mariotte’s bottle was 111 mm (Figure 11).
It can be seen from Figure 12 that the water content of medium fine sand and silty fine sand changed obviously after the water release, whereas the water content of clay soil and silty soil changed inconspicuously, with the clay soil especially keeping high water content. The water content below the silty fine sand dropped significantly in 11 h (from 21:30 17th to 08:30 18th), and in the following 24 h, the decreasing rate of water content below the silty fine sand slowed. The moisture distribution of the soil profile was tested respectively again, at 09:00 on the 21st, 22nd and 24th. Since the water content distribution of three times had no noticeable change, the water-release test was finished. According to test results by Trime, the water-holding capacity was 340.21 mm (including the amount of water supplied by the Mariotte’s bottle and the previously held by the previous soil column) in the soil column before water releasing, and 305.09 mm after water releasing (Table 3). Therefore, the total water-releasing capacity was 35.12 mm (the cumulative drainage recorded by the Mariotte’s bottle was 111 mm, including 75.88 mm of water in the filter and the aquifer), in which the medium fine sand and silty fine sand layer contributed 35.06 mm, accounting for 99.83% in the drainage.

3.6. Available Water Content of Soil Column after the Water-Releasing Experiment

According to the experiment results in this paper, the water content of the whole soil column was still 305.09 mm (Table 3) according to the water content of the soil profile that was measured at 09:00 on 24 October. However, the available water content of the soil column after water release still needs to be discussed. It is well known that not all of the water in the soil can be used by vegetation, and the available water content is determined by the field capacity and the permanent wilting point [35,36,37,38]. According to the method in reference [39], when the soil pressure head is −1.5 MPa, the corresponding water content is at permanent wilting point. In this study, hydrus-1D software was used to fit the soil water characteristic curves of four different lithologies according to the neural network prediction model, as shown from Figure 13. According to the soil water characteristic curves, the corresponding permanent wilting point was determined, and the permanent wilting points of clay, silty fine sand and medium fine sand are 12.47%, 12.37, 6.51% and 4.88% (Table 4), respectively. The water content of the soil-column section after releasing water subtracts the wilting water content of the corresponding soil, and the remaining water is available for vegetation. The available water content of each testing depth in soil column after the water-releasing experiment are shown from Table 5, and the available water storage is 182.54 mm, which can be used by vegetation.

4. Discussion

4.1. Water Consumption of Typical Vegetation in Arid Area

The arid area in northwest China is characterized by sparse precipitation, large potential evapotranspiration and low soil-moisture content, so the xerophytic vegetation depends on groundwater growth. There are many studies on water consumption of typical plants in different watersheds in arid areas of northwest China. Generally speaking, the transpiration water consumption of different plants is quite different in the growing season. The transpiration of aquatic plants such as Phragmites australis and Tamarix ramosissima in the riparian zone is higher than that of xerophytes such as Apocynum venetum, Nitraria, Populus euphratica and Haloxylon. Moreover, the water consumption of the same vegetation in the growing season is different under different groundwater depths. With the increase of groundwater depth, the water available to the roots of desert plants decreases. Plants can reduce water consumption by regulating physiological characteristics such as stomatal conductance. The growth period of vegetation in arid areas is from April to September. Moreover, in the same growing season, the plant evapotranspiration can be reduced by about 50% when the suitable ecological groundwater depth decreases to the limit of ecological groundwater depth. Using the Tardieu–Davies model, Zhang et al. [40] calculated water consumption during the growing season under different ecological groundwater depths in the Shiyang River Basin. Under suitable ecological groundwater depths, the transpiration of Phragmites australis and Tamarix chinensis were 1292 mm and 1147 mm, respectively, during the growing season; and the transpiration of Haloxylon ammodendron was 279 mm; and Halocnemum, Apocynum venetum, Nitraria tangutorum and Populus euphratica were 940 mm, 913 mm, 534 mm and 448 mm, respectively. However, under the extreme ecological groundwater depth, the transpiration of all vegetation decreased to some extent. The transpiration of Phragmites australis and Tamarix chinensis were 1050 mm and 999 mm, respectively; the transpiration of Haloxylon ammodendron was 130 mm; and that of Halocnemum, Apocynum venetum, Nitraria tangutorum and Populus euphratica were 706 mm, 686 mm, 345 mm and 298 mm, respectively.

4.2. Influence of Multi-Layer Structure of Vadose Zone on Maintaining Ecological Effect of Groundwater

After the groundwater level drops, the moisture content of the soil above the water table can still maintain the vegetation to a certain extent. Taking Haloxylon ammodendron as an example, the water-supply time of the layer structure vadose zone was analyzed after the groundwater level dropped. Under the condition of suitable ecological groundwater depth, the water consumption of Haloxylon ammodendron in the growing season (from April to September) was 279 mm, which is a typical arid vegetation in the arid area of northwest China. In contrast, the water consumption of vegetation decreased by about 50% (that is 130.0 mm) under the extreme ecological groundwater depth. Therefore, according to the experimental results in this paper (Table 5), it can be concluded that under ideal conditions, the water stored in the stratum (182.54 mm) can maintain a growing season (6 months) of Haloxylon ammodendron after the groundwater level drops. In particular, the existence of a clay soil structure makes the soil profile hold higher water content after the water level drops, which has a significant positive effect on maintaining the benign development of land surface ecology.

5. Conclusions

(1)
According to the field survey, it was found that the structure of the vadose zone of natural vegetation in arid areas is a multi-layered structure, which has a positive effect on maintaining ecology. Moreover, the clay soil layer in the vadose zone will make the groundwater capillary rise higher (140 cm), making it easier for the roots of the vegetation on the ground to obtain water;
(2)
When the groundwater level drops, the water-holding capacity of the multi-layered vadose zone will continue to maintain high water content for the survival of surface vegetation in a certain time and range. Combined with the ecological water consumption of Haloxylon ammodendron in the growing season and indoor soil-column test, the available water held by multi-layered structure of the vadose zone can maintain Haloxylon ammodendron to survive in a growing season (6 months) after the groundwater level drops;
(3)
The research results can provide a scientific basis for the fine management of groundwater in arid regions. As we all know, artificial water delivery is an effective means to control ecological problems, but the period and the amount of water delivery have been debated by a wide range of scholars. The research calculated the available water held by the multi-layered vadose zone and the survival time of the typical vegetation when the groundwater depth drops below the ecological groundwater level. It will provide a scientific basis for the optimization of an ecological water-delivery period and water quantity;
(4)
Based on the above research, it is recommended to investigate the structure of the vadose zone in natural vegetation areas in arid areas to carry out the district management of ecological protection [41,42,43,44], realize refined management of groundwater and improve water-resource utilization efficiency.

Author Contributions

Conceptualization, G.Z.; methodology, H.C.; software, H.C.; validation, J.W., Q.W. and X.L.; formal analysis, H.C.; investigation, J.W., Q.W. and H.C.; data curation, G.Z.; writing—original draft preparation, H.C.; writing—review and editing, X.L.; visualization, H.C.; supervision, X.L.; project administration, G.Z.; funding acquisition, G.Z., Q.W. and H.C. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by National Key R&D Program of China (Grant No. 2017YFC0406106), the National Natural Science Foundation of China (Grant No. 41807217), and the Key Innovative Research Project of National Special Fundamental Research Funds (Grant No. JYYWF20180403).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The data presented in this study are available on request from the corresponding author.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Liu, M.; Liu, P.; Guo, Y.; Wang, Y.; Geng, X.; Nie, Z.; Yu, Y. Change-Point Analysis of Precipitation and Drought Extremes in China over the Past 50 Years. Atmosphere 2019, 11, 11. [Google Scholar] [CrossRef] [Green Version]
  2. Liu, M.; Shen, Y.; Qi, Y.; Wang, Y.; Geng, X. Changes in Precipitation and Drought Extremes over the Past Half Century in China. Atmosphere 2019, 10, 203. [Google Scholar] [CrossRef] [Green Version]
  3. Chen, M.; Zhang, Q.; Wang, Y.; Yan, L.; Deng, W. Critical depth of recharge of the vegetation by groundwater in the West Liaohe Plain. Adv. Water Sci. 2019, 30, 24–33. [Google Scholar] [CrossRef]
  4. Wang, W.; Gong, C.; Zhang, Z.; Chen, L. Research Status and Prospect of the Subsurface Hydrology and Ecological Effect in Arid Regions. Adv. Earth Sci. 2018, 33, 702–718. [Google Scholar] [CrossRef]
  5. Li, F.; Wang, Y.; Zhao, Y.; Chen, X. Modelling the response of vegetation restoration to changes in groundwater level, based on ecologically suitable groundwater depth. Hydrogeol. J. 2018, 26, 2189–2204. [Google Scholar] [CrossRef]
  6. Johansen, O.M.; Andersen, D.K.; Ejrnæs, R.; Pedersen, M.L. Relations between vegetation and water level in groundwater dependent terrestrial ecosystems (GWDTEs). Limnologica 2018, 68, 130–141. [Google Scholar] [CrossRef] [Green Version]
  7. Cao, L.; Nie, Z.; Liu, M.; Wang, L.; Wang, J.; Wang, Q. The Ecological Relationship of Groundwater–Soil–Vegetation in the Oasis–Desert Transition Zone of the Shiyang River Basin. Water 2021, 13, 1642. [Google Scholar] [CrossRef]
  8. Hu, S.; Ma, R.; Sun, Z.; Ge, M.; Zeng, L.; Huang, F.; Bu, J.; Wang, Z. Determination of the optimal ecological water conveyance volume for vegetation restoration in an arid inland river basin, northwestern China. Sci. Total Environ. 2021, 788, 147775. [Google Scholar] [CrossRef]
  9. Huang, F.; Zhang, D.; Chen, X. Vegetation Response to Groundwater Variation in Arid Environments: Visualization of Research Evolution, Synthesis of Response Types, and Estimation of Groundwater Threshold. Int. J. Environ. Res. Public Health 2019, 16, 1849. [Google Scholar] [CrossRef] [Green Version]
  10. Cheng, G.; Liu, D.; Wen, L.; Liu, H.; Dong, S.; Wang, L. Effects of Groundwater and Soil Chemical Properties on Distribution of Natural Vegetation around Playa. Arid Zone Res. 2019, 36, 85–94. [Google Scholar] [CrossRef]
  11. Jia, L.; Guo, Z.; Long, Y.; Guo, K.; Liao, Z. Rearch advances in ecological groundwater level in arid areas. Ecol. Sci. 2015, 34, 187–193. [Google Scholar] [CrossRef]
  12. Taufik, M.; Veldhuizen, A.A.; Wösten, J.H.M.; Lanen, H.A.J. Exploration of the importance of physical properties of Indonesian peatlands to assess critical groundwater table depths, associated drought and fire hazard. Geoderma 2019, 347, 160–169. [Google Scholar] [CrossRef]
  13. Zhang, C.; Shao, J.; Li, C.; Cui, Y. A study on the ecological groundwater table in the North China plain. J. Jilin Univ. Earth Sci. Ed. 2003, 33, 323–326. [Google Scholar] [CrossRef]
  14. Eamus, D.; Zolfaghar, S.; Villalobos-Vega, R.; Cleverly, J.; Huete, A. Groundwater-dependent ecosystems: Recent insights from satellite and field-based studies. Hydrol. Earth Syst. Sci. 2015, 19, 4229–4256. [Google Scholar] [CrossRef] [Green Version]
  15. Liao, Z.; Ma, Z.; Cheng, S.; Xie, X.; He, H. Dominant critical water level of groundwater and its determing method. Water Resour. Hydropower Eng. 2018, 49, 26–32. [Google Scholar] [CrossRef]
  16. Huang, F.; Zhang, Y.; Zhang, D.; Chen, X. Environmental groundwater depth for groundwater-dependent terrestrial ecosystems in arid/semiarid regions: A review. Int. J. Environ. Res. Public Health 2019, 16, 763. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  17. Zhu, L.; Xu, G.; Li, Y.; Tang, L.; Niu, Z. Relationships among plant species diversity, biomass, and the groundwater table in the Hailiutu River basin. Acta Ecol. Sin. 2017, 37, 1912–1921. [Google Scholar] [CrossRef] [Green Version]
  18. Wang, S.; Wu, B.; Yang, P.; Li, X.; Dong, X. Determination of the Ecological Groundwater Depth Considering Ecological Integrity over Oasis Irrigation Areas in the Yanqi Basin. Resour. Sci. 2011, 33, 422–430. [Google Scholar]
  19. Chui, T.F.M.; Low, S.Y.; Liong, S. An ecohydrological model for studying groundwater–vegetation interactions in wetlands. J. Hydrol. 2011, 409, 291–304. [Google Scholar] [CrossRef]
  20. Jiang, L.; Shang, S.; Mao, X. Simulation on the impact of the groundwater table on evapotranspiration of a riparian forest in arid region. J. Tsinghua Univ. (Sci. Tech.) 2013, 53, 601–605, 612. [Google Scholar] [CrossRef]
  21. Zhang, Y.; Chen, X.; Gao, M.; Liu, X. Meta-analysis of ecological depth to groundwater table and its influencing factors in arid region of northwest China. South-to-North Water Transf. Water Sci. Technol. 2020, 18, 57–65. [Google Scholar] [CrossRef]
  22. Sheng, J.; Qiao, Y.; Liu, H.; Zhai, Z.; Guo, Y. A Study on the Root System of Haloxylon Aammodendron (C.A.Mey.) Bunge. Acta Agrestia Sin. 2004, 12, 91–94. [Google Scholar] [CrossRef]
  23. Zhou, H. Review of studies on the relationship between soil water movement and energyand their driving forces in the vadose zone of arid regions. Acta Ecol. Sin. 2019, 39, 6586–6597. [Google Scholar] [CrossRef]
  24. Zhou, H.; Zhao, W. Soil physical characteristics of shallow vadose zone and modeling its effects on upward capillary rise of groundwater in an arid-desert area. Chin. J. Appl. Ecol. 2019, 30, 2999–3009. [Google Scholar] [CrossRef]
  25. Huang, M.; Barbour, S.L.; Elshorbagy, A.; Zettl, J.D.; Si, B.C. Infiltration and drainage processes in multi-layered coarse soils. Can. J. Soil Sci. 2011, 91, 169–183. [Google Scholar] [CrossRef]
  26. Zornberg, J.G.; Bouazza, A.; McCartney, J.S. Geosynthetic capillary barriers: Current state of knowledge. Geosynth. Int. 2010, 17, 273–300. [Google Scholar] [CrossRef] [Green Version]
  27. Shi, W.; Shen, B.; Wang, Z.; Wang, W. Maximum height of upward capillary water movement in layered soil. Agric. Res. Arid Areas 2007, 25, 94–97. [Google Scholar] [CrossRef]
  28. Cui, H.; Zhang, B.; Feng, X.; Gao, Y. Soil water-holding properties of different soil body configuration. Agric. Res. Arid Areas 2016, 34, 1–5. [Google Scholar] [CrossRef]
  29. Ge, J.; Huang, D.; Gao, X.; Tang, J.; Shen, H. Water Retention Capacity of Drained Soil Columns with Grained Layers. Southwest China J. Agric. Sci. 2019, 32, 2126–2132. [Google Scholar] [CrossRef]
  30. Chunyu, X.; Huang, F.; Xia, Z.; Zhang, D.; Chen, X.; Xie, Y. Assessing the Ecological Effects of Water Transport to a Lake in Arid Regions: A Case Study of Qingtu Lake in Shiyang River Basin, Northwest China. Int. J. Environ. Res. Public Health 2019, 16, 145. [Google Scholar] [CrossRef] [Green Version]
  31. Huang, F.; Chunyu, X.; Zhang, D.; Chen, X.; Ochoa, C.G. A framework to assess the impact of ecological water conveyance on groundwater-dependent terrestrial ecosystems in arid inland river basins. Sci. Total Environ. 2020, 709, 136155. [Google Scholar] [CrossRef]
  32. Huang, F.; Ochoa, C.G.; Chen, X.; Zhang, D. Modeling oasis dynamics driven by ecological water diversion and implications for oasis restoration in arid endorheic basins. J. Hydrol. 2021, 593, 125774. [Google Scholar] [CrossRef]
  33. Qian, J.; Chang, J.; Ma, H.; Ran, X.; Wang, Y. Comprehensive Water Resources Planning of Minqin County; Lanzhou University: Lanzhou, China, 2013; Unpublished work. [Google Scholar]
  34. Shao, M.; Wang, Q.; Huang, M. Soil Physics; Higher Education Press: Beijing, China, 2006. [Google Scholar]
  35. Chen, C.; Zhou, H.; Shang, J.; Hu, K.; Ren, T. Estimation of soil water content at permanent wilting point using hygroscopic water content. Eur. J. Soil Sci. 2020, 71, 392–398. [Google Scholar] [CrossRef]
  36. Taşan, S.; Demir, Y. Comparative Analysis of MLR, ANN, and ANFIS Models for Prediction of Field Capacity and Permanent Wilting Point for Bafra Plain Soils. Commun. Soil Sci. Plant Anal. 2020, 51, 604–621. [Google Scholar] [CrossRef]
  37. Ostovari, Y.; Asgari, K.; Cornelis, W. Performance evaluation of pedotransfer functions to predict field capacity and perma- nent wilting point using UNSODA and HYPRES Dataset. Arid Land Res. Manage. 2015, 29, 383–398. [Google Scholar] [CrossRef]
  38. Assi, A.T.; Blake, J.; Mohtar, R.H.; Braudeau, E. Soil aggregates structure-based approach for quantifying the field capacity, permanent wilting point and available water capacity. Irrig. Sci. 2019, 19, 3805–3816. [Google Scholar] [CrossRef]
  39. Wiecheteck, L.H.; Giarola, N.F.B.; de Lima, R.P.; Tormena, C.A.; Torres, L.C.; de Paula, A.L. Comparing the classical permanent wilting point concept of soil (−15,000 hPa) to biological wilting of wheat and barley plants under contrasting soil textures. Agric. Water Manag. 2020, 230, 105965. [Google Scholar] [CrossRef]
  40. Zhang, Y.; Chen, X.; Gao, M.; Liu, X. Simulation of transpiration for typical xeromorphic plants in inland arid region of Northwestern China. Acta Ecol. Sin. 2021, 41, 7751–7762. [Google Scholar] [CrossRef]
  41. Wang, X.; Shao, J.; Cui, Y.; Zhang, Q. Application of a Surrogate Model for a Groundwater Numerical Simulation Model for Determination of the Annual Control Index of the Groundwater Table in China. Sustainability 2020, 12, 5752. [Google Scholar] [CrossRef]
  42. Wang, J.; Zhang, G.; Yan, M.; Tian, Y.; Wang, Q. Index weight analysis of groundwater function evaluation and zoning system in arid areas. Trans. CSAE 2020, 36, 133–143. [Google Scholar] [CrossRef]
  43. Wang, J.; Zhang, G.; Wang, Q.; Cui, H.; Liu, P. Construction and application of evaluation index system of groundwater ecological function in northwest arid area. Acta Geol. Sin. 2021, 95, 1573–1581. [Google Scholar] [CrossRef]
  44. Wang, J.; Zhang, G.; Cui, H.; Wang, Q.; Dong, H.; Hao, J. System index attribute and application of groundwater function zoning in northwest inland area of China. J. Hydraul. Eng. 2020, 51, 796–804. [Google Scholar] [CrossRef]
Figure 1. Mechanism of maintaining natural vegetation by supportive capillary water.
Figure 1. Mechanism of maintaining natural vegetation by supportive capillary water.
Water 14 00059 g001
Figure 2. Remote-sensing image of the Shiyang River Basin and the study area of Minqin Basin.
Figure 2. Remote-sensing image of the Shiyang River Basin and the study area of Minqin Basin.
Water 14 00059 g002
Figure 3. The location map and photos of the survey pit. (a) location map of the survey pit; (b) pit section sampling; (c) cover the shading net after sampling; (d) pit profile under the shading net.
Figure 3. The location map and photos of the survey pit. (a) location map of the survey pit; (b) pit section sampling; (c) cover the shading net after sampling; (d) pit profile under the shading net.
Water 14 00059 g003
Figure 4. The structure and photograph of the test equipment.
Figure 4. The structure and photograph of the test equipment.
Water 14 00059 g004
Figure 5. Groundwater depths of the boreholes.
Figure 5. Groundwater depths of the boreholes.
Water 14 00059 g005
Figure 6. Calibration curve of different lithological soil filled in the column.
Figure 6. Calibration curve of different lithological soil filled in the column.
Water 14 00059 g006aWater 14 00059 g006b
Figure 7. Capillary rise height of groundwater measured by in situ pit test.
Figure 7. Capillary rise height of groundwater measured by in situ pit test.
Water 14 00059 g007
Figure 8. Cumulative height of wetting front and the cumulative water supply increasing with time.
Figure 8. Cumulative height of wetting front and the cumulative water supply increasing with time.
Water 14 00059 g008
Figure 9. Comparison of initial water content and final water content.
Figure 9. Comparison of initial water content and final water content.
Water 14 00059 g009
Figure 10. The accumulated height curve of wetting front increases with time in the early stage of the experiment.
Figure 10. The accumulated height curve of wetting front increases with time in the early stage of the experiment.
Water 14 00059 g010
Figure 11. The cumulative drainage of soil column changes with time.
Figure 11. The cumulative drainage of soil column changes with time.
Water 14 00059 g011
Figure 12. Variation curve of water content at different times after water release from soil column.
Figure 12. Variation curve of water content at different times after water release from soil column.
Water 14 00059 g012
Figure 13. Soil water characteristic curve fitted by Hydrus-1D.
Figure 13. Soil water characteristic curve fitted by Hydrus-1D.
Water 14 00059 g013
Table 1. Particle composition and dry bulk density of soil samples.
Table 1. Particle composition and dry bulk density of soil samples.
LithologyPercentage of Different Soil Particles (%)Dry Bulk Density (g/cm3)
>1 mm1–0.5 mm0.5–0.25 mm0.25–0.1 mm0.1–0.075 mm<0.075 mm
silty soil000.0447.7820.8231.361.35
clay soil0.234.653.1429.1728.9833.841.44
silty
fine sand
00.2417.8853.9015.3212.661.51
medium fine sand0.0116.8526.9333.4014.857.961.60
Table 2. The total amount of water supplied to the soil column during the capillary-water rising experiment.
Table 2. The total amount of water supplied to the soil column during the capillary-water rising experiment.
Testing Depth
(cm)
Lithologyθfinal (%)θinitial (%)θfinal – θinitial (%)ΔW (mm)
5Silty soil6.416.69
15Silty soil7.476.690.780.78
25Silty soil7.216.690.510.51
35Silty soil17.526.6910.8310.83
45Clay soil29.761.5428.2228.22
55Clay soil32.171.5430.6330.63
65Clay soil32.951.5431.4131.41
75Clay soil34.161.5432.6232.62
85Silty fine sand26.090.2625.8325.83
105Medium fine sand18.710.5818.1318.13
115Medium fine sand21.960.5821.3821.38
125Medium fine sand28.590.5828.0128.01
135Medium fine sand32.350.5831.7731.77
145Medium fine sand31.740.5831.1731.17
Cumulative water supply (mm)291.28
Testing Depth: distance from position of test point to top of soil column; θinitial: volumetric moisture content before the experiment; θfinal: volumetric moisture content after the experiment; ΔW: water storage at various depths, ΔW = (testing depthn − testing depthn-1) × (θfinal − θinitial); “—“: eliminated value which is negative because of the measurement error.
Table 3. Water holding capacity of soil column before and after the water-releasing experiment.
Table 3. Water holding capacity of soil column before and after the water-releasing experiment.
Testing Depth
(cm)
Lithologyθinitial (%)ΔW1 (mm)θfinal (%)ΔW2 (mm)
10Silty soil6.946.946.606.60
15Silty soil7.473.736.743.37
20Silty soil7.343.677.423.71
25Silty soil7.213.6010.935.47
30Silty soil12.376.1813.856.92
35Silty soil17.528.7614.597.29
40Silty soil23.6411.8223.6411.82
45Clay soil29.7614.8830.0715.04
50Clay soil30.9615.4830.9815.49
55Clay soil32.1716.0832.1916.09
60Clay soil32.5616.2833.0516.52
65Clay soil32.9516.4733.4316.72
70Clay soil33.5516.7833.6416.82
75Clay soil34.1617.0834.4217.21
80Clay soil30.1215.0627.4613.73
85Silty fine sand26.0913.0520.3010.15
90Silty fine sand24.2512.1216.108.05
95Silty fine sand22.4011.2012.266.13
100Silty fine sand20.5510.2810.385.19
105Medium fine sand18.719.3514.757.38
110Medium fine sand20.3310.1717.518.76
115Medium fine sand21.9610.9819.449.72
120Medium fine sand25.2712.6421.9610.98
125Medium fine sand28.5914.2925.0512.53
130Medium fine sand30.4715.2328.3014.15
135Medium fine sand32.3516.1729.1314.57
140Medium fine sand32.0416.0227.9613.98
145Medium fine sand31.7415.8721.4410.72
Water storage 340.21 305.09
Testing depth: distance from position of test point to top of soil column; θinitial: volumetric moisture content before the experiment; θfinal: volumetric moisture content after the experiment; ΔW1: initial water storage at various depths, ΔW1 = (testing depthn − testing depthn-1) × θinitial-n; ΔW2: final water storage at various depths, ΔW2 = (testing depthn − testing depthn-1) × θfinal-n.
Table 4. Fitting parameters and permanent wilting point of the different soils.
Table 4. Fitting parameters and permanent wilting point of the different soils.
LithologyFitting ParametersPermanent Wilting Point (%)
QrQsAlpha (1/cm)nKs (cm/days)l
Clay Soil0.08130.4310.01411.38958.650.512.47
Silty Soil0.08040.45360.01711.388420.210.512.37
Silty fine sand0.04960.39040.03131.503249.330.56.51
Medium Fine sand0.04250.36030.04011.613160.660.54.88
Qr: residual water content; Qs: saturated water content; Ks: saturated hydraulic conductivity; Alpha, n and l: empirical coefficients affecting the shape of the hydraulic functions in HYDRUS.
Table 5. Available water content of soil column after the water-releasing experiment.
Table 5. Available water content of soil column after the water-releasing experiment.
Testing Depth
(cm)
Lithologyθfinal (%)ΔW2 (mm)θpwp (%)θavailable (%)ΔW3 (mm)
10Silty soil6.606.6012.37
15Silty soil6.743.3712.37
20Silty soil7.423.7112.37
25Silty soil10.935.4712.370.93
30Silty soil13.856.9212.373.850.74
35Silty soil14.597.2912.374.591.11
40Silty soil23.6411.8212.3713.645.64
45Clay soil30.0715.0412.4720.078.80
50Clay soil30.9815.4912.4720.989.26
55Clay soil32.1916.0912.4722.199.86
60Clay soil33.0516.5212.4723.0510.29
65Clay soil33.4316.7212.4723.4310.48
70Clay soil33.6416.8212.4723.6410.59
75Clay soil34.4217.2112.4724.4210.98
80Clay soil27.4613.7312.4717.467.50
85Silty fine sand20.3010.156.5110.306.90
90Silty fine sand16.108.056.516.104.80
95Silty fine sand12.266.136.512.262.88
100Silty fine sand10.385.196.510.381.94
105Medium fine sand14.757.384.884.754.94
110Medium fine sand17.518.764.887.516.32
115Medium fine sand19.449.724.889.447.28
120Medium fine sand21.9610.984.8811.968.54
125Medium fine sand25.0512.534.8815.0510.09
130Medium fine sand28.3014.154.8818.3011.71
135Medium fine sand29.1314.574.8819.1312.13
140Medium fine sand27.9613.984.8817.9611.54
145Medium fine sand21.4410.724.8811.448.28
Water storage 305.09 182.54
Testing Depth: distance from position of test point to top of soil column; θfinal: volumetric moisture content after the experiment; θpwp: permanent wilting point, when the soil moisture content is lower than θpwp, the water will not be absorbed by the vegetation; θavailable: available volumetric moisture content, θfinal − θpwp; ΔW2: final water storage at various depths, ΔW2 = (testing depthn − testing depthn-1) × θfinal-n; ΔW3: available water storage at various depths, ΔW3 = (testing depthn − testing depthn–1)*θavailable-n; ”—”: eliminated value which is negative.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Cui, H.; Zhang, G.; Wang, J.; Wang, Q.; Lang, X. Influence of Multi-Layered Structure of Vadose Zone on Ecological Effect of Groundwater in Arid Area: A Case Study of Shiyang River Basin, Northwest China. Water 2022, 14, 59. https://doi.org/10.3390/w14010059

AMA Style

Cui H, Zhang G, Wang J, Wang Q, Lang X. Influence of Multi-Layered Structure of Vadose Zone on Ecological Effect of Groundwater in Arid Area: A Case Study of Shiyang River Basin, Northwest China. Water. 2022; 14(1):59. https://doi.org/10.3390/w14010059

Chicago/Turabian Style

Cui, Haohao, Guanghui Zhang, Jinzhe Wang, Qian Wang, and Xujuan Lang. 2022. "Influence of Multi-Layered Structure of Vadose Zone on Ecological Effect of Groundwater in Arid Area: A Case Study of Shiyang River Basin, Northwest China" Water 14, no. 1: 59. https://doi.org/10.3390/w14010059

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

Article Metrics

Back to TopTop