A hydropedological approach to simulate streamflow and soil water contents with SWAT+

Reflecting internal catchment hydrological processes in hydrological models is important for accurate predictions of the impact of climate and land-use change on water resources. Characterizing these processes is however difficult and expensive due to their dynamic nature and spatio-temporal variability. Hydropedology is a relatively new discipline focusing on the synergistic integration of hydrology, soil physics and pedology. Hydropedological interpretations of soils and soil distribution can be used to characterize key hydrological processes, especially in areas with no or limited hydrometric measurements. Here we applied a hydropedological approach to reflect flowpaths through detailed routing in SWAT+ for a 157 ha catchment (Weatherley) in South Africa. We compared the hydropedological approach and a standard (no routing) approach against measured streamflow (two weirs) and soil water contents (13 locations). The catchment was treated as ‘ungauged’ and the model was not calibrated against hydrometric measurements in order to determine the direct contribution of hydropedology on modelling efficiency. Streamflow was predicted well without calibration (NSE > 0.8; R² > 0.82) for both approaches at both weirs. The standard approach yielded slightly better streamflow predictions. The hydropedological approach resulted in considerable improvements in the simulation of soil water contents (R² increased from 0.40 to 0.49 and PBIAS decreased from 40% to 20%). The routing capacity of SWAT+ as employed in the hydropedological approach reduced the underestimation of wetland water regimes drastically and resulted in a more accurate representation of the dominant hydrological processes in this catchment. We concluded that hydropedology can be a valuable source of ‘soft data’ to reflect internal catchment structure and processes and, potentially, for realistic calibrations in other studies, especially those conducted in areas with limited hydrometric measurements.     

Rainfall partitioning through a mixed cedar swamp and associated C and N fluxes in Southern Ontario,Canada

Partitioning of rainfall through a forest canopy into throughfall, stemflow, and canopy interception is a critical process in the water cycle, and the contact of precipitation with vegetated surfaces leads to increased delivery of solutes to the forest floor. This study investigates the rainfall partitioning over a growing season through a temperate, riparian, mixed coniferous‐deciduous cedar swamp, an ecosystem not well studied with respect to this process. Seasonal throughfall, stemflow, and interception were 69.2%, 1.5%, and 29.3% of recorded above‐canopy precipitation, respectively. Event throughfall ranged from a low of 31.5 ± 6.8% for a small 0.8‐mm event to a high of 82.9 ± 2.4% for a large 42.7‐mm event. Rain fluxes of at least 8 mm were needed to generate stemflow from all instrumented trees. Most trees had funnelling ratios <1.0, with an exponential decrease in funnelling ratio with increasing tree size. Despite this, stand‐scale funnelling ratios averaged 2.81 ± 1.73, indicating equivalent depth of water delivered across the swamp floor by stemflow was greater than incident precipitation. Throughfall dissolved organic carbon (DOC) and total dissolved nitrogen (TDN) averaged 26.60 ± 2.96 and 2.02 ± 0.16 mg L^?1, respectively, which were ~11 and three times above‐canopy rain levels. Stemflow DOC averaged 73.33 ± 7.43 mg L^?1, 35 times higher than precipitation, and TDN was 4.45 ± 0.56 mg L^?1, 7.5 times higher than rain. Stemflow DOC concentration was highest from Populus balsamifera and TDN greatest from Thuja occidentalis trees. Although total below‐canopy flux of TDN increased with increasing event size, DOC flux was greatest for events 20–30 mm, suggesting a canopy storage threshold of DOC was readily diluted. In addition to documenting rainfall partitioning in a novel ecosystem, this study demonstrates the excess carbon and nitrogen delivered to riparian swamps, suggesting the assimilative capacity of these zones may be underestimated.     

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.     

Changes in dissolved organic matter (DOM) amount and composition along nested headwater stream locations during baseflow and stormflow

Amount and composition of dissolved organic matter (DOM) were evaluated for multiple, nested stream locations in a forested watershed to investigate the role of hydrologic flow paths, wetlands and drainage scale. Sampling was performed over a 4‐year period (2008–2011) for five locations with drainage areas of 0.62, 3.5, 4.5, 12 and 79 ha. Hydrologic flow paths were characterized using an end‐member mixing model. DOM composition was determined using a suite of spectrofluorometric indices and a site‐specific parallel factor analysis model. Dissolved organic carbon (DOC), humic‐like DOM and fluorescence index were most sensitive to changes with drainage scale, whereas dissolved organic nitrogen, specific UV absorbance, Sr and protein‐like DOM were least sensitive. DOM concentrations and humic‐like DOM constituents were highest during both baseflow and stormflow for a 3.5‐ha catchment with a wetland near the catchment outlet. Whereas storm‐event concentrations of DOC and humic DOM constituents declined, the mass exports of DOC increased with increasing catchment scale. A pronounced dilution in storm‐event DOC concentration was observed at peak stream discharge for the 12‐ha drainage location, which was not as apparent at the 79‐ha scale, suggesting key differences in supply and transport of DOM. Our observations indicate that hydrologic flow paths, especially during storms, and the location and extent of wetlands in the catchment are key determinants of DOM concentration and composition. This study furthers our understanding of changes in DOM with drainage scale and the controls on DOM in headwater, forested catchments. Copyright © 2014 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.     

River evaporation,condensation and heat fluxes within a first‐order tributary of Catamaran Brook (New Brunswick,Canada)

Stream temperature plays an important role in many biotic and abiotic processes, as it influences many physical, chemical and biological properties in rivers. As such, a good understanding of the thermal regime of rivers is essential for effective fisheries management and the protection aquatic habitats. Moreover, a thorough understanding of underlying physical processes and river heat fluxes is essential in the development of better and more adaptive water temperature models. Very few studies have measured river evaporation and condensation and subsequently calculated corresponding heat fluxes in small tributary streams, mainly because microclimate data (data collected within the stream environment) are essential and rarely available. As such, the present study will address these issues by measuring river evaporation and condensation in tributary 1 (Trib 1, a small tributary within Catamaran Brook) using floating minipans. The latent heat flux and other important fluxes were calculated. Results showed that evaporation was low within the small Trib 1 of Catamaran Brook, less than 0.07 mm day^?1. Results showed that condensation played an important role in the latent heat flux. In fact, condensation was present during 34 of 92 days (37%) during the summer, which occurred when air temperature was greater than water temperature by 4–6 °C. Heat fluxes within this small stream showed that solar radiation dominated the heat gains and long‐wave radiation dominated the heat losses. © 2015 Her Majesty the Queen in Right of Canada. Hydrological Processes. © 2015 John Wiley & Sons, Ltd.     

Modelling groundwater/surface water interaction in a managed riparian chalk valley wetland

Understanding hydrological processes in wetlands may be complicated by management practices and complex groundwater/surface water interactions. This is especially true for wetlands underlain by permeable geology, such as chalk. In this study, the physically based, distributed model MIKE SHE is used to simulate hydrological processes at the Centre for Ecology and Hydrology River Lambourn Observatory, Boxford, Berkshire, UK. This comprises a 10‐ha lowland, chalk valley bottom, riparian wetland designated for its conservation value and scientific interest. Channel management and a compound geology exert important, but to date not completely understood, influences upon hydrological conditions. Model calibration and validation were based upon comparisons of observed and simulated groundwater heads and channel stages over an equally split 20‐month period. Model results are generally consistent with field observations and include short‐term responses to events as well as longer‐term seasonal trends. An intrinsic difficulty in representing compressible, anisotropic soils limited otherwise excellent performance in some areas. Hydrological processes in the wetland are dominated by the interaction between groundwater and surface water. Channel stage provides head boundaries for broad water levels across the wetland, whilst areas of groundwater upwelling control discrete head elevations. A relic surface drainage network confines flooding extents and routes seepage to the main channels. In‐channel macrophyte growth and its management have an acute effect on water levels and the proportional contribution of groundwater and surface water. The implications of model results for management of conservation species and their associated habitats are discussed. Copyright © 2015 John Wiley & Sons, Ltd.