Soil moisture causes dynamic adjustments to root reinforcement that reduce slope stability

In steep soil‐mantled landscapes, the initiation of shallow landslides is strongly controlled by the distribution of vegetation, whose roots reinforce the soil. The magnitude of root reinforcement depends on the number, diameter distribution, orientation and the mechanical properties of roots that cross potential failure planes. Understanding how these properties vary in space and time in forests remains a significant challenge. Here we test the hypothesis that spatio‐temporal variations in root reinforcement along a hillslope occur as a function of topographic soil moisture gradients. To test this hypothesis we compared root reinforcement measurements from relatively dry, divergent noses to relatively wet, convergent hollows in the southern Appalachian Mountains, North Carolina, USA. Our initial results showed that root reinforcement decreased in areas of higher soil moisture because the tensile strength of roots decreased. A post hoc laboratory experiment further demonstrated that root tensile strength decreased as root moisture content increased. This effect is consistent with other experiments on stem woods showing that increased water content in the cell wall decreases tensile strength. Our experimental data demonstrated that roots can adjust to changes in the external root moisture conditions within hours, suggesting that root moisture content will change over the timescale of large storm events (hours–days). We assessed the effects of the dynamic changes in root tensile strength to the magnitude of apparent cohesion within the infinite slope stability model. Slopes can be considerably less stable when precipitation‐driven increases in saturated soil depth both increase pore pressures and decrease root reinforcement. Copyright © 2016 John Wiley & Sons, Ltd.     

Spatial and temporal analysis of hillslope–channel coupling and implications for the longitudinal profile in a dryland basin

The long‐term evolution of channel longitudinal profiles within drainage basins is partly determined by the relative balance of hillslope sediment supply to channels and the evacuation of channel sediment. However, the lack of theoretical understanding of the physical processes of hillslope–channel coupling makes it challenging to determine whether hillslope sediment supply or channel sediment evacuation dominates over different timescales and how this balance affects bed elevation locally along the longitudinal profile. In this paper, we develop a framework for inferring the relative dominance of hillslope sediment supply to the channel versus channel sediment evacuation, over a range of temporal and spatial scales. The framework combines distinct local flow distributions on hillslopes and in the channel with surface grain‐size distributions. We use these to compute local hydraulic stresses at various hillslope‐channel coupling locations within the Walnut Gulch Experimental Watershed (WGEW) in southeast Arizona, USA. These stresses are then assessed as a local net balance of geomorphic work between hillslopes and channel for a range of flow conditions generalizing decadal historical records. Our analysis reveals that, although the magnitude of hydraulic stress in the channel is consistently higher than that on hillslopes, the product of stress magnitude and frequency results in a close balance between hillslope supply and channel evacuation for the most frequent flows. Only at less frequent, high‐magnitude flows do channel hydraulic stresses exceed those on hillslopes, and channel evacuation dominates the net balance. This result suggests that WGEW exists mostly (~50% of the time) in an equilibrium condition of sediment balance between hillslopes and channels, which helps to explain the observed straight longitudinal profile. We illustrate how this balance can be upset by climate changes that differentially affect relative flow regimes on slopes and in channels. Such changes can push the long profile into a convex or concave condition. © 2018 The Authors. Earth Surface Processes and Landforms published by John Wiley & Sons Ltd.     

Effects of gradient,distance, curvature and aspect on steep burned and unburned hillslope soil erosion and deposition

Wildfire denudes vegetation and impacts chemical and physical soil properties, which can alter hillslope erosion rates. Post‐wildfire erosion can also contribute disproportionately to long‐term erosion rates and landscape evolution. Post‐fire hillslope erosion rates remain difficult to predict and document at the hillslope scale. Here we use ²¹⁰Pb_aex (lead‐210 mineral‐adsorbed excess) inventories to describe net sediment erosion on steep, convex hillslopes in three basins (unburned, moderately and severely burned) in mountainous central Idaho. We analyzed nearly 300 soil samples for ²¹⁰Pb_aex content with alpha spectrometry and related net sediment erosion to burn severity, aspect, gradient, curvature and distance from ridgetop. We also tested our data against models for advective, linear and non‐linear diffusive erosion. Statistically lower net soil losses on north‐ versus south‐facing unburned hillslopes suggest that greater vegetative cover and soil cohesion on north‐facing slopes decrease erosion. On burned hillslopes, erosion differences between aspects were less apparent and net erosion was more variable, indicating that vegetation influences erosion magnitude and fire drives erosion variability. We estimated net soil losses throughout the length of unburned hillslopes, including through a footslope transition to concave form. In contrast, on burned hillslopes, the subtle shift from convex to concave form was associated with deposition of a post‐fire erosion pulse. Such overall patterns of erosion and deposition are consistent with predictions from a non‐linear diffusion equation. This finding also suggests that concave sections of overall convex hillslopes affect post‐disturbance soil erosion and deposition. Despite these patterns, no strong relationships were evident between local net soil losses and gradient, curvature, distance from ridgetop, or erosion predicted with advection or diffusion equations. The observed relationship between gradient and erosion is therefore likely more complex or stochastic than often described theoretically, especially over relatively short timescales (60–100 years). Copyright © 2016 John Wiley & Sons, Ltd.     

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.     

Without power? Landslide inventories in the face of climate change

Projected scenarios of climate change involve general predictions about the likely changes to the magnitude and frequency of landslides, particularly as a consequence of altered precipitation and temperature regimes. Whether such landslide response to contemporary or past climate change may be captured in differing scaling statistics of landslide size distributions and the erosion rates derived thereof remains debated. We test this notion with simple Monte Carlo and bootstrap simulations of statistical models commonly used to characterize empirical landslide size distributions. Our results show that significant changes to total volumes contained in such inventories may be masked by statistically indistinguishable scaling parameters, critically depending on, among others, the size of the largest of landslides recorded. Conversely, comparable model parameter values may obscure significant, i.e. more than twofold, changes to landslide occurrence, and thus inferred rates of hillslope denudation and sediment delivery to drainage networks. A time series of some of Earth's largest mass movements reveals clustering near and partly before the last glacial‐interglacial transition and a distinct step‐over from white noise to temporal clustering around this period. However, elucidating whether this is a distinct signal of first‐order climate‐change impact on slope stability or simply coincides with a transition from short‐term statistical noise to long‐term steady‐state conditions remains an important research challenge. Copyright © 2011 John Wiley & Sons, Ltd.     

Effects of spatially structured vegetation patterns on hillslope erosion in a semiarid Mediterranean environment: a simulation study

A general trend of decreasing soil loss rates with increasing vegetation cover fraction is widely accepted. Field observations and experimental work, however, show that the form of the cover‐erosion function can vary considerably, in particular for low cover conditions that prevail on arid and semiarid hillslopes. In this paper the structured spatial distribution of the vegetation cover and associated soil attributes is proposed as one of the possible causes of variation in cover–erosion relationships, in particular in dryland environments where patchy vegetation covers are common. A simulation approach was used to test the hypothesis that hillslope discharge and soil loss could be affected by variation in the spatial correlation structure of coupled vegetation cover and soil patterns alone. The Limburg Soil Erosion Model (LISEM) was parameterized and verified for a small catchment with discontinuous vegetation cover at Rambla Honda, SE Spain. Using the same parameter sets LISEM was subsequently used to simulate water and sediment fluxes on 1 ha hypothetical hillslopes with simulated spatial distributions of vegetation and soil parameters. Storms of constant rainfall intensity in the range of 30–70 mm h^?1 and 10–30 min duration were applied. To quantify the effect of the spatial correlation structure of the vegetation and soil patterns, predicted discharge and soil loss rates from hillslopes with spatially structured distributions of vegetation and soil parameters were compared with those from hillslopes with spatially uniform distributions. The results showed that the spatial organization of bare and vegetated surfaces alone can have a substantial impact on predicted storm discharge and erosion. In general, water and sediment yields from hillslopes with spatially structured distributions of vegetation and soil parameters were greater than from identical hillslopes with spatially uniform distributions. Within a storm the effect of spatially structured vegetation and soil patterns was observed to be highly dynamic, and to depend on rainfall intensity and slope gradient. Copyright © 2005 John Wiley & Sons, Ltd.     

Changes in hillslope hydrology in a perched,shallow soil system due to clearcutting and residual biomass removal

Sustainable fuels legislation and volatility in energy prices have put additional pressures on the forestry sector to intensify the harvesting of biomass for “advanced biofuel” production. To better understand how residual biomass removal after harvest affects forest hydrology in relatively low slope terrain, a Before-After-Control-Impact (BACI) study was conducted in the USDA Forest Service's Marcell Experimental Forest, Minnesota, USA. Hydrological measurements were made from 2010–2013 on a forested hillslope that was divided into three treatment blocks, where one block was harvested and residual biomass removed (Biomass Removed), the second was harvested and residual biomass left (Biomass Left), and the last block was left as an Unharvested Control. The pre-harvest period (2 years) was 2010–11 and post-harvest (2 years) was 2012–13. Water table elevation at the upslope and downslope position, subsurface runoff, and soil moisture were measured between May–November. Mixed effect statistical models were used to compare both the before-after and “control” treatment ratios (ratios between harvested hillslopes and the Unharvested Control hillslope). Subsurface runoff significantly increased (p < .05) at both harvested hillslopes but to a greater degree on the Biomass Left hillslope. Greater subsurface runoff volumes at both harvested hillslopes were driven by substantial increases during fall, with additional significant increases during summer on the Biomass Left hillslope. The hydrological connectivity, inferred from event runoff ratios, increased due to harvesting at both hillslopes but only significantly on the Biomass Left hillslope. The winter harvest minimized soil disturbance, resulting in no change to the effective hydraulic conductivity distribution with depth. Thus, the observed hydrological changes were driven by increased effective precipitation and decreased evapotranspiration, increasing the duration that both harvested hillslopes were hydrologically active. The harvesting of residual biomass appears to lessen hydrological connectivity relative to leaving residual biomass on the hillslope, potentially decreasing downstream hydrological impacts of similar forestry operations.     

Interactions and connectivity between runoff generation processes of different spatial scales

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

Including root reinforcement variability in a probabilistic 3D stability model

Forests play a significant role in protecting people, settlements in mountainous terrains from hydrogeomorphic hazards, including shallow landslides. Although several studies have investigated the interactions between forests and slope instabilities, a full understanding of them has not yet been obtained. Additionally, models that incorporate forest stand properties into slope failure probability analyses have not been developed. In principle, physical‐based models, which are powerful tools for landslide hazard analyses, represent an appropriate approach to linking stand properties and slope stability. However, the reliability of these models depends on numerous parameters that describe highly complex geotechnical and hydrological processes (e.g. potential failure depth, saturation ratio, root reinforcement, etc.) that are difficult to measure and model. In particular, the spatial heterogeneity of root reinforcement remains a problem, and the use of physically based models from a forest management perspective has been limited. This paper presents a procedure for assessing slope stability in terms of the Factor of Safety that accounts for forest stand characteristics such as tree density, average diameter at breast height and minimum distance between trees. The procedure combines a three‐dimensional (3D) slope stability model with an evaluation of the variability of root reinforcement in terms of a probability distribution, according to forest characteristics. Monte Carlo simulation is used to account for the residual uncertainties in both stand characteristics and 3D stability model parameters. The proposed method was applied in a subalpine catchment in the Italian Alps, mainly covered by coniferous forest and characterized by steep slopes and high landslide risk. The results suggest that the procedure is highly reliable, according to landslide inventory maps [area under the ROC curve (AUC) is 0.82 and modified success rate (MSR) is 0.70]. Thus, it represents a promising tool for studying the role of root reinforcement in landslide hazard mapping and guiding forest management from a slope stability perspective. Copyright © 2017 John Wiley & Sons, Ltd.     

The effect of accelerated soil erosion on hillslope morphology

Intensive agricultural land use can have detrimental effects on landscape properties, greatly accelerating soil erosion, with consequent fertility loss and reduced agricultural potential. To quantify the effects of such erosional processes on hillslope morphology and gain insight into the underlying dynamics, we use a twofold approach. First, a statistical analysis of topographical features is conducted, with a focus on slope and gradient distributions. The accelerated soil erosion is shown to be fingerprinted in the distribution tails, which provide a clear statistical signature of this human-induced land modification. Theoretical solutions are then derived for the hillslope morphology and the associated creep and runoff erosion fluxes, allowing us to distinguish between the main erosional mechanisms operating in disturbed and undisturbed areas. We focus our application on the landscape at the Calhoun Critical Zone Observatory in the US Southern Piedmont, where severe soil erosion followed intensive cotton cultivation, resulting in highly eroded and gullied hillslopes. The observed differences in hillslope morphologies in disturbed and undisturbed areas are shown to be related to the disruption of the natural balance between soil creep and runoff erosion. The relaxation time required for the disturbed hillslopes to reach a quasi-equilibrium condition is also investigated. © 2019 John Wiley & Sons, Ltd.     

Downvalley fining of hillslope sediment in an alpine catchment: implications for downstream fining of sediment flux in mountain rivers

The size distributions of sediment delivered from hillslopes to rivers profoundly influence river morphodynamics, including river incision into bedrock and the quality of aquatic habitat. Yet little is known about the factors that influence size distributions of sediment produced by weathering on hillslopes. We present results of a field study of hillslope sediment size distributions at Inyo Creek, a steep catchment in granitic bedrock of the Sierra Nevada, USA. Particles sampled near the base of hillslopes, adjacent to the trunk stream, show a pronounced decrease in sediment size with decreasing sample elevation across all but the coarsest size classes. Measured size distributions become increasingly bimodal with decreasing elevation, exhibiting a coarse, bouldery mode that does not change with elevation and a more abundant finer mode that shifts from cobbles at the highest elevations to gravel at mid elevations and finally to sand at low elevations. We interpret these altitudinal variations in hillslope sediment size to reflect changes in physical, chemical, and biological weathering that can be explained by the catchment's strong altitudinal gradients in topography, climate, and vegetation cover. Because elevation and travel distance to the outlet are closely coupled, the altitudinal trends in sediment size produce a systematic decrease in sediment size along hillslopes parallel to the trunk stream. We refer to this phenomenon as ‘downvalley fining.’ Forward modeling shows that downvalley fining of hillslope sediment is necessary for downstream fining of the long-term average flux of coarse sediment in mountain landscapes where hillslopes and channels are coupled and long-term net sediment deposition is negligible. The model also shows that abrasion plays a secondary role in downstream fining of coarse sediment flux but plays a dominant role in partitioning between the bedload and suspended load. Patterns observed at Inyo Creek may be widespread in mountain ranges around the world. © 2020 The Authors. Earth Surface Processes and Landforms published by 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.     

The influence of climate and microclimate (aspect) on soil creep efficiency: Cinder cone morphology and evolution along the eastern Mediterranean Golan Heights

Although hillslope evolution has been subject to much investigation for more than a century, the effect of climate on the morphology of soil‐mantled hillslopes remains poorly constrained. In this study, we perform numerical simulations of volcanic cinder cones in the Golan Heights (eastern Mediterranean) to estimate soil transport efficiency across a significant north–south gradient in mean annual precipitation (1100 to 500 mm). We use the initial cinder cone morphology (constrained by stratigraphy), the modern hillslope form (surveyed with sub‐meter accuracy) and the eruption age (based on ⁴⁰Ar–³⁹Ar chronology) to predict the best‐fit value of the soil transport coefficient (‘diffusivity’) based on a nonlinear transport model. Our results indicate that the best‐fit diffusivity (K ) varies from 1 to 6 m² ka^?1 among the five cinder cones in our field area. Diffusivity (K ) values vary systematically with precipitation and hillslope aspect; specifically, K is higher on south‐facing (drier) hillslopes and decreases with mean annual precipitation. We interpret this climate dependency to reflect vegetation‐driven variations in apparent soil cohesion, which increases with root network density, and attenuation of rain splash and overland flow erosion, which increases with vegetative ground cover. To assess how vegetative root mass and ground cover vary with precipitation and aspect, we quantified the spatial distribution of NDVI (normalized difference vegetation index) from ASTER satellite images and observed spatial variations that correlate with our calibrated values of K . Analysis of previously studied cinder cones in the USA can be used to extend our framework to arid domains. This endeavor suggests a humped relationship between K and precipitation with maximum diffusivity at mean annual precipitation of 400–600 mm. Copyright © 2017 John Wiley & Sons, Ltd.     

Quantifying lateral root reinforcement in steep slopes – from a bundle of roots to tree stands

A review of present modelling approaches for root reinforcement in vegetated steep hillslopes reveals critical gaps in consideration of plant–soil interactions at various scales of interest for shallow landslide prediction. A new framework is proposed for systematic quantification of root reinforcement at scales ranging from single root to tree root system, to a stand of trees. In addition to standard basal reinforcement considered in most approaches, the critical role of roots in stabilizing slopes through lateral reinforcement is highlighted. Primary geometrical and mechanical properties of root systems and their function in stabilizing the soil mass are reviewed. Stress–strain relationships are considered for a bundle of roots using the formalism of the fiber bundle model (FBM) that offers a natural means for upscaling mechanical behavior of root systems. An extension of the FBM is proposed, considering key root and soil parameters such as root diameter distribution, tortuosity, soil type, soil moisture and friction between soil and root surface. The spatial distribution of root mechanical reinforcement around a single tree is computed from root diameter and density distributions based on easy to measure properties. The distribution of root reinforcement for a stand of trees was obtained from spatial and mechanical superposition of individual tree values with regard to their positions on a hillslope. Potential applications of the proposed approach are illustrated in a numerical experiment of spatial strength distribution in a hypothetical slope with 1000 trees randomly distributed. The analyses result in spatial distribution of weak and strong zones within the soil where landslide triggering is expected in large and continuous zones with low reinforcement values. Mapping such zones would enhance the quality of landslide susceptibility maps and optimization of silvicultural measures in protection forests. Copyright © 2010 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.     

Rainfall thresholds for shallow landsliding derived from pressure‐head monitoring: cases with permeable and impermeable bedrocks in Boso Peninsula,Japan

Rainfall thresholds for shallow landslide initiation were determined for hillslopes with two types of bedrock, permeable sandstone and impermeable mudstone, in the Boso Peninsula, Japan. The pressure‐head response to rainfall was monitored above a slip scarp due to earlier landslides. Multiple regression analysis estimated the rainfall thresholds for landsliding from the relation between the magnitude of the rainfall event and slope instability caused by the increased pressure heads. The thresholds were expressed as critical combinations of rainfall intensity and duration, incorporating the geotechnical properties of the hillslope materials and also the slope hydrological processes. The permeable sandstone hillslope has a greater critical rainfall and hence a longer recurrence interval than the impermeable mudstone hillslope. This implies a lower potential for landsliding in sandstone hillslopes, corresponding to lower landslide activity. Copyright © 2007 John Wiley & Sons, Ltd.     

MATHEMATICAL MODEL OF OVERLAND FLOW AND MECHANISM OF SOIL CONSERVATION FOR FORESTED STEEP HILLSLOPE

ＭＡＴＨＥＭＡＴＩＣＡＬＭＯＤＥＬＯＦＯＶＥＲＬＡＮＤＦＬＯＷＡＮＤＭＥＣＨＡＮＩＳＭＯＦＳＯＩＬＣＯＮＳＥＲＶＡＴＩＯＮＦＯＲＦＯＲＥＳＴＥＤＳＴＥＥＰＨＩＬＬＳＬＯＰＥ（１）ＡｎａｌｙｔｉｃａｌＳｏｌｕｔｉｏｎｔｏｔｈｅＯｖｅｒｌａｎｄＦｌｏｗｏ...     

The sensitivity of hillslope bedrock erosion to precipitation

Decoupling the impacts of climate and tectonics on hillslope erosion rates is a challenging problem. Hillslope erosion rates are well known to respond to changes in hillslope boundary conditions (e.g. channel incision rates) through their dependence on soil thickness, and precipitation is an important control on soil formation. Surprisingly though, compilations of hillslope denudation rates suggest little precipitation sensitivity. To isolate the effects of precipitation and boundary condition, we measured rates of soil production from bedrock and described soils on hillslopes along a semi‐arid to hyperarid precipitation gradient in northern Chile. In each climate zone, hillslopes with contrasting boundary conditions (actively incising channels versus non‐eroding landforms) were studied. Channel incision rates, which ultimately drive hillslope erosion, varied with precipitation rather than tectonic setting throughout the study area. These precipitation‐dependent incision rates are mirrored on the hillslopes, where erosion shifts from relatively fast and biologically‐driven to extremely slow and salt‐driven as precipitation decreases. Contrary to studies in humid regions, bedrock erosion rates increase with precipitation following a power law, from ～1 m Ma^?1 in the hyperarid region to ～40 m Ma^?1 in the semi‐arid region. The effect of boundary condition on soil thickness was observed in all climate zones (thicker soils on hillslopes with stable boundaries compared to hillslopes bounded by active channels), but the difference in bedrock erosion rates between the hillslopes within a climate region (slower erosion rates on hillslopes with stable boundaries) decreased as precipitation decreased. The biotic‐abiotic threshold also marks the precipitation rate below which bedrock erosion rates are no longer a function of soil thickness. Our work shows that hillslope processes become sensitive to precipitation as life disappears and the ability of the landscape to respond to tectonics decreases. Copyright © 2010 John Wiley & Sons, Ltd.     

How deep can forest vegetation cover extend their hydrological reinforcing contribution?

An experimental campaign was set up to quantify the contribution of evapotranspiration fluxes on hillslope hydrology and stability for different forest vegetation cover types. Three adjacent hillslopes, respectively, covered by hardwood, softwood, and grass were instrumented with nine access tubes each to monitor soil water dynamics at the three depths of 30, 60, and 100 cm, using a PR2/6 profile probe (Delta‐T Devices Ltd) for about 6 months including wet periods. Soil was drier under softwood and wetter under grass at all the three depths during most of the monitoring period. Matric suction derived via the soil moisture measurements was more responsive to changes in the atmospheric conditions and also recovered faster at the 30 cm depth. Results showed no significant differences between mean matric suction under hardwood (101.6 kPa) with that under either softwood or grass cover. However, a significant difference was found between mean matric suction under softwood (137.5 kPa) and grass (84.3 kPa). Results revealed that, during the wettest period, the hydrological effects from all three vegetation covers were substantial at the 30 cm depth, whereas the contribution from grass cover at 60 cm (2.0 kPa) and 100 cm (1.1 kPa) depths and from hardwood trees at 100 cm depth (1.2 kPa) was negligible. It is surmised that potential instability would have occurred at these larger depths along hillslopes where shallow hillslope failures are most likely to occur in the region. The hydrological effects from softwood trees, 8.1 and 3.9 kPa, were significant as the corresponding factor of safety values showed stable conditions at both depths of 60 and 100 cm, respectively. Therefore, the considerable hydrological reinforcing effects from softwood trees to the 100 cm depth suggest that a hillslope stability analysis would show that hillslopes with softwood trees will be stable even during the wet season.     

A numerical investigation of bedrock groundwater recharge and exfiltration on soil mantled hillslopes

Here we use Richards Equation models of variably saturated soil and bedrock groundwater flow to investigate first-order patterns of the coupling between soil and bedrock flow systems. We utilize a Monte Carlo sensitivity analysis to identify important hillslope parameters controlling bedrock recharge and then model the transient response of bedrock and soil flow to seasonal precipitation. Our results suggest that hillslopes can be divided into three conceptual zones of groundwater interaction, (a) the zone of lateral unsaturated soil moisture accumulation (upper portion of hillslope), (b) the zone of soil saturation and bedrock recharge (middle of hillslope) and (c) the zone of saturated-soil lateral flow and bedrock groundwater exfiltration (bottom of hillslope). Zones of groundwater interaction expand upslope during periods of precipitation and drain downslope during dry periods. The amount of water partitioned to the bedrock groundwater system a can be predicted by the ratio of bedrock to soil saturated hydraulic conductivity across a variety of hillslope configurations. Our modelled processes are qualitatively consistent with observations of shallow subsurface saturation and groundwater fluctuation on hillslopes studied in our two experimental watersheds and support a conceptual model of tightly coupled shallow and deep subsurface circulation where groundwater recharge and discharge continuously stores and releases water from longer residence time storage.