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.     

A long‐term study of stable isotopes as tracers of processes governing water flow and quality in a lowland river basin: the upper Thames,UK

A long‐term study of O, H and C stable isotopes has been undertaken on river waters across the 7000‐km² upper Thames lowland river basin in the southern UK. During the period, flow conditions ranged from drought to flood. A 10‐year monthly record (2003–2012) of the main River Thames showed a maximum variation of 3‰ (δ¹⁸O) and 20‰ (δ²H), although interannual average values varied little around a mean of –6.5‰ (δ¹⁸O) and –44‰ (δ²H). A δ²H/δ¹⁸O slope of 5.3 suggested a degree of evaporative enrichment, consistent with derivation from local rainfall with a weighted mean of –7.2‰ (δ¹⁸O) and –48‰ (δ²H) for the period. A tendency towards isotopic depletion of the river with increasing flow rate was noted, but at very high flows (>100 m³/s), a reversion to the mean was interpreted as the displacement of bank storage by rising groundwater levels (corroborated by measurements of specific electrical conductivity). A shorter quarterly study (October 2011–April 2013) of isotope variations in 15 tributaries with varying geology revealed different responses to evaporation, with a well‐correlated inverse relationship between Δ¹⁸O and baseflow index for most of the rivers. A comparison with aquifer waters in the basin showed that even at low flow, rivers rarely consist solely of isotopically unmodified groundwater. Long‐term monitoring (2003–2007) of carbon stable isotopes in dissolved inorganic carbon (DIC) in the Thames revealed a complex interplay between respiration, photosynthesis and evasion, but with a mean interannual δ¹³C‐DIC value of –14.8 ± 0.5‰, exchange with atmospheric carbon could be ruled out. Quarterly monitoring of the tributaries (October 2011–April 2013) indicated that in addition to the aforementioned factors, river flow variations and catchment characteristics were likely to affect δ¹³C‐DIC. Comparison with basin groundwaters of different alkalinity and δ¹³C‐DIC values showed that the origin of river baseflow is usually obscured. The findings show that long‐term monitoring of environmental tracers can help to improve the understanding of how lowland river catchments function. Copyright © NERC 2015. Hydrological Processes © 2015 John Wiley & Sons, Ltd.     

Using hysteretic behaviour and hydrograph classification to identify hydrological function across the “hillslope–depression–stream” continuum in a karst catchment

In cockpit karst landscapes, fluxes from upland areas contribute large volumes of water to low-lying depressions and stream flow. Hydrograph hysteresis and similarity between monitoring sites is important for understanding the space–time variability of hydrologic responses across the “hillslope–depression–stream” continuum. In this study, the hysteretic feature of hydrographs was assessed by characterizing the loop-like relationships between responses at upstream sites relative to subsurface discharge at the outlet of a small karst catchment. A classification of hydrograph responses based on the multi-scale smoothing Kernel -derived distance classifies the hydrograph responses on the basis of similarities between hillslope and depression sites, and those at the catchment outlet. Results demonstrate that the temporal and spatial variability of hydrograph hysteresis and similarity between hillslope flow and outlet stream flow can be explained by the local heterogeneity of depression aquifer. Large depression storage deficits emerging in the highly heterogeneous aquifer produce strong hysteresis and multiple relationships of upstream hydrographs relative to the outlet subsurface discharge. In contrast, when depression storage deficits are filled during consecutive rainfall events, depression hydrographs at the high permeability sites are almost synchronous or exhibit a monotonous function with the hydrographs at the outlet. This reduced hydrograph hysteresis enhances preferential flow paths in fractured rocks and conduits that can accelerate the hillslope flow to the outlet. Therefore, classification of hydrograph similarities between any upstream sites and the catchment outlet can help to identify the dominant hydrological functions in the heterogeneous karst catchment.     

Investigation of the elevation of saltwater wedge due to subsurface dams

Subsurface dams are rather effective and used for the prevention of saltwater intrusion in coastal regions around the world. We carried out the laboratory experiments to investigate the elevation of saltwater wedge after the construction of subsurface dams. The elevation of saltwater wedge refers to the upward movement of the downstream saltwater wedge because the subsurface dams obstruct the regional groundwater flow and reduce the freshwater discharge. Consequently, the saltwater wedge cannot further extend in the longitudinal direction but rises in the vertical profile resulting in significant downstream aquifer salinization. In order to quantitatively address this issue, field-scale numerical simulations were conducted to explore the influence of various dam heights, distances, and hydraulic gradients on the elevation of saltwater wedge. Our investigation shows that the upward movement of the saltwater wedge and its areal extension in the vertical domain of the downstream aquifer become more severe with a higher dam and performed a great dependence on the freshwater discharge. Furthermore, the increase of the hydraulic gradient and the dam distance from the sea boundary leads to a more pronounced wedge elevation. This phenomenon comes from the variation of the freshwater discharge due to the modification of dam height, location, and hydraulic gradient. Large freshwater discharge can generate greater repulsive force to restrain the elevation of saltwater wedge. These conclusions provide theoretical references for the behaviour of the freshwater–seawater interface after the construction of subsurface dams and help optimize the design strategy to better utilize the coastal groundwater resources.     

A nonstationary extreme value distribution for analysing the cessation of karst spring discharge

The effects of climate change and population growth in recent decades are leading us to consider their combined and potentially extreme consequences, particularly regarding hydrological processes, which can be modeled using a generalized extreme value (GEV) distribution. Most of the GEV models were based on a stationary assumption for hydrological processes, in contrast to the nonstationary reality due to climate change and human activities. In this paper, we present the nonstationary generalized extreme value (NSGEV) distribution and use it to investigate the risk of Niangziguan Springs discharge decreasing to zero. Rather than assuming the location, scale, and shape parameters to be constant as one might do for a stationary GEV distribution analysis, the NSGEV approach can reflect the dynamic processes by defining the GEV parameters as functions of time. Because most of the GEV model is designed to evaluate maxima (e.g. flooding, represented by positive numbers), and spring discharge cessation is a ?minima’, we deduced an NSGEV model for minima by applying opposite numbers, i.e. negative instead of positive numbers. The results of the model application to Niangziguan Springs showed that the probability of zero discharge at Niangziguan Springs will be 1/80 in 2025, and 1/10 in 2030. After 2025, the rate of decrease in spring discharge will accelerate, and the probability that Niangziguan Springs will cease flowing will dramatically increase. The NSGEV model is a robust method for analysing karst spring discharge. Copyright © 2013 John Wiley & Sons, Ltd.     

Divergence and flow direction as indicators of subsurface heterogeneity and stage‐dependent storage in alluvial floodplains

Assuming homogeneity in alluvial aquifers is convenient, but limits our ability to accurately predict stream‐aquifer interactions. Research is needed on (i) identifying the presence of focused, as opposed to diffuse, groundwater discharge/recharge to streams and (ii) the magnitude and role of large‐scale bank and transient storage in alluvial floodplains relative to changes in stream stage. The objective of this research was to document and quantify the effect of stage‐dependent aquifer heterogeneity and bank storage relative to changes in stream stage using groundwater flow divergence and direction. Monitoring was performed in alluvial floodplains adjacent to the Barren Fork Creek and Honey Creek in northeastern Oklahoma. Based on results from subsurface electrical resistivity mapping, observation wells were installed in high and low electrical resistivity subsoils. Water levels in the wells were recorded real time using pressure transducers (August to October 2009). Divergence was used to quantify heterogeneity (i.e. variation in hydraulic conductivity, porosity, and/or aquifer thickness), and flow direction was used to assess the potential for large‐scale (100 m) bank or transient storage. Areas of localized heterogeneity appeared to act as divergence zones allowing stream water to quickly enter the groundwater system, or as flow convergence zones draining a large groundwater area. Maximum divergence or convergence occurred with maximum rates of change in flow rates or stream stage. Flow directions in the groundwater changed considerably between base and high flows, suggesting that the floodplains acted as large‐scale bank storage zones, rapidly storing and releasing water during passage of a storm hydrograph. During storm events at both sites, the average groundwater direction changed by at least 90° from the average groundwater direction during baseflow. Aquifer heterogeneity in floodplains yields hyporheic flows that are more responsive and spatially and temporally complex than would be expected compared to more common assumptions of homogeneity. Copyright © 2012 John Wiley & Sons, Ltd.