Contributions to storm quickflow in a small headwater catchment—the role of natural pipes and soil macropores

Analysis of hydrographs from a 4·3 hectare stream head catchment indicates that storm runoff is generated from dynamic source areas. The volume and timing of contributions from different parts of the catchment show, when compared with the extent of surface saturation, that pipeflow generated from areas not saturated at the soil surface is a significant component of the quickflow hydrograph. A simple model of pipeflow generation and contribution is discussed in the light of field results.     

Testing the influence of topography and material properties on catchment‐scale soil moisture patterns using remotely sensed vegetation patterns in a humid temperate catchment,northern Britain

In order to evaluate the relationship between the apparent complexity of hillslope soil moisture and the emergent patterns of catchment hydrological behaviour and water quality, we need fine‐resolution catchment‐wide data on soil moisture characteristics. This study proposes a methodology whereby vegetation patterns obtained from high‐resolution orthorectified aerial photographs are used as an indicator of soil moisture characteristics. This enables us to examine a set of hypotheses regarding what drives the spatial patterns of soil moisture at the catchment scale (material properties or topography). We find that the pattern of Juncus effusus vegetation is controlled largely by topography and mediated by the catchment's material properties. Characterizing topography using the topographic index adds value to the soil moisture predictions relative to slope or upslope contributing area (UCA). However, these predictions depart from the observed soil moisture patterns at very steep slopes or low UCAs. Copyright © 2012 John Wiley & Sons, Ltd.     

THE INFLUENCE OF TOPOGRAPHY ON TIME AND SPACE DISTRIBUTION OF SOIL SURFACE WATER CONTENT

The distribution of water content in time and space at the soil surface has been investigated on a small farmland catchment (1.3 km² ) from four field surveys corresponding to different moisture statuses. For each survey, about 400 samples were collected at the soil surface at a depth of 5 cm along ten axes parallel to the greatest slope. The relationship between the measurements and the topography has been analysed. The structure of the data is well explained by a topographic index referring to the downslope conditions and defined as the elevation difference between the sample point and the stream point corresponding to the outlet of the water pathway derived from the digital elevation model (DEM). This index can be considered as an hydraulic head, at least for saturated conditions. A threshold for this index allows two domains within the catchment to be distinguished; an upper domain where the water content is nearly constant and varies slowly, and a lower domain where moisture status increases and is highly variable. The spatial distribution of these two domains is well correlated to the spatial distribution of the soils. Thus, both topography and the spatial distribution of soil appear to control the spatial distribution of surface water content at the 1-km² scale. © 1997 by John Wiley & Sons Ltd.     

Spatial variation of effective porosity and its implications for discharge in an upland headwater catchment in Scotland

Spatial and temporal measurements of shallow sub-surface soil physical properties were made within a 1 km² upland catchment. The surface soil layer of the catchment was organic rich (>70% organic matter) with a corresponding total porosity of 81%. Monthly point observations of volumetric water content (θ) were combined with point estimates of total porosity () and the porosity <50 μm (_residual), to define the ratio of water filled pore volume:pore volume in pores <50 μm (=θ/_residual). Values of θ/_residual were compared with discharge to test whether mass flow occurred when θ/_residual>1. A correlation between water content and discharge was found, with discharge increasing rapidly when θ/_residual approached unity. Similar relationships between water content and catchment discharge were identified for soil units adjacent to the stream when θ/_residual approached unity. These data suggest that soil pores >50 μm are of crucial importance in determining catchment discharge. Spatial and temporal variations in soil properties related to moisture content of the soil were also observed. Under dry conditions, a clear division based on aspect was noted, the west-facing side of the catchment being wettest. In wetter months, total porosity and soil water content were significantly affected by soil type and the spatial pattern of soil water content was more variable than in the dryer months. The physical quantification of soil properties in the shallow sub-surface layer proved important in explaining different initial changes in discharge from the catchment in response to a rainfall event.     

A hydrogeochemical model for stream chemistry,Cascade range,Washington, U.S.A.

A simple mixing model demonstrates that chemical variations in Cascade surface waters reflect flow from three general zones: alpine areas, forested colluvial slopes, and seasonally saturated areas. The chemistry of weathering solutions in alpine portions of the Williamson Creek catchment (North Cascade Range) results from alteration of plagioclase, hornblende, and biotite to kaolinitic material and vermiculite. Surface and shallow groundwater in forested portions of the catchment reflect these reactions, dissolution of small quantities of carbonate, and biologic activity. Both at-a-point and downstream chemical variations are explained quantitatively by the volume of water that originates in each of the hydrogeochemical source areas. Water from the forested colluvial slopes is most significant on an annual basis. However, summer low-flow is a mixture of colluvial waters and dilute solutions from the alpine zone, whereas 10 to 30 per cent of peak flow in snowmelt and rainstorms is produced from seasonally saturated areas. Poor concentration/discharge (C/Q) correlations, typical of Cascade rivers, result from mixing of significant C/Q relations for water leaving each source area. Model predictions could be substantially improved by better data for the effects of temperature, water-contact time, and biologic cycling on the chemistry of soil water from forested zones.     

Poro-mechanical coupling influences on potential for rainfall-induced shallow landslides in unsaturated soils

Rainfall-induced landslides are a common occurrence in terrain with steep topography and soils that have degradable strength. Rainfall infiltration into a partially saturated slope of infinite extent can lead to either a decrease or complete elimination of soil suction, compromising the slopes' stability. In this research the rainfall infiltration coupled with deformation of a partially saturated soil slope during rainfall infiltration is analyzed. The limit equilibrium conditions and the shear strength relationship of a partially saturated soil are employed to develop an analytical solution for calculating the stability of an infinite partially saturated slope due to rainfall infiltration. The analytical solutions are able to consider the influence of the coupled effects on the stability of the slope. The factors that affect the safety of a partially saturated slope of infinite extent are discussed. The results indicate that the poro-mechanical coupling of water infiltration and deformation has an important effect on the stability of the infinite unsaturated slope.     

Which catchment characteristics control the temporal dependence structure of daily river flows?

Hydrological classification systems seek to provide information about the dominant processes in the catchment to enable information to be transferred between catchments. Currently, there is no widely agreed‐upon system for classifying river catchments. This paper develops a novel approach to classifying catchments based on the temporal dependence structure of daily mean river flow time series, applied to 116 near‐natural ‘benchmark’ catchments in the UK. The classification system is validated using 49 independent catchments. Temporal dependence in river flow data is driven by the flow pathways, connectivity and storage within the catchment and can thus be used to assess the influence catchment characteristics have on moderating the precipitation‐to‐flow relationship. Semi‐variograms were computed for the 116 benchmark catchments to provide a robust and efficient way of characterising temporal dependence. Cluster analysis was performed on the semi‐variograms, resulting in four distinct clusters. The influence of a wide range of catchment characteristics on the semi‐variogram shape was investigated, including: elevation, land cover, physiographic characteristics, soil type and geology. Geology, depth to gleyed layer in soils, slope of the catchment and the percentage of arable land were significantly different between the clusters. These characteristics drive the temporal dependence structure by influencing the rate at which water moves through the catchment and/or the storage in the catchment. Quadratic discriminant analysis was used to show that a model with five catchment characteristics is able to predict the temporal dependence structure for un‐gauged catchments. This method could form the basis for future regionalisation strategies, as a way of transferring information on the precipitation‐to‐flow relationship between gauged and un‐gauged catchments. © 2014 The Authors. Hydrological Processes by published by John Wiley & Sons, Ltd.     

Dynamics of suspended sediment transport and yield in a large agricultural catchment,southwest France

The dynamics of suspended sediment transport were monitored continuously in a large agricultural catchment in southwest France from January 2007 to March 2009. The objective of this paper is to analyse the temporal variability in suspended sediment transport and yield in that catchment. Analyses were also undertaken to assess the relationships between precipitation, discharge and suspended sediment transport, and to interpret sediment delivery processes using suspended sediment‐discharge hysteresis patterns. During the study period, we analysed 17 flood events, with high resolution suspended sediment data derived from continuous turbidity and automatic sampling. The results revealed strong seasonal, annual and inter‐annual variability in suspended sediment transport. Sediment was strongly transported during spring, when frequent flood events of high magnitude and intensity occurred. Annual sediment transport in 2007 yielded 16 614 tonnes, representing 15 t km^?2 (85% of annual load transport during floods for 16% of annual duration), while the 2008 sediment yield was 77 960 tonnes, representing 70 t km^?2 (95% of annual load transport during floods for 20% of annual duration). Analysis of the relationships between precipitation, discharge and suspended sediment transport showed that there were significant correlations between total precipitation, peak discharge, total water yield, flood intensity and sediment variables during the flood events, but no relationship with antecedent conditions. Flood events were classified in relation to suspended sediment concentration (SSC)–discharge hysteretic loops, complemented with temporal dynamics of SSC–discharge ranges during rising and falling flow. The hysteretic shapes obtained for all flood events reflected the distribution of probable sediment sources throughout the catchment. Regarding the sediment transport during all flood events, clockwise hysteretic loops represented 68% from river deposited sediments and nearby source areas, anticlockwise 29% from distant source areas, and simultaneity of SSC and discharge 3%. Copyright © 2010 John Wiley & Sons, Ltd.     

A comparison of wetness indices for the prediction of observed connected saturated areas under contrasting conditions

For lack of other widely available spatial information, topography is often used to predict water fluxes and water quality in mesoscale watersheds. Such data have however proven to be misleading in many environments where large and flat valley bottoms and/or highly conducive soil covers determine water storage and water transport mechanisms. Also, the focus is generally on the prediction of saturation areas regardless of whether they are connected to the catchment hydrographic network or rather present in isolated topographic depressions. Here soil information was coupled with terrain data towards the targeted prediction of connected saturated areas. The focus was on the 30 km² Girnock catchment (Cairngorm Mountains, northeast Scotland) and its 3 km² sub‐catchment, Bruntland Burn in which seven field surveys were done to capture actual maps of connected saturated areas in both dry and humid conditions. The 1 km² resolution UK Hydrology of Soil Types (HOST) classification was used to extract relevant, spatially variable, soil parameters. Results show that connected saturated areas were fairly well predicted by wetness indices but only in wet conditions when they covered more than 30% of the whole catchment area. Geomorphic indices including information on terrain shape, steepness, aspect, soil texture and soil depth showed potential but generally performed poorly. Indices based on soil and topographic data did not have more predictive power than those based on topographic information only: this was attributed to the coarse resolution of the HOST classification. Nevertheless, analyses provided interesting insights into the scale‐dependent water storage and transport mechanisms in both study catchments. Copyright © 2013 John Wiley & Sons, Ltd.     

Saturated area dynamics and streamflow generation from coupled surface–subsurface simulations and field observations

A distributed physically-based model describing coupled surface–subsurface flows is applied to an instrumented catchment to investigate the links between runoff generation processes and the dynamics of saturated areas. The spatial characterization of the system is obtained through geophysical measurements and in situ observations. The model is able to reproduce the dynamics of the system through the calibration of only few parameters with a clear physical interpretation, providing a solid basis for our numerical investigations. Such investigations demonstrate the important control exerted by surface topography on the time evolution of saturated area patterns, mainly mediated by topographic curvature, that dictates both the dominant streamflow generation process at the local scale and the connection-disconnection dynamics of saturated areas. The relation between hillslope water storage and streamflow, Q = f(V), is shown to be highly hysteretical and dependent on the mean saturation of the catchment: higher degrees of saturation tend to yield one-to-one relationships between streamflow and water storage. On the contrary, streamflow-water storage relations are importantly affected by the specific configuration of saturated areas connected to the outlet when the system is far from complete saturation. This observation contradicts common assumptions of a one-to-one relationship Q = f(V) often used to justify widely observed power-law Q vs. dQ/dt recession curves. Furthermore, even when Q = f(V) becomes unique at high degrees of saturation, no power-law form emerged in our runs, speculatively because of the small size of the catchment formed by a single incision and the corresponding hillslope.     

Assessment of overland flow variation and blue water production in a farmed semi-arid water harvesting catchment

Upgrading agriculture in semi-arid areas and ensuring its sustainability require an optimal management of rainfall partition between blue and green waters in the farmed water harvesting catchment. The main objective of this study is to analyze the influence of heterogeneous land use on the spatial and temporal variation of rainfall partitioning and blue water production within a typical farmed catchment located in north-eastern Tunisia. The catchment has an area of 2.6 km² and comprises at its outlet a dam, which retains the runoff water in a reservoir. Overland flow and soil water balance components were monitored during two cropping seasons (2000/2001 and 2001/2002) on a network of eleven plots of 2 m² each with different land use and soil characteristics. The hydrological balances of both the catchment and reservoir have been monitored since 1994.Observed data showed a very large temporal and spatial variability of overland flow within the catchment reflecting the great importance of total rainfall as well as land use. During the 2001/2002 season the results showed a large variation of the number of observed runoff events, from 27 to 39, and of the annual overland flow depths, from 8 mm (under vineyard on calcaric cambisols) up to 43 mm (under shrubs-pasture on haplic regosols), between the plots. The annual runoff amounts were moderate; they always corresponded to less than 15% of the annual rainfall amount whatever the observation scale. It was also observed that changes in land use in years with similar rainfall could lead to significant differences in blue water flow. An attempt for predicting the overland flow by the general linear regression approach showed an r² of 31%, the predictors used are the class of soil infiltration capacity, the initial moisture saturation ratio of the soil surface layer and the total rainfall amounts.These experimental results indicate that the variation in land use in a semi-arid catchment is a main factor of variation in soil surface conditions and explain the major role played by the former on hydrological behavior of the upstream area and on rainfall partition between overland flow and infiltration. Therefore, to predict the water harvesting capacities in terms of blue water production of a farmed catchment in semi-arid areas it seems essential to consider precisely its land use and its temporal evolution related to management practices.     

Distributed and integrated response of a geographic information system‐based hydrologic model in the eastern Palouse region,Idaho

Verification of distributed hydrologic models is rare owing to the lack of spatially detailed field measurements and a common mismatch between the scale at which soil hydraulic properties are measured and the scale of a single modelling unit. In this study, two of the most commonly calibrated parameters, i.e. soil depth and the vertical distribution of lateral saturated hydraulic conductivity K_s, were eliminated by a spatially detailed soil characterization and results of a hillslope‐scale field experiment. The soil moisture routing (SMR) model, a geographic information system‐based hydrologic model, was modified to represent the dominant hydrologic processes for the Palouse region of northern Idaho. Modifications included K_s as a double exponential function of depth in a single soil layer, a snow accumulation and melt algorithm, and a simple relationship between storage and perched water depth (PWD) using the drainable porosity. The model was applied to a 2 ha catchment without calibration to measured data. Distributed responses were compared with observed PWD over a 3‐year period on a 10 m × 15 m grid. Integrated responses were compared with observed surface runoff at the catchment outlet. The modified SMR model simulated the PWD fluctuations remarkably well, especially considering the shallow soils in this catchment: a 0·20 m error in PWD is equivalent to only a 1·6% error in predicted soil moisture content. Simulations also captured PWD fluctuations during a year with high spatial variability of snow accumulation and snowmelt rates at upslope, mid‐slope, and toe slope positions with errors as low as 0·09 m, 0·12 m, and 0·12 m respectively. Errors in distributed and integrated model simulations were attributed mostly to misrepresentation of rain events and snowmelt timing problems. In one location in the catchment, simulated PWD was consistently greater than observed PWD, indicating a localized recharge zone, which was not identified by the soil morphological survey. Copyright © 2006 John Wiley & Sons, Ltd.     

Soil moisture variation in relation to topography and land use in a hillslope catchment of the Loess Plateau, China

The profile characteristics and the temporal dynamics of soil moisture variation were studied at 26 locations in Da Nangou catchment (3.5 km²) in the loess area of China. Soil moisture measurements were performed biweekly at five depths in the soil profile (0–5, 10–15, 20–25, 40–45 and 70–75 cm) from May to October 1998 using Delta-T theta probe. Soil moisture profile type and temporal variation type and their relationship to topography and land use were identified by detrended canonical correspondence analysis (DCCA) and correlation analysis. The profile distribution of time-averaged soil moisture content can be classified into three types i.e. decreasing-type, waving-type and increasing-type. The profile features of soil moisture (e.g. profile gradient and profile variability) are influenced by different environmental factors. The profile type of soil moisture is only attributed to land use while profile gradient and profile variability of soil moisture is mainly related to land use and topography (e.g. landform type and slope). The temporal dynamics of layer-averaged soil moisture content is grouped into three types including three-peak type, synchro-four-peak type and lagged-four-peak type. These types are controlled by topography rather than by land use. The temporal dynamic type of soil moisture shows significant correlation with relative elevation, slope, aspect, while temporal variance displays significant relation with slope shape. The mean soil moisture is related to both the profile and dynamics features of soil moisture and is controlled by both land use and topography (e.g. aspect, position, slope and relative elevation). The spatial variability of soil moisture across landscape varies with both soil depths and temporal evolution.     

Mechanisms responsible for streamflow generation on a small,salt-affected and deeply weathered hillslope

This paper considers the contributions of overland flow, throughflow and deep seepage to the generation of streamflow in a salt-affected, deeply weathered landscape. Runoff mechanisms on a small hillslope in south-western Australia were dependent on the extent and development of variable source areas. In winter, streamflow generation was controlled by returnflow, saturation overland flow and throughflow. In summer, post-ponding, infiltration-excess and saturation overland flow dominated. The extent of the variable source area and the magnitude of streamflow were due to antecedent soil moisture, rainfall and slope morphology. Concave hillslope sections accumulated soil moisture due to both saturated and unsaturated lateral flow processes. Throughflow provided the mechanism and vehicle for solute movement from the groundwater discharge area to the stream. However, discharge from the deep aquifer was the primary mechanism responsible for soil salinity and maintaining the core of the variable source area. Estimates of throughflow which only take account of soil-water movement and disregard returnflow, will underestimate the magnitude of throughflow.     

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.     

Simulation of syn-eruptive floods in the circumvesuvian plain (southern Italy)

During explosive eruptions the deposition of fine-grained volcanic ash fallout reduces soil permeability, favouring runoff of meteoric water and thus increasing the occurrence of catastrophic floods. A fully dynamic, two-dimensional model was used to simulate flooding scenarios in the Vesuvian area following an explosive volcanic eruption. The highest risk occurs in the catchment area of the Acerra-Nola Plain N and NE of Vesuvius. This plain has a population of 70,000 living in low-lying areas. This catchment area is vulnerable to ash fall because it lies downwind of the dominant synoptic circulation and it lacks a natural outflow toward the sea. Our numerical simulations predict dangerous scenarios, even in quiescent periods, during extreme rain events (return periods of 200 years have been considered), and a significant increase in the extent of the flooded areas due to renewed volcanic activity. Based on these simulations a hazard zonation has been proposed. Editorial responsibility: A Woods     

Relationship between tree height and landslide characteristics obtained by GIS assessment

Landslides in forested landscapes have far-reaching implications, beyond that of just destroying the forest itself, sometimes initiating large-scale sediment disasters. Although vegetation increases slope stability through its root network, it is hard to evaluate its contribution to slope stability over a wide area. In this study, the relationship between tree height and landslide characteristics in the Ikawa catchment, central Japan, was investigated to develop a method for evaluating the effects of forest cover on slope stability over a regional extent. Catchment-wide tree height was obtained using airborne LiDAR point cloud data and used in conjunction with the root depth profile, measured for trees of various height by digging trenches. Root tensile strength per unit area of soil was calculated from individual root diameters and empirical power law equations on the relationship between root diameter and root tensile force in order to better understand the effect that tree height has on slope stability. Landslide density in the Ikawa catchment shows that landslides occur more frequently in forests with shorter trees, with occurrence decreasing as tree height increases. This is likely due to the stabilizing features of larger trees having a greater network of roots, which is supported by the general increase in total root area and the deeper penetration of root biomass into the soil as the height of trees surveyed increases. Landslide density was not solely affected by tree height, but also by slope gradient and plane curvature. Decreasing landslide occurrence and landslide area as tree height increases suggests that slope stability increases with tree height, while the random distribution of results when comparing landslide depth to tree height suggests that while tree height has an impact on relative slope stability, the landslide failure depth is independent of tree height, and thus controlled by other factors. © 2020 John Wiley & Sons, Ltd.     

Hillslope‐scale probabilistic characterization of soil moisture dynamics and average water balance

Probabilistic water balance modelling provides a useful framework for investigating the interactions between soil, vegetation, and the atmosphere. It has been used to estimate temporal soil moisture dynamics and ecohydrological responses at a point. This study combines a nonlinear rainfall–runoff theory with probabilistic water balance model to represent varied source area runoff as a function of rainfall depth and a runoff coefficient at hillslope scale. Analytical solutions of the soil‐moisture probability density function and average water balance model are then developed. Based on a sensitivity analysis of soil moisture dynamics, we show that when varied source area runoff is incorporated, mean soil moisture is always lower and total runoff higher, compared with the original probabilistic water balance model. The increased runoff from the inclusion of varied source area runoff is mainly because of a reduction in leakage when the index of dryness is less than one and evapotranspiration when the index of dryness is greater than one. Inclusion of varied source area runoff in the model means that the actual evapotranspiration is limited by less available water (i.e. water limit), which is stricter than Budyko’s and Milly’s water limit. Application of the model to a catchment located in Western Australia showed that the method can predict annual value of actual evapotranspiration and streamflow accurately. Copyright © 2012 John Wiley & Sons, Ltd.