Next Article in Journal
Prediction and Management of the Groundwater Environmental Pollution Impact in Anning Refinery in Southern China
Previous Article in Journal
Changes in the Characteristics of Zooplankton Communities in Response to Shifts in the Aquatic Environment in the Shallow Waters of Northern Liaodong Bay, China
Previous Article in Special Issue
Effects of Patch Properties of Submerged Vegetation on Sediment Scouring and Deposition
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Exploring the Effects of Fissures on Hydraulic Parameters in Subsurface Flows from the Perspective of Energy Changes

1
Hunan Academy of Agricultural Science, Changsha 410125, China
2
Key Laboratory of Arable Land Conservation (Middle and Lower Reaches of Yangtze River) of the Ministry of Agriculture, Huazhong Agricultural University, Wuhan 430070, China
3
State Key Laboratory of Water Resources and Hydropower Engineering Science, Wuhan University, Wuhan 430072, China
*
Author to whom correspondence should be addressed.
Yu Tao and Siwen Peng are co-first authors.
Water 2024, 16(19), 2712; https://doi.org/10.3390/w16192712
Submission received: 27 July 2024 / Revised: 29 August 2024 / Accepted: 30 August 2024 / Published: 24 September 2024
(This article belongs to the Special Issue Monitoring and Control of Soil and Water Erosion)

Abstract

:
Reynolds number (Re), pore water pressure (P), and water flow shear force (τ) are primary indicators reflecting the characteristics of subsurface flow. Exploring the calculation of these parameters will facilitate the understanding of the hydrodynamic characteristics in different subsurface flows and quantify their differences. Hence, we conducted a study to monitor soil water content, matrix potential, and pore water pressure in two typical soil profiles (with and without fissures). The distribution of Re, P, and τ in both matrix flow (MF) and preferential flow (PF) were calculated with an improved calculation method, focusing on their energy changes. Results showed that these hydrologic parameters are quite different between MF and PF. Re values in MF remained below 0.1, indicating lower water flow velocities, while the Re values ranged from 0.8 to 2 in PF, indicating higher flow velocities. The P values in PF was tens to hundreds of times higher than that in MF, which is mainly due to the rapid accumulation and leakage of water within soil fissures. Additionally, the larger hydraulic radius and gradient in PF also resulted in higher τ values in PF (2~6 N m−2) than in MF (0~1.5 N m−2). In PF, the pressure potential was the significant factor for τ, while τ in MF was dominated by the matrix potential and varies with the magnitude of the matrix potential gradient. This study suggests that Re, P, and τ could be considered as the major indexes to reflect dynamic characteristics of subsurface flow.

1. Introduction

Subsurface flow refers to water movement within the soil, either laterally along soil horizons or vertically through percolation. The factors influencing subsurface flow include soil physical properties, soil horizon stratification, and water supply conditions [1,2]. Additionally, features like macropores, animal burrows, and plant roots can facilitate the formation of subsurface flow [2]. Compared to surface runoff, subsurface flow is generally slower and can persist for longer durations, spanning days, weeks, or even longer [3]. The dynamics of subsurface flow play a crucial role in understanding its behavior and characteristics. Understanding subsurface flow dynamics allows researchers to calculate important parameters, such as the Reynolds number and water flow shear force, for different flow patterns. By identifying the key factors that influence these forces, researchers can better understand the mechanisms driving subsurface flow. Meanwhile, these parameters provide valuable information about the flow regime and its impact on the surrounding environment, such as contributing to the prevention of soil erosion induced by subsurface flow.
Matrix flow (MF) and preferential flow (PF) are two typical types of subsurface flow [4,5]. When MF occurs, soil moisture moves slowly in the capillary pores of soil, when PF occurs, soil moisture infiltrates rapidly directly through the large pores in the soil, resulting in a faster water flow than MF [4,5]. Hence, MF follows Darcy’s law, in which the flow rate is positively correlated with the pressure difference and the permeability of the soil [6]. However, the soil macropores enable rapid water movement, where PF occurs, and its flow characteristics do not follow Darcy’s law [7]. Both MF and PF are influenced by gravitational potential, matrix potential, and pore water pressure potentials [8]. In MF, the difference in matrix potential among adjacent soils drives subsurface flow, and the water always flows from the high-potential region to the low-potential region [9]. Nevertheless, in PF, water movement through macropores is not affected by matrix potentials when the macropores are saturated.
Pore water pressure is another driving force for subsurface flow, pushing soil water flow from areas of high pore water pressure to low, and this process is often accompanied by rapid and dramatic changes in pore water pressure [10]. Furthermore, changes in soil matrix potential and pore water pressure not only drive subsurface flow but also impact the physical properties and turbulent structure of water flow [11]. These changes result in alterations in energy loss and hydraulic characteristics. Therefore, understanding the dynamics of subsurface flow and its different flow patterns is essential for effective management and protection of groundwater resources.
Soil water potential, pore water pressure, water velocity, the Reynolds number (Re), and water flow shear force (τ) are primary parameters used to characterize the hydraulic characteristics of water flow. Among these, soil water potential and pore water pressure can be directly measured while the other three need to be calculated. The calculation of the water flow rate is straightforward, but Re and τ are slightly more complicated. Re is a dimensionless parameter that represents the ratio of inertial force to viscous force in fluid flow [12]. In 1883, Reynolds suggested that the flow pattern of a fluid is related to its flow velocity, the diameter of the tube, and the viscosity and density of the fluid (ρ) [12]. Thus, the Re of the subsurface flow can be calculated as Re = vd/4υ, where v (m s−1) is the linear velocity of the water in the soil, υ (N s−1 m−2) is the viscosity coefficient of water at a certain temperature, and d (m) is the diameter of the soil pores. When Re is small, the flow is mainly influenced by viscous forces, resulting in a stable flow field [13]. Conversely, when Re is large, inertial forces dominate, leading to an unstable flow field with turbulent flow [14]. In the case of water flowing vertically downward in a soil pipe, the flow is affected by matrix potential tension in addition to pressure, friction, inertial, and viscous forces [15]. The influence of the inertial force is negligible in this scenario, and the flow is primarily determined by the tension generated by matrix potential and pressure, with gravity acting as a constant term [16]. The soil water flow state and the friction and viscous forces vary with different Re values [17]. While friction and viscous forces cannot be directly observed, their resultant force can be numerically represented by τ, according to the Law of Conservation of Energy and the Principle of Balance of Forces.
The τ of subsurface flow is the force generated by the flow of water moving in the direction of the hydraulic gradient [18], which is mainly frictional and viscous force. These forces acting in the direction of the hydraulic gradient will scour the soil, disrupt the soil structure, and result in soil particles transporting out of the slope body along with the water flow [19,20]. According to the Law of Conservation of Energy, in uniform pipeline flow, the work performed by the water flow overcomes frictional resistance (f), also known as head loss (hf), caused by τ [21]. The formula for shear stress of thin-layer water flow under ideal conditions is given as τ = ρwgRJ, where τ represents water flow shear force, ρw denotes the density of water, g is the gravitational acceleration constant, R represents the hydraulic radius, and J signifies the hydraulic gradient [22]. The hydraulic gradient is crucial for calculating τ and is the biggest difference between calculating the τ of the subsurface flow and the pipeline flow [13].
Currently, limited research has been conducted on the calculation and application of Re and τ in subsurface flow. We argue that it is possible to draw on the principles and methods of pipeline flow hydraulics and improve the corresponding formulas for use in calculating the hydraulic factors of subsurface flow [23]. The Reynolds equation is the most typical calculation for Re, and τ is usually calculated with the water head loss equation which is followed by energy conservation. We attempt to apply the basic rules of Re and τ of pipeline flow to calculate τ in subsurface flow. Thus, the objectives of this study are the following: (1) to calculate Re for two typical subsurface flow patterns (matrix flow and preferential flow) and examine the effect of pore conditions on the size and distribution of Re and (2) to investigate the variations in τ under matrix flow and preferential flow, and explore the key factors affecting τ under these two subsurface flow patterns. This knowledge can contribute to the development of effective management and mitigation strategies for subsurface flow-related problems, ultimately leading to improved water resource management and decreased environmental degradation risk.

2. Materials and Methods

2.1. The Calculation of Hydraulic Factors

2.1.1. Reynolds Number (Re)

For the two typical subsurface flows, the preferential flow (PF) and the matrix flow (MF), the water flow velocity and Re can be expressed as the following equations:
υ = as/t
Re = vd/4υ
where v (m s−1) represents the water flow linear velocity in the soil, while υ (N s−1 m−2) is the viscosity coefficient of water at a specific temperature (in this study, it is 0.8 × 10−6 m2 s−1 at an air temperature of 30 °C). d (m) is the soil-equivalent pore diameter at a given matrix suction. The variables s (m) and t (s) denote the length of the water flow infiltration path (the distance between different soil layers) and the time taken for rainwater to infiltrate into different soil layers, respectively. The parameter a represents the curvature of the soil water infiltration channel (the ratio of the actual length of the path through which water flows to the distance between soil layers). According to the Brooks–Corey equation [24], a was valued as 2 in the calculation for MF, while it is the distance between different soil layers for PF and valued as 1 in this study.
The most important aspect of the calculation process was the determination of the infiltration time and the infiltration rate. The average infiltration rate is computed by considering the length of the water flow infiltration path and the time taken by rainwater to infiltrate into different soil layers. When calculating Re for MF, the starting time point selected is when rainwater begins to maintain a stable and uniform infiltration between the two soil layers. Conversely, for PF, the starting time point is when PF begins to occur between the two soil layers during the infiltration process. The equivalent pore diameter is determined by the matrix suction, which is calculated using the soil water content when the rainwater reaches a specific soil layer. The relationship between soil water content and matrix suction is obtained from the soil water characteristic curve (SWCC).

2.1.2. Water Flow Shear Force (τ)

According to the Law of Conservation of Energy, the energy balance equation of the water flow movement process in a pipeline flow (Figure 1) can be described as follows [25]:
h m + P m ρ w g + a m v m 2 2 g = h n + P n   ρ w g + a n v n   2 2 g + h f
Since there is a uniform flow, acceleration is zero and thus vm = vn. Equation (3) can be simplified as follows:
  h f = h m h n + P m ρ w g P n ρ w g
The equation for all external forces in the direction of motion is as follow:
  P m A P n A + ρ w g A l c o s θ = τ χ l  
where
c o s θ = h m h n l  
R = A χ = d 4
J = h f l  
Per the above equations, τ is calculated as follows:
τ = ρ w g d 4 h f l
where h is the potential energy, α (≈1) is the kinetic energy correction factor, P (Pa) is the dynamic water pressure at both ends, A (m2) is the area of the water flow section, l (m) is the linear distance of the water flow path, and χ (m) is the wetted perimeter.
For subsurface flow where water moves vertically downwards in a soil pipe, the water flow shear force of subsurface flow can be expressed as in Figure 2. Due to the influence of the matrix potential and gravity, based on the Law of Conservation of Energy, the energy balance of the water flow movement process follows Equations (10) and (11).
h m + P m ρ w g + α m v m   2 2 g + φ m = h n + P n ρ w g + α n v n   2 2 g + φ n + h s
h s = h m h n + P m ρ w g P n ρ w g + φ m φ n + α m v m   2 α n v n   2 2 g
In PF (Figure 3), flow velocity changes and the acceleration is (αnvn2αmvm2)/2l. The energy required for acceleration is negligible and can be disregarded during calculations. The soil surrounding the water pathway is fully saturated, resulting in a zero-matrix potential. The equilibrium of all external forces in the direction of motion follows Equation (12).
P m A P n A + ρ w g A l ρ w g A α m v m   2 a n v n   2 2 l = τ χ l
At this time, τ is calculated as follows:
τ = ρ w g d 4 h s l
For MF, the hydraulic gradient is as follows:
J = ψ l
In this case, τ is calculated as follows:
τ = ρ w g d 4   ψ l
According to Darcy’s law, i is derived as follows:
i = k d ψ d y
The calculation of τ for MF is further simplified as follows:
τ = i k d 4 ρ w g
where ψ (m) is the soil water potential, k (m s−1) is the saturated hydraulic conductivity, and i (m s−1) is the water flow flux.

2.2. Field Validation

2.2.1. Study Location Description

The study area is situated in Tongcheng County (29.33° N, 113.77° E), Hubei Province, China. It exhibits a typical subtropical monsoon climate characterized by an annual rainfall of 1512.8 mm and an average annual temperature of 16.1 °C (China Meteorological Administration). The soil in this region predominantly originates from granite parent material, which has undergone significant weathering due to consistent rainfall and high temperatures. Consequently, the soil in the study area is loose and deep, ranging from 10 to 20 m in depth. According to the U.S. Soil Taxonomy System, it falls under the classification of Ultisols.
The study site is an active Benggang, characterized by vertical collapse due to water and gravity [26]. It consists of a catchment slope, gully head (Benggang wall), avalanche pile, flow channel, and alluvial fan (Figure 4). Our study focused on a representative Benggang in an active area where three collapses occurred during observation. The Benggang wall is steep, with a height difference of 5.2 m and a slope angle of over 75°. Soil layering varies with slope position: an upper shallow surface soil layer (layer A), a middle and lower red clay layer (layer B) exposed by water erosion, and even a sandy layer (layer BC) in some areas. The depth of 20–40 cm is consistently red clay, followed by a 40–80 cm transition layer, and a 60–80 cm sandy soil layer on the Benggang wall. The detritus layer (layer C) is deep and not visible in the soil profile.

2.2.2. Monitoring of Soil Water Content, Pore Water Pressure, and Matrix Suction

To investigate the changes in soil water dynamics during rainfall infiltration, two soil profiles were selected on the Benggang wall (Figure 4). These profiles were formed naturally as a result of the Benggang collapse. Profile A represents a homogeneous soil profile, while profile B exhibits visible fissures on the exposed side of the Benggang (Figure 5).
At specific depths of 20, 40, 60, and 80 cm on both profiles, we embedded soil water sensors (SMEC 300 SM/EC/Temp Sensor, Spectrum Technologies, Inc., Chicago, IL, USA), pore pressure sensors (Navigation Technology Company, Xi’an, China), and micro-tensiometers (Umwelt-Geräte-Technik GmbH, Berlin, Germany). These sensors were positioned consistently between the profiles. A 2 min resolution was used for the monitoring of soil water content, pore water pressure, and soil matrix potential. In the cracked profile (profile B), the fissure widths at the sensor positions were recorded as 2.11 mm, 2.34 mm, 1.98 mm, and 2.02 mm. Table 1 provides an overview of the soil physical properties for the four soil layers present in profiles A and B.

2.2.3. Measurements of Soil Macropores

According to the saturated volumetric soil water content and residual soil water content (drying method) of the different soil layers in the two profiles, and the soil water characteristic curve (SWCC), the soil macropores distribution was obtained. Among them, the soil water characteristic curve is drawn using RETC software 6.02. The soil water content under a series of matrix suction (under the suction of 15 cm-H2O, 9 cm-H2O, and 3 cm-H2O) was measured with a micro tensiometer, and then the obtained soil water characteristics data were input into the RETC computer program (American Salinity Laboratory, 4500 Glenwood Drive, Riverside, CA 92051, USA). In RETC, the VG model [27] is used to fit and draw the SWCC. According to Brooks’ research, the equivalent pore diameter is calculated as d = 3/h, where h (cm-H2O) is the height of the water applied to the matrix suction, d (mm) is the diameter of the soil pore equivalent to the suction [28]. Different matrix suctions represent different soil pore diameters. Finally, according to the change in soil water content under different suction conditions, the specific distribution of soil macropores with different diameters is calculated.

2.3. Data Statistics and Analysis

All statistical analyses were performed with the software Origin 2021 (Origin Lab Corporation, Northampton, MA, USA). An alpha of 0.05% was used to identify significant results.

3. Results

3.1. The Reynolds Number under Different Subsurface Flows

Table 2 presents the relevant calculation parameters and Re values during the initial stage of infiltration for three rainfall events, considering pore pressure. When subsurface flow is driven by matrix flow (MF), Re values are generally small. Among the nine scenarios (three soil layers × three rainfall events), profile A consistently showed Re values for MF below 0.1. The highest Re value (0.091) was observed in the 20–40 cm soil layer during the second artificial rainfall. Re values decrease from the surface soil to deeper layers, with the lowest values appearing at 60–80 cm. The velocity of water flow in MF is also very low, usually below 0.1 mm s−1, and there is a time lag of over an hour between the soil water content probes in different layers (Figure 6). This suggests that the small Re values during MF can be attributed to the extremely slow infiltration velocity of subsurface flow.
In profile B, subsurface flow occurs as preferential flow (PF) in fissures, and the velocity of water flow in PF is much higher than that in the soil matrix (Figure 7). When PF drives subsurface flow, the Re values due to pore pressure are tens-of-times greater than that in MF, which ranged from 0.8 to 2. However, due to the limitations of the soil layer thickness (20 cm) and monitoring frequency (every 5 min), the maximum Re value during PF cannot be measured. Under the same rainfall intensity, the infiltration velocity of rainwater in PF is significantly higher than that in MF, generally exceeding 1 mm s−1. This fast infiltration velocity likely leads to higher Re values. Additionally, in profile B, the deeper soil layer is predominantly sandy with numerous macropores (Table 1). Consequently, Re values tend to increase from the surface to deeper layers in most rainfall events (Table 2).

3.2. The Distribution of Soil Pore Pressure under Different Subsurface Flows

Pore water pressures were compared between profiles A and B in three rainfall events to examine the impact of soil water content on pore water pressure distribution in MF and PF (Figure 8). High soil water content promotes pore water pressure. During MF, when water content is below 0.3 m3 m−3, pore water pressures in different soil layers are close to 0 kPa. Nevertheless, pore water pressures become polarized when soil water content exceeds 0.3 m3 m−3. However, certain soil water contents correspond to multiple pore water pressure values, especially at the 20 cm and 40 cm soil layers, where maximum and minimum values of 2 kPa and −2 kPa appear, respectively, in the 40 cm soil layer. Due to the thickness of the soil layer and its depth, positive pore water pressure is more likely, particularly in the deeper layer (such as 80 cm) (Figure 8a).
In PF, as soil water content increases, pore water pressure in each soil layer gradually transitions to positive pressure and continues to rise (Figure 8b). Generally, pore water pressure throughout the entire profile is positive in PF, except for the 20 cm soil layer where soil water content is below 0.3 m3 m−3 and an isolated area in the 40 cm soil layer. In deeper soil layers, pore water pressures are positive during infiltration, except a minimum value of −1 kPa observed at low soil water content in the 20 cm soil layer. At the 80 cm soil layer, pore water pressure increases linearly with water content, reaching a maximum value of 1.5 kPa. In PF, the rapid pooling of soil water may be the reason for the linear increase in pore water pressure.

3.3. The Water Flow Shear Force under Different Subsurface Flows

Table 3 presents the distribution of shear stress (τ) in profiles A and B at the initial stage of rainwater infiltration across three rainfall events. In MF, τ values are generally below 2 N m−2, with the maximum value (1.729 N m−2) occurring in the 40~60 cm soil layer. Noteworthy, the τ generated via matrix potential is tens to hundreds of times higher than that via pore water pressure, with the magnitude of τ determined by the matrix potential gradient. In PF, τ is approximately ten times higher than in MF. The value of τ in PF is typically above 2 N m−2, with the maximum value (6.949 N m−2) occurring in the 60~80 cm soil layer. Moreover, the τ value generated via gravity is fixed and only related to the vertical height difference and soil pore diameter. However, the pore water pressure in the fissure of PF is tens to hundreds of times higher than that in MF, and it is the key factor determining τ in PF.
Due to the long equilibrium time and lag in the determination of soil matrix potentials, we calculated τ for MF according to Darcy’s law. Table 4 shows that the calculated τ for MF ranges from 0 to 2 N m−2, with a maximum value of 1.606 N m−2. The water flow flux at this time is 0.083 mm s−1. This result is consistent with the calculation based on the Law of Energy Conservation. However, the size of τ calculated with the two methods differs in the same soil layer. When using Darcy’s law, water flow flux is the determining factor for τ. In our three rainfall events, water flow flux gradually decreased from the surface downward, while saturated hydraulic conductivity increased. Two of the rainfall events indicated that the Re for MF decreased from the surface downward.
Figure 6 and Figure 7 provide insights into the continuous changes in τ during rainfall for MF and PF, respectively. These values were calculated using Darcy’s law and the Law of Energy Conservation. In MF, τ remains below 0.5 N m−2 throughout the rainfall process. The influence of τ becomes evident only when the soil water content begins to change, and the response of τ in a soil layer is positively correlated with the change in water content. However, the water content level does not determine the magnitude of τ between different soil layers. In the case of PF, τ undergoes significant changes during rainfall and when MF occurs. Notably, τ is particularly high at the beginning of infiltration. In Figure 8, changes in τ often occur when there are substantial changes in soil water content or pore water pressure. However, an increase in soil water content does not necessarily lead to an increase in τ. When there is a steady increase or decrease in soil water content, τ remains relatively stable. The same holds for changes in τ caused by pore water pressure, where a sharp change in pore water pressure leads to changes in τ depending on the difference in pore water pressure at both ends.

4. Discussion

4.1. The Distribution of Re in Different Profiles

Re is a crucial parameter in analyzing subsurface flow dynamics, with its value depending on the type of flow. Understanding the distribution of Re values helps to comprehend the behavior of subsurface flow, which is essential for various applications like groundwater management, material transport, and energy conversion. In this study, we focused on quantifying Re in two drastically different subsurface flows.
When MF dominated the subsurface flow, Re is small and the maximum value is below 0.1. During the infiltration process, MF follows Darcy’s law, and the soil pore diameter used to transport MF is small, both of which result in a slow flow velocity of rainwater infiltration [29] and a lower Re. Meanwhile, Re with dynamic pressure gradually decreases with an increase in profile depth, which is mainly caused by the gradual decrease in water velocity in a profile [30]. Additionally, the number of macropores gradually decreases in the deeper profile. Since macropores can promote the movement of water flow, Re with MF in the upper part of a profile, with more macropores, also will be larger.
Conversely, for PF, the water movement does not follow Darcy’s law, but bypasses the soil matrix and moves downward in a disorderly and inhomogeneous manner through the cracks [31]. Due to the larger width of fissures, water flow velocity is faster, and the Re for PF is much larger than MF. We found that when PF occurs, the Re for PF is nearly always greater than 1. In addition, during downward infiltration with PF, the soil pipe wall is saturated and not affected by the matrix potential [32]. In this situation, gravitational potential is constantly converted into kinetic energy and the velocity of water flow in the soil is constantly accelerating under the action of gravity. As the profile deepens, the Re for PF will continue to increase.

4.2. The Effects of Flow Regimes on Pore Water Pressure

The distribution of pore water pressure during subsurface flow is intricately linked to the continuous change in soil water content in different soil layers [8]. As soil water flows, it aggregates or leaks, resulting in varying water pressure states in the soil pores [33]. This interplay between soil water content and pore water pressure in different subsurface flows is complex and fascinating. When MF dominates subsurface flow, the pore water pressure rapidly polarizes as the water content increases. In different soil layers, the uniform distribution of soil pores and soil texture results in a different hydraulic conductivity, which plays a significant role in determining the variation in pore water pressure [34]. When the hydraulic conductivity of the lower soil layer is greater than that of the corresponding upper layer, a water scarcity zone is easy to form between adjacent soil layers, resulting in cavities and negative pore water pressure in the soil body. Conversely, when the lower layer’s hydraulic conductivity is lower, rainwater will easily accumulate between adjacent soil layers, resulting in a stagnant layer and positive pore water pressure [35]. Noteworthy, with the high soil water condition, the soil body is more prone to produce stagnant layers or cavities, which makes the change in positive and negative pore water pressure especially obvious in this stage [36]. Meanwhile, during rainfall events, the hydraulic conductivity of each soil layer will also change with the variation in soil water content. This suggests that the changes in soil water content within any of the two adjacent soil layers are accompanied by a change in hydraulic conductivity. Our results show that in general, the pore water pressure can vary even with constant soil water content, and different stages of rainfall can result in different pore water pressures for the same water content.
The distribution of pore water pressure in subsurface flow is influenced by the change in soil water content within different soil layers. As soil water content increases, pore water pressure transitions to positive pressure and continues to rise. In PF, the water flow bypasses the soil matrix and moves quickly through fissures, and the soil structure has little influence on water flow velocity [37]. Positive pressure is generated on the fissure wall as soil water fills the fissures [38]. However, PF can also cause blockages and dredging of fissures, leading to changes in positive and negative pore water pressures [39]. When fissures are blocked, positive pore water pressure increases, and when fissures are drained, negative pore water pressure occurs.

4.3. Factors Controlling Water Flow Shear Force

Soil water content, pore water pressure, and subsurface flow are interconnected factors that play a crucial role in various applications. τ is a function of the hydraulic radius and the hydraulic gradient, with its magnitude depending on the combination of the hydraulic radius and the hydraulic gradient [40]. While the hydraulic radius is constant, the difference in hydraulic gradient is key to controlling τ under different flow conditions. The hydraulic gradient is influenced by gravity potential, pressure potential, matrix potential, and kinetic energy of water flow [13]. Under different conditions, one or more of these energies may dominate the hydraulic gradient [41]. In addition, since soil water flow velocities are usually small, the variations in the velocity of subsurface flow have less influence on the changes in water flow energy. Therefore, τ of subsurface flow is usually calculated using a constant flow velocity [42].
When MF occurs, gravity potential and matrix potential differences become the dominant forces affecting τ [43]. As the gradient value of the gravitational potential energy is fixed, matrix potential plays the most critical factor. In our research, τ caused by matrix potential difference is tens to hundreds of times greater than τ generated by pressure potential difference when MF occurs (Table 3). The magnitude of the matrix potential difference largely controls τ in different soil layers, and it correlates positively with soil water content. Additionally, due to the continuous downward movement of water flow and its decreasing hydraulic gradient, the direction of τ is always downward [44,45]. However, since the measurement of soil matrix potential requires a long equilibrium time, it is difficult to accurately measure and also makes the calculation of τ during MF complex. It is a better choice to calculate τ using Darcy’s law; in this process, the water flow flux can be obtained by the change in soil water content per unit time.
In PF, the wall of the fissure quickly saturates, rendering the matrix potential irrelevant, and then the pressure potential becomes the key factor in controlling τ [46]. The accumulation or leakage of soil water is a key cause of changes in pressure potential and streamflow shear [47]. In our study, the water flow rate in PF is fast, and the rapid accumulation and seepage of soil water within the fissure result in a pressure potential energy that is tens to hundreds of times higher than that of MF (Table 3). Furthermore, the hydraulic radius in PF is much larger than in MF, resulting in τ in PF being tens to hundreds of times greater.

4.4. Understanding Subsurface Flow Dynamics and Energy Changes

Our results show that the changes in hydrodynamic factors driving subsurface flow can be calculated by analyzing the changes in potential and kinetic energy during different subsurface flow motions. In MF, potential energy is the key factor controlling the state of subsurface flow. In PF, the kinetic energy state of water flow and the change in pressure potential are the key factors controlling the state of subsurface flow. This provides us with a new perspective of thinking in regard to subsurface flow. Under field conditions, we can infer and predict the specific state of subsurface flow by monitoring energy changes in the soil profile. For example, some measurable factors such as soil temperature and soil pore pressure are used to determine the form of subsurface flow. These parameters provide a new way to calculate the hydraulics driving subsurface flow in different soils (e.g., rate of water flow, water flow shear force, Reynolds number, etc.). The Law of Conservation of Energy can be used to solve hydraulic equations and provide us with new parameters for the simulation of subsurface flow and the establishment of subsurface flow models. When other parameters cannot be determined due to the state of subsurface flow, energy changes can be used as a reference. Thus, in the absence of other suitable methods, justifying the movement of subsurface flow from the perspective of energy exchanges may result in new ways of thinking about subsurface flow.

5. Conclusions

In our study, we made an improvement in the calculation method for hydraulic factors in uniform pipe flow. By applying this improved method, we were able to analyze the distribution of Reynolds number (Re), pore water pressure, and τ in two profiles under two different flow regimes: preferential flow (PF) and matrix flow (MF). In the MF regime, the water flow velocity and Re values were generally below 0.1 mm s−1 and 0.1, respectively. However, under the PF regime, Re values were significantly higher, ranging from 0.8 to nearly 2. We identified that high water content plays a crucial role in the generation of pore water pressure across different flow regimes. Specifically, greater water content resulted in higher pore water pressure. These findings attempt to analyze the characteristics of different subsurface flow from the perspective of energy changes and will provide some guidelines for exploring soil erosion. Moreover, this approach also holds promise for gaining a deeper understanding of the mechanisms and behavior of groundwater flow.

Author Contributions

Conceptualization, Y.T. and J.C.; methodology, Y.T.; software, S.P.; validation, S.L. and Y.T.; formal analysis, B.L.; investigation, Y.T. and B.L.; resources, S.P. and S.L.; data curation, Y.T. and B.L.; writing—original draft preparation, Y.T. and B.L.; writing—review and editing, B.L.; visualization, Y.T. and B.L.; supervision, J.C.; project administration, J.C. and Y.T.; funding acquisition, J.C. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the National Natural Science Foundation of China, grant number 42177307 and 41571258, and the National Key R&D Program of China, grant number 2021YFD1500703.

Data Availability Statement

The original contributions presented in this study are included in the article. Further inquiries can be directed to the corresponding authors.

Conflicts of Interest

The authors declare no conflicts of interest.

References

  1. Renard, P.; Allard, D. Connectivity metrics for subsurface flow and transport. Adv. Water Resour. 2013, 51, 168–196. [Google Scholar] [CrossRef]
  2. Guo, L.; Lin, H. Addressing two bottlenecks to advance the understanding of preferential flow in soils. Adv. Agron. 2018, 147, 61–117. [Google Scholar]
  3. Lin, H. Temporal stability of soil moisture spatial pattern and subsurface preferential flow pathways in the shale hills catchment. Soil Sci. Soc. Am. 2006, 5, 317–340. [Google Scholar] [CrossRef]
  4. Preziosi, L.; Farina, A. On Darcy’s law for growing porous media. Int. J. Non-Linear Mech. 2002, 37, 485–491. [Google Scholar] [CrossRef]
  5. Clothier, B.; Green, S.; Deurer, M. Preferential flow and transport in soil: Progress and prognosis. Eur. J. Soil Sci. 2007, 59, 2–13. [Google Scholar] [CrossRef]
  6. Glass, R.; Nicholl, M.; Pringle, S.; Wood, T. Unsaturated flow through a fracture–matrix network: Dynamic preferential pathways in mesoscale laboratory experiments. Water Resour. Res. 2002, 38, 17-1–17-17. [Google Scholar] [CrossRef]
  7. Tyner, J.; Wright, W.; Yoder, R. Identifying long-term preferential and matrix flow recharge at the field scale. Trans. ASABE 2007, 50, 2001–2006. [Google Scholar] [CrossRef]
  8. Rahardjo, H.; Leong, E.; Rezaur, R. Effect of antecedent rainfall on pore-water pressure distribution characteristics in residual soil slopes under tropical rainfall. Hydrol. Process. 2008, 22, 506–523. [Google Scholar] [CrossRef]
  9. Campbell, G. Soil water potential measurement: An overview. Irrig. Sci. 1988, 9, 265–273. [Google Scholar] [CrossRef]
  10. Uchida, T.; Asano, Y.; Mizuyama, T.; McDonnell, J. Role of upslope soil pore pressure on lateral subsurface storm flow dynamics. Water Resour. Res. 2004, 40, W12401. [Google Scholar] [CrossRef]
  11. Zhang, W.; Tang, X.; Xian, Q.; Weisbrod, N.; Yang, J.; Wang, H. A field study of colloid transport in surface and subsurface flows. J. Hydrol. 2016, 542, 101–114. [Google Scholar] [CrossRef]
  12. Reynolds, O. An experimental investigation of the circumstances which determine whether the motion of water shall be direct or sinuous, and of the law of resistance in parallel channels. Proc. R. Soc. Lond. 1883, 35, 84–99. [Google Scholar]
  13. Atwood, D.; Gorelick, S. Hydraulic gradient control for groundwater contaminant removal. J. Hydrol. 1985, 76, 85–106. [Google Scholar] [CrossRef]
  14. Reynolds, O. On the dynamical theory of incompressible viscous fluids and the determination of the criterion. Proc. R. Soc. Lond. 1894, 56, 40–45. [Google Scholar]
  15. Towner, G. The application of classical physics transport theory to water movement in soil: Development and deficiencies. J. Soil Sci. 1989, 40, 251–260. [Google Scholar] [CrossRef]
  16. Lu, T.; Biggar, J.; Nielsen, D. Water movement in glass bead porous media: 1. Experiments of capillary rise and hysteresis. Water Resour. Res. 1994, 30, 3275–3281. [Google Scholar] [CrossRef]
  17. Vasil’ev, M. On the force acting on a cylinder in a steady stream of viscous fluid at low reynolds number. J. Appl. Math. Mech. 1981, 45, 632–635. [Google Scholar] [CrossRef]
  18. Govan, A.; Hewitt, G.; Owen, D.; Burnett, G. Wall shear stress measurements in vertical air-water annular two-phase flow. Int. J. Multiph. Flow 1989, 15, 307–325. [Google Scholar] [CrossRef]
  19. Iverson, R.; Reid, M. Gravity-driven groundwater flow and slope failure potential: 1. Elastic Effective-Stress Model. Water Resour. Res. 1992, 28, 925–938. [Google Scholar] [CrossRef]
  20. Liu, Y.; Yang, J.; Hu, J.; Tang, C.; Zheng, H. Characteristics of the surface–subsurface flow generation and sediment yield to the rainfall regime and land-cover by long-term in-situ observation in the red soil region. South. China J. Hydrol. 2016, 539, 457–467. [Google Scholar] [CrossRef]
  21. Goulter, I.; Lussier, B.; Morgan, D. Implications of head loss path choice in the optimization of water distribution networks. Water Resour. Res. 1986, 22, 819–822. [Google Scholar] [CrossRef]
  22. Lilleleht, L.; Hanratty, T. Relation of interfacial shear stress to the wave height for concurrent air-water flow. AIChE J. 1961, 7, 548–550. [Google Scholar] [CrossRef]
  23. Dong, Q.; Xu, D.; Zhang, S.; Bai, M.; Li, Y. A hybrid coupled model of surface and subsurface flow for surface irrigation. J. Hydrol. 2013, 500, 62–74. [Google Scholar] [CrossRef]
  24. Brooks, R.; Corey, A. Hydraulic Properties of Porous Media: Hydrology Papers 3; Colorado State University: Fort Collins, CO, USA, 1964. [Google Scholar]
  25. Reichardt, K.; Timm, L. The Movement of Water in the Systems. In Soil, Plant and Atmosphere; Springer: Berlin/Heidelberg, Germany, 2020. [Google Scholar]
  26. Shi, D. Soil erosion and its control in the granite region of Southern China. J. Soil Water Conserv. 1991, 3, 63–72. (In Chinese) [Google Scholar]
  27. Van-Genuchten, M. A closed-form equation for predicting the hydraulic conductivity of unsaturated soils. Soil Sci. Soc. Am. J. 1980, 44, 892–898. [Google Scholar] [CrossRef]
  28. Brooks, R.; Corey, A. Properties of porous media affecting fluid flow. J. Irrig. Drain. Div. 1966, 92, 61–88. [Google Scholar] [CrossRef]
  29. Nimmo, J.; Rubin, J.; Hammermeister, D. Unsaturated flow in a centrifugal field: Measurement of hydraulic conductivity and testing of Darcy’s law. Water Resour. Res. 1987, 23, 124–134. [Google Scholar] [CrossRef]
  30. Dakshanamurthy, V.; Fredlund, D. A mathematical model for predicting moisture flow in an unsaturated soil under hydraulic and temperature gradients. Water Resour. Res. 1981, 17, 714–722. [Google Scholar] [CrossRef]
  31. Greco, R. Preferential flow in macroporous swelling soil with internal catchment: Model development and applications. J. Hydrol. 2002, 269, 150–168. [Google Scholar] [CrossRef]
  32. Ludwig, R.; Gerke, H.; Wendroth, O. Describing water flow in macroporous field soils using the modified macro model. J. Hydrol. 1999, 215, 135–152. [Google Scholar] [CrossRef]
  33. Revol, P.; Vauclin, M.; Vachaud, G.; Clothier, B. Infiltration from a surface point source and drip irrigation: 1. The midpoint soil water pressure. Water Resour. Res. 1997, 33, 1861–1867. [Google Scholar] [CrossRef]
  34. Chu-Agor, M.; Fox, G.; Cancienne, R.; Wilson, G. Seepage caused tension failures and erosion undercutting of hillslopes. J. Hydrol. 2008, 359, 247–259. [Google Scholar] [CrossRef]
  35. Shoushtari, S.; Cartwright, N.; Perrochet, P.; Nielsen, P. Two-dimensional vertical moisture-pressure dynamics above groundwater waves: Sand flume experiments and modelling. J. Hydrol. 2017, 544, 467–478. [Google Scholar] [CrossRef]
  36. Anochikwa, C.; Kamp, G.; Barbour, S. Interpreting pore-water pressure changes induced by water table fluctuations and mechanical loading due to soil moisture changes. Can. Geotech. J. 2012, 49, 357–366. [Google Scholar] [CrossRef]
  37. Tatard, L.; Planchon, O.; Wainwright, J.; Nord, G.; Favis-Mortlock, D.; Silvera, N.; Ribolzi, O.; Esteves, M.; Huang, C. Measurement and modelling of high-resolution flow-velocity data under simulated rainfall on a low-slope sandy soil. J. Hydrol. 2008, 348, 1–12. [Google Scholar] [CrossRef]
  38. Shoushtari, S.; Nielsen, P.; Cartwright, N.; Perrochet, P. Periodic seepage face formation and water pressure distribution along a vertical boundary of an aquifer. J. Hydrol. 2015, 523, 24–33. [Google Scholar] [CrossRef]
  39. Cuss, R.; Harrington, J.; Graham, C.; Noy, D. Observations of pore pressure in clay-rich materials; implications for the concept of effective stress applied to unconventional hydrocarbons. Energy Procedia 2014, 59, 59–66. [Google Scholar] [CrossRef]
  40. Shu, A.; Zhou, X.; Yu, M.; Duan, G.; Zhu, F. Characteristics for circulating currents and water-flow shear stress under the condition of bank slope collapse. J. Hydraul. Eng. 2018, 49, 271–281. [Google Scholar]
  41. Lu, J.; Lee, J.; Lu, T.; Hong, J. Experimental study of extreme shear stress for shallow flow under simulated rainfall. Hydrol. Process. 2009, 23, 1660–1667. [Google Scholar] [CrossRef]
  42. Onitsuka, K.; Akiyama, J.; Matsuoka, S. Prediction of velocity profiles and Reynolds stress distributions in turbulent open-channel flows with adverse pressure gradient. J. Hydraul. Res. 2009, 47, 58–65. [Google Scholar] [CrossRef]
  43. Zheng, F.; Huang, C.; Norton, L. Vertical hydraulic gradient and run-on water and sediment effects on erosion processes and sediment regimes. Soil Sci. Soc. Am. J. 2000, 64, 4–11. [Google Scholar] [CrossRef]
  44. Green, C.; Stonestrom, D.; Bekins, B.; Akstin, K.; Schulz, M. Percolation and transport in a sandy soil under a natural hydraulic gradient. Water Resour. Res. 2005, 41, W10414. [Google Scholar] [CrossRef]
  45. Liu, H.; Birkholzer, J. On the relationship between water flux and hydraulic gradient for unsaturated and saturated clay. J. Hydrol. 2012, 475, 242–247. [Google Scholar] [CrossRef]
  46. Newman, B.; Campbell, A.; Wilcox, B. Tracer-based studies of soil water movement in semi-arid forests of New Mexico. J. Hydrol. 1997, 196, 251–270. [Google Scholar] [CrossRef]
  47. Harr, R. Water flux in soil and subsoil on a steep forested slope. J. Hydrol. 1977, 33, 37–58. [Google Scholar] [CrossRef]
Figure 1. The energy distribution of the moving water in pipeline flow.
Figure 1. The energy distribution of the moving water in pipeline flow.
Water 16 02712 g001
Figure 2. The energy distribution when the subsurface flow moves downward in the soil pipe.
Figure 2. The energy distribution when the subsurface flow moves downward in the soil pipe.
Water 16 02712 g002
Figure 3. Force analysis of soil water in soil pipeline when preferential flow occurs.
Figure 3. Force analysis of soil water in soil pipeline when preferential flow occurs.
Water 16 02712 g003
Figure 4. The five parts of typical Benggang erosion and the location of the two research soil profiles (A and B) in the hillslope.
Figure 4. The five parts of typical Benggang erosion and the location of the two research soil profiles (A and B) in the hillslope.
Water 16 02712 g004
Figure 5. The soil profile of soil moisture and pore water pressure monitoring sites (naturally exposed in the outside). The (A) profile has no crack; the (B) profile has an obvious crack.
Figure 5. The soil profile of soil moisture and pore water pressure monitoring sites (naturally exposed in the outside). The (A) profile has no crack; the (B) profile has an obvious crack.
Water 16 02712 g005
Figure 6. The water flow shear force changes (df)with soil water content (ac) when matrix flow (MF) occurs in each soil layer.
Figure 6. The water flow shear force changes (df)with soil water content (ac) when matrix flow (MF) occurs in each soil layer.
Water 16 02712 g006
Figure 7. The water flow shear force (df) changes and soil pore water pressure (gi) changes with soil water content (ac) when preferential flow (PF) occurs in each soil layer.
Figure 7. The water flow shear force (df) changes and soil pore water pressure (gi) changes with soil water content (ac) when preferential flow (PF) occurs in each soil layer.
Water 16 02712 g007
Figure 8. The variation in pore water pressure with soil water content in matrix flow (a) and preferential flow (b).
Figure 8. The variation in pore water pressure with soil water content in matrix flow (a) and preferential flow (b).
Water 16 02712 g008
Table 1. Soil physical properties of four soil layers at two profiles.
Table 1. Soil physical properties of four soil layers at two profiles.
ProfileSoil Layer (cm)Soil Particle Size Distribution (%)Soil Porosity (%)Bulk Density
(g cm−3)
Sand *SiltClay>1 mm1~0.5 mm0.5~0.2 mm<0.2 mm
A0~2043.6632.1924.151.521.161.7638.591.58
20~4043.3332.3824.291.513.61.5236.771.54
40~6044.235.8719.931.263.360.9334.21.55
60~8044.6535.5819.771.32.980.7434.61.54
B0~2047.6732.6919.641.63.061.1135.721.55
20~4050.5931.8917.522.956.273.2728.611.56
40~6056.9729.0713.965.455.91.0129.941.53
60~8058.9426.914.169.7210.170.5422.171.52
Notes: * The size of sand, silt, and clay were 2~0.05 mm, 0.05~0.002 mm, and less than 0.002 mm, respectively (according to USAD).
Table 2. The equivalent pore diameter (the equivalent pore diameter for PF refers to the width of the fracture of the soil layer where the PF occurs), infiltration rate, and the Reynolds number due to pore pressure when MF and PF occurred during the initial stage of infiltration at each profile and each soil layer.
Table 2. The equivalent pore diameter (the equivalent pore diameter for PF refers to the width of the fracture of the soil layer where the PF occurs), infiltration rate, and the Reynolds number due to pore pressure when MF and PF occurred during the initial stage of infiltration at each profile and each soil layer.
ProfileRainfall AmountRainfall IntensitySoil LayerInitial Water
Content
Equivalent
Pore Diameter
(PF/MF)
Infiltration Rate
(PF/MF)
Reynolds
Number
(MF)
Reynolds
Number
(PF)
mmmm h−1cmm3 m−3mmmm s−1
A60, artificial rainfall100~200.2220.20
20~400.3360.320.0330.033
40~600.2670.290.0430.029
60~800.2780.310.0130.014
60, artificial rainfall100~200.2540.25
20~400.3540.340.0460.049
40~600.3010.340.0180.021
60~800.3020.330.0160.016
34, natural rainfall100~200.2430.23
20~400.3730.350.0830.091
40~600.2910.320.0330.034
60~800.3190.340.0150.016
B30, artificial rainfall100~200.2842.11
20~400.2942.341.334 0.975
40~600.1961.981.334 0.825
60~800.1872.022.667 1.684
30, artificial rainfall100~200.2152.11
20~400.2672.341.334 0.975
40~600.1741.982.667 1.652
60~800.2662.022.667 1.684
27, natural rainfall100~200.3812.11
20~400.3732.341.334 0.975
40~600.3441.981.334 0.825
60~800.2832.022.667 1.684
Table 3. The water flow shear force when MF and PF occurred during the initial stage of infiltration at each profile and each soil layer.
Table 3. The water flow shear force when MF and PF occurred during the initial stage of infiltration at each profile and each soil layer.
ProfileSoil Layer
(cm)
Pore Water Pressure
(KPa)
Soil Matrix Potential (MF)
(KPa)
Water Flow Shear Force
(N m−2)
Water Flow Shear Force (N m−2)
via Pore
Pressure
via Gravityvia Matrix
Potential
A, artificial rainfall0–200.016−2.3
20–400.014−2.10.5470.0010.591−0.045
40–600.011−5.21.5790.0010.5641.014
60–800.052−4.60.343−0.0120.588−0.233
A, artificial rainfall0–200.346−2.2
20–400.008−1.60.4610.1640.711−0.414
40–600.011−3.11.216−0.0010.7110.506
60–800.076−3.00.648−0.0240.709−0.039
A, natural rainfall0–200.014−2.1
20–40−0.279−1.20.5930.1260.761−0.294
40–60−0.292−4.21.7290.0040.6621.063
60–800.077−2.80.076−0.1430.760−0.543
B, artificial rainfall0–200.842
20–401.0536 4.439−1.3355.737
40–600.161 6.9022.0474.855
60–800.277 4.576−0.3844.960
B, artificial rainfall0–20−0.561
20–400.257 2.406−3.3325.737
40–600.341 3.623−1.2324.855
60–80−0.357 6.9491.9894.960
B, natural rainfall0–20−0.059
20–400.231 5.621−0.1165.737
40–600.115 5.1210.2664.855
60–800.091 5.046−0.0864.960
Table 4. The water flow shear force calculated via Darcy’s law when MF occurred during the initial stage of infiltration at each profile and each soil layer.
Table 4. The water flow shear force calculated via Darcy’s law when MF occurred during the initial stage of infiltration at each profile and each soil layer.
ProfileSoil Layer
(cm)
Equivalent Pore Diameter (MF)
(mm)
Saturated Hydraulic Conductivity
(cm d−1)
Water Flow Flux
(mm s−1)
Water Flow Shear Force Due to Pore Pressure
(N m−2)
B, artificial rainfall0–200.30.163
20–400.320.3530.0330.907
40–600.230.3490.0430.252
60–800.240.5770.0130.147
B, artificial rainfall0–200.210.163
20–400.340.3530.0460.631
40–600.290.3490.0180.766
60–800.290.5770.0161.419
B, natural rainfall0–200.210.163
20–400.350.3530.0831.606
40–600.270.3490.0330.465
60–800.310.5770.0150.132
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Tao, Y.; Peng, S.; Chen, J.; Long, S.; Liao, B. Exploring the Effects of Fissures on Hydraulic Parameters in Subsurface Flows from the Perspective of Energy Changes. Water 2024, 16, 2712. https://doi.org/10.3390/w16192712

AMA Style

Tao Y, Peng S, Chen J, Long S, Liao B. Exploring the Effects of Fissures on Hydraulic Parameters in Subsurface Flows from the Perspective of Energy Changes. Water. 2024; 16(19):2712. https://doi.org/10.3390/w16192712

Chicago/Turabian Style

Tao, Yu, Siwen Peng, Jiazhou Chen, Shiping Long, and Bin Liao. 2024. "Exploring the Effects of Fissures on Hydraulic Parameters in Subsurface Flows from the Perspective of Energy Changes" Water 16, no. 19: 2712. https://doi.org/10.3390/w16192712

APA Style

Tao, Y., Peng, S., Chen, J., Long, S., & Liao, B. (2024). Exploring the Effects of Fissures on Hydraulic Parameters in Subsurface Flows from the Perspective of Energy Changes. Water, 16(19), 2712. https://doi.org/10.3390/w16192712

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