The impact of glaciation cycles on groundwater flow was studied within the framework of nuclear waste storage in underground geological formations. The eastern section of the Paris Basin (a layered aquifer with impervious/pervious alternations) in France was considered for the last 120 ka. Cold periods corresponded with arid climates. The issue of talik development below water bodies was addressed. These unfrozen zones can maintain open pathways for aquifer recharge. Transient thermal evolution was simulated on a small-scale generic unit of the landscape including a “river” and “plain”. Coupled thermo-hydraulic modeling and simplified conductive heat transfer were considered for a broad range of scenarios. The results showed that when considering the current limited river dimensions and purely conductive heat transfer, taliks are expected to close within a few centuries. However, including coupled advection for flows from the river to the plain (probably pertinent for the eastern Paris Basin aquifer recharge zones) strongly delays talik closure (millennium scale). The impact on regional underground flows is expected to vary from a complete stop of recharge to a reduced recharge, corresponding to the talik zones. Consequences for future modeling approaches of the Paris Basin are discussed.

The recent and rapid warming of the Arctic leads to thawing of permafrost, which influences and changes subsurface water-flow systems in such landscapes. This study explores the utility of catchments as “sentinels of change” by considering long-term discharge data from 17 stations on unregulated rivers in northern Sweden and analyzing trends in annual minimum discharge and recession flow characteristics. For the catchments considered, the annual minimum discharge has increased significantly (based on the Mann Kendall test at a 95 % confidence level) in nine of the catchments and decreased significantly in one catchment. Considering changes in recession-flow characteristics, seven catchments showed significant trends consistent with permafrost thawing while two catchments showed significant trends in the opposite direction. These results are mechanistically consistent with generic physically based modeling studies and the geological setting, as the catchments considered span the spatial limit of permafrost extent. This study illuminates the potential for using hydrologic observations to monitor changes in catchment-scale permafrost. Further, this opens the door for research to isolate the mechanisms behind the different trends observed and to gauge their ability to reflect actual permafrost conditions at the catchment scale.

Ranchi urban area (257 km (super 2) ) depends on aquifers for 30% of its total drinking-water supply of 17 million m (super 3) year (super -1) . Local hydrostratigraphy is represented by a heterogeneous, weathered and fractured aquifer system, typical of the Precambrian suite of rocks in the Indian subcontinent. Intensive development of the fractured aquifers, up to 200 m below ground, has lowered the hydraulic head and resulted in dwindling yields from fractures during the summer. To understand the groundwater flow regime and aquifer recharge mechanism, the present study examines delta (super 18) O and delta D variation in aquifer-specific samples along with water levels, yield of the fractures, EC and Cl (super -) . Three types of groundwater have been identified based on isotopic composition and d-excess values, each representing different recharge source-water and pathways. The major source of recharge for the aquifers is infiltration from rainfall. Two large reservoirs and an excavated lake within the study area contribute to the recharge process but insignificantly. Isotopic compositions and the relatively low EC and low Cl (super -) concentrations of high-yielding bore wells in some places indicate the presence of fast-conducting fracture zones receiving copious recharge from rainfall. Such fractures can be developed further through bore wells for drinking supply with due provision for artificial recharge.

Groundwater models simulating flow in buried valleys interacting with regional aquifers are often based on hydrogeological models interpreted from dense geophysical datasets and scarce borehole data. For three simple synthetic cases, it is demonstrated that alternative methods of inversion of transient electro-magnetic (TEM) data can lead to very different interpretations of the hydrogeology inside and surrounding a buried valley. The alternative interpreted hydrogeological models are used in numerical modelling of groundwater flow to a pumping well. It is demonstrated that the alternative models result in quite different groundwater-model predictions of capture zone, recharge area, and groundwater age for the pumping well. It is briefly demonstrated that model calibration against hydraulic head data is not likely to improve the predictions or to identify the structural error of the interpreted hydrogeological models. It is therefore concluded that when TEM-based resistivity models are interpreted to construct the hydrogeological framework of a groundwater model, it must be done cautiously with support from deep borehole information. Too much reliance on geophysical mapping can lead to seriously wrong hydrogeological models and correspondingly wrong groundwater-model predictions.

Water resources in agriculture-dominated basins of the arid western United States are stressed due to long-term impacts from pumping. A review of 88 regional groundwater-flow modeling applications from seven intensively irrigated western states (Arizona, California, Colorado, Idaho, Kansas, Nebraska and Texas) was conducted to provide hydrogeologists, modelers, water managers, and decision makers insight about past modeling studies that will aid future model development. Groundwater models were classified into three types: resource evaluation models (39%), which quantify water budgets and act as preliminary models intended to be updated later, or constitute re-calibrations of older models; management/planning models (55%), used to explore and identify management plans based on the response of the groundwater system to water-development or climate scenarios, sometimes under water-use constraints; and water rights models (7%), used to make water administration decisions based on model output and to quantify water shortages incurred by water users or climate changes. Results for 27 model characteristics are summarized by state and model type, and important comparisons and contrasts are highlighted. Consideration of modeling uncertainty and the management focus toward sustainability, adaptive management and resilience are discussed, and future modeling recommendations, in light of the reviewed models and other published works, are presented.

In Canada’s western Arctic, perennial discharge from permafrost watersheds is the surface manifestation of active groundwater flow systems with features including the occurrence of year-round open water and the formation of icings, yet understanding the mechanisms of groundwater recharge and flow in periglacial environments remains enigmatic. Stable isotopes (δ¹⁸O, δD, δ¹³CDIC), and noble gases have proved useful to study groundwater recharge and flow of groundwater which discharges along rivers in Canada’s western Arctic. In these studies of six catchments, groundwater recharge was determined to be a mix of snowmelt and precipitation. All systems investigated show that groundwater has recharged through organic soils with elevated PCO₂, which suggests that recharge occurs largely during summer when biological activity is high. Noble gas concentrations show that the recharge temperature was between 0 and 5 °C, which when considered in the context of discharge temperatures, suggests that there is no significant imbalance of energy flux into the subsurface. Groundwater circulation times were found to be up to 31 years for non-thermal waters using the ³ H-³He method.

Management of groundwater resources can be enhanced by using numerical models to improve development strategies. However, the lack of basic data often limits the implementation of these models. The Kharga Oasis in the western desert of Egypt is an arid area that mainly depends on groundwater from the Nubian Sandstone Aquifer System (NSAS), for which the hydrogeological data needed for groundwater simulation are lacking, thereby introducing a problem for model calibration and validation. The Grey Model (GM) was adopted to analyze groundwater flow. This model combines a finite element method (FEM) with a linear regression model to try to obtain the best-fit piezometric-level trends compared to observations. The GM simulation results clearly show that the future water table in the northeastern part of the study area will face a severe drawdown compared with that in the southwestern part and that the hydraulic head difference between these parts will reach 140 m by 2060. Given the uncertainty and limitation of available data, the GM produced more realistic results compared with those obtained from a FEM alone. The GM could be applied to other cases with similar data limitations.

The Pleistocene Kimitsu aquifer was selected for examination of the relationship between groundwater age and chemical evolution of Ca(HCO₃)₂-type groundwater. For the most part, the aquifer is confined and composed mainly of quartz and feldspar with a small amount of calcite. The groundwater ages calculated by ¹⁴C were adjusted by using a carbon mass-balance method and corrected for effects of ¹⁴C diffusion. Groundwater ages in the Kimitsu aquifer vary from modern (upgradient) to approximately 2,400 years at 4.4 km from the edge of the recharge area. The ¹⁴C age was verified by groundwater velocity calculated from the hydraulic gradient and hydraulic conductivity. The confined groundwater evolved to Ca(HCO₃)₂-type around 50 years after recharge and this has been maintained for more than 8,300 years due to low chemical reactivity, derived from equilibrium with calcite, kaolinite and Ca-montmorillonite. In addition, high pH prevents the dissolution of Fe and Mn. Consequently, the rate of increase in electrical conductivity ranges from 10 to 30 μS/cm per 1,000 years. On the other hand, leakage from the deep region, which is recognized from high Cl– levels, causes remarkable increases in CH₄ and HCO₃ – concentrations, resulting in an apparent sulfidic zone at 500-m depth in most downgradient regions.

On-going geological disposal programs for spent nuclear fuel have generated strong demands for investigation and characterization of deep-lying groundwater systems. Because of the long time scales for which radiological safety needs to be demonstrated in safety assessment applications, an analysis of the hydrogeological performance of the geosphere system during glacial climate conditions is needed. Groundwater flow at depth in crystalline rock during the passage of an ice-sheet margin is discussed based on performed groundwater-flow modeling of two bedrock sites, Forsmark and Laxemar, in the Fennoscandian Shield, Sweden. The modeled ice sheet mimics the Weichselian ice sheet during its last major advance and retreat over northern Europe. The paper elaborates and analyzes different choices of top boundary conditions at the ice sheet–subsurface interface (e.g. ice-sheet thickness and ice-margin velocity) and in the proglacial area (presence or lack of permafrost) and relates these choices to available groundwater-flow-model hydraulic output and prevailing conceptual hydrogeochemical models of the salinity evolution at the two sites. It is concluded that the choice of boundary conditions has a strong impact on results and that the studied sites behave differently for identical boundary conditions due to differences in their structural-hydraulic properties.

In cold regions, hydrologic systems possess seasonal and perennial ice-free zones (taliks) within areas of permafrost that control and are enhanced by groundwater flow. Simulation of talik development that follows lake formation in watersheds modeled after those in the Yukon Flats of interior Alaska (USA) provides insight on the coupled interaction between groundwater flow and ice distribution. The SUTRA groundwater simulator with freeze–thaw physics is used to examine the effect of climate, lake size, and lake–groundwater relations on talik formation. Considering a range of these factors, simulated times for a through-going sub-lake talik to form through 90 m of permafrost range from ∼200 to > 1,000 years (vertical thaw rates < 0.1–0.5 m yr⁻¹). Seasonal temperature cycles along lake margins impact supra-permafrost flow and late-stage cryologic processes. Warmer climate accelerates complete permafrost thaw and enhances seasonal flow within the supra-permafrost layer. Prior to open talik formation, sub-lake permafrost thaw is dominated by heat conduction. When hydraulic conditions induce upward or downward flow between the lake and sub-permafrost aquifer, thaw rates are greatly increased. The complexity of ground-ice and water-flow interplay, together with anticipated warming in the arctic, underscores the utility of coupled groundwater-energy transport models in evaluating hydrologic systems impacted by permafrost.