Characterizing rainfall–groundwater dynamics in a hard‐rock aquifer system using time series,geographic information system and geostatistical modelling

The aim of this study was to investigate rainfall–groundwater dynamics over space and annual time scales in a hard‐rock aquifer system of India by employing time series, geographic information system and geostatistical modelling techniques. Trends in 43‐year (1965–2007) annual rainfall time series of ten rainfall stations and 16‐year (1991–2006) pre‐monsoon and post‐monsoon groundwater levels of 140 sites were identified by using Mann–Kendall, Spearman rank order correlation and Kendall rank correlation tests. Trends were quantified by Kendall slope method. Furthermore, the study involves novelty of examining homogeneity of pre‐monsoon and post‐monsoon groundwater levels, for the first time, by applying seven tests. Regression analysis between rainfall and post‐monsoon groundwater levels was performed. The pre‐monsoon and post‐monsoon groundwater levels for four periods – 1991–1994, 1995–1998, 1999–2002 and 2003–2006 – were subjected to geographic information system‐based geostatistical modelling. The rainfall showed considerable spatiotemporal variations, with a declining trend at the Mavli rainfall station (p‐value < 0.05). The Levene's tests revealed spatial homogeneity of rainfall at α = 0.05. Regression analyses indicated significant relationships (r² > 0.5) between groundwater level and rainfall for eight rainfall stations. Non‐homogeneity and declining trends in the groundwater level, attributed to anthropogenic and hydrologic factors, were found at 5–61 more sites in pre‐monsoon compared with post‐monsoon season. The groundwater declining rates in phyllite–schist, gneiss, schist and granite formations were found to be 0.18, 0.26, 0.21 and 0.14 m year^?1 and 0.13, 0.19, 0.16 and 0.02 m year^?1 during the pre‐monsoon and post‐monsoon seasons, respectively. The geostatistical analyses for four time periods revealed linkages between the rainfall and groundwater levels. Copyright © 2013 John Wiley & Sons, Ltd.     

Dynamic evolution of nutrient discharge under stormflow and baseflow conditions in a coastal agricultural watershed in Ishigaki Island,Okinawa, Japan

Excessive terrestrial nutrient loadings adversely impact coral reefs by primarily enhancing growth of macroalgae, potentially leading to a phase‐shift phenomenon. Hydrological processes and other spatial and temporal factors affecting nutrient discharge must be examined to be able to formulate effective measures for reducing nutrient export to adjacent reefs. During storm events and baseflow periods, water samples were obtained from the tropical Todoroki River, which drains an intensively agricultural watershed into Shiraho coral reef. In situ nutrient analyzers were deployed for 6 months to hourly measure dissolved nutrient (NO₃⁻‐N and PO₄³⁻‐P) concentrations. Total phosphorus (TP) and suspended solid concentration (TSS) were increased by higher rainfall intensity (r = 0·94, p < 0·01) and river discharge Q (r = 0·88, p < 0·01). In contrast, NO₃⁻‐N concentration tends to decrease drastically (e.g. from 3 to 1 mg l⁻¹) during flood events. When base flow starts to dominate afterwards, NO₃⁻‐N manifested an increasing trend, but decreases when baseflow discharge becomes low. This counter‐clockwise hysteresis for NO₃⁻‐N highlights the significant influence of groundwater discharge. N delivery can therefore be considered a persistent process compared to sediment and P discharge, which are highly episodic in nature. Based on GIS analysis, nutrient concentration along the Todoroki River was largely affected by the percentage of sugarcane/bare areas and bedrock type. The spatial distribution of N concentration in the river reflects the considerable influence of subsurface geology—higher N levels in limestone‐dominated areas. P concentrations were directly related to the total length of artificial drainage, which enhances sediment transport. The use of high‐resolution monitoring data coupled with GIS‐based spatial analysis therefore enabled the clarification of control factors and the difference in the spatio‐temporal discharge characteristics between N and P. Thus, although erosion‐reduction schemes would reduce P discharge, other approaches (e.g. minimize fertilizer) are needed to reduce N discharge. Copyright © 2010 John Wiley & Sons, Ltd.     

Groundwater flow modelling of Kwa Ibo River watershed,southeastern Nigeria

Groundwater flow modelling of the Kwa Ibo River watershed in Abia State of Nigeria is presented in this paper with the aim of assessing the degree of interaction between the Kwa Ibo River and the groundwater regime of the thick sandy aquifer. The local geology of the area comprises the Quaternary to recent Benin Formation. Potential aquifer zones that were delineated earlier using geoelectrical resistivity soundings and borehole data for the area formed the basis for groundwater flow modelling. The watershed has been modelled with a grid of 65 rows by 43 columns and with two layers. Lateral inflow from the north has been simulated with constant heads at the Government College, Umuahia, and outflow at Usaka Elegu in the south. The Kwa Ibo River traverses the middle of the watershed from north to south. The river‐stage data at Umudike, Amawom, Ntalakwu and Usaka Elegu have been used for assigning surface water levels and riverbed elevations in the model. Permeability distribution was found to vary from 3 to 14·5 m day^?1. Natural recharge due to rainfall formed the main input to the aquifer system, and abstraction from wells was the main output. A steady‐state groundwater flow simulation was carried out and calibrated against the May 1980 water levels using 26 observation wells. The model computations have converged after 123 iterations. Under the transient‐state calibration, the highest rainfall (and hence groundwater recharge) over the 10‐year study period was recorded in 1996, whereas the lowest was recorded in 1991. The computed groundwater balance of 55 274 m³ day^?1 was comparable to that estimated from field investigations. Results from the modelling show that abstraction is much less than groundwater recharge. Hence there is the possibility for additional groundwater exploitation in the watershed through drilling of boreholes. Copyright © 2008 John Wiley & Sons, Ltd.     

Isotopic evidence for widespread cold‐season‐biased groundwater recharge and young streamflow across central Canada

Transformations of precipitation into groundwater and streamflow are fundamental hydrological processes, critical to irrigated agriculture, hydroelectric power generation, and ecosystem health. Our understanding of the timing of groundwater recharge and streamflow generation remains incomplete, limiting our ability to predict fresh water, nutrient, and contaminant fluxes, especially in large basins. Here, we analyze thousands of rain, snow, groundwater, and streamflow δ¹⁸O and δ²H values in the Nelson River basin, which covers 1.2 million km² of central Canada. We show that the fraction of precipitation that recharges aquifers is ~1.3–5 times higher for precipitation falling during cold months with subzero mean monthly temperatures than for precipitation falling during warmer months. The near‐ubiquity of cold‐season‐biased groundwater recharge implies that changes to winter water balances may have disproportionate impacts on annual groundwater recharge rates. We also show that young streamflow—defined as precipitation that enters a river in less than ~2.3 months—comprises ~27% of annual streamflow but varies widely among tributaries in the Nelson River basin (1–59%). Young streamflow fractions are lower in steep catchments and higher in flatter catchments such as the transboundary Red River basin. Our findings imply that flat, lower permeability, heavily tiled landscapes favor more rapid transmission of precipitation into rivers, possibly mobilizing excess soluble fertilizers and exacerbating eutrophication events in Lake Winnipeg.     

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.     

Spatial distribution and content of soil organic matter in an agricultural field in eastern Canada,as estimated from geostatistical tools

Soil erosion induces soil redistribution within the landscape and thus contributes to the spatial variability of soil quality. This study complements a previous experimentation initiated by the authors focusing on soil redistribution as a result of soil erosion, as indicated by caesium‐137 (¹³⁷Cs) measurements, in a small agricultural field in Canada. The spatial variability of soil organic matter (SOM) was characterized using geostatistics, which consider the randomized and structured nature of spatial variables and the spatial distribution of the samples. The spatial correlation of SOM (in percentages) patterns in the topsoil was established taking into account the spatial structure present in the data. A significant autocorrelation and reliable variograms were found with a R² ≥ 0·9, thus demonstrating a strong spatial dependence. Ordinary Kriging (OK) interpolation provided the best cross validation (r² = 0·35). OK and inverse distance weighting power two (IDW2) interpolation approaches produced similar estimates of the total SOM content of the topsoil (0–20 cm) of the experimental field, i.e. 211 and 213 tonnes, respectively. However, the two approaches produced differences in the spatial distribution patterns and the relative magnitude of some SOM content classes. The spatialization of SOM and soil redistribution variability – as evidenced by ¹³⁷Cs measurements – is a first step towards the assessment of the impact of soil erosion on SOM losses to recommend conservation measures. Copyright © 2009 John Wiley & Sons, Ltd.     

Assimilation of historical head data to estimate spatial distributions of stream bed and aquifer hydraulic conductivity fields

Management of water resources in alluvial aquifers relies mainly on understanding interactions between hydraulically connected streams and aquifers. Numerical models that simulate this interaction often are used as decision support tools for water resource management. However, the accuracy of numerical predictions relies heavily on unknown system parameters (e.g., streambed conductivity and aquifer hydraulic conductivity), which are spatially heterogeneous and difficult to measure directly. This paper employs an ensemble smoother to invert groundwater level measurements to jointly estimate spatially varying streambed and alluvial aquifer hydraulic conductivity along a 35.6‐km segment of the South Platte River in Northeastern Colorado. The accuracy of the inversion procedure is evaluated using a synthetic experiment and historical groundwater level measurements, with the latter constituting the novelty of this study in the inversion and validation of high‐resolution fields of streambed and aquifer conductivities. Results show that the estimated streambed conductivity field and aquifer conductivity field produce an acceptable agreement between observed and simulated groundwater levels and stream flow rates. The estimated parameter fields are also used to simulate the spatially varying flow exchange between the alluvial aquifer and the stream, which exhibits high spatial variability along the river reach with a maximum average monthly aquifer gain of about 2.3 m³/day and a maximum average monthly aquifer loss of 2.8 m³/day, per unit area of streambed (m²). These results demonstrate that data assimilation inversion provides a reliable and computationally affordable tool to estimate the spatial variability of streambed and aquifer conductivities at high resolution in real‐world systems.     

Spatial evolution of nitrogen and phosphorus loads along a small Mediterranean river: implication of bed sediments

This study addresses the spatial variations in water quality along the River Vène (France). The Vène drains a 67 km² rural basin, with a large karstic area, located in a Mediterranean context. A 1 day sampling campaign was conducted along the river, in winter low‐flow conditions (February 2003). Physico‐chemical parameters and water flow discharge were measured in situ during the sampling campaign. Water quality was evaluated by determining the concentrations of nitrogen and phosphorus in water and bed‐sediment samples. Nitrogen and phosphorus loads were evaluated taking into account the measured concentrations and discharge. The campaign included 18 sampling points and concerned the whole river from the spring to the outlet, plus the main inputs, i.e. sewage treatment works, main tributaries and karstic springs. The spatial evolution of nitrogen and phosphorus loads along the river allowed the significant role of point‐source inputs to be demonstrated. The decrease in nutrient loads along the river occurred mainly in specific reaches where fine sediments had accumulated. In these zones, phosphorus is trapped in the bed sediments in calcium‐bound phosphates due to precipitation processes. Copyright © 2005 John Wiley & Sons, Ltd.     

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.     

Temporal and spatial variation of water stable isotopes (18O and 2H) in the Kaidu River basin,Northwestern China

Water resources are the most critical factors to ecology and society in arid basins, such as Kaidu River basin. Isotope technique was convenient to trace this process and reveal the influence from the environment. In this paper, we try to investigate the temporal and spatial characteristics in stable isotope (¹⁸O and ²H) of surface water and groundwater in Kaidu River. Through the water stable isotope composition measurement, spatial and temporal characteristics of deuterium (δ²H) and oxygen 18 (δ¹⁸O) were analysed. It is revealed that (1) comparing the stream water line with the groundwater line and local meteorological water line of Urumqi City, it is found that the contribution of precipitation to surface water in stream runoff is the main source, whereas the surface water is the main source of groundwater. Groundwater is mainly drainage of surface runoff in the river; (2) in the main stream of Kaidu River, the spatial variability of river water showed a ‘heavier‐lighter‐heavier’ change along with the main stream for δ¹⁸O, and temporal variability showed higher in summer and lower in winter; (3) the δ¹⁸O and δ²H values of groundwater samples ranged from ?11.36 to ?7.97‰ and ?73.45 to ?60.05‰, respectively. There is an increasing trend of isotopic values along the groundwater flow path. The seasonal fluctuation of δ¹⁸O is not clear in most samples. Copyright © 2012 John Wiley & Sons, Ltd.     

Two‐dimensional flow response to tidal fluctuation in a heterogeneous aquifer‐aquitard system

In alluvial coastal aquifers, finer sediments are preferentially deposited along the downstream direction, so the hydraulic conductivity is generally heterogeneous and changes with distance from the coastline. To investigate the influence of aquifer heterogeneity on seawater‐groundwater interaction, a new two‐dimensional model characterising groundwater flow in an aquifer‐aquitard system was developed assuming that the hydraulic conductivity of the aquifer linearly increases with the distance from the coastline along the inland direction. A closed‐form analytical solution was derived using the separation‐of‐variables method. Comparing the new solution with the numerical solution by comsol Multiphysics (Sweden) based on the finite‐element method, one can see that the new solution agreed with the numerical solution very well except at the early time. We found that both aquitard leakance and the heterogeneity factor (b) could result in the propagation bias. The propagation bias represents the inconsistency between the theoretical calculation and the observed strong attenuation and small time lag between the head and tide fluctuations. The attenuation decreased with perpendicular distance from the coastline (x‐axis), whereas the time lag increased with distance along the x‐axis. The relationship between the time lag and the distance along the x‐axis seemed to be linear when b was 0.001 m^?1, whereas it obeyed a power function when b was greater than 0.01 m^?1. Copyright © 2014 John Wiley & Sons, Ltd.     

Modelling the effects of land‐use change on water and salt delivery from a catchment affected by dryland salinity in south‐east Australia

A comprehensive framework for the assessment of water and salt balance for large catchments affected by dryland salinity is applied to the Boorowa River catchment (1550 km²), located in south‐eastern Australia. The framework comprised two models, each focusing on a different aspect and operating on a different scale. A quasi‐physical semi‐distributed model CATSALT was used to estimate runoff and salt fluxes from different source areas within the catchment. The effects of land use, climate, topography, soils and geology are included. A groundwater model FLOWTUBE was used to estimate the long‐term effects of land‐use change on groundwater discharge. Unlike conventional salinity studies that focus on groundwater alone, this study makes use of a new approach to explore surface and groundwater interactions with salt stores and the stream. Land‐use change scenarios based on increased perennial pasture and tree‐cover content of the vegetation, aimed at high leakage and saline discharge areas, are investigated. Likely downstream impacts of the reduction in flow and salt export are estimated. The water balance model was able to simulate both the daily observed stream flow and salt load at the catchment outlet for high and low flow conditions satisfactorily. Mean leakage rate of about 23·2 mm year^?1 under current land use for the Boorowa catchment was estimated. The corresponding mean runoff and salt export from the catchment were 89 382 ML year^?1 and 38 938 t year^?1, respectively. Investigation of various land‐use change scenarios indicates that changing annual pastures and cropping areas to perennial pastures is not likely to result in substantial improvement of water quality in the Boorowa River. A land‐use change of about 20% tree‐cover, specifically targeting high recharge and the saline discharge areas, would be needed to decrease stream salinity by 150 µS cm^?1 from its current level. Stream salinity reductions of about 20 µS cm^?1 in the main Lachlan River downstream of the confluence of the Boorowa River is predicted. The FLOWTUBE modelling within the Boorowa River catchment indicated that discharge areas under increased recharge conditions could re‐equilibrate in around 20 years for the catchment, and around 15 years for individual hillslopes. Copyright © 2004 John Wiley & Sons, Ltd.     

Arsenic Control During Aquifer Storage Recovery Cycle Tests in the Floridan Aquifer

Implementation of aquifer storage recovery (ASR) for water resource management in Florida is impeded by arsenic mobilization. Arsenic, released by pyrite oxidation during the recharge phase, sometimes results in groundwater concentrations that exceed the 10 µg/L criterion defined in the Safe Drinking Water Act. ASR was proposed as a major storage component for the Comprehensive Everglades Restoration Plan (CERP), in which excess surface water is stored during the wet season, and then distributed during the dry season for ecosystem restoration. To evaluate ASR system performance for CERP goals, three cycle tests were conducted, with extensive water‐quality monitoring in the Upper Floridan Aquifer (UFA) at the Kissimmee River ASR (KRASR) pilot system. During each cycle test, redox evolution from sub‐oxic to sulfate‐reducing conditions occurs in the UFA storage zone, as indicated by decreasing Fe²⁺/H₂S mass ratios. Arsenic, released by pyrite oxidation during recharge, is sequestered during storage and recovery by co‐precipitation with iron sulfide. Mineral saturation indices indicate that amorphous iron oxide (a sorption surface for arsenic) is stable only during oxic and sub‐oxic conditions of the recharge phase, but iron sulfide (which co‐precipitates arsenic) is stable during the sulfate‐reducing conditions of the storage and recovery phases. Resultant arsenic concentrations in recovered water are below the 10 µg/L regulatory criterion during cycle tests 2 and 3. The arsenic sequestration process is appropriate for other ASR systems that recharge treated surface water into a sulfate‐reducing aquifer.     

Coupled surface water–groundwater model and its application in the arid Shiyang River basin,China

The planning and management of water resources in the Shiyang River basin, China require a tool for assessing the impact of groundwater and stream use on water supply reliabilities and improving many environment‐related problems such as soil desertification induced by recent water‐related human activities. A coupled model, integrating rule‐based lumped surface water model and distributed three‐dimensional groundwater flow model, has been established to investigate surface water and groundwater management scenarios that may be designed to restore the deteriorated ecological environment of the downstream portion of the Shiyang River basin. More than 66% of the water level among 24 observation wells have simulation error less than 1·0 m. The overall trend of the temporal changes of simulated and observed surface runoff at the Caiqi gauging station remains almost the same. The calibration was considered satisfactory. Initial frameworks for water allocation, including agricultural water‐saving projects, water diversion within the basin and inter‐basin water transfer, reducing agricultural irrigation area and surface water use instead of groundwater exploitation at the downstream were figured out that would provide a rational use of water resources throughout the whole basin. Sixteen scenarios were modelled to find out the most appropriate management strategies. The results showed that in the two selected management options, the groundwater budget at the Minqin basin was about 1·4 × 10⁸ m³/a and the ecological environment would be improved significantly, but the deficit existed at the Wuwei basin and the number was about 0·8 × 10⁸ m³/a. Water demand for domestic, industry and urban green area would be met in the next 30 years, but the water shortage for meeting the demand of agricultural water use in the Shiyang River basin was about 2·2 × 10⁸ m³/a. It is suggested that more inter‐basin water transfer should be required to obtain sustainable water resource use in the Shiyang River basin. Copyright © 2009 John Wiley & Sons, Ltd.     

Multiple‐indicator study of groundwater flow and chemistry and the impacts of river and paddy water on groundwater in the alluvial fan of the Tedori River,Japan

To investigate the source, flow paths, and chemistry of rich resources of high‐quality, shallow groundwater in the alluvial fan between the Tedori and Sai rivers in central Japan, we analysed stable isotope ratios of H, O, and Sr and concentrations of major dissolved ions and trace elements in groundwater, river water, and paddy water. The ⁸⁷Sr/⁸⁶Sr ratios of the groundwater are related to near‐surface geology: groundwater in sediment from the Tedori River has high ⁸⁷Sr/⁸⁶Sr ratios (>0.711), whereas that from the Sai River in the north of the fan has low ⁸⁷Sr/⁸⁶Sr ratios (<0.711). δ²H and δ¹⁸O values and ⁸⁷Sr/⁸⁶Sr ratios indicate that groundwater in the central and southern fans is recharged by the Tedori River, whereas recharge in the north is from the Sai River. Mg²⁺, Ca²⁺, Sr²⁺, HCO₃^?, and SO₄^2? concentrations and δ²H and δ¹⁸O values in the groundwater are high in the central fan and, except for the northern area, tend to increase with distance from the Tedori River. There are linear relationships between ⁸⁷Sr/⁸⁶Sr ratio and the reciprocal concentrations of Sr²⁺, Mg²⁺, and Ca²⁺. These geochemical characteristics suggest that as groundwater recharged from the Tedori River flows towards the central fan, it mixes with waters derived from precipitation and paddy water that have become enriched in these components during downward infiltration. These results are consistent with our hydrological analysis and numerical simulation of groundwater flow, thus verifying the validity of the model we used in our simulation of groundwater flow. Copyright © 2016 John Wiley & Sons, Ltd.     

Characteristics of preferential flow and groundwater discharge to Shingobee Lake,Minnesota, USA

Small‐scale heterogeneities and large changes in hydraulic gradient over short distances can create preferential groundwater flow paths that discharge to lakes. A 170 m² grid within an area of springs and seeps along the shore of Shingobee Lake, Minnesota, was intensively instrumented to characterize groundwater‐lake interaction within underlying organic‐rich soil and sandy glacial sediments. Seepage meters in the lake and piezometer nests, installed at depths of 0·5 and 1·0 m below the ground surface and lakebed, were used to estimate groundwater flow. Statistical analysis of hydraulic conductivity estimated from slug tests indicated a range from 21 to 4·8 × 10^?3 m day^?1 and small spatial correlation. Although hydraulic gradients are overall upward and toward the lake, surface water that flows onto an area about 2 m onshore results in downward flow and localized recharge. Most flow occurred within 3 m of the shore through more permeable pathways. Seepage meter and Darcy law estimates of groundwater discharge agreed well within error limits. In the small area examined, discharge decreases irregularly with distance into the lake, indicating that sediment heterogeneity plays an important role in the distribution of groundwater discharge. Temperature gradients showed some relationship to discharge, but neither temperature profiles nor specific electrical conductance could provide a more convenient method to map groundwater–lake interaction. These results suggest that site‐specific data may be needed to evaluate local water budget and to protect the water quality and quantity of discharge‐dominated lakes. Copyright © 2002 John Wiley & Sons, Ltd.     

Assessing regional‐scale spatio‐temporal patterns of groundwater–surface water interactions using a coupled SWAT‐MODFLOW model

Interaction between groundwater and surface water in watersheds has significant impacts on water management and water rights, nutrient loading from aquifers to streams, and in‐stream flow requirements for aquatic species. Of particular importance are the spatial patterns of these interactions. This study explores the spatio‐temporal patterns of groundwater discharge to a river system in a semi‐arid region, with methods applied to the Sprague River Watershed (4100 km²) within the Upper Klamath Basin in Oregon, USA. Patterns of groundwater–surface water interaction are explored throughout the watershed during the 1970–2003 time period using a coupled SWAT‐MODFLOW model tested against streamflow, groundwater level and field‐estimated reach‐specific groundwater discharge rates. Daily time steps and coupling are used, with groundwater discharge rates calculated for each model computational point along the stream. Model results also are averaged by month and by year to determine seasonal and decadal trends in groundwater discharge rates. Results show high spatial variability in groundwater discharge, with several locations showing no groundwater/surface water interaction. Average annual groundwater discharge is 20.5 m³/s, with maximum and minimum rates occurring in September–October and March–April, respectively. Annual average rates increase by approximately 0.02 m³/s per year over the 34‐year period, negligible compared with the average annual rate, although 70% of the stream network experiences an increase in groundwater discharge rate between 1970 and 2003. Results can assist with water management, identifying potential locations of heavy nutrient mass loading from the aquifer to streams and ecological assessment and planning focused on locations of high groundwater discharge. Copyright © 2016 John Wiley & Sons, Ltd.     

Estimating groundwater recharge from water isotope (δ2H, δ18O) depth profiles in the Densu River basin,Ghana

Abstract

Accurate estimation of groundwater recharge is essential for the proper management of aquifers. A study of water isotope (δ²H, δ¹⁸O) depth profiles was carried out to estimate groundwater recharge in the Densu River basin in Ghana, at three chosen observation sites that differ in their altitude, geology, climate and vegetation. Water isotopes and water contents were analysed with depth to determine water flow in the unsaturated zone. The measured data showed isotope enrichment in the pore water near the soil surface due to evaporation. Seasonal variations in the isotope signal of the pore water were also observed to a depth of 2.75 m. Below that depth, the seasonal variation of the isotope signal was attenuated due to diffusion/dispersion and low water flow velocities. Groundwater recharge rates were determined by numerical modelling of the unsaturated water flow and water isotope transport. Different groundwater recharge rates were computed at the three observation sites and were found to vary between 94 and 182 mm/year (± max. 7%). Further, the approximate peak-shift method was applied to give information about groundwater recharge rates. Although this simple method neglects variations in flow conditions and only considers advective transport, it yielded mean groundwater recharge rates of 110–250 mm/year (± max. 30%), which were in the same order of magnitude as computed numerical modelling values. Integrating these site-specific groundwater recharge rates to the whole catchment indicates that more water is potentially renewed than consumed nowadays. With increases in population and irrigation, more clean water is required, and knowledge about groundwater recharge rates – essential for improving the groundwater management in the Densu River basin – can be easily obtained by measuring water isotope depth profiles and applying a simple peak-shift approach.

Citation Adomako, D., Maloszewski, P., Stumpp, C., Osae, S. & Akiti, T. T. (2010) Estimating groundwater recharge from water isotope (δ²H, δ¹⁸O) depth profiles in the Densu River basin, Ghana. Hydrol. Sci. J. 55(8), 1405–1416.     

Relationships between riparian evapotranspiration and groundwater depth along a semiarid irrigated river valley

Evapotranspiration (ET) from riparian vegetation can be difficult to estimate due to relatively abundant water supply, spatial vegetation heterogeneity, and interactions with anthropogenic influences such as shallower groundwater tables, increased salinity, and nonpoint source pollution induced by irrigation. In semiarid south-eastern Colorado, reliable ET estimates are scarce for the riparian corridor that borders the Arkansas River. This work investigates relationships between the riparian ecosystem along the Arkansas River and an underlying alluvial aquifer using ET estimates from remotely sensed data and modelled water table depths. Results from a calibrated, finite-difference groundwater model are used to estimate weekly water table fluctuations in the riparian ecosystem from 1999 to 2009, and estimates of ET are calculated using the Operational Simplified Surface Energy Balance (SSEBop) model with over 200 Landsat scenes covering over 30 km² of riparian ecosystem along a 70-km stretch of the river. Comparison of calculated monthly SSEBop ET to estimated alfalfa reference ET from local micrometeorological station data indicated statistically significant high linear correspondence (R² = .87). Daily calculated SSEBop ET showed statistically significant moderate linear correspondence with data from a local weighing lysimeter (R² = .59). Simulated monthly SSEBop ET values were larger in drier years compared with wetter years, and ET variability was also larger in drier years. Peak ET most commonly occurred during the month of June for all 11 years of analysis. Relationships between ET and water table depth showed that peak monthly ET was highest when groundwater depths were less than about 3 m, and ET values were significantly lower for groundwater depths greater than 3 m. Negative sample Spearman correlation highlighted riparian corridor locations where ET increased as a result of decreased groundwater depths across years with different hydroclimatic conditions. This study shows how a combination of remotely sensed riparian ET estimates and a regional groundwater model can improve our understanding of linkages between riparian consumptive use and near-river groundwater conditions influenced by irrigation return flow and different climatic drivers.     

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.