Understanding natural flows in porous media with low Reynolds number (Re) has significant implications for both science and engineering. However, knowledge and experimental gaps remain regarding such natural flows. In this context, we designed a sand column-based laboratory filtration experiment to investigate flow characteristics in porous media with low Re. Four media were considered including two silts (silt-I and silt-II), one silty sand, and one medium sand. Results show that constant hydraulic head is presumed to be an important factor that affects flow regime in porous media. In general, the flow approaches Darcian at a constant hydraulic head of ~170 cm, whereas it becomes non-Darcian at a constant hydraulic head of ~230 cm. The type of media determines the Re range that delimitates between Darcy and non-Darcy flows. Specifically, the transition appears at 0.015 < Re < 0.020 for medium sand and 0.000027 < Re < 0.000029 for silt-II, respectively. In the condition of constant hydraulic heads, the breakdown of Darcy’s law may occur at very low Re values (Re→0). Media dependent Re ranges are probably needed to describe the beginning of non-Darcy flows, rather than 1 to 10 or other value for all media. Findings in this study can offer insights into calculation and simulation of flows in low-permeability reservoirs, pumping process of foundation pit excavation, and other non-Darcy flows in low-permeability media.

Knowledge of dam construction in floodplain systems and its hydrodynamic effects plays a critical role in managing various kinds of floodplains. This study uses 3D floodplain hydrodynamic modeling to explore the possible effects of a proposed hydraulic project in Poyang Lake (PLHP) on the hydrodynamics, exemplified by a large floodplain system. Simulations showed that the water levels across most lake regions presented more significant changes than in the floodplain areas during the study period. The increased water levels upstream from the PLHP (~1.0 m) were distinctly higher than that downstream (~0.1 m). The PLHP may decrease the magnitude of the water velocities in the main channels of the lake, whereas velocities may experience mostly minor changes in the floodplains, depending upon the altered flow dynamics and transport. On average, the water temperature may exhibit mostly minor changes (~<1.0 °C) for both the horizontal and vertical scales within the flood-pulse-influenced lake system. Additionally, the model results indicated that the outflow process caused by the PLHP may be altered from the natural discharge into the Yangtze River to frequent backflow events during the storage period, demonstrating the non-negligible effect of the PLHP on the water supply for the downstream Yangtze River in the future.

The baseflow of the Yellow River is vital and important for water resource management and for understanding the hydrological cycle and ecohydrology setting in this arid and semi-arid basin. This study uses a Lyne and Hollick digital filtering technique to investigate the behaviors of the baseflow and the baseflow index in the upper reaches of the Yellow River Basin (China). The observed streamflow discharges along the river were used to analyze the baseflow trend, persistence, and periodic characteristics during the period of 1950–2000. The results show that the average baseflow and BFI in the upper reaches of the Yellow River exhibit a decreasing trend and will continue to decline in the future. Generally, the annual average baseflow and BFI for the most upstream areas of the Yellow River show little difference, while the baseflow and BFI exhibit significant differences for the downstream areas. The filtered annual baseflow varied between 128 × 10⁸ m³/year and 193 × 10⁸ m³/year for the Yellow River. The BFI ranged from 0.54 to 0.65, with an average of 0.60. This indicates that on average, 60% of the long-term streamflow is likely controlled by groundwater discharge and shallow subsurface flow. Statistics show that two periodic variations were observed in the baseflow evolution process. The results indicate that on average, the first and second main cycles of baseflow behaviors occur at 28 years and 12–17 years, respectively. Correspondingly, the estimation indicates that the abrupt change points tend to appear in the 1960s, the 1980s, and the 1990s. An improved understanding of baseflow behaviors can help guide future strategies to manage the river regime, its water resources, and water quality.

Light non-aqueous phase liquids (LNAPLs), which include various petroleum products, are a significant source of groundwater contamination globally. Once introduced into the subsurface, these contaminants tend to accumulate in the vadose zone, causing chronic soil and water pollution. The vadose zone often contains lens-shaped bodies with diverse properties that can significantly influence the migration and distribution of LNAPLs. Understanding the interaction between LNAPLs and these lens-shaped bodies is crucial for developing effective environmental management and remediation strategies. Prior research has primarily focused on LNAPL behavior in homogeneous media, with less emphasis on the impact of heterogeneous conditions introduced by lens-shaped bodies. To investigate the impact of lens-shaped structures on the migration of LNAPLs and to assess the specific effects of different types of lens-shaped structures on the distribution characteristics of LNAPL migration, this study simulates the LNAPL leakage process using an indoor two-dimensional sandbox. Three distinct test groups were conducted: one with no lens-shaped aquifer, one with a low-permeability lens, and one with a high-permeability lens. This study employs a combination of oil front curve mapping and high-density resistivity imaging techniques to systematically evaluate how the presence of lens-shaped structures affects the migration behavior, distribution patterns, and corresponding resistivity anomalies of LNAPLs. The results indicate that the migration rate and distribution characteristics of LNAPLs are influenced by the presence of a lens in the gas band of the envelope. The maximum vertical migration distances of the LNAPL are as follows: high-permeability lens (45 cm), no lens-shaped aquifer (40 cm), and low-permeability lens (35 cm). Horizontally, the maximum migration distances of the LNAPL to the upper part of the lens body decreases in the order of low-permeability lens, high-permeability lens, and no lens-shaped aquifer. The low-permeability lens impedes the vertical migration of the LNAPL, significantly affecting its migration path. It creates a flow around effect, hindering the downward migration of the LNAPL. In contrast, the high-permeability lens has a weaker retention effect and creates preferential flow paths, promoting the downward migration of the LNAPL. Under conditions with no lens-shaped aquifer and a high-permeability lens, the region of positive resistivity change rate is symmetrical around the axis where the injection point is located. Future research should explore the impact of various LNAPL types, lens geometries, and water table fluctuations on migration patterns. Incorporating numerical simulations could provide deeper insights into the mechanisms controlling LNAPL migration in heterogeneous subsurface environments.

Contamination from light nonaqueous phase fluids (LNAPLs) and their derivatives during mining, production, and transportation has become a concern. Scholars have extensively studied LNAPL contamination, but the role of water content variation on its migration process in the unsaturated zone has not been sufficiently researched. The specific issue addressed in this study is the impact of water content on the migration of light nonaqueous phase liquids (LNAPLs) in sandy soils, a critical yet under-researched aspect of subsurface contamination. To tackle this, we employed indoor simulated vertical, one-dimensional, multiphase flow soil column experiments, utilizing the orthogonal experimental method to systematically evaluate the effects of varying water contents on the occurrence state and migration rate of LNAPLs. The experimental results indicate the following: (1) The migration rate of LNAPL exhibits an L-shaped trend during subsurface imbibition and a nonlinear relationship with migration time. The migration rate and migration time of surface infiltration have a linear growth relationship. (2) The residual rate of LNAPL is negatively correlated with water content and positively correlated with oil content in the homogeneous non-saturated state. With the increase in the amount of leaked oil, 40% of the leaked LNAPL is sorbed within the soil. (3) When the water content of the test medium is below 14%, and the oil content is below 11%, LNAPL appears in the unsaturated zone in a solid phase. As the water content increases, the adsorption rate of the oil phase gradually decreases and eventually reaches the oil saturation point. (4) When the water content of the medium exceeds 8%, over time, LNAPL will be subject to oil–water interfacial tension, and the rate of LNAPL movement first decreases and then increases, displaying nonlinear growth. The innovation of this work lies in the comprehensive analysis of LNAPL migration under controlled laboratory conditions, providing results that enhance the understanding of LNAPL behavior in sandy soils. These quantitative insights are crucial for developing targeted remediation strategies for LNAPL-induced pollution in the unsaturated zone.