1. Introduction
The main source of groundwater recharge in arid and semi-arid regions is usually surface water from a river [
1]. As the supply of river surface water gradually declines, the hydraulic connection between the river and groundwater weakens, and eventually a disconnect occurs. Studies of the interaction between rivers and groundwater over the past decade have mainly focused on their hydraulic connection under a saturated state. The earliest study used two-dimensional integral equations to simulate groundwater-surface water interactions, demonstrating that the presence of an unsaturated zone below the riverbed results in the rapid flow of water into the aquifer [
2]. This flow preferentially travels laterally through the high-infiltration zone of the riverbank. Understanding the relationship between river water and groundwater and its evolution can facilitate the sustainable utilization, management, and protection of regional groundwater resources, particularly in arid and semi-arid regions in which the main source of groundwater recharge is river seepage [
3,
4]. Research on the relationship between river water and groundwater has for a long time mainly focused on their hydraulic connection under saturated conditions [
5]. Human demand for water is increasing with rapid human development. Consequently, there is increasing focus on understanding the relationship between river water and groundwater. There has been a gradual realization that groundwater levels will continue to decline as groundwater extraction increases. This decline in groundwater levels will lead to the loss of the hydraulic connection between river water and the groundwater under a saturated condition [
6,
7]. A disconnection between the river and groundwater results in the formation of an inverted saturated zone and an unsaturated zone under the river and between the inverted groundwater table and the groundwater table, respectively [
6,
8,
9].
The interaction between river water and groundwater is complex and is affected by many factors, including river morphology, riverbed composition, river hydrodynamic conditions, and human activities [
2,
10,
11,
12,
13]. The interactions between river water and groundwater have been extensively studied using analytical, numerical, and experimental methods [
1,
5,
8,
9,
14,
15,
16,
17,
18]. Using sandbox experiments, Riviere et al. [
19] studied the dynamic evolution of the relationship between river water and groundwater under conditions ranging from saturated to a complete disconnection. The results showed that the groundwater infiltration rate is regulated by the hydraulic gradient during the river and groundwater truncation process and tends to stabilize at a maximum value as the truncation stage completes. The results by Rivière [
19] are consistent with the relationship between the groundwater infiltration rate and river water level proposed by Fox and Durnford [
6] and Brunner et al. [
10]. A disconnect between river water and groundwater can only occur in partially connected rivers since fully connected rivers remain in a state of hydraulic connection under saturated conditions [
4]. In addition, the low permeability of a riverbed has a significant impact on the river water infiltration. Fox and Durnford [
6] and Brunner et al. [
10] studied the law of seepage between river water and aquifers and proposed that the development of a low-permeability layer in the riverbed is a necessary precondition for a disconnection between river water and groundwater. However, a study based on laboratory sandbox experiments by Wang et al. [
9] concluded that regardless of the presence of an aquitard at the bottom of the riverbed, the hydraulic connection between the river and groundwater will be broken under certain conditions. Wang et al. [
17] concluded that the thickness of the inverted saturation zone that develops under the riverbed is always equal to the depth of the river. However, this was contradictory to the results of Xie et al. [
18]. In addition, Xie et al. [
18] found that a large groundwater depth would develop in the absence of a barrier layer and only under a very shallow river water depth and narrow river channel.
Recent research has provided an increased understanding of the conditions under which the hydraulic connection between river water and groundwater under a saturated condition breaks as well as the corresponding laws of river infiltration. Xi’an et al. [
20] found that lateral flow resulting from capillary action or heterogeneity plays an important role in the relationship between river water and aquifers. They also found that the lowest point of the inverted water table can develop at the bottom of the riverbed. Xian et al. [
21] developed a bio-clogging model integrating saturated flow, microbial growth, and biogeochemical reactions to study the feedback relationship between nutrient transport and microbial growth during river and groundwater evolution and its control on soil infiltration.
The relationship between river water and groundwater is complex and affected by many factors. Therefore, the conditions under which the connection between river water and groundwater breaks remain unclear. Identifying the hydraulic connection between river water and groundwater can facilitate the sustainable utilization of regional groundwater resources and the protection of water resources. The predecessors put forward the theory that the recharge capacity of rivers to groundwater is limited, and rivers “recharge groundwater by vertical leaching during the flat and dry seasons and only recharge groundwater in a short period of time by means of uniform infiltration surface infiltration. When there is no obvious aquifer under the riverbed, under certain conditions, the groundwater level and the river will be “disconnected”. There is no doubt that the inverted saturated zone relates to the river water. Both the river water and the water in the inverted saturated zone show continuity and conduction of water pressure. Therefore, when exposed, the water in the inverted saturated zone can be discharged naturally and continuously under the action of gravity. In addition, the water in the inverted saturated zone is dynamic. The continuous downward leakage of water in the river results in the formation of an inverted saturated zone under the riverbed. In other words, the inverted saturated zone disappears when the river is not cut off or downward leakage no longer occurs.
The present study designed a laboratory sandbox experiment using three types of sand with different particle sizes. The aim of the present study was to identify the hydraulic connection between river water and groundwater under different aquifer mediums.
3. Results
The results of the sandbox experiment conducted in the present study were similar to those of previous studies. The sandbox experiments investigated three water bearing media, sand A, sand B, and sand C, and in each experiment, the river water depth was regulated at 18.5 cm and 10 cm.
Figure 3,
Figure 4,
Figure 5,
Figure 6,
Figure 7,
Figure 8,
Figure 9,
Figure 10 and
Figure 11 show the results of the change in the relationship between river and groundwater. In general, the hydraulic connection between river water and groundwater was gradually lost with the continuous decrease in discharged water level. This resulted in the formation of the river, inverted saturation zone, unsaturated zone, and saturation zone from top to bottom.
3.1. Sand C
Figure 3 shows the results of the sandbox experiment under sand C and a river water depth of 18.5 cm. Under a high discharged water level, the river and groundwater showed a uniform water table curve and a hydraulic connection under saturated conditions was maintained between them, with river water directly recharging groundwater. A decline in the level of discharge water to 140 cm resulted in the water table curve shrinking to the lower part of the river, with a tendency towards breaking the hydraulic connection. This result indicates that the relationship between the river and groundwater entered a disconnected transition stage. The hydraulic connection between groundwater and river water was lost at a discharge level of 130 cm, resulting in the formation of an inverted saturated zone with a thickness of 15–20 cm. However, the thickness of the inverted saturated zone stabilized with continuing reduction in the discharged water level.
The river level controller was adjusted to reduce the river water depth to 10 cm, after which the experiment was repeated. The results are shown in
Figure 4. In general, the hydraulic connection between river water and groundwater gradually weakened with a declining level of discharged water. There was a clear bend in the water table curve when the level of water discharged reached 160 cm, indicating a trend of disconnection of the hydraulic connection. There was a complete disconnect between river water and groundwater at a level of discharged water of 150 cm. An inverted water-saturated zone with a thickness of 5–10 cm was formed directly below the riverbed.
As shown in
Figure 5, discharge flow during stabilization was measured on both sides of the sandbox. The discharge flow under a river depth of 18.5 cm clearly exceeded that under a depth of 10 cm. This indicates that discharged flow showed a positive relationship with river depth. The discharged flow on both sides of the sandbox increased continuously with a decreasing level of discharged water under a hydraulic connection between the river and groundwater under saturated conditions. However, the discharged flow showed a slight decline followed by stabilization when this hydraulic connection was broken. This result indicates that recharge of groundwater by river water reaches a maximum after the hydraulic connection between the river and groundwater is broken.
3.2. Sand B
Figure 6 shows the results of the sandbox experiments under the aquifer medium sand B. When controlling the river water level at 18.5 cm, the hydraulic connection between the river and groundwater was maintained under a high level of discharged water, and the unsaturated zone only appeared on the left- and right-hand sides of the river at the upper part of the sandbox. The range of the unsaturated zone gradually expanded with a decline in the level of discharged water, and the water curve gradually extended to the lower part of the riverbed. The water table curve disconnected and the hydraulic connection between groundwater and the river was broken when the discharged water level was reduced to 150 cm. The thickness of the inverted saturated zone under the riverbed was 5 cm–10 cm. The size of the inverted saturated zone decreased slightly with the decline in the discharged water level to 140 cm. However, the maximum thickness of the inverted, saturated zone remained stabilized at 5–10 cm.
Reducing the depth of the river water to 10 cm resulted in the gradual loss of the connection between the river water and groundwater with a decreasing level of water discharge, and the connection finally breaking (
Figure 7) at a discharge level of 160 cm. At the period of hydraulic disconnection, the thickness of the inverted saturated zone was 5–10 cm. However, the size of the inverted saturation zone gradually decreased and disappeared within the continuing decrease in the level of discharged water. The size of the suspension-saturated belt could not be determined by the sandbox device (thickness of the suspension saturated belt < 5 cm) when the discharged water level reached 140 cm.
Figure 8 shows the relationship between discharge flow and discharge water level on both sides of the sand tank. The result indicated a gradual increase in discharged flow on both sides with decreasing discharged water level. The discharged flow eventually reached an extreme value, at which point the hydraulic connection between the river water and groundwater was broken. After this, the level of discharge continued to decline, whereas the discharged flow decreased and stabilized. In addition, under the same discharged water level, there was a positive correlation between the depth of river water and the discharge.
3.3. Sand A
Figure 9 shows the results of the experiment under water-bearing medium sand A and a river water depth of 18.5 cm. The results were consistent with those under sand B and sand C. Under a high level of discharged water, the hydraulic connection between the river and groundwater was maintained and a unified water table curve was noted. The range of the unsaturated zone gradually increased with a decreasing level of discharged water, and gradually expanded directly below the riverbed. The water table curve in the lower part of the riverbed showed clear shrinkage and narrowing when the drainage water level reached 80 cm and the hydraulic connection showed a weakening tendency. The hydraulic connection between the river and groundwater was in a state of transition. The river and groundwater no longer showed a uniform water table curve as the level of discharged water continued to decline to 60 cm, and the hydraulic connection was broken. An inverted saturated zone was then formed under the riverbed with a thickness of 15–20 cm.
Under a river water depth of 10 cm and a level of discharged water exceeding 140 cm, the hydraulic connection between the river and groundwater was maintained and the unsaturated zone only appeared on the left- and right-hand sides of the river (
Figure 10). The scope of the unsaturated zone gradually increased with a continuing decline in the level of discharged water. The water table curve under the riverbed began to contract when the level of discharged water dropped to 140 cm. This result indicated that there was a weakening of the hydraulic connection between the river water and groundwater. This hydraulic connection was broken when the level of discharged water reached 90 cm; the water table curve was disconnected, and an inverted saturated zone was formed under the riverbed. The thickness of the inverted saturated zone was then 20 cm–25 cm (
Figure 10). However, the size of the inverted saturated zone changed with a continued decline in the level of discharged water, decreasing to 10 cm–15 cm. The level of discharged water continued to decrease, and the groundwater table gradually decreased. However, the size of the inverted water saturation zone stabilized.
Explain that the three different sand discharge water levels are different and irregular. During the experiment, under the steady state of stable river water depth and zero flow, the authors lowered the water tank with a constant water head on the east and west sides of the discharge water level by 10 cm each time and observed the suspended saturated water zone of the river and groundwater surface, until the suspended saturated zone of the river becomes “disconnected” from the groundwater surface, i.e., loses hydraulic connection. Next, the discharge water level was increased to speed up the experiment, dropping 40 cm each time and observed the connection between the saturated water zone of the river and the groundwater surface. Take Sand B as an example; when the river water depth is 18.5 cm and the stabilization time exceeds 1 h, the discharge water level (movable water tanks on the east and west sides) stabilizes at 200 cm for the first time, and the pressure measurement and discharge data connected to the sand box are read. The water level of the second discharge is stable at 190 cm, and the pressure measurement and discharge data are read again, and so on, until the movable water tanks on the east and west sides are stable at 160 cm, and the river suspension saturated zone and the groundwater surface are about to be disconnected, that is, the connection surface between the river and the groundwater begins to shrink. When the discharge water level stabilizes at 150 cm, the river and groundwater are disconnected. In order to speed up the progress of the experiment again, the water tank controlling the discharge water level was lowered to 130 cm, and the pressure measurement discharge was observed again in the shape of the saturated zone. In order to express the experimental results concisely, all the data before the disconnection between the river and the groundwater are not drawn, and only the three sets of data in the middle are selected. Similarly, when the water depth of the river is 10 cm, the contact surface between the river and groundwater tends to decrease when the discharge water level of Sand B is 170 cm, and the disconnection occurs when the discharge water level is 140 cm. This set of experiments is the disconnection between the river and groundwater. For the group with denser former data, these groups of data are selected and drawn into the state diagram before disjointedness. The state of the last set of discharge water levels of 80 cm is selected to be compared with the state of the same discharge water level of a river water depth of 18.5 cm, and the last set of data of these experiments is when it drops to 80 cm.
Figure 11 shows the relationship between discharged flow and discharge water level as well as river water depth under sand A. Under the same discharge water level, the depth of the river water was proportional to the discharged water. Under the same river and groundwater water table curve, the level of discharged water was inversely proportional to discharged flow. However, the discharged flow stabilized after the hydraulic connection between the river and groundwater was broken. This result was consistent with the overall trend in changing discharged flow with changing discharged water levels in the sand B and sand C experiments.
As shown in
Figure 5, under river water depths of 18.5 cm and 10 cm and while the hydraulic connection between groundwater and river water was maintained, the discharge flow increased with a decreasing level of the discharged water level. The flow in the river showed a decreasing trend with a decreasing level of discharged water, which it stabilized with the same regularity on both sides of the river.
3.4. Discharge Flow
Figure 12 shows the relationship between discharged flow and discharge water level under the three types of aqueous media. Under the same river water depth and level of discharged water, the rank of the three sand treatments according to discharged flow was: sand C > sand B > sand A. This result corresponded with the permeability coefficients of the three sands (
KC > KB > KA) (
Table 1). The overall relationship between discharged flow and discharged water level noted in the present study was similar to that noted in previous studies, although some obvious differences were noted. Previous studies have shown that discharged flow reaches a maximum before the adjacent water level is reached, and the hydraulic connection between the river and groundwater is broken, and then slightly decreases and stabilizes. In addition, the results of the present study indicated that discharged flow reaches a maximum after the hydraulic connection between the river and groundwater is broken. However, previous studies [
10] have shown that the maximum discharged flow occurs before the hydraulic connection between the river and groundwater is broken, at which point the flow decreases slightly and then stabilizes. The authors have suggested that this observation can be attributed to the area of the unsaturated zone during the period of no hydraulic connection exceeding that of the unsaturated zone during the period of connection. However, this explanation seems weak since the area of the unsaturated zone increases with the decreasing level of discharged water during the period of the hydraulic connection.
Figure 12 shows the relationship between the discharged flow and discharged level of water at river water depths of 18.5 cm and 10 cm. The river was disconnected from the critical discharged water level of groundwater and the maximum discharged flow was reached. The discharged flow increased slightly after the hydraulic connection was broken, which then stabilized. The results under a river water depth of 18.5 cm were consistent with those under a river water depth of 10 cm.
In
Figure 12, the X-axis represents the drainage water level, and the value range is 40–200 cm. The left Y-axis represents the discharge flow of coarse sand at different river depths. When the water depth of the river is large, the corresponding discharge flow is also large, so the value range of the discharge flow in the above figure is larger than that in the below graph; the Y-axis on the right represents the discharge flow of medium sand and fine sand, the discharge flows are much smaller than that of coarse sand, so the value range is correspondingly smaller. The experimental results of the three water-containing media are that the discharge flow reaches the maximum value at the critical discharge water level where the river and groundwater are disconnected. After the disconnection, the discharge flow decreases slightly and remains basically stable. There are some differences in the experimental results (
Figure 12). In addition, in the experimental results of A sand and B sand, when the river water depth is 18.5 cm, the discharge flow reaches its maximum.
The degree of decline after the extreme value is greater than that when the river water depth is 10 cm, while the results of sand C are similar to the experimental results of sand A and sand B. When the water-containing medium is C sand, the discharge flow will not change after reaching the maximum value. Comparing the discharge flow of the three kinds of sand, it can be seen that under the same river water depth or discharge water level, the discharge flow of sand A is the smallest, followed by B sand, and then C sand, which has the largest discharge flow; that is, the greater the permeability coefficient, the greater the discharge flow.
4. Discussion
The present study assessed three kinds of sand with different particle sizes using a sandbox device. Although the results of the present study were similar to those of previous studies, some results were significantly different [
10]. The results of the present study that were consistent with those of previous studies included: (1) The unified infiltration curve for river and groundwater gradually lost its unity with the continuous decline in discharge water level. This resulted in the hydraulic connection being lost. The characteristic river→inverted saturated zone→unsaturated zone→saturated zone structure was then formed from top to bottom. (2) Under the same water depth, the level of discharge increased with decreasing level of discharge, finally reaching an extreme value. After the breakdown in the hydraulic connection between the river and the groundwater, the volume of discharge decreased slightly, following which it stabilized. (3) Under the same level of discharged water, the depth of the river water was proportional to the discharged flow and the effort needed to break the hydraulic connection.
The novel results obtained in the present study included:
- (1).
Discharge flow. A comparison of the discharge flow (the quantity of river water leakage) for the three sand types of different particle sizes under the same river water depth, and the discharged water level clearly showed that discharged flow was proportional to the diameter of the sand particles. The rank of the three sand types in terms of discharge flow was: coarse sand > medium sand > fine sand.
- (2).
Inverted saturated zone. The discharged water level resulting in the breaking of the hydraulic connection between river water and groundwater under a river depth of 18.5 cm was lower than that under a river depth of 10 cm. In other words, the thickness of the inverted saturated zone was inversely proportional to the depth of river water. There were different heights of discharge water level and different thicknesses of inverted saturated zone among different river water depths.
Previous studies have shown that the thickness of the inverted saturated zone is equal to the depth of the river water during periods of separation of the hydraulic connection between the river and groundwater [
10]. Under the same river water depth conditions, it seems unlikely that the inverted saturated zone formed in a gravel bed will have the same thickness that formed in a fine sand bed. In addition, the thickness of the inverted saturated zone under the river observed in the current study was not always equal to the depth of the river. For example, under a river water depth of 18.5 cm and the water-bearing medium sand C, the thickness of the inverted saturated zone was 10–15 cm, whereas that under water-bearing medium sand B was 5–10 cm and clearly smaller than the river water depth of 18.5 cm. The rank of the sand types according to the thicknesses of the inverted saturated zones under the same river water depth was sand A > sand C > sand B. Therefore, it is proposed that the thickness of the inverted saturated zone formed by the disconnection of the hydraulic link between the river and groundwater is not only affected by the depth of river water, but also related to the physical characteristics of the water-bearing medium.
The results showed that under the same conditions, the thickness of the inverted saturated zone changes with changing depth of the river water. The depth of the river water is proportional to the thickness of the suspended saturated zone. The inverted saturated zone under the river is unable to form when the depth of the river tends to 0. In addition, there is a large quantity of groundwater recharged by the river when the river depth is large. It results in the formation of the hydraulic connection between the river and the groundwater and will prevent the formation of the inverted saturated zone. Therefore, the thickness of the inverted saturated zone is related to river water depth, and within the critical range of river water depth, the thickness of the inverted saturated zone is directly proportion to river water depth.
The results show that under the same conditions, the thickness of the inverted saturation zone increases with the depth of the river. The depth of the river water is proportional to the thickness of the suspended saturation zone. When the river depth tends to 0, the river hardly infiltrates to recharge the groundwater, and the inverted saturation zone at the bottom of the river cannot be formed. In addition, when the river water depth is relatively large, a large amount of river water infiltrates to recharge groundwater, resulting in the formation of a hydraulic connection between the river and groundwater, and the formation of a suspended saturated zone. Therefore, the thickness of the inverted saturation zone is related to the river water depth, and within the critical range of the river water depth, the thickness of the inverted saturation zone is proportional to the river water depth (
Figure 13).
When the particles are very small, the permeable medium will become a weak permeable layer, thereby hindering the formation of a suspended saturated zone. When the particles are very small, the permeable medium will become a water-resisting layer, thereby hindering the formation of an inverted saturated zone. In contrast, when the particles are large, the water flow will directly pass through the water-bearing medium, thereby preventing the formation of an inverted saturated zone. Under certain conditions, an inverted saturated zone will be formed under the riverbed with increasing sand particle size. The inverted saturated zone will gradually increase with increasing sand particle size until the hydraulic connection between groundwater and the river is formed, at which point the inverted saturated zone will disappear. Therefore, within a certain critical condition, the thickness of the inverted saturated zone formed after the disconnection of the hydraulic connection between the river and groundwater is inversely proportional to porosity, while it is directly proportional to particle size (hydraulic conductivity), as shown in
Figure 14.
The present study determined the range of the inverted saturated zone by observing whether water flowed freely from the pressure tap after removing the pressure tube. Therefore, the inverted water table was located between two pressure taps that flow freely and contained no water. The range of thickness of the inverted saturated zone after the hydraulic connection between the river and the groundwater was broken, and the different water-bearing media conditions are shown in
Table 2. The present study assumed that the inverted water table obtained from the experiment was located in the middle of the two freely flowing pressure taps containing no water. The results suggested that the thickness of the inverted saturated zone was significantly correlated with porosity, as shown in
Figure 15. The amount of data on the thickness of the suspended saturated zone of the river obtained in this experiment is too small (only three). If it is expressed by a polynomial, it cannot more intuitively reflect the relationship between the suspended saturated zone and the river water depth and the porosity of the water-containing medium. In addition to the high fitting degree, the power function is selected for fitting, and it is also considered that the definition domain of the power function is generally in the first quadrant, that is, the porosity is a number greater than 0, which is more in line with the actual conditions.
The results of the present study indicated that the thickness of the inverted saturated zone formed below the riverbed after the breaking of the hydraulic link between the river and groundwater was significantly correlated with both the river depth and porosity. An equation is given to estimate the thickness of the inverted saturated zone:
In Equation (1), M is the thickness of the inverted saturated zone, α is a constant (and different for each aquifer medium), H is the depth of river water, and n is the porosity.
Table 3 shows the thicknesses of the inverted saturated zone calculated using this formula. All of the calculated thicknesses were within the range of the experimental observations. It should be noted that the provided formula is based only on the results of the current study and may not reflect field conditions.
Future studies should obtain more data to further verify and improve the estimation formula.
The evolution of the hydraulic relationship between rivers and groundwater is very complex. There has not yet been a comprehensive study of the interaction between rivers and groundwater, and there remains no effective method to study the relationship in practice, and this should be a focus of future research. The experiments conducted in the present study are based on the theory of the disconnect of the hydraulic connection between rivers and groundwater by Liu et al. [
22]. This theory states that the hydraulic connection between rivers and groundwater weakens and breaks during the dry season, and a suspended saturated zone is then formed. The results of a study by Feng et al., based on this theory, concluded that subsequent to the breaking of the hydraulic connection, the thickness of the suspended saturated zone is equal to the thickness of the river water layer [
23]. The present study tested this theory using three kinds of sand with different particle sizes [
24]. The results of the present study examined the distribution of negative pressure during the experiment and other conditions influencing the disconnect of the hydraulic connection between rivers and groundwater and the suspended saturated zone. The estimation formula proposed by the present study represents an advancement of understanding of these processes, although it requires further improvement.