Intense rainstorms in 2008 resulted in wide-spread flooding across the Midwestern United States. In Wisconsin, floodwater inundated a 17.7-km²area on an outwash terrace, 7.5 m above the mapped floodplain of the Wisconsin River. Surface-water runoff initiated the flooding, but results of field investigation and modeling indicate that rapid water-table rise and groundwater inundation caused the long-lasting flood far from the riparian floodplain. Local geologic and geomorphic features of the landscape lead to spatial variability in runoff and recharge to the unconfined sand and gravel aquifer, and regional hydrogeologic conditions increased groundwater discharge from the deep bedrock aquifer to the river valley. Although reports of extreme cases of groundwater flooding are uncommon, this occurrence had significant economic and social costs. Local, state and federal officials required hydrologic analysis to support emergency management and long-term flood mitigation strategies. Rapid, sustained water-table rise and the resultant flooding of this high-permeability aquifer illustrate a significant aspect of groundwater system response to an extreme precipitation event. Comprehensive land-use planning should encompass the potential for water-table rise and groundwater flooding in a variety of hydrogeologic settings, as future changes in climate may impact recharge and the water-table elevation.

Natural mineral waters (NMW), often used to produce bottled water, are of high socio-economic interest and need appropriate management to ensure the sustainability of the resource. A complex sparkling NMW system at La Salvetat, southern France, was investigated using a multidisciplinary approach. Geological and geophysical investigations, pumping test analyses, time-series signal processing, hydrogeochemical and isotopic data (both stable and radiogenic), and numerical modelling provided complementary information on the geometry, hydrodynamic characteristics and functioning of this mineral system. The conceptual model consists of a compartmentalized reservoir characterized by two subvertical, parallel deeply rooted hydraulically independent permeable structures that are fed by deep CO₂-rich crustal fluids. The non-mineralized shallow aquifer system corresponds to a fissured layer within the weathered zone that is recharged by leakage from the overlying saprolite. This surficial aquifer responds rapidly to recharge (40–80 days), whereas the deep system’s response to recharge is much longer (up to 120 days). This research demonstrates the need for multidisciplinary approaches and modelling (quantity, hydrochemistry) for understanding complex NMW systems. This knowledge is already being applied by the bottling company that manages the resource at La Salvetat, and would be useful for conceptualising other NMW sites.

Groundwater in mountainous karst regions is vital for regional water budgets and freshwater supply. Owing to increasing water demand and climate change, detailed knowledge of the highly heterogeneous alpine aquifer systems is required. Multi-tracer analyses have been conducted in the steep karstic Wetterstein Mountains, which includes Germany’s highest summit, Zugspitze (2,962 m asl). Results of artificial tracer tests demonstrate well-developed flow paths through the unsaturated zone (up to 1,000 m thickness). Flow paths cross topographic divides and contribute to deep drainage systems underneath alpine valleys. Cross-formational flow has been identified. Quantitative analysis of tailing-dominated breakthrough curves and stable isotopes (¹⁸O) has enabled determination of the mean transit-time distribution. A fast-flow component with transit times between 3 and 13 days was found in karst conduits and open fissures, dependent on flow conditions. An intermediate-flow component, showing mean transit times of about 2.9–4.9 months, was found in well-drained fissures and fractures. A slow-flow component with mean transit times greater than 1 year is attributable to slow flow and low storage in the poorly drained fissures and rock matrix. The conceptual model enables a better understanding of drainage, water resources and vulnerability of the high-alpine karst system.

A prerequisite for minimizing contamination risk whilst conducting managed aquifer recharge (MAR) with recycled water is estimating the residence time in the zone where pathogen inactivation and biodegradation processes occur. MAR in Western Australia’s coastal aquifers is a potential major water source. As MAR with recycled water becomes increasingly considered in this region, better knowledge of applied and incidental tracer-based options from case studies is needed. Tracer data were collected at a MAR site in Floreat, Western Australia, under a controlled pumping regime over a distance of 50 m. Travel times for bromide-spiked groundwater were compared with two incidental tracers in recycled water: chloride and water temperature. The average travel time using bromide was 87 ± 6 days, whereas the estimates were longer based on water temperature (102 ± 17 days) and chloride (98 ± 60 days). The estimate of average flow velocity based on water temperature data was identical to the estimate based on bromide within a 25-m section of the aquifer (0.57 ± 0.04 m day⁻¹). This case study offers insights into the advantages, challenges and limitations of using incidental tracers in recycled water as a supplement to a controlled tracer test for estimating aquifer residence times.

Parameter structure identification is formulated in terms of solving an inverse problem, which allows for a determination of an appropriate level of parameter structure complexity, and the identification of its pattern and the associated parameter values. With the increasing complexity of parameter structure identification in groundwater modeling, demand for robust, fast, and accurate optimizers is on the rise among researchers from groundwater hydrology fields. A novel global optimizer, differential evolution (DE), has been proposed to solve the parameter-structure-identification problem. The Voronoi tessellation is adopted for the automatic parameterization. The stepwise regression method and the error covariance matrix are used to determine the optimal structure complexity. Numerical experiments with a continuous hydraulic conductivity distribution are conducted to demonstrate the proposed methodology. The results indicate that the DE can identify the global optimum effectively and efficiently. A sensitivity analysis of the control parameters and mutation schemes implemented in the DE is employed to examine their influence on the objective function. The comparison between DE and genetic algorithm shows the advantage of DE in terms of robustness and efficiency. The proposed methodology is also applied to a real groundwater system, Pingtung Plain in Taiwan, and the properties of aquifers are successfully identified.

Influences of hydraulic conductivity (K) heterogeneities on bedrock groundwater (BG) flow systems in mountainous topography are investigated using a conceptual 2D numerical modelling approach. A conceptual model for K heterogeneity in crystalline bedrock mountainous environments is developed based on a review of previous research, and represents heterogeneities due to weathering profile, bedrock fracture characteristics, and catchment-scale (∼0.1–1 km) structural features. Numerical groundwater modelling of K scenarios for hypothetical mountain catchment topography indicates that general characteristics of the BG flow directions are dominated by prominent topographic features. Within the modelled saturated BG flow system, ∼90 % or more of total BG flux is focussed within a fractured bedrock zone, extending to depths of ∼100–200 m below the ground surface, overlying lower-K bedrock. Structural features and heterogeneities, represented as discrete zones of higher or lower K relative to surrounding bedrock, locally influence BG flow, but do not influence general BG flow patterns or general positions of BG flow divides. This result is supported by similar BG transit-time distribution shapes and statistics for systems with and without structural features. The results support the development of topography-based methods for predicting general locations of BG flow-system boundaries in mountain regions.

Temporal changes in the quantity and chemical status of groundwater resources must be accurately quantified to aid sustainable management of aquifers. Monitoring data show that the groundwater level in Shahrood alluvial aquifer, northeastern Iran, continuously declined from 1993 to 2009, falling 11.4 m in 16 years. This constitutes a loss of 216 million m³from the aquifer’s stored groundwater reserve. Overexploitation and reduction in rainfall intensified the declining trend. In contrast, the reduced abstraction rate, the result of reduced borehole productivity (related to the reduction in saturated-zone thickness over time), slowed down the declining trend. Groundwater salinity varied substantially showing a minor rising trend. For the same 16-year period, increases were recorded in the order of 24% for electrical conductivity, 12.4% for major ions, and 9.9% for pH. This research shows that the groundwater-level declining trend was not interrupted by fluctuation in rainfall and it does not necessarily lead to water-quality deterioration. Water-level drop is greater near the aquifer’s recharging boundary, while greater rates of salinity rise occur around the end of groundwater flow lines. Also, fresher groundwater experiences a greater rate of salinity increase. These findings are of significance for predicting the groundwater level and salinity of exhausted aquifers.

Numerical groundwater models were used to assess groundwater sustainability on Jeju Island, South Korea, for various climate and groundwater withdrawal scenarios. Sustainability criteria included groundwater-level elevation, spring flows, and salinity. The latter was studied for the eastern sector of the island where saltwater intrusion is significant. Model results suggest that there is a need to revise the current estimate of sustainable yield of 1.77 × 10⁶m³/day. At the maximum extraction of 84 % of the sustainable yield, a 10-year drought scenario would decrease spring flows by 28 %, dry up 27 % of springs, and decrease hydraulic head by an island-wide average of 7 m. Head values are particularly sensitive to changes in recharge in the western parts of the island, due to the relatively low hydraulic conductivity of fractured volcanic aquifers and increased groundwater extraction for irrigation. Increases in salinity are highest under drought conditions around the current 2-m head contour line, with an estimated increase of up to 9 g/L under 100 % sustainable-yield use. The study lists recommendations towards improving the island’s management of potable groundwater resources. However, results should be treated with caution given the available data limitations and the simplifying assumptions of the numerical modeling approaches.

Assessing factors that influence groundwater levels such as land use and pumping strategy, is essential to adequately manage groundwater resources. A transient numerical model for groundwater flow with infiltration was developed for the Tedori River alluvial fan (140 km²), Japan. The main water input into the groundwater body in this area is irrigation water, which is significantly influenced by land use, namely paddy and upland fields. The proposed model consists of two models, a one-dimensional (1-D) unsaturated-zone water flow model (HYDRUS-1D) for estimating groundwater recharge and a 3-D groundwater flow model (MODFLOW). Numerical simulation of groundwater flow from October 1975 to November 2009 was performed to validate the model. Simulation revealed seasonal groundwater level fluctuations, affected by paddy irrigation management. However, computational accuracy was limited by the spatiotemporal data resolution of the groundwater use. Both annual groundwater levels and recharge during the irrigation periods from 1975 to 2009 showed long-term decreasing trends. With the decline in rice-planted paddy field area, groundwater recharge cumulatively decreased to 61 % of the peak in 1977. A paddy-upland crop-rotation system could decrease groundwater recharge to 73–98 % relative to no crop rotation.

Numerical investigations and a thermohydraulic evaluation are presented for two-well models of an aquifer thermal energy storage (ATES) system operating under a continuous flow regime. A three-dimensional numerical model for groundwater flow and heat transport is used to analyze the thermal energy storage in the aquifer. This study emphasizes the influence of regional groundwater flow on the heat transfer and storage of the system under various operation scenarios. For different parameters of the system, performances were compared in terms of the temperature of recovered water and the temperature field in the aquifer. The calculated temperature at the producing well varies within a certain range throughout the year, reflecting the seasonal (quarterly) temperature variation of the injected water. The pressure gradient across the system, which determines the direction and velocity of regional groundwater flow, has a substantial influence on the convective heat transport and performance of aquifer thermal storage. Injection/production rate and geometrical size of the aquifer used in the model also impact the predicted temperature distribution at each stage and the recovery water temperature. The hydrogeological-thermal simulation is shown to play an integral part in the prediction of performance of processes as complicated as those in ATES systems.