How well can the subsurface storage–discharge relation be interpreted and predicted using the geometric factors in headwater areas?

Headwater storage–discharge (S–Q) remains one of the least understood processes, and there is renewed interest in the S–Q relation. How well can the S–Q relation be interpreted mechanistically using geometric factors? In this paper, the hillslope storage Boussinesq and hillslope storage kinematic wave equation were adopted to guide the theoretical derivations. Analytical solutions were derived based on the hsKW equation for nine idealized hillslope aquifers, which were subdivided into two groups, i.e. hillslope aquifers with exponential hillslope width function (C1) and hillslope aquifers with Gaussian hillslope width function (C2). We found that analytical expressions of the S–Q relation can be derived for C1 hillslope aquifers. For more compound hillslope aquifers, i.e. C2, no explicit S–Q relation can be obtained. The whole subsurface recession after a rainstorm is simulated by applying the initial saturation condition. We found that the simulated S–Q processes can be characterized by a two‐phase recession, i.e. quick and slow recession. The time (t_b) at the dividing point of the quick and slow recessions depends on the geometric factors, such as the plan and profile curvature. In the quick recession for C1, many of the S–Q curves can be described as linear or quasi‐linear functions, which indicate that linear reservoir models can be applied approximately for recession simulations. However, during the slow recession phase of C1 and during the whole recession of C2, the S–Q relations are highly non‐linear. Finally, we compared the hillslope storage kinematic wave and hillslope storage Boussinesq models for simulating subsurface water recession after a rainstorm event in a real‐world headwater catchment (G5) in China. Through comparison of the recession slope curves, we found that the simulated results of the models employing the Gaussian hillslope width function match the observed hydrograph. The results indicated that appropriate organization of the hillslope geometric factors enhances our ability to make S–Q predictions. Copyright © 2016 John Wiley & Sons, Ltd.     

The introduction of runoff routing into large-scale hydrological models

Two alternative schemes are presented that are appropriate for the representation of runoff routing in large-scale grid-based hydrological models and atmospheric general circulation models (AGCMs). The first scheme characterizes routing processes as a single conceptual store. The second scheme, developed by Naden (1992), uses the normalized network width function to characterize the channel network form and a linear solution to the convective diffusion equation of one-dimensional flow to characterize the routing effect of a single channel. Both schemes are applied to the Severn catchment at the daily time-scale for the period 1981 to 1990 using a grid resolution of 40 km. Comparable results were obtained using both schemes (efficiencies were of the order of 80% in both cases). A combined model using a conceptual reservoir to represent hillslope routing and the network-based scheme to represent channel routing was developed to investigate the relative roles of hillslope and channel routing at the catchment scale. The application of this model demonstrated the important role of hillslope routing in reproducing the low frequency component of the catchment response. However, in terms of goodness-of-fit there was little to choose between the three schemes. Consequently, it is recommended that additional a priori knowledge of the routing processes should be used to condition the choice of model structure. © 1997 John Wiley & Sons, Ltd.     

Hydrological flow paths at the allt A' mharcaidh catchment: An analysis of plot and catchment scale observations

An integrated programme of hydrological monitoring at the 10 km² Allt a' Mharcaidh catchment in north-east Scotland has been based on observations at plot, hillslope and catchment scale. The resonse of the principal soil types has been characterized from a combination of throughflow and three-dimensional tensiometer data at plot scale, and plot sequences have been used to investigate hillslope scale effects. Seep emergence is associated with downslope drainage and local topographic convergence; in parallel preferential pathways generate a highly dynamic throughflow response. Catchment and subcatchment hydrographs mirror the twin dynamic observed at hillslope scale, and a unified hypothesis of response is presented which is consistent with all scales of observations.     

A watershed rainfall data recovery approach with application to distributed hydrological models

Watershed areal rainfall estimation, which is one of the most important and fundamental aspects in hydrological forecasting and various kinds of catchment‐scale hydrological models, is widely used in the analysis of hydrological regime change, and its precision has a direct influence on the accuracy of hydrological forecasting and hydrological simulation. In China, it is difficult to obtain the watershed areal rainfall estimate with reliable precision and avoid the phenomenon of ‘the same effect of different parameters’ because of the low density of the rain gauge network. Therefore, a watershed rainfall data recovery approach of improving the precision of watershed areal rainfall estimation is proposed here. This approach is to build new observatories, establish the time–space relations of rainfall between newly built observatories and previously built observatories in a relatively short interval and then recover the rainfall data of newly built observatories prior to their construction through simulating the relations over a longer time. As a result, watershed rainfall information could be elaborated to improve the precision of watershed areal rainfall estimate and avoid the phenomenon of ‘the same effect of different parameters’ to a certain degree in the process of hydrological simulation. The approach is used in the hydrological simulation of Hali River catchment. In combination with the Soil Water Assessment Tool model, a better result can be obtained in the hydrological simulation. Therefore, the approach can be used in other similar catchments. Copyright © 2011 John Wiley & Sons, Ltd.     

Modelling of sediment dynamics in a laboratory‐scale experimental catchment

The variability of hillslope form and function is examined experimentally using a simple model catchment in which most landscape development parameters are either known or controlled. It is demonstrated that there is considerable variability in sediment output from similar catchments, subjected to the same hydrological processes, and for which the initial hillslope profiles are the same. The results demonstrate that, in the case of catchments with a linear initial hillslope profile, the sediment output is initially high but reduces through time, whereas for a concave initial profile the sediment output was smaller and relatively constant. Concave hillslope profiles also displayed reduced sediment output when compared with linear slopes with the same overall slope. Using this experimental model catchment data, the SIBERIA landscape evolution model was tested for its ability to predict temporal sediment transport. When calibrated for the rainfall and erodible material, SIBERIA is able to simulate mean temporal sediment output for the experimental catchment over a range of hillslope profiles and rainfall intensities. SIBERIA is also able to match the hillslope profile of the experimental catchments. The results of the study provide confidence in the ability of SIBERIA to predict temporal sediment output. The experimental and modelling data also demonstrate that, even with all geomorphic and hydrological variables being known and/or controlled, there is still a need for long‐term stream gauging to obtain reliable assessments of field catchment hydrology and sediment transport. Copyright © 2005 John Wiley & Sons, Ltd.     

A hillslope-based hydrological model using catchment area and width functions

Abstract

A flow-interval hillslope discretization scheme is proposed for catchment hydrological modelling. By this scheme, a two-dimensional catchment is simplified into a one-dimensional cascade of flow intervals linked by the main stream. Each flow interval comprises a set of parallel hillslopes. The hillslope is the fundamental computational unit in the hydrological model providing lateral inflow to the main stream. The size of hillslope is determined by the catchment area and width functions. Catchment runoff is the total of hillslope responses through the river routing. Tests in four Japanese catchments showed that the model performed well on simulating the overall water balance, general flow pattern, and daily and hourly hydrographs of a whole catchment, as well as simultaneous simulation in different subcatchments. Characteristics of catchment hydrological responses and model applicability are discussed.     

Interactions and connectivity between runoff generation processes of different spatial scales

Monitoring runoff generation processes in the field is a prerequisite for developing conceptual hydrological models and theories. At the same time, our perception of hydrological processes strongly depends on the spatial and temporal scale of observation. Therefore, the aim of this study is to investigate interactions between runoff generation processes of different spatial scales (plot scale, hillslope scale, and headwater scale). Different runoff generation processes of three hillslopes with similar topography, geology and soil properties, but differences in vegetation cover (grassland, coniferous forest, and mixed forest) within a small v‐shaped headwater were measured: water table dynamics in wells with high spatial and temporal resolution, subsurface flow (SSF) of three 10 m wide trenches at the bottom of the hillslopes subdivided into two trench sections each, overland flow at the plot scale, and catchment runoff. Bachmair et al. ( 2012 ) found a high spatial variability of water table dynamics at the plot scale. In this study, we investigate the representativity of SSF observations at the plot scale versus the hillslope scale and vice versa, and the linkage between hillslope dynamics (SSF and overland flow) and streamflow. Distinct differences in total SSF within each 10 m wide trench confirm the high spatial variability of the water table dynamics. The representativity of plot scale observations for hillslope scale SSF strongly depends on whether or not wells capture spatially variable flowpaths. At the grassland hillslope, subsurface flowpaths are not captured by our relatively densely spaced wells (3 m), despite a similar trench flow response to the coniferous forest hillslope. Regarding the linkage between hillslope dynamics and catchment runoff, we found an intermediate to high correlation between streamflow and hillslope hydrological dynamics (trench flow and overland flow), which highlights the importance of hillslope processes in this small watershed. Although the total contribution of SSF to total event catchment runoff is rather small, the contribution during peak flow is moderate to substantial. Additionally, there is process synchronicity between spatially discontiguous measurement points across scales, potentially indicating subsurface flowpath connectivity. Our findings stress the need for (i) a combination of observations at different spatial scales, and (ii) a consideration of the high spatial variability of SSF at the plot and hillslope scale when designing monitoring networks and assessing hydrological connectivity. Copyright © 2013 John Wiley & Sons, Ltd.     

Sensitivity of hydrological models to uncertainty in rainfall input

Abstract

Different approaches used in hydrological modelling are compared in terms of the way each one takes the rainfall data into account. We examine the errors associated with accounting for rainfall variability, whether in hydrological modelling (distributed vs lumped models) or in computing catchment rainfall, as well as the impact of each approach on the representativeness of the parameters it uses. The database consists of 1859 rainfall events, distributed on 500 basins, located in the southeast of France with areas ranging from 6.2 to 2851 km². The study uses as reference the hydrographs computed by a distributed hydrological model from radar rainfall. This allows us to compare and to test the effects of various simplifications to the process when taking rainfall information (complete rain field vs sampled rainfall) and rainfall–runoff modelling (lumped vs distributed) into account. The results appear to show that, in general, the sampling effect can lead to errors in discharge at the outlet that are as great as, or even greater than, those one would get with a fully lumped approach. We found that small catchments are more sensitive to the uncertainties in catchment rainfall input generated by sampling rainfall data as seen through a raingauge network. Conversely, the larger catchments are more sensitive to uncertainties generated when the spatial variability of rainfall events is not taken into account. These uncertainties can be compensated for relatively easily by recalibrating the parameters of the hydrological model, although such recalibrations cause the parameter in question to completely lose physical meaning.

Citation Arnaud, P., Lavabre, J., Fouchier, C., Diss, S. & Javelle, P. (2011) Sensitivity of hydrological models to uncertainty of rainfall input. Hydrol. Sci. J. 56(3), 397–410.     

Capabilities and limitations of detailed hillslope hydrological modelling

Hillslope hydrological modelling is considered to be of great importance for the understanding and quantification of hydrological processes in hilly or mountainous landscapes. In recent years a few comprehensive hydrological models have been developed at the hillslope scale which have resulted in an advanced representation of hillslope hydrological processes (including their interactions), and in some operational applications, such as in runoff and erosion studies at the field scale or lateral flow simulation in environmental and geotechnical engineering. An overview of the objectives of hillslope hydrological modelling is given, followed by a brief introduction of an exemplary comprehensive hillslope model, which stimulates a series of hydrological processes such as interception, evapotranspiration, infiltration into the soil matrix and into macropores, lateral and vertical subsurface soil water flow both in the matrix and preferential flow paths, surface runoff and channel discharge. Several examples of this model are presented and discussed in order to determine the model's capabilities and limitations. Finally, conclusions about the limitations of detailed hillslope modelling are drawn and an outlook on the future prospects of hydrological models on the hillslope scale is given.The model presented performed reasonable calculations of Hortonian surface runoff and subsequent erosion processes, given detailed information of initial soil water content and soil hydraulic conditions. The vertical and lateral soil moisture dynamics were also represented quite well. However, the given examples of model applications show that quite detailed climatic and soil data are required to obtain satisfactory results. The limitations of detailed hillslope hydrological modelling arise from different points: difficulties in the representations of certain processes (e.g. surface crusting, unsaturated–saturated soil moisture flow, macropore flow), problems of small‐scale variability, a general scarcity of detailed soil data, incomplete process parametrization and problems with the interdependent linkage of several hillslopes and channel–hillslope interactions. Copyright © 1999 John Wiley & Sons, Ltd.     

A study on the impact of climate change on streamflow at the watershed scale in the humid tropics

Abstract

Climate change is recognized to be one of the most serious challenges facing mankind today. Driven by anthropogenic activities, it is known to be a direct threat to our food and water supplies and an indirect threat to world security. Increase in the concentration of carbon dioxide and other greenhouse gases in the atmosphere will certainly affect hydrological regimes. The consequent global warming is expected to have major implications on water resources management. The objective of this research is to present a general approach for evaluating the impacts of potential climate change on streamflow in a river basin in the humid tropical zone of India. Large-scale global climate models (GCMs) are the best available tools to provide estimates of the effect of rising greenhouse gases on rainfall and temperature. However the spatial resolution of these models (250 km?×?250 km) is not compatible with that of watershed hydrological models. Hence the outputs from GCMs have to be downscaled using regional climate models (RCMs), so as to project the output of a GCM to a finer resolution (50 km?×?50 km). In the present work, the projections of a GCM for two scenarios, A2 and B2 are downscaled by a RCM to project future climate in a watershed. Projections for two important climate variables, viz. rainfall and temperature are made. These are then used as inputs for a physically-based hydrological model, SWAT, in order to evaluate the effect of climate change on streamflow and vegetative growth in a humid tropical watershed.

Citation Raneesh, K. Y. & Santosh, G. T. (2011) A study on the impact of climate change on streamflow at the watershed scale in the humid tropics. Hydrol. Sci. J. 56(6), 946–965.     

Assessment on the function of reservoirs for flood control during typhoon seasons based on a distributed hydrological model

A grid‐based distributed hydrological model, PDTank model, is used to simulate hydrological processes in the upper Tone River catchment. The Tone River catchment often suffers from heavy rainfall events during the typhoon seasons. The reservoirs located in the catchment play an important role in flood regulation. Through the coupling of the PDTank model and a reservoir module that combines the storage function and operation function, the PDTank model is used for flood forecasting in this study. By comparing the hydrographs simulated using gauging and radar rainfall data, it is found that the spatial variability of rainfall is an important factor for flood simulation and the accuracy of the hydrographs simulated using radar rainfall data is slightly improved. The simulation of the typhoon flood event numbered No. 9 shows that the reservoirs in the catchment attenuate the peak flood discharge by 423·3 m³/s and validates the potential applicability of the distributed hydrological model on the assessment of function of reservoirs for flood control during typhoon seasons. Copyright © 2011 John Wiley & Sons, Ltd.     

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.     

Impact of climate change on hydrological processes over a basin scale in northern Taiwan

This study investigates the impact of climate change on rainfall, evapotranspiration, and discharge in northern Taiwan. The upstream catchment of the Shihmen reservoir in northern Taiwan was chosen as the study area. Both observed discharge and soil moisture were simultaneously adopted to optimize the HBV‐based hydrological model, clearly improving the simulation of the soil moisture. The delta change of monthly temperature and precipitation from the grid cell of GCMs (General Circulation Models) that is closest to the study area were utilized to generate the daily rainfall and temperature series based on a weather generating model. The daily rainfall and temperature series were further inputted into the calibrated hydrological model to project the hydrological variables. The studies show that rainfall and discharge will be increased during the wet season (May to October) and decreased during the dry season (November to April of the following year). Evapotranspiration will be increased in the whole year except in November and December. Copyright © 2009 John Wiley & Sons, Ltd.     

An analysis of the impact of spatial variability in rainfall on runoff and sediment predictions from a distributed model

This paper investigates the effect of introducing spatially varying rainfall fields to a hydrological model simulating runoff and erosion. Pairs of model simulations were run using either spatially uniform (i.e. spatially averaged) or spatially varying rainfall fields on a 500‐m grid. The hydrological model used was a simplified version of Thales which enabled runoff generation processes to be isolated from hillslope averaging processes. Both saturation excess and infiltration excess generation mechanisms were considered, as simplifications of actual hillslope processes. A 5‐year average recurrence interval synthetic rainfall event typical of temperate climates (Melbourne, Australia) was used. The erosion model was based on the WEPP interrill equation, modified to allow nonlinear terms relating the erosion rate to rainfall or runoff‐squared. The model results were extracted at different scales to investigate whether the effects of spatially varying rainfall were scale dependent. A series of statistical metrics were developed to assess the variability due to introducing the spatially varying rainfall field. At the catchment (approximately 150 km²) scale, it was found that particularly for saturation excess runoff, model predictions of runoff were insensitive to the spatial resolution of the rainfall data. Generally, erosion processes at smaller sub‐catchment scales, particularly when the sediment generation equation had non linearity, were more sensitive to spatial rainfall variability. Introducing runon infiltration reduced the total runoff and sediment yield at all scales, and this process was also most sensitive to the rainfall resolution. Copyright © 2011 John Wiley & Sons, Ltd.     

Effects of differential hillslope‐scale water retention characteristics on rainfall–runoff response at the Landscape Evolution Observatory

Hillslopes turn precipitation into runoff and thus exert important controls on various Earth system processes. It remains difficult to collect reliable data necessary for understanding and modeling these Earth system processes in real catchments. To overcome this problem, controlled experiments are being conducted at the Landscape Evolution Observatory at Biosphere 2, The University of Arizona. Previous experiments have revealed differences in hydrological response between 2 landscapes within Landscape Evolution Observatory, even though both landscapes were designed to be identical. In an attempt to discover where the observed differences stem from, we use a fully 3‐dimensional hydrological model (CATchment HYdrology) to show the effect of soil water retention characteristics and saturated hydraulic conductivity on the hydrological response of these 2 hillslopes. We also show that soil water retention characteristics can be derived at hillslope scale from experimental observations of soil moisture and matric potential. It is found that differences in soil packing between the 2 landscapes may be responsible for the observed differences in hydrological response. This modeling study also suggests that soil water retention characteristics and saturated hydraulic conductivity have a profound effect on rainfall–runoff processes at hillslope scale and that parametrization of a single hillslope may be a promising step in modeling rainfall–runoff response in real catchments.     

Evaluation of remote‐sensing‐based rainfall products through predictive capability in hydrological runoff modelling

The emergence of regional and global satellite‐based rainfall products with high spatial and temporal resolution has opened up new large‐scale hydrological applications in data‐sparse or ungauged catchments. Particularly, distributed hydrological models can benefit from the good spatial coverage and distributed nature of satellite‐based rainfall estimates (SRFE). In this study, five SRFEs with temporal resolution of 24 h and spatial resolution between 8 and 27 km have been evaluated through their predictive capability in a distributed hydrological model of the Senegal River basin in West Africa. The main advantage of this evaluation methodology is the integration of the rainfall model input in time and space when evaluated at the sub‐catchment scale. An initial data analysis revealed significant biases in the SRFE products and large variations in rainfall amounts between SRFEs, although the spatial patterns were similar. The results showed that the Climate Prediction Center/Famine Early Warning System (CPC‐FEWS) and cold cloud duration (CCD) products, which are partly based on rain gauge data and produced specifically for the African continent, performed better in the modelling context than the global SRFEs, Climate Prediction Center MORPHing technique (CMORPH), Tropical Rainfall Measuring Mission (TRMM) and Precipitation Estimation from Remotely Sensed Information using Artificial Neural Networks (PERSIANN). The best performing SRFE, CPC‐FEWS, produced good results with values of R²_NS between 0·84 and 0·87 after bias correction and model recalibration. This was comparable to model simulations based on traditional rain gauge data. The study highlights the need for input specific calibration of hydrological models, since major differences were observed in model performances even when all SRFEs were scaled to the same mean rainfall amounts. This is mainly attributed to differences in temporal dynamics between products. Copyright © 2009 John Wiley & Sons, Ltd.     

Spatial and temporal variability of the hydrological response in a small Mediterranean research catchment (Vallcebre,Eastern Pyrenees)

This paper analyses the spatial and temporal variability of the hydrological response in a small Mediterranean catchment (Cal Rodó). The first part of the analysis focuses on the rainfall–runoff relationship at seasonal and monthly scale, using an 8‐year data set. Then, using storm‐flow volume and coefficient, the temporal variability of the rainfall–runoff relationship and its relationship with several hydrological variables are analysed at the event scale from hydrographs observed over a 3‐year period. Finally, the spatial non‐linearity of the hydrological response is examined by comparing the Cal Rodó hydrological response with the Can Vila sub‐catchment response at the event scale. Results show that, on a seasonal and monthly scale, there is no simple relationship between rainfall and runoff depths, and that evapotranspiration is a factor that introduced some non‐linearity in the rainfall–runoff relationship. The analysis of monthly values also reveals the existence of a threshold in the relationship between rainfall and runoff depths, denoting a more contrasted hydrological response than the one usually observed in humid catchments. At the event scale, the storm‐flow coefficient has a clear seasonal pattern with an alternance between a wet period, when the catchment is hydrologically responsive, and a dry summer period, when the catchment is much less reactive to any rainfall. The relationship between the storm‐flow coefficient and rainfall depth, rainfall maximum intensity and base‐flow shows that observed correlations are the same as those observed for humid conditions, even if correlation coefficients are notably lower. Comparison with the Can Vila sub‐catchment highlights the spatial heterogeneity of the rainfall‐runoff relationship at the small catchment scale. Although interpretation in terms of runoff processes remains delicate, heterogeneities between the two catchments seem to be related to changes in the ratio between infiltration excess and saturation processes in runoff formation. Copyright © 2007 John Wiley & Sons, Ltd.     

Spatial scale effect on seasonal streamflows in permafrost catchments on the Qinghai–Tibet Plateau

The scale issue is of central concern in hydrological processes to understand the potential upscaling or downscaling methodologies, and to develop models for scaling the dominant processes at different scales and in different environments. In this study, a typical permafrost watershed in the Qinghai‐Tibet Plateau was selected. Its hydrological processes were monitored for 4 years from 2004 to 2008, measuring the effects of freezing and thawing depth of active soil layers on runoff processes. To identify the nature and cause of variation in the runoff response in different size catchments, catchments ranging from 1·07 to 112 km² were identified in the watershed. The results indicated that the variation of runoff coefficients showed a ‘V’ shape with increasing catchment size during the spring and autumn seasons, when the active soil was subjected to thawing or freezing processes. A two‐stage method was proposed to create runoff scaling models to indicate the effects of scale on runoff processes. In summer, the scaling transition model followed an exponential function for mean daily discharge, whereas the scaling model for flood flow exhibited a linear function. In autumn, the runoff process transition across multiple scales followed an exponential function with air temperature as the driving factor. These scaling models demonstrate relatively high simulation efficiency and precision, and provide a practical way for upscaling or downscaling runoff processes in a medium‐size permafrost watershed. For permafrost catchments of this scale, the results show that the synergistic effect of scale and vegetation cover is an important driving factor in the runoff response. Copyright © 2011 John Wiley & Sons, Ltd.     

Study of dynamic behaviour of recession curves

Since Brutsaert and Neiber (1977), recession curves are widely used to analyse subsurface systems of river basins by expressing ? dQ/dt as a function of Q, which typically take a power law form: ? dQ/dt = kQ^α, where Q is the discharge at a basin outlet at time t. Traditionally recession flows are modelled by single reservoir models that assume a unique relationship between ? dQ/dt and Q for a basin. However, recent observations indicate that ? dQ/dt ‐ Q relationship of a basin varies greatly across recession events, indicating the limitation of such models. In this study, the dynamic relationship between ? dQ/dt and Q of a basin is investigated through the geomorphological recession flow model which models recession flows by considering the temporal evolution of its active drainage network (the part of the stream network of the basin draining water at time t). Two primary factors responsible for the dynamic relationship are identified: (i) degree of aquifer recharge (ii) spatial variation of rainfall. Degree of aquifer recharge, which is likely to be controlled by (effective) rainfall patterns, influences the power law coefficient, k. It is found that k has correlation with past average streamflow, which confirms the notion that dynamic ? dQ/dt ‐ Q relationship is caused by the degree of aquifer recharge. Spatial variation of rainfall is found to have control on both the exponent, α, and the power law coefficient, k. It is noticed that that even with same α and k, recession curves can be different, possibly due to their different (recession) peak values. This may also happen due to spatial variation of rainfall. Copyright © 2012 John Wiley & Sons, Ltd.