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.     

Strong slope-aspect control of regolith thickness by bedrock foliation

The porous near-surface layer of the Earth's crust – the critical zone – constitutes a vital reservoir of water for ecosystems, provides baseflow to streams, guides recharge to deep aquifers, filters contaminants from groundwater, and regulates the long-term evolution of landscapes. Recent work suggests that the controls on regolith thickness include climate, tectonics, lithology, and vegetation. However, the relative paucity of observations of regolith structure and properties at landscape scales means that theoretical models of critical zone structure are incompletely tested. Here we present seismic refraction and electrical resistivity surveys that thoroughly characterize subsurface structure in a small catchment in the Santa Catalina Mountains, Arizona, USA, where slope-aspect effects on regolith structure are expected based on differences in vegetation. Our results show a stark contrast in physical properties and inferred regolith thickness on opposing slopes, but in the opposite sense of that expected from environmental models and observed vegetation patterns. Although vegetation (as expressed by normalized difference vegetation index [NDVI]) is denser on the north-facing slope, regolith on the south-facing slope is four times thicker (as indicated by lower seismic velocities and resistivities). This contrast cannot be explained by variations in topographic stress or conventional hillslope morphology models. Instead, regolith thickness appears to be controlled by metamorphic foliation: regolith is thicker where foliation dips into the topography, and thinner where foliation is nearly parallel to the surface. We hypothesize that, in this catchment, hydraulic conductivity and infiltration capacity control weathering: infiltration is hindered and regolith is thin where foliation is parallel to the surface topography, whereas water infiltrates deeper and regolith is thicker where foliation intersects topography at a substantial angle. These results suggest that bedrock foliation, and perhaps by extension sedimentary layering, can control regolith thickness and must be accounted for in models of critical zone development. © 2020 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.     

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.     

The measurement and modelling of rill erosion at angle of repose slopes in mine spoil

The process of rill erosion causes significant amounts of sediment to be moved in both undisturbed and disturbed environments and can be a significant issue for agriculture as well as mining lands. Rills also often develop very quickly (from a single rainfall event to a season) and can develop into gullies if sufficient runoff is available to continue their development. This study examines the ability of a terrestrial laser scanner to quantify rills that have developed on fresh and homogeneous mine spoil on an angle of repose slope. It also examines the ability of the SIBERIA erosion model to simulate the rill's spatial and temporal behaviour. While there has been considerable work done examining rill erosion on rehabilitated mine sites and agricultural fields, little work has been done to examine rill development at angle of repose sites. Results show that while the overall hillslope morphology was captured by the laser scanner, with the morphology of the rills being broadly captured, the characteristics of the rills were not well defined. The digital elevation model created by the laser scanner failed to capture the rill thalwegs and tops of the banks, therefore delineating a series of ill defined longitudinal downslope depressions. These results demonstrate that an even greater density of points is needed to capture sufficient rill morphology. Nevertheless, SIBERIA simulations of the hillslope demonstrated that the model was able to capture rill behaviour in both space and time when correct model parameters were used. This result provides confidence in the SIBERIA model and its parameterization. The results demonstrate the sensitivity of the model to changes in parameters and the importance of the calibration process. Copyright © 2007 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.     

The effect of landform variation on vegetation patterning and related sediment dynamics

Semi‐arid ecosystems are often spatially self‐organized in typical patterns of vegetation bands with high plant cover interspersed with bare soil areas, also known as ‘tiger bush’. In modelling studies, most often, straight planar slopes were used to analyse vegetation patterning. The effect of slope steepness has been investigated widely, and some studies investigated the effects of microtopography and hillslope orientation. However, at the larger catchment scale, the overall form of the landscape may affect vegetation patterning and these more complex landscapes are much more prevalent than straight slopes. Hence, our objective was to determine the effect of landform variation on vegetation patterning and sediment dynamics. We linked two well‐established models that simulate (a) plant growth, death and dispersal of vegetation, and (b) erosion and sedimentation dynamics. The model was tested on a straight planar hillslope and then applied to (i) a set of simple synthetic topographies with varying curvature and (ii) three more complex, real‐world landscapes of distinct morphology. Results show banded vegetation patterning on all synthetic topographies, always perpendicular to the slope gradient. Interestingly, we also found that movement of bands – a debated phenomenon – seems to be dependent on curvature. Vegetation banding was simulated on the slopes of the alluvial fan and along the valley slopes of the dissected and rolling landscapes. In all landscapes, local valleys developed a full vegetation cover induced by water concentration, which is consistent with observations worldwide. Finally, banded vegetation patterns were found to reduce erosion significantly as compared to other vegetation configurations. © 2018 The Authors. Earth Surface Processes and Landforms published by John Wiley & Sons Ltd.     

Hydro-geomorphic behavior of contrasting tropical landscapes and critical zone response to changing climate

The potentially important influence of climate change on landscape evolution and on critical zone processes is not sufficiently understood. The relative contribution of hydro-climatic factors on hillslope erosion rates may significantly vary with topography at the watershed scale. The objective of this study is to quantify the hydro-geomorphic behavior of two contrasting landscapes in response to different climate change scenarios in the Luquillo Critical Zone Observatory, a site of particular geomorphological interest, in terms of hillslope erosion and rainfall-triggered landslides. We investigate the extent to which hillslope erosion and landslide occurrence remain relatively invariant with future hydro-climatic perturbations. The adjacent Mameyes and Icacos watersheds are studied, which are underlain by contrasting lithologies. A high resolution coupled hydro-geomorphic model based on tRIBS (Triangulated Irregular Network-based Real-time Integrated Basin Simulator) is used. Observations of landslide activity and hillslope erosion are used to evaluate the model performance. The process-based model quantifies feedbacks among different hydrologic processes, landslide occurrence, and topsoil erosion and deposition. Simulations suggest that the propensity for landslide occurrence in the Luquillo Mountains is controlled by tropical storms, subsurface water flow, and by non-climatic factors, and is expected to remain significant through 2099. The Icacos watershed, which is underlain by quartz diorite, is dominated by relatively large landslides. The relative frequency of smaller landslides is higher at the Mameyes watershed, which is underlain by volcaniclastic rock. While projections of precipitation decrease at the study site may lead to moderate decline in hillslope erosion rates, the simulated erosional potential of the two diverse landscapes likely remains significant. © 2018 John Wiley & Sons, Ltd.     

Subsurface lateral flow generation in aspen and conifer‐dominated hillslopes of a first order catchment in northern Utah

Mountain headwater catchments in the semi‐arid Intermountain West are important sources of surface water because these high elevations receive more precipitation than neighboring lowlands. This study examined subsurface runoff in two hillslopes, one aspen dominated, the other conifer dominated, adjacent to a first order stream in snow‐driven northern Utah. Snow accumulation, soil moisture, trenchflow and streamflow were examined in hillslopes and their adjacent stream. Snow water equivalents (SWEs) were greater under aspen stands compared to conifer, the difference increasing with higher annual precipitation. Semi‐variograms of shallow spatial soil moisture patterns and transects of continuous soil moisture showed no increase in soil moisture downslope, suggesting the absence of subsurface flow in shallow (～12 cm) soil layers of either vegetation type. However, a clear threshold relationship between soil moisture and streamflow indicated hillslope–stream connectivity, deeper within the soil profile. Subsurface flow was detected at ～50 cm depth, which was sustained for longer in the conifer hillslope. Soil profiles under the two vegetation types varied, with deep aspen soils having greater water storage capacity than shallow rocky conifer soils. Though SWEs were less under the conifers, the soil profile had less water storage capacity and produced more subsurface lateral flow during the spring snowmelt. Copyright © 2010 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.     

Controls of preferential orientation of earthquake- and rainfall-triggered landslides in Taiwan's orogenic mountain belt

Landsliding usually occurs on specific hillslope aspect, which may reflect the control of specific geo-environmental factors, triggering factors, or their interaction. To explore this notion, this study used island-wide landslide inventories of the Chi-Chi earthquake in 1999 (M_W = 7.6) and Typhoon Morakot in 2009 in Taiwan to investigate the preferential orientation of landslides and the controls of landslide triggers and geological settings. The results showed two patterns. The orientations of earthquake-triggered landslides were toward the aspect facing away from the epicenter in areas with peak ground acceleration (PGA) ≥ 0.6 g and landslide ratio ≥ 1%, suggesting that the orientations were controlled by seismic wave propagation. Rainfall-triggered landslides tended to occur on dip slopes, instead of the windward slopes, suggesting that geological settings were a more effective control of the mass wasting processes on hillslope scale than the rainfall condition. This study highlights the importance of the endogenic processes, namely seismic wave and geological settings, on the predesigned orientation of landslides triggered by either earthquake or rainfall, which can in turn improve our knowledge of landscape evolution and landslide prediction. © 2019 John Wiley & Sons, Ltd.     

Spatial structures of stream and hillslope drainage networks following gully erosion after wildfire

The drainage networks of catchment areas burned by wildfire were analysed at several scales. The smallest scale (1–1000 m²) representative of hillslopes, and the small scale (1000 m² to 1 km²), representative of small catchments, were characterized by the analysis of field measurements. The large scale (1–1000 km²), representative of perennial stream networks, was derived from a 30‐m digital elevation model and analysed by computer analysis. Scaling laws used to describe large‐scale drainage networks could be extrapolated to the small scale but could not describe the smallest scale of drainage structures observed in the hillslope region. The hillslope drainage network appears to have a second‐order effect that reduces the number of order 1 and order 2 streams predicted by the large‐scale channel structure. This network comprises two spatial patterns of rills with width‐to‐depth ratios typically less than 10. One pattern is parallel rills draining nearly planar hillslope surfaces, and the other pattern is three to six converging rills draining the critical source area uphill from an order 1 channel head. The magnitude of this critical area depends on infiltration, hillslope roughness and critical shear stress for erosion of sediment, all of which can be substantially altered by wildfire. Order 1 and 2 streams were found to constitute the interface region, which is altered by a disturbance, like wildfire, from subtle unchannelized drainages in unburned catchments to incised drainages. These drainages are characterized by gullies also with width‐to‐depth ratios typically less than 10 in burned catchments. The regions (hillslope, interface and channel) had different drainage network structures to collect and transfer water and sediment. Copyright © 2005 John Wiley & Sons, Ltd.     

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

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

Linking landslides,hillslope erosion,and landscape evolution

This Virtual Issue highlights 10 recent innovative, unconventional, or otherwise significant contributions to Earth Surface Processes and Landforms that help advance the state‐of‐the‐art in research on linkages between landslides, hillslope erosion, and landscape evolution. The selected studies address this feedback within a temporal spectrum that ranges from the event to the millennial scale, thus underscoring the importance of detailed field observations, high‐resolution digital topographic data, and geochronological methods for increasing our capability of quantifying landslide processes and hillslope erosion. Copyright © 2009 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.     

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.     

DRAINAGE NETWORKS AFTER WILDFIRE

Predicting runoff and erosion from watersheds burned by wildfires requires an understanding of the three-dimensional structure of both hillslope and channel drainage networks. We investigate the small- and large-scale structures of drainage networks using field studies and computer analysis of 30- m digital elevation model. Topologic variables were derived from a composite 30-m DEM, which included 14 order 6 watersheds within the Pikes Peak batholith. Both topologic and hydraulic variables were measured in the field in two smaller burned watersheds （3.7 and 7.0 hectares） located within one of the order 6 watersheds burned by the 1996 Buffalo Creek Fire in Central Colorado. Horton ratios of topologic variables （stream number, drainage area, stream length, and stream slope） for small-scale and large-scale watersheds are shown to scale geometrically with stream order （i.e., to be scale invariant）. However, the ratios derived for the large-scale drainage networks could not be used to predict the rill and gully drainage network structure. Hydraulic variables （width, depth, cross- sectional area, and bed roughness） for small-scale drainage networks were found to be scale invariant across 3 to 4 stream orders. The relation between hydraulic radius and cross-sectional area is similar for fills and gullies, suggesting that their geometry can be treated similarly in hydraulic modeling. Additionally, the rills and gullies have relatively small width-to-depth ratios, implying sidewall friction may be important to the erosion and evolutionary process relative to main stem channels.     

Soil water dynamics under different forest vegetation cover: Implications for hillslope stability

Though it is well known that vegetation affects the water balance of soils through canopy interception and evapotranspiration, its hydrological contribution to soil hydrology and stability is yet to be fully quantified. To improve understanding of this hydrological process, soil water dynamics have been monitored at three adjacent hillslopes with different vegetation covers (deciduous tree cover, coniferous tree cover, and grass cover), for nine months from December 2014 to September 2015. The monitored soil moisture values were translated into soil matric suction (SMS) values to facilitate the analysis of hillslope stability. Our observations showed significant seasonal variations in SMS for each vegetation cover condition. However, a significant difference between different vegetation covers was only evident during the winter season where the mean SMS under coniferous tree cover (83.6 kPa) was significantly greater than that under grass cover (41 kPa). The hydrological reinforcing contribution due to matric suction was highest for the hillslope with coniferous tree cover, while the hillslope with deciduous tree cover was second and the hillslope with grass cover was third. The greatest contributions for all cover types were during the summer season. During the winter season, the wettest period of the monitoring study, the additional hydrological reinforcing contributions provided by the deciduous tree cover (1.5 to 6.5 kPa) or the grass cover (0.9 to 5.4 kPa) were insufficient to avoid potential slope failure conditions. However, the additional hydrological reinforcing contribution from the coniferous tree cover (5.8 to 10.4 kPa) was sufficient to provide potentially stable hillslope conditions during the winter season. Our study clearly suggests that during the winter season the hydrological effects from both deciduous tree and grass covers are insufficient to promote slope stability, while the hydrological reinforcing effects from the coniferous tree cover are sufficient even during the winter season. Copyright © 2018 John Wiley & Sons, Ltd.     

Using stable isotopes to identify major flow pathways in a permafrost influenced alpine meadow hillslope during summer rainfall period

Global warming has leaded to permafrost degradation, with potential impacts on the runoff generation processes of permafrost influenced alpine meadow hillslope. Stable isotopes have the potential to trace the complex runoff generation processes. In this study, precipitation, hillslope surface and subsurface runoff, stream water, and mobile soil water (MSW) at different hillslope positions and depths were collected during the summer rainfall period to analyse the major flow pathway based on stable isotopic signatures. The results indicated that (a) compared with precipitation, the δ²H values of MSW showed little temporal variation but strong heterogeneity with enriched isotopic ratios at lower hillslope positions and in deeper soil layers. (b) The δ²H values of middle-slope surface runoff and shallow subsurface flow were similar to those of precipitation and MSW of the same soil layer, respectively. (c) Middle-slope shallow subsurface flow was the major flow pathway of the permafrost influenced alpine meadow hillslope, which turned into surface runoff at the riparian zone before contributing to the streamflow. (d) The slight variation of δ²H values in stream water was shown to be related to mixing processes of new water (precipitation, 2%) and old water (middle-slope shallow subsurface flow, 98%) in the highly transmissive shallow thawed soil layers. It was inferred that supra-permafrost water levels would be lowered to a less conductive, deeper soil layer under further warming and thawing permafrost, which would result in a declined streamflow and delayed runoff peak. This study explained the “rapid mobilization of old water” paradox in permafrost influenced alpine meadow hillslope and improved our understanding of permafrost hillslope hydrology in alpine regions.     

Variable hillslope-channel coupling and channel characteristics of forested mountain streams in glaciated landscapes

Channel morphology of forested, mountain streams in glaciated landscapes is regulated by a complex suite of processes, and remains difficult to predict. Here, we analyze models of channel geometry against a comprehensive field dataset collected in two previously glaciated basins in Haida Gwaii, B.C., to explore the influence of variable hillslope–channel coupling imposed by the glacial legacy on channel form. Our objective is to better understand the relation between hillslope–channel coupling and stream character within glaciated basins. We find that the glacial legacy on landscape structure is characterized by relatively large spatial variation in hillslope–channel coupling. Spatial differences in coupling influence the frequency and magnitude of coarse sediment and woody material delivery to the channel network. Analyses using a model for channel gradient and multiple models for width and depth show that hillslope–channel coupling and high wood loading induce deviations from standard downstream predictions for all three variables in the study basins. Examination of model residuals using Boosted Regression Trees and nine additional channel variables indicates that ~10 to ~40% of residual variance can be explained by logjam variables, ~15–40% by the degree of hillslope–channel coupling, and 10–20% by proximity to slope failures. These results indicate that channel classification systems incorporating hillslope–channel coupling, and, indirectly, the catchment glacial legacy, may present a more complete understanding of mountain channels. From these results, we propose a conceptual framework which describes the linkages between landscape history, hillslope–channel coupling, and channel form. © 2018 John Wiley & Sons, Ltd.