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.     

Characterizing spatial and temporal variation in 18O and 2H content of New Zealand river water for better understanding of hydrologic processes

Time series of hydrogen and oxygen stable isotope ratios (δ²H and δ¹⁸O) in rivers can be used to quantify groundwater contributions to streamflow, and timescales of catchment storage. However, these isotope hydrology techniques rely on distinct spatial or temporal patterns of δ²H and δ¹⁸O within the hydrologic cycle. In New Zealand, lack of understanding of spatial and temporal patterns of δ²H and δ¹⁸O of river water hinders development of regional and national-scale hydrological models. We measured δ²H and δ¹⁸O monthly, together with river flow rates at 58 locations across New Zealand over a two-year period. Results show: (a) general patterns of decreasing δ²H and δ¹⁸O with increasing latitude were altered by New Zealand's major mountain ranges; δ²H and δ¹⁸O were distinctly lower in rivers fed from higher elevation catchments, and in eastern rain-shadow areas of both islands; (b) river water δ²H and δ¹⁸O values were partly controlled by local catchment characteristics (catchment slope, PET, catchment elevation, and upstream lake area) that influence evaporation processes; (c) regional differences in evaporation caused the slope of the river water line (i.e., the relationship between δ²H and δ¹⁸O in river water) for the (warmer) North Island to be lower than that of the (cooler, mountain-dominated) South Island; (d) δ²H seasonal offsets (i.e., the difference between seasonal peak and mean values) for individual sites ranged from 0.50‰ to 5.07‰. Peak values of δ¹⁸O and δ²H were in late summer, but values peaked 1 month later at the South Island sites, likely due to greater snow-melt contributions to streamflow. Strong spatial differences in river water δ²H and δ¹⁸O caused by orographic rainfall effects and evaporation may inform studies of water mixing across landscapes. Generally distinct seasonal isotope cycles, despite the large catchment sizes of rivers studied, are encouraging for transit time analysis applications.     

Water uptake of riparian plants in the lower Lhasa River Basin,South Tibetan Plateau using stable water isotopes

Riparian plants can adapt their water uptake strategies based on climatic and hydrological conditions within a river basin. The response of cold-alpine riparian trees to changes in water availability is poorly understood. The Lhasa River is a representative cold-alpine river in South Tibet and an under-studied environment. Therefore, a 96 km section of the lower Lhasa River was selected for a study on the water-use patterns of riparian plants. Plant water, soil water, groundwater and river water were measured at three sites for δ¹⁸O and δ²H values during the warm-wet and cold-dry periods in 2018. Soil profiles differed in isotope values between seasons and with the distance along the river. During the cold-dry period, the upper parts of the soil profiles were significantly affected by evaporation. During the warm-wet period, the soil profile was influenced by precipitation infiltration in the upper reaches of the study area and by various water sources in the lower reaches. Calculations using the IsoSource model indicated that the mature salix and birch trees (Salix cheilophila Schneid. and Betula platyphylla Suk.) accessed water from multiple sources during the cold-dry period, whereas they sourced more than 70% of their requirement from the upper 60–80 cm of the soil profile during the warm-wet period. The model indicated that the immature rose willow tree (Tamarix ramosissima Ledeb) accessed 66% of its water from the surface soil during the cold-dry period, but used the deeper layers during the warm-wet period. The plant type was not the dominant factor driving water uptake patterns in mature plants. Our findings can contribute to strategies for the sustainable development of cold-alpine riparian ecosystems. It is recommended that reducing plantation density and collocating plants with different rooting depths would be conducive to optimal plant growth in this environment.     

Application of a confluence‐based sediment‐fingerprinting approach to a dynamic sedimentary catchment,New Zealand

Fine sediment is a dynamic component of the fluvial system, contributing to the physical form, chemistry and ecological health of a river. It is important to understand rates and patterns of sediment delivery, transport and deposition. Sediment fingerprinting is a means of directly determining sediment sources via their geochemical properties, but it faces challenges in discriminating sources within larger catchments. In this research, sediment fingerprinting was applied to major river confluences in the Manawatu catchment as a broad‐scale application to characterizing sub‐catchment sediment contributions for a sedimentary catchment dominated by agriculture. Stepwise discriminant function analysis and principal component analysis of bulk geochemical concentrations and geochemical indicators were used to investigate sub‐catchment geochemical signatures. Each confluence displayed a unique array of geochemical variables suited for discrimination. Geochemical variation in upstream sediment samples was likely a result of the varying geological source compositions. The Tiraumea sub‐catchment provided the dominant signature at the major confluence with the Upper Manawatu and Mangatainoka sub‐catchments. Subsequent downstream confluences are dominated by the upstream geochemical signatures from the main stem of Manawatu River. Variability in the downstream geochemical signature is likely due to incomplete mixing caused in part by channel configuration. Results from this exploratory investigation indicate that numerous geochemical elements have the ability to differentiate fine sediment sources using a broad‐scale confluence‐based approach and suggest there is enough geochemical variation throughout a large sedimentary catchment for a full sediment fingerprint model. Combining powerful statistical procedures with other geochemical analyses is critical to understanding the processes or spatial patterns responsible for sediment signature variation within this type of catchment. Copyright © 2015 John Wiley & Sons, Ltd.     

Water source dynamics of high Arctic river basins

Arctic river basins are amongst the most vulnerable to climate change. However, there is currently limited knowledge of the hydrological processes that govern flow dynamics in Arctic river basins. We address this research gap using natural hydrochemical and isotopic tracers to identify water sources that contributed to runoff in river basins spanning a gradient of glacierization (0–61%) in Svalbard during summer 2010 and 2011. Spatially distinct hydrological processes operating over diurnal, weekly and seasonal timescales were characterized by river hydrochemistry and isotopic composition. Two conceptual water sources (‘meltwater’ and ‘groundwater’) were identified and used as a basis for end‐member mixing analyses to assess seasonal and year‐to‐year variability in water source dynamics. In glacier‐fed rivers, meltwater dominated flows at all sites (typically >80%) with the highest contributions observed at the beginning of each study period in early July when snow cover was most extensive. Rivers in non‐glacierized basins were sourced initially from snowmelt but became increasingly dependent on groundwater inputs (up to 100% of total flow volume) by late summer. These hydrological changes were attributed to the depletion of snowpacks and enhanced soil water storage capacity as the active layer expanded throughout each melt season. These findings provide insight into the processes that underpin water source dynamics in Arctic river systems and potential future changes in Arctic hydrology that might be expected under a changing climate. Copyright © 2013 John Wiley & Sons, Ltd.     

Quantification of river water infiltration in shallow aquifers using acesulfame and anthropogenic gadolinium

This study has investigated the use of the artificial sweetener acesulfame and the magnetic resonance imaging contrast agent gadolinium as quantitative tracers for river water infiltration into shallow groundwater. The influence of a river on alluvial groundwater in a subalpine catchment in western Europe has been assessed using the ‘classical’ hydrochemical tracer chloride and the trace contaminants acesulfame and anthropogenic gadolinium. Mixing ratios for riverine bank filtrate with ambient groundwater and the uncertainties associated with the temporal and spatial tracer variability were calculated using acesulfame and gadolinium and compared with those obtained using chloride. The temporal variability of tracer concentrations in river water of gadolinium (standard deviation SD: 63%) and acesulfame (SD: 71%) both exceeded that of chloride (SD: 27%), and this was identified as the main source of uncertainty in the mixing analysis. Similar spatial distributions were detected in the groundwater for chloride and gadolinium, but not for acesulfame. Mixing analyses using acesulfame resulted in calculated mixing ratios that differed from those obtained using gadolinium and chloride by up to 83% and 92%, respectively. At the investigated site, which had oxic conditions and moderate temperatures, acesulfame was found to be a less reliable tracer than either gadolinium or chloride, probably because of natural attenuation and input from other sources. There was no statistically significant difference between the mixing ratios obtained using chloride or gadolinium, the mixing ratios obtained using gadolinium were 40–50% lower than those obtained using chloride. This is mainly due to a bias of the mean gadolinium concentration in river water towards higher values. In view of the uncertainties of the two tracers, neither could be preferred over the other for the quantification of bank filtrate in groundwater. At this specific site gadolinium was able to reliably identify river water infiltration and was a more precise tracer than chloride at low mixing ratios (<20%), because of the exclusive occurrence of gadolinium in river water and its high dynamic range. Copyright © 2015 John Wiley & Sons, Ltd.     

Scales and magnitude of hyporheic,river–aquifer and bank storage exchange fluxes

Many studies have investigated the exchange processes that occur between rivers and groundwater systems and have successfully quantified the water fluxes involved. Specifically, these exchange processes include hyporheic exchange, river–aquifer exchange (groundwater discharge and river loss) and bank storage exchange. Remarkably, there are relatively few examples of field studies where more than one exchange process is quantified, and as a consequence, the relationships between them are not well understood. To compare the relative magnitudes of these common exchange processes, we have collected data from 54 studies that have quantified one or more of these exchange flux types. Each flux value is plotted against river discharge at the time of measurement to allow the different exchange flux types to be compared. We show that there are positive relationships between the magnitude of each exchange flux type and increasing river discharge across the different studies. For every one order of magnitude increase in river discharge, the hyporheic, river–aquifer and bank storage exchange fluxes increase by factors of 2.7, 2.9 and 2.5, respectively. On average, hyporheic exchange fluxes are almost an order of magnitude greater than river–aquifer exchange fluxes, which are, in turn, approximately four times greater than bank storage exchange fluxes for the same river discharge. Unless measurement approaches that can distinguish between different types of exchange flux are used, there is potential for hyporheic exchange fluxes to be misinterpreted as river–aquifer exchange fluxes, with possible implications for water resource management decisions. Copyright © 2015 John Wiley & Sons, Ltd.     

Reach‐scale contributions of road‐surface sediment to the Honna River,Haida Gwaii,BC

Gravel road surfaces can be a major source of fine sediment to streams, yet their contribution to channel reach sediment balances remains poorly documented. To quantify the input of road surface material and to compare this input with natural sediment sources at the reach scale, suspended sediment dynamics was examined and a 16‐month sediment balance was developed for a ~35 channel‐width (approx. 425 m) reach of the Honna River, a medium‐size, road‐affected stream located in coastal British Columbia. Of the 105 ± 33 t of suspended material passing through the reach, 18 ± 6% was attributed to the road surface. The high availability of sediment on the road surface appears to limit hysteresis in road run‐off. During rainstorms that increase streamflow, road surface material composed 0.5–15% of sediment inputs during relatively dry conditions from April to the end of September and 5–70% through wetter conditions from October to the end of March, but our data do not show evidence of major sediment accumulation on the riverbed in the reach. A comparison of modelled sediment production on the road surface with observed yields from drainage channels suggests that (1) during low intensity rainfall, ditches and drainage channels may trap sediment from road run‐off, which is subsequently released during events of greater intensity, and/or (2) production models do not effectively describe processes, such as deposition or erosion of sediment in ditches, which control sediment transport and delivery. Our findings further emphasize the risk of unpaved roads in polluting river systems and highlight the continued need for careful road design and location away from sensitive aquatic environments. Copyright © 2016 John Wiley & Sons, Ltd.     

Watershed delineation using hydrographic features and a DEM in plain river network region

Watershed delineation is a required step when conducting any spatially distributed hydrological modelling. Automated approaches are often proposed to delineate a watershed based on a river network extracted from the digital elevation model (DEM) using the deterministic eight‐neighbour (D8) method. However, a realistic river network cannot be derived from conventional DEM processing methods for a large flat area with a complex network of rivers, lakes, reservoirs, and polders, referred to as a plain river network region (PRNR). In this study, a new approach, which uses both hydrographic features and DEM, has been developed to address the problems of watershed delineation in PRNR. It extracts the river nodes and determines the flow directions of the river network based on a vector‐based hydrographic feature data model. The river network, lakes, reservoirs, and polders are then used to modify the flow directions of grid cells determined by D8 approach. The watershed is eventually delineated into four types of catchments including lakes, reservoirs, polders, and overland catchments based on the flow direction matrix and the location of river nodes. Multiple flow directions of grid cells are represented using a multi‐direction encoding method, and multiple outflows of catchments are also reflected in the topology of catchments. The proposed approach is applied to the western Taihu watershed in China. Comparisons between the results obtained from the D8 approach, the ‘stream burning’ approach, and those from the proposed approach clearly demonstrate an improvement of the new approach over the conventional approaches. This approach will benefit the development of distributed hydrological models in PRNR for the consideration of different types and multiple inlets and outlets of catchments. Copyright © 2015 John Wiley & Sons, Ltd.