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]

• 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]

• 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]

Summary: Multi-scale heterogeneity of geological formations is a rule, which consists of random (local-scale) and systematic (large-scale) variability of hydraulic conductivity. The random variability and depth–decaying trend, a systematic variability, have different effects on subsurface flow, thus on groundwater discharge into tunnels. Little research has examined this problem in the past. Using Monte Carlo simulation and information of statistics of heterogeneity, we evaluate the most likely (ensemble average) discharge rate into a tunnel in geologic media with the multi-scale heterogeneity and uncertainty associated with this estimate. We find that the ensemble average discharge rate is larger than the discharge rate predicted by geometric mean of hydraulic conductivity, and smaller than the discharge rate predicted by arithmetic mean of hydraulic conductivity. Moreover, the ensemble average discharge rate decreases with the decay exponent of the depth–decaying trend, and increases with the standard deviation as well as the correlation scale of the stationary log-conductivity fields. We also find that the discharge rate of a tunnel is highly influenced by the hydraulic conductivity near the tunnel. Furthermore, deviation of the true discharge rate from the ensemble average can be large and increases with the decay exponent, standard deviation and correlation scale of log-conductivity fields. The largest uncertainty of discharge rate prediction in the shallow subsurface is controlled by the variability of conductivity fields and the uncertainty at the deep subsurface is by the depth–decaying trend of hydraulic conductivity. Therefore, accurate prediction of groundwater discharge into tunnels requires detailed characterization of multi-scale heterogeneity. [Copyright &y& Elsevier]