A low‐dimensional physically based model of hydrologic control of shallow landsliding on complex hillslopes

Hillslopes have complex three‐dimensional shapes that are characterized by their plan shape, profile curvature of surface and bedrock, and soil depth. To investigate the stability of complex hillslopes (with different slope curvatures and plan shapes), we combine the hillslope‐storage Boussinesq (HSB) model with the infinite slope stability method. The HSB model is based on the continuity and Darcy equations expressed in terms of storage along the hillslope. Solutions of the HSB equation account explicitly for plan shape by introducing the hillslope width function and for profile curvature through the bedrock slope angle and the hillslope soil depth function. The presented model is composed of three parts: a topography model conceptualizing three‐dimensional soil mantled landscapes, a dynamic hydrology model for shallow subsurface flow and water table depth (HSB model) and an infinite slope stability method based on the Mohr–Coulomb failure law. The resulting hillslope‐storage Boussinesq stability model (HSB‐SM) is able to simulate rain‐induced shallow landsliding on hillslopes with non‐constant bedrock slope and non‐parallel plan shape. We apply the model to nine characteristic hillslope types with three different profile curvatures (concave, straight, convex) and three different plan shapes (convergent, parallel, divergent). In the presented model, the unsaturated storage has been calculated based on the unit head gradient assumption. To relax this assumption and to investigate the effect of neglecting the variations of unsaturated storage on the assessment of slope stability in the transient case, we also combine a coupled model of saturated and unsaturated storage and the infinite slope stability method. The results show that the variations of the unsaturated zone storage do not play a critical role in hillslope stability. Therefore, it can be concluded that the presented dynamic slope stability model (HSB‐SM) can be used safely for slope stability analysis on complex hillslopes. Our results show that after a certain period of rainfall the convergent hillslopes with concave and straight profiles become unstable more quickly than others, whilst divergent convex hillslopes remain stable (even after intense rainfall). In addition, the relation between subsurface flow and hillslope stability has been investigated. Our analyses show that the minimum safety factor (FS) occurs when the rate of subsurface flow is a maximum. In fact, by increasing the subsurface flow, stability decreases for all hillslope shapes. Copyright © 2008 John Wiley & Sons, Ltd.     

Effects of hillslope topography on hydrological responses in a weathered granite mountain,Japan: comparison of the runoff response between the valley‐head and the side slope

To evaluate the effects of hillslope topography on storm runoff in a weathered granite mountain, discharge rate, soil pore water pressures, and water chemistry were observed on two types of hillslope: a valley‐head (a concave hillslope) and a side slope (a planar hillslope). Hydrological responses on the valley‐head and side slope reflected their respective topographic characteristics and varied with the rainfall magnitude. During small rainfall events (<35 mm), runoff from the side slope occurred rapidly relative to the valley‐head. The valley‐head showed little response in storm runoff. As rainfall amounts increased (35–60 mm), the valley‐head yielded a higher flow relative to the side slope. For large rainfall events (>60 mm), runoff from both hillslopes increased with rainfall, although that from the valley‐head was larger than that from the side slope. The differences in the runoff responses were caused by differences in the roles of lower‐slope soils and the convergence of the hillslope. During small rainfall events, the side slope could store little water; in contrast, all rainwater could be stored in the soils at the valley‐head hollow. As the amount of rainfall increased, the subsurface saturated area of the valley‐head extended from the bottom to the upper portion of the slope, with the contributions of transient groundwater via lateral preferential flowpaths due to the high concentration of subsurface water. Conversely, saturated subsurface flow did not contribute to runoff responses, and the subsurface saturated area at the side slope did not extend to the upper slope for the same storm size. During large rainfall events, expansion of the subsurface saturated area was observed in both hillslopes. Thus, differences in the concentration of subsurface water, reflecting hillslope topography, may create differences in the extension of the subsurface saturated area, as well as variability in runoff responses. Copyright © 2007 John Wiley & Sons, Ltd.     

Rainfall characteristics for shallow landsliding in Seattle,Washington, USA

Shallow landsliding in the Seattle, Washington, area, has caused the occasional loss of human life and millions of dollars in damage to property. The effective management of the hazard requires an understanding of the rainfall conditions that result in landslides. We present an empirical approach to quantify the antecedent moisture conditions and rainstorm intensity and duration that have triggered shallow landsliding using 25 years of hourly rainfall data and a complementary record of landslide occurrence. Our approach combines a simple water balance to estimate the antecedent moisture conditions of hillslope materials and a rainfall intensity–duration threshold to identify periods when shallow landsliding can be expected. The water balance is calibrated with field‐monitoring data and combined with the rainfall intensity–duration threshold using a decision tree. Results are cast in terms of a hypothetical landslide warning system. Two widespread landslide events are correctly identified by the warning scheme; however, it is less accurate for more isolated landsliding. Copyright © 2005 John Wiley & Sons, Ltd.     

Hillslopes in humid-tropical climates aren't always wet: Implications for hydrologic response and landslide initiation in Puerto Rico

The devastating impacts of the widespread flooding and landsliding in Puerto Rico following the September 2017 landfall of Hurricane Maria highlight the increasingly extreme atmospheric disturbances and enhanced hazard potential in mountainous humid-tropical climate zones. Long-standing conceptual models for hydrologically driven hazards in Puerto Rico posit that hillslope soils remain wet throughout the year, and therefore, that antecedent soil wetness imposes a negligible effect on hazard potential. Our post-Maria in situ hillslope hydrologic observations, however, indicate that while some slopes remain wet throughout the year, others exhibit appreciable seasonal and intra-storm subsurface drainage. Therefore, we evaluated the performance of hydro-meteorological (soil wetness and rainfall) versus intensity-duration (rainfall only) hillslope hydrologic response thresholds that identify the onset of positive pore-water pressure, a predisposing factor for widespread slope instability in this region. Our analyses also consider the role of soil-water storage and infiltration rates on runoff generation, which are relevant factors for flooding hazards. We found that the hydro-meteorological thresholds outperformed intensity-duration thresholds for a seasonally wet, coarse-grained soil, although they did not outperform intensity-duration thresholds for a perennially wet, fine-grained soil. These end-member soils types may also produce radically different stormflow responses, with subsurface flow being more common for the coarse-grained soils underlain by intrusive rocks versus infiltration excess and/or saturation excess for the fine-grained soils underlain by volcaniclastic rocks. We conclude that variability in soil-hydraulic properties, as opposed to climate zone, is the dominant factor that controls runoff generation mechanisms and modulates the relative importance of antecedent soil wetness for our hillslope hydrologic response thresholds.     

Coupling between landslides and eroding stream channels reconstructed from spruce tree rings (examples from the Carpathians and Sudetes – Central Europe)

The analysis of the positive feedback between landslides and erosion requires determination of the precise temporal and spatial relations between events of colluvium delivery and fluvial erosion. In our study we use decennial datasets on the occurrence of landsliding and erosion achieved through dendrochronological methods. Four sites covering areas of landslide slopes and adjacent valley floors with stream channels were studied. Landsliding on slopes was dated from the tree‐ring eccentricity developed in stems tilted due to bedrock instability. Erosion in channels was dated using the wood anatomy of roots exposed by erosion of the soil cover. Analysis of the temporal relations between dated landsliding, erosion and precipitation record has revealed that two types of repeating sequences can be observed: (1) rainfall → landsliding → erosion; (2) rainfall → erosion → landsliding. These sequences are an indication of the occurrence of slope‐channel positive feedback in the sites studied. In the first type, landsliding triggered by rainfall delivers colluvia into the valley floor and causes its narrowing, which in turn causes increased erosion. In the second type erosion triggered by rainfall disturbs the slope equilibrium and causes landsliding. Landsliding and erosion, once triggered by precipitation, can occur alternately in years with average precipitation and reinforce one another. Bidirectional coupling between landsliding and channel erosion was shown notably through the effects of channel shifting and forced sinuosity and by increased erosion of the slopes opposite the active landslides. Observations also suggest that the repetition of sequences described over longer periods of time can lead to a general widening of the valley floor at the expense of slopes and to a gradual change of the valley cross‐profile from narrow, V‐shaped into a wide flat‐bottomed. Thus landsliding–erosion coupling/positive feedback was recognized as an important factor shaping hillslope–valley topography of the mid‐mountain areas studied. Copyright © 2014 John Wiley & Sons, Ltd.     

Post-deforestation soil loss from steepland hillslopes in Taranaki,New Zealand

Soil erosion on steepland hillslopes in Taranaki, New Zealand, where landsliding is the dominant erosion form, was investigated by comparing mean regolith depths between first-order basins that have had their forest cover removed for different periods of time. Regolith depth and slope angle data were collected along 19 profile lines and 30 profile lines from steepland basins that had been deforested for 10 and 85 years, respectively. These profile lines were subdivided into a total of 236 profile segments of relatively linear slope angle and uniform regolith depth, that averaged 17·5 m in length. The depth of pre-existing regolith on post-deforestation landslide sites is estimated from a regression of regolith depth on slope angle for undisturbed (non-landslide) profile segments. Regolith depletion on landslide sites is in turn estimated by subtracting the depth of regolith on landslide sites from the estimate of pre-existing regolith depth. Regolith depletion by post-deforestation landslides, averaged over the entire length of profile lines, gives an estimate of average surface lowering. For the area deforested for 85 years, average surface lowering by post-deforestation landslides is 0·15 ± 0·04 m, and is the same as the difference in mean depth of 0·15 ± 0·11 m between this area and the area deforested for 10 years. Erosion of regolith from hillslopes by processes other than landsliding appears to be minimal. The 0·15 m average surface lowering represents a regolith depletion rate of 1·8 ± 0±5 mm yr^?1. For hillslopes steeper than 28°, where all post-deforestation landslides occur, average surface lowering is 0·20 ± 0·05 m, and the regolith depletion rate is 2±4 · 0±6 mm yr^?1. Average surface lowering is greatest at 0·23 ± 0·07 m on hillslopes steeper than 32° where most post-deforestation landslides occur. Here, the regolith depletion rate is 2·7 ± 0·8 mm yr^?1. A large-magnitude, low-frequency storm in March 1990, produced an average surface lowering of 0·041 m. There were proportionately more landslides in the area deforested for 10 years, illustrating the importance of previous erosion history of hillslopes on the spatial distribution of landslides. There were also comparatively few landslides on steeper hillslopes because previous lower magnitude storms had already removed much of the deeper regolith.     

Shallow groundwater response to rainfall on a forested headwater catchment in northern coastal California: implications of topography,rainfall, and throughfall intensities on peak pressure head generation

The pore water pressure head that builds in the soil during storms is a critical factor for the prediction of potential slope instability. We report findings from a 3‐year study of pressure head in 83 piezometers distributed within a 13‐ha forested catchment on the northern coast of California. The study's primary objective was to observe the seasonal and storm‐based dynamics of pressure head at a catchment scale in relation to observed rainfall characteristics and in situ topography to better understand landscape patterns of pressure head. An additional goal was to determine the influence of the interaction between rainfall and forest canopy in altering delivery of water and pressure head during the large storms necessary to induce landsliding. We found that pressure head was highly variable in space and time at the catchment scale. Pore pressures peaked close to maximum rainfall intensity during the largest storms measured. The difference between rainfall and throughfall delivered through the canopy was negligible during the critical landslide‐producing peak rainfall periods. Pore pressure was spatially variable within the catchment and did not strongly correlate with surficial topographic features. Only 23% of the piezometers located in a variety of slope positions were found to be highly responsive to rainfall. Topographic index statistically explained peak pressure head at responsive locations during common storms, but not during the larger storms with potential to produce landslides. Drainage efficiency throughout the catchment increased significantly in storms exceeding 2 to 7 months peak pressure head return period indicated by slowing or cessation of the rate of increase of pressure head with increasing storm magnitude. This asymptotic piezometric pattern persisted through the largest storm measured during the study. Faster soil drainage suppressed pressure head response in larger storms with important process implications for pore pressure development and landslide hazard modelling. Copyright © 2012 John Wiley & Sons, Ltd.     

Effects of hillslope topography on runoff response in a small catchment in the Fudoji Experimental Watershed,central Japan

We investigated the role of different hillslope units with different topographic characteristics on runoff generation processes based on field observations at two types of hillslopes (0·1 ha): a valley‐head (a convergent hillslope) and a side slope (a planar hillslope), as well as at three small catchments having two types of slopes with different drainage areas ranging from 1·9 to 49·7 ha in the Tanakami Mountains, central Japan. We found that the contribution of the hillslope unit type to small catchment runoff varied with the magnitude of rainfall. When the total amount of rainfall for a single storm event was < 35 mm, runoff in the small catchment was predominantly generated from the side slope. As the amount of rainfall increased (>35 mm), the valley‐head also began to contribute to the catchment runoff, adding to runoff from the side slope. Although the direct runoff from the valley‐head was greater than that from the side slope, the contribution from the side slope was quantitatively greater than that from the valley‐head due to the proportionally larger area occupied by the side slope in the small catchment. The storm runoff responses of the small catchments reflected the change in the runoff components of each hillslope unit as the amount of rainfall increased and rainfall patterns changed. However, similar runoff responses were found for the small catchments with different areas. The similarity of the runoff responses is attributable to overlay effects of different hillslope units and the similar composition ratios of the valley‐head and side slope in the catchments. This study suggests that the relative roles of the valley‐head and side slope are important in runoff generation and solute transport as the catchment size increases from a hillslope/headwater to a small catchment. Copyright © 2011 John Wiley & Sons, Ltd.     

The role of the seismic trigger in the Calitri landslide (Italy): historical reconstruction and dynamic analysis

The study covered by this paper was focused on the historical case of the Calitri landslide, which was repeatedly reactivated by earthquakes, as reported since 1694. The town of Calitri (Southern Italy) is located on a ridge whose southern slope, from its top to the Ofanto river valley floor, has been historically affected by major landsliding. The last record of recurrence of the Calitri landslide leads back to the 1980 Irpinia earthquake, which caused significant damage to the town and had pervasive and visible ground effects. Based on a detailed historical reconstruction of landsliding and seismicity at Calitri, the study analysed the current static and dynamic stability of the landsliding slope by means of a finite-difference numerical analysis, taking into account the various factors of landslide initiation and kinematics. The results of the numerical analysis: (i) were consistent with the roto-translational mechanism observed upon the latest reactivation of the landslide; (ii) demonstrated that excess pore pressure redistribution caused a lag between the seismic trigger and the initiation of landsliding; and (iii) showed the impact of seismic input frequency on propagation and depth of slope instability.     

Lithological controls on hillslope sediment supply: insights from landslide activity and grain size distributions

The volumes, rates and grain size distributions of sediment supplied from hillslopes represent the initial input of sediment delivered from upland areas and propagated through sediment routing systems. Moreover, hillslope sediment supply has a significant impact on landscape response time to tectonic and climatic perturbations. However, there are very few detailed field studies characterizing hillslope sediment supply as a function of lithology and delivery process. Here, we present new empirical data from tectonically‐active areas in southern Italy that quantifies how lithology and rock strength control the landslide fluxes and grain size distributions supplied from hillslopes. Landslides are the major source of hillslope sediment supply in this area, and our inventory of ~2800 landslides reveals that landslide sediment flux is dominated by small, shallow landslides. We find that lithology and rock strength modulate the abundance of steep slopes and landslides, and the distribution of landslide sizes. Outcrop‐scale rock strength also controls the grain sizes supplied by bedrock weathering, and influences the degree of coarsening of landslide supply with respect to weathering supply. Finally, we show that hillslope sediment supply largely determines the grain sizes of fluvial export, from catchments and that catchments with greater long‐term landslide rates deliver coarser material. Therefore, our results demonstrate a dual control of lithology on hillslope sediment supply, by modulating both the sediment fluxes from landslides and the grain sizes supplied by hillslopes to the fluvial system. Copyright © 2017 John Wiley & Sons, Ltd.     

A spatiotemporal probabilistic modelling of storm‐induced shallow landsliding using aerial photographs and logistic regression

Landslides constitute one of the major natural hazards that could cause significant losses of life and property. Mapping or delineating areas prone to landsliding is therefore essential for land‐use activities and management decision making in hilly or mountainous regions. A landslide hazard map can be constructed by a qualitative combination of maps of site conditions, including geology, topography and geomorphology, by statistical methods through correlating landslide occurrence with geologic and geomorphic factors, or by using safety factors from stability analysis. A landslide hazard map should provide information on both the spatial and temporal probabilities of landsliding in a certain area. However, most previous studies have focused on susceptibility mapping, rather than on hazard mapping in a spatiotemporal context. This study aims at developing a predictive model, based on both quasi‐static and dynamic variables, to determine the probability of landsliding in terms of space and time. The study area selected is about 13 km² in North Lantau, Hong Kong. The source areas of the landslides caused by the rainstorms of 18 July 1992 and 4–5 November 1993 were interpreted from multi‐temporal aerial photographs. Landslide data, lithology, digital elevation model data, land cover, and rainfall data were digitized into a geographic information system database. A logistic regression model was developed using lithology, slope gradient, slope aspect, elevation, slope shape, land cover, and rolling 24 h rainfall as independent variables, since the dependent variable could be expressed in a dichotomous way. This model achieved an overall accuracy of 87·2%, with 89·5% of landslide grid cells correctly classified and found to be performing satisfactorily. The model was then applied to rainfalls of a variety of periods of return, to predict the probability of landsliding on natural slopes in space and time. It is observed that the modelling techniques described here are useful for predicting the spatiotemporal probability of landsliding and can be used by land‐use planners to develop effective management strategies. Copyright © 2003 John Wiley & Sons, Ltd.     

Parametric study of a physically-based,plot-scale hillslope hydrological model through virtual experiments

Abstract

A physically-based hillslope hydrological model with shallow overland flow and rapid subsurface stormflow components was developed and calibrated using field experiments conducted on a preferential path nested hillslope in northeast India. Virtual experiments were carried out to perform sensitivity analysis of the model using the automated parameter estimation (PEST) algorithm. Different physical parameters of the model were varied to study the resulting effects on overland flow and subsurface stormflow responses from the theoretical hillslopes. It was observed that topographical shapes had significant effects on overland flow hydrographs. The slope profiles, surface storage, relief, rainfall intensity and infiltration rates primarily controlled the overland flow response of the hillslopes. Prompt subsurface stormflow responses were mainly dominated by lateral preferential flow, as soil matrix flow rates were very slow. Rainfall intensity and soil macropore structures were the most influential parameters on subsurface stormflow. The number of connected soil macropores was a more sensitive parameter than the size of macropores. In hillslopes with highly active vertical and lateral preferential pathways, saturation excess overland flow was not evident. However, saturation excess overland flow was generated if the lateral macropores were disconnected. Under such conditions, rainfall intensity, duration and preferential flow rate governed the process of saturation excess overland flow generation from hillslopes.

Editor D. Koutsoyiannis; Associate editor C. Perrin     

The design of post‐mining landscapes using geomorphic principles

Nature can provide analogues for post‐mining landscapes in terms of landscape stability and also in terms of the rehabilitated structure ‘blending in’ with the surrounding undisturbed landscape. In soil‐mantled landscapes, hillslopes typically have a characteristic pro?le that has a convex upper hillslope pro?le with a concave pro?le lower down the slope. In this paper hillslope characteristic form is derived using the area–slope relationship from pre‐mining topography at two sites in Western Australia. Using this relationship, concave hillslope pro?les are constructed and compared to linear hillslopes in terms of sediment loss using the SIBERIA erosion model. It is found that concave hillslopes can reduce sediment loss by up to ?ve times that of linear slopes. Concave slopes can therefore provide an alternative method for the construction of post‐mining landscapes. An understanding of landscape geomorphological properties and the use of erosion models can greatly assist in the design of post‐mining landscapes. Copyright © 2003 John Wiley & Sons, Ltd.     

Applying Probability Determination to Refine Landslide-triggering Rainfall Thresholds Using an Empirical “Antecedent Daily Rainfall Model”

—Rainfall-triggered landslides constitute a serious hazard and an important geomorphic process in many parts of the world. Attempts have been made at various scales in a number of countries to investigate triggering conditions in order to identify patterns in behaviour and, ultimately, to define or calculate landslide-triggering rainfall thresholds. This study was carried out in three landslide-prone regions in the North Island of New Zealand. Regional landslide-triggering rainfall thresholds were calculated using an empirical “Antecedent Daily Rainfall Model.” In this model, first introduced by, triggering rainfall conditions are represented by a combination of rainfall occurring in a period before the event (antecedent rainfall) and rainfall on the day of the event. A physically-based decay coefficient is derived for each region from the recessional behaviour of storm hydrographs and is used to produce an index for antecedent rainfall. Statistical techniques are employed to obtain the thresholds which best separate the rainfall conditions associated with landslide occurrence from those of non-occurrence or a given probability of occurrence.The resultant regional models are able to represent the probability of occurrence of landsliding events on the basis of rainfall conditions. The calculated thresholds show regional differences in susceptibility of a given landscape to landslide-triggering rainfall. These differences relate to both the landslide database and the difference of existing physical conditions between the regions.     

Hillslope response to knickpoint migration in the Southern Appalachians: implications for the evolution of post‐orogenic landscapes

The southern Appalachians represent a landscape characterized by locally high topographic relief, steep slopes, and frequent mass movement in the absence of significant tectonic forcing for at least the last 200 Ma. The fundamental processes responsible for landscape evolution in a post‐orogenic landscape remain enigmatic. The non‐glaciated Cullasaja River basin of south‐western North Carolina, with uniform lithology, frequent debris flows, and the availability of high‐resolution airborne lidar DEMs, is an ideal natural setting to study landscape evolution in a post‐orogenic landscape through the lens of hillslope–channel coupling. This investigation is limited to channels with upslope contributing areas >2.7 km², a conservative estimate of the transition from fluvial to debris‐flow dominated channel processes. Values of normalized hypsometry, hypsometric integral, and mean slope vs elevation are used for 14 tributary basins and the Cullasaja basin as a whole to characterize landscape evolution following upstream knickpoint migration. Results highlight the existence of a transient spatial relationship between knickpoints present along the fluvial network of the Cullasaja basin and adjacent hillslopes. Metrics of topography (relief, slope gradient) and hillslope activity (landslide frequency) exhibit significant downstream increases below the current position of major knickpoints. The transient effect of knickpoint‐driven channel incision on basin hillslopes is captured by measuring the relief, mean slope steepness, and mass movement frequency of tributary basins and comparing these results with the distance from major knickpoints along the Cullasaja River. A conceptual model of area–elevation and slope distributions is presented that may be representative of post‐orogenic landscape evolution in analogous geologic settings. Importantly, the model explains how knickpoint migration and channel–hillslope coupling is an important factor in tectonically‐inactive (i.e. post‐orogenic) orogens for the maintenance of significant relief, steep slopes, and weathering‐limited hillslopes. Copyright © 2011 John Wiley & Sons, Ltd.     

Physically based modelling of sheet erosion (detachment and deposition processes) in complex hillslopes

A one‐dimensional uncoupled model governed by this research is a physics‐based modelling of the rainfall‐runoff induced erosion process. The presented model is composed of three parts of a three‐dimensional (3D) hillslope geometry, a nonlinear storage (kinematic wave) model for hillslope hydrological response, and an unsteady physically based surface erosion model. The 3D hillslope geometry model allows describing of the hillslope morphology by defining their plan shape and profile curvature. By changing these two topographic parameters, nine basic hillslope types are derived. The modelling of hillslope hydrological response is based on a flow continuity equation as the relation of discharge and flow depth is passed on kinematic wave approximation. The erosion model is based on a mass conservation equation for unsteady flow. The model assumes that suspended sediment does not affect flow dynamics. The model also accounts for the effect of flow depth plus loose soil depth on soil detachment. The presented model was run for two different precipitations, slope content, and length, and results were plotted for sediment detachment/deposition rate. Based on the obtained results, in hillslopes with convex and straight profile curvatures, sediment detachment only occurred in the whole length of the hillslope. However, in concave ones, sediment detachment and deposition only occurred together in hillslope. The hillslopes with straight profiles and convergent plans have the highest rate of detachment. Also, results show that most detachment rates occur in convex profile curvatures, which are about 15 times more than in straight profiles. Copyright © 2015 John Wiley & Sons, Ltd.     

Spatio‐temporal variability of piezometric response on two steep alpine hillslopes

Information on the main drivers of subsurface flow generation on hillslopes of alpine headwater catchments is still missing. Therefore, the dominant factors controlling the water table response to precipitation at the hillslope scale in the alpine Bridge Creek Catchment, Northern Italy, were investigated. Two steep hillslopes of similar size, soil properties and vegetation cover but contrasting topography were instrumented with 24 piezometric wells. Sixty‐three (63) rainfall‐runoff events were selected over three years in the snow‐free months to analyse the influence of rainfall depth, antecedent moisture conditions, hillslope topographic characteristics and soil depth on shallow water table dynamics. Piezometric response, expressed as percentage of well activation and water peak magnitude, was strongly correlated with soil moisture status, as described by an index combining antecedent soil moisture and rainfall depth. Hillslope topography was found to be a dominant control only for the convex‐divergent hillslope and during wet conditions. Timing of water table response depended primarily on soil depth and topographic position, with piezometric peak response occurring later and showing a greater temporal variability at the hillslope bottom, characterized by thicker soil. The relationship between mean hillslope water table level and standard deviation for all wells reflected the timing of the water table response at the different locations along the hillslopes. The outcomes of this research contribute to a better understanding of the controls on piezometric response at the hillslope scale in steep terrain and its role on the hydrological functioning of the study catchment and of other sites with similar physiographic characteristics. Copyright © 2014 John Wiley & Sons, Ltd.     

Dynamic domain kinematic modelling for predicting interflow over leaky impeding layers

Traditional Boussinesq or kinematic simulations of interflow (i.e., lateral subsurface flow) assume no leakage through the impeding layer and require a no-flow boundary condition at the ridge top. However, recent analyses of many interflow-producing landscapes indicate that leaky impeding layers are common, that most interflow percolates well before reaching the toe slope, and therefore, the downslope contributing length is shorter than the hillslope length. In watersheds characterised by perched interflow over a low conductivity layer through permeable topsoil, interflow with percolation may be modelled with a kinematic wave model using a mobile upslope boundary condition defining the hillslope portion contributing interflow to valleys. Here, we developed and applied a dynamic interflow model to simulate interflow using a downslope travel distance concept such that only the active contributing length is modelled at any time. The model defines a variable active area based on the depth of the perched layer, the topographic slope and the ratio of the hydraulic conductivity of topsoil to that of the impeding layer. It incorporates a two-layer soil moisture accounting water balance analysis, a pedo-transfer function, and percolation and evaporation routines to predict interflow rates in continuous and event-based scenarios. We tested the modelling concept on two sets of data (2-year dataset of rainfall observations for the continuous simulation and a multi-day irrigation experiment for the event simulation) from a 121-m-long open interflow collection trench on an experimental hillslope at the Savannah River Site, South Carolina. The continuous model simulation partially represented the observed interflow hydrograph and perched water depth in the experimental hillslope with correlation coefficients of 0.85 and 0.35, respectively. Model performance improved significantly at event-scale analysis. The modelling approach realistically represents interflow dynamics in hillslopes with leaky impeding layers and can be integrated into catchment-scale hydrology models for more detailed hillslope process modelling.     

How long is a hillslope?

Hillslope length is a fundamental attribute of landscapes, intrinsically linked to drainage density, landslide hazard, biogeochemical cycling and hillslope sediment transport. Existing methods to estimate catchment average hillslope lengths include inversion of drainage density or identification of a break in slope–area scaling, where the hillslope domain transitions into the fluvial domain. Here we implement a technique which models flow from point sources on hilltops across pixels in a digital elevation model (DEM), based on flow directions calculated using pixel aspect, until reaching the channel network, defined using recently developed channel extraction algorithms. Through comparisons between these measurement techniques, we show that estimating hillslope length from plots of topographic slope versus drainage area, or by inverting measures of drainage density, systematically underestimates hillslope length. In addition, hillslope lengths estimated by slope–area scaling breaks show large variations between catchments of similar morphology and area. We then use hillslope length–relief structure of landscapes to explore nature of sediment flux operating on a landscape. Distinct topographic forms are predicted for end‐member sediment flux laws which constrain sediment transport on hillslopes as being linearly or nonlinearly dependent on hillslope gradient. Because our method extracts hillslope profiles originating from every ridgetop pixel in a DEM, we show that the resulting population of hillslope length–relief measurements can be used to differentiate between linear and nonlinear sediment transport laws in soil mantled landscapes. We find that across a broad range of sites across the continental United States, topography is consistent with a sediment flux law in which transport is nonlinearly proportional to topographic gradient. © 2016 The Authors. Earth Surface Processes and Landforms published by John Wiley & Sons Ltd.