Estimating the sensitivity of mean annual runoff to climate change using selected hydrological models

Hydrological model sensitivity to climate change can be defined as the response of a particular hydrological model to a known quantum of climate change. This paper estimates the hydrological sensitivity, measured as the percentage change in mean annual runoff, of two lumped parameter rainfall-runoff models, SIMHYD and AWBM and an empirical model, Zhang01, to changes in rainfall and potential evaporation. These changes are estimated for 22 Australian catchments covering a range of climates, from cool temperate to tropical and moist to arid. The results show that the models display different sensitivities to both rainfall and potential evaporation changes. The SIMHYD, AWBM and Zhang01 models show mean sensitivities of 2.4%, 2.5% and 2.1% change in mean annual flow for every 1% change in mean annual rainfall, respectively. All rainfall sensitivities have a lower limit of 1.8% and show upper limits of 4.1%, 3.4% and 2.5%, respectively. The results for potential evaporation change are −0.5%, −0.8% and −1.0% for every 1% increase in mean annual potential evaporation, respectively, with changes rainfall being approximately 3–5 times more sensitive than changes in potential evaporation for each 1% change in climate. Despite these differences, the results show similar correlations for several catchment characteristics. The most significant relationship is between percent change in annual rainfall and potential evaporation to the catchment runoff coefficient. The sensitivity of both A and B factors decreases with an increasing runoff coefficient, as does the uncertainty in this relationship. The results suggest that a first-order relationship can be used to give a rough estimate of changes in runoff using estimates of change in rainfall and potential evaporation representing small to modest changes in climate. Further work will develop these methods further, by investigating other regions and changes on the subannual scale.     

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

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

The relationship between runoff rate and lag time and the effects of surface treatments at the plot scale

Abstract

The manner in which both the seasonal and regional variations in storm duration, intensity and inter-storm period manifest in the runoff response of agricultural water supply catchments is investigated. High-resolution rainfall data were analysed for a network of 17 raingauges located across the semiarid (200–500 mm year^?1) agricultural districts of southwest Western Australia. Seasonal variations in mean storm duration, mean rainfall intensity and mean inter-storm period were modelled using simple periodic functions whose parameters were then also regressed with geographic and climatic indices to create spatial fields for each of these statistics. Based on these mean values, a continuous rainfall time series can be synthesized for any location within the region, with the rainfall depth within each storm being downscaled to 5-min time steps using a bounded random cascade model. Runoff from six different catchment surface treatments (“engineered” catchments) was simulated using a conceptual water-balance model, validated using rainfall—runoff data from an experimental field site. The expected yield of the various catchment types at any other location within the study region is then simulated using the above rainfall—runoff model and synthetic rainfall and potential evaporation time series under a range of climatic settings representative of regional climate variation. The resulting coupled model can be used to estimate the catchment area required to yield an acceptable volume of runoff for any location and dam capacity, at a specified reliability level, thus providing a tool for water resource managers to design engineered catchments for water supply. Although the model presented is specific for Western Australia's southwest region, the methodology itself is applicable to other locations.     

Controls on event runoff coefficients in the eastern Italian Alps

Analyses of event runoff coefficients provide essential insight on catchment response, particularly if a range of catchments and a range of events are compared by a single indicator. In this study we examine the effect of climate, geology, land use, flood types and initial soil moisture conditions on the distribution functions of the event runoff coefficients for a set of 14 mountainous catchments located in the eastern Italian Alps, ranging in size from 7.3 to 608.4 km². Runoff coefficients were computed from hourly precipitation, runoff data and estimates of snowmelt. A total of 535 events were analysed over the period 1989–2004. We classified each basin using a “permeability index” which was inferred from a geologic map and ranged from “low” to “high permeability”. A continuous soil moisture accounting model was applied to each catchment to classify ‘wet’ and ‘dry’ initial soil moisture conditions. The results indicate that the spatial distribution of runoff coefficients is highly correlated with mean annual precipitation, with the mean runoff coefficient increasing with mean annual precipitation. Geology, through the ‘permeability index’, is another important control on runoff coefficients for catchments with mean annual precipitation less than 1200 mm. Land use, as indexed by the SCS curve number, influences runoff coefficient distribution to a lesser degree. An analysis of the runoff coefficients by flood type indicates that runoff coefficients increase with event snowmelt. Results show that there exists an intermediate region of subsurface water storage capacity, as indexed by a flow–duration curve-based index, which maximises the impact of initial wetness conditions on the runoff coefficient. This means that the difference between runoff coefficients characterised by wet and dry initial conditions is negligible both for basins with very large storage capacity and for basins with small storage capacity. For basins with intermediate storage capacities, the impact of the initial wetness conditions may be relatively large.     

Estimation of actual evaporation using precipitation and soil moisture records in arid climates

There is little knowledge available about infiltration and evaporation processes in wadi channels in arid regions. This work was conducted to determine the actual evaporation from bare soils in wadi channels in the south-western region of the Kingdom of Saudi Arabia. The estimation of soil evaporation is highly dependent on the availability of moisture in the upper layers of alluvial wadis, in which the areal rainfall, flood hydrograph and soil properties play a significant part. The study was conducted by estimating the actual evaporation using soil moisture data, precipitation and runoff depths in a representative basin. The results are compared with potential rates. The actual rates were 1.5 mm/day immediately after a rainy day and then decreased to 0.42 mm/day. The minimum rate was about 0.1–0.2 mm/day during the dry season. The potential rates were about 9.5 mm/day in June and July, decreasing to 3.5 mm/day in December and January.     

Comparison of varied complexity models simulating recharge at the field scale

Numerical models are frequently used for the regional quantification of groundwater recharge. However there is a wide range of potential models available that represent the land surface with varying degrees of complexity, but which are rarely tested against observations at the field scale. We compared four models that simulate potential recharge at four intensively monitored sites with different vegetation and soil types in two adjacent catchments. These models were: Penman–Grindley, UN Food and Agricultural Organization, SPAtial Distributed Evaporation and Joint UK Land Environment Simulator. Standardized, unoptimized land surface datasets and pertinent literature were used for parameterization to reflect practice in regional water resource management and planning in the UK. The models were validated against soil moisture observations at all sites, as well as observed transpiration and interception and calculated total evaporation over a year at a woodland site. Soil moisture observations were generally reproduced well, but there were significant differences in how the models apportioned precipitation through the hydrological cycle. This demonstrates that soil moisture data alone are not a good diagnostic for groundwater recharge models. Significant differences in potential recharge were produced by models at both grassland sites, although simulated average annual potential recharge varied by only 15% at the grassland site on permeable soil. At the woodland sites, soil moisture contents were reproduced least accurately, and there were large differences in potential recharge at both woodland sites. This predominantly resulted from varied and inaccurate simulation of evaporation, particularly in the form of interception losses where this was explicitly represented in models. Differences in model structure, such as runoff representation, and parameter selection also influenced all results. Hydrological Processes © 2013 John Wiley & Sons, Ltd.     

Observation and estimation of evaporation from the ground surface of the cryosphere in eastern Asia

The characteristics of evaporation from the ground surface of Asian cryosphere sites are presented, as estimated by the lysimeter method, a profile method, and a heat budget method. The observation sites were located on the eastern Tibetan Plateau, in the Qilian and Tianshan Mountains of China, and in eastern Siberia. The lysimeter method has been demonstrated to be a reliable observation technique for estimating daily evaporation from the land surface, given suitable experiment design and operation. Daily mean evaporation varied within the range of 0·3 to 3·5 mm on the permafrost surface, and regional differences in evaporation were strongly related to surface soil moisture. Locally, topography, by way of its influence on surface soil moisture, was found to control evaporation systematically. Seasonality of ground evaporation in permafrost regions is dominated by thaw–freeze cycles at the surface; evaporation from the melting permafrost surface is up to four to seven times greater than that from frozen ground. In forested terrain, the interception of precipitation can reduce daily evaporation by 60 to 70%. Sublimation from the snow surface was observed at some sites in the range of 0·2 to 1·0 mm daily; atmospheric conditions, such as wind speed and saturation deficit, were dominant factors in determining snow sublimation. Copyright © 2003 John Wiley & Sons, Ltd.     

The influence of antecedent rainfall on Sahelian surface evaporation

This study quantifies the influence of rainfall on surface evaporation in the Sahel. A numerical model of the surface is used to extend the observations taken during the HAPEX–Sahel project, and is forced by 2 years of rain‐gauge data. The model is applied to the Southern Super Site (SSS), which covers an area of less than 100 km². The effects of rainfall variability (spatial and temporal) on soil moisture, vegetation growth and evaporation are explored. Contrasting rainfall conditions between the two years produce observed differences in the timing of the seasonal growth cycle. This correlates well with modelled root‐zone moisture deficits, and exerts a modest influence on transpiration rates. The evolution of surface evaporation is dominated, however, by the bare soil contribution in the day or two after a storm. This component also exerts a strong influence on the spatial variability of fluxes across the SSS, particularly when rain falls only in part of the area. In these cases, differences in evaporation between recently wetted and dry areas can reach 3\5 mm day⁻¹. Observations indicate that during a period of persistent rainfall gradients across the SSS, the lower atmosphere maintained a ‘memory’ of past rainfall patterns through humidity contrasts. These were the result of gradients in surface soil moisture, and therefore evaporation. The model results therefore support the possibility of a positive surface feedback mechanism affecting rainfall patterns in the region. Copyright © 2000 John Wiley & Sons, Ltd.     

Development of a one-parameter variable source area runoff model for ungauged basins

This research develops a one-parameter model of saturated source area dynamics and the spatial distribution of soil moisture. The single required parameter is the maximum soil moisture deficit within the catchment. The concept behind the development of the model comes from the fact that the complexity of topographically-driven runoff generation can be reduced through the use of geomorphological scaling relations. The scaling formulation allows the prediction of the dynamics of saturated source areas as a function of basin-wide soil moisture state. This model offers a number of potential advantages. Firstly, the model parameter is independent of topographic index distribution and its associated scale effects. Secondly, it may be possible to measure this single parameter using field measurements or perhaps remote sensing, which gives the model significant potential for application in ungauged basins. Finally, the fact that this parameter is a physical characteristic of the basin, estimation of this parameter avoids regionalization and parameter transferability problems. The model is tested using rainfall–runoff data from the 10.4 ha experimental catchment known as Tarrawara in Australia, the 37 km² Town Creek catchment in U.S.A., and the 620 km² Balaphi and the 850 km² Likhu sub-catchments of the Koshi river in Nepal. In sub-catchments of Koshi river, the simulation results compare favorably against the calibrated TOPMODEL both in terms of direct runoff and the spatial distribution of soil moisture state. In the Tarrawara and Town Brook catchments, simulation results compare favorably against observed storm runoff using all observed data, without calibration.     

Mapping of erosion rates in marly badlands based on a coupling of anatomical changes in exposed roots with slope maps derived from LiDAR data

Black marls form very extensive outcrops in the Alps and constitute some of the most eroded terrains, thus causing major problems of sedimentation in artificial storage systems (e.g. reservoirs) and river systems. In the experimental catchments near Draix (France), soil erosion rates have been measured in the past at the plot scale through a detailed monitoring of surface elevation changes and at the catchment scale through continuous monitoring of sediment yield in traps at basin outlets. More recently, erosion rates have been determined by means of dendrogeomorphic techniques in three monitored catchments of the Draix basin. A total of 48 exposed roots of Scots pine have been sampled and anatomical variations in annual growth rings resulting from denudation analysed. At the plot scale, average medium‐term soil erosion rates derived from exposed roots vary between 1·8 and 13·8 mm yr^?1 (average: 5·9 mm yr^?1) and values are significantly correlated with slope angle. The dendrogeomorphic record of point‐scale soil erosion rates matches very well with soil erosion rates measured in the Draix basins. Based on the point‐scale measurements and dendrogeomorphic results obtained at the point scale, a linear regression model involving slope angle was derived and coupled to high‐resolution slope maps obtained from a LiDAR‐generated digital elevation model so as to generate high‐resolution soil erosion maps. The resulting regression model is statistically significant and average soil erosion rates obtained from the areal erosion map (5·8, 5·2 and 6·2 mm yr^?1 for the Roubine, Moulin and Laval catchments, respectively) prove to be well in concert with average annual erosion rates measured in traps at the outlet of these catchments since 1985 (6·3, 4·1 and 6·4 mm yr^?1). This contribution demonstrates that dendrogeomorphic analyses of roots clearly have significant potential and that they are a powerful tool for the quantification and mapping of soil erosion rates in areas where measurements of past erosion is lacking. Copyright © 2011 John Wiley & Sons, Ltd.     

Hydroclimatology of the Volta River Basin in West Africa: Trends and variability from 1901 to 2002

Long-term historical records of rainfall (P), runoff (Q) and other climatic factors were used to investigate hydrological variability and trends in the Volta River Basin over the period 1901-2002. Potential (E_p) and actual evaporation (E), rainfall variability index (δ), Budyko’s aridity index (I_A), evaporation ratio (C_E) and runoff ratio (C_Q) were estimated from the available hydroclimatological records. Mann-Kendall trend analysis and non-parametric Sen’s slope estimates were performed on the respective time series variables to detect monotonic trend direction and magnitude of change over time.Rainfall variability index showed that 1968 was the wettest year (δ = +1.75) while 1983 was the driest (δ = −3.03), with the last three decades being drier than any other comparable period in the hydrological history of the Volta. An increase of 0.2 mm/yr² (P < 0.05) was observed in E_p for the 1901-1969 sub-series while an increased of 1.8 mm/yr² (P < 0.01) was recorded since 1970. Rainfall increased at the rate of 0.7 mm/yr² or 49 mm/yr between 1901 and 1969, whereas a decrease of 0.2 mm/yr² (6 mm/yr) was estimated for 1970-2002 sub-series. Runoff increased significantly at the rate of 0.8 mm/yr (23 mm/yr) since 1970. Runoff before dam construction was higher (87.5 mm/yr) and more varied (CV = 41.5%) than the post-dam period with value of 73.5 mm/yr (CV = 23.9%). A 10% relative decrease in P resulted in a 16% decrease in Q between 1936 and 1998. Since 1970, all the months showed increasing runoff trends with significant slopes (P < 0.05) in 9 out of the 12 months. Possible causes, such as climate change and land cover change, on the detected changes in hydroclimatology are briefly discussed.     

Estimates of spatial variation in evaporation using satellite‐derived surface temperature and a water balance model

Evaporation dominates the water balance in arid and semi‐arid areas. The estimation of evaporation by land‐cover type is important for proper management of scarce water resources. Here, we present a method to assess spatial and temporal patterns of actual evaporation by relating water balance evaporation estimates to satellite‐derived radiometric surface temperature. The method is applied to a heterogeneous landscape in the Krishna River basin in south India using 10‐day composites of NOAA advanced very high‐resolution radiometer satellite imagery. The surface temperature predicts the difference between reference evaporation and modelled actual evaporation well in the four catchments (r² = 0·85 to r² = 0·88). Spatial and temporal variations in evaporation are linked to vegetation type and irrigation. During the monsoon season (June–September), evaporation occurs quite uniformly over the case‐study area (1·7–2·1 mm day^?1), since precipitation is in excess of soil moisture holding capacity, but it is higher in irrigated areas (2·2–2·7 mm day^?1). In the post‐monsoon season (December–March) evaporation is highest in irrigated areas (2·4 mm day^?1). A seemingly reasonable estimate of temporal and spatial patterns of evaporation can be made without the use of more complex and data‐intensive methods; the method also constrains satellite estimates of evaporation by the annual water balance, thereby assuring accuracy at the seasonal and annual time‐scales. Copyright © 2007 John Wiley & Sons, Ltd.     

Factors controlling sediment yield in a major South American drainage basin: the Magdalena River,Colombia

The Magdalena River, a major fluvial system draining most of the Colombian Andes, has the highest sediment yield of any medium-sized or large river in South America. We examined sediment yield and its response to control variables in the Magdalena drainage basin based on a multi-year dataset of sediment loads from 32 tributary catchments. Various morphometric, hydrologic, and climatic variables were estimated in order to understand and predict the variation in sediment yield. Sediment yield varies from 128 to 2200 t km⁻² yr⁻¹ for catchments ranging from 320 to 59,600 km². The mean sediment yield for 32 sub-basins within the Magdalena basin is ∼690 t km⁻² yr⁻¹. Mean annual runoff is the dominant control and explains 51% of the observed variance in sediment yield. A multiple regression model, including two control variables, runoff and maximum water discharge, explains 58% of the variance. This model is efficient (ME=0.89) and is a valuable tool for predicting total sediment yield from tributary catchments in the Magdalena basin. Multiple correlations for those basins corresponding to the upper Magdalena, middle basin, Eastern Cordillera, and catchment areas greater than 2000 km², explain 75, 77, 89, and 78% of the variance in sediment yield, respectively. Although more variance is explained when dataset are grouped into categories, the models are less efficient (ME<0.72). Within the spatially distributed models, six catchment variables predict sediment yield, including runoff, precipitation, precipitation peakedness, mean elevation, mean water discharge, and relief. These estimators are related to the relative importance of climate and weathering, hillslope erosion, and fluvial transport processes. Time series analysis indicates that significant increases in sediment load have occurred over 68% of the catchment area, while 31% have experienced a decreasing trend in sediment load and thus yield. Land use analysis and increasing sediment load trends indicate that erosion within the catchment has increased over the last 10–20 years.     

Remotely sensed ET for streamflow modelling in catchments with contrasting flow characteristics: an attempt to improve efficiency

An efficient calibration with remotely sensed (RS) data is important for accurate predictions at ungauged catchments. This study investigates the advantages of streamflow-sensitive regionalization on calibration with RS evapotranspiration (ET). Regionalization experiments are performed at 28 catchments in Australia. The catchments are classified into three groups based on annual rainfall and runoff coefficients. Streamflow, RS ET, and a multi-objective RS ET-streamflow calibration are performed using the DiffeRential Evolution Adaptive Metropolis algorithm in each catchment. Simplified Australian Water Resource Assessment-Landscape model is calibrated for a selection of five parameters. Posterior probability distributions of parameters from three calibrations performed at donor catchments in each group are inspected to find the parameter for regionalization in the individual group. In group 1 of wetter catchments, regionalization of parameter FsoilEmax (soil evaporation scaling factor) helps to simplify the calibration without any deterioration in ET, soil moisture (SM) and streamflow predictions. Regionalization of parameter Beta (coefficient describing rate of hydraulic conductivity increase with water content) in group 2 assists to improve the streamflow predictions with no decrement in ET and SM predictions. However, regionalization is not able to provide satisfactory results in group 3. Group 3 includes low-yielding catchments, with average annual rainfall below 1000 mm/year and runoff coefficient less than 0.1, where traditional streamflow calibration also fails to produce accurate results. This study concludes that streamflow-sensitive regionalization is effective for improving the efficacy of RS ET calibration in wetter catchments.     

Effects of land‐cover changes on the hydrological response of interior Columbia River basin forested catchments

The topographically explicit distributed hydrology–soil–vegetation model (DHSVM) is used to simulate hydrological effects of changes in land cover for four catchments, ranging from 27 to 1033 km², within the Columbia River basin. Surface fluxes (stream flow and evapotranspiration) and state variables (soil moisture and snow water equivalent) corresponding to historical (1900) and current (1990) vegetation are compared. In addition a sensitivity analysis, where the catchments are covered entirely by conifers at different maturity stages, was conducted. In general, lower leaf‐area index (LAI) resulted in higher snow water equivalent, more stream flow and less evapotranspiration. Comparisons with the macroscale variable infiltration capacity (VIC) model, which parameterizes, rather than explicitly represents, topographic effects, show that runoff predicted by DHSVM is more sensitive to land‐cover changes than is runoff predicted by VIC. This is explained by model differences in soil parameters and evapotranspiration calculations, and by the more explicit representation of saturation excess in DHSVM and its higher sensitivity to LAI changes in the calculation of evapotranspiration. Copyright © 2002 John Wiley & Sons, Ltd.     

River water and nutrient discharges in the Northern Adriatic Sea: Current importance and long term changes

Runoff and nutrient transport by rivers were analysed in the Northern Adriatic continental shelf, in order to evaluate their interannual and multidecal variability, as well as their current contribution to determine freshwater and nutrient budgets in this marine region. During the years 2004-2007, the runoff in the basin (34.1-64.6 km³ yr⁻¹) was highly imbalanced, being 84% of freshwater discharged along the western coast, because of the contributions of Po, Adige and Brenta rivers. In the northern and eastern sections of the coast, freshwater discharge by rivers was less important (10 and 6%, respectively), but not negligible in determining the oceanographic properties at sub-regional scales. The oscillations of the transport of biogenic elements (124-262×10³ t N yr⁻¹ for TN, 72-136×10³ t N yr⁻¹ for DIN, 4.5-11.1×10³t P yr⁻¹ for TP, 2.2-3.5×10³ t P yr⁻¹ for PO₄ and 104-196×10³ t Si yr⁻¹ for SiO₂) were strictly dependant to the differences in the annual runoff. A strong excess of N load in comparison to P load characterised all rivers, both in inorganic nutrient (DIN/PO4=37-418) and total (TN/TP=48-208) pools, particularly in the northern and eastern areas of the basin.The annual runoff showed significant oscillations for Po on multidecadal time scale, whereas a general decrease (−33%) was observed for the other N Adriatic rivers as the recent discharges were compared to those before the 1980s. During the dry years 2005-2007, a strong reduction of river water flows and nutrient loads was experienced by the N Adriatic ecosystem with respect to years characterised by medium-high regimes. An increased frequency of similar drought periods, due to ongoing climate changes or to a larger human usage of continental waters, would be easily able to significantly change the biogeochemistry of this basin.     

A temperature‐precipitation‐based model of thirty‐year mean snowpack accumulation and melt in Oregon,USA

High‐resolution, spatially extensive climate grids can be useful in regional hydrologic applications. However, in regions where precipitation is dominated by snow, snowmelt models are often used to account for timing and magnitude of water delivery. We developed an empirical, nonlinear model to estimate 30‐year means of monthly snowpack and snowmelt throughout Oregon. Precipitation and temperature for the period 1971–2000, derived from 400‐m resolution PRISM data, and potential evapotranspiration (estimated from temperature and day length) drive the model. The model was calibrated using mean monthly data from 45 SNOTEL sites and accurately estimated snowpack at 25 validation sites: R² = 0·76, Nash‐Sutcliffe Efficiency (NSE) = 0·80. Calibrating it with data from all 70 SNOTEL sites gave somewhat better results (R² = 0·84, NSE = 0·85). We separately applied the model to SNOTEL stations located < 200 and ≥ 200 km from the Oregon coast, since they have different climatic conditions. The model performed equally well for both areas. We used the model to modify moisture surplus (precipitation minus potential evapotranspiration) to account for snowpack accumulation and snowmelt. The resulting values accurately reflect the shape and magnitude of runoff at a snow‐dominated basin, with low winter values and a June peak. Our findings suggest that the model is robust with respect to different climatic conditions, and that it can be used to estimate potential runoff in snow‐dominated basins. The model may allow high‐resolution, regional hydrologic comparisons to be made across basins that are differentially affected by snowpack, and may prove useful for investigating regional hydrologic response to climate change. Published in 2011 by John Wiley & Sons, Ltd.     

The role of overland flow and subsurface flow on the spatial distribution of soil moisture in the topsoil

Many investigations show relationships between topographical factors and the spatial distribution of soil moisture in catchments. However, few quantitative analyses have been carried out to elucidate the role of different hydrological processes in the spatial distribution of topsoil moisture in catchments. A spatially distributed rainfall—runoff model was used to investigate contributions of subsurface matric flow, macropore flow and surface runoff to the spatial distribution of soil moisture in a cultivated catchment. The model results show that lateral subsurface flow in the soil matrix or in macropores has a minor effect on the spatial distribution of soil moisture. Only when a perched groundwater table is maintained long enough, which is only possible if the subsurface is completely impermeable, may a spatial distribution in moisture content occur along the slope. Surface runoff, producing accumulations of soil moisture in flat flow paths of agricultural origin (field boundaries), was demonstrated to cause significant spatial variations in soil moisture within a short period after rainfall (<2 days). When significant amounts of surface runoff are produced, wetter moisture conditions will be generated at locations with larger upstream contributing areas. Copyright © 2001 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.