Search

Hierarchically nested groundwater flow systems have been widely investigated by numerical modeling and laboratory experiments but seldom recognized from in situ observations and testing. Groundwater age and geochemical profiles were obtained by the authors using a packer system along two wellbores drilled in the Ordos Plateau, China. Groundwater age profiles were constructed from 3H, 85Kr, 14C, and 81Kr data, with the apparent age generally increasing with depth, while exhibiting several inflection points. These inflection points were compared with geochemical profiles and numerical modeling results to delineate groundwater flow systems. Inflection points with significant increases in the gradients of age and Cl− concentration versus depth may indicate interfaces among flow systems. A concurrent decrease‐increase turning of groundwater age and Cl− concentration is influenced by a stagnation point. The groundwater age profiles are quite different from previous studies, indicating that the method can greatly reduce the uncertainty in model construction. Plain Language Summary: The provision of renewable groundwater, a crucial source of drinking water worldwide, relies on a circulation path from recharge to discharge zones with a subsurface flow in reservoirs (aquifers). Basin‐scale groundwater circulation is speculated to form hierarchically nested flow systems with differing capabilities for renewal. However, flow systems are invisible and difficult to detect through in situ observations on groundwater. In this study, nested groundwater flow systems were identified from groundwater age and geochemical profiles along wellbores. Inflection points on these profiles provide significant insights into the interfaces and stagnation points among nested groundwater flow systems. Uncertainties could be brought by the aquifer heterogeneity in the hydraulic conductivity and localized soluble salts. This adopted method is expected to be improved by incorporating geophysical survey and high‐resolution geochemical analyses to eliminate uncertainties in model construction. Key Points: Groundwater age and geochemical profiles were obtained to investigate groundwater flow systemsInflection points on profiles are used to infer interfaces or stagnation points among different groundwater flow systemsGroundwater flow models constrained by inflection points greatly reduce the uncertainty of flow system characterization [ABSTRACT FROM AUTHOR]

Freezing-induced groundwater-level decline is widely observed in regions with a shallow water table, but many existing studies on freezing-induced groundwater migration do not account for freezing-induced water-level fluctuations. Here, by combining detailed field observations of liquid soil water content and groundwater-level fluctuations at a site in the Ordos Plateau, China, and numerical modeling, we showed that the interaction of soil water and groundwater dynamics was controlled by wintertime atmospheric conditions and topographically driven lateral groundwater inflow. With an initial water table depth of 120 cm and a lateral groundwater inflow rate of 1.03 mmd-1 , the observed freezing and thawing-induced fluctuations of soil water content and groundwater level are well reproduced. By calculating the budget of groundwater, the mean upward flux of freezing-induced groundwater loss is 1.46 mmd-1 for 93 d , while the mean flux of thawing-induced groundwater recharge is as high as 3.94 mmd-1 for 32 d. These results could be useful for local water resources management when encountering seasonally frozen soils and for future studies on two- or three-dimensional transient groundwater flow in semi-arid and seasonally frozen regions. By comparing models under a series of conditions, we found the magnitude of freezing-induced groundwater loss decreases with initial water table depth and increases with the rate of groundwater inflow. We also found a fixed-head lower boundary condition would overestimate freezing-induced groundwater migration when the water table depth is shallow. Therefore, an accurate characterization of freezing-induced water table decline is critical to quantifying the contribution of groundwater to hydrological and ecological processes in cold regions. [ABSTRACT FROM AUTHOR]

The groundwater divide is a key feature of river basins and significantly influenced by subsurface hydrological processes. For an unconfined aquifer between two parallel rivers or ditches, it has long been defined as the top of the water table based on the Dupuit–Forchheimer approximation. However, the exact groundwater divide is subject to the interface between two local flow systems transporting groundwater to rivers from the infiltration recharge. This study contributes a new analytical model for two‐dimensional groundwater flow between rivers of different water levels. The flownet is delineated in the model to identify groundwater flow systems and the exact groundwater divide. Formulas with two dimensionless parameters are derived to determine the distributed hydraulic head, the top of the water table and the groundwater divide. The locations of the groundwater divide and the top of the water table are not the same. The distance between them in horizontal can reach up to 8.9% of the distance between rivers. Numerical verifications indicate that simplifications in the analytical model do not significantly cause misestimates in the location of the groundwater divide. In contrast, the Dupuit–Forchheimer approximation yields an incorrect water table shape. The new analytical model is applied to investigate groundwater divides in the Loess Plateau, China, with a Monte Carlo simulation process taking into account the uncertainties in the parameters. [ABSTRACT FROM AUTHOR]

The Badain Jaran Desert (BJD) in China is a desert with impressive sand dunes and a groundwater situation that has attracted numerous researchers. This paper gives an overview of the mysteries of groundwater in the BJD that are exhibited as five key problems identified in previous studies. These problems relate to the origin of the groundwater, the hydrological connection between the BJD and the Heihe River Basin (HRB), the infiltration recharge, the lake-groundwater interactions, and the features of stable isotope analyses. The existing controversial analyses and hypotheses have caused debate and have hindered effective water resources management in the region. In recent years, these problems have been partly addressed by additional surveys. It has been revealed that the Quaternary sandy sediments and Neogene-Cretaceous sandstones form a thick aquifer system in the BJD. Groundwater flow at the regional scale is dominated by a significant difference in water levels between the surrounding mountains and lowlands at the western and northern edges. Discharge of groundwater from the BJD to the downstream HRB occurs according to the regional flow. Seasonal fluctuations of the water level in lakes are less than 0.5 m due to the quasi-steady groundwater discharge. The magnitude of infiltration recharge is still highly uncertain because significant limitations existed in previous studies. The evaporation effect may be the key to interpreting the anomalous negative deuterium-excess in the BJD groundwater. Further investigations are expected to reveal the hydrogeological conditions in more detail. [ABSTRACT FROM AUTHOR]

A new method based on a graphics processing unit ( GPU) library is proposed in the paper to parallelize MODFLOW. Two programs, GetAb_CG and CG_GPU, have been developed to reorganize the equations in MODFLOW and solve them with the GPU library. Experimental tests using the NVIDIA Tesla C1060 show that a 1.6- to 10.6-fold speedup can be achieved for models with more than 105 cells. The efficiency can be further improved by using up-to-date GPU devices. [ABSTRACT FROM AUTHOR]

Unconfined aquifers beneath piedmont pluvial fans are widely distributed in front of mountains and proper for water supply with pumping wells. However, the catchment zone and capture zones of a pumping well in such an unconfined aquifer is not well known. We develop a preliminary simplified model where groundwater flows between a segmental inflow boundary and a discharge boundary of constant head. The catchment zone is delineated from numerical simulation via MODFLOW and MODPATH. Results are expressed with dimensionless variables and lumped parameters to show general behaviors. Sensitive analyses indicate that there are 4 types of the catchment zone according to different connections to the boundaries. The shape of the catchment zone is quantitatively analyzed with typical shape factors. Capture zones with respect to special travel times are identified from travel time distribution in the catchment zone. The modeling results can be applied in the design of water supply wells and delineation of protection zones at a site with similar hydrogeological conditions. [ABSTRACT FROM AUTHOR]

• We summarize field observations of 14C age showing obvious GW aging in two basins. • Pumping causes not only rejuvenation, but also aging of GW in layered basins. • Although GW rejuvenation is larger in area, aging is greater in magnitude. • We identify two new mechanisms of pumping-induced GW aging in layered basins. • Topography-driven flow is the main cause of GW aging in the discharge area. As an indicator of the renewability of a groundwater basin, the spatial distribution of groundwater age has attracted intensive research. As reported in the literature, 3H measurements in 74 US aquifers showed pumping-induced groundwater rejuvenation, whereas 14C ages in two heavily developed aquifers in France and China indicated pumping-induced groundwater aging. Unfortunately, the mechanisms leading to the seemingly opposite behavior of pumping-induced aging and rejuvenation remain unclear. In this study, we conducted numerical simulations of groundwater age in layered aquifer-aquitard systems under different pumping scenarios and aquitard permeability to show that pumping-induced groundwater rejuvenation and aging always occur simultaneously. Groundwater rejuvenation widely occurs throughout the basin wherever the well is located, whereas groundwater aging mainly occurs to the downgradient parts of the basin. Although groundwater rejuvenation is larger in area, groundwater aging is greater in magnitude, and the magnitude of aging increases with pumping intensity and duration. Moreover, when pumping occurs in the discharge area, aging could occur in almost the entire aquifer to the downgradient of the pumping well, and the magnitude of aging is larger. The large area and large magnitude of aging in the discharge area is a direct result of upward leakage of old groundwater from underlying aquitards/aquifers, which is controlled by the topography-driven groundwater flow. By demonstrating the importance of topography-driven groundwater flow and underlying formations in groundwater age redistribution, this study has implications for determining the bottom of a groundwater flow and solute transport model and interpretation of various tracers. [ABSTRACT FROM AUTHOR]

• Field data show the control of soil texture on freezing-induced water table recession. • Ku below the frozen zone is the limiting factor of freezing-induced water flow. • We obtain extinction depths of freezing-induced GW migration for different soils. • The extinction depths of freezing-induced and evaporation-induced GW migration. Although water table recessions associated with seasonal freezing have been widely reported, the conditions that groundwater could respond to seasonal freezing are poorly known. Based on 30 records of winter-time water table recessions reported in the literature, we find the mean water table recession increases from sand to loam to silt to clay. To obtain a quantitative understanding of the conditions of freezing-induced groundwater migration and the role of soil texture, we construct a one-dimensional SHAW model with a dynamic water table that could respond to winter-time atmospheric conditions. By simulating freezing-induced groundwater migration under a series of pre-freezing water table depths, soil textures and freezing indices, we find the extinction depth of freezing-induced groundwater migration, i.e., the limit water table depth beyond which there is negligible freezing-induced groundwater flow, increases from coarse to fine soils and increases with freezing index. Moreover, we find there is a good linear relationship between the distance from the freezing front to the water table corresponding to the extinction depth of freezing-induced groundwater migration and the distance from the land surface to the water table corresponding to the extinction depth of evaporation-induced groundwater migration for different soils. To demonstrate the importance of defining extinction depth, we also show the overestimation of freezing-induced water gain in the frozen zone by the widely used fixed water table condition when the water table depth is shallower than the extinction depth. This study has implications for future hydrologic and engineering studies in critical zones with a cold climate. [ABSTRACT FROM AUTHOR]

• We separate the roles of upward discharge and aquitards in flowing artesian zone. • Basinal flow are coupled with well hydraulics to examine hydraulics of flowing wells. • A low K aquitard causes a large flowing artesian zone but low stable discharge rate. • Baseflow reduction is the main source of groundwater derived in flowing wells. Causes and water sources of flowing artesian wells attracted the interest of many hydrogeologists throughout history, however, a quantitative model that satisfactorily considers the roles of topography, groundwater recharge/discharge and aquitards on hydraulics of flowing wells is still lacking. In this study, a three-layer river-valley basin with a recharge boundary is used to obtain the basinal flow field, spatial distribution and transient discharge rates of flowing wells. Even if there is an aquitard separating the unconfined and the confined aquifers, groundwater discharge to the river still plays a critical role in occurrence of flowing wells. An aquitard with lower permeability would enlarge the zone with flowing wells, implying that flowing wells are more likely to occur in basins with continuous aquitards. By comparing flowpaths and inter-aquifer leakages before and after numerically installing a flowing well in a three-layer basin, we find decreased upward leakage in the discharge limb plays a much larger role than increased downward leakage in the recharge limb in the stable discharge rate of a flowing well. This recognition is different from previous studies where the stable discharge rate is supported only by downward leakage from the overlying aquifer. By considering well hydraulics in a basinal background flow field, this study improves understanding of the mechanisms of flowing artesian wells and the interaction of surface water and deep groundwater, indicating that flowing wells are valuable for sampling deep groundwater representative of deep and old groundwater in the discharge area of a basin. [ABSTRACT FROM AUTHOR]