Exploitation of groundwater has greatly increased since the 1970s to meet the increased water demand due to fast economic development in China. Correspondingly, the regional groundwater level has declined substantially in many areas of China. Water sources are scarce in northern and northwestern China, and the anthropogenic pollution of groundwater has worsened the situation. Groundwater containing high concentrations of geogenic arsenic, fluoride, iodine, and salinity is widely distributed across China, which has negatively affected safe supply of water for drinking and other purposes. In addition to anthropogenic contamination, the interactions between surface water and groundwater, including seawater intrusion, have caused deterioration of groundwater quality. The ecosystem and geo-environment have been severely affected by the depletion of groundwater resources. Land subsidence due to excessive groundwater withdrawal has been observed in more than 50 cities in China, with a maximum accumulated subsidence of 2–3 m. Groundwater-dependent ecosystems are being degraded due to changes in the water table or poor groundwater quality. This paper reviews these changes in China, which have occurred under the impact of rapid economic development. The effects of economic growth on groundwater systems should be monitored, understood and predicted to better protect and manage groundwater resources for the future.

The close link between land use and groundwater has long been recognised, but not widely translated into integrated policy and management practices. Common understanding is needed to facilitate cross-sector dialogue on governance. The process of land-use planning advocated by the United Nations Sustainable Development Goal (UN-SDG) 15 for 2030, coupled with the launch of an independent global land-use monitoring initiative known as Land Matrix, appear to provide windows of opportunity for hydrogeologists to make specific proposals for the inclusion of groundwater protection needs in national land-use plans and the consideration of groundwater sustainability threats from major land deals and contracts. Ignoring the groundwater dimension in land-use management can result in high long-run costs for drinking-water supply and aquatic ecosystems. Thus, coordinated governance based on a coherent set of land-use sustainability criteria, aimed at enhancing both the food and groundwater harvest, is crucial for the future.

A three-dimensional groundwater model was used to improve water resource management for a study area in north-west Switzerland, where drinking-water production is close to former landfills and industrial areas. To avoid drinking-water contamination, artificial groundwater recharge with surface water is used to create a hydraulic barrier between the contaminated sites and drinking-water extraction wells. The model was used for simulating existing and proposed water management strategies as a tool to ensure the utmost security for drinking water. A systematic evaluation of the flow direction between existing observation points using a developed three-point estimation method for a large number of scenarios was carried out. It is demonstrated that systematically applying the developed methodology helps to identify vulnerable locations which are sensitive to changing boundary conditions such as those arising from changes to artificial groundwater recharge rates. At these locations, additional investigations and protection are required. The presented integrated approach, using the groundwater flow direction between observation points, can be easily transferred to a variety of hydrological settings to systematically evaluate groundwater modelling scenarios.

Sharp-interface (or interface) flow models with Dupuit-Forchheimer approximation are widely used to assess, to first order, an aquifer’s vulnerability to seawater intrusion (SWI) and to evaluate sustainable management options for coastal groundwater resources at the screening level. Recognising that interface flow models overestimate SWI, corrections have been proposed to account for the neglected mixing and also for the outflow through a finite gap. These corrections, however, were introduced in the context of specific studies and may not be generally applicable as proposed. The interface model is revisited, placing its corrections in the context of variable-density flow (VDF) theory, by expressing them in terms of the dimensionless parameters governing VDF in schematised (aspect ratio = thickness/length) homogeneous confined coastal aquifers: the coupling parameter (α), a Péclet number (Pe), and the dispersivities ratio (rα). Interfaces are compared to the 50%-salinity lines of VDF numerical solutions and regression equations are developed for estimating the outflow gap and for correcting the length of the interface (terminating with a blunted edge); the dispersion correction, which modifies the interface curvature, is restated with a variable exponent. The corrections for dispersion and for the interface length appear to be the most effective; an outflow gap is important only at small α values (strong advection relative to vertical flow due to density differences). These concepts are applied successfully to calculate the interface position in the lowermost confined sub-unit of the Coastal Plain aquifer of Israel, as an estimate of SWI.

Many numerical methods that simulate groundwater flow, particularly the continuous Galerkin finite element method, do not produce velocity information directly. Many algorithms have been proposed to improve the accuracy of velocity fields computed from hydraulic potentials. The differences in the streamlines generated from velocity fields obtained using different algorithms are presented in this report. The superconvergence method employed by FEFLOW, a popular commercial code, and some dual-mesh methods proposed in recent years are selected for comparison. The applications to depict hydrogeologic conditions using streamlines are used, and errors in streamlines are shown to lead to notable errors in boundary conditions, the locations of material interfaces, fluxes and conductivities. Furthermore, the effects of the procedures used in these two types of methods, including velocity integration and local conservation, are analyzed. The method of interpolating velocities across edges using fluxes is shown to be able to eliminate errors associated with refraction points that are not located along material interfaces and streamline ends at no-flow boundaries. Local conservation is shown to be a crucial property of velocity fields and can result in more accurate streamline densities. A case study involving both three-dimensional and two-dimensional cross-sectional models of a coal mine in Inner Mongolia, China, are used to support the conclusions presented.

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.

Many environmental problems are rooted in human behavior. This study aimed to explore the causal effect of cultural environmental bias on ‘sustainable behavior’ among agricultural groundwater users in Fars province, Iran, according to Klockner’s comprehensive model. A survey-based research project was conducted to gathering data on the paradigm of environmental psychology. The sample included agricultural groundwater users (n = 296) who were selected at random within a structured sampling regime involving study areas that represent three (higher, medium and lower) bounds of the agricultural-groundwater-vulnerability spectrum. Results showed that the “environment as ductile (EnAD)” variable was a strong determinant of sustainable behavior as it related to groundwater use, and that EnAE had the highest causal effect on the behavior of agricultural groundwater users. The adjusted model explained 41% variance of “groundwater sustainable behavior”. Based on the results, the groundwater sustainable behaviors of agricultural groundwater users were found to be affected by personal and subjective norm variables and that they are influenced by casual effects of the “environment as ductile (EnAD)” variable. The conclusions reflect the Fars agricultural groundwater users’ attitude or worldview on groundwater as an unrecoverable resource; thus, it is necessary that scientific disciplines like hydrogeology and psycho-sociology be considered together in a comprehensive approach for every groundwater study.

Hydrogeochemical analysis and multivariate statistics were applied to identify flow patterns and major processes controlling the hydrogeochemistry of groundwater in the Jianghan Plain, which is located in central Yangtze River Basin (central China) and characterized by intensive surface-water/groundwater interaction. Although HCO₃-Ca-(Mg) type water predominated in the study area, the 457 (21 surface water and 436 groundwater) samples were effectively classified into five clusters by hierarchical cluster analysis. The hydrochemical variations among these clusters were governed by three factors from factor analysis. Major components (e.g., Ca, Mg and HCO₃) in surface water and groundwater originated from carbonate and silicate weathering (factor 1). Redox conditions (factor 2) influenced the geogenic Fe and As contamination in shallow confined groundwater. Anthropogenic activities (factor 3) primarily caused high levels of Cl and SO₄ in surface water and phreatic groundwater. Furthermore, the factor score 1 of samples in the shallow confined aquifer gradually increased along the flow paths. This study demonstrates that enhanced information on hydrochemistry in complex groundwater flow systems, by multivariate statistical methods, improves the understanding of groundwater flow and hydrogeochemical evolution due to natural and anthropogenic impacts.

Submarine groundwater discharge (SGD) is an important pathway for groundwater and associated chemicals to discharge to the sea. Groundwater levels monitored along a transect perpendicular to the shoreline are used to calculate SGD flux from the nearshore aquifer to Tolo Harbor, Hong Kong (China). The calculated SGD flux—recharge/discharge measured with Darcy’s Law methods—agrees well with estimates based on geo-tracer techniques and seepage meter in Tolo Harbor during previous studies. The estimated freshwater SGD is 1.69–2.0 m²/d at the study site and 0.3 ± 0.04 cm/d for the whole of Tolo Harbor, which is comparable to the river discharge (0.25 ± 0.07 cm/d) and precipitation (0.45 ± 0.15 cm/d). The tide-driven SGD in the intertidal zone is 13.98–17.59 m²/d at the study site and 2.42 ± 0.56 cm/d for the whole of Tolo Harbor. The SGD occurring in the subtidal zone and the bottom of Tolo Harbor is 3.12 ± 4.63 cm/d. Fresh SGD accounts for ~5% of the total SGD, while the rest (~95%) is contributed by saline SGD driven by various forces. About 96% of the tide-driven SGD in the intertidal zone occurs in the ebbing tide period because the head difference between the groundwater level and sea level is great during this period. Tide-driven SGD in the spring tide is ~1.2 times that during neap tide. The tidal fluctuation amplitude and tide-driven SGD in the intertidal zone are positively correlated to each other; thus, a spring neap variation of the tide-driven SGD is observed.

Changes in groundwater level during earthquakes have been reported worldwide. In this study, field observations of co-seismic groundwater-level changes in wells under different aquifer conditions and sampling intervals due to near-field earthquake events in Taiwan are presented. Sustained changes, usually observed immediately after earthquakes, are found in the confined aquifer. Oscillatory changes due to the dynamic strain triggered by passing earthquake waves can only be recorded by a high-frequency data logger. While co-seismic changes recover rapidly in an unconfined aquifer, they can sustain for months or longer in a confined aquifer. Three monitoring wells with long-term groundwater-level data were examined to understand the association of co-seismic changes with local hydrogeological conditions. The finite element software ABAQUS is used to simulate the pore-pressure changes induced by the displacements due to fault rupture. The calculated co-seismic change in pore pressure is related to the compressibility of the formation. The recovery rate of the change is rapid in the unconfined aquifer due to the hydrostatic condition at the water table, but slow in the confined aquifer due to the less permeable confining layer. Fracturing of the confining layer during earthquakes may enhance the dissipation of pore pressure and induce the discharge of the confined aquifer. The study results indicated that aquifer characteristics play an important role in determining groundwater-level changes during and after earthquakes.