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.     

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.     

Spatial and temporal variability of pore‐water pressures in residual soil slopes in a tropical climate

It is critical to understand and quantify the temporal and spatial variability in hillslope hydrological data in order to advance hillslope hydrological studies, evaluate distributed parameter hydrological models, analyse variability in hydrological response of slopes and design efficient field data sampling networks. The spatial and temporal variability of field‐measured pore‐water pressures in three residual soil slopes in Singapore was investigated using geostatistical methods. Parameters of the semivariograms, namely the range, sill and nugget effect, revealed interesting insights into the spatial structure of the temporal situation of pore‐water pressures in the slopes. While informative, mean estimates have been shown to be inadequate for modelling purposes, indicator semivariograms together with mean prediction by kriging provide a better form of model input. Results also indicate that significant temporal and spatial variability in pore‐water pressures exists in the slope profile and thereby induces variability in hydrological response of the slope. Spatial and temporal variability in pore‐water pressure decreases with increasing soil depth. The variability decreases during wet conditions as the slope approaches near saturation and the variability increases with high matric suction development following rainfall periods. Variability in pore‐water pressures is greatest at shallow depths and near the slope crest and is strongly influenced by the combined action of microclimate, vegetation and soil properties. Copyright © 2002 John Wiley & Sons, Ltd.     

Exploring geochemical controls on weathering and erosion of convex hillslopes: beyond the empirical regolith production function

Landscape curvature evolves in response to physical, chemical, and biological influences that cannot yet be quantified in models. Nonetheless, the simplest models predict the existence of equilibrium hillslope profiles. Here, we develop a model describing steady‐state regolith production caused by mineral dissolution on hillslopes which have attained an equilibrium parabolic profile. When the hillslope lowers at a constant rate, the rate of chemical weathering is highest at the ridgetop where curvature is highest and the ridge develops the thickest regolith. This result derives from inclusion of all the terms in the mathematical definition of curvature. Including these terms shows that the curvature of a parabolic hillslope profile varies with distance from the ridge. The hillslope model (meter‐scale) is similar to models of weathering rind formation (centimeter‐scale) where curvature‐driven solute transport causes development of the thickest rinds at highly curved clast corners. At the clast scale, models fit observations. Here, we similarly explore model predictions of the effect of curvature at the hillslope scale. The hillslope model shows that when erosion rates are small and vertical porefluid infiltration is moderate, the hill weathers at both ridge and valley in the erosive transport‐limited regime. For this regime, the reacting mineral is weathered away before it reaches the land surface: in other words, the model predicts completely developed element‐depth profiles at both ridge and valley. In contrast, when the erosion rate increases or porefluid velocity decreases, denudation occurs in the weathering‐limited regime. In this regime, the reacting mineral does not weather away before it reaches the land surface and simulations predict incompletely developed profiles at both ridge and valley. These predictions are broadly consistent with observations of completely developed element‐depth profiles along hillslopes denuding under erosive transport‐limitation but incompletely developed profiles along hillslopes denuding under weathering limitation in some field settings. Copyright © 2013 John Wiley & Sons, Ltd.     

PLAN FORM AND PROFILE FORM OF HILLSLOPES

For 27 hillslope profiles located within four first-order drainage basins on the dip slope of the South Downs in East Sussex, contour curvature (P), measured in degrees per 100 m, is highly variable within any one hillslope. Hillslopes may be described as wholly concave in plan where P is less than ?40, wholly convex in plan where P is greater than +40 and convexo-concave in plan where P lies between ?40 and +40. One measure of P taken at the profile sampling line is found to give a better estimate of P than measuring P at the steepest point on the hillslope. P is found to be significantly negatively correlated with slope in 38 per cent of cases, correlated with distance from the divide in 55 per cent of cases (positively where P <0 and negatively where P <0) and positively correlated with profile curvature in 70 per cent of cases. Implications of these correlations for hillslope process-response models are discussed.     

Modelling of sediment dynamics in a laboratory‐scale experimental catchment

The variability of hillslope form and function is examined experimentally using a simple model catchment in which most landscape development parameters are either known or controlled. It is demonstrated that there is considerable variability in sediment output from similar catchments, subjected to the same hydrological processes, and for which the initial hillslope profiles are the same. The results demonstrate that, in the case of catchments with a linear initial hillslope profile, the sediment output is initially high but reduces through time, whereas for a concave initial profile the sediment output was smaller and relatively constant. Concave hillslope profiles also displayed reduced sediment output when compared with linear slopes with the same overall slope. Using this experimental model catchment data, the SIBERIA landscape evolution model was tested for its ability to predict temporal sediment transport. When calibrated for the rainfall and erodible material, SIBERIA is able to simulate mean temporal sediment output for the experimental catchment over a range of hillslope profiles and rainfall intensities. SIBERIA is also able to match the hillslope profile of the experimental catchments. The results of the study provide confidence in the ability of SIBERIA to predict temporal sediment output. The experimental and modelling data also demonstrate that, even with all geomorphic and hydrological variables being known and/or controlled, there is still a need for long‐term stream gauging to obtain reliable assessments of field catchment hydrology and sediment transport. Copyright © 2005 John Wiley & Sons, Ltd.     

Paraglacial modification of slope form

The morphological consequences of paraglacial modification of valley-side drift slopes are investigated at six sites in Norway. Here, paraglacial slope adjustment operates primarily through the development of gully systems, whereby glacigenic sediment is stripped from the upper drift slope and redeposited in debris cones downslope. This results in an overall lowering of average gradient by up to 4·5° along gully axes. In general, slope profile adjustment appears to be characterized by a convergence of slope profiles towards an ‘equilibrium form’ with an upper rectilinear slope gradient at 29°± 4° and a range of concavities of approximately 0·0 to 0·4. After initial rapid incision, further gully deepening is limited, but gullies become progressively wider as sidewall gradients decline to c. 25°, after which parallel retreat appears to predominate. The final form of mature paraglacial gully systems consists of an upper bedrock-floored source area, a mid-slope area of broad gullies whose sidewalls rest at stable, moderate gradients, and a lower slope zone where gullies discharge onto the surfaces of debris cones and fans. Some gullies appear to have attained this final form and have stabilized following exhaustion of readily entrainable sediment within decades of gully initiation. At most sites, paraglacial activity has transformed steep drift-mantled valley sides into gullied slopes where an average of c. 2–3 m of surface lowering has taken place. At the most active sites, these average amounts imply minimum erosion rates averaging c. 90 mm a⁻¹ since gully initiation, which highlights the extreme rapidity of paraglacial erosion of deglaciated drift-mantled slopes. Copyright © 1999 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.     

The topographic signature of a major typhoon

In August 2009, the typhoon Morakot, characterized by a cumulative rainfall up to 2884 mm in about three days, triggered thousands of landslides in Taiwan. The availability of LiDAR surveys before (2005) and after (2010) this event offers a unique opportunity to investigate the topographic signatures of a major typhoon. The analysis considers the comparison of slope–area relationships derived by LiDAR digital terrain models (DTMs). This approach has been successfully used to distinguish hillslope from channelized processes, as a basis to develop landscape evolution models and theories, and understand the linkages between landscape morphology and tectonics, climate, and geology. We considered six catchments affected by a different degree of erosion: three affected by shallow and deep‐seated landslides, and three not affected by erosion. For each of these catchments, 2 m DTMs were derived from LiDAR data. The scaling regimes of local slope versus drainage area suggested that for the catchments affected by landslides: (i) the hillslope‐to‐valley transitions morphology, for a given value of drainage area, is shifted towards higher value of slopes, thus indicating a likely migration of the channelized processes and erosion toward the catchment boundary (the catchment head becomes steeper because of erosion); (ii) the topographic gradient along valley profiles tends to decrease progressively (the valley profile becomes gentler because of sediment deposition after the typhoon). The catchments without any landslides present a statistically indistinguishable slope–area scaling regime. These results are interesting since for the first time, using multi‐temporal high‐resolution topography derived by LiDAR, we demonstrated that a single climate event is able to cause significant major geomorphic changes on the landscape, detectable using slope–area scaling analysis. This provides new insights about landscape evolution under major climate forcing. Copyright © 2015 John Wiley & Sons, Ltd.     

Landscape position,surface hydraulic gradients and erosion processes

Different hydraulic gradients, especially due to seepage or drainage, at different locations on a hillslope profile may have a profound effect on the dominant erosion processes. A laboratory study was designed to simulate hillslope processes and quantify effects of surface hydraulic gradients on erosion for a Glynwood clay loam soil (fine, illitic, mesic Aquic Hapludalf). A 5 m long, 1·2 m wide soil pan was used at 5 and 10 per cent slopes with an external watering tube to vary the soil bed's hydrological conditions. Different combinations of slope steepness with seepage or drainage gradients were used to simulate the hydrologic conditions on a 5 m segment of a hillslope profile. Runoff samples were taken during rainfall-only and rainfall with added inflow. Results showed that, under drainage conditions, interrill processes dominated and rilling was limited. The surface contained scattered crescent-shaped pits after the run. Under seepage conditions, rilling processes dominated and the inflow introduced at the top of the soil pan further accelerated the headward erosion of the rills. Erosion rates increased by as much as 60 times under seepage conditions representative of the lower backslope when compared to drainage conditions that generally occur at the upper backslope. This indicated that rills and gullies on backslopes and footslopes may be catalysed or enhanced by seepage conditions rather than form from flow hydraulic shear stress alone. An understanding of spatial and temporal changes that affect both hillslope hydrology and erosional processes is needed to develop accurate process-based erosion prediction models. This knowledge may lead to different management practices on landscape positions where seepage occurs. © 1998 John Wiley & Sons, Ltd.     

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.     

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.     

Quantifying river incision into low‐relief surfaces using local and catchment‐wide 10Be denudation rates

Relief generation in non‐glaciated regions is largely controlled by river incision into bedrock but datable fluvial terraces that allow quantifying incision rates are not always present. Here we suggest a new method to determine river incision rates in regions where low‐relief surfaces are dissected by streams. The approach consists of three steps and requires the ¹⁰Be concentrations of a stream sediment sample and a regolith sample from the low‐relief surface. In the first step, the spatial distribution of ¹⁰Be surface concentrations in the given catchment is modelled by assuming that denudation rates are controlled by the local hillslope angles. The slope–denudation rate relation for this catchment is then quantified by adjusting the relation between slope angle and denudation rate until the average ¹⁰Be concentration in the model is equal to the one measured in the stream sediment sample. In the second step, curved swath profiles are used to measure hillslope angles adjacent to the main river channel. Third, the mean slope angle derived from these swath profiles and the slope–denudation relation are used to quantify the river incision rate (assuming that the incision rate equals the denudation rate on adjacent hillslopes). We apply our approach to two study areas in southern Tibet and central Europe (Black Forest). In both regions, local ¹⁰Be denudation rates on flat parts of the incised low‐relief surface are lower than catchment‐wide denudation rates. As the latter integrate across the entire landscape, river incision rates must exceed these spatially averaged denudation rates. Our approach yields river incision rates between ~15 and ~30 m/Ma for the Tibetan study area and incision rates of ~70 to ~100 m/Ma in the Black Forest. Taking the lowering of the low‐relief surfaces into account suggests that relief in the two study areas increases at rates of 10–20 and 40–70 m/Ma, respectively. Copyright © 2018 John Wiley & Sons, Ltd.     

Spatial patterns of sediment delivery to valley floors: sensitivity to sediment transport capacity and hillslope hydrology relations

There is considerable interest in large‐scale spatial patterns of sediment transport in catchments, and this topic is often approached using terrain‐based modelling. In such models topography influences the discharge of overland flow and its sediment transport capacity. The sediment transport capacity of overland flow is commonly expressed as a power function of slope and discharge (i.e. q_s=k₁q^βS^γ). The relationship between discharge and contributing area can also be expressed as a power function. Several reviews reveal a limited range of values for the two exponents β and γ. In this paper we examine the sensitivity of catchment‐scale patterns of sediment delivery to valley floors to a range of sediment transport capacity and hillslope hydrology parameterizations, using two catchments on the southern tablelands of New South Wales. The results indicate that, over the limited range of β and γ identified within the literature, sediment deliveries to valley floors across the two catchments are similar for all but one of five sediment transport capacity relationships. The patterns are dominated by the trend in slope through each catchment. The sensitivity to hillslope hydrology of predicted sediment delivery patterns is strong in the catchment with systematic variation in unit hillslope area, and weak in the catchment for which there are no systematic trends in unit hillslope area. We believe there is less experimental evidence to restrict choice of hillslope hydrology parameters than there is for sediment transport capacity. Copyright © 2001 John Wiley & Sons, Ltd.     

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.     

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.     

Soil saturation patterns in steep,convergent hill-slopes under forest and pasture vegetation

The relationships between hillslope form and the development of zones of hillslope saturation are central to studies of process hydrology. These relationships are examined in two small, steep (30°) first-order basins under forest and pasture vegetation. Hillslope saturation patterns are described using a probability analysis of observations of peak water table elevations made at individual sites in each basin. In both first-order basins, the effects of topographic convergence on the spatial variability of hillslope saturation zones are most pronounced on shallow slopes and in steeper areas where hillslope form is strongly concave. The implications of vegetation type on the formation and variation of hillslope saturation are also examined. The differences in hillslope saturation patterns in each vegetation type are attributed primarily to differences in net rainfall and not soil physical properties or other features of topography.     

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.     

Analysis of the contributions of topographic,soil, and vegetation features on the spatial distributions of surface soil moisture in a steep natural forested headwater catchment

Surface soil moisture has been extensively studied for various land uses and landforms. Although many studies have reported potential factors that control surface soil moisture over space or time, the findings have not always been consistent, indicating a need for identification of the main factors. This study focused on the static controls of topographic, soil, and vegetation features on surface soil moisture in a steep natural forested headwater catchment consisting of three hillslope units of a gully area, side slope, and valley‐head slope. Using a simple correlation analysis to investigate the effects of the static factors on surface soil moisture at depths of 0–20 cm at 470 points in 13 surveys, we addressed the characteristics of surface soil moisture and its main controlling factors. The results indicated that the mean of surface soil moisture was in the decreasing order of gully area > valley‐head slope > side slope. The relationship between the mean and standard deviation of surface soil moisture showed a convex‐upward shape in the headwater catchment, a negative curvilinear shape in the gully area, and positive curvilinear shapes at the side and valley‐head slopes. At the headwater catchment and valley‐head slope, positive contributions of soil porosity and negative contributions of slope gradient and saturated hydraulic conductivity were the main controlling factors of surface soil moisture under wetter conditions, whereas positive contributions of topographic wetness index and negative contributions of vegetation density were the main controlling factors of surface soil moisture under drier conditions. At the side slope underlain by fractured bedrocks, only saturated hydraulic conductivity and vegetation density were observed to be the controlling factors. Surface soil moisture in the gully area was mainly affected by runoff rather than were static features. Thus, using hillslope units is effective for approximately estimating the hydrological behaviours of surface moisture on a larger scale, whereas dependency between the main static factors and moisture conditions is helpful for estimating the spatial distributions of surface moisture on a smaller scale.     

Erosional unloading,hillslope geometry,and the height of the Cascade Range,Washington State,USA

Peaks in the Cascade Range in northern Washington State are on average ～800 m higher than in southern Washington. The influences of differential valley excavation and variations in hillslope length and average slope on these altitudinal trends were tested using a 3‐dimensional model for isostatic rock uplift and calculations of hillslope length and slope respectively. The magnitude of isostatic peak uplift calculated by the model is highly dependent on the flexural rigidity (D) and the related effective elastic thickness (T_e) of the crust of this region. Crustal rigidity was constrained using published estimates and by estimating the depth of the seismogenic zone in the area (D > 1 × 10²³ Nm and T_e > 24 km). With these constraints, isostatic compensation due to differential erosion added < 700 m and 300 m, or < 25% overall, of height to peaks in the northern and southern Cascades, respectively. Deeper valley incision in the northern Cascades accounts for < 300 m of the 800 m difference in peak altitudes between north and south. Similarly, variation in valley spacing and slope account for < 350 m of the difference in mean altitude between northern and southern regions. Hence, at least several hundred m difference in altitude between the northern and southern regions of the Cascades in Washington must be due to tectonic, geologic, or geophysical factors rather than surficial and geomorphic effects like isostatic response to valley incision and hillslope geometry. Copyright © 2009 John Wiley & Sons, Ltd.