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.

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.

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.

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.

Rapid population growth in sub-Saharan West Africa and related cropland expansion were shown in some places to have increased focused recharge through ponds, raising the water table. To estimate changes in diffuse recharge, the water content and matric potential were monitored during 2009 and 2010, and modeling was performed using the Hydrus-1D code for two field sites in southwest Niger: (1) fallow land and (2) rainfed millet cropland. Monitoring results of the upper 10 m showed increased water content and matric potential to greater depth under rainfed cropland (>2.5 m) than under fallow land (≤1.0 m). Model simulations indicate that conversion from fallow land to rainfed cropland (1) increases vadose-zone water storage and (2) should increase drainage flux (∼25 mm year⁻¹) at 10-m depth after a 30–60 year lag. Therefore, observed regional increases in groundwater storage may increasingly result from diffuse recharge, which could compensate, at least in part, groundwater withdrawal due to observed expansion in irrigated surfaces; and hence, contribute to mitigate food crises in the Sahel.

Groundwater recharge sets a constraint on aquifer water balance in the context of water management. Historical data on groundwater and other relevant hydrological processes can be used to understand the effects of climatic variability on recharge, but such data sets are rare. The climate of the Canadian prairies is characterized by large inter-annual and inter-decadal variability in precipitation, which provides opportunities to examine the response of groundwater recharge to changes in meteorological conditions. A decadal study was conducted in a small (250 km²) prairie watershed in Alberta, Canada. Relative magnitude of annual recharge, indicated by water-level rise, was significantly correlated with a combination of growing-season precipitation and snowmelt runoff, which drives depression-focussed infiltration of meltwater. Annual precipitation was greater than vapour flux at an experimental site in some years and smaller in other years. On average precipitation minus vapour flux was 10 mm y⁻¹, which was comparable to the magnitude of watershed-scale groundwater recharge estimated from creek baseflow. Average baseflow showed a distinct shift from a low value (4 mm y⁻¹) in 1982–1995 to a high value (15 mm y⁻¹) in 2003–2013, indicating the sensitivity of groundwater recharge to a decadal-scale variability of meteorological conditions.

Geochemical and isotopic characterization of groundwater and lake-water samples were combined with water and total dissolved solids balances to evaluate sources of groundwater quality deterioration in eastern Hetao Basin, Inner Mongolia, China. Groundwater quality is poor; 11 of 13 wells exceed drinking-water guidelines for at least one health-based parameter and all wells exceed aesthetic guidelines. The well water is largely derived from Yellow River irrigation water. Notably high uranium concentrations in the Yellow River, relative to world rivers, suggest groundwater uranium and other trace elements may originate in the river-derived irrigation water. Complex hydrostratigraphy and spatial variation in groundwater recharge result in spatially complex groundwater flow and geochemistry. Evapotranspiration of irrigation water causes chloride concentration increases of up to two orders of magnitude in the basin, notably in shallow groundwater around Wuliangsuhai Lake. In addition to evapotranspiration, groundwater quality is affected by mineral precipitation and dissolution, silicate weathering, and redox processes. The lake-water and TDS balances suggest that a small amount of discharge to groundwater (but associated with very high solute concentrations) contributes to groundwater salinization in this region. Increasing salinity in the groundwater and Wuliangsuhai Lake will continue to deteriorate water quality unless irrigation management practices improve.

Cost-benefit analysis can be used to identify the optimum remedial strategy for contaminated groundwater, but requires that the benefits of remedial action be quantified. When the groundwater is located in an unused part of an aquifer, it can be difficult to estimate the economic loss of resource. ‘Option values’, which represent the value placed on groundwater for possible future use, may be used. Calculation of the overall loss requires consideration of the spatial (physical) extent of an aquifer where groundwater use is constrained by the presence of contamination and economic elements. An approach is presented for estimating the spatial element, represented as the annualised abstraction loss from the groundwater resource. Two distinct components have been identified based on the potential restriction on location of a new abstraction well: (1) the direct resource loss, which relates to the contaminated zone of the aquifer within which new abstraction is no longer an option, and, (2) the indirect location loss, which relates to the loss of an opportunity to site a new abstraction well in the surrounding uncontaminated aquifer due to the presence of the plume. Methods are presented for their calculation and the approach is illustrated with a case study.

A hydrogeological study was conducted in northwestern New Brunswick, Canada, to improve the predictability of fracture-dominated groundwater flow within folded bedrock composed of fine-grained turbidites. Borehole televiewer logging and outcrop mapping, integrated with hydraulic packer tests revealed enhanced hydraulic conductivity associated with northeasterly striking bedding-plane fractures formed during folding and flexural slip. These fractures impart azimuthal anisotropy to the aquifer because of moderately dipping fold limbs. High-angle fractures form a well-developed non-stratabound network, comprising two open fracture sets striking NNE parallel to the current direction of principal stress, and WNW parallel to the direction of principal stress that dominated during the Acadian orogeny. The subset of fractures showing significant oxidation, deemed most important to the groundwater flow system, is dominated by bedding-plane and high-angle fractures striking near-parallel to the maximum principal stress direction, resulting in extensional opening and enhanced hydraulic conductivities. An equivalent porous media model, incorporating anisotropy and varying hydraulic conductivity with depth, indicates that horizontal flow dominates the aquifer with relatively minor exchange between different model layers. These findings have implications for understanding flow directions in the Black Brook Watershed and elsewhere in the Matapédia Basin where fractures formed under similar stress conditions.