Simulated response of an intermittent stream to rainfall frequency patterns

The response of intermittent catchments to rainfall is complex and difficult to model. This study uses the spatially distributed CATchment HYdrology (CATHY) model to explore how the frequency of daily rainfall (λ) can affect the hydrologic regime of intermittent catchments. After a multi-objective calibration and validation of CATHY against experimental measurements of streamflow and groundwater levels in a catchment used as a pasture, the role of λ in affecting streamflow characteristics was explored using different scenarios. With different values of λ for the dry and wet periods of the year, CATHY showed that a series of frequent rainfall events was often associated with incipient streamflow, independent of the season. Activation of streamflow during the wet season was related to multiple factors and was not often associated with the shallow groundwater levels near the outlet of the catchment. The interplay between rainfall depth and intensity acted as the most important factor for the generation of streamflow. Using the difference between accumulated rainfall and evapotranspiration as a measure of wetness, saturated subsurface flow mechanism generated streamflow in simulations with wetness at least three times larger than mean wetness of other simulations. Although groundwater uprise near the outlet did not effectively contribute to streamflow in the initial days of flow, it strongly correlated with the magnitude of the runoff coefficient. Values of λ close or equal to the maximum value in the wet season can sustain the connectivity between groundwater and streamflow in the riparian zone. This connectivity increases the catchment wetness, which consequently results in an increase of the generated streamflow. Our study showed that rainfall regimes characterized by different λ were able to identify distinct flow regimes typical of either intermittent, ephemeral, or nonflowing catchments. Decrease of λ in the wet season is likely associated with a reduction of streamflow, with a shift of flow regime from intermittent to ephemeral or no-flow.     

RELATIONSHIPS BETWEEN CATCHMENT ATTRIBUTES AND HYDROLOGICAL RESPONSE CHARACTERISTICS IN SMALL AUSTRALIAN MOUNTAIN ASH CATCHMENTS

Sixteen small catchments in the Maroondah region of Victoria, Australia were analysed using rainfall, temperature and streamflow time series with a rainfall–runoff model whose parameters efficiently characterize the hydrological response of a catchment. A set of catchment attributes for each of these catchments was then compared with the associated set of hydrological response characteristics of the catchments as estimated by the model. The time constant governing quickflow recession of streamflow (τ_q) was related to the drainage network and catchment area. The time constant governing slowflow recession of streamflow (τ_s) was related to the slope and shape of the catchment. The parameter governing evapotranspirative losses ( f ) was related to catchment gradient and vegetative water use. Forestry activities in the catchments changed evapotranspirative losses and thus total volume of streamflow, but did not affect the rate of streamflow recession.     

MODELLING RAINFALL-RUNOFF FROM LARGE CATCHMENT TO BASIN SCALE: THE GOULBURN VALLEY,VICTORIA

A lumped parameter dynamic rainfall-runoff model, IHACRES, is applied to the large upland area (more than 4500 km²) of the Goulburn Valley Basin, Victoria, Australia to predict streamflow under different climatic conditions. This paper presents the first evaluation of a rainfall–runoff model at large catchment scale, which is comprehensive in terms of the number of catchments investigated and the number of calibration and simulation periods used. The basin is subdivided into 12 catchments (from 100 to 700 km²), each of which is calibrated separately. High values of model efficiency and low bias are consistently obtained for different calibration sub-periods for all catchments in the basin. Simulation or so-called validation tests are used to select the best models for each catchment. This allows simulation of the water regime during long historical (approximately 90 year) periods when only climatological (rainfall and temperature) data were available. This procedure is extremely important for the estimation of the effect of climate variability and of the possible impact of climate change on the hydrological regime in the region and, in particular, for supporting irrigation management of the basin. Analysis of a composite catchment (2417 km²) and its five separate subcatchments indicates that the information content in the rainfall–streamflow data is independent of catchment size. Dynamic modelling of the daily water balance at the macroscale is limited principally by the adequacy of the precipitation gauging network. When a good estimate of areal precipitation is available for a catchment, it is not necessary to consider subcatchment-scale variability for modelling if the only interest is the daily discharge and evaporation losses from the catchment.     

Spatial GR4J conceptualization of the Tamor glaciated alpine catchment in Eastern Nepal: evaluation of GR4JSG against streamflow and MODIS snow extent

Snow and glacial melt processes are an important part of the Himalayan water balance. Correct quantification of melt runoff processes is necessary to understand the region's vulnerability to climate change. This paper describes in detail an application of conceptual GR4J hydrological model in the Tamor catchment in Eastern Nepal using typical elevation band and degree‐day factor approaches to model Himalayan snow and glacial melt processes. The model aims to provide a simple model that meets most water planning applications. The paper contributes a model conceptualization (GR4JSG) that enables coarse evaluation of modelled snow extents against remotely sensed Moderate Resolution Imaging Spectroradiometer snow extent. Novel aspects include the glacial store in GR4JSG and examination of how the parameters controlling snow and glacial stores correlate with existing parameters of GR4J. The model is calibrated using a Bayesian Monte Carlo Markov Chain method against observed streamflow for one glaciated catchment with reliable data. Evaluation of the modelled streamflow with observed streamflow gave Nash Sutcliffe Efficiency of 0.88 and Percent Bias of <4%. Comparison of the modelled snow extents with Moderate Resolution Imaging Spectroradiometer gave R² of 0.46, with calibration against streamflow only. The contribution of melt runoff to total discharge from the catchment is 14–16% across different experiments. The model is highly sensitive to rainfall and temperature data, which suffer from known problems and biases, for example because of stations being located predominantly in valleys and at lower elevations. Testing of the model in other Himalayan catchments may reveal additional limitations. © 2016 The Authors. Hydrological Processes published by John Wiley & Sons Ltd.     

Estimation of sub-annual inter-catchment groundwater flow using short-term water balance method

Predicting inter-catchment groundwater flow (IGF) is essential because IGF greatly affects stream water discharge and water chemistry. However, methods for estimating sub-annual IGF and clarifying its mechanisms using minimal data are limited. Thus, we quantified the sub-annual IGF and elucidated its driving factors using the short-term water balance method (STWB) for three forest headwater catchments in Japan (named here catchment A, B and As). Our previous study using the chloride mass balance indicated that annual IGF of catchment A (49.0 ha) can be negligible. Therefore, we calculated the daily evapotranspiration (ET) rate using the Priestley–Taylor expression and the 5-year water balance in catchment A (2010–2014). The sub-annual IGF of the three catchments was then calculated by subtracting the ET rate from the difference between rainfall and stream discharge during the sub-annual water balance periods selected using the STWB. The IGF rates of catchment B (7.0 ha), which is adjacent to catchment A, were positive in most cases, indicating that more groundwater flowed out of the catchment than into it, and exhibited positive linear relationships with rainfall and stream discharge. This suggested that as the catchments became wetter, more groundwater flowed out of catchment B. Conversely, the IGF rates of catchment As (5.3 ha), included in catchment A, were negative in most cases, indicating that more groundwater flowed into the catchment than out from it, and exhibited negative linear relationships with rainfall and stream discharge. Given the topography of the catchments studied, infiltration into the bedrock was the probable reason for the IGF outflow from catchment B. We hypothesized that in catchment As, the discrepancy between the actual hydrological boundary and the surface topographic boundary could have caused an IGF inflow. This study provides a useful tool for determining an IGF model structure to be incorporated into rainfall-runoff models.     

Estimation of rainfall elasticity of streamflow in Australia

Abstract

Estimates of rainfall elasticity of streamflow in 219 catchments across Australia are presented. The rainfall elasticity of streamflow is defined here as the proportional change in mean annual streamflow divided by the proportional change in mean annual rainfall. The elasticity is therefore a simple estimate of the sensitivity of long-term streamflow to changes in long-term rainfall, and is particularly useful as an initial estimate of climate change impact in land and water resources projects. The rainfall elasticity of streamflow is estimated here using a hydrological modelling approach and a nonparametric estimator. The results indicate that the rainfall elasticity of streamflow (? _P ) in Australia is about 2.0–3.5 (observed in about 70% of the catchments), that is, a 1% change in mean annual rainfall results in a 2.0–3.5% change in mean annual streamflow. The rainfall elasticity of streamflow is strongly correlated to runoff coefficient and mean annual rainfall and streamflow, where streamflow is more sensitive to rainfall in drier catchments, and those with low runoff coefficients. There is a clear relation-ship between the ? _P values estimated using the hydrological modelling approach and those estimated using the nonparametric estimator for the 219 catchments, although the values estimated by the hydrological modelling approach are, on average, slightly higher. The modelling approach is useful where a detailed study is required and where there are sufficient data to reliably develop and calibrate a hydrological model. The nonparametric estimator is useful where consistent estimates of the sensitivity of long-term streamflow to climate are required, because it is simple to use and estimates the elasticity directly from the historical data. The nonparametric method, being model independent, can also be easily applied in comparative studies to data sets from many catchments across large regions.     

Evaluating the SWAT's snow hydrology over a Northern Quebec watershed

The snow treatment becomes an important component of Soil and Water Assessment Tool (SWAT)’s hydrology when spring flows are dominated by snow melting. However, little is known about SWAT's snow hydrology performance because most studies using SWAT were conducted in rainfall‐driven catchments. To fill this gap, the present study aims to evaluate the ability of SWAT in simulating snow‐melting‐dominated streamflow in the Outardes Basin in Northern Quebec. SWAT performance in simulating snowmelt is evaluated against observed streamflow data and compared to simulations from the operationally used Streamflow Synthesis and Reservoir Regulation (SSARR) model over that catchment. The SWAT 5‐year calibration showed a satisfactory performance at the daily and seasonal time scales with low volume biases. The SWAT validation was conducted over two (17‐year and 15‐year) periods. Performances were similar to the calibration period in simulating the daily and seasonal streamflows again with low model biases. The spring‐snowmelt‐generated peak flow was accurately simulated by SWAT both in magnitude and timing. When SWAT's results are compared to SSARR, similar performances in simulating the daily discharges were observed. SSARR simulates more accurately streamflow generated at the snowmelt onset whereas SWAT better predicts streamflow in summer, fall and winter. SWAT provided reasonable streamflow simulations for our snow‐covered catchment, but refinement of the process‐driven baseflow during the snowmelt onset could improve spring performances. Therefore, SWAT becomes an attractive tool for evaluating water resources management in Nordic environments when a distributed model is preferred or when water quality information (e.g. temperature) is required. Copyright © 2013 John Wiley & Sons, Ltd.     

A comparison of streamflow,salt and water balances in adjacent farmland and forest catchments in south‐western Victoria,Australia

A study of the hydrologic effects of catchment change from pasture to plantation was carried out in Gatum, south‐western Victoria, Australia. This study describes the hydrologic characteristics of two adjacent catchments: one with 97% grassland and the other one with 62% Eucalyptus globulus plantations. Streamflow from both catchments was intermittent during the 20‐month study period. Monthly streamflow was always greater in the pasture‐dominated catchment compared with the plantation catchment because of lower evapotranspiration in the pasture‐based catchment. This difference in streamflow was also observed even during summer 2010/2011 when precipitation was 74% above average (1954–2012) summer rainfall. Streamflow peaks in the plantation‐based catchment were smaller than in the pasture‐dominated system. Flow duration curves show differences between the pasture and plantation‐dominated catchments and affect both high‐flow and low‐flow periods. Groundwater levels fell (up to 4.4 m) in the plantation catchment during the study period but rose (up to 3.2 m) in the pasture catchment. Higher evapotranspiration in the plantation catchment resulted in falling groundwater levels and greater disconnection of the groundwater system from the stream, resulting in lower baseflow contribution to streamflow. Salt export from each catchment increases with increasing flow and is higher at the pasture catchment, mainly because of the higher flow. Reduced salt loading to streams due to tree planting is generally considered environmentally beneficial in saline areas of south‐eastern Australia, but this benefit is offset by reduced total streamflow. Copyright © 2014 John Wiley & Sons, Ltd.     

Change in low flows due to catchment management dynamics—Application of a comparative modelling approach

Understanding the natural low flow of a catchment is critical for effective water management policy in semi-arid and arid lands. The Geba catchment in Ethiopia, forming the headwaters of Tekeze-Atbara basin was known for its severe land degradation before the recent large scale Soil and Water conservation (SWC) programs. Such interventions can modify the hydrological processes by changing the partitioning of the incoming rainfall on the land surface. However, the literature lacks studies to quantify the hydrological impacts of these interventions in the semi-arid catchments of the Nile basin. Statistical test and Indicators of Hydrological Alteration (IHA) were used to identify the trends of streamflow in two comparatives adjacent (one treated with intensive SWC intervention and control with fewer interventions) catchments. A distributed hydrological model was developed to understand the differences in hydrological processes of the two catchments. The statistical and IHA tools showed that the low flow in the treated catchment has significantly increased while considerably decreased in the control catchment. Comparative analysis confirmed that the low flow in the catchment with intensive SWC works was greater than that of the control by >30% while the direct runoff was lower by >120%. This implies a large proportion of the rainfall in the treated catchment is infiltrated and recharge aquifers which subsequently contribute to streamflow during the dry season. The proportion of soil storage was more than double compared to the control catchment. Moreover, hydrological response comparison from pre- and post-intervention showed that a drastic reduction in direct runoff (>84%) has improved the low flow by >55%. This strongly suggests that the ongoing intensive SWC works have significantly improved the low flows while it contributed to the reduction of total streamflow in the catchment.     

Understanding changes in annual runoff following land use changes: a systematic data‐based approach

A simple modelling framework for assessing the response of ungauged catchments to land use change in South‐Western Australia is presented. The framework uses knowledge of transpiration losses from native vegetation and pasture and then partitions the ‘excess’ water (resulting from reduced transpiration after land use change) between runoff and deep storage. The simple partitioning is achieved by using soft information (satellite imagery, previous mapping and field assessment) to delimit the spread of the permanently saturated area close to the stream. Runoff is then assumed to increase in proportion to the saturated area, with the residual difference going to deep storage. The model parameters to describe the annual water yield are obtained a priori and no calibration to streamflow is required. We tested the model using gauged records over 25 years from paired catchment experiments in South‐Western Australia. Very good estimates of runoff were obtained from high rainfall (>1100 mm yr⁻¹) catchments (R² > 0·9) and for low rainfall (<900 mm yr⁻¹) catchments after clearing (R² = 0·96) but results were poorer (R² = 0·55) for an uncleared low rainfall catchment, although overall balances were reasonable. In the drier uncleared catchments, the within‐year distributions of rainfall may exert a substantial influence on runoff response that is not completely captured by the presented model. Copyright © 2005 John Wiley & Sons, Ltd.     

Water and salt balance modelling to predict the effects of land-use changes in forested catchments. 1. Small catchment water balance model

A long-term water balance model has been developed to predict the hydrological effects of land-use change (especially forest clearing) in small experimental catchments in the south-west of Western Australia. This small catchment model has been used as the building block for the development of a large catchment-scale model, and has also formed the basis for a coupled water and salt balance model, developed to predict the changes in stream salinity resulting from land-use and climate change. The application of the coupled salt and water balance model to predict stream salinities in two small experimental catchments, and the application of the large catchment-scale model to predict changes in water yield in a medium-sized catchment that is being mined for bauxite, are presented in Parts 2 and 3, respectively, of this series of papers. The small catchment model has been designed as a simple, robust, conceptually based model of the basic daily water balance fluxes in forested catchments. The responses of the catchment to rainfall and pan evaporation are conceptualized in terms of three interdependent subsurface stores A, B and F. Store A depicts a near-stream perched aquifer system; B represents a deeper, permanent groundwater system; and F is an intermediate, unsaturated infiltration store. The responses of these stores are characterized by a set of constitutive relations which involves a number of conceptual parameters. These parameters are estimated by calibration by comparing observed and predicted runoff. The model has performed very well in simulations carried out on Salmon and Wights, two small experimental catchments in the Collie River basin in south-west Western Australia. The results from the application of the model to these small catchments are presented in this paper.     

Modelling daily streamflow at ungauged catchments: what information is necessary?

Rainfall–runoff modelling at ungauged catchments often involves the transfer of calibrated model parameters from ‘donor’ gauged catchments. However, in any rainfall–runoff model, some parameters tend to be more sensitive to the objective function, whereas others are insensitive over their entire feasible range. In this paper, we analyse the effect of selectively transferring sensitive versus insensitive parameters on streamflow predictability at ungauged catchments. We develop a simple daily time‐step rainfall–runoff model [exponential bucket hydrologic model (EXP‐HYDRO)] and calibrate it at 756 catchments within the continental USA. Nash–Sutcliffe efficiency of (NS) is used as the objective function. The model simulates satisfactorily at 323 catchments (NS > 0.6), most of which are located in the eastern part of the USA, along the Rocky Mountain Range, and near the western Pacific coast. Of the six calibration parameters, only three parameters are found to be sensitive to NS. Two of these parameters control the hydrograph recession behaviour of a catchment, and the third parameter controls the snowmelt rate. We find that when only sensitive parameters are transferred, model performance at ungauged catchments is almost at par with that of transferring all six parameters. Conversely, the transfer of only insensitive parameters results in a significant deterioration in model performance. Results suggest that streamflow predictability at ungauged catchments using rainfall–runoff models is largely dependent on the transfer of a small subset of parameters. We recommend that, in any modelling framework, such parameters should be identified and further characterized to better understand the information controlling streamflow predictability at ungauged catchments. Copyright © 2012 John Wiley & Sons, Ltd.     

Multi‐variable and multi‐site calibration and validation of SWAT in a large mountainous catchment with high spatial variability

Many methods developed for calibration and validation of physically based distributed hydrological models are time consuming and computationally intensive. Only a small set of input parameters can be optimized, and the optimization often results in unrealistic values. In this study we adopted a multi‐variable and multi‐site approach to calibration and validation of the Soil Water Assessment Tool (SWAT) model for the Motueka catchment, making use of extensive field measurements. Not only were a number of hydrological processes (model components) in a catchment evaluated, but also a number of subcatchments were used in the calibration. The internal variables used were PET, annual water yield, daily streamflow, baseflow, and soil moisture. The study was conducted using an 11‐year historical flow record (1990–2000); 1990–94 was used for calibration and 1995–2000 for validation. SWAT generally predicted well the PET, water yield and daily streamflow. The predicted daily streamflow matched the observed values, with a Nash–Sutcliffe coefficient of 0·78 during calibration and 0·72 during validation. However, values for subcatchments ranged from 0·31 to 0·67 during calibration, and 0·36 to 0·52 during validation. The predicted soil moisture remained wet compared with the measurement. About 50% of the extra soil water storage predicted by the model can be ascribed to overprediction of precipitation; the remaining 50% discrepancy was likely to be a result of poor representation of soil properties. Hydrological compensations in the modelling results are derived from water balances in the various pathways and storage (evaporation, streamflow, surface runoff, soil moisture and groundwater) and the contributions to streamflow from different geographic areas (hill slopes, variable source areas, sub‐basins, and subcatchments). The use of an integrated multi‐variable and multi‐site method improved the model calibration and validation and highlighted the areas and hydrological processes requiring greater calibration effort. Copyright © 2005 John Wiley & Sons, Ltd.     

The use of simple process-based models in the estimate of water balances for mixed land use catchments in East Africa

Soil moisture measurements by the neutron probe method were analysed to provide the parameters required for a daily model of actual evaporation from three land uses—grassland, indigenous bamboo and plantation softwood—in the Aberdare range of hills, Kenya. These daily estimates of evaporation were summed to provide annual totals and used, on a percentage land cover basis, in water balance calculations for three experimental mixed land use catchments, two of which were undergoing land use change. The annual water use, given by the difference between rainfall inputs and streamflow outputs, of the undisturbed catchment could normally be predicted to within 10%, whereas differences in the predicted and measured water use of the other two catchments were related to the changes in vegetation.     

Streamflow response to Pinus plantation harvesting: Canobolas State forest,southeastern Australia

Streamflows were measured in two Pinus radiata plantation catchments and one native eucalypt forest catchment in Canobolas State forest from 1999 to 2007. In 2002/2003, clearfall harvesting of 43·2 and 40·3% of two plantation catchments occurred, respectively. Water yields increased by 54 mm (52%), 71 mm (35%) and 50 mm (19%) in the first three years post‐harvest in treated catchment A and by 103 mm (118%), 157 mm (82%) and 119 mm (48%) in treated catchment B relative to the native forest control catchment. In the fourth post‐harvest water year annual rainfall was only 488 mm, which resulted in negligible run‐off in all catchments, regardless of forest cover. In both plantation catchments, monthly streamflows increased significantly (p = 0·01, p < 0·001) due to a significant increase in baseflows (p < 0·001) after harvesting. Monthly stormflows were not significantly affected by harvesting. Flow duration curve analyses indicated a variable response between the two plantation catchments. Treated catchment A was converted from an ephemeral stream flowing 42% of the time pre‐harvest to a temporary stream flowing 82% of the time post‐harvest. These changes occurred throughout all seasons of the year but were most pronounced during summer and autumn when baseflows were maintained post‐harvest but were not observed under native forest or mature pine plantations. By contrast, flow duration increased in treated catchment B from 12% of the time pre‐harvest to 38% of the time post‐harvest with the greatest changes measured during the winter and spring months when streamflow would normally occur under native forest conditions. These observations have important implications for the development of models of plantation water use to be utilized in water resource planning in Australia. Copyright © 2009 John Wiley & Sons, Ltd.     

Trend analysis of nutrient loadings in a large prairie catchment

Nutrient loadings in many river catchments continue to increase due to rapid expansion of agriculture, urban and industrial development, and population growth. Nutrient enrichment of water bodies has intensified eutrophication which degrades water quality and ecosystem health. In this study, we carried out a trend analysis of total phosphorus and total nitrogen loads in the South Saskatchewan River (SSR) catchment using a novel approach to analyse nutrient time series. Seasonal analysis of trends at each of the water quality stations was performed to determine the relationships between annual flow regimes and nutrient loads in the catchment, in particular, the influence of the high spring runoff on nutrient export. Decadal analysis was also performed to determine the long-term relationships of nutrients with anthropogenic changes in the catchment. Although it was found that seasonal and historical variability of nutrient load trends is mainly determined by streamflow regime changes, there is evidence that increases in nitrogen concentration can also be attributed to anthropogenic changes.     

Nutrient source attribution: Quantitative typology distinction of active and legacy source contributions to waterborne loads

Distinction between active and legacy sources of nutrients is needed for effective reduction of waterborne nutrient loads and associated eutrophication. This study quantifies main typological differences in nutrient load behaviour versus water discharge for active and legacy sources. This quantitative typology is used for source attribution based on monitoring data for water discharge and concentrations of total nitrogen (TN) and total phosphorous (TP) from 37 catchments draining into the Baltic Sea along the coastline of Sweden over the period 2003–2013. Results indicate dominant legacy source contributions to the monitored loads of TN and TP in most (33 of the total 37) study catchments. Dominant active sources are indicated in 1 catchment for TN, and mixed sources are indicated in 3 catchments for TN, and 4 catchments for TP. The TN and TP concentration contributions are quantified to be overall higher from the legacy than the active sources. Legacy concentrations also correlate well with key indicators of human activity in the catchments, agricultural land share for TN (R² = 0.65) and population density for TP (R² = 0.56). Legacy-dominated nutrient concentrations also change more slowly than in catchments with dominant active or mixed sources. Various data-based results and indications converge in indicating legacy source contributions as largely dominant, mainly anthropogenic, and with near-zero average change trends in the present study of catchments draining into the Baltic Sea along the coastline of Sweden, as in other parts of the world. These convergent indications emphasize needs to identify and map the different types of sources in each catchment, and differentiate strategies and measures to target each source type for possible achievement of shorter- and longer-term goals of water quality improvement.