The response of flow duration curves to afforestation

The hydrologic effect of replacing pasture or other short crops with trees is reasonably well understood on a mean annual basis. The impact on flow regime, as described by the annual flow duration curve (FDC) is less certain. A method to assess the impact of plantation establishment on FDCs was developed. The starting point for the analyses was the assumption that rainfall and vegetation age are the principal drivers of evapotranspiration. A key objective was to remove the variability in the rainfall signal, leaving changes in streamflow solely attributable to the evapotranspiration of the plantation. A method was developed to (1) fit a model to the observed annual time series of FDC percentiles; i.e. 10th percentile for each year of record with annual rainfall and plantation age as parameters, (2) replace the annual rainfall variation with the long term mean to obtain climate adjusted FDCs, and (3) quantify changes in FDC percentiles as plantations age. Data from 10 catchments from Australia, South Africa and New Zealand were used. The model was able to represent flow variation for the majority of percentiles at eight of the 10 catchments, particularly for the 10–50th percentiles. The adjusted FDCs revealed variable patterns in flow reductions with two types of responses (groups) being identified. Group 1 catchments show a substantial increase in the number of zero flow days, with low flows being more affected than high flows. Group 2 catchments show a more uniform reduction in flows across all percentiles. The differences may be partly explained by storage characteristics. The modelled flow reductions were in accord with published results of paired catchment experiments. An additional analysis was performed to characterise the impact of afforestation on the number of zero flow days (N_zero) for the catchments in group 1. This model performed particularly well, and when adjusted for climate, indicated a significant increase in N_zero. The zero flow day method could be used to determine change in the occurrence of any given flow in response to afforestation. The methods used in this study proved satisfactory in removing the rainfall variability, and have added useful insight into the hydrologic impacts of plantation establishment. This approach provides a methodology for understanding catchment response to afforestation, where paired catchment data is not available.     

Use of Gene Expression Programming in regionalization of flow duration curve

In this paper, a recently introduced artificial intelligence technique known as Gene Expression Programming (GEP) has been employed to perform symbolic regression for developing a parametric scheme of flow duration curve (FDC) regionalization, to relate selected FDC characteristics to catchment characteristics. Stream flow records of selected catchments located in the Auckland Region of New Zealand were used. FDCs of the selected catchments were normalised by dividing the ordinates by their median value. Input for the symbolic regression analysis using GEP was (a) selected characteristics of normalised FDCs; and (b) 26 catchment characteristics related to climate, morphology, soil properties and land cover properties obtained using the observed data and GIS analysis. Our study showed that application of this artificial intelligence technique expedites the selection of a set of the most relevant independent variables out of a large set, because these are automatically selected through the GEP process. Values of the FDC characteristics obtained from the developed relationships have high correlations with the observed values.     

A new method for modelling flow duration curves and predicting streamflow regimes under altered land-use conditions / Une nouvelle méthode de modélisation des courbes de débits classés et de prévision des régimes d'écoulement sous conditions modifiées d'occupation du sol

Abstract

The impact of changes in land-use/cover on streamflow at short time scales is evaluated by examining the changes in the flow duration curve (FDC) before and after land-use/cover change. The FDC characterizes the relationship between the magnitude and frequency and hence provides the complete range of streamflow over time. Two issues need to be considered in predicting the FDC due to land-use changes: (a) the appropriate parametric form of the FDC that enables application of the same expression of FDCs before and after the changes; and (b) the ability of parameters to capture and characterize the effect of land-use changes. In this paper, we propose a model which is a four-parameter double power form as a function of the FDC, where the two hydrological parameters represent the mean annual flow ([Qbar]) and the cease-to-flow point (τ expressed as a percentage), while the other two parameters (α and β) determine the shape of the FDC. The properties of this function are investigated in order to assign parameters to cope with the land-use changes. The model is used for several typical catchments in Australia for demonstration.     

Prediction of streamflow regimes over large geographical areas: interpolated flow–duration curves for the Danube region

Flow–duration curves (FDCs) are essential to support decisions on water resources management, and their regionalization is fundamental for the assessment of ungauged basins. In comparison with calibrated rainfall–runoff models, statistical methods provide data-driven estimates representing a useful benchmark. The objective of this work is the interpolation of FDCs from ~500 discharge gauging stations in the Danube. To this aim we use total negative deviation top-kriging (TNDTK), as multi-regression models are shown to be unsuitable for representing FDCs across all durations and sites. TNDTK shows a high accuracy for the entire Danube region, with overall Nash-Sutcliffe efficiency values computed in a leave-p-out cross-validation scheme (p equal to one site, one-third and half of the sites), all above 0.88. A reliability measure based on kriging variance is attached to each interpolated FDC at ~4000 prediction nodes. The GIS layer of regionalized FDCs is made available for broader use in the region.     

Prediction of streamflow regimes over large geographical areas: interpolated flow–duration curves for the Danube region

ABSTRACT

Flow–duration curves (FDCs) are essential to support decisions on water resources management, and their regionalization is fundamental for the assessment of ungauged basins. In comparison with calibrated rainfall–runoff models, statistical methods provide data-driven estimates representing a useful benchmark. The objective of this work is the interpolation of FDCs from ~500 discharge gauging stations in the Danube. To this aim we use total negative deviation top-kriging (TNDTK), as multi-regression models are shown to be unsuitable for representing FDCs across all durations and sites. TNDTK shows a high accuracy for the entire Danube region, with overall Nash-Sutcliffe efficiency values computed in a leave-p-out cross-validation scheme (p equal to one site, one-third and half of the sites), all above 0.88. A reliability measure based on kriging variance is attached to each interpolated FDC at ~4000 prediction nodes. The GIS layer of regionalized FDCs is made available for broader use in the region.     

Combining a rainfall–runoff model and a regionalization approach for flood and water resource assessment in the western Po Valley,Italy

ABSTRACT

Selecting the best structure and parameterization of rainfall–runoff models is not straightforward and depends on a broad number of factors. In this study, the “Modello Idrologico Semi-Distribuito in continuo” (MISDc) was tested on 63 mountainous catchments in the western Po Valley (Italy) and the optimal model parameters were regionalized using different strategies. The model performance was evaluated through several indexes analysing hydrological regime, high-flow condition and flow–duration curve (FDC). In general, MISDc provides a good fit behaviour with a Kling-Gupta Efficiency index greater than 0.5 for 100% and 84% of cases for calibration and validation, respectively. Concerning the regionalization, spatial proximity approach is the most accurate solution obtaining satisfactory performance. Lastly, the predicted FDCs showed an excellent similarity with the observed ones. Results encourage to apply MISDc over the study area for flood forecasting and for assessing water resources availability thanks to the modest computational efforts and data requirements.     

THE EFFECTS OF INACCURATE INPUT PARAMETERS ON DEPOSITION OF SUSPENDED SEDIMENT

I.~OntjcnONEngineeringPredictionsofsedimenttransPOrtinopen-Channelsystems,especiallyifunsteadyflowconditionmustbeconsidered,isfraughtwithambigUityduetospatialandtemporalvariability,measurementerrors,lumtedsamplingoftheParameters,bounceandscallconditionS,andsink/sourcetenns.Duringrecentyearsmanysophisticatedmathematicalmodelshavebecomeavailabletocalculateunsteadyflowinone,tWoorthreedimensionS.ThesemodelswillbecoupledwithtranSportequationSforbedloadorsuSpendedloadtocalculateforexamplelongt…     

Regional models of flow-duration curves of perennial and intermittent streams and their use for calibrating the parameters of a rainfall–runoff model

Abstract

One of the main challenges faced by hydrologists and water engineers is the estimation of variables needed for water resources planning and management in ungauged river basins. To this end, techniques for transposing information, such as hydrological regional analyses, are widely employed. A method is presented for regionalizing flow-duration curves (FDCs) in perennial, intermittent and ephemeral rivers, based on the extended Burr XII probability distribution. This distribution shows great flexibility to fit data, with accurate reproduction of flow extremes. The performance analysis showed that, in general, the regional models are able to synthesize FDCs in ungauged basins, with a few possible drawbacks in the application of the method to intermittent and ephemeral rivers. In addition to the regional models, we summarize the experience of using synthetic FDCs for the indirect calibration of the Rio Grande rainfall–runoff model parameters in ungauged basins.

Editor D. Koutsoyiannis

Citation Costa, V., Fernandes, W., and Naghettini, M., 2013. Regional models of flow-duration curves of perennial and intermittent streams and their use for calibrating the parameters of a rainfall–runoff model. Hydrological Sciences Journal, 59 (2), 262–277.     

River water quality response under hypothetical climate change scenarios in Tunga‐Bhadra river,India

Analysis of climate change impacts on streamflow by perturbing the climate inputs has been a concern for many authors in the past few years, but there are few analyses for the impacts on water quality. To examine the impact of change in climate variables on the water quality parameters, the water quality input variables have to be perturbed. The primary input variables that can be considered for such an analysis are streamflow and water temperature, which are affected by changes in precipitation and air temperature, respectively. Using hypothetical scenarios to represent both greenhouse warming and streamflow changes, the sensitivity of the water quality parameters has been evaluated under conditions of altered river flow and river temperature in this article. Historical data analysis of hydroclimatic variables is carried out, which includes flow duration exceedance percentage (e.g. Q90), single low‐flow indices (e.g. 7Q10, 30Q10) and relationships between climatic variables and surface variables. For the study region of Tunga‐Bhadra river in India, low flows are found to be decreasing and water temperatures are found to be increasing. As a result, there is a reduction in dissolved oxygen (DO) levels found in recent years. Water quality responses of six hypothetical climate change scenarios were simulated by the water quality model, QUAL2K. A simple linear regression relation between air and water temperature is used to generate the scenarios for river water temperature. The results suggest that all the hypothetical climate change scenarios would cause impairment in water quality. It was found that there is a significant decrease in DO levels due to the impact of climate change on temperature and flows, even when the discharges were at safe permissible levels set by pollution control agencies (PCAs). The necessity to improve the standards of PCA and develop adaptation policies for the dischargers to account for climate change is examined through a fuzzy waste load allocation model developed earlier. Copyright © 2011 John Wiley & Sons, Ltd.     

Rating curve uncertainty and change detection in discharge time series: case study with 44‐year historic data from the Nyangores River,Kenya

The intersection of the developing topic of rating curve and discharge series uncertainty with the topic of hydrological change detection (e.g., in response to land cover or climatic change) has not yet been well studied. The work herein explores this intersection, with consideration of a long‐term discharge response (1964–2007) for a ~650‐km² headwater basin of the Mara River in west Kenya, starting with stream rating and daily gauge height data. A rating model was calibrated using Bayesian methods to quantify uncertainty intervals in model parameters and predictions. There was an unknown balance of random and systemic error in rating data scatter (a scenario not likely unique to this basin), which led to an unknown balance of noise and information in the calibrated statistical error model. This had implications on testing for hydrological change. Overall, indications were that shifts in basin's discharge response were rather subtle over the 44‐year period. A null hypothesis for change using flow duration curves (FDCs) from four different 8‐year data intervals could be either accepted or rejected over much of the net flow domain depending on different applications of the statistical error model (each with precedence in the literature). The only unambiguous indication of change in FDC comparisons appeared to be a reduction in lowest baseflow in recent years (flows with >98% exceedance probability). We defined a subjective uncertainty interval based on an intermediate balance of random and systematic error in the rating model that suggested a possibility of more prevalent impacts. These results have relevance to management in the Mara basin and to future studies that might establish linkages to historic land use and climatic factors. The concern about uncertain uncertainty intervals (uncertainty²) extends beyond the Mara and is relevant to testing change where non‐random rating errors may be important and subtle responses are investigated. Copyright © 2013 John Wiley & Sons, Ltd.     

A Method to Represent a Well in a Three-Dimensional Discrete Fracture Network Model

In discrete fracture network (DFN) modeling, fractures are randomly generated and placed in the model domain. The rock matrix is considered impermeable. Small fractures and isolated fractures are often ignored to reduce computational expense. As a result, the rock matrix between fractures could be large and intersections may not be found between a well introduced in the model and the hydraulically connected fracture networks (fracture backbones). To overcome this issue, this study developed a method to conceptualize a well in a three-dimensional (3D) DFN using two orthogonal rectangular fractures oriented along the well's axis. Six parameters were introduced to parameterize the well screen and skin zone, and to control the connectivity between the well and the fracture backbones. The two orthogonal fractures were discretized using a high-resolution mesh to improve the quality of flow and transport simulations around and along the well. The method was successfully implemented within dfnWorks 2.0 (Hyman et al. 2015) to incorporate a well in a 3D DFN and to track particles leaving an injection well and migrating to a pumping well. Verification of the method against MODFLOW/MODPATH found a perfect match in simulated hydraulic head and particle tracking. Using three examples, the study showed that the method ensured the connectivity between wells and fracture backbones, and honored the physical processes of flow and transport along and around wells in DFNs. Recommendations are given for estimating the values of the six introduced well parameters in a real-world case study.     

Using synthetic flow duration curves for rainfall–runoff model calibration at ungauged sites

The primary objective of the study is to propose a strategy for rainfall–runoff model calibration at ungauged sites. This strategy comprises two main components: (1) development of the regional analysis method to synthesize the flow duration curves at ungauged sites; and (2) utilization of the synthetic flow duration curves for model calibration. Since the regional analysis method can synthesize the flow duration curves at ungauged sites, the continuous rainfall–runoff model coupled with a global optimization method were applied in southern Taiwan using the synthetic flow duration curve as an objective for model calibration. The results reveal that the regional flow duration curve and the strategy for model calibration at ungauged sites have good performances in the study area. Copyright © 2000 John Wiley & Sons, Ltd.     

Low and high flow analyses and wavelet application for characterization of the Blue Nile River system

The low and high flow characteristic of the Blue Nile River (BNR) basin is presented. The study discusses low and high flow, flow duration curve (FDC) and trend analysis of the BNR and its major tributaries. Different probability density functions were fitted to better describe the low and high flows of the BNR and major tributaries in the basin. Wavelet analysis was used in understanding the variance and frequency‐time localization and detection of dominant oscillations in rainfall and flow. FDCs were developed, and low flow (below 50% exceedance) and high flow (over 75% exceedance) of the curves were analysed and compared. The Gravity Recovery and Climate Experiment (GRACE) satellite‐based maps of monthly changes in gravity converted to water equivalents from 2003 to 2006 for February, May and September showed an increase in the moisture influx in the BNR basin for the month of September, and loss of moisture in February and May. It was also shown that 2004 and 2005 were drier with less moisture influx compared to 2003 and 2006. On the basis of the Kolmogorov‐Smirnov, Anderson‐Darling and Chi‐square tests, Gen. Pareto, Frechet 3P, Log‐normal, Log‐logistics, Fatigue Life and Phased Bi‐Weibull distributions best describe the low and high flows within the BNR basin. This will be beneficial in developing flow hydrographs for similar ungauged watersheds within the BNR basin. The below 50% and above 75% exceedance on the FDC for five major rivers in addition to the BNR showed different characteristics depending on size, land cover, topography and other factors. The low flow frequency analysis of the BNR at Bahir Dar showed 0·55 m³/s as the monthly low flow with recurrence interval of 10 years. The wavelet analysis of the rainfall (at Bahir Dar and basin‐wide) and flows at three selected stations shows inter‐ and intra‐annual variability of rainfall and flows at various scales. Copyright © 2009 John Wiley & Sons, Ltd.     

Analysis of nonstationarity in low flow in the Loess Plateau of China

An approach for nonstationary low‐flow frequency analysis is developed and demonstrated on a dataset from the rivers on the Loess Plateau of China. Nonstationary low‐flow frequency analysis has drawn significant attention in recent years by establishing relationships between low‐flow series and explanatory variables series, but few studies have tested whether the time‐varying moments of low flow can be fully described by the time‐varying moments of the explanatory variables. In this research, the low‐flow distributions are analytically derived from the 2 basic explanatory variables—the recession duration and the recession coefficient—with the assumption that the recession duration and recession coefficient variables follow exponential and gamma distributions, respectively; the derived low‐flow distributions are applied to test whether the time‐varying moments of explanatory variables can explain the nonstationarities found in the low‐flow variable. The effects of ecosystem construction measures, that is, check dam, terrace, forest, and grassland, on the recession duration and recession coefficient are further discussed. Daily flow series from 11 hydrological stations from the Loess Plateau are used and processed with a moving average technique. Low‐flow data are extracted following the pit under threshold approach. Six of the 11 low‐flow series show significant nonstationarities at the 5% significance level, and the trend curves of the moments of low flow are in close agreement with the curves estimated from the derived distribution with time‐dependent moments of the recession duration and time‐constant moments of the recession coefficient. It is indicated that the nonstationarity in the low‐flow distribution results from the nonstationarity in the recession duration in all 6 cases, and the increase in the recession duration is resulted from large‐scale ecosystem constructions rather than climate change. The large‐scale ecosystem constructions are found to have more influence on the decrease in streamflow than on the increase in watershed storage, thus resulting in the reduction of low flow. A high return period for the initial fixed design value decreases dramatically with an increasing recession duration.     

Synthetic regional flow duration curve for southern Taiwan

The primary purpose of this study is to develop the regional flow duration curves for southern Taiwan. To define homogeneous regions for developing regional flow duration curves, multivariate statistical analysis (principal component and cluster analysis) was applied to daily flow data from 34 stream-gauged stations in southern Taiwan. Two kinds of clustering variables, the dimensionless flow duration curve and specific flow duration curve, were compared in this study. It was found that three homogeneous regions delineated by specific flow duration curves as clustering variables have more reasonable results. The three homogeneous regions not only have well-defined geographical boundaries, but also correspond to the rainfall and geology characteristics of the regions. It seems that the technique of cluster analysis can reasonably define the homogeneous regions. In each homogeneous region, the synthetic regional flow duration curves were developed by a family of parametric duration curves. This approach has the advantage of being simple and needing only the basin area as an index. The performance of the regional flow duration curve was verified by the comparison of areas under the actual and synthetic flow duration curves; the latter were generated from the regional flow duration curve. Almost all the 34 stream-gauged stations had less than 25% absolute error.     

A statistical analysis of flank eruptions on Etna volcano

A singularly complete record exists for the eruptive activity of Etna volcano. The time series of occurrence of flank eruptions in the period 1600–1980, in which the record is presumably complete, is found to follow a stationary Poisson process. A revision of the available data shows that eruption durations are rather well correlated with the estimates of the volume of lava flows. This implies that the magnitude of an eruption can be defined directly by its duration. Extreme value statistics are then applied to the time series, using duration as a dependent variable. The probability of occurrence of a very long (300 days) eruption is greater than 50% only in time intervals of the order of 50 years. The correlation found between duration and total output also allows estimation of the probability of occurrence of a major event which exceeds a given duration and total flow of lava. The composite probabilities do not differ considerably from the pure ones. Paralleling a well established application to seismic events, extreme value theory can be profitably used in volcanic risk estimates, provided that appropriate account is also taken of all other variables.     

An approach to measure parameter sensitivity in watershed hydrological modelling

Hydrological responses vary spatially and temporally according to watershed characteristics. In this study, the hydrological models that we developed earlier for the Little Miami River (LMR) and Las Vegas Wash (LVW) watersheds in the USA were used for detailed sensitivity analyses. To compare the relative sensitivities of the hydrological parameters of these two models, we used normalized root mean square error (NRMSE). By combining the NRMSE index with the flow duration curve analysis, we derived an approach to measure parameter sensitivities under different flow regimes. Results show that the parameters related to groundwater are highly sensitive in the LMR watershed, whereas the LVW watershed is primarily sensitive to near-surface and impervious parameters. The high and medium flows are more impacted by most of the parameters. The low flow regime was highly sensitive to groundwater-related parameters. Moreover, our approach is found to be useful in facilitating model development and calibration.

EDITOR D. Koutsoyiannis

ASSOCIATE EDITOR S. Huang     

Stochastic analysis of flood series

A stochastic model based on the renewal process was developed and used to analyse the characteristics of floods: the volume exceedence, the duration of the flood and the maximum annual flow. The model contains a method for determination of total annual volume exceedence and total annual duration of floods, as well as a method for calculation of maximum annual exceedence, maximum flood duration and maximum flow. The subset of the flood occurrence number in a given time interval is common for all analysed phenomena (volume exceedence, flood duration, maximum flow). The subset of given exceedences is common for total annual volume exceedence, as well as for maximum annual volume exceedence. The same holds for durations of individual floods. The model was then applied to analyse the floods on the Drina River at the Paunci hydrological station and on the Danube River at the Bezdan station.     

Changes in streamflow regime following vegetation changes from paired catchments

Vegetation changes can significantly affect catchment water balance. It is important to evaluate the effects of vegetation cover change on streamflow as changes in streamflow relate to water security. This study focuses on the use of statistical methods to determine responses in streamflow at seven paired catchments in Australia, New Zealand, and South Africa to vegetation change. The non‐parametric Mann–Kendall test and Pettitt's test were used to identify trends and change points in the annual streamflow records. Statistically significant trends in annual streamflow were detected for most of the treated catchments. It took between 3 and 10 years for a change in vegetation cover to result in significant change in annual streamflow. Presence of the change points in streamflow was associated with changes in the mean, variance, and distribution of annual streamflow. The streamflow in the deforestation catchments increased after the change points, whereas reduction in streamflow was observed in the afforestation catchments. The streamflow response is mainly affected by the climate and underlying vegetation change. Daily flow duration curves (FDCs) for the whole period and pre‐change and post‐change point periods also were analysed to investigate the changes in flow regime. Three types of vegetation change effects on the flow regime have been identified. The relative reductions in most percentile flows are constant in the afforestation catchments. The comparison of trend, change point, and FDC in the annual streamflow from the paired experiments reflects the important role of the vegetation change. Copyright © 2011 John Wiley & Sons, Ltd.