Patterns of groundwater hydrochemistry in Apan-Tochac sub-basin,Mexico

Abstract

On the basis of the degree of mineralization, the groundwater of Apan-Tochac sub-basin may be considered as fresh (TDS < 500 ppm). However, chlorination is necessary to make it fit for human consumption. Major ion analyses of over 235 water samples reveal a striking relationship between hydrochemical evolution and the groundwater flow system. A high content of total dissolved solids, and low values of the Ca:Mg ratio are present in wells located on the plain (discharge zone), whereas opposite conditions are associated with wells located in higher regions (recharge zone). Statistical data analysis using the method of principal components allowed to differentiation of two hydrochemical families: (a) low mineralization corresponding to the recharge zone, and (b) high mineralization corresponding to the discharge zone. Waters of the Ca + Mg + HCO₃, and Na + Mg + HCO₃ hydrochemical fades are present and the former is dominant. The water is slightly alkaline, having slight problems of salinity during the year owing mainly to Ca²⁺HCO₃ ^? and Na⁺Cl^? salts. The hydrochemistry of the groundwater reflects the pattern of local groundwater flow for this sub-basin.     

Estimation of groundwater recharge in sandy till with two different methods using groundwater level fluctuations

Two methods for estimation of groundwater recharge, both based on groundwater level fluctuations, were applied in a moraine area in southeastern Sweden. The first method utilized a onedimensional soil water model which was tested against observed groundwater levels. The boundary conditions were defined by using standard meteorological data and submodels for precipitation, snow dynamics, interception, evapotranspiration and horizontal groundwater outflow. The second method directly transformed groundwater level fluctuations to equivalent amounts of water from a constructed recession curve and the specific yield concept. Conceptually the two methods could be characterized as inflow and response methods respectively.

A good fit between observed and simulated groundwater levels was obtained by the soil water modelling. The results were, however, shown to be rather insensitive to displacement in the water balance between evapotranspiration and groundwater outflow, giving a good fit for a simulated net groundwater recharge ranging between 134 and 197 mm. The results from the attempts to use groundwater level fluctuations directly were discouraging. Compared to the soil water simulations the results were unstable and quite different for different years. It was impossible to use a constant specific yield or even different specific yields depending on depth.

The conclusion was that the possibilities to use groundwater level data for quantitative water balance studies are limited under the studied climatical and hydrogeological conditions. The modelling effort clearly demonstrated the need for a better quantitative knowledge on soil properties if water balance information is to be deduced. The soil water model though, could be a valuable tool studying variations within and between different years as well as processes and single events.     

Principles of zoning a territory by the hazard and risks of groundwater pollution

Principles of zoning a territory by the pollution hazard to groundwater depending on the pollution sources and their position in the environment are considered. Pollution hazard classifications are proposed for groundwater, pollutants (by their chemical and hydrochemical properties), and groundwater pollution risks, depending on the degree of hazard and the type of pollutant. A procedure is proposed for mapping the hazard and risks of groundwater pollution by various pollutants on the basis of the performed zoning.     

Assessment of seasonal groundwater quality and potential saltwater intrusion: a study case in Urmia coastal aquifer (NW Iran) using the groundwater quality index (GQI) and hydrochemical facies evolution diagram (HFE-D)

Coastal aquifers are at threat of salinization in most parts of the world. This work investigated the seasonal hydrochemical evolution of coastal groundwater resources in Urmia plain, NW Iran. Two recently proposed methods have been used to comparison, recognize and understand the temporal and spatial evolution of saltwater intrusion in a coastal alluvial aquifer. The study takes into account that saltwater intrusion is a dynamic process, and that seasonal variations in the balance of the aquifer cause changes in groundwater chemistry. Pattern diagrams, which constitute the outcome of several hydrochemical processes, have traditionally been used to characterize vulnerability to sea/saltwater intrusion. However, the formats of such diagrams do not facilitate the geospatial analysis of groundwater quality, thus limiting the ability of spatio-temporal mapping and monitoring. This deficiency calls for methodologies which can translate information from some diagrams such Piper diagram into a format that can be mapped spatially. Distribution of groundwater chemistry types in Urmia plain based on modified Piper diagram using GQI_Piper(mix) and GQI_Piper(dom) indices that Mixed Ca–Mg–Cl and Ca-HCO₃ are the dominant water types in the wet and dry seasons, respectively. In this study, a groundwater quality index specific to seawater intrusion (GQI_SWI) was used to check its efficiency for the groundwater samples affected by Urmia hypersaline Lake, Iran. Analysis of the main processes, by means of the Hydrochemical Facies Evolution Diagram (HFE-Diagram), provides essential knowledge about the main hydrochemical processes. Subsequently, analysis of the spatial distribution of hydrochemical facies using heatmaps helps to identify the general state of the aquifer with respect to saltwater intrusion during different sampling periods. The HFE-D results appear to be very successful for differentiating variations through time in the salinization processes caused by saltwater intrusion into the aquifer, distinguishing the phase of saltwater intrusion from the phase of recovery, and their respective evolutions. Both GQI and HFE-D methods show that hydrochemical variations can be read in terms of the pattern of saltwater intrusion and groundwater quality status. But generally, in this case (i.e. saltwater and not seawater intrusion) the HFE-D method was presented better efficiency than GQI method (including GQI_Piper and GQI_SWI).     

Isotopic and hydrochemical insights into the groundwater characteristics along an arid to semi-humid climate gradient in China

In arid to semi-arid regions, groundwater is a critical water resource heavily relied upon, with the recharge sources and patterns being predominantly shaped by climate change and regional disparities. To compare the characteristics of groundwater in the endorheic and exorheic river basins with the climate transition zone of Gansu Province, this study uses isotopic hydrochemical analyses. This study summarizes the differences in regional groundwater recharge and evolutionary patterns. The results shows that the distribution patterns of precipitation isotopes in endorheic and exorheic river basins are opposite to those of groundwater isotopes. Specifically, the precipitation in the endorheic areas is more depleted in heavy isotopes, whereas the groundwater is more enriched. Both endorheic areas and exorheic areas exhibit similar characteristics of groundwater hydrochemical evolution, evolving from low-mineralization Mg²⁺ HC ${\mathrm{O}}_3^{-}$ recharge water to Na⁺ Cl⁻ type water with saline characteristics. The former is primarily replenished by surface water, whereas the latter is primarily replenished by precipitation. Variations in recharge patterns along with the differences in climatic conditions lead to distinct groundwater conditions in the two regions.     

Environmental isotopes and hydrochemical study applied to surface water and groundwater interaction in the Awash River basin

An environmental isotope and hydrochemical study was carried out to conceptualize the surface water and groundwater interaction and to explore the groundwater flow pattern in relation to the geological setting. More emphasis is given to the Afar Depression where groundwater is a vital source of water supply. Conventional field hydrogeological study and river discharge records support the isotope and hydrochemical analysis. The region is tectonically active, comprising rift volcanic terrain bordered by highlands. The result revealed that recent meteoric water is the major source of recharge. Three distinct groundwater zones were identified associated with the highlands, transitional escarpment and the rift. Towards the rift, the ionic concentration and isotopic enrichment (δ²H and δ^18degO) increases following the groundwater flow paths, which is strongly controlled by axial rift faults. The groundwater flow converges to the seismically active volcano–tectonic depressions with internal drainage and to the Awash River. Within the Afar Depression, at least four groundwater regimen are identified: (1) fresh and shallow groundwater associated with alluvial deposits ultimately recharged by isotopically depleted recent highland rainfall and the evaporated Awash River; (2) cold and relatively younger groundwater within localized fractured volcanics showing mixed origin in axial fault zones; (3) old groundwater with very high ionic concentration and low isotopic signature localized in deep volcanic aquifers; and (4) old and hot saline groundwaters connected to geothermal systems. The study demonstrated that dependable groundwater can only be obtained from the first two aquifer types in aerially restricted zones in flat plains following river courses, local wadis and volcano–tectonic depressions. The conventional hydrogeological survey and discharge records indicate substantial channel losses from the Awash River, which becomes a more dominant source of recharge in central and lower Awash valleys. Copyright © 2007 John Wiley & Sons, Ltd.     

Using hydrochemical,stable isotope,and river water recharge data to identify groundwater flow paths in a deeply buried karst system

The deeply buried river‐connected Xishan karst aquifer (XKA) in western Beijing, China, has been suffering from diminishing recharge for several decades, which in turn leads to the disappearing of spring water outflows and continuously lowering of groundwater level in the area. Thus, it is important to correctly recognize the groundwater recharge and flow paths for the sustainable development of the XKA. To investigate these issues, the hydrochemical and isotopic compositions are analysed for both surface water and groundwater samples collected over an area of about 280 km². Results show that (a) the river water is characterized by high Na contents; (b) the δ²H and δ¹⁸O values in the river water are distinctively higher than those of groundwater samples, after experiencing the long‐time evaporative enrichment in the upstream reservoir; (c) the Sr concentrations and ⁸⁷Sr/⁸⁶Sr ratios of groundwater clearly indicated the interaction between water and carbonate minerals but excluded the water–silicate interaction; and (d) the groundwater samples in the direct recharge area of the XKA have the lowest Na concentrations and the δ²H and δ¹⁸O values. Based on the large differences in the Na contents and ¹⁸O values of groundwater and surface water, a simple two‐component mixing model is developed for the study area and the fractions of the river water are estimated for groundwater samples. We find that the distribution pattern of the river water fractions in the XKA clearly shows a change of directions in the preferential flow path of the groundwater from its source zone to the discharge area. Overall, our results suggest that the recharged surface water can be a useful evidence for delineating the groundwater flow path in river‐connected karst aquifer. This study improves our understanding of the heterogeneity in karst groundwater systems.     

Dispersion of naturally occurring ions in groundwater from various rock types in a portion of the San Pedro river basin, Arizona

The groundwater in an alluvial basin in southern Arizona was analyzed for concentrations of Ca⁺⁺, Mg⁺⁺, Na⁺, and Cl⁻, ions.

The variety of rock types in the area, plus the undisturbed state of the groundwater basin, make comparative rock mineralization-groundwater ionization interpretations possible. Ionic dispersion in groundwater eminating from source areas composed of differing rock types is plotted as isogram maps. These isolated areas of differing mineral composition each exhibit a unique ionic contribution to groundwater. The ion concentrations in groundwater were then used as naturally occurring tracers to determine source areas of recharge and to delineate subsurface barriers to the normal basin flow net. Ion dispersion plots reveal the carbonates of the Dragoon Mountains to be a major contributor of Ca⁺⁺ and Mg⁺⁺ to the deep alluvial portion of the basin. Cl⁻ dispersion patterns show the granitic intrusives of the Tombstone Hills produce a barrier effect in the normal flow pattern of the basin as well as being a contributor of Cl⁻ to groundwater.     

Modeling surface water–groundwater interaction in New Zealand: Model development and application

Most rivers worldwide have a strong interaction with groundwater when they leave the mountains and flow over alluvial plains before flowing into the seas or disappearing in the deserts, and in New Zealand, typically, rivers lose water to the groundwater in the upper plains and generally gain water from the groundwater in the lower plains. Aiming at simulating surface water–groundwater interaction nationally in New Zealand, we developed a conceptual groundwater module for the national hydrologic model TopNet to simulate surface water–groundwater interaction, groundwater flow, and intercatchment groundwater flow. The developed model was applied to the Pareora catchment in South Island of New Zealand, where there are concurrent spot gauged flows. Results show that the model simulations not only fit quite well to flow measurement but also to concurrent spot gauged flows, and compared to the original TopNet, it has a significant improvement in the low flows. Sensitivity analysis shows river flow is sensitive to the river losing/gaining rate instead of groundwater characteristic, while groundwater storage is sensitive to both river losing/gaining rate and groundwater characteristic. This indicates our conceptual approach is promising for nationwide modeling without the large amount of geology and aquifer data typically required by physically‐based modeling approaches.     

Proglacial groundwater storage dynamics under climate change and glacier retreat

Proglacial aquifers are an important water store in glacierised mountain catchments that supplement meltwater-fed river flows and support freshwater ecosystems. Climate change and glacier retreat will perturb water storage in these aquifers, yet the climate-glacier-groundwater response cascade has rarely been studied and remains poorly understood. This study implements an integrated modelling approach that combines distributed glacio-hydrological and groundwater models with climate change projections to evaluate the evolution of groundwater storage dynamics and surface-groundwater exchanges in a temperate, glacierised catchment in Iceland. Focused infiltration along the meltwater-fed Virkisá River channel is found to be an important source of groundwater recharge and is projected to provide 14%–20% of total groundwater recharge by the 2080s. The simulations highlight a mechanism by which glacier retreat could inhibit river recharge in the future due to the loss of diurnal melt cycling in the runoff hydrograph. However, the evolution of proglacial groundwater level dynamics show considerable resilience to changes in river recharge and, instead, are driven by changes in the magnitude and seasonal timing of diffuse recharge from year-round rainfall. The majority of scenarios simulate an overall reduction in groundwater levels with a maximum 30-day average groundwater level reduction of 1 m. The simulations replicate observational studies of baseflow to the river, where up to 15% of the 30-day average river flow comes from groundwater outside of the melt season. This is forecast to reduce to 3%–8% by the 2080s due to increased contributions from rainfall and meltwater runoff. During the melt season, groundwater will continue to contribute 1%–3% of river flow despite significant reductions in meltwater runoff inputs. Therefore it is concluded that, in the proglacial region, groundwater will continue to provide only limited buffering of river flows as the glacier retreats.     

The impact of well drawdowns on the mixing process of river water and groundwater and water quality in a riverside well field,Northeast China

Riverbank filtration (RBF) has been widely used throughout the world as an effective means to regulate surface water and groundwater resources and pretreat raw water for municipal water supply. The quality of the water from a riverside well field and the mixing ratios of surface water and groundwater is primarily impacted by the hydrodynamic processes in the RBF system. The RBF system is largely controlled by the water exploiting system in addition to the natural hydrologic condition of the river–aquifer system. As one of the most important design parameters of the riverside well field, the drawdown of groundwater level greatly determines the water head differences between the river water and groundwater as well as the field flow net, which subsequently impacts the mixing of river water and groundwater and water quality significantly. This study aimed to improve the understanding of the mixing process between the surface water and groundwater and estimate the impact of the RBF on the mixing ratio of surface water–groundwater and water quality quantitatively. A set of field pumping tests with various groundwater level drawdowns were carried out independently and successively at a riverside field with a single pumping well near the Songhua River in Northeast China in August 2017. During these tests, the water levels and hydrochemical parameters of the Songhua River, the adjacent aquifer, and the pumping well were monitored. The dynamic mixing process of the river water and groundwater and water quality under various drawdown conditions were analysed systematically using analytical methods. The results obtained from Dupuit method and the Mirror Image method in conjunction with the Hydrochemical Tracing method showed that the pumping water directly from the river water reached 60% ± 10% after a steady flow net was established. The larger the proportion of the pumping water captured from the river, the better quality of the pumping water was, because the quality of the river water (only limited to some water quality parameters monitored which were minority) was better than that of the groundwater. The results also showed that total Fe, TDS, total hardness, COD_Mn, and K⁺ were relatively sensitive to the changes of groundwater drawdown, and their concentrations decreased with an increase in the groundwater drawdown. It can be concluded that both the mixing ratio of the surface water and the groundwater and the water quality of the riverside well field can be regulated through adjusting the designed drawdown of the groundwater level, which is helpful for the design and the optimization of the riverside well water intake engineering.     

Stable isotopic and geochemical identification of groundwater evolution and recharge sources in the arid Shule River Basin of Northwestern China

Stable isotopic (δD_VSMOW and δ¹⁸O_VSMOW) and geochemical signatures were employed to constrain the geochemical evolution and sources of groundwater recharge in the arid Shule River Basin, Northwestern China, where extensive groundwater extraction occurs for agricultural and domestic supply. Springs in the mountain front of the Qilian Mountains, the Yumen‐Tashi groundwater (YTG), and the Guazhou groundwater (GZG) were Ca‐HCO₃, Ca‐Mg‐HCO₃‐SO₄ and Na‐Mg‐SO₄‐Cl type waters, respectively. Total dissolved solids (TDS) and major ion (Mg²⁺, Na⁺, Ca²⁺, K⁺, SO₄^2?, Cl^? and NO₃^?) concentrations of groundwater gradually increase from the mountain front to the lower reaches of the Guazhou Basin. Geochemical evolution in groundwater was possibly due to a combination of mineral dissolution, mixing processes and evapotranspiration along groundwater flow paths. The isotopic and geochemical variations in melt water, springs, river water, YTG and GZG, together with the end‐member mixing analysis (EMMA) indicate that the springs in the mountain front mainly originate from precipitation, the infiltration of melt water and river in the upper reaches; the lateral groundwater from the mountain front and river water in the middle reaches are probably effective recharge sources for the YTG, while contribution of precipitation to YTG is extremely limited; the GZG is mainly recharged by lateral groundwater flow from the Yumen‐Tashi Basin and irrigation return flow. The general characteristics of groundwater in the Shule River Basin have been initially identified, and the results should facilitate integrated management of groundwater and surface water resources in the study area. Copyright © 2015 John Wiley & Sons, Ltd.     

Multi isotopic approach to study temporal variation of groundwater quality in coastal aquifer of Saijo plain, Shikoku Island, Japan

In order to delineate geochemical processes and its relation with groundwater quality evolution, hydrochemical and isotopic analysis were carried out in coastal groundwater of Saijo plain, western Japan. From analytical results, even within a small distance from the coast; ionic concentration of water samples varies a lot which infers non-homogenous and patchy distribution of different aquifer system (complex geology). From stable isotopic results, it was found that most of sample points plotted near the local meteoric water line (LMWL) i.e. origin of ground water is meteoric in principle; however point away from the LMWL might favors exchange with rock minerals.     

Multi isotopic approach to study temporal variation of groundwater quality in coastal aquifer of Saijo plain, Shikoku Island, Japan

In order to delineate geochemical processes and its relation with groundwater quality evolution, hydrochemical and isotopic analysis were carried out in coastal groundwater of Saijo plain, western Japan. From analytical results, even within a small distance from the coast; ionic concentration of water samples varies a lot which infers non-homogenous and patchy distribution of different aquifer system (complex geology). From stable isotopic results, it was found that most of sample points plotted near the local meteoric water line (LMWL) i.e. origin of ground water is meteoric in principle; however point away from the LMWL might favors exchange with rock minerals.     

Evaluation of the soil water balance in an alluvial flood plain with a shallow groundwater table

Abstract

Many of the hydrological and ecological functions of alluvial flood plains within watersheds depend on the water flow exchanges between the vadoze soil zone and the shallow groundwater. The water balance of the soil in the flood plain is investigated, in order to evaluate the main hydrological processes that underlie the temporal dynamics of soil moisture and groundwater levels. The soil moisture and the groundwater level in the flood plain were monitored continuously for a three-year period. These data were integrated with the results derived from applying a physically-based numerical model which simulated the variably-saturated vertical water flow in the soil. The analysis indicated that the simultaneous processes of lateral groundwater flow and the vertical recharge from the unsaturated zone caused the observed water table fluctuations. The importance of these flows in determining the rises in the water table varied, depending on soil moisture and groundwater depth before precipitation. The monitoring period included two hydrological years (September 2009–September 2011). About 13% of the precipitation vertically recharged the groundwater in the first year and about 50% in the second. The difference in the two recharge coefficients was in part due to the lower groundwater levels in the recharge season of the first hydrological year, compared to those observed in the second. In the latter year, the shallow groundwater increased the soil moisture in the unsaturated zone due to capillary rise, and so the mean hydraulic conductivity of the unsaturated soil was high. This moisture state of soil favoured a more efficient conversion of infiltrated precipitation into vertical groundwater recharge. The results show that groundwater dynamics in the flood plain are an important source of temporal variability in soil moisture and vertical recharge processes, and this variability must be properly taken into account when the water balance is investigated in shallow groundwater environments.

Citation Pirastru, M. and Niedda, M., 2013. Evaluation of the soil water balance in an alluvial flood plain with a shallow groundwater table. Hydrological Sciences Journal, 58 (4), 898–911.     

Simulated groundwater interaction with rivers and springs in the Heihe river basin

Surface water and groundwater in the Heihe river basin of China are interconnected and the pattern of water resources exploitation has a direct effect on the interaction of groundwater and surface water, especially on a downstream oasis. A three‐dimensional groundwater flow simulation model with eight model layers was established to simulate the regional groundwater flow in the multilayered aquifer system and the interaction among the rivers, springs, and groundwater. The model was calibrated not only with historical water levels but also with the investigated baseflow and spring flux. The simulation results of the numerical model match reasonably well with the observed groundwater levels, baseflow to rivers, and spring flux. The numerical simulation also demonstrates that the hydraulic connection between the river and the aquifers has transferred from the coupling to decoupling at some reaches. It is suggested that there is a vital need to reduce groundwater withdrawal and to rationalize the use of both groundwater and surface water in order to maintain sustainable development in the study area. Copyright © 2007 John Wiley & Sons, Ltd.     

Comparison of environmental tracers including organic micropollutants as groundwater exfiltration indicators into a small river of a karstic catchment

Understanding groundwater–surface water (GW–SW) interactions is vital for water management in karstic catchments due to its impact on water quality. The objective of this study was to evaluate and compare the applicability of seven environmental tracers to quantify and localize groundwater exfiltration into a small, human-impacted karstic river system. Tracers were selected based on their emission source to the surface water either as (a) dissolved, predominantly geogenic compounds (radon-222, sulphate and electrical conductivity) or (b) anthropogenic compounds (predominantly) originating from wastewater treatment plant (WWTP) effluents (carbamazepine, tramadol, sodium, chloride). Two contrasting sampling approaches were compared (a) assuming steady-state flow conditions and (b) considering the travel time of the water parcels (Lagrangian sampling) through the catchment to account for diurnal changes in inflow from the WWTP. Spatial variability of the concentrations of all tracers indicated sections of preferential groundwater inflow. Lagrangian sampling techniques seem highly relevant for capturing dynamic concentration patterns of WWTP-derived compounds. Quantification of GW inflow with the finite element model FINIFLUX, based on observed in-stream Rn activities led to plausible fluxes along the investigated river reaches (0.265 m³ s⁻¹), while observations of other natural or anthropogenic environmental tracers produced less plausible water fluxes. Important point sources of groundwater exfiltration can be ascribed to locations where the river crosses geological fault lines. This indicates that commonly applied concepts describing groundwater–surface water interactions assuming diffuse flow in porous media are difficult to transfer to karstic river systems whereas concepts from fractured aquifers may be more applicable. In general, this study helps selecting the best suited hydrological tracer for GW exfiltration and leads to a better understanding of processes controlling groundwater inflow into karstic river systems.     

Improved simulation of river water and groundwater exchange in an alluvial plain using the SWAT model

Hydrological interaction between surface and subsurface water systems has a significant impact on water quality, ecosystems and biogeochemistry cycling of both systems. Distributed models have been developed to simulate this function, but they require detailed spatial inputs and extensive computation time. The soil and water assessment tool (SWAT) model is a semi‐distributed model that has been successfully applied around the world. However, it has not been able to simulate the two‐way exchanges between surface water and groundwater. In this study, the SWAT‐landscape unit (LU) model – based on a catena method that routes flow across three LUs (the divide, the hillslope and the valley) – was modified and applied in the floodplain of the Garonne River. The modified model was called SWAT‐LUD. Darcy's equation was applied to simulate groundwater flow. The algorithm for surface water‐level simulation during flooding periods was modified, and the influence of flooding on groundwater levels was added to the model. Chloride was chosen as a conservative tracer to test simulated water exchanges. The simulated water exchange quantity from SWAT‐LUD was compared with the output of a two‐dimensional distributed model, surface–subsurface water exchange model. The results showed that simulated groundwater levels in the LU adjoining the river matched the observed data very well. Additionally, SWAT‐LUD model was able to reflect the actual water exchange between the river and the aquifer. It showed that river water discharge has a significant influence on the surface–groundwater exchanges. The main water flow direction in the river/groundwater interface was from groundwater to river; water that flowed in this direction accounted for 65% of the total exchanged water volume. The water mixing occurs mainly during high hydraulic periods. Flooded water was important for the surface–subsurface water exchange process; it accounted for 69% of total water that flowed from the river to the aquifer. The new module also provides the option of simulating pollution transfer occurring at the river/groundwater interface at the catchment scale. Copyright © 2015 John Wiley & Sons, Ltd.     

River–aquifer interactions in a semi‐arid environment stressed by groundwater abstraction

Rivers and aquifers are, in many cases, a connected resource and as such the interactions between them need to be understood and quantified for the resource to be managed appropriately. The objective of this paper is to advance the understanding of river–aquifer interactions processes in semi‐arid environments stressed by groundwater abstraction. This is performed using data from a specific catchment where records of precipitation, evapotranspiration, river flow, groundwater levels and groundwater abstraction are analysed using basic statistics, hydrograph analysis and a simple mathematical model to determine the processes causing the spatial and temporal changes in river–aquifer interactions. This combined approach provides a novel but simple methodology to analyse river–aquifer interactions, which can be applied to catchments worldwide. The analysis revealed that the groundwater levels have declined (~ 3 m) since the onset of groundwater abstraction. The decline is predominantly due to the abstraction rather than climatic changes (r = 0.84 for the relationship between groundwater abstraction and groundwater levels; r = 0.92 for the relationship between decline in groundwater levels and magnitude of seasonal drawdown). It is then demonstrated that, since the onset of abstraction, the river has changed from being gaining to losing during low‐flow periods, defined as periods with flow less than 0.5, 1.0 or 1.5 GL/day (1 GL/day = 1 × 10⁶ m³/day). If defined as < 1.0 GL/day, low‐flow periods constitute approximately 65% of the river flows; the periods where the river is losing at low‐flow conditions are thus significant. Importantly, there was a significant delay (> 10 years) between the onset of groundwater abstraction and the changeover from gaining to losing conditions. Finally, a relationship between the groundwater gradient towards the river and the river flow at low‐flow is demonstrated. The results have important implications for water management as well as water ecology and quality. Copyright © 2012 John Wiley & Sons, Ltd.     

Analytical investigation of the exact groundwater divide between rivers beyond the Dupuit–Forchheimer approximation

The groundwater divide is a key feature of river basins and significantly influenced by subsurface hydrological processes. For an unconfined aquifer between two parallel rivers or ditches, it has long been defined as the top of the water table based on the Dupuit–Forchheimer approximation. However, the exact groundwater divide is subject to the interface between two local flow systems transporting groundwater to rivers from the infiltration recharge. This study contributes a new analytical model for two-dimensional groundwater flow between rivers of different water levels. The flownet is delineated in the model to identify groundwater flow systems and the exact groundwater divide. Formulas with two dimensionless parameters are derived to determine the distributed hydraulic head, the top of the water table and the groundwater divide. The locations of the groundwater divide and the top of the water table are not the same. The distance between them in horizontal can reach up to 8.9% of the distance between rivers. Numerical verifications indicate that simplifications in the analytical model do not significantly cause misestimates in the location of the groundwater divide. In contrast, the Dupuit–Forchheimer approximation yields an incorrect water table shape. The new analytical model is applied to investigate groundwater divides in the Loess Plateau, China, with a Monte Carlo simulation process taking into account the uncertainties in the parameters.