The geormorphology of subalpine snow avalanche runout zones: San juan mountains,colorado

A comparison was made between the distal ends of twenty-two avalanche and fifteen non-avalanche slopes in the San Juan Mountains of Colorado, U.S.A. All slopes occurred in the subalpine zone. Six characteristics were used for analysis: type of slope, surface material, longitudinal profile, perched debris or debris tails, avalanche impact on opposite valley wall, and transverse profile. Both fan and roadbank avalanche slope types were found along with the non-avalanche slopes. Almost all slopes were turf covered rather than talus since the work was done below treeline. Twenty avalanche slopes had a distinctive concave longitudinal profile. Little debris of any kind was found since the slopes were in an area of insignificant amounts of detritus. Many of the larger and two of the smaller avalanche slopes showed evidence of impact upon the opposite slope. Eighteen of the avalanche slopes had convex transverse profiles.     

Talus slope geomorphology investigated at multiple time scales from high-resolution topographic surveys and historical aerial photographs (Sanetsch Pass,Switzerland)

Talus slopes are common places for debris storage in high-mountain environments and form an important step in the alpine sediment cascade. To understand slope instabilities and sediment transfers, detailed investigations of talus slope geomorphology are needed. Therefore, this study presents a detailed analysis of a talus slope on Col du Sanetsch (Swiss Alps), which is investigated at multiple time scales using high-resolution topographic (HRT) surveys and historical aerial photographs. HRT surveys were collected during three consecutive summers (2017–2019), using uncrewed aerial vehicle (UAV) and terrestrial laser scanning (TLS) measurements. To date, very few studies exist that use HRT methods on talus slopes, especially to the extent of our study area (2 km²). Data acquisition from ground control and in situ field observations is challenging on a talus slope due to the steep terrain (30–37°) and high surface roughness. This results in a poor spatial distribution of ground control points (GCPs), causing unwanted deformation of up to 2 m in the gathered UAV-derived HRT data. The co-alignment of UAV imagery from different survey dates improved this deformation significantly, as validated by the TLS data. Sediment transfer is dominated by small-scale but widespread snow push processes. Pre-existing debris flow channels are prone to erosion and redeposition of material within the channel. A debris flow event of high magnitude occurred in the summer of 2019, as a result of several convective thunderstorms. While low-magnitude (<5,000 m³) debris flow events are frequent throughout the historical record with a return period of 10–20 years, this 2019 event exceeded all historical debris flow events since 1946 in both extent and volume. Future climate predictions show an increase of such intense precipitation events in the region, potentially altering the frequency of debris flows in the study area and changing the dominant geomorphic process which are active on such talus slopes. © 2020 John Wiley & Sons, Ltd.     

Physical modelling of rainfall‐ and snowmelt‐induced erosion of stony slope underlain by permafrost

Physical modelling experiments have been carried out in a cold room to test on a small scale, the effects of water supply during the thaw of an experimental slope with permafrost. Permafrost was maintained at depth and a thin active layer was frozen and thawed from the surface. Data from the experiments relate to two different conditions, first with moderate rainfall, and second with heavy rainfall during the thaw period. When moderate rainfall is applied during thaw phases, the experimental slope is slightly degraded. At the scale of the experiment, erosion processes involve frost jacking of the coarse blocks, frost creep and gelifluction that induce slow and gradual down slope displacements of the active layer, but also small landslides leading to large but slow mass movements with short displacements. Changes in experimental slope morphology are marked by the initiation of a small‐scale drainage network and the development of a little crest line which shows a progressive upslope migration. With such boundary conditions, there is not enough water supply to evacuate downslope the whole of the eroded material and a topographic smoothing is observed. When heavy rainfall is applied during thaw periods, rapid mass wasting (small mud‐flows and debris flows) become prominent. Slope failures are largely controlled by the water saturation of the active layer and by the occurrence of steeper slopes. At the scale of the experiment, rates of erosion and maximum incision increase by about 100% leading to significant slope degradation with marked and specific scars comparable to gullying. These morphological changes are dependant on both the size and the frequency of catastrophic events. These experiments provide detailed data that could improve the knowledge of the physical parameters that control the initiation, at a small‐scale, of erosion processes on periglacial slopes with a thin active layer and/or with thin cover of mobilizable slope deposits. Copyright © 2010 John Wiley & Sons, Ltd.     

Offshore aeolian sediment transport across a human‐modified foredune

Concepts derived from previous studies of offshore winds on natural dunes are evaluated on a dune maintained for shore protection during three offshore wind events. The potential for offshore winds to form a lee‐side eddy on the backshore or transfer sediment from the dune and berm crest to the water are evaluated, as are differences in wind speed and sediment transport on the dune crest, berm crest and a pedestrian access gap. The dune is 18–20 m wide near the base and has a crest 4.5 m above backshore elevation. Two sand‐trapping fences facilitate accretion. Data were obtained from wind vanes on the crest and lee of the dune and anemometers and sand traps placed across the dune, on the beach berm crest and in the access gap. Mean wind direction above the dune crest varied from 11 to 3 deg from shore normal. No persistent recirculation eddy occurred on the 12 deg seaward slope. Wind speed on the berm crest was 85–89% of speed at the dune crest, but rates of sediment transport were 2.27 times greater during the strongest winds, indicating that a wide beach overcomes the transport limitation of a dune barrier. Limited transport on the seaward dune ramp indicates that losses to the water are mostly from the backshore, not the dune. The seaward slope gains sand from the landward slope and dune crest. Sand fences causing accretion on the dune ramp during onshore winds lower the seaward slope and reduce the likelihood of detached flows during offshore winds. Transport rates are higher in access gaps than on the dune crest despite lower wind speeds because of flatter slopes and absence of vegetation. Transport rates across dunes and through gaps can be reduced using vegetation and raised walkover structures. Copyright © 2017 John Wiley & Sons, Ltd.     

Stability of supraglacial debris

Rock debris on the surface of ablating glaciers is not static, and is often transported across the ice surface as relief evolves during melt. This supraglacial debris transport has a strong influence on the spatial distribution of melt, and is implicated in the formation of hummocky glacial topography in deglaciated terrain. Furthermore, as ice‐dammed lakes and ice‐cored slopes become increasingly common in deglaciating watersheds, there is rising concern about hazards to humans and infrastructure posed by mass‐wasting of ice‐cored debris. The existing quantitative framework for describing these debris transport processes is limited, making it difficult to account for transport in mass balance, hazard assessment, and landscape development models. This paper develops a theoretical framework for assessing slope stability and gravitational mass transport in a debris‐covered ice setting. Excess water pressure at the interface between ablating ice and lowering debris is computed by combining Darcy's law with a meltwater balance. A limit‐equilibrium slope stability analysis is then applied to hypothetical debris layers with end‐member moisture conditions derived from a downslope meltwater balance that includes production and seepage. The resulting model system constrains maximum stable slope angles and lengths that vary with debris texture, thickness, and the rate of meltwater production. Model predictions are compared with field observations and with digital elevation model (DEM)‐derived terrain metrics from two modern debris‐covered glaciers on Mount Rainier, USA. Copyright © 2017 John Wiley & Sons, Ltd.     

Seepage erosion and its implication to the formation of amphitheatre valley heads: A case study at Obara,Japan

Seepage erosion was investigated in an amphitheatre with a semicircular valley head, steep slopes, and a flat bottom developed in granodiorite hills at Obara, Aichi prefecture, Japan. A high sediment yield occurred where the measuring sites were located at the base of the landslide debris in the base of the convex slopes, whereas sediment outflows were small where the measuring sites were located at the base of the strong convex slopes. This implies that the seepage erosion was an effective agent for removal of debris deposited at the base of the slope. Small landslides can be found at the lower slopes within the area of the observed amphitheatre. The slope stability analysis and subsurface water observation of the lower slope suggest that the small landslides in this amphitheatre are due to over-steepened slopes, and relatively insensitive to subsurface water status. Colluvium in the flat valley bottom thinly covers the bedrock surface. Therefore the topography of the amphitheatre was found to be formed by parallel retreat of slopes by the repetition of basal seepage erosion and subsequent small landslides.     

The effects of slope shape and polyacrylamide application on runoff,erosion and nutrient loss from hillslopes under simulated rainfall

Natural hillslopes are mostly composed of complex slope shapes, which significantly affect soil erosion. However, existing studies have mainly focused on uniform slopes to simplify complex hillslopes, and the mechanisms responsible for the influence of slope shape on soil and nutrient losses are still not well understood, especially in the application of soil improvers to reduce soil loss. To investigate the effects of slope shape and polyacrylamide (PAM) application on runoff, soil erosion and nutrient loss, this study conducted artificial field rainfall experiments involving two PAM application rates and nine slope shapes. The results indicate that the average amount of soil loss from convex slopes was 1.5 and 1.3 times greater than that from concave and uniform slopes, respectively, and the average amount of ammonia nitrogen loss and phosphate loss increased by 24.0%–58.6%. Soil and nutrient losses increased as the convexity of the convex slopes increased. For runoff, there was little difference between concave and convex slopes, but the runoff amount for both slopes was greater than that for uniform slopes. After PAM application, the soil loss decreased by more than 90%, and the nutrient loss decreased by 28.2%–68.1%. The application of PAM was most effective in reducing soil erosion and nutrient loss from convex slopes, and it is recommended to appropriately increase the PAM application rate for convex slopes. A strong linear relationship between ammonia nitrogen and phosphate concentrations and sediment concentrations was found in the runoff on slopes with no PAM application. However, this linear relationship weakened for slopes with PAM application. The findings of this study may be valuable for optimizing nonpoint source pollution management in basins.     

The morphology and evolution of the Stromboli 2002–2003 lava flow field: an example of a basaltic flow field emplaced on a steep slope

The use of a hand-held thermal camera during the 2002–2003 Stromboli effusive eruption proved essential in tracking the development of flow field structures and in measuring related eruption parameters, such as the number of active vents and flow lengths. The steep underlying slope on which the flow field was emplaced resulted in a characteristic flow field morphology. This comprised a proximal shield, where flow stacking and inflation caused piling up of lava on the relatively flat ground of the vent zone, that fed a medial–distal lava flow field. This zone was characterized by the formation of lava tubes and tumuli forming a complex network of tumuli and flows linked by tubes. Most of the flow field was emplaced on extremely steep slopes and this had two effects. It caused flows to slide, as well as flow, and flow fronts to fail frequently, persistent flow front crumbling resulted in the production of an extensive debris field. Channel-fed flows were also characterized by development of excavated debris levees in this zone (Calvari et al. 2005). Collapse of lava flow fronts and inflation of the upper proximal lava shield made volume calculation very difficult. Comparison of the final field volume with that expecta by integrating the lava effusion rates through time suggests a loss of ~70% erupted lava by flow front crumbling and accumulation as debris flows below sea level. Derived relationships between effusion rate, flow length, and number of active vents showed systematic and correlated variations with time where spreading of volume between numerous flows caused an otherwise good correlation between effusion rate, flow length to break down. Observations collected during this eruption are useful in helping to understand lava flow processes on steep slopes, as well as in interpreting old lava–debris sequences found in other steep-sided volcanoes subject to effusive activity.     

Grain size-dependent 10Be concentrations in alluvial stream sediment of the Huasco Valley,a semi-arid Andes region

Terrestrial cosmogenic nuclide concentrations in sediment are used to quantify mean denudation rates in catchments. This article explores the differences between the ¹⁰Be concentration in fine (sand) and in coarse (1–3 or 5–10 cm pebbles) river sediment. Sand and pebbles were sampled at four locations in the Huasco Valley, in the arid Chilean Andes. Sand has ¹⁰Be concentrations between 4.8 and 8.3·10⁵ at g⁻¹, while pebbles have smaller concentrations between 2.2 and 3.3·10⁵ at g⁻¹. It appears that the different concentrations, systematically measured between sand and pebbles, are the result of different denudation rates, linked with the geomorphologic processes that originated them. We propose that the ¹⁰Be concentrations in sand are determined by the mean denudation rate of all of the geomorphologic processes taking place in the catchment, including debris flow processes as well as slower processes such as hill slope diffusion. In contrast, the concentrations in pebbles are probably related to debris flows occurring in steep slopes. The mean denudation rates calculated in the catchment are between 30 and 50 m/Myr, while the denudation rates associated with debris flow are between 59 and 81 m/Myr. These denudation rates are consistent with those calculated using different methods, such as measuring eroded volumes.     

An enhanced temperature index model for debris-covered glaciers accounting for thickness effect

Debris-covered glaciers are increasingly studied because it is assumed that debris cover extent and thickness could increase in a warming climate, with more regular rockfalls from the surrounding slopes and more englacial melt-out material. Debris energy-balance models have been developed to account for the melt rate enhancement/reduction due to a thin/thick debris layer, respectively. However, such models require a large amount of input data that are not often available, especially in remote mountain areas such as the Himalaya, and can be difficult to extrapolate. Due to their lower data requirements, empirical models have been used extensively in clean glacier melt modelling. For debris-covered glaciers, however, they generally simplify the debris effect by using a single melt-reduction factor which does not account for the influence of varying debris thickness on melt and prescribe a constant reduction for the entire melt across a glacier.In this paper, we present a new temperature-index model that accounts for debris thickness in the computation of melt rates at the debris-ice interface. The model empirical parameters are optimized at the point scale for varying debris thicknesses against melt rates simulated by a physically-based debris energy balance model. The latter is validated against ablation stake readings and surface temperature measurements. Each parameter is then related to a plausible set of debris thickness values to provide a general and transferable parameterization. We develop the model on Miage Glacier, Italy, and then test its transferability on Haut Glacier d’Arolla, Switzerland.The performance of the new debris temperature-index (DETI) model in simulating the glacier melt rate at the point scale is comparable to the one of the physically based approach, and the definition of model parameters as a function of debris thickness allows the simulation of the nonlinear relationship of melt rate to debris thickness, summarised by the Østrem curve. Its large number of parameters might be a limitation, but we show that the model is transferable in time and space to a second glacier with little loss of performance. We thus suggest that the new DETI model can be included in continuous mass balance models of debris-covered glaciers, because of its limited data requirements. As such, we expect its application to lead to an improvement in simulations of the debris-covered glacier response to climate in comparison with models that simply recalibrate empirical parameters to prescribe a constant across glacier reduction in melt.     

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.     

Arctic rockwall retreat rates estimated using laboratory‐calibrated ERT measurements of talus cones in Longyeardalen,Svalbard

Holocene rockwall retreat rates quantify integral values of rock slope erosion and talus cone evolution. Here we investigate Holocene rockwall retreat of exposed arctic sandstone cliffs in Longyeardalen, central Svalbard and apply laboratory‐calibrated electrical resistivity tomography (ERT) to determine talus sediment thickness. Temperature–resistivity functions of two sandstone samples are measured in the laboratory and compared with borehole temperatures from the talus slope. The resistivity of the higher and lower‐porosity sandstone at relevant borehole permafrost temperatures defines a threshold range that accounts for the lithological variability of the dominant bedrock and debris material. This helps to estimate the depth of the transition from higher resistivities of ice‐rich debris to lower resistivities of frozen bedrock in the six ERT transects. The depth of the debris–bedrock transition in ERT profiles is confirmed by a pronounced apparent resistivity gradient in the raw data plotted versus depth of investigation. High‐resolution LiDAR‐scanning and ERT subsurface information were collated in a GIS to interpolate the bedrock surface and to calculate the sediment volume of the talus cones. The resulting volumes were referenced to source areas to calculate rockwall retreat rates. The rock mass strength was estimated for the source areas. The integral rockwall retreat rates range from 0.33 to 1.96 mm yr^–1, and are among the highest rockwall retreat rates measured in arctic environments, presumably modulated by harsh environmental forcing on a porous sandstone rock cliff with a comparatively low rock mass strength. Here, we show the potential of laboratory‐calibrated ERT to provide accurate estimates of rockwall retreat rates even in ice‐rich permafrost talus slopes. Copyright © 2012 John Wiley & Sons, Ltd.     

Processes at the margins of supraglacial debris cover: Quantifying dirty ice ablation and debris redistribution

Current glacier ablation models have difficulty simulating the high-melt transition zone between clean and debris-covered ice. In this zone, thin debris cover is thought to increase ablation compared to clean ice, but often this cover is patchy rather than continuous. There is a need to understand ablation and debris dynamics in this transition zone to improve the accuracy of ablation models and the predictions of future debris cover extent. To quantify the ablation of partially debris-covered ice (or ‘dirty ice’), a high-resolution, spatially continuous ablation map was created from repeat unmanned aerial systems surveys, corrected for glacier flow in a novel way using on-glacier ablation stakes. Surprisingly, ablation is similar (range ~ 5 mm w.e. per day) across a wide range of percentage debris covers (~ 30–80%) due to the opposing effects of a positive correlation between percentage debris cover and clast size, countered by a negative correlation with albedo. Once debris cover becomes continuous, ablation is significantly reduced (by 61.6% compared to a partial debris cover), and there is some evidence that the cleanest ice (<~ 15% debris cover) has a lower ablation than dirty ice (by 3.7%). High-resolution feature tracking of clast movement revealed a strong modal clast velocity where debris was continuous, indicating that debris moves by creep down moraine slopes, in turn promoting debris cover growth at the slope toe. However, not all slope margins gain debris due to the removal of clasts by supraglacial streams. Clast velocities in the dirty ice area were twice as fast as clasts within the continuously debris-covered area, as clasts moved by sliding off their boulder tables. These new quantitative insights into the interplay between debris cover characteristics and ablation can be used to improve the treatment of dirty ice in ablation models, in turn improving estimates of glacial meltwater production. © 2020 The Authors. Earth Surface Processes and Landforms published by John Wiley & Sons Ltd     

Effect of basin physical characteristics on solute fluxes in nine alpine/subalpine basins,Colorado, USA

Alpine/subalpine basins may exhibit substantial variability in solute fluxes despite many apparent similarities in basin characteristics. An evaluation of controls on spatial patterns in solute fluxes may allow development of predictive tools for assessing basin sensitivity to outside perturbations such as climate change or deposition of atmospheric pollutants. Relationships between basin physical characteristics, determined from geographical information system (GIS) tools, and solute fluxes and mineral weathering rates were explored for nine alpine/subalpine basins in Rocky Mountain National Park, Colorado, using correlation analyses for 1993 and 1994 data. Stream‐water nitrate fluxes were correlated positively with basin characteristics associated with the talus environment; i.e., the fractional amounts of steep slopes (≥ 30°), unvegetated terrain and young debris (primarily Holocene till) in the basins, and were correlated negatively with fractional amounts of subalpine meadow terrain. Correlations with nitrate indicate the importance of the talus environment in promoting nitrate flux and the mitigating effect of areas with established vegetation, such as subalpine meadows. Total mineral weathering rates for the basins ranged from about 300 to 600 mol ha^?1 year^?1. Oligoclase weathering accounted for 30 to 73% of the total mineral weathering flux, and was positively correlated with the amount of old debris (primarily Pleistocene glacial till) in the basins. Although calcite is found in trace amounts in bedrock, calcite weathering accounted for up to 44% of the total mineral weathering flux. Calcite was strongly correlated with steep slope, unvegetated terrain, and young debris—probably because physical weathering in steep‐gradient areas exposes fresh mineral surfaces that contain calcite for chemical weathering. Oligoclase and calcite weathering are the dominant sources of alkalinity in the basins. However, atmospherically deposited acids consume much of the alkalinity generated by weathering of calcite and other minerals in the talus environment. Published in 2001 by John Wiley & Sons, Ltd.     

Slope dynamics of Lake Albano (Rome,Italy): insights from high resolution bathymetry

New detailed data about the morphology of the submerged slopes of Lake Albano (Rome, Italy) have been collected by a sonar multibeam survey financed by the Italian Department of Civil Protection. These data allow for investigation of the subaqueous slope dynamics of the lake, which partially fills a volcanic depression, and the elucidation of the relationships between subaqueous and subaerial slope processes. Subaerial, submerged and combined subaerial/submerged landslide‐related morphologies were detected around the inner slopes of the lake. In the submerged slopes, several gravity‐induced landforms were recognized: landslide scar areas, landslide accumulations, erosional chutes and channels, block fields, isolated blocks, scarps and slope breaks. An attempt to evaluate the state of activity of the submerged slopes was carried out by taking into consideration the relative freshness of some selected landforms. Interpretation of bathymetric data, as well as direct surveys of the subaerial slopes, was used to assess the morphometric features and interpret the type of movement of the landslides. We propose a comprehensive classification based on the landslide's size and type of movement. We recognized rock fall/topples, debris flows, rock slides and slump, complex rock slides/channelled flows and debris slide and slump. The volume of the main landslides ranged between 10¹ and 10³ m³, while a few rock and debris slides have volumes ranging between 10³ and 10⁵ m³. Two large palaeo‐landslides with volumes on the order of 10⁶ m³ were identified in the southern and northern part of the lake, respectively. Velocities of the recognized landslides range from rapid to extremely rapid. Two main landslide hazard scenarios have been depicted from the results of the integrated analysis of both subaerial and submerged gravity‐induced landforms. The most hazardous scenario involves extremely rapid large volume events (>10⁶ m³) that could, if they interacted with water, induce catastrophic tsunamis. Copyright © 2009 John Wiley & Sons, Ltd.     

The structure and sedimentology of relict talus,Trotternish, northern Skye,Scotland

Sections up to 3·5 m deep cut through the upper rectilinear segment of relict, vegetated talus slopes at the foot of the Trotternish escarpment reveal stacked debris-flow deposits intercalated with occasional slopewash horizons and buried organic soils. Radiocarbon dating of buried soil horizons indicates that reworking of sediment by debris flows predates 5·9–5·6 Cal ka BP , and has been intermittently active throughout the late Holocene. Particle size analyses of 18 bulk samples from these units indicates that c. 27–30 per cent of the talus deposit is composed of fine (<2 mm) sediment. Sedimentological comparison with tills excludes a glacigenic origin for the talus debris, and the angularity of constituent clasts suggests that in situ weathering has been insignificant in generating fine material. We conclude that the fine sediment within the talus is derived primarily by granular weathering of the rockwall, with syndepositional accumulation of both fine and coarse debris, implying that c. 27–30 per cent of rockwall retreat since deglaciation reflects granular weathering rather than rockfall. The abundance of fines within the talus deposits is inferred to have been of critical importance in facilitating build-up of porewater pressures during rainstorms, leading to episodic failure and flow of debris on the upper parts of the slope. A wider implication of these findings is that the mechanical properties of talus slopes cannot be regarded as those of free-draining accumulations of coarse clasts, and that models that treat talus slopes as such have limited value in explaining their form and evolution. Our findings lend support to models that envisage the upper straight slope on talus accumulations as the product of mass-transport as well as rockfall, and indicate that episodic debris flow has been the primary agent of mass-transport at this site. © 1998 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.     

Geomorphological evolution of a debris-covered glacier surface

There exists a need to advance our understanding of debris-covered glacier surfaces over relatively short timescales due to rapid, climatically induced areal expansion of debris cover at the global scale, and the impact debris has on mass balance. We applied unpiloted aerial vehicle structure-from-motion (UAV-SfM) and digital elevation model (DEM) differencing with debris thickness and debris stability modelling to unravel the evolution of a 0.15 km² region of the debris-covered Miage Glacier, Italy, between June 2015 and July 2018. DEM differencing revealed widespread surface lowering (mean 4.1 ± 1.0 m a^-1; maximum 13.3 m a^-1). We combined elevation change data with local meteorological data and a sub-debris melt model, and used these relationships to produce high resolution, spatially distributed maps of debris thickness. These maps were differenced to explore patterns and mechanisms of debris redistribution. Median debris thicknesses ranged from 0.12 to 0.17 m and were spatially variable. We observed localized debris thinning across ice cliff faces, except those which were decaying, where debris thickened. We observed pervasive debris thinning across larger, backwasting slopes, including those bordered by supraglacial streams, as well as ingestion of debris by a newly exposed englacial conduit. Debris stability mapping showed that 18.2–26.4% of the survey area was theoretically subject to debris remobilization. By linking changes in stability to changes in debris thickness, we observed that slopes that remain stable, stabilize, or remain unstable between periods almost exclusively show net debris thickening (mean 0.07 m a^-1) whilst those which become newly unstable exhibit both debris thinning and thickening. We observe a systematic downslope increase in the rate at which debris cover thickens which can be described as a function of the topographic position index and slope gradient. Our data provide quantifiable insights into mechanisms of debris remobilization on glacier surfaces over sub-decadal timescales, and open avenues for future research to explore glacier-scale spatiotemporal patterns of debris remobilization. © 2020 The Authors. Earth Surface Processes and Landforms published by John Wiley & Sons Ltd     

Rates of rockwall slope erosion in the upper Bhagirathi catchment,Garhwal Himalaya

Rockwall slope erosion is defined for the upper Bhagirathi catchment using cosmogenic Beryllium-10 (¹⁰Be) concentrations in sediment from medial moraines on Gangotri glacier. Beryllium-10 concentrations range from 1.1 ± 0.2 to 2.7 ± 0.3 × 10⁴ at/g SiO₂, yielding rockwall slope erosion rates from 2.4 ± 0.4 to 6.9 ± 1.9 mm/a. Slope erosion rates are likely to have varied over space and time and responded to shifts in climate, geomorphic and/or tectonic regime throughout the late Quaternary. Geomorphic and sedimentological analyses confirm that the moraines are predominately composed of rockfall and avalanche debris mobilized from steep relief rockwall slopes via periglacial weathering processes. The glacial rockwall slope erosion affects sediment flux and storage of snow and ice at the catchment head on diurnal to millennial timescales, and more broadly influences catchment configuration and relief, glacier dynamics and microclimates. The slope erosion rates exceed the averaged catchment-wide and exhumation rates of Bhagirathi and the Garhwal region on geomorphic timescales (10³−10⁵ years), supporting the view that erosion at the headwaters can outpace the wider catchment. The ¹⁰Be concentrations of medial moraine sediment for the upper Bhagirathi catchment and the catchments of Chhota Shigri in Lahul, northern India and Baltoro glacier in Central Karakoram, Pakistan show a tentative relationship between ¹⁰Be concentration and precipitation. As such there is more rapid glacial rockwall slope erosion in the monsoon-influenced Lesser and Greater Himalaya compared to the semi-arid interior of the orogen. Rockwall slope erosion in the three study areas, and more broadly across the northwest Himalaya is likely governed by individual catchment dynamics that vary across space and time. © 2019 The Authors. Earth Surface Processes and Landforms Published by John Wiley & Sons, Ltd.     

Runoff generation and routing on artificial slopes in a Mediterranean–continental environment: the Teruel coalfield,Spain

The aim of this study was to identify the mechanisms of runoff generation and routing and their controlling factors at the hillslope scale, on artificial slopes derived from surface coal mining reclamation in a Mediterranean–continental area. Rainfall and runoff at interrill and microcatchment scales were recorded for a year on two slopes with different substrata: topsoil cover and overburden cover. Runoff coefficient and runoff routing from interrill areas to microcatchment outlets were higher in the overburden substratum than in topsoil, and greater in the most developed rill network. Rainfall volume is the major parameter responsible for runoff response on overburden, suggesting that this substratum is very impermeable—at least during the main rainfall periods of the year (late spring and autumn) when the soil surface is sealed. In such conditions, most rainfall input is converted into runoff, regardless of its intensity. Results from artificial rainfall experiments, conducted 3 and 7 years after seeding, confirm the low infiltration capacity of overburden when sealed. The hydrological response shows great seasonal variability on the overburden slope in accordance with soil surface changes over the year. Rainfall volume and intensities (I₃₀, I₆₀) explain runoff at the interrill scale on the topsoil slope, where rainfall experiments demonstrated a typical Hortonian infiltration curve. However, no correlation was found at the microcatchment level, probably because of the loss of functionality of the only rill as ecological succession proceeded. The runoff generation mechanism on the topsoil slope is more homogeneous throughout the year. Runoff connectivity, defined as the ratio between runoff rates recorded at the rill network scale and those recorded at the interrill area scale in every rainfall event, was also greater on the rilled overburden slope, and in the most developed rill network. The dense rill networks of the overburden slope guarantee very effective runoff drainage, regardless of rainfall magnitude. Rills drain overland flow from interrill‐sealed areas, reducing the opportunity of reinfiltration in areas not affected by siltation. Runoff generation and routing on topsoil slopes are controlled by grass cover and soil moisture content, whereas on overburden slopes rill network density and soil moisture content are the main controlling factors. Copyright © 2002 John Wiley & Sons, Ltd.