Relating hydrogeomorphic properties to stream buffering chemistry in the Neversink River watershed,New York State,USA

Monitoring the effects of acidic deposition on aquatic ecosystems in the Northeastern US has generally required regular measurements of stream buffering chemistry (i.e. acid‐neutralizing capacity (ANC) and calcium Ca²⁺), which can be expensive and time consuming. The goal of this paper was to develop a simple method for predicting baseflow buffering chemistry based on the hydrogeomorphic properties of ten nested watersheds in the Neversink River basin (2·0–176·0 km²), an acid‐sensitive basin in the Catskill Mountains, New York State. The tributaries and main reach watersheds have strongly contrasting mean baseflow ANC values and Ca²⁺ concentrations, despite rather homogeneous vegetation, bedrock geology, and soils. A stepwise regression was applied to relate 13 hydrogeomorphic properties to the mean baseflow ANC values and Ca²⁺ concentrations. The regression analysis showed that watersheds with lower ANC values had a higher mean ratio of ‘quickflow’ runoff to precipitation during 20 non‐snowmelt runoff events (referred to as mean runoff ratio). The mean runoff ratio could explain at least 80% of the variability in mean baseflow ANC values and Ca²⁺ concentrations among the ten watersheds. Greater mean runoff ratios also correlated with steeper slopes and greater drainage densities, thus allowing the prediction of baseflow ANC values (r² = 0·75) and Ca²⁺ concentrations (r² = 0·77) with widely available spatial data alone. These results indicate that hydrogeomorphic properties can predict a watershed's sensitivity to acid deposition in regions where the spatial sources of stream buffering chemistry from the bedrock mineralogy and soils are fairly uniform. Copyright © 2010 John Wiley & Sons, Ltd.     

Catchment storage and residence time in a periodically irrigated watershed

Understanding anthropogenic impacts on water storage and water flow pathways in catchments is an ongoing challenge in hydrology. Here, we study the dynamics of subsurface storage and residence time of water in a catchment in Berkeley, California, that is within a regional park but contains diverse land use within its perimeter, including a periodically irrigated golf course. Our study combines several isotopic tracers with water budget data to examine sources of water in a stream draining the site. Irrigation water, applied to a small area of the watershed, is a minor component of the water budget. However, geochemical tracers reveal that irrigation water is a significant fraction of stream flow downstream of the golf course during baseflow and during precipitation events. Isotopic tracers indicate that the watershed has a preference to release young water for stream flow generation, resulting in contrasting tritium ages for stream water and groundwater of 1.3 ± 0.5 year and 8.2 ± 1.7 year, respectively. We determined that the older water is a very small component (0.7%) of the stream water in the tail of an assumed exponential distribution. We used the seasonal variation of stable water isotopes in precipitation and stream water over two water years to explain the damping of the isotopic signature of stream water, which yields information about the catchment's response to the input signal. The methods described here may be applicable to other urban or suburban headwater catchments in areas with a component of non-natural recharge from, for example, leaky infrastructure, storm water routing or dry season irrigation.     

High‐frequency precipitation and stream water quality time series from Plynlimon,Wales: an openly accessible data resource spanning the periodic table

This scientific briefing announces the availability of a new multi‐element high‐frequency water quality data set that is openly accessible to the research community. The data set comprises up to 2 years of 7‐hourly water quality data for two streams and one rainfall site in the Upper Severn catchment at Plynlimon in Mid‐Wales. The measurements cover 50 analytes ranging from H⁺ to U and spanning six orders of magnitude in concentration, including major, minor and trace elements as well as nutrients, and they complement decades of weekly measurements of the same analytes at the Upper Severn. Together, the weekly and 7‐hourly time series provide a unique data set for studying both long‐term trends and short‐term dynamics. The data show complex behaviour over a wide range of timescales, challenging our understanding of catchment processes and informing future modelling efforts. Copyright © 2013 John Wiley & Sons, Ltd.     

Characterizing ponds in a watershed simulation and evaluating their influence on streamflow in a Mississippi watershed

Small water bodies, such as ponds and wetlands, are common landscape features, but often are not simulated within a watershed modeling framework. The wetland modeling tool, AgWET, uses a GIS framework to characterize the features of ponds and wetlands so that they can be incorporated into watershed simulations using the Annualized Agricultural Non-Point Source (AnnAGNPS) pollution model. AgWET was used to characterize farm ponds on the Goodwin Creek Experimental Watershed in northwest Mississippi and AnnAGNPS simulated watershed hydrology. Monthly streamflow was validated at four watershed gauges with Nash-Sutcliffe efficiency values between 0.91 and 0.94. Ponds influenced watershed hydrology at various scales, with a decrease in average streamflow by 4% at the watershed outlet, 8% at the sub-watershed scale, and an average of 56% immediately downstream of the ponds. AgWET can be used to simulate ponds in watershed assessments for improved results and evaluation of future scenarios.     

The sensitivity of simulated flow and water quality response to spatial heterogeneity on a hillslope in the Tarrawarra catchment,Australia

Models are widely used to simulate hydrological response and the generation and transport of constituents such as salt, phosphorus, and nitrogen from catchments to streams. Several models use a spatial representation with catchments divided into subcatchments. Variations in land use and other characteristics within subcatchments are represented by spatially lumped hydrologic response units (HRUs) or functional units instead of using fully distributed models. This approach disregards any spatial interaction between HRUs, including their connectivity to each other and to the stream and the influence of these interactions on water and constituent export. A spatially explicit hydrological model (Thales) was used to simulate a variety of theoretical catchments with soils dominated by combinations of infiltration excess, saturation excess, and subsurface stormflow processes and different soil constituent concentrations that were spatially interacting (i.e. located along a hillslope sequence). The modelling results show that the response of both runoff and concentration is sensitive to varying spatial arrangements due to interactions of runoff, infiltration, and chemical processes between the different soil types in many but not all situations. Results highlight the importance of considering connectivity of pathways when modelling hydrological response and constituents export. This is achieved by comparing pairs of simulations and the corresponding differences in the exported loads. Copyright © 2009 John Wiley & Sons, Ltd.     

Water and salt balance modelling to predict the effects of land-use changes in forested catchments. 3. The large catchment model

This paper presents an application of a long-term, large catchment-scale, water balance model developed to predict the effects of forest clearing in the south-west of Western Australia. The conceptual model simulates the basic daily water balance fluxes in forested catchments before and after clearing. The large catchment is divided into a number of sub-catchments (1–5 km² in area), which are taken as the fundamental building blocks of the large catchment model. The responses of the individual subcatchments to rainfall and pan evaporation are conceptualized in terms of three inter-dependent subsurface stores A, B and F, which are considered to represent the moisture states of the subcatchments. Details of the subcatchment-scale water balance model have been presented earlier in Part 1 of this series of papers. The response of any subcatchment is a function of its local moisture state, as measured by the local values of the stores. The variations of the initial values of the stores among the subcatchments are described in the large catchment model through simple, linear equations involving a number of similarity indices representing topography, mean annual rainfall and level of forest clearing. The model is applied to the Conjurunup catchment, a medium-sized (39·6 km²) catchment in the south-west of Western Australia. The catchment has been heterogeneously (in space and time) cleared for bauxite mining and subsequently rehabilitated. For this application, the catchment is divided into 11 subcatchments. The model parameters are estimated by calibration, by comparing observed and predicted runoff values, over a 18 year period, for the large catchment and two of the subcatchments. Excellent fits are obtained.     

Assessing sensitivity of hydrologic responses to climate change from forested watershed in Mississippi

The climatic processes such as changes in precipitation, temperature and atmospheric CO₂ concentration can intensify the effects on water resources. An assessment of the effects of long‐term climate change on water resources is essential to the development of water quality improvement programs. This study was conducted in the Upper Pearl River Watershed (UPRW) in east‐central Mississippi to assess the effects of long‐term potential future climate change on average mean monthly stream flow from the five spatially distributed U. S. Geological Survey (USGS) gage stations in the UPRW using the Soil and Water Assessment Tool. The model was calibrated (January 1981 to December 1994) and validated (January 1995 to September 2008) using monthly measured stream flow data. The calibrated and validated model determined good to very good performance for stream flow prediction (R² and E from 0·60 to 0·86) between measured and predicted stream flow values. The root mean square error values (from 14 to 37 m³ s^?1) were estimated at similar levels of errors during model calibration and validation. The results showed that long‐term (50 years) average monthly stream flow sensitivity due to climate change effects was found the greatest as a result of percentage change in the precipitation followed by carbon dioxide (CO₂) concentration and temperature. The long‐term model simulation scenarios as compared with the base scenario for all five spatially distributed USGS gage stations in the UPRW estimated an average monthly stream flow decrease (from 54 to 67%) and average monthly stream flow increase (from 67 to 79%) depending on the spatial characteristics of the USGS gage stations. Overall, the results indicate that the UPRW hydrology is very sensitive to potential future climate changes and that these changes could stimulate increased streamflow generation from the watershed. Copyright © 2010 John Wiley & Sons, Ltd.     

Water and salt balance modelling to predict the effects of land-use changes in forested catchments. 2. coupled model of water and salt balances

A long-term salt balance model is coupled with the small catchment water balance model presented in Part 1 of this series of papers. The salt balance model was designed as a simple robust, conceptually based model of the fundamental salt fluxes and stores in forested and cleared catchments. The model has four interdependent stores representing salt storage in the unsaturated zone, the deep permanent saturated groundwater system, the near-stream perched groundwater system and in a ‘salt bulge’ just above the permanent water-table. The model has performed well in simulations carried out on Salmon and Wights, two small experimental catchments in south-west Western Australia. When applied to Wights catchment the salt balance model was able to predict the stream salinities prior to clearing of native forests, and the increased salinities after the clearing.     

Consistency between hydrological models and field observations: linking processes at the hillslope scale to hydrological responses at the watershed scale

The purpose of this paper is to identify simple connections between observations of hydrological processes at the hillslope scale and observations of the response of watersheds following rainfall, with a view to building a parsimonious model of catchment processes. The focus is on the well‐studied Panola Mountain Research Watershed (PMRW), Georgia, USA. Recession analysis of discharge Q shows that while the relationship between dQ/dt and Q is approximately consistent with a linear reservoir for the hillslope, there is a deviation from linearity that becomes progressively larger with increasing spatial scale. To account for these scale differences conceptual models of streamflow recession are defined at both the hillslope scale and the watershed scale, and an assessment made as to whether models at the hillslope scale can be aggregated to be consistent with models at the watershed scale. Results from this study show that a model with parallel linear reservoirs provides the most plausible explanation (of those tested) for both the linear hillslope response to rainfall and non‐linear recession behaviour observed at the watershed outlet. In this model each linear reservoir is associated with a landscape type. The parallel reservoir model is consistent with both geochemical analyses of hydrological flow paths and water balance estimates of bedrock recharge. Overall, this study demonstrates that standard approaches of using recession analysis to identify the functional form of storage–discharge relationships identify model structures that are inconsistent with field evidence, and that recession analysis at multiple spatial scales can provide useful insights into catchment behaviour. Copyright © 2008 John Wiley & Sons, Ltd.     

Quantitative evaluation of different hydrological modelling approaches in a partly glacierized Swiss watershed

Mountain water resources management often requires hydrological models that need to handle both snow and ice melt. In this study, we compared two different model types for a partly glacierized watershed in central Switzerland: (1) an energy‐balance model primarily designed for snow simulations; and (2) a temperature‐index model developed for glacier simulations. The models were forced with data extrapolated from long‐term measurement records to mimic the typical input data situation for climate change assessments. By using different methods to distribute precipitation, we also assessed how various snow cover patterns influenced the modelled runoff. The energy‐balance model provided accurate discharge estimations during periods dominated by snow melt, but dropped in performance during the glacier ablation season. The glacier melt rates were sensitive to the modelled snow cover patterns and to the parameterization of turbulent heat fluxes. In contrast, the temperature‐index model poorly reproduced snow melt runoff, but provided accurate discharge estimations during the periods dominated by glacier ablation, almost independently of the method used to distribute precipitation. Apparently, the calibration of this model compensated for the inaccurate precipitation input with biased parameters. Our results show that accurate estimates of snow cover patterns are needed either to correctly constrain the melt parameters of the temperature‐index model or to ensure appropriate glacier surface albedos required by the energy‐balance model. Thus, particularly when only distant meteorological stations are available, carefully selected input data and efficient extrapolation methods of meteorological variables improve the reliability of runoff simulations in high alpine watersheds. Copyright © 2011 John Wiley & Sons, Ltd.     

LISEM: A SINGLE-EVENT PHYSICALLY BASED HYDROLOGICAL AND SOIL EROSION MODEL FOR DRAINAGE BASINS. I: THEORY,INPUT AND OUTPUT

A new physically based hydrological and soil erosion model has been developed, which can be used for planning and conservation purposes: the LImburg Soil Erosion Model (LISEM). The LISEM model is one of the first examples of a physically based model that is completely incorporated in a raster Geographical Information System. This incorporation facilitates easy application in larger catchments, improves the user friendliness by avoiding conversion routines and allows remotely sensed data to be used. Processes incorporated in the model are rainfall, interception, surface storage in micro-depressions, infiltration and vertical movement of water in the soil, overland flow, channel flow, detachment by rainfall and throughfall, detachment by overland flow and transport capacity of the flow. Special attention has been given to the influence of tractor wheelings, small roads and surface sealing. Vertical movement of water in the soil is simulated using the Richard's equation. Optionally, the user can choose the Holtan or the Green–Ampt infiltration model. For the distribution flow routing, a four-point finite-difference solution of the kinematic wave is used together with Manning's equation.     

Snow hydrology of a constructed watershed in the Athabasca oil sands region,Alberta, Canada

Wetlands are now being integrated into oil sands mining landscape closure design plans. These wetland ecosystems will be constructed within a regional sub‐humid climate where snowfall represents ~25% of annual precipitation. However, few studies focus on the distribution of snow and, hence, the storage of winter precipitation in reclaimed oil sands landscapes. In this study, the distribution, ablation and fate of snowmelt waters are quantified within a constructed watershed in a post‐mining oil sands environment. Basin‐averaged peak SWE was 106 mm, with no significant difference between reclaimed slopes with vegetation and those that were sparsely vegetated or bare. Snow depth was greatest and more variable near the toe of slopes and became progressively shallower towards the crest. Snow ablation started first on the vegetated slope, which also exhibited the maximum observed ablation rates. This enhanced melt was attributed to increased absorption of short‐wave radiation by vegetation stems and branches. Recharge to reclaimed slopes and a constructed aquifer during the snowmelt period was minimal, as the presence of ground frost minimized infiltration. Accordingly, substantial surface run‐off was observed from all reclaimed slopes, despite being designed to reduce run‐off and increase water storage. This could result in increased flashiness of downstream watercourses during the spring freshet that receive run‐off from post‐mining landscapes where large reclaimed slopes are prolific. Run‐off ratios for the reclaimed slopes were between 0.7 and 0.9. Thus, it is essential to consider snow dynamics when designing landscape‐scale constructed ecosystems. This research demonstrates that the snowmelt period hydrology within reclaimed landscapes is fundamentally different from that reported for natural settings and represents one of the first studies on snow dynamics in constructed watershed systems in the post‐mined oil sands landscape. Copyright © 2016 John Wiley & Sons, Ltd.     

Does the incorporation of process conceptualization and tracer data improve the structure and performance of a simple rainfall‐runoff model in a Scottish mesoscale catchment?

A geomorphological instantaneous unit hydrograph (GIUH) rainfall‐runoff model was applied in a 31 km² montane catchment in Scotland. Modelling was based on flow path length distributions derived from a digital terrain model (DTM). The model was applied in two ways; a single landscape unit response based on the DTM alone, and a two‐landscape unit response, which incorporated the distribution of saturated areas derived from field‐validated geographic information system (GIS) analysis based on a DTM and soil maps. This was to test the hypothesis that incorporation of process‐information would enhance the model performance. The model was applied with limited multiple event calibration to produce parameter sets which could be applied to a spectrum of events with contrasting characteristics and antecedent conditions. Gran alkalinity was used as a tracer to provide an additional objective measure for assessing model performance. The models captured the hydrological response dynamics of the catchment reasonably well. In general, the single landscape unit approach produced the best individual model performance statistics, though the two‐landscape unit approach provided a range of models, which bracketed the storm hydrograph response more realistically. There was a tendency to over‐predict the rising limb of the hydrograph, underestimate large storm event peaks and anticipate the hydrograph recession too rapidly. Most of these limitations could be explained by the simplistic assumptions embedded within the GIUH approach. The modelling also gave feasible predictions of stream water chemistry, though these could not be used as a basis for model rejection. Nevertheless, the study suggested that the approach has potential for prediction of hydrological response in ungauged montane headwater basins. Copyright © 2007 John Wiley & Sons, Ltd.     

The Glendhu experimental catchment study,upland east Otago,New Zealand: 34 years of hydrological observations on the afforestation of tussock grasslands

This paper presents results from 34 years of the Glendhu Experimental Catchment Study, established in 1979 by the former New Zealand Forest Service in upland east Otago in New Zealand's South Island to determine the hydrological consequences of converting indigenous tussock grassland to plantation forestry. A traditional paired catchment approach was adopted; after a 2.5‐year pretreatment period, one catchment (GH2, 310 ha) was planted over two thirds of its area in Pinus radiata, and an adjacent catchment (GH1, 216 ha) was left in tussock as a control. The average annual reduction in water yield from the planted catchment between canopy closure in 1991 and 2013, compared with that in tussock, was 273 mm (33%). Annual water yields from the planted catchment continued to decline relative to the tussock catchment until 2010. Since then, the difference in annual water yields between the two catchments has narrowed. Ripping before planting caused some redistribution of the total streamflow from stormflow to baseflow. Following canopy closure, afforestation has reduced the low flow (Q₉₅) by an average of 26% compared with the tussock catchment. Average peak flows for small events (2–5 L/s/ha) were reduced by 78%, but only by 37% for larger, less frequent storms (>15 L/s/ha), suggesting that peak flows during high magnitude storms are less dependent on the prevailing land cover.     

IHMS—Integrated Hydrological Modelling System. Part 1. Hydrological processes and general structure

A newly Integrated Hydrological Modelling System (IHMS) has been developed to study the impact of changes in climate, land use and water management on groundwater and seawater intrusion (SWI) into coastal areas. The system represents the combination of three models, which can, if required, be run separately. It has been designed to assess the combined impact of climate, land use and groundwater abstraction changes on river, drainage and groundwater flows, groundwater levels and, where appropriate, SWI. The approach is interdisciplinary and reflects an integrated water management approach. The system comprises three packages: the Distributed Catchment Scale Model (DiCaSM), MODFLOW (96 and 2000) and SWI models. In addition to estimating all water balance components, DiCaSM, produces the recharge data that are used as input to the groundwater flow model of the US Geological Survey, MODFLOW. The latter subsequently generates the head distribution and groundwater flows that are used as input to the SWI model, SWI. Thus, any changes in land use, rainfall, water management, abstraction, etc. at the surface are first handled by DiCaSM, then by MODFLOW and finally by the SWI. The three models operate at different spatial and temporal scales and a facility (interface utilities between models) to aggregate/disaggregate input/output data to meet a desired spatial and temporal scale was developed allowing smooth and easy communication between the three models. As MODFLOW and SWI are published and in the public domain, this article focuses on DiCaSM, the newly developed unsaturated zone DiCaSM and equally important the interfacing utilities between the three models. DiCaSM simulates a number of hydrological processes: rainfall interception, evapotranspiration, surface runoff, infiltration, soil water movement in the root zone, plant water uptake, crop growth, stream flow and groundwater recharge. Input requirements include distributed data sets of rainfall, land use, soil types and digital terrain; climate data input can be either distributed or non‐distributed. The model produces distributed and time series output of all water balance components including potential evapotranspiration, actual evapotranspiration, rainfall interception, infiltration, plant water uptake, transpiration, soil water content, soil moisture (SM) deficit, groundwater recharge rate, stream flow and surface runoff. This article focuses on details of the hydrological processes and the various equations used in DiCaSM, as well as the nature of the interface to the MODFLOW and SWI models. Furthermore, the results of preliminary tests of DiCaSM are reported; these include tests related to the ability of the model to predict the SM content of surface and subsurface soil layers, as well as groundwater levels. The latter demonstrates how the groundwater recharge calculated from DiCaSM can be used as input into the groundwater model MODFLOW using aggregation and disaggregation algorithms (built into the interface utility). SWI has also been run successfully with hypothetical examples and was able to reproduce the results of some of the original examples of Bakker and Schaars ( 2005 ). In the subsequent articles, the results of applications to different catchments will be reported. Copyright © 2010 John Wiley & Sons, Ltd.     

CANOPEX: A Canadian hydrometeorological watershed database

Over the past few years, many international initiatives and collaborations were launched to improve and share knowledge in hydrology research. Large databases allowed finding patterns and relationships across regions and scales. This paper introduces the Canadian model parameter experiment (CANOPEX) database, which is adapted from the US MOPEX project data and methods. The CANOPEX database includes meteorological and hydrometric data as well as watershed boundaries for 698 basins. Two sets of basin‐averaged meteorological data (Maximum and minimum temperature and precipitation) are provided. The first dataset is directly taken from Environment Canada's weather stations whereas the second is extracted from the Natural Resources Canada gridded climate data product. Data are provided in MOPEX and Matlab formats. CANOPEX watersheds are well distributed over Canada, which allows investigating a variety of physiological and climatological conditions. The CANOPEX database can be used in a variety of hydrologic research projects such as climate change impact studies, model comparisons, multi‐modelling, ensemble streamflow prediction and model parameter estimation. CANOPEX could be used to generalize findings to other cold climate catchments as well as assess the robustness of research methodologies and procedures. Copyright © 2016 John Wiley & Sons, Ltd.     

An overview of hydrometeorological datasets from a small agricultural catchment (Nučice) in the Czech Republic

We introduce the freely available web-based Water in an Agricultural Landscape—NUčice Database (WALNUD) dataset that includes both hydrological and meteorological records at the Nučice experimental catchment (0.53 km²), which is representative of an intensively farmed landscape in the Czech Republic. The Nučice experimental catchment was established in 2011 for the observation of rainfall–runoff processes, soil erosion processes, and water balance of a cultivated landscape. The average altitude is 401 m a.s.l., the mean land slope is 3.9%, and the climate is humid continental (mean annual temperature 7.9°C, annual precipitation 630 mm). The catchment is drained by an artificially straightened stream and consists of three fields covering over 95% of the area which are managed by two different farmers. The typical crops are winter wheat, rapeseed, and alfalfa. The installed equipment includes a standard meteorological station, several rain gauges distributed across the basin, and a flume with an H-type facing that is used to monitor stream discharge, water turbidity, and basic water quality indicators. Additionally, the groundwater level and soil water content at various depths near the stream are recorded. Recently, large-scale soil moisture monitoring efforts have been introduced with the installation of two cosmic-ray neutron sensors for soil moisture monitoring. The datasets consist of observed variables (e.g. measured precipitation, air temperature, stream discharge, and soil moisture) and are available online for public use. The cross-seasonal, open access datasets at this small-scale agricultural catchment will benefit not only hydrologists but also local farmers.