Impacts of post‐glacial rebound on landslide spatial distribution at a regional scale in northern Iceland (Skagafjörður)

In de‐glaciated areas, para‐glaciation (i.e. the conditioning of landscapes by prior glaciation) has often been considered a major predisposing factor in landslide occurrence; its consequences have been particularly well identified at a fine scale (especially on bedrock jointing). Hitherto, the relative impacts of para‐glaciation on hillslope dynamics at a regional scale had nevertheless not been quantified statistically. We examine Skagafjörður area (northern Iceland) where landslides are widespread (at least 108 were mapped in an area of c. 3000 km²). We compare the role of para‐glaciation (debuttressing, influence of post‐glacial rebound) with that of classic factors (topography, lithology, etc.) in landslide occurrence and location, using a spatial analysis based on a chi‐square test. On the one hand, the results highlight that landslides are over‐represented in areas where post‐glacial rebound was at its maximum, with a stronger concentration of landslides in the northern part of the fjord. On the other hand, the distribution of landslides did not show any clear relationship with the pattern of glacial debuttressing. Tschuprow coefficient highlights that the influence of post‐glacial rebound on landslide location is higher than the combined influence of slope gradient, curvature or geological structure. This result is supported by our initial evidence for the timing of landslides in the area: most landslides occurred during the first half of the Holocene, and a period of hillslope instability was initiated when the post‐glacial uplift was at its maximum. Finally, the mechanisms that link post‐glacial rebound and landsliding as well as the geomorphic impacts of landslides, are discussed. Copyright © 2013 John Wiley & Sons, Ltd.     

Modelling catchment‐scale shallow landslide occurrence and sediment yield as a function of rainfall return period

A model‐based method is proposed for improving upon existing threshold relationships which define the rainfall conditions for triggering shallow landslides but do not allow the magnitude of landsliding (i.e. the number of landslides) to be determined. The SHETRAN catchment‐scale shallow landslide model is used to quantify the magnitude of landsliding as a function of rainfall return period, for focus sites of 180 and 45 km² in the Italian Southern Alps and the central Spanish Pyrenees. Rainfall events with intensities of different return period are generated for a range of durations (1‐day to 5‐day) and applied to the model to give the number of landslides triggered and the resulting sediment yield for each event. For a given event duration, simulated numbers of landslides become progressively less sensitive to return period as return period increases. Similarly, for an event of given return period, landslide magnitude becomes less sensitive to event duration as duration increases. The temporal distribution of rainfall within an event is shown to have a significant impact on the number of landslides and the timing of their occurrence. The contribution of shallow landsliding to catchment sediment yield is similarly quantified as a function of the rainfall characteristics. Rainfall intensity–duration curves are presented which define different levels of landsliding magnitude and which advance our predictive capability beyond, but are generally consistent with, published threshold curves. The magnitude curves are relevant to the development of guidelines for landslide hazard assessment and forecasting. Copyright © 2011 John Wiley & Sons, Ltd.     

Supra‐glacial deposition and flux of catastrophic rock–slope failure debris,south‐central Alaska

The ongoing debate over the effects of global environmental change on Earth's cryosphere calls for detailed knowledge about process rates and their variability in cold environments. In this context, appraisals of the coupling between glacier dynamics and para‐glacial erosion rates in tectonically active mountains remain rare. We contribute to filling this knowledge gap and present an unprecedented regional‐scale inventory of supra‐glacial sediment flux and hillslope erosion rates inferred from an analysis of 123 large (> 0·1 km²) catastrophic bedrock landslides that fell onto glaciers in the Chugach Mountains, Alaska, as documented by satellite images obtained between 1972 to 2008. Assuming these supra‐glacial landslide deposits to be passive strain markers we infer minimum decadal‐scale sediment yields of 190 to 7400 t km^–2 yr^–1 for a given glacier‐surface cross‐section impacted by episodic rock–slope failure. These rates compare to reported fluvial sediment yields in many mountain rivers, but are an order of magnitude below the extreme sediment yields measured at the snouts of Alaskan glaciers, indicating that the bulk of debris discharged derives from en‐glacial, sub‐glacial or ice‐proximal sources. We estimate an average minimum para‐glacial erosion rate by large, episodic rock–slope failures at 0·5–0·7 mm yr^–1 in the Chugach Mountains over a 50‐yr period, with earthquakes likely being responsible for up to 73% of this rate. Though ranking amongst the highest decadal landslide erosion rates for this size of study area worldwide, our inferred rates of hillslope erosion in the Chugach Mountains remain an order of magnitude below the pace of extremely rapid glacial sediment export and glacio‐isostatic surface uplift previously reported from the region. Copyright © 2012 John Wiley & Sons, Ltd.     

Late-stage volcano geomorphic evolution of the Pleistocene San Francisco Mountain, Arizona (USA), based on high-resolution DEM analysis and 40Ar/39Ar chronology

The cone-building volcanic activity and subsequent erosion of San Francisco Mountain, AZ, USA, were studied by using high-resolution digital elevation model (DEM) analysis and new ⁴⁰Ar/³⁹Ar dating. By defining remnants or planèzes of the volcano flanks in DEM-derived images, the original edifice can be reconstructed. We propose a two-cone model with adjacent summit vents which were active in different times. The reconstructed cones were 4,460 and 4,350 m high a.s.l., corresponding to ∼2,160 and 2,050 m relative height, respectively. New ⁴⁰Ar/³⁹Ar data allow us to decipher the chronological details of the cone-building activity. We dated the Older and Younger Andesites of the volcano that, according to previous mapping, built the stage 2 and stage 3 stratocones, respectively. The new ⁴⁰Ar/³⁹Ar plateau ages yielded 589–556 ka for the Older and 514–505 ka for the Younger Andesites, supporting their distinct nature with a possible dormant period between. The obtained ages imply an intense final (≤100 ka long) cone-building activity, terminating ∼100 ka earlier than indicated by previous K-Ar ages. Moreover, ⁴⁰Ar/³⁹Ar dating constrains the formation of the Inner Basin, an elliptical depression in the center of the volcano initially created by flank collapse. A 530 ka age (with a ±58.4 ka 2σ error) for a post-depression dacite suggests that the collapse event is geochronologically indistinguishable from the termination of the andesitic cone-building activity. According to our DEM analysis, the original cone of San Francisco Mountain had a volume of about 80 km³. Of this volume, ∼7.5 km³ was removed by the flank collapse and subsequent glacial erosion, creating the present-day enlarged Inner Basin, and ∼2 km³ was removed from the outer valleys by erosion. Based on volumetric analysis and previous and new radiometric ages, the average long-term eruption rate of San Francisco Mountain was ∼0.2 km³/ka, which is a medium rate for long-lived stratovolcanoes. However, according to the new ⁴⁰Ar/³⁹Ar dates for the last ≤100 ka period, the final stratovolcanic activity was characterized by a greater ∼0.3 km³/ka rate.     

OSL dating of past lake levels for a large dammed lake in southern Tibet and determination of possible controls on lake evolution

Fluvio‐lacustrine terraces along Phung Chu (river) on the central southern Tibetan Plateau indicate that a large palaeo‐dammed‐lake formerly existed in this area. Based on landscape survey, optically stimulated luminescence (OSL) dating and sedimentary analyses, this research shows that the Phung Chu was blocked and a dammed‐lake over 2500 km² in size formed before 30 ka ago. OSL dating analysis suggests the fluvio‐lacustrine sediments were well bleached and yield accurate age estimates for two lake drainage events. The first drainage event took place after 30 ka, resulted in river incision and formed a high terrace at 50 m height from the present river level. The second drainage happened after 3.7 ka, resulted in further river incision and formed the second terrace at 25 m height from the present river level. According to the distribution of the fluvio‐lacustrine sediments, active normal faults (particularly the Kharta Fault) in this region and the high gradient slopes after Phung Chu enters the Yö Ri gorge, seismically‐induced landsliding is regarded as highly likely to have been the cause of river blockage and associated formation of a dammed‐lake, although glacial damming is also a possible cause. The volume of drainages from this dammed‐lake may have led to catastrophic flooding and analogous modern lakes represent significant geo‐hazard risks to down‐river human settlements. As dammed‐lakes are special phases in fluvial evolution, often involving river blockage, breakthrough and drastic catchment change, these processes can reveal how tectonic or climatic events modify landforms. However, such tectonic‐derived landform changes can also impact palaeo‐climate of the region. Thus this study has added new evidence regarding the evolutionary history of a dammed lake including its formation, duration, extent and final drainage, which is crucial for understanding its general landscape process mechanisms and for better assessing geo‐hazard risks. Copyright © 2016 John Wiley & Sons, Ltd.     

Luminescence chronology from the Verkhoyansk Mountains (North-Eastern Siberia)

This study presents new dating results from the Verkhoyansk Mountains in northeastern Siberia. Pleistocene sediments of aeolian and glacial origin have been studied and dated by infrared optically stimulated luminescence (IRSL). The chronostratigraphy of this vast area is mainly based on radiocarbon dating up to now. Aeolian sediments are widespread in the foreland of the Verkhoyansk Mountains. IRSL-dating results indicate two major periods of increased accumulation between 33 and 24 ka and between 13 and 9 ka. A new stratigraphy of Pleistocene glaciations in this area has been developed. Up to five end moraines have been identified in two catchments areas. IRSL was applied to date the sandy and silty sediments covering the glacial deposits in order to get minimum ages. Furthermore, glacial and glacio-fluvial deposits were also sampled from a few locations. According to the IRSL dating results, the uppermost end moraine was deposited prior to 50 ka. The three outermost moraines might have been formed during Early Weicheslian to Saalian times. It is very likely that no major glaciers reached the foreland of the Verkhoyansk Mountains during the Last Glacial Maximum (LGM, MIS 2).     

Enhanced rock‐slope failure following ice‐sheet deglaciation: timing and causes

The temporal pattern of rock‐slope failures (RSFs) following Late Pleistocene deglaciation on tectonically stable terrains is controversial: previous studies variously suggest (1) a rapid response due to removal of supporting ice (‘debuttressing’), (2) a progressive decline in RSF frequency, and (3) a millennial‐scale delay before peak RSF activity. We test these competing models through beryllium‐10 (¹⁰Be) exposure dating of five closely‐spaced quartzite RSFs on the Isle of Jura, Scotland, to establish the relationship between timing of failure and those of deglaciation, episodes of rapid warming and periods of rapid glacio‐isostatic uplift. All five dated RSFs occurred at least 720–2240 years after deglaciation, with the probability of failure peaking ~2 ka after deglaciation, consistent with millennial‐scale delay model (3). This excludes debuttressing as an immediate cause of failure, though it is likely that time‐dependent stress release due to deglacial unloading resulted in progressive development of failure planes within the rock. Thaw of permafrost ice in joints is unlikely to have been a prime trigger of failure as some RSFs occurred several centuries after the onset of interstadial warming. Conversely, the timespan of the RSFs coincides with the period of maximum glacio‐isostatic crustal uplift, suggesting that failure was triggered by uplift‐driven seismic events acting on fractured rock masses. Implications of this and related research are: (1) that retreat of the last Pleistocene ice sheets across tectonically‐stable mountainous terrains was succeeded by a period of enhanced rock‐slope failure due to deglacial unloading and probably uplift‐driven seismicity; (2) that the great majority of RSFs in the British Isles outside the limits of Loch Lomond Stadial (= Younger Dryas) glaciation are of Lateglacial (pre‐Holocene) age; and (3) numerous RSFs must also have occurred inside Loch Lomond Stadial (LLS) glacial limits, but that runout debris was removed by LLS glaciers. Copyright © 2013 John Wiley & Sons, Ltd.     

Large landslides and their effect on sediment flux in South Westland,New Zealand

Landslides and runoff are dominant erosional agents in the tectonically active alpine South Westland area of New Zealand, characterized by high uplift rates and extreme orographic precipitation. Despite a high density of shallow debris slides and flows, the geomorphic imprints of deep‐seated bedrock failures are dominant and persistent. Over 50 large (>1 km²) landslides comprising rock slide[sol ]avalanches, complex rotational and rock‐block slides, wedge failures, and deep‐seated gravitational slope deformation were detected on air photos and shaded‐relief images. Major long‐term impacts on alpine rivers include (1) forced alluviation upstream of landslide dams, (2) occlusion of gorges and triggering of secondary riparian landslides, and (3) diversion of channels around deposits to form incised meandering gorges. Remnants of large prehistoric (i.e. pre‐1840) landslide deposits possibly represent the low‐frequency (in terms of total area affected yet dominant) end of the spectrum of mass wasting in the western Southern Alps. This is at odds with high erosion rates in an active erosional landscape. Large landslides appear to have dual roles of supplying and retaining sediment. The implications of these roles are that (1) previous models of (shallow) landslide‐derived sediment flux need to be recalibrated, and (2) geomorphic effects of earthquake‐induced landsliding may persist for at least 10² years. Copyright © 2005 John Wiley & Sons, Ltd.     

Postglacial sediment budget of Chilliwack Valley,British Columbia

The paraglacial reworking of glacial sediments by rivers and mass wasting is an important conditioning factor for modern sediment yields in mountainous catchments in formerly glaciated regions. Catchment scale and patterns of sediment storage are important influences in the rate of postglacial adjustment. We develop a quantitative framework to estimate the volume, sediment type, and fractional size distribution of legacy glacial materials in a large (1230 km²) watershed in the North Cascade Mountains in south‐western British Columbia, Canada. Chilliwack Valley is exceptional because of the well‐dated bounds of deglaciation. Interpolation of paleo‐surfaces from partially eroded deposits in the valley allows us to estimate the total evacuated sediment volume. We present a chronology of sediment evacuation from the valley and deposition in the outlet fan, based on infrared stimulated luminescence (IRSL) and ¹⁴ C dating of river terraces and fan strata, respectively. The effects of paraglacial sedimentation in Chilliwack Valley were intensified through a major fall in valley base‐level following ice retreat. The steepened mainstem valley gradient led to deep incision of valley fills and fan deposits in the lower valley network. The results of this integrated study provide a postglacial chronology and detailed sediment budget, accounting for long‐term sorting of the original sediments, lag deposit formation in the mainstem, deposition in the outlet fan, and approximate downstream losses of suspended sediment and wash load. The mass balance indicates that a bulk volume of approximately 3.2 km³ of glacial material has been evacuated from the valley. Copyright © 2012 John Wiley & Sons, Ltd.     

Application of the SHETRAN basin‐scale,landslide sediment yield model to the Llobregat basin,Spanish Pyrenees

The first application of the SHETRAN basin‐scale, landslide erosion and sediment yield model is carried out for a major landsliding event in the upper 505 km² of the Llobregat basin, in the eastern Spanish Pyrenees, in November 1982. The model simulates the spatial distribution of shallow landslides and their sediment yield. Acknowledging uncertainty in the model parameter evaluation, the aim of the application was not to reproduce the observed occurrence of landslides as accurately as possible with one simulation, but to bracket the observed pattern with several simulations representing uncertainty in the key input conditions. Bounds on the landslide simulations were thus determined as a function of uncertainty in the vegetation root cohesion (used in the model factor of safety calculations). The resulting upper bound considerably overestimates the observed pattern (17 000 landslides compared with an observation of around 700), but it reproduces several of the principal clusters in the observed pattern. The lower bound contains around 500 landslides. The sediment yield estimates (2670–14 630 t km^?2) are comparable to measurements elsewhere in the Pyrenees for extreme events. The results demonstrate an ability to simulate the basin‐scale landslide response to a rainfall event and the resulting sediment yield. They also highlight the need for further research in setting the uncertainty bounds and in avoiding large overestimates of landslide occurrence arising in part from a current inability to model small‐scale controls for a basin of the given size. Copyright © 2006 John Wiley & Sons, Ltd.     

Postglacial topographic evolution of glaciated valleys: a stochastic landscape evolution model

The retreat of valley glaciers has a dramatic effect on the stability of glaciated valleys and exerts a prolonged influence on the subsequent fluvial sediment transport regime. We have studied the evolution of an idealized glaciated valley during the period following retreat of ice using a numerical model. The model incorporates a stochastic process to represent deep‐seated landsliding, non‐linear diffusion to represent shallow landsliding and an approximation of the Bagnold relation to represent fluvial sediment transport. It was calibrated using field data from several recent surveys within British Columbia, Canada. We present ensemble model results and compare them with results from a deterministic linear‐diffusion model to show that explicit representation of large landslides is necessary to reproduce the morphology and channel network structure of a typical postglacial valley. Our model predicts a rapid rate of fluvial sediment transport following deglaciation with a subsequent gradual decline, similar to that inferred for Holocene time. We also describe how changes in the model parameters affect the estimated magnitude and duration of the paraglacial sediment pulse. Copyright © 2005 John Wiley & Sons, Ltd.     

Age of glacial and periglacial landforms in northern Hokkaido, Japan, using OSL dating of fine grain quartz

The timing of glacial advances, periglacial phenomena, and the ages of two marker tephras in northern Hokkaido were estimated by OSL dating. It appears that the glacier of Yamunai 2 stage on Rishiri Island expanded between 24 and 15 ka. In northern Hokkaido, OSL ages indicate ice wedge formation during the period 24–18 ka. These results indicate that both the glacial advance and the development of ice wedges were synchronous phenomena relating to the Last Glacial Maximum.     

Estimation of regional material yield from coastal landslides based on historical digital terrain modelling

High‐resolution historical (1942) and recent (1994) digital terrain models were derived from aerial photographs along the Big Sur coastline in central California to measure the long‐term volume of material that enters the nearshore environment. During the 52‐year measurement time period, an average of 21 000 ± 3100 m³ km^?1 a^?1 of material was eroded from nine study sections distributed along the coast, with a low yield of 1000 ± 240 m³ km^?1 a^?1 and a high of 46 700 ± 7300 m³ km^?1 a^?1. The results compare well with known volumes from several deep‐seated landslides in the area and suggest that the processes by which material is delivered to the coast are episodic in nature. In addition, a number of parameters are investigated to determine what influences the substantial variation in yield along the coast. It is found that the magnitude of regional coastal landslide sediment yield is primarily related to the physical strength of the slope‐forming material. Coastal Highway 1 runs along the lower portion of the slope along this stretch of coastline, and winter storms frequently damage the highway. The California Department of Transportation is responsible for maintaining this scenic highway while minimizing the impacts to the coastal ecosystems that are part of the Monterey Bay National Marine Sanctuary. This study provides environmental managers with critical background data on the volumes of material that historically enter the nearshore from landslides, as well as demonstrating the application of deriving historical digital terrain data to model landscape evolution. Published in 2005 by John Wiley & Sons, Ltd.     

10Be‐derived assessment of accelerated erosion in a glacially conditioned inner gorge,Entlebuch, Central Alps of Switzerland

Inner gorges often result from the propagation of erosional waves related to glacial/interglacial climate shifts. However, only few studies have quantified the modern erosional response to this glacial conditioning. Here, we report in situ ¹⁰Be data from the 64 km² Entlen catchment (Swiss Alps). This basin hosts a 7 km long central inner gorge with two tributaries that are >100 m‐deeply incised into thick glacial till and bedrock. The ¹⁰Be concentrations measured at the downstream end of the gorge yield a catchment‐wide erosion rate of 0.42 ± 0.04 mm yr^‐1, while erosion rates are consistently lower upstream of the inner gorge, ranging from 0.14 ± 0.01 mm yr^‐1 to 0.23 ± 0.02 mm yr^‐1. However, ¹⁰Be‐based sediment budget calculations yield rates of ~1.3 mm yr^‐1 for the inner gorge of the trunk stream. Likewise, in the two incised tributary reaches, erosion rates are ~2.0 mm yr^‐1 and ~1.9 mm yr^‐1. Moreover, at the erosional front of the gorge, we measured bedrock incision rates ranging from ~2.5 mm yr^‐1 to ~3.8 mm yr^‐1. These rates, however, are too low to infer a post‐glacial age (15–20 ka) for the gorge initiation. This would require erosion rates that are between 2 and 6 times higher than present‐day estimates. However, the downcutting into unconsolidated glacial till favored high erosion rates through knickzone propagation immediately after the retreat of the LGM glaciers, and subsequent hillslope relaxation led to a progressive decrease in erosion rates. This hypothesis of a two‐ to sixfold decrease in erosion rates does not conflict with the ¹⁰Be‐based erosion rate budgets, because the modern erosional time scale recorded by ¹⁰Be cover the past 2–3 ka only. These results point to the acceleration of Holocene erosion in response to the glacial overprint of the landscape. Copyright © 2012 John Wiley & Sons, Ltd.     

Ice cap erosion patterns from bedrock 10Be and 26Al,southeastern Tibetan Plateau

Quantifying glacial erosion contributes to our understanding of landscape evolution and topographic relief production in high altitude and high latitude areas. Combining in situ ¹⁰Be and ²⁶Al analysis of bedrock, boulder, and river sand samples, geomorphological mapping, and field investigations, we examine glacial erosion patterns of former ice caps in the Shaluli Shan of the southeastern Tibetan Plateau. The general landform pattern shows a zonal pattern of landscape modification produced by ice caps of up to 4000 km² during pre-LGM (Last Glacial Maximum) glaciations, while the dating results and landforms on the plateau surface imply that the LGM ice cap further modified the scoured terrain into different zones. Modeled glacial erosion depth of 0–0.38 m per 100 ka bedrock sample located close to the western margin of the LGM ice cap, indicates limited erosion prior to LGM and Late Glacial moraine deposition. A strong erosion zone exists proximal to the LGM ice cap marginal zone, indicated by modeled glacial erosion depth >2.23 m per 100 ka from bedrock samples. Modeled glacial erosion depths of 0–1.77 m per 100 ka from samples collected along the edge of a central upland, confirm the presence of a zone of intermediate erosion in-between the central upland and the strong erosion zone. Significant nuclide inheritance in river sand samples from basins on the scoured plateau surface also indicate restricted glacial erosion during the last glaciation. Our study, for the first time, shows clear evidence for preservation of glacial landforms formed during previous glaciations under non-erosive ice on the Tibetan Plateau. As patterns of glacial erosion intensity are largely driven by the basal thermal regime, our results confirm earlier inferences from geomorphology for a concentric basal thermal pattern for the Haizishan ice cap during the LGM. © 2018 John Wiley & Sons, Ltd.     

Use of LIDAR‐derived images for mapping old landslides under forest

Large, deep‐seated landslides are common features in the Flemish Ardennes (Belgium). As most of these old (>100 years) landslides are located under forest in this hilly region, aerial photograph interpretation is not an appropriate landslide mapping method. This study tested the potential of LIDAR (Light Detection and Ranging) images for mapping old landslides under forest. Landslide inventory maps were created for a 125 km² area by applying the expert knowledge of seven geomorphologists to LIDAR‐derived hillshade, slope and contour line maps in a GIS environment. Each of the seven LIDAR‐based landslide inventories was compared (i) with the other six, (ii) with a detailed field survey‐based inventory, and (iii) with a comparable study in which topographic data were extracted from a topographical map. The combination of the percentage of field landslides indicated by an expert and the percentage of positional discrepancies (expressed in terms of positional mismatch) were used to evaluate the quality of the LIDAR‐based inventory maps. High‐quality LIDAR‐derived landslide inventory maps contain more than 70 per cent of the landslides mapped during the field survey, and have positional discrepancies smaller than 70 per cent when compared with the field survey‐based inventory map. Four experts and the combination map of all experts satisfied these criteria. Together the seven experts indicated all landslides mapped in the field. Importantly, LIDAR enabled the experts to find ten new landslides and to correct the boundaries of eleven (of the 77) landslides mapped during the field survey. Hence, this study showed that large‐scale LIDAR‐derived maps analysed by experienced geomorphologists can significantly improve field survey‐based inventories of landslides with a subdued morphology in hilly regions. Copyright © 2006 John Wiley & Sons, Ltd.     

Cosmogenic 36Cl glacial chronologies of the Late Quaternary glaciers on Mount Geyikdağ in the Eastern Mediterranean

We report the timing of glaciations during the Late Quaternary in the central Taurus Mountains of Turkey in the Eastern Mediterranean. Forty moraine samples from three glacial valleys on Mount Geyikdağ (36.53°N, 32.10^°E, 2877 m), near the Eastern Mediterranean coast of Turkey, were dated with in-situ cosmogenic ³⁶Cl. These glacial valleys are located on the southern flank of the mountain and were filled with few km long glaciers that terminated at elevations of about 1750 m above sea level. Three glacial retreats/advances were determined in this study. During the Last Glacial Maximum (LGM), glaciers reached their maximum positions at 20.6 ± 0.6 ka ago (±1σ). This date is in accordance with the timing of local maximum ice extent, represented by piedmont glaciers in the northern side of the mountain. Glaciers started to retreat after the LGM and shortly stabilized or re-advanced two times before they completely vanished out. The first stage ended before 13.7 ± 0.8 ka ago during the Late-glacial. The last glaciation occurred during the Holocene and ended between 9.6 ± 1.4 ka and 5.9 ± 0.5 ka ago. Later, glaciers mostly vanished from the study area, but a few rock glaciers developed during the Late Holocene. Glacial chronologies of Mount Geyikdağ are mostly comparable with the globally observed advances elsewhere.     

Digital landscapes of deglaciation: identifying Late Quaternary glacial lake outburst floods using LiDAR

High resolution DEMs obtained from LiDAR topographic data have led to improved landform inventories (e.g. landslides and fault scarps) and understanding of geomorphic event frequency. Here we use airborne LiDAR mapping to investigate meltwater pathways associated with the Tweed Valley palaeo ice‐stream (UK). In particular we focus on a gorge downstream of Palaeolake Milfield, previously mapped as a sub‐glacial meltwater channel, where the identification of abandoned headcut channels, run‐up bars, rock‐cut terrace surfaces and eddy flow features attest to formation by a sub‐aerial glacial lake outburst flood (GLOF) caused by breaching of a sediment dam, likely an esker ridge. Mapping of these landforms combined with analysis of the gorge rim elevations and cross‐section variability revealed a two phase event with another breach site downstream following flow blockage by higher elevation drumlin topography. We estimate the magnitude of peak flow to be 1–3 × 10³ m³/s, duration of the event to range from 16–155 days, and a specific sediment yield of 10⁷–10⁹ m³/km²/yr. We identified other outburst pathways in the lower Tweed basin that help delineate an ice margin position of the retreating Tweed Valley ice stream. The results suggest that low magnitude outburst floods are under‐represented in Quaternary geomorphological maps. We therefore recommend regional LiDAR mapping of meltwater pathways to identify other GLOFs in order to better quantify the pattern of freshwater and sediment fluxes from melting ice sheets to oceans. Despite the relatively low magnitude of the Till outburst event, it had a significant impact on the landscape development of the lower Tweed Valley through the creation of a new tributary pathway and triggering of rapid knickpoint retreat encouraging new regional models of post‐glacial fluvial landscape response. Copyright © 2015 John Wiley & Sons, Ltd.     

10Be surface-exposure age dating of the Last Glacial Maximum in the northern Pamir (Tajikistan)

Knowledge of the spatial and temporal variations in Alpine glaciations is essential for reconstructing the regional and global timing of ice ages. This study investigates glacial deposits at the mouth of the Muksu catchment in the northern Pamir using ¹⁰Be surface-exposure age dating. We sampled boulders from the furthest downstream recessional moraine (20 samples) and five lateral moraines (41 samples) near the former terminus of the Fedchenko Glacier, the longest (∼72 km) present-day Alpine glacier of the Pamir. After the identification of outliers, the boulder population of the recessional moraine yielded a mean exposure age of 17.5 ± 1.9 ka. The maximum exposure age of the lateral moraines, collected ∼5 km up-valley of the recessional moraine, is 18.2 ± 1.7 ka. The boulder ages reflect glacial deposition during the Last Glacial Maximum (Marine Isotope Stage 2) in the region; they are in accordance with published glacial deposition ages in the western Tian Shan.