Modelling surface runoff to evaluate the effects of wildfires in multiple semi‐arid,shrubland‐dominated catchments

Wildfires change the infiltration properties of soil, reduce the amount of interception and result in increased runoff. A wildfire at Northeast Attica, Central Greece, in August 2009, destroyed approximately one third of a study area consisting of a mixture of shrublands, pastures and pines. The present study simultaneously models multiple semi‐arid, shrubland‐dominated Mediterranean catchments and assesses the hydrological response (mean annual and monthly runoff and runoff coefficients) during the first few years following wildfires. A physically based, hydrological model (MIKE SHE) was chosen. Calibration and validation results of mean monthly discharge presented very good agreement with the observed data for the pre‐wildfire and post‐wildfire period for two subcatchments (Nash–Sutcliffe Efficiency coefficient of 79.7%). The model was then used to assess the pre‐wildfire and post‐wildfire runoff responses for each of seven catchments in the study area. Mean annual surface runoff increased for the first year and after the second year following the wildfires increased by 112% and 166%, respectively. These values are within the range observed in similar cases of monitored sites. This modelling approach may provide a way of prioritizing catchment selection with respect to post‐fire remediation activities. Additionally, this modelling assessment methodology would be valuable to other semi‐arid areas because it provides an important means for comprehensively assessing post‐wildfire response over large regions and therefore attempts to address some of the scaled issues in the specific literature field of research. Copyright © 2015 John Wiley & Sons, Ltd.     

The importance of surface controls on overland flow connectivity in semi‐arid environments: results from a numerical experimental approach

In semi‐arid environments, the characteristics of the land surface determine how rainfall is transformed into surface runoff and influences how this runoff moves from the hillslopes into river channels. Whether or not water reaches the river channel is determined by the hydrological connectivity. This paper uses a numerical experiment‐based approach to systematically assess the effects of slope length, gradient, flow path convergence, infiltration rates and vegetation patterns on the generation and connectivity of runoff. The experiments were performed with the Connectivity of Runoff Model, 2D version distributed, physically based, hydrological model. The experiments presented are set within a semi‐arid environment, characteristic of south‐eastern Spain, which is subject to low frequency high rainfall intensity storm events. As a result, the dominant hydrological processes are infiltration excess runoff generation and surface flow dynamics. The results from the modelling experiments demonstrate that three surface factors are important in determining the form of the discharge hydrograph: the slope length, the slope gradient and the infiltration characteristics at the hillslope‐channel connection. These factors are all related to the time required for generated runoff to reach an efficient flow channel, because once in this channel, the transmission losses significantly decrease. Because these factors are distributed across the landscape, they have a fundamental role in controlling the landscape hydrological response to storm events. Copyright © 2013 John Wiley & Sons, Ltd.     

Erosion and terrace failure due to agricultural land abandonment in a semi‐arid environment

Agricultural land abandonment is currently widely spread in Mediterranean countries and a further increase is expected. Previous research has shown that abandoned fields in semi‐arid areas are more vulnerable to gully erosion. The absence of ploughing and slow vegetation recovery cause the formation of soil crusts with low infiltration rates, resulting in increased runoff and gully erosion risk. The objective of our study was to assess the extent and causes of erosion and terrace failure on abandoned fields and to discuss options for mitigation. The study was carried out in the Carcavo basin, a semi‐arid catchment in southeast Spain. At catchment scale all abandoned fields were surveyed and characteristics of each field were described. Additionally we surveyed abandoned and cultivated terraces and used statistical analyses to determine the factors that induce terrace failure. At field scale we constructed a detailed digital elevation model (DEM) for an abandoned terrace field in order to calculate sediment losses since time of abandonment. The results revealed that more than half the abandoned fields had moderate to severe erosion and the statistical analysis showed that these fields had significantly steeper slopes, were terraced and had cereals as previous land use. Factors that increase the risk of terrace failure were land abandonment, steeper terrace slope, loam texture, valley‐bottom position and shrubs on the terrace wall. The reconstructed erosion rate (87 ton ha^?1 year^?1) confirmed the importance of gully erosion on these abandoned terrace fields. Potential soil and water conservation practices to mitigate soil erosion after abandonment are: (1) maintenance of terrace walls, as a result more water is retained, which increases vegetation cover and consequently decreases erosion. (2) Revegetation with indigenous grass species on spots with concentrated flow, especially near terrace walls. Copyright © 2008 John Wiley & Sons, Ltd.     

Using high resolution tracer data to constrain water storage,flux and age estimates in a spatially distributed rainfall‐runoff model

Models simulating stream flow and conservative tracers can provide a representation of flow paths, storage distributions and mixing processes that is advantageous for many predictive purposes. Compared with models that only simulate stream flow, tracer data can be used to investigate the internal consistency of model behaviour and to gain insight into model performance. Here, we examine the strengths and weaknesses of a data‐driven, spatially distributed tracer‐aided rainfall‐runoff model. The model structure allowed us to assess the influence of landscape characteristics on the routing and mixing of water and tracers. The model was applied to a site in the Scottish Highlands with a unique tracer data set; ~4 years of daily isotope ratios in stream water and precipitation were available, as well as 2 years of weekly soil and ground water isotopes. The model structure was based on an empirically based, lumped tracer‐aided model previously developed for the catchment. The best model runs were selected from Monte Carlo simulations based on dual calibration criteria using objective functions for both stream isotopes and discharge at the outlet. Model performance for these criteria was reasonable (Nash–Sutcliffe efficiencies for discharge and isotope ratios were ~0.4–0.6). The model could generally reproduce the variable isotope signals in the soils of the steeper hill slopes where storage was low, and damped isotope responses in valley bottom cells with high storage. The model also allowed us to estimate the age distributions of internal stores, water fluxes and stream flow. Average stream water age was ~1.6 years, integrating older groundwater in the valley bottom and dynamic younger soil waters. By tracking water ages and simulating isotopes, the model captured the changes in connectivity driven by distributed storage dynamics. This has substantially improved the representation of spatio‐temporal process dynamics and gives a more robust framework for projecting environmental change impacts. Copyright © 2016 The Authors Hydrological Processes Published by John Wiley & Sons Ltd.     

Hydrological response of a Brazilian semi‐arid catchment to different land use and climate change scenarios: a modelling study

Brazilian semi‐arid regions are characterized by water scarcity, vulnerability to desertification, and climate variability. The investigation of hydrological processes in this region is of major interest not only for water planning strategies but also to address the possible impact of future climate and land‐use changes on water resources. A hydrological distributed catchment‐scale model (DiCaSM) has been applied to simulate hydrological processes in a small representative catchment of the Brazilian northeast semi‐arid region, and also to investigate the impact of climate and land‐use changes, as well as changes associated with biofuel/energy crops production. The catchment is part of the Brazilian network for semi‐arid hydrology, established by the Brazilian Federal Government. Estimating and modelling streamflow (STF) and recharge in semi‐arid areas is a challenging task, mainly because of limitation in in situ measurements, and also due to the local nature of some processes. Direct recharge measurements are very difficult in semi‐arid catchments and contain a high level of uncertainty. The latter is usually addressed by short‐ and long‐time‐scale calibration and validation at catchment scale, as well as by examining the model sensitivity to the physical parameters responsible for the recharge. The DiCaSM model was run from 2000 to 2008, and streamflow was successfully simulated, with a Nash–Sutcliffe (NS) efficiency coefficient of 0·73, and R² of 0·79. On the basis of a range of climate change scenarios for the region, the DiCaSM model forecasted a reduction by 35%, 68%, and 77%, in groundwater recharge (GWR), and by 34%, 65%, and 72%, in streamflow, for the time spans 2010–2039, 2040–2069, and 2070–2099, respectively, could take place for a dry future climate scenario. These reductions would produce severe impact on water availability in the region. Introducing castor beans to the catchment would increase the GWR and streamflow, mainly if the caatinga areas would be converted into castor beans production. Changing an area of 1000 ha from caatinga to castor beans would increase the GWR by 46% and streamflow by 3%. If the same area of pasture is converted into castor beans, there would be an increase in GWR and streamflow by 24% and 5%, respectively. Such results are expected to contribute towards environmental policies for north‐east Brazil (NEB), and to biofuel production perspectives in the region. Copyright © 2010 John Wiley & Sons, Ltd.     

Simulating saturation‐excess surface run‐off in a semi‐distributed hydrological model

The level of complexity, and the number of parameters, to include in a hydrological model is a relatively contentious issue in hydrological modelling. However, it can be argued that explicitly representing important run‐off generation processes can improve the practical value of a model's outputs. This paper explores the benefits of including a new function into an existing semi‐distributed hydrological model (the Pitman model) that is widely used in the sub‐Saharan Africa region. The new function was designed to represent saturation‐excess surface run‐off processes at subcatchment scales and was motivated by the evidence of dambo (low topography riparian areas) type features in many sub‐Saharan river basins. The results for uncertainty versions of the model, with and without the new function, were compared for 25 catchments, which were divided up into those where evidence of dambos exists and those where there is no such evidence. The results suggest that the new function certainly improves the model results for the catchments where dambos exist, but not in situations where saturation‐excess surface run‐off is not expected to occur. The overall conclusion is therefore that the addition of the new function is justified.     

Process‐based hydrological modelling: The potential of a bottom‐up approach for runoff predictions in ungauged catchments

Conceptual rainfall–runoff models are a valuable tool for predictions in ungauged catchments. However, most of them rely on calibration to determine parameter values. Improving the representation of runoff processes in models is an attractive alternative to calibration. Such an approach requires a straightforward, a priori parameter allocation procedure applicable on a wide range of spatial scales. However, such a procedure has not been developed yet. In this paper, we introduce a process‐based runoff generation module (RGM‐PRO) as a spin‐off of the traditional runoff generation module of the PREVAH hydrological modelling system. RGM‐PRO is able to exploit information from maps of runoff types, which are developed on the basis of field investigations and expert knowledge. It is grid based, and within each grid cell, the process heterogeneity is considered to avoid information loss due to grid resolution. The new module is event based, and initial conditions are assimilated and downscaled from continuous simulations of PREVAH, which are also available for real‐time applications. Four parameter allocation strategies were developed, on the basis of the results of sprinkling experiments on 60‐m² hillslope plots at several grassland locations in Switzerland, and were tested on five catchments on the Swiss Plateau and Prealps. For the same catchments, simulation results obtained with the best parameter allocation strategy were compared with those obtained with different configurations of the traditional runoff generation module of PREVAH, which was also applied as an event‐based module here. These configurations include a version that avoids calibration, one that transfers calibrated parameters, and one that uses regionalised parameter values. RGM‐PRO simulated heavy events in a more realistic way than the uncalibrated traditional runoff generation module of PREVAH, and, in some instances, it even exceeded the performance of the calibrated traditional one. The use of information on the spatial distribution of runoff types additionally proved to be valuable as a regionalisation technique and showed advantages over the other regionalisation approaches, also in terms of robustness and transferability.     

Changes in ecosystem structure,function and hydrological connectivity control water,soil and carbon losses in semi‐arid grass to woody vegetation transitions

Connectivity has recently emerged as a key concept for understanding hydrological response to vegetation change in semi‐arid environments, providing an explanatory link between abiotic and biotic, structure and function. Reduced vegetation cover following woody encroachment, generally promotes longer, more connected overland flow pathways, which has the potential to result in an accentuated rainfall‐runoff response and fluxes of both soil erosion and carbon. This paper investigates changing hydrological connectivity as an emergent property of changing ecosystem structure over two contrasting semi‐arid grass to woody vegetation transitions in New Mexico, USA. Vegetation structure is quantified to evaluate if it can be used to explain observed variations in water, sediment and carbon fluxes. Hydrological connectivity is quantified using a flow length metric, combining topographic and vegetation cover data. Results demonstrate that the two woody‐dominated sites have significantly longer mean flowpath lengths (4 · 3 m), than the grass‐dominated sites (2 · 4 m). Mean flowpath lengths illustrate a significant positive relationship with the functional response. The woody‐dominated sites lost more water, soil and carbon than their grassland counterparts. Woody sites erode more, with mean event‐based sediment yields of 1203 g, compared to 295 g from grasslands. In addition, the woody sites lost more organic carbon, with mean event yields of 39 g compared to 5 g from grassland sites. Finally, hydrological connectivity (expressed as mean flowpath length) is discussed as a meaningful measure of the interaction between structure and function and how this manifests under the extreme rainfall that occurs in semi‐arid deserts. In combination with rainfall characteristics, connectivity emerges as a useful tool to explain the impact of vegetation change on water, soil and carbon losses across semi‐arid environments. Copyright © 2013 John Wiley & Sons, Ltd.     

The use of agent based modelling techniques in hydrology: determining the spatial and temporal origin of channel flow in semi‐arid catchments

Information on the spatial and temporal origin of runoff entering the channel during a storm event would be valuable in understanding the physical dynamics of catchment hydrology; this knowledge could be used to help design flood defences and diffuse pollution mitigation strategies. The majority of distributed hydrological models give information on the amount of flow leaving a catchment and the pattern of fluxes within the catchment. However, these models do not give any precise information on the origin of runoff within the catchment. The spatial and temporal distribution of runoff sources is particularly intricate in semi‐arid catchments, where there are complex interactions between runoff generation, transmission and re‐infiltration over short temporal scales. Agents are software components that are capable of moving through and responding to their local environment. In this application, the agents trace the path taken by water through the catchment. They have information on their local environment and on the basis of this information make decisions on where to move. Within a given model iteration, the agents are able to stay in the current cell, infiltrate into the soil or flow into a neighbouring cell. The information on the current state of the hydrological environment is provided by the environment generator. In this application, the Connectivity of Runoff Model (CRUM) has been used to generate the environment. CRUM is a physically based, distributed, dynamic hydrology model, which considers the hydrological processes relevant for a semi‐arid environment at the temporal scale of a single storm event. During the storm event, agents are introduced into the model across the catchment; they trace the flows of water and store information on the flow pathways. Therefore, this modelling approach is capable of giving a novel picture of the temporal and spatial dynamics of flow generation and transmission during a storm event. This is possible by extracting the pathways taken by the agents at different time slices during the storm. The agent based modelling approach has been applied to two small catchments in South East Spain. The modelling approach showed that the two catchments responded differently to the same rainfall event due to the differences in the runoff generation and overland flow connectivity between the two catchments. The model also showed that the time of travel to the nearest flow concentration is extremely important for determining the connectivity of a point in the landscape with the catchment outflow. Copyright © 2007 John Wiley & Sons, Ltd.     

Semi‐distributed flood runoff model at the subcontinental scale for southwestern Iran

Rainfall–runoff modelling, as a surface hydrological process, on large‐scale data‐poor basins is currently a major topic of investigation that requires the model parameters be identified by using basin physical characteristics rather than calibration. This paper describes the application of the TOPMODEL framework accompanied by a kinematic wave model to the Karun River sub‐basins in southwestern Iran with just one conceptual parameter for calibration. ISLSCP1, HYDRO1K and Reynolds data sets are presented in a geographical information system and used as data sources for meteorological information, hydrological features and soil characteristics of the study area respectively. The results show that although the model developed can adequately predict flood runoff in the catchment with only one calibrated parameter, it is suggested that the effect of surface reservoirs be considered in the proposed model. Copyright © 2007 John Wiley & Sons, Ltd.     

Use of the Connectivity of Runoff Model (CRUM) to investigate the influence of storm characteristics on runoff generation and connectivity in semi‐arid areas

Much attention has been given to the surface controls on the generation and transmission of runoff in semi‐arid areas. However, the surface controls form only one part of the system; hence, it is important to consider the effect that the characteristics of the storm event have on the generation of runoff and the transmission of flow across the slope. The impact of storm characteristics has been investigated using the Connectivity of Runoff Model (CRUM). This is a distributed, dynamic hydrology model that considers the hydrological processes relevant to semi‐arid environments at the temporal scale of a single storm event. The key storm characteristics that have been investigated are the storm duration, rainfall intensity, rainfall variability and temporal structure. This has been achieved through the use of a series of defined storm hydrographs and stochastic rainfall. Results show that the temporal fragmentation of high‐intensity rainfall is important for determining the travel distances of overland flow and, hence, the amount of runoff that leaves the slope as discharge. If the high‐intensity rainfall is fragmented, then the runoff infiltrates a short distance downslope. Longer periods of high‐intensity rainfall allow the runoff to travel further and, hence, become discharge. Therefore, storms with similar amounts of high‐intensity rainfall can produce very different amounts of discharge depending on the storm characteristics. The response of the hydrological system to changes in the rainfall characteristics can be explained using a four‐stage model of the runoff generation process. These stages are: (1) all water infiltrating, (2) the surface depression store filling or emptying without runoff occurring, (3) the generation and transmission of runoff and (4) the transmission of runoff without new runoff being generated. The storm event will move the system between the four stages and the nature of the rainfall required to move between the stages is determined by the surface characteristics. This research shows the importance of the variable‐intensity rainfall when modelling semi‐arid runoff generation. The amount of discharge may be greater or less than the amount that would have been produced if constant rainfall intensity is used in the model. Copyright © 2006 John Wiley & Sons, Ltd.     

Soil moisture states,lateral flow,and streamflow generation in a semi‐arid,snowmelt‐driven catchment

Hydraulic connectivity on hillslopes and the existence of preferred soil moisture states in a catchment have important controls on runoff generation. In this study we investigate the relationships between soil moisture patterns, lateral hillslope flow, and streamflow generation in a semi‐arid, snowmelt‐driven catchment. We identify five soil moisture conditions that occur during a year and present a conceptual model based on field studies and computer simulations of how streamflow is generated with respect to the soil moisture conditions. The five soil moisture conditions are (1) a summer dry period, (2) a transitional fall wetting period, (3) a winter wet, low‐flux period, (4) a spring wet, high‐flux period, and (5) a transitional late‐spring drying period. Transitions between the periods are driven by changes in the water balance between rain, snow, snowmelt and evapotranspiration. Low rates of water input to the soil during the winter allow dry soil regions to persist at the soil–bedrock interface, which act as barriers to lateral flow. Once the dry‐soil flow barriers are wetted, whole‐slope hydraulic connectivity is established, lateral flow can occur, and upland soils are in direct connection with the near‐stream soil moisture. This whole‐slope connectivity can alter near‐stream hydraulics and modify the delivery of water, pressure, and solutes to the stream. Copyright © 2005 John Wiley & Sons, Ltd.     

Hydrograph separation using tracers and digital filters to quantify runoff components in a semi‐arid mesoscale catchment

Chemical hydrograph separation using electrical conductivity and digital filters is applied to quantify runoff components in the 1,640 km² semi‐arid Kaap River catchment and its subcatchments in South Africa. A rich data set of weekly to monthly water quality data ranging from 1978 to 2012 (450 to 940 samples per site) was analysed at 4 sampling locations in the catchment. The data were routinely collected by South Africa's national Department of Water and Sanitation, using standard sampling procedures. Chemical hydrograph separation using electrical conductivity (EC) as a tracer was used as reference and a recursive digital filter was then calibrated for the catchment. Results of the two‐component hydrograph separation indicate the dominance of baseflow in the low flow regime, with a contribution of about 90% of total flow; however, during the wet season, baseflow accounts for 50% of total flow. The digital filter parameters were very sensitive and required calibration, using chemical hydrograph separation as a reference. Calibrated baseflow estimates ranged from 40% of total flow at the catchment outlet to 70% in the tributaries. The study demonstrates that routinely monitored water quality data, especially EC, can be used as a meaningful tracer, which could also aid in the calibration of a digital filter method and reduce uncertainty of estimated flow components. This information enhances our understanding of how baseflow is generated and contributed to streamflow throughout the year, which can aid in quantification of environmental flows, as well as to better parameterize hydrological models used for water resources planning and management. Baseflow estimates can also be useful for groundwater and water quality management.     

Hydrological processes of a closed catchment‐lake system in a semi‐arid Mediterranean environment

Data collected in 4 years of field observations were used in conjunction with continuous simulation models to study, at the small‐basin scale, the water balance of a closed catchment‐lake system in a semi‐arid Mediterranean environment. The open water evaporation was computed with the Penman equation, using the data set collected in the middle of the lake. The surface runoff was partly measured at the main tributary and partly simulated using a distributed, catchment, hydrological model, calibrated with the observed discharge. The simplified structure of the developed modelling mainly concerns soil moisture dynamics and bedrock hydraulics, whereas the flow components are physically based. The calibration produced high efficiency coefficients and showed that surface runoff is greatly affected by soil water percolation into fractured bedrock. The bedrock reduces the storm‐flow peaks and the interflow and has important multi‐year effects on the annual runoff coefficients. The net subsurface outflow from the lake was calculated as the residual of the lake water balance. It was almost constant in the dry seasons and increased in the wet seasons, because of the moistening of the unsaturated soil. During the years of observation, rainfall 30% higher than average caused abundant runoff and a continuous rise in the lake water levels. The analysis allows to predict that, in years with lower than the average rainfall, runoff will be drastically reduced and will not be able to compensate for negative balance between precipitation and lake evaporation. Such highly unsteady situations, with great fluctuations in lake levels, are typical of closed catchment‐lake systems in the semi‐arid Mediterranean environment. Copyright © 2012 John Wiley & Sons, Ltd.     

Assessment of rainfall‐runoff models based upon wavelet analysis

A basic hypothesis is proposed: given that wavelet‐based analysis has been used to interpret runoff time‐series, it may be extended to evaluation of rainfall‐runoff model results. Conventional objective functions make certain assumptions about the data series to which they are applied (e.g. uncorrelated error, homoscedasticity). The difficulty that objective functions have in distinguishing between different realizations of the same model, or different models of the same system, is that they may have contributed in part to the occurrence of model equifinality. Of particular concern is the fact that the error present in a rainfall‐runoff model may be time dependent, requiring some form of time localization in both identification of error and derivation of global objective functions. We explore the use of a complex Gaussian (order 2) wavelet to describe: (1) a measured hydrograph; (2) the same hydrograph with different simulated errors introduced; and (3) model predictions of the same hydrograph based upon a modified form of TOPMODEL. The analysis of results was based upon: (a) differences in wavelet power (the wavelet power error) between the measured hydrograph and both the simulated error and modelled hydrographs; and (b) the wavelet phase. Power difference and wavelet phase were used to develop two objective functions, RMSE(power) and RMS(phase), which were shown to distinguish between simulated errors and model predictions with similar values of the commonly adopted Nash‐Sutcliffe efficiency index. These objective functions suffer because they do not retain time, frequency or time‐frequency localization. Consideration of wavelet power spectra and time‐ and frequency‐integrated power spectra shows that the impacts of different types of simulated error can be seen through retention of some localization, especially in relation to when and the scale over which error was manifest. Theoretical objections to the use of wavelet analysis for this type of application are noted, especially in relation to the dependence of findings upon the wavelet chosen. However, it is argued that the benefits of localization and the qualitatively low sensitivity of wavelet power and phase to wavelet choice are sufficient to warrant further exploration of wavelet‐based approaches to rainfall‐runoff model evaluation. Copyright © 2006 John Wiley & Sons, Ltd.     

Application of a data‐based mechanistic modelling (DBM) approach for predicting runoff generation in semi‐arid regions

This paper addresses the application of a data‐based mechanistic (DBM) modelling approach using transfer function models (TFMs) with non‐linear rainfall filtering to predict runoff generation from a semi‐arid catchment (795 km²) in Tanzania. With DBM modelling, time series of rainfall and streamflow were allowed to suggest an appropriate model structure compatible with the data available. The model structures were evaluated by looking at how well the model fitted the data, and how well the parameters of the model were estimated. The results indicated that a parallel model structure is appropriate with a proportion of the runoff being routed through a fast flow pathway and the remainder through a slow flow pathway. Finally, the study employed a Generalized Likelihood Uncertainty Estimation (GLUE) methodology to evaluate the parameter sensitivity and predictive uncertainty based on the feasible parameter ranges chosen from the initial analysis of recession curves and calibration of the TFM. Results showed that parameters that control the slow flow pathway are relatively more sensitive than those that control the fast flow pathway of the hydrograph. Within the GLUE framework, it was found that multiple acceptable parameter sets give a range of predictions. This was found to be an advantage, since it allows the possibility of assessing the uncertainty in predictions as conditioned on the calibration data and then using that uncertainty as part of the decision‐making process arising from any rainfall‐runoff modelling project. Copyright © 2001 John Wiley & Sons, Ltd.     

Simulating the effects of spatial configurations of agricultural ditch drainage networks on surface runoff from agricultural catchments

The study of runoff is a crucial issue because it is closely related to flooding, water quality and erosion. In cultivated catchments, agricultural ditch drainage networks are known to influence runoff. As anthropogenic elements, agricultural ditch drainage networks can therefore be altered to better manage surface runoff in cultivated catchments. However, the relationship between the spatial configuration, i.e. the density and the topology, of agricultural ditch drainage networks and surface runoff in cultivated catchments is not understood. We studied this relationship by using a random network simulator that was coupled to a distributed hydrological model. The simulations explored a large variety of spatial configurations corresponding to a thousand stochastic agricultural ditch drainage networks on a 6.4 km² Mediterranean cultivated catchment. Next, several distributed hydrological functions were used to compute water flow paths and runoff for each simulation. The results showed that (i) denser networks increased the drained volume and the peak discharge and decreased hillslopes runoff, (ii) greater network density did not affect the surface runoff any further above a given network density, (iii) the correlation between network density and runoff was weaker for small subcatchments (< 2 km²) where the variability in the drained area that resulted from changes in agricultural ditch drainage networks increased the variability of runoff and (iv) the actual agricultural ditch drainage network appeared to be well optimized for managing runoff as compared with the simulated networks. Finally, our results highlighted the role of agricultural ditch drainage networks in intercepting and decreasing overland flow on hillslopes and increasing runoff in drainage networks. Copyright © 2011 John Wiley & Sons, Ltd.     

Effects of conditional parameterization on performance of rainfall‐runoff model regarding hydrologic non‐stationarity

Better parameterization of a hydrological model can lead to improved streamflow prediction. This is particularly important for seasonal streamflow forecasting with the use of hydrological modelling. Considering the possible effects of hydrologic non‐stationarity, this paper examined ten parameterization schemes at 12 catchments located in three different climatic zones in east Australia. These schemes are grouped into four categories according to the period when the data are used for model calibration, i.e. calibration using data: (1) from a fixed period in the historical records; (2) from different lengths of historical records prior to prediction year; (3) from different climatic analogue years in the past; and (4) data from the individual months. Parameterization schemes were evaluated according to model efficiency in both the calibration and verification period. The results show that the calibration skill changes with the different historic periods when data are used at all catchments. Comparison of model performance between the calibration schemes indicates that it is worth calibrating the model with the use of data from each individual month for the purpose of seasonal streamflow forecasting. For the catchments in the winter‐dominant rainfall region of south‐east Australia, a more significant shift in rainfall‐runoff relationships at different periods was found. For those catchments, model calibration with the use of 20 years of data prior to the prediction year leads to a more consistent performance. Copyright © 2011 John Wiley & Sons, Ltd.     

A semi‐distributed parallel‐type linear reservoir rainfall‐runoff model and its application in Taiwan

The main purpose of this paper is to introduce a semi‐distributed parallel surface rainfall‐runoff conceptual model. In this paper, a general solution of the instantaneous unit hydrograph (IUH) has been derived successfully for N linearly connected reservoirs, each having a different storage constant. The solution is a function of geomorphologic parameters, meteorologic factors and roughness coefficients. The model also takes into account the hydrologic response which is influenced by outflow downstream of a reservoir. For calibration, the shuffled complex evolution (SCE) algorithm is used to search for the global optimal parameters of the model. Because of the parallel structure, the mean roughness parameter of the channel becomes a “conceptual” parameter without a real physical meaning. To evaluate the adaptability of the model adopted, three watersheds around the city of Taipei in Taiwan were chosen to test the effectiveness of the model. The study provides an appropriate rainfall‐runoff model for planning flood mitigation in Taiwan. Copyright © 1999 John Wiley & Sons, Ltd.