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.     

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.     

Impacts of groundwater flow on the evolution of a thermokarst lake in the permafrost–dominated region on the Qinghai–Tibet plateau

As a result of global warming induced permafrost degradation in recent decades, thermokarst lakes in the Qinghai–Tibet plateau (QTP) have been regulating local hydrological and ecological processes. Simulations with coupled moisture–heat numerical models in the Beiluhe basin (located in the hinterland of permafrost regions on the QTP) have provided insights into the interaction between groundwater flow and the freeze–thaw process. A total of 30 modified SUTRA scenarios were established to examine the effects of hydrodynamic forces, permeability, and climate on thermokarst lakes. The results indicate that the hydrodynamic condition variables regulate the permafrost degradation around the lakes. In case groundwater recharges to the lake, a low–temperature groundwater flow stimulates the expansion of the surrounding thawing regions through thermal convection. The thawing rate of the permafrost underlying the lake intensifies when groundwater is discharged from the lake. Under different permeability conditions, spatiotemporal variations in the active layer thickness significantly influence the occurrence of an open talik at the lake bottom. A warmer and wetter climate will inevitably lead to a sharp decrease in the upper limit of the surrounding permafrost, with a continual decrease in the duration of open talik events. Overall, our results underscore that comprehensive consideration of the relevant hydrologic processes is critical for improving the understanding of environmental and ecological changes in cold environments.     

Spatiotemporal estimation of snow depth using point data from snow stakes,digital terrain models,and satellite data

Snow availability in Alpine catchments plays an important role in water resources management. In this paper, we propose a method for an optimal estimation of snow depth (areal extension and thickness) in Alpine systems from point data and satellite observations by using significant explanatory variables deduced from a digital terrain model. It is intended to be a parsimonious approach that may complement physical‐based methodologies. Different techniques (multiple regression, multicriteria analysis, and kriging) are integrated to address the following issues: We identify the explanatory variables that could be helpful on the basis of a critical review of the scientific literature. We study the relationship between ground observations and explanatory variables using a systematic procedure for a complete multiple regression analysis. Multiple regression models are calibrated combining all suggested model structures and explanatory variables. We also propose an evaluation of the models (using indices to analyze the goodness of fit) and select the best approaches (models and variables) on the basis of multicriteria analysis. Estimation of the snow depth is performed with the selected regression models. The residual estimation is improved by applying kriging in cases with spatial correlation. The final estimate is obtained by combining regression and kriging results, and constraining the snow domain in accordance with satellite data. The method is illustrated using the case study of the Sierra Nevada mountain range (Southern Spain). A cross‐validation experiment has confirmed the efficiency of the proposed procedure. Finally, although it is not the scope of this work, the snow depth is used to asses a first estimation of snow water equivalent resources.     

Groundwater – the disregarded component in lake water and nutrient budgets. Part 2: effects of groundwater on nutrients

Lacustrine groundwater discharge (LGD) transports nutrients from a catchment to a lake, which may fuel eutrophication, one of the major threats to our fresh waters. Unfortunately, LGD has often been disregarded in lake nutrient studies. Most measurement techniques are based on separate determinations of volume and nutrient concentration of LGD: Loads are calculated by multiplying seepage volumes by concentrations of exfiltrating water. Typically low phosphorus (P) concentrations of pristine groundwater often are increased due to anthropogenic sources such as fertilizer, manure or sewage. Mineralization of naturally present organic matter might also increase groundwater P. Reducing redox conditions favour P transport through the aquifer to the reactive aquifer‐lake interface. In some cases, large decreases of P concentrations may occur at the interface, for example, due to increased oxygen availability, while in other cases, there is nearly no decrease in P. The high reactivity of the interface complicates quantification of groundwater‐borne P loads to the lake, making difficult clear differentiation of internal and external P loads to surface water. Anthropogenic sources of nitrogen (N) in groundwater are similar to those of phosphate. However, the environmental fate of N differs fundamentally from P because N occurs in several different redox states, each with different mobility. While nitrate behaves essentially conservatively in most oxic aquifers, ammonium's mobility is similar to that of phosphate. Nitrate may be transformed to gaseous N₂ in reducing conditions and permanently removed from the system. Biogeochemical turnover of N is common at the reactive aquifer‐lake interface. Nutrient loads from LGD were compiled from the literature. Groundwater‐borne P loads vary from 0.74 to 2900 mg PO₄‐P m^?2 year^?1; for N, these loads vary from 0.001 to 640 g m^?2 year^?1. Even small amounts of seepage can carry large nutrient loads due to often high nutrient concentrations in groundwater. Large spatial heterogeneity, uncertain areal extent of the interface and difficult accessibility make every determination of LGD a challenge. However, determinations of LGD are essential to effective lake management. Copyright © 2014 John Wiley & Sons, Ltd.     

Varying water utilization of Haloxylon ammodendron plantations in a desert‐oasis ecotone

Haloxylon ammodendron is a desert shrub used extensively in China for restoring degraded dry lands. An understanding of the water source used by H. ammodendron plantations is critical achieving sustainable vegetation restoration. We measured mortality, shoot size, and rooting depth in 5‐, 10‐, 20‐, and 40‐year‐old H. ammodendron plantations. We examined stable isotopic ratios of oxygen (δ¹⁸O) in precipitation, groundwater, and soil water in different soil layers and seasons, and in plant stem water to determine water sources at different shrub ages. We found that water acquisition patterns in H. ammodendron plantations differed with plantation age and season. Thus, the main water source for 5‐year‐old shrubs was shallow soil water. Water sources of 10‐year‐old shrubs shifted depending on the soil water conditions during the season. Although their tap roots could absorb deep soil water, the plantation main water sources were from soil water, and about 50% of water originated from shallow and mid soil. This pattern might occur because main water sources in these plantations were changeable over time. The 20‐ and 40‐year‐old shrubs acquired water mainly from permanent groundwater. We conclude that the main water source of a young H. ammodendron plantation was soil water recharged by precipitation. However, when roots reached sufficient depth, water originated mainly from the deep soil water, especially in the dry season. The deeply rooted 20‐ and 40‐year‐old shrubs have the ability to exploit a deep and reliable water source. To achieve sustainability in these plantations, we recommend a reduction in the initial density of H. ammodendron in the desert‐oasis ecotone to decelerate the consumption of shallow soil water during plantation establishment.     

Likely effects of climate change on groundwater availability in a Mediterranean region of Southeastern Spain

Groundwater resources are typically the main fresh water source in arid and semi‐arid regions. Natural recharge of aquifers is mainly based on precipitation; however, only heavy precipitation events (HPEs) are expected to produce appreciable aquifer recharge in these environments. In this work, we used daily precipitation and monthly water level time series from different locations over a Mediterranean region of Southeastern Spain to identify the critical threshold value to define HPEs that lead to appreciable aquifer recharge in this region. Wavelet and trend analyses were used to study the changes in the temporal distribution of the chosen HPEs (≥20 mm day^?1) over the observed period 1953–2012 and its projected evolution by using 18 downscaled climate projections over the projected period 2040–2099. The used precipitation time series were grouped in 10 clusters according to similarities between them assessed by using Pearson correlations. Results showed that the critical HPE threshold for the study area is 20 mm day^?1. Wavelet analysis showed that observed significant seasonal and annual peaks in global wavelet spectrum in the first sub‐period (1953–1982) are no longer significant in the second sub‐period (1983–2012) in the major part of the ten clusters. This change is because of the reduction of the mean HPEs number, which showed a negative trend over the observed period in nine clusters and was significant in five of them. However, the mean size of HPEs showed a positive trend in six clusters. A similar tendency of change is expected over the projected period. The expected reduction of the mean HPEs number is two times higher under the high climate scenario (RCP8.5) than under the moderate scenario (RCP4.5). The mean size of these events is expected to increase under the two scenarios. The groundwater availability will be affected by the reduction of HPE number which will increase the length of no aquifer recharge periods (NARP) accentuating the groundwater drought in the region. Copyright © 2016 John Wiley & Sons, Ltd.     

Integration of aquifer geology,groundwater flow and arsenic distribution in deltaic aquifers – A unifying concept

Groundwater arsenic (As) presents a public health risk of great magnitude in densely populated Asian delta regions, most acutely in the Bengal Basin (West Bengal, India and Bangladesh). Research has focused on the sources, mobilisation, and heterogeneity of groundwater As, but a consistent explanation of As distribution from local to basin scale remains elusive. We show for the Bengal Aquifer System that the numerous, discontinuous silt‐clay layers together with surface topography impose a hierarchical pattern of groundwater flow, which constrains As penetration into the aquifer and controls its redistribution towards discharge zones, where it is re‐sequestered to solid phases. This is particularly so for the discrete periods of As release to groundwater in the shallow subsurface associated with sea level high‐stand conditions of Quaternary inter‐glacial periods. We propose a hypothesis concerning groundwater flow ( S ilt‐clay layers I mpose H ierarchical groundwater flow patterns constraining A rsenic progression [SIHA]), which links consensus views on the As source and history of sedimentation in the basin to the variety of spatial and depth distributions of groundwater As reported in the literature. SIHA reconciles apparent inconsistencies between independent, in some cases contrasting, field observations. We infer that lithological and topographic controls on groundwater flow, inherent to SIHA, apply more generally to deltaic aquifers elsewhere. The analysis suggests that groundwater As may persist in the aquifers of Asian deltas over thousands of years, but in certain regions, particularly at deeper levels, As will not exceed low background concentrations unless groundwater flow systems are short‐circuited by excessive pumping.     

The role of baseflow in dissolved solids delivery to streams in the Upper Colorado River Basin

Salinity has a major effect on water users in the Colorado River Basin, estimated to cause almost $300 million per year in economic damages. The Colorado River Basin Salinity Control Program implements and manages projects to reduce salinity loads, investing millions of dollars per year in irrigation upgrades, canal projects, and other mitigation strategies. To inform and improve mitigation efforts, there is a need to better understand sources of salinity to streams and how salinity has changed over time. This study explores salinity in the baseflow fraction of streamflow, assessing whether groundwater is a significant contributor of dissolved solids to streams in the Upper Colorado River Basin (UCRB). Chemical hydrograph separation was used to estimate baseflow discharge and baseflow dissolved solids loads at stream gages (n = 69) across the UCRB. On average, it is estimated that 89% of dissolved solids loads originate from the baseflow fraction of streamflow, indicating that subsurface transport processes play a dominant role in delivering dissolved solids to streams in the UCRB. A statistical trend analysis using weighted regressions on time, discharge, and season was used to evaluate changes in baseflow dissolved solids loads in streams (n = 27) from 1986 to 2011. Decreasing trends in baseflow dissolved solids loads were observed at 63% of streams. At the three most downstream sites, Green River at Green River, UT, Colorado River at Cisco, UT, and the San Juan River near Bluff, UT, baseflow dissolved solids loads decreased by a combined 823,000 metric tons (mT), which is approximately 69% of projected basin‐scale decreases in total dissolved solids loads as a result of salinity control efforts. Decreasing trends in baseflow dissolved solids loads suggest that salinity mitigation projects, landscape changes, and/or climate are reducing dissolved solids transported to streams through the subsurface. Notably, the pace and extent of decreases in baseflow dissolved solids loads declined during the most recent decade; average decreasing loads during the 2000s (28,200 mT) were only 54% of average decreasing loads in the 1990s (51,700 mT).     

Application of Distributed Temperature Sensing for coupled mapping of sedimentation processes and spatio‐temporal variability of groundwater discharge in soft‐bedded streams

The delineation of groundwater discharge areas based on Distributed Temperature Sensing (DTS) data of the streambed can be difficult in soft‐bedded streams where sedimentation and scouring processes constantly change the position of the fibre optic cable relative to the streambed. Deposition‐induced temperature anomalies resemble the signal of groundwater discharge while scouring will cause the cable to float in the water column and measure stream water temperatures. DTS applied in a looped layout with nine fibre optic cable rows in a 70 × 5 m section of a soft‐bedded stream made it possible to detect variability in streambed temperatures between October 2011 and January 2012. Detailed monthly streambed elevation surveys were carried out to monitor the position of the fibre optic cable relative to the streambed and to quantify the effect of sedimentation processes on streambed temperatures. Based on the simultaneous interpretation of streambed temperature and elevation data, a method is proposed to delineate potential high‐groundwater discharge areas and identify deposition‐induced temperature anomalies in soft‐bedded streams. Potential high‐discharge sites were detected using as metrics the daily minimum, maximum and mean streambed temperatures as well as the daily amplitude and standard deviation of temperatures. The identified potential high‐discharge areas were mostly located near the channel banks, also showing temporal variability because of the scouring and redistribution of streambed sediments, leading to the relocation of pool‐riffle sequences. This study also shows that sediment deposits of 0.1 m thickness already resulted in an increase in daily minimum streambed temperatures and decrease in daily amplitude and standard deviation. Scouring sites showed lower daily minimum streambed temperatures and higher daily amplitude and standard deviation compared with areas without sedimentation and scouring. As a limitation of the approach, groundwater discharge occurring at depositional and scouring areas cannot be identified by the metrics applied. Copyright © 2015 John Wiley & Sons, Ltd.