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.     

Downslope soil detachment–transport on steep slopes via rain splash

This study developed a one‐dimensional model of downslope rain splash transport based on field experiments and previous studies. The developed model considers soil detachment processes, ground cover, probability densities, and the effect of overland run‐off in preventing detachment. Field monitoring was conducted to observe precipitation run‐off, ground cover, and sediment production on steep hillslopes. Field‐observed data were used to develop the splash detachment rate equation, probability densities for splash transport, and the maximum splash transport distance. Observed and estimated splash transport showed overall agreement, with some differences for small storm events or events with relatively low intensity, probably caused by variation of overland run‐off depth and connectivity as well as differences in soil surface cohesion at various degrees of wetness. Our model can provide insights on the interactions among rainfall intensity, soil surface condition, soil wetness, and splash transport on forested hillslopes. Copyright © 2011 John Wiley & Sons, Ltd.     

Gopher bioturbation: field evidence for non‐linear hillslope diffusion

It has generally been assumed that diffusive sediment transport on soil‐mantled hillslopes is linearly dependent on hillslope gradient. Fieldwork was done near Santa Barbara, California, to develop a sediment transport equation for bioturbation by the pocket gopher (Thomomys bottae) and to determine whether it supports linear diffusion. The route taken by the sediment is divided into two parts, a subsurface path followed by a surface path. The first is the transport of soil through the burrow to the burrow opening. The second is the discharge of sediment from the burrow opening onto the hillslope surface. The total volumetric sediment flux, as a function of hillslope gradient, is found to be: q_s (cm³ cm⁻¹ a⁻¹) = 176(dz/dx)³ − 189(dz/dx)² + 68(dz/dx) + 34(dz/dx)^0·4. This result does not support the use of linear diffusion for hillslopes where gopher bioturbation is the dominant mode of sediment transport. A one‐dimensional hillslope evolution program was used to evolve hillslope profiles according to non‐linear and linear diffusion and to compare them to a typical hillslope. The non‐linear case more closely resembles the actual profile with a convex cap at the divide leading into a straight midslope section. Copyright © 2000 John Wiley & Sons, Ltd.     

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.     

Mining soil databases for landscape‐scale patterns in the abundance and size distribution of hillslope rock fragments

Landscape‐scale variation in rock fragments on soil‐mantled hillslopes is poorly understood, despite the potential importance of rock fragments in soil weathering and coarse sediment supply to river networks. We explored the utility of soil survey databases for data mining, with the goals of identifying landscape‐scale patterns in the abundance and size distribution of rock fragments (diameter D > 2 mm) and potential controls on grain size production. We focus on data from three regions: the Hawaiian Islands, and the Sierra Nevada and Cascade Mountains, where elevation transects span a range of environmental conditions. We selected pedons from pits dug on hillslopes with active soil production and transport. For the 27 pedons selected, we constructed depth‐averaged grain size distributions and calculated the mass fraction of rock fragments (F_RF) and the median rock fragment grain size (D_50RF). We also categorized as bimodal, size distributions with a clear ‘breakpoint’ between fine and coarse modes. Several strong patterns emerge from the data. We find rock fragments in 85% of the pedons, primarily in distinct coarse modes within bimodal size distributions. Values of F_RF and D_50RF are strongly correlated, although the best‐fit power law scaling between F_RF and D_50RF differs between the warmer Hawaiian, and colder Sierra Nevada and Cascade Mountain sites. We also find a regional contrast in the variation in F_RF with elevation; F_RF declines with elevation in Hawaii, but increases in the mainland sites. Although this contrast could be an artifact of variable lithology, precipitation may influence many patterns in the data. Lower mean‐annual precipitation correlates with higher F_RF, dominantly bimodal distributions and surface enrichment in the vertical distribution of rock fragments. These observations may be useful in refining models of coarse sediment supply to rivers, and suggest opportunities for future work to test mechanistic hypotheses for rock fragment production on soil‐mantled hillslopes. Copyright © 2011 John Wiley & Sons, Ltd.     

Hillslope response to tectonic forcing in threshold landscapes

Hillslopes are thought to poorly record tectonic signals in threshold landscapes. Numerous previous studies of steep landscapes suggest that large changes in long‐term erosion rate lead to little change in mean hillslope angle, measured at coarse resolution. New LiDAR‐derived topography data enables a finer examination of threshold hillslopes. Here we quantify hillslope response to tectonic forcing in a threshold landscape. To do so, we use an extensive cosmogenic beryllium‐10 (¹⁰Be)‐based dataset of catchment‐averaged erosion rates combined with a 500 km² LiDAR‐derived 1 m digital elevation model to exploit a gradient of tectonic forcing and topographic relief in the San Gabriel Mountains, California. We also calibrate a new method of quantifying rock exposure from LiDAR‐derived slope measurements using high‐resolution panoramic photographs. Two distinct trends in hillslope behavior emerge: below catchment‐mean slopes of 30°, modal slopes increase with mean slopes, slope distribution skewness decreases with increasing mean slope, and bedrock exposure is limited; above mean slopes of 30°, our rock exposure index increases strongly with mean slope, and the prevalence of angle‐of‐repose debris wedges keeps modal slopes near 37°, resulting in a positive relationship between slope distribution skewness and mean slope. We find that both mean slopes and rock exposure increase with erosion rate up to 1 mm/a, in contrast to previous work based on coarser topographic data. We also find that as erosion rates increase, the extent of the fluvial network decreases, while colluvial channels extend downstream, keeping the total drainage density similar across the range. Our results reveal important textural details lost in 10 or 30 m resolution digital elevation models of steep landscapes, and highlight the need for process‐based studies of threshold hillslopes and colluvial channels. Copyright © 2012 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.     

Modelling investigation of water partitioning at a semiarid ponderosa pine hillslope

The effects of vegetation root distribution on near‐surface water partitioning can be two‐fold. On the one hand, the roots facilitate deep percolation by root‐induced macropore flow; on the other hand, they reduce the potential for deep percolation by root‐water‐uptake processes. Whether the roots impede or facilitate deep percolation depends on various conditions, including climate, soil, and vegetation characteristics. This paper examines the effects of root distribution on deep percolation into the underlying permeable bedrock for a given soil profile and climate condition using HYDRUS modelling. The simulations were based on previously field experiments on a semiarid ponderosa pine (Pinus ponderosa) hillslope. An equivalent single continuum model for simulating root macropore flow on hillslopes is presented, with root macropore hydraulic parameterization estimated based on observed root distribution. The sensitivity analysis results indicate that the root macropore effect dominates saturated soil water flow in low conductivity soils (K_matrix below 10^?7 m/s), while it is insignificant in soils with a K_matrix larger than 10^?5 m/s, consistent with observations in this and other studies. At the ponderosa pine site, the model with simple root‐macropore parameterization reasonably well reproduces soil moisture distribution and some major runoff events. The results indicate that the clay‐rich soil layer without root‐induced macropores acts as an impeding layer for potential groundwater recharge. This impeding layer results in a bedrock percolation of less than 1% of the annual precipitation. Without this impeding layer, percolation into the underlying permeable bedrock could be as much as 20% of the annual precipitation. This suggests that at a surface with low‐permeability soil overlying permeable bedrock, the root penetration depth in the soil is critical condition for whether or not significant percolation occurs. Copyright © 2010 John Wiley & Sons, Ltd.     

Strong variation in weathering of layered rock maintains hillslope‐scale strength under high precipitation

The evolution of volcanic landscapes and their landslide potential are both dependent upon the weathering of layered volcanic rock sequences. We characterize critical zone structure using shallow seismic V_p and V_s profiles and vertical exposures of rock across a basaltic climosequence on Kohala peninsula, Hawai’i, and exploit the dramatic gradient in mean annual precipitation (MAP) across the peninsula as a proxy for weathering intensity. Seismic velocity increases rapidly with depth and the velocity–depth gradient is uniform across three sites with 500–600 mm/yr MAP, where the transition to unaltered bedrock occurs at a depth of 4 to 10 m. In contrast, velocity increases with depth less rapidly at wetter sites, but this gradient remains constant across increasing MAP from 1000 to 3000 mm/yr and the transition to unaltered bedrock is near the maximum depth of investigation (15–25 m). In detail, the profiles of seismic velocity and of weathering at wet sites are nowhere monotonic functions of depth. The uniform average velocity gradient and the greater depths of low velocities may be explained by the averaging of velocities over intercalated highly weathered sites with less weathered layers at sites where MAP > 1000 mm/yr. Hence, the main effect of climate is not the progressive deepening of a near‐surface altered layer, but rather the rapid weathering of high permeability zones within rock subjected to precipitation greater than ~1000 mm/yr. Although weathering suggests mechanical weakening, the nearly horizontal orientation of alternating weathered and unweathered horizons with respect to topography also plays a role in the slope stability of these heterogeneous rock masses. We speculate that where steep, rapidly evolving hillslopes exist, the sub‐horizontal orientation of weak/strong horizons allows such sites to remain nearly as strong as their less weathered counterparts at drier sites, as is exemplified by the 50°–60° slopes maintained in the amphitheater canyons on the northwest flank of the island. Copyright © 2017 John Wiley & Sons, Ltd.     

A one‐dimensional model for simulating armouring and erosion on hillslopes: 2. Long term erosion and armouring predictions for two contrasting mine spoils

This paper investigates the dynamics of soil armouring as a result of fluvial erosion for a non‐cohesive sandy gravel spoil from the Ranger Mine, Australia, and a cohesive silt loam spoil from the Northparkes Mine, Australia, using a model for hillslope soil armouring. These long term predictions concentrate on the temporal and spatial changes of the spoil grading and erosion over 100–200 years for the flat cap regions (1–2%) and steep batter edges (10–30%) typically encountered on waste rock dumps. The existence of a significant rock fragment fraction in the Ranger spoil means that it armours readily, while Northparkes does not. For Ranger the waste rock showed reductions in (1) cumulative erosion of up to 81% from that obtained by extrapolating the initial erosion rate out 100 years and (2) the erosion/year by more than 10‐fold. For Northparkes reductions were less marked, with the maximum reduction in erosion/year being 37% after 200 years. For Ranger the reductions were greatest and fastest for intermediate gradient hillslopes. For the steepest hillslopes the armouring decreased because the flow shear stresses were large enough to mobilize all material in the armour layer. Model uncertainty was assessed with probabilistic confidence limits demonstrating that these erodibility reductions were statistically significant. A commonly used hillslope erosion model (sediment flux = β₁ discharge m₁ slope n₁) was fitted to these predictions. The erodibility, β₁, and m₁ decreased with time, which was consistent with our physical intuition about armouring. At Ranger the parameter m₁ asymptoted to 1·5–1·6 while at Northparkes it asymptoted to 1·2–1·3. At Ranger transient spatial trends in armouring led to a short term (50–200 years in the future) reduction in n₁, to below zero under certain circumstances, recovering to an asymptote of about 0·5–1. At Northparkes n₁ asymptoted to about 0·6, with no negative transients predicted. The m₁ and n₁ parameters predicted for Ranger were shown to be consistent with field data from a 10‐year‐old armoured hillslope and consistent with published relationships between erodibility and rock content for natural hillslopes. Copyright © 2007 John Wiley & Sons, Ltd.     

The roles of channels and hillslopes in rainfall/run‐off lag times during intense storms in a steep catchment

The response time (lag time) between rainfall input and run‐off output in headwater catchments is a key parameter for flood prediction. Lag times are expected to be controlled by run‐off processes, both on hillslopes and in channels. To demonstrate these effects on peak lag times within a 4.5‐km² catchment, we measured stream water levels at up to 16 channel locations at 1‐min intervals and compared the lag times with topographic indices describing the length and gradient of the hillslope and channel flow path. We captured storm events with a total precipitation of 38–198 mm and maximum hourly precipitation intensity of 9–90 mm/hr. There were positive relationships between lag time and flow path length as well as the ratio of the flow path length and the square root of the gradient of channels for the most intense storms, demonstrating that channel flow paths generally defined the variation in lag times. Topographic analysis showed that hillslope flow path lengths were similar among locations, whereas channel flow path length increased almost one order of magnitude with a 100‐fold increase in catchment area. Thus, the relative importance of hillslope flow path decreased with increasing catchment area. Our results indicate that the variation in lag times is small when hillslopes are sufficiently wet; thus, catchment‐scale variation in lag times can be explained almost entirely by channel processes. Detailed topographic channel information can improve prediction of flood peak timing, whereas hillslopes can be treated as homogeneous during large flood events.     

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.     

Lateral flow thresholds for aspen forested hillslopes on the Western Boreal Plain,Alberta, Canada

To predict the long‐term sustainability of water resources on the Boreal Plain region of northern Alberta, it is critical to understand when hillslopes generate runoff and connect with surface waters. The sub‐humid climate (P ≤ ET) and deep glacial sediments of this region result in large available soil storage capacity relative to moisture surpluses or deficits, leading to threshold‐dependent rainfall‐runoff relationships. Rainfall simulation experiments were conducted using large magnitude and high intensity applications to examine the thresholds in precipitation and soil moisture that are necessary to generate lateral flow from hillslope runoff plots representative of Luvisolic soils and an aspen canopy. Two adjacent plots (areas of 2·95 and 3·4 m²) of contrasting antecedent moisture conditions were examined; one had tree root uptake excluded for two months to increase soil moisture content, while the second plot allowed tree uptake over the growing season resulting in drier soils. Vertical flow as drainage and soil moisture storage dominated the water balances of both plots. Greater lateral flow occurred from the plot with higher antecedent moisture content. Results indicate that a minimum of 15–20 mm of rainfall is required to generate lateral flow, and only after the soils have been wetted to a depth of 0·75 m (C‐horizon). The depth and intensity of rainfall events that generated runoff > 1 mm have return periods of 25 years or greater and, when combined with the need for wet antecendent conditions, indicate that lateral flow generation on these hillslopes will occur infrequently. Copyright © 2008 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.     

Evaluating young fluvial terrace riser degradation using a nonlinear transport model: Application to the Kongur Normal Fault in the Pamir,northwest China

Selecting an appropriate hillslope transport model and calibrating model parameters are essential for morphological dating of fault and fluvial scarps. In this paper, we refine the method of profile-based morphologic dating by updating the representation of nonlinearity in sediment flux dependence on the hillslope gradient. We apply this revised method to fluvial scarps bounding fluvial terraces offset along the Kongur Normal Fault in the semi-arid high-altitude Pamir mountains, northwestern China. One of these terraces, the T3 surface, is dated to 7.0^+1.9/_−1.6 kyr using ¹⁰Be cosmogenic depth profile analysis. Well-preserved, dated terraces make this an ideal site to test the utility of morphological modelling in constraining ages of the young terrace risers. To do this, 35 topographic swath profiles across the terrace risers are extracted from a 0.2 m-resolution digital elevation model produced using structure from motion from photos collected by an unmanned aerial vehicle. The best estimates of morphological age are 13.9 ± 1.3 m² for the riser T3/T4 and 11.9 ± 1.3 m² for T2/T3 using a linear diffusion approach. These two morphological ages overlap within uncertainty and fail to distinguish between two young terrace risers. Alternatively, we employed a nonlinear diffusion model, calibrated with transport constant k = 1 m² kyr⁻¹, nonlinearity n = 2, and critical gradient S_c = tan(33°). This nonlinear model produces ages of 7.3 ± 0.5 kyr for T3/T4 and 4.0 ± 0.2 kyr for T2/T3; these ages are consistent with terrace surface ages deduced by using vertical offset divided by independently determining average throw rate. This comparison shows the advantage of a nonlinear model in defining ages of young scarps. Furthermore, we explored the minor effect of heterogeneous degradation along steep sections of the scarp profiles. The nonlinear scarp modelling scheme we develop in this paper is suitable for studying scarp degradation in other regions. © 2020 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.     

Stone movement on hillslopes in the Mojave Desert,California: A 16-year record

Previous research has shown that either hydraulic action or creep may be the dominant process transporting coarse debris down hillslopes in the American Southwest. This study analyses the movement over 16 years of painted stones on two hillslopes in the central Mojave Desert to ascertain which of these two processes dominate in this region. The distance moved (M) is found to be directly related to length of overland flow (X) and hillslope gradient (S), and inversely related to particle size (D). The fact that M is more highly correlated with X than with S suggests that hydraulic action rather than creep is the dominant process. It is concluded that this is probably the case over most of the Mojave and Sonoran Deserts on slopes up to at least 24°, and that it is only at higher elevations where winters are more severe that creep may become dominant.     

Mechanisms and pathways of lateral flow on aspen‐forested,Luvisolic soils,Western Boreal Plains,Alberta, Canada

Rainfall simulation experiments by Redding and Devito ( 2008 , Hydrological Processes 23: 4287–4300) on two adjacent plots of contrasting antecedent soil moisture storage on an aspen‐forested hillslope on the Boreal Plain showed that lateral flow generation occurred only once large soil storage capacity was saturated combined with a minimum event precipitation of 15–20 mm. This paper extends the results of Redding and Devito ( 2008 , Hydrological Processes 23: 4287–4300) with detailed analysis of pore pressure, soil moisture and tracer data from the rainfall simulation experiments, which is used to identify lateral flow generation mechanisms and flow pathways. Lateral flow was not generated until soils were wet into the fine textured C horizon. Lateral flow occurred dominantly through the clay‐rich Bt horizon by way of root channels. Lateral flow during the largest event was dominated by event water, and precipitation intensity was critical in lateral flow generation. Lateral flow was initiated as preferential flow near the soil surface into root channels, followed by development of a perched water table at depth, which also interacted with preferential flow pathways to move water laterally by the transmissivity feedback mechanism. The results indicate that lateral flow generated by rainfall on these hillslopes is uncommon because of the generally high available soil moisture storage capacity and the low probability of rainfall events of sufficient magnitude and intensity. Copyright © 2010 John Wiley & Sons, Ltd.     

Measurements of rain splash on bench terraces in a humid tropical steepland environment

Soil loss continues to threaten Java's predominantly bench‐terraced volcanic uplands. Sediment transport processes on back‐sloping terraces with well‐aggregated clay‐rich oxisols in West Java were studied using two different techniques. Splash on bare, cropped, or mulched sub‐horizontal (2–3°) terrace beds was studied using splash cups of different sizes, whereas transport of sediment on the predominantly bare and steep (30–40/deg ) terrace risers was measured using a novel device combining a Gerlach‐type trough with a splash box to enable the separate measurement of transport by wash and splash processes. Measurements were made during two consecutive rainy seasons. The results were interpreted using a recently developed splash distribution theory and related to effective rainfall erosive energy. Splash transportability (i.e. transport per unit contour length and unit erosive energy) on the terrace risers was more than an order of magnitude greater than on bare terrace beds (0·39–0·57 versus 0·013–0·016 g m J^?1). This was caused primarily by a greater average splash distance on the short, steep risers (>11 cm versus c. 1 cm on the beds). Splashed amounts were reduced by the gradual formation of a protective ‘pavement’ of coarser aggregates, in particular on the terrace beds. Soil aggregate size exhibited an inverse relationship with detachability (i.e. detachment per unit area and unit erosive energy) and average splash length, and therefore also with transportability, as did the degree of canopy and mulch cover. On the terrace risers, splash‐creep and gravitational processes transported an additional 6–50% of measured rain splash, whereas transport by wash played a marginal role. Copyright © 2002 John Wiley & Sons, Ltd.     

A methodology to study rain splash and wash processes under natural rainfall

The complex interactions between rainfall‐driven erosion processes and rainfall characteristics, slope gradient, soil treatment and soil surface processes are not very well understood. A combination of experiments under natural rainfall and a consistent physical theory for their interpretation is needed to shed more light on the underlying processes. The present study demonstrates such a methodology. An experimental device employed earlier in laboratory studies was used to measure downslope rain splash and ‘splash‐creep’, lateral splash, upslope splash and rainfall‐driven runoff transport (wash) from a highly aggregated clay‐rich oxisol exposed to natural rainfall in West Java, Indonesia. Two series of measurements were made: the first with the soil surface at angles of 0°, 5°, 15° and 40°; and the second all at an angle of 5° but with different tillage and mulching treatments. A number of rainfall erosivity indices were calculated from rainfall intensity measurements and compared with measured transport components. Overall storm kinetic energy correlated reasonably well with sediment transport, but much better agreement was obtained when a threshold rainfall intensity (20 mm h^?1) was introduced. Rain splash transport measurements were interpreted using a recently developed theory relating detachment to sediment transport. Furthermore, a conceptually sound yet simple wash transport model is advanced that satisfactorily predicted observed washed sediment concentrations. The lack of replication precluded rigorous assessment of the effect of slope and soil treatment on erosion processes, but some general conclusions could still be drawn. The results stress the importance of experiments under conditions of natural rainfall. Copyright © 2002 John Wiley & Sons, Ltd.