Search

• 3D pathways of IGF are identified with an assessment on flow rates. • IGF contributes 36%∼59% to the baseflow and 32%∼52% to the streamflow in the studied river. • Infiltration recharge and aquifer properties affect the development of IGF. Inter-basin groundwater flow (IGF), defined as groundwater flows across the topographic divides, occurs in nature as a common phenomenon, and seriously affect water balance and solute transport in a basin. Although the existence of IGF has been checked in drainage basins with indirect methods, such as using information of hydrogeochemistry, few studies are conducted to directly identify and quantify the IGF in the three-dimensional (3D) space using numerical modeling of groundwater flow. In this study, we identify, partition and quantify IGF based on a regional-scale 3D groundwater flow model covering several traditional surface water basins. The investigation is focused on the Hailiutu River, a groundwater-fed stream in the Ordos Plateau, China. It is the first time that the source zones and 3D paths of IGF as well as its contribution to the streamflow are fully captured for a groundwater-fed stream. Results indicate that the IGF contributes 36%∼59% to the baseflow and 32%∼52% to the streamflow of the Hailiutu River, through different groups of circulation paths with the maximum depth up to ∼ 1000 m. Uncertainty analyses of the model show that the contribution of IGF to the streamflow is positively correlated to the infiltration recharge and hydraulic conductivities of the aquifer media. Several indirect methods are also checked for comparison, but they could only be used to assess the potential contribution of IGF, not the IGF paths. This study will facilitate in deeply understanding the evolution mechanism of the inter-basin groundwater circulation and the interaction between surface water and groundwater. [ABSTRACT FROM AUTHOR]

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]

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 coupled regional and local model is required when groundwater flow and solute transport are to be simulated in local areas of interest with a finer grid while regional aquifer boundary and major stresses should be retained with a coarser grid. The coupled model should also maintain interactions between the regional and local flow systems. In the Beijing Plain (China), assessment of managed aquifer recharge (MAR), groundwater pollution caused by rivers, capture zone of well fields, and land subsidence at the cone of depression requires a coupled regional and local model. This study evaluates three methods for coupling regional and local flow models for simulating MAR in the Chaobai River catchment in the Beijing Plain. These methods are the conventional grid refinement (CGR) method, the local grid refinement (LGR) method and the unstructured grid (USG) method. The assessment included the comparison of the complexity of the coupled model construction, the goodness of fit of the computed and observed groundwater heads, the consistency of regional and local groundwater budgets, and the capture zone of a well filed influenced by the MAR site. The results indicated that the CGR method based on MODFLOW-2005 is the easiest to implement the coupled model, capable of reproducing regional and local groundwater heads and budget, and already coupled with density and viscosity dependent model codes for transport simulation. However, the CGR method inherits shortcomings of finite difference grids to create multiple local models with inefficient computing efforts. The USG method based on MODFLOW-USG has the advantage of creating multi-scale models and is flexible to simulate rivers, wells, irregular boundaries, heterogeneities and the MAR site. However, it is more difficult to construct the coupled models with the unstructured grids, therefore, a good graphic user interface is necessary for efficient model construction. The LGR method based on MODFLOW-LGR can be used to create multiple local models in uniform aquifer systems. So far, little effort has been devoted to upgrade the LGR method for complex aquifer structures and develop the coupled transport models. [ABSTRACT FROM AUTHOR]

• Sustainable groundwater development is achievable in Beijing plain. • Substituting groundwater abstraction with transferred surface water has stopped depletion. • Managed aquifer recharge is very effective for sustaining groundwater abstraction. • Climate variation has larger impact on groundwater resources than the anticipated future climate change. To cope with the groundwater depletion problem and achieve sustainable groundwater development, groundwater conservation measures and managed aquifer recharge (MAR) have been implemented worldwide. However, knowledge gaps exit how does the aquifer system respond to these interventions differently and if these interventions are adequate to lead to long-term sustainable groundwater development under future climate change. In Beijing Plain, two measures have been implemented: reduction of groundwater abstraction by substituting groundwater abstraction with transferred surface water and implementation of managed aquifer recharge (MAR) in two major rivers. This study aims to assess how do the shallow and deep aquifers respond to these measures and if these measures can lead to long-term sustainable groundwater development in Beijing Plain under future climate change. A 3-D transient groundwater flow model was calibrated and used to simulate groundwater level and budget changes from 2021 to 2050. The monthly groundwater recharge was estimated using the projected monthly precipitation from three downscaled regional climate models under two scenarios (RCP4.5 and RCP8.5). The results show that declines in groundwater head and storage can be reversed with the combined two measures, thereby contributing to achieve sustainable groundwater development. The reduction of abstractions is a deciding measure to reverse the trend of groundwater depletion, especially in the deep confined aquifers, while large scale MAR schemes can restore the cones of depressions in shallow aquifers and maintain the groundwater abstraction. Climate variation has large impacts on groundwater resources, especially, consecutive dry years can cause rapid groundwater storage depletion. The projected monthly precipitation from 2021 to 2050 is not significantly different from the past. Therefore, the projected future precipitation has minor impacts on groundwater resources in the next 30 years. The findings from the study will support the Beijing municipality to maintain the tight control on groundwater abstraction and to implement large-scale MAR schemes in two rivers. This successful example will encourage managers of other heavily exploited aquifers to take similar measures to achieve sustainable groundwater development. [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]

• Groundwater mounds resulting from nearshore recharge is identified in a sandy hillslope. • The thickness of vadose zone determines fluxes across water table along the hillslope. • Nearshore recharge temporarily changes streamlines and increases tortuosity of flowpaths. • Flowpaths near the wetland are mainly controlled by regional GW discharge and strong ET. Groundwater mounds and depressions induced by rainfall and evapotranspiration near wetlands/lakes have been widely observed, but their effects on groundwater flow field are poorly known. Here, a sandy hillslope near a wetland in the Ordos Plateau, China, which has a semi-arid climate, is selected for investigation. After identifying temporary groundwater mounds/depressions based on three water level observation wells, a HYDRUS-2D model is setup to obtain the highly transient daily flow field. The variable thickness of unsaturated zone near the wetland, which controls soil water content and unsaturated hydraulic conductivity, is found to be responsible for groundwater mounds resulting from focused recharge and groundwater depressions resulting from nonuniform evapotranspiration. On nonrainy days, evapotranspiration in the lowland is supported by flowpaths from regional groundwater discharge, and in the slope is supported by a flow cell from the upland to the upper part of the slope, which causes water table recession in the upland. On rainy days, focused recharge in the lowland leads to two temporary local flow cells, which break the regional groundwater discharge flowpaths into two temporary flow cells. The two temporary flow cells toward the upland contribute to water table rises in the uplandMoreover, we find focused recharge in rainy days does not change the overall direction of flowpaths controlled by regional groundwater discharge and strong evapotranspiration, but only increase the tortuosity of flowpaths. By coupling daily atmospheric conditions and regional groundwater discharge, this study has implications for groundwater recharge estimation and future groundwater modelling in arid regions. [ABSTRACT FROM AUTHOR]

• Three alternative conceptual models were constructed. • Each model was calibrated automatically and evaluated statistically. • Hydrogeological interpretations are the basis for the conceptualisation. • Modelling objectives determine the complexity of the conceptual models. Three alternative groundwater flow models were evaluated for Beijing Plain, China. The first model (AM1) was constructed with the "thin layer approach" in which all 9 model layers, including five aquifers separated by four aquitards, are continuously present in the same model area. The second model (AM2) was constructed with the "quasi-3D approach" in which the hydrogeological formations were classified into five aquifer units consisting of mixed permeable and semi-permeable layers at different depth ranges. The third model (AM3) was constructed with the "true layer approach" in which aquifers and aquitards were defined according to hydrostratigraphic properties, and model layers are absent in the area where corresponding hydrogeological formations intersect bedrocks. All 3 models were calibrated with the parameter optimization method under the steady state flow condition with the same hydrological stresses and observation data. All three models fit to observations well with the similar calibration criteria values. Furthermore, AIC and BIC information criteria could not distinguish three alternative models. Only KIC could identify AM3 as the best model. Major differences of the three alternative models were identified from a hydrogeological perspective. The AM1 model depicted an illusion through contour maps that groundwater was present everywhere in the deep aquifers. The model computed larger vertical leakages because more abstraction rates were assigned improperly in deep aquifers. The AM2 model was able to compute regional groundwater balances and depicted spatial groundwater level variations. However, the AM2 model computed longer groundwater travel times around the wellfield and should not be used for the delineation of the well field protection zones and contaminant transport simulation. The AM3 model could not only compute the regional groundwater balances and describe spatial groundwater distribution in deep confined aquifers, but also delineate the capture zone of the wellfield. It can therefore be used for simulating contaminant transport. Furthermore, the AM3 model is suitable to construct a coupled regional and local flow model for simulating a managed aquifer recharge scheme. [ABSTRACT FROM AUTHOR]