Refine
Year of publication
Document Type
Is part of the Bibliography
- yes (80)
Keywords
- Himalayas (5)
- Landslide (5)
- landslide (5)
- Earthquake (4)
- Erosion (4)
- hydropower (4)
- Germany (3)
- Himalaya (3)
- Patagonia (3)
- Tibetan Plateau (3)
Glacial advances constrained by Be-10 exposure dating of bedrock landslides, Kyrgyz Tien Shan
(2011)
Numerous large landslide deposits occur in the Tien Shan, a tectonically active intraplate orogen in Central Asia. Yet their significance in Quaternary landscape evolution and natural hazard assessment remains unresolved due to the lack of "absolute" age constraints. Here we present the first Be-10 exposure ages for three prominent (>10(7) m(3)) bedrock landslides that blocked major rivers and formed lakes, two of which subsequently breached, in the northern Kyrgyz Tien Shan. Three Be-10 ages reveal that one landslide in the Alamyedin River occurred at 11-15 ka, which is consistent with two C-14 ages of gastropod shells from reworked loess capping the landslide. One large landslide in Aksu River is among the oldest documented in semi-arid continental interiors, with a Be-10 age of 63-67 ka. The Ukok River landslide deposit(s) yielded variable Be-10 ages, which may result from multiple landslides, and inheritance of Be-10. Two Be-10 ages of 8.2 and 5.9 ka suggest that one major landslide occurred in the early to mid-Holocene, followed by at least one other event between 1.5 and 0.4 ka. Judging from the regional glacial chronology, all three landslides have occurred between major regional glacial advances. Whereas Alamyedin and Ukok can be considered as postglacial in this context, Aksu is of interglacial age. None of the landslide deposits show traces of glacial erosion, hence their locations and I Be ages mark maximum extents and minimum ages of glacial advances, respectively. Using toe-to-headwall altitude ratios of 0.4-0.5, we reconstruct minimum equilibrium-line altitudes that exceed previous estimates by as much as 400 m along the moister northern fringe of the Tien Shan. Our data show that deposits from large landslides can provide valuable spatio-temporal constraints for glacial advances in landscapes where moraines and glacial deposits have low preservation potential. (C) 2011 University of Washington. Published by Elsevier Inc. All rights reserved.
Earthquake-triggered landslide dams are potentially dangerous disrupters of water and sediment flux in mountain rivers, and capable of releasing catastrophic outburst flows to downstream areas. We analyze an inventory of 828 landslide dams in the Longmen Shan mountains, China, triggered by the M-w 7.9 2008 Wenchuan earthquake. This database is unique in that it is the largest of its kind attributable to a single regional-scale triggering event: 501 of the spatially clustered landslides fully blocked rivers, while the remainder only partially obstructed or diverted channels in steep watersheds of the hanging wall of the Yingxiu-Beichuan Fault Zone. The size distributions of the earthquake-triggered landslides, landslide dams, and associated lakes (a) can be modeled by an inverse gamma distribution; (b) show that moderate-size slope failures caused the majority of blockages; and (c) allow a detailed assessment of seismically induced river-blockage effects on regional water and sediment storage. Monte Carlo simulations based on volumetric scaling relationships for soil and bedrock failures respectively indicate that 14% (18%) of the estimated total coseismic landslide volume of 6.4 (14.6) x 10(9) m(3) was contained in landslide dams, representing only 1.4% of the >60,000 slope failures attributed to the earthquake. These dams have created storage capacity of similar to 0.6x 10(9) m(3) for incoming water and sediment. About 25% of the dams containing 2% of the total river-blocking debris volume failed one week after the earthquake; these figures had risen to 60% (similar to 20%), and >90% (>90%) within one month, and one:year, respectively, thus also emptying similar to 92% of the total potential water and sediment storage behind these, dams within one year following the earthquake. Currently only similar to 0.08 x 10(9) m(3) remain available as natural reservoirs for storing water and sediment, while similar to 0.19 x 10(9) m(3), i.e. about a third of the total river-blocking debris volume, has been eroded by rivers. Dam volume and upstream catchment area control to first order the longevity of the barriers, and bivariate domain plots are consistent with the observation that most earthquake-triggered landslide dams were ephemeral. We conclude that the river-blocking portion of coseismic slope failures disproportionately modulates the post-seismic sediment flux in the Longmen Shan on annual to decadal timescales.
Quantitative estimates of sediment flux and the global cycling of sediments from hillslopes to rivers, estuaries, deltas, continental shelves, and deep-sea basins have a long research tradition. In this context, extremely large and commensurately rare sediment transport events have so far eluded a systematic analysis. To start filling this knowledge gap I review some of the highest reported sediment yields in mountain rivers impacted by volcanic eruptions, earthquake- and storm-triggered landslide episodes, and catastrophic dam breaks. Extreme specific yields, defined here as those exceeding the 95th percentile of compiled data, are similar to 10(4) t km(-2) yr(-1) if averaged over 1 yr. These extreme yields vary by eight orders of magnitude, but systematically decay with reference intervals from minutes to millennia such that yields vary by three orders of magnitude for a given reference interval. Sediment delivery from natural dam breaks and pyroclastic eruptions dominate these yields for a given reference interval. Even if averaged over 10(2)-10(3) yr, the contribution of individual disturbances may remain elevated above corresponding catchment denudation rates. I further estimate rates of sediment (re-)mobilisation by individual giant terrestrial and submarine mass movements. Less than 50 postglacial submarine mass movements have involved an equivalent of similar to 10% of the contemporary annual global flux of fluvial sediment to Earth's oceans, while mobilisation rates by individual events rival the decadal-scale sediment discharge from tectonically active orogens such as Taiwan or New Zealand. Sediment flushing associated with catastrophic natural dam breaks is non-stationary and shows a distinct kink at the last glacial-interglacial transition, owing to the drainage of very large late Pleistocene ice-marginal lakes. Besides emphasising the contribution of high-magnitude and low-frequency events to the global sediment cascade, these findings stress the importance of sediment storage for fuelling rather than buffering high sediment transport rates.
Climate change, manifested by an increase in mean, minimum, and maximum temperatures and by more intense rainstorms, is becoming more evident in many regions. An important consequence of these changes may be an increase in landslides in high mountains. More research, however, is necessary to detect changes in landslide magnitude and frequency related to contemporary climate, particularly in alpine regions hosting glaciers, permafrost, and snow. These regions not only are sensitive to changes in both temperature and precipitation, but are also areas in which landslides are ubiquitous even under a stable climate. We analyze a series of catastrophic slope failures that occurred in the mountains of Europe, the Americas, and the Caucasus since the end of the 1990s. We distinguish between rock and ice avalanches, debris flows from de-glaciated areas, and landslides that involve dynamic interactions with glacial and river processes. Analysis of these events indicates several important controls on slope stability in high mountains, including: the non-linear response of firn and ice to warming; three-dimensional warming of subsurface bedrock and its relation to site geology; de-glaciation accompanied by exposure of new sediment; and combined short-term effects of precipitation and temperature. Based on several case studies, we propose that the following mechanisms can significantly alter landslide magnitude and frequency, and thus hazard, under warming conditions: (1) positive feedbacks acting on mass movement processes that after an initial climatic stimulus may evolve independently of climate change; (2) threshold behavior and tipping points in geomorphic systems; (3) storage of sediment and ice involving important lag-time effects.
Projected scenarios of climate change involve general predictions about the likely changes to the magnitude and frequency of landslides, particularly as a consequence of altered precipitation and temperature regimes. Whether such landslide response to contemporary or past climate change may be captured in differing scaling statistics of landslide size distributions and the erosion rates derived thereof remains debated. We test this notion with simple Monte Carlo and bootstrap simulations of statistical models commonly used to characterize empirical landslide size distributions. Our results show that significant changes to total volumes contained in such inventories may be masked by statistically indistinguishable scaling parameters, critically depending on, among others, the size of the largest of landslides recorded. Conversely, comparable model parameter values may obscure significant, i.e. more than twofold, changes to landslide occurrence, and thus inferred rates of hillslope denudation and sediment delivery to drainage networks. A time series of some of Earth's largest mass movements reveals clustering near and partly before the last glacial-interglacial transition and a distinct step-over from white noise to temporal clustering around this period. However, elucidating whether this is a distinct signal of first-order climate-change impact on slope stability or simply coincides with a transition from short-term statistical noise to long-term steady-state conditions remains an important research challenge.
Supra-glacial deposition and flux of catastrophic rock-slope failure debris, south-central Alaska
(2013)
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 center dot 1km2) 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 7400tkm-2yr-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 center dot 5-0 center dot 7mmyr-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.
The 12 January 2010 M-w 7.0 Haiti earthquake occurred in a complex deformation zone at the boundary between the North American and Caribbean plates. Combined geodetic, geological and seismological data posited that surface deformation was driven by rupture on the Leogane blind thrust fault, while part of the rupture occurred as deep lateral slip on the Enriquillo-Plantain Garden Fault (EPGF). The earthquake triggered >4490 landslides, mainly shallow, disrupted rock falls, debris-soil falls and slides, and a few lateral spreads, over an area of similar to 2150 km(2). The regional distribution of these slope failures defies those of most similar earthquake-triggered landslide episodes reported previously. Most of the coseismic landslides did not proliferate in the hanging wall of the main rupture, but clustered instead at the junction of the blind Leogane and EPGF ruptures, where topographic relief and hillslope steepness are above average. Also, low-relief areas subjected to high coseismic uplift were prone to lesser hanging wall slope instability than previous studies would suggest. We argue that a combined effect of complex rupture dynamics and topography primarily control this previously rarely documented landslide pattern. Compared to recent thrust fault-earthquakes of similar magnitudes elsewhere, we conclude that lower static stress drop, mean fault displacement, and blind ruptures of the 2010 Haiti earthquake resulted in fewer, smaller, and more symmetrically distributed landslides than previous studies would suggest. Our findings caution against overly relying on across-the-board models of slope stability response to seismic ground shaking. (C) 2012 Elsevier B.V. All rights reserved.
Any understanding of sediment routing from mountain belts to their forelands and offshore sinks remains incomplete without estimates of intermediate storage that decisively buffers sediment yields from erosion rates, attenuates water and sediment fluxes, and protects underlying bedrock from incision. We quantify for the first time the sediment stored in > 38000 mainly postglacial Himalayan valley fills, based on an empirical volume-area scaling of valley-fill outlines automatically extracted from digital topographic data. The estimated total volume of 690(+452/-242) km(3) is mostly contained in few large valley fills > 1 km(3), while catastrophic mass wasting adds another 177(31) km(3). Sediment storage volumes are highly disparate along the strike of the orogen. Much of the Himalaya's stock of sediment is sequestered in glacially scoured valleys that provide accommodation space for similar to 44% of the total volume upstream of the rapidly exhuming and incising syntaxes. Conversely, the step-like long-wave topography of the central Himalayas limits glacier extent, and thus any significant glacier-derived storage of sediment away from tectonic basins. We show that exclusive removal of Himalayan valley fills could nourish contemporary sediment flux from the Indus and Brahmaputra basins for > 1 kyr, though individual fills may attain residence times of > 100 kyr. These millennial lag times in the Himalayan sediment routing system may sufficiently buffer signals of short-term seismic as well as climatic disturbances, thus complicating simple correlation and interpretation of sedimentary archives from the Himalayan orogen, its foreland, and its submarine fan systems. (C) 2013 Elsevier B.V. All rights reserved.
Natural and human-induced erosion supplies high amounts of soil organic carbon (OC) to terrestrial drainage networks. Yet OC fluxes in rivers were considered in global budgets only recently. Modern estimates of annual carbon burial in inland river sediments of 0.6 Gt C, or 22% of C transferred from terrestrial ecosystems to river channels, consider only lakes and reservoirs and disregard any long-term carbon burial in hillslope or floodplain sediments. Here we present the first assessment of sediment-bound OC storage in Central Europe from a synthesis of similar to 1500 Holocene hillslope and floodplain sedimentary archives. We show that sediment storage increases with drainage-basin size due to more extensive floodplains in larger river basins. However, hillslopes retain hitherto unrecognized high amounts of eroded soils at the scale of large river basins such that average agricultural erosion rates during the Holocene would have been at least twice as high as reported previously. This anthropogenic hillslope sediment storage exceeds floodplain storage in drainage basins <10(5) km(2), challenging the notion that floodplains are the dominant sedimentary sinks. In terms of carbon burial, OC concentrations in floodplains exceed those on hillslopes, and net OC accumulation rates in floodplains (0.70.2 g C m(-2)a(-1)) surpass those on hillslopes (0.40.1 g C m(-2)a(-1)) over the last 7500 years. We conclude that carbon burial in floodplains and on hillslopes in Central Europe exceeds terrestrial carbon storage in lakes and reservoirs by at least 2 orders of magnitude and should thus be considered in continental carbon budgets.
Quantifying volumes and rates of delivery of terrestrial sediment from island arcs to subduction zones is indispensable for refining estimates of the thickness of trench fills that may eventually control the location and timing of submarine landslides and tsunami-generating mega-earthquakes. Despite these motivating insights, knowledge about the rates of erosion and sediment export from the Japanese islands to their Pacific subduction zones remains patchy regardless of the increasing availability of highly resolved data on surface deformation, climate, geology, and topography. Traditionally, natural erosion rates across the island arc have been estimated from regression of topographic catchment metrics and reservoir sedimentation rates that were recorded over several years to decades. We review current research in this context, correct for a systematic bias in one of the most widely used predictions, and present new estimates of decadal to millennial-scale erosion rates of Japan's terrestrial inner forearc. We draw on several independent and unprecedented inventories of mass wasting, reservoir sedimentation, and concentrations of cosmogenic Be-10 in river sands. We find that natural Be-10-derived denudation rates of several mm yr(-1) in the Japanese Alps have been sustained over several centuries to millennia, and are, within error, roughly consistent with sediment yields inferred from artificial reservoir sedimentation. Local exceptions may likely result from release of sediment storage or regional landsliding episodes that trigger transient sediment pulses. Our synopsis further reveals that catchments draining Japan's eastern seaboard differ distinctly in their tectonic, lithological, topographic, and climatic characteristics between the Tohoku, Japanese Alps, and Nankai inner forearc segments, which is underscored by a marked asymmetric pattern of erosion rates along the island arc. Erosion rates are highest (up to at least 3 mm yr(-1)) in the Japanese Alps that mark the collision of two subduction zones, where high topographic relief, hillslope and bedrock-channel steepness foster rapid denudation by mass wasting. Comparable, if slightly lower, erosion rates characterise the Nankai inner forearc in southwest Japan, most likely due to higher typhoon-driven rainfall totals and variability rather than its high topographic relief. In contrast, our estimated erosion and flux rates are lowest in the Tohoku inner forearc catchments that feed sediment into the Japan Trench. We conclude that collisional mountain building of the Japanese Alps drives some of the highest erosion rates in the island arc despite similar uplift and precipitation controls in southwest Japan. We infer that, prior to extensive river damming, reservoir construction, and coastal works, the gross of Japan's total sediment export to the Pacific Ocean entered the accretionary margin of the Nankai Trough as opposed to the comparatively sediment-starved Japan Trench. Compared to documented contemporary rates of sediment flux from mountainous catchments elsewhere in the Pacific, the rivers draining Japan's inner forearc take an intermediate position despite high relief, steep slopes, very high seismicity, and frequent rainstorms. However, the average rates of millennial-scale denudation in the Japanese Alps particularly are amongst the highest reported worldwide.
Local mismatches between these late Holocene and modern rates emphasise the anthropogenic fingerprint on sediment retention that may have significantly reduced the island arc's mass flux to its subduction zones, as is the case elsewhere in east and southeast Asia. (C) 2014 Elsevier B.V. All rights reserved.
Late Quaternary valley infill and dissection in the Indus River, western Tibetan Plateau margin
(2014)
The Indus, one of Earth's major rivers, drains large parts of the NW Himalaya and the Transhimalayan ranges that form part of the western Tibetan Plateau margin. In the western Himalayan syntaxis, where local topographic relief exceeds 7 km, the Indus has incised a steep bedrock gorge at rates of several mm yr(-1). Upstream, however, the upper Indus and its tributaries alternate between bedrock gorges and broad alluvial flats flanked by the Ladakh and Zanskar ranges. We review the late Quaternary valley history in this region with a focus on the confluence of the Indus and Zanskar Rivers, where vast alluvial terrace staircases and lake sediments record major episodes of aggradation and incision. New absolute dating of high-level fluvial terrace remnants using cosmogenic Be-10, optically and infrared stimulated luminescence (OSL, IRSL) indicates at least two phases of late Quaternary valley infilling. These phases commenced before similar to 200 ka and similar to 50-20 ka, judging from terrace treads stranded >150 m and similar to 30-40 m above modern river levels, respectively. Numerous stacks of lacustrine sediments that straddle the Indus River >200 km between the city of Leh and the confluence with the Shyok River share a distinct horizontal alignment. Constraints from IRSL samples of lacustrine sequences from the Leh-Spituk area reveal a protracted lake phase from >177 ka to 72 ka, locally accumulating >50-m thick deposits. In the absence of tectonic faulting, major lithological differences, and stream capture, we attribute the formation of this and other large lakes in the region to natural damming by large landslides, glaciers, and alluvial fans. The overall patchy landform age constraints from earlier studies can be reconciled by postulating a major deglacial control on sediment flux, valley infilling, and subsequent incision that has been modulated locally by backwater effects of natural damming. While comparison with Pleistocene monsoon proxies reveals no obvious correlation, a lateor post-glacial sediment pulse seems a more likely source of this widespread sedimentation that has partly buried the dissected bedrock topography. Overall, the long residence times of fluvial, alluvial and lacustrine deposits in the region (>500 ka) support previous studies, but remain striking given the dominantly steep slopes and deeply carved valleys that characterise this high-altitude mountain desert. Recalculated late Quaternary rates of fluvial bedrock incision in the Indus and Zanskar of 1.5 +/- 0.2 mm yr(-1) are at odds with the longevity of juxtaposed valley-fill deposits, unless a lack of decisive lateral fluvial erosion helps to preserve these late Pleistocene sedimentary archives. We conclude that alternating, similar to 10(4)-yr long, phases of massive infilling and incision have dominated the late Quaternary history of the Indus valley below the western Tibetan Plateau margin. (C) 2014 Elsevier Ltd. All rights reserved.
The 2002 M-w 7.9 Denali Fault earthquake, Alaska, provides an unparalleled opportunity to investigate in quantitative detail the regional hillslope mass-wasting response to strong seismic shaking in glacierized terrain. We present the first detailed inventory of similar to 1580 coseismic slope failures, out of which some 20% occurred above large valley glaciers, based on mapping from multi-temporal remote sensing data. We find that the Denali earthquake produced at least one order of magnitude fewer landslides in a much narrower corridor along the fault ruptures than empirical predictions for an M 8 earthquake would suggest, despite the availability of sufficiently steep and dissected mountainous topography prone to frequent slope failure. In order to explore potential controls on the reduced extent of regional coseismic landsliding we compare our data with inventories that we compiled for two recent earthquakes in periglacial and formerly glaciated terrain, i.e. at Yushu, Tibet (M-w 6.9, 2010), and Aysen Fjord, Chile (2007 M-w 6.2). Fault movement during these events was, similarly to that of the Denali earthquake, dominated by strike-slip offsets along near-vertical faults. Our comparison returns very similar coseismic landslide patterns that are consistent with the idea that fault type, geometry, and dynamic rupture process rather than widespread glacier cover were among the first-order controls on regional hillslope erosional response in these earthquakes. We conclude that estimating the amount of coseismic hillslope sediment input to the sediment cascade from earthquake magnitude alone remains highly problematic, particularly if glacierized terrain is involved. (C) 2014 Elsevier Ltd. All rights reserved.
Dealing with predicted increases in extreme weather conditions due to climate change requires robust knowledge about controls on rainfall-triggered landslides. We explore relationships between rainfall and landslide size throughout the Japanese archipelago. We test whether the total volume of landslides can be predicted directly from rainfall totals, intensity, and duration using a nationwide inventory of 4744 rainfall-triggered landslides recorded from A.D. 2001 to 2011. We find that larger landslides were more abundant at the expense of smaller ones when total, maximum, and mean rainfall intensity exceeded similar to 250 mm, similar to 35 mm/h, and similar to 4 mm/h, respectively. Frequency distributions of these rainfall parameters are peaked and heavily skewed. Yet neither the most frequent nor the most extreme values of these rainfall metrics coincide consistently with the maximum landslide volumes. A striking decrease of landslide volumes at both mean and maximum rainfall intensity, as well as duration, points to an exhaustion in hillslope geomorphic response regardless of sample size, landslide type, mobilized volume, dominant lithology, or reporting bias. Our results underscore substantial offsets between the peaks of rainfall metrics and maximum associated landslide volumes, thus complicating straightforward estimates of geomorphic work from metrics of rainstorm magnitude or frequency. Only the rainfall total appears to be a suitable monotonic predictor of landslide volumes mobilized during typhoons and frontal storms.
The growing body of research on large-scale mass wasting events so far has only scarcely investigated the sedimentology of chaotic deposits from non-volcanic terrestrial landslides such that any overarching and systematic terminological framework remains elusive. Yet recent work has emphasized the need for better understanding the internal structure and composition of rockslide deposits as a means to characterise the mechanics during the final stages of runout and emplacement. We offer a comprehensive overview on the occurrence of rock fragmentation and frictional melt both at different geographic locations, and different sections within large (>10(6) m(3)) rockslide masses. We argue that exposures of pervasively fragmented and interlocked jigsaw-cracked rock masses; basal melange containing rip-up clasts and phantom blocks; micro-breccia; and thin bands of basal frictionite are indispensable clues for identifying deposits from giant rockslides that may remain morphologically inconspicuous otherwise. These sedimentary assemblages are diagnostic tools for distinguishing large rockslide debris from macro and microscopically similar glacial deposits, tectonic fault-zone breccias, and impact breccias, and thus help avoid palaeoclimatic and tectonic misinterpretations, let alone misestimates of the hazard from giant rockslides. Moreover, experimental results from Mossbauer spectroscopy of frictionite samples support visual interpretations of thin sections, and demonstrate that short-lived (<10 s) friction-induced partial melting at temperatures >1500 degrees C in the absence of water occurred at the base of several giant moving rockslides. This finding supports previous theories of dry excess runout accompanied by comminution of rock masses down to gm-scale, and indicates that catastrophic motion of large fragmenting rock masses does not require water as a potential lubricant.
The Norwegian traffic network is impacted by about 2000 landslides, avalanches, and debris flows each year that incur high economic losses. Despite the urgent need to mitigate future losses, efforts to locate potential debris flow source areas have been rare at the regional scale. We tackle this research gap by exploring a minimal set of possible topographic predictors of debris flow initiation that we input to a Weights-of-Evidence (WofE) model for mapping the regional susceptibility to debris flows in western Norway. We use an inventory of 429 debris flows that were recorded between 1979 and 2008, and use the terrain variables of slope, total curvature, and contributing area (flow accumulation) to compute the posterior probabilities of local debris flow occurrence. The novelty of our approach is that we quantify the uncertainties in the WofE approach arising from different predictor classification schemes and data input, while estimating model accuracy and predictive performance from independent test data. Our results show that a percentile-based classification scheme excels over a manual classification of the predictor variables because differing abundances in manually defined bins reduce the reliability of the conditional independence tests, a key, and often neglected, prerequisite for the WofE method. The conditional dependence between total curvature and flow accumulation precludes their joint use in the model. Slope gradient has the highest true positive rate (88%), although the fraction of area classified as susceptible is very large (37%). The predictive performance, i.e. the reduction of false positives, is improved when combined with either total curvature or flow accumulation. Bootstrapping shows that the combination of slope and flow accumulation provides more reliable predictions than the combination of slope and total curvature, and helps refining the use of slope-area plots for identifying morphometric fingerprints of debris flow source areas, an approach used outside the field of landslide susceptibility assessments.
Lacustrine sediments have been widely used to investigate past climatic and environmental changes on millennial to seasonal time scales. Sedimentary archives of lakes in mountainous regions may also record non-climatic events such as earthquakes. We argue herein that a set of 64 annual laminae couplets reconciles a stratigraphically inconsistent accelerator mass spectrometry (AMS) C-14 chronology in a similar to 4-m-long sediment core from Lake Mengda, in the north-eastern Tibetan Plateau. The laminations suggest the lake was formed by a large landslide, triggered by the 1927 Gulang earthquake (M = 8.0). The lake sediment sequence can be separated into three units based on lithologic, sedimentary, and isotopic characteristics. Starting from the bottom of the sequence, these are: (1) unweathered, coarse, sandy valley-floor deposits or landslide debris that pre-date the lake, (2) landslide-induced, fine-grained soil or reworked landslide debris with a high organic content, and (3) lacustrine sediments with low organic content and laminations. These annual laminations provide a high-resolution record of anthropogenic and environmental changes during the twentieth century, recording enhanced sediment input associated with two phases of construction activities. The high mean sedimentation rates of up to 4.8 mm year(-1) underscore the potential for reconstructing such distinct sediment pulses in remote, forested, and seemingly undisturbed mountain catchments.
Many mountain belts sustain prolonged snow cover for parts of the year, although enquiries into rates of erosion in these landscapes have focused almost exclusively on the snow-free periods. This raises the question of whether annual snow cover contributes significantly to modulating rates of erosion in high-relief terrain. In this context, the sudden release of snow avalanches is a frequent and potentially relevant process, judging from the physical damage to subalpine forest ecosystems, and the amount of debris contained in avalanche deposits. To quantitatively constrain this visual impression and to expand the sparse literature, we sampled sediment concentrations of n = 28 river-spanning snow-avalanche deposits (snow bridges) in the area around Davos, eastern Swiss Alps, and inferred an orders-of-magnitude variability in specific fine sediment and organic carbon yields (1.8 to 830 t km(-2) yr(-1), and 0.04 to 131 tC km(-2) yr(-1), respectively). A Monte Carlo simulation demonstrates that, with a minimum of free parameters, such variability is inherent to the geometric scaling used for computing specific yields. Moreover, the widely applied method of linearly extrapolating plot scale sample data may be prone to substantial under- or overestimates. A comparison of our inferred yields with previously published work demonstrates the relevance of wet snow avalanches as prominent agents of soil erosion and transporters of biogeochemical constituents to mountain rivers. Given that a number of snow bridges persisted below the insulating debris cover well into the summer months, snow-avalanche deposits also contribute to regulating in-channel sediment and organic debris storage on seasonal timescales. Finally, our results underline the potential shortcomings of neglecting erosional processes in the winter and spring months in mountainous terrain subjected to prominent snow cover.
Modern natural hazards research requires dealing with several uncertainties that arise from limited process knowledge, measurement errors, censored and incomplete observations, and the intrinsic randomness of the governing processes. Nevertheless, deterministic analyses are still widely used in quantitative hazard assessments despite the pitfall of misestimating the hazard and any ensuing risks.
In this paper we show that Bayesian networks offer a flexible framework for capturing and expressing a broad range of uncertainties encountered in natural hazard assessments. Although Bayesian networks are well studied in theory, their application to real-world data is far from straightforward, and requires specific tailoring and adaptation of existing algorithms. We offer suggestions as how to tackle frequently arising problems in this context and mainly concentrate on the handling of continuous variables, incomplete data sets, and the interaction of both. By way of three case studies from earthquake, flood, and landslide research, we demonstrate the method of data-driven Bayesian network learning, and showcase the flexibility, applicability, and benefits of this approach.
Our results offer fresh and partly counterintuitive insights into well-studied multivariate problems of earthquake-induced ground motion prediction, accurate flood damage quantification, and spatially explicit landslide prediction at the regional scale. In particular, we highlight how Bayesian networks help to express information flow and independence assumptions between candidate predictors. Such knowledge is pivotal in providing scientists and decision makers with well-informed strategies for selecting adequate predictor variables for quantitative natural hazard assessments.
The Indus River, one of Asia's premier rivers, drains the western Tibetan Plateau and the Nanga Parbat syntaxis. These two areas juxtapose some of the lowest and highest topographic relief and commensurate denudation rates in the Himalaya-Tibet orogen, respectively, yet the spatial pattern of denudation rates upstream of the syntaxis remains largely unclear, as does the way in which major rivers drive headward incision into the Tibetan Plateau. We report a new inventory of Be-10-based basinwide denudation rates from 33 tributaries flanking the Indus River along a 320 km reach across the western Tibetan Plateau margin. We find that denudation rates of up to 110 mm k.y.(-1) in the Ladakh and Zanskar Ranges systematically decrease eastward to 10 mm k.y.(-1) toward the Tibetan Plateau. Independent results from bulk petrographic and heavy mineral analyses support this denudation gradient. Assuming that incision along the Indus exerts the base-level control on tributary denudation rates, our data show a systematic eastward decrease of landscape downwearing, reaching its minimum on the Tibetan Plateau. In contrast, denudation rates increase rapidly 150-200 km downstream of a distinct knick-point that marks the Tibetan Plateau margin in the Indus River longitudinal profile. We infer that any vigorous headward incision and any accompanying erosional waves into the interior of the plateau mostly concerned reaches well below this plateau margin. Moreover, reported long-term (>10(6) yr) exhumation rates from low-temperature chronometry of 0.1-0.75 mm yr(-1) consistently exceed our Be-10-derived denudation rates. With averaging time scales of 10(3)-10(4) yr for our denudation data, we report postglacial rates of downwearing in a tectonically idle landscape. To counterbalance this apparent mismatch, denudation rates must have been higher in the Quaternary during glacial-interglacial intervals.
Deforestation is a prominent anthropogenic cause of erosive overland flow and slope instability, boosting rates of soil erosion and concomitant sediment flux. Conventional methods of gauging or estimating post-logging sediment flux often focus on annual timescales but overlook potentially important process response on shorter intervals immediately following timber harvest. We resolve such dynamics with non-parametric quantile regression forests (QRF) based on high-frequency (3 min) discharge measurements and sediment concentration data sampled every 30-60 min in similar-sized (similar to 0.1 km(2)) forested Chilean catchments that were logged during either the rainy or the dry season. The method of QRF builds on the random forest algorithm, and combines quantile regression with repeated random sub-sampling of both cases and predictors. The algorithm belongs to the family of decision-tree classifiers, which allow quantifying relevant predictors in high-dimensional parameter space. We find that, where no logging occurred, similar to 80% of the total sediment load was transported during extremely variable runoff events during only 5% of the monitoring period. In particular, dry-season logging dampened the relative role of these rare, extreme sediment-transport events by increasing load efficiency during more efficient moderate events. We show that QRFs outperform traditional sediment rating curves (SRCs) in terms of accurately simulating short-term dynamics of sediment flux, and conclude that QRF may reliably support forest management recommendations by providing robust simulations of post-logging response of water and sediment fluxes at high temporal resolution.
Limits to lichenometry
(2015)
Lichenometry is a straightforward and inexpensive method for dating Holocene rock surfaces. The rationale is that the diameter of the largest lichen scales with the age of the originally fresh rock surface that it colonised. The success of the method depends on finding the largest lichen diameters, a suitable lichen-growth model, and a robust calibration curve. Recent critique of the method motivates us to revisit the accuracy and uncertainties of lichenometry. Specifically, we test how well lichenometry is capable of resolving the ages of different lobes of large active rock glaciers in the Kyrgyz Tien Shan. We use a bootstrapped quantile regression to calibrate local growth curves of Xanthoria elegans, Aspicilia tianshanica, and Rhizocarpon geographicum, and report a nonlinear decrease in dating accuracy with increasing lichen diameter. A Bayesian type of an analysis of variance demonstrates that our calibration allows discriminating credibly between rock-glacier lobes of different ages despite the uncertainties tied to sample size and correctly identifying the largest lichen thalli. Our results also show that calibration error grows with lichen size, so that the separability of rock-glacier lobes of different ages decreases, while the tendency to assign coeval ages increases. The abundant young (<200 yr) specimen of fast-growing X elegans are in contrast with the fewer, slow-growing, but older (200-1500 yr) R. geographicum and A. tianshanica, and record either a regional reactivation of lobes in the past 200 years, or simply a censoring effect of lichen mortality during early phases of colonisation. The high variance of lichen sizes captures the activity of rock-glacier lobes, which is difficult to explain by regional climatic cooling or earthquake triggers alone. Therefore, we caution against inferring palaeoclimatic conditions from the topographic position of rock-glacier lobes. We conclude that lichenometry works better as a tool for establishing a relative, rather than an absolute, chronology of rock-glacier lobes in the northern Tien Shan. (C) 2015 Elsevier Ltd. All rights reserved.
The 2008 eruption of Chaiten volcano in southern Chile severely impacted several densely forested river catchments by supplying excess pyroclastic sediment to the channel networks. Our aim is to substantiate whether and how channel geometry and forest stands changed in the Rayas River following the sudden input of pyroclastic sediment. We measured the resulting changes to channel geometry and riparian forest stands along 17.6 km of the impacted gravel-bed Rayas River (294 km(2)) from multiple high-resolution satellite images, aerial photographs, and fieldwork to quantify yield volume characteristics of the forest stands. Limited channel changes during the last 60 years before the eruption reflect a dynamic equilibrium condition of the river corridor, despite the high annual precipitation and the sediment supply from Chaiten and Michinmahuida volcanoes in the headwaters. Images taken in 1945, 2004, and 2005 show that total size of the vegetated channel islands nearly doubled between 1945 and 2004 and remained unchanged between 2004 and 2005. Pyroclastic sediment entering the Rayas River after the 2008 eruption caused only minor average channel widening (7%), but killed all island vegetation in the study reach. Substantial shifts in the size distribution of in-channel vegetation patches reflect losses in total island area of 46% from 2005 to 2009 and an additional 34% from 2009 to 2012. The estimated pulsed release of organic carbon into the channel, mainly in the form of large wood from obliterated island and floodplain forests, was 78-400 tC/km/y and surpasses most documented yields from small mountainous catchments with similar rainfall, forest cover, and disturbance history, while making up between 20% and 60% of the annual carbon burial rate of fluvial sediments in the northern Patagonian fjords. We conclude that the carbon footprint of the 2008 Chaiten eruption on the Rayas River was more significant than the measured geomorphic impacts on channel geometry for the first five years following disturbance. The modest post-eruptive geomorphic response in this river is a poor indicator of its biogeochemical response. (C) 2015 Elsevier B.V. All rights reserved.
Mass wasting is an important process for denuding hillslopes and lowering ridge crests in active mountain belts such as the Himalaya-Karakoram ranges (HKR). Such a high-relief landscape is likely to be at its mechanical threshold, maintained by competing rapid rock uplift, river incision, and pervasive slope failure. We introduce excess topography, Z(E), for quantifying potentially unstable rock-mass volumes inclined at angles greater than a specified threshold angle. We find that Z(E) peaks along major fluvial and glacial inner gorges, which is also where the majority of 492 large (>0.1 km(2)) rock-slope failures occur in the Himalaya's largest cluster of documented Pleistocene to Holocene bedrock landslides. Our data reveal that bedrock landslides in the HKR chiefly detached from near or below the median elevation, whereas glaciers and rock glaciers occupy higher-elevation bands almost exclusively. Less than 10% of the area of the HKR is upslope of glaciers, such that possible censoring of evidence of large bedrock landslides above the permanent snow line barely affects this finding. Bedrock landslides appear to preferentially undermine topographic relief in response to fluvial and glacial incision along inner gorges, unless more frequent and smaller undetected failures, or rigorous (peri-)glacial erosion, compensate for this role at higher elevation. Either way, the distinct patterns of excess topography and large bedrock landsliding in the HKR juxtapose two stacked domains of landslide and (peri-)glacial erosion that may respond to different time scales of perturbation. Our findings call for more detailed analysis of vertical erosional domains and their geomorphic coupling in active mountain belts.
In this study, we investigate how immersive 3D geovisualization can be used in higher education. Based on MacEachren and Kraak's geovisualization cube, we examine the usage of immersive 3D geovisualization and its usefulness in a research-based learning module on flood risk, called GEOSimulator. Results of a survey among participating students reveal benefits, such as better orientation in the study area, higher interactivity with the data, improved discourse among students and enhanced motivation through immersive 3D geovisualization. This suggests that immersive 3D visualization can effectively be used in higher education and that 3D CAVE settings enhance interactive learning between students.
Roads at risk
(2015)
Globalisation and interregional exchange of people, goods, and services has boosted the importance of and reliance on all kinds of transport networks. The linear structure of road networks is especially sensitive to natural hazards. In southern Norway, steep topography and extreme weather events promote frequent traffic disruption caused by debris flows. Topographic susceptibility and trigger frequency maps serve as input into a hazard appraisal at the scale of first-order catchments to quantify the impact of debris flows on the road network in terms of a failure likelihood of each link connecting two network vertices, e.g. road junctions. We compute total additional traffic loads as a function of traffic volume and excess distance, i.e. the extra length of an alternative path connecting two previously disrupted network vertices using a shortest-path algorithm. Our risk metric of link failure is the total additional annual traffic load, expressed as vehicle kilometres, because of debris-flow-related road closures. We present two scenarios demonstrating the impact of debris flows on the road network and quantify the associated path-failure likelihood between major cities in southern Norway. The scenarios indicate that major routes crossing the central and north-western part of the study area are associated with high link-failure risk. Yet options for detours on major routes are manifold and incur only little additional costs provided that drivers are sufficiently well informed about road closures. Our risk estimates may be of importance to road network managers and transport companies relying on speedy delivery of services and goods.
Roads at risk
(2015)
Globalisation and interregional exchange of people, goods, and services has boosted the importance of and reliance on all kinds of transport networks. The linear structure of road networks is especially sensitive to natural hazards. In southern Norway, steep topography and extreme weather events promote frequent traffic disruption caused by debris flows. Topographic susceptibility and trigger frequency maps serve as input into a hazard appraisal at the scale of first-order catchments to quantify the impact of debris flows on the road network in terms of a failure likelihood of each link connecting two network vertices, e.g. road junctions. We compute total additional traffic loads as a function of traffic volume and excess distance, i.e. the extra length of an alternative path connecting two previously disrupted network vertices using a shortest-path algorithm. Our risk metric of link failure is the total additional annual traffic load, expressed as vehicle kilometres, because of debris-flow-related road closures. We present two scenarios demonstrating the impact of debris flows on the road network and quantify the associated path-failure likelihood between major cities in southern Norway. The scenarios indicate that major routes crossing the central and north-western part of the study area are associated with high link-failure risk. Yet options for detours on major routes are manifold and incur only little additional costs provided that drivers are sufficiently well informed about road closures. Our risk estimates may be of importance to road network managers and transport companies relying on speedy delivery of services and goods.
Sediments in rivers record the dynamics of erosion processes. While bulk sediment fluxes are easily and routinely obtained, sediment caliber remains underexplored when inferring erosion mechanisms. Yet sediment grain size distributions may be the key to discriminating their origin. We have studied grain size-specific suspended sediment fluxes in the Kali Gandaki, a major trans-Himalayan river. Two strategically located gauging stations enable tracing of sediment caliber on either side of the Himalayan orographic barrier. The data show that fine sediment input into the northern headwaters is persistent, while coarse sediment comes from the High Himalayas during the summer monsoon. A temporally matching landslide inventory similarly indicates the prominence of monsoon-driven hillslope mass wasting. Thus, mechanisms of sediment supply can leave strong traces in the fluvial caliber, which could project well beyond the mountain front and add to the variability of the sedimentary record of orogen erosion.
The use of topographic metrics for estimating the susceptibility to, and reconstructing the characteristics of, debris flows has a long research tradition, although largely devoted to humid mountainous terrain. The exceptional 2010 monsoonal rainstorms in the high-altitude mountain desert of Ladakh and Zanskar, NW India, were a painful reminder of how susceptible arid regions are to rainfall-triggered flash floods, landslides, and debris flows. The rainstorms of August 4-6 triggered numerous debris flows, killing 182 people, devastating 607 houses, and more than 10 bridges around Ladakh's capital of Leh. The lessons from this disaster motivated us to revisit methods of predicting (a) flow parameters such as peak discharge and maximum velocity from field and remote sensing data, and (b) the susceptibility to debris flows from catchment morphometry. We focus on quantifying uncertainties tied to these approaches. Comparison of high-resolution satellite images pre- and post-dating the 2010 rainstorm reveals the extent of damage and catastrophic channel widening. Computations based on these geomorphic markers indicate maximum flow velocities of 1.6-6.7 m s(-1) with runout of up to similar to 10 km on several alluvial fans that sustain most of the region's settlements. We estimate median peak discharges of 310-610 m(3) s(-1), which are largely consistent with previous estimates. Monte Carlo-based error propagation for a single given flow-reconstruction method returns a variance in discharge similar to one derived from juxtaposing several different flow reconstruction methods. We further compare discriminant analysis, classification tree modelling, and Bayesian logistic regression to predict debris-flow susceptibility from morphometric variables of 171 catchments in the Ladakh Range. These methods distinguish between fluvial and debris flow-prone catchments at similar success rates, but Bayesian logistic regression allows quantifying uncertainties and relationships between potential predictors. We conclude that, in order to be robust and reliable, morphometric reconstruction of debris-flow properties and susceptibility requires careful assessment and reporting of errors and uncertainties. (C) 2015 Elsevier B.V. All rights reserved.
Himalayan water resources attract a rapidly growing number of hydroelectric power projects (HPP) to satisfy Asia's soaring energy demands. Yet HPP operating or planned in steep, glacier-fed mountain rivers face hazards of glacial lake outburst floods (GLOFs) that can damage hydropower infrastructure, alter water and sediment yields, and compromise livelihoods downstream. Detailed appraisals of such GLOF hazards are limited to case studies, however, and a more comprehensive, systematic analysis remains elusive. To this end we estimate the regional exposure of 257 Himalayan HPP to GLOFs, using a flood-wave propagation model fed by Monte Carlo-derived outburst volumes of >2300 glacial lakes. We interpret the spread of thus modeled peak discharges as a predictive uncertainty that arises mainly from outburst volumes and dam-breach rates that are difficult to assess before dams fail. With 66% of sampled HPP are on potential GLOF tracks, up to one third of these HPP could experience GLOF discharges well above local design floods, as hydropower development continues to seek higher sites closer to glacial lakes. We compute that this systematic push of HPP into headwaters effectively doubles the uncertainty about GLOF peak discharge in these locations. Peak discharges farther downstream, in contrast, are easier to predict because GLOF waves attenuate rapidly. Considering this systematic pattern of regional GLOF exposure might aid the site selection of future Himalayan HPP. Our method can augment, and help to regularly update, current hazard assessments, given that global warming is likely changing the number and size of Himalayan meltwater lakes.
Geomorphic footprints of past large Himalayan earthquakes are elusive, although they are urgently needed for gauging and predicting recovery times of seismically perturbed mountain landscapes. We present evidence of catastrophic valley infill following at least three medieval earthquakes in the Nepal Himalaya. Radiocarbon dates from peat beds, plant macrofossils, and humic silts in fine-grained tributary sediments near Pokhara, Nepal’s second-largest city, match the timing of nearby M > 8 earthquakes in ~1100, 1255, and 1344 C.E. The upstream dip of tributary valley fills and x-ray fluorescence spectrometry of their provenance rule out local sources. Instead, geomorphic and sedimentary evidence is consistent with catastrophic fluvial aggradation and debris flows that had plugged several tributaries with tens of meters of calcareous sediment from a Higher Himalayan source >60 kilometers away.
Himalayan water resources attract a rapidly growing number of hydroelectric power projects (HPP) to satisfy Asia's soaring energy demands. Yet HPP operating or planned in steep, glacier-fed mountain rivers face hazards of glacial lake outburst floods (GLOFs) that can damage hydropower infrastructure, alter water and sediment yields, and compromise livelihoods downstream. Detailed appraisals of such GLOF hazards are limited to case studies, however, and a more comprehensive, systematic analysis remains elusive. To this end we estimate the regional exposure of 257 Himalayan HPP to GLOFs, using a flood-wave propagation model fed by Monte Carlo-derived outburst volumes of >2300 glacial lakes. We interpret the spread of thus modeled peak discharges as a predictive uncertainty that arises mainly from outburst volumes and dam-breach rates that are difficult to assess before dams fail. With 66% of sampled HPP are on potential GLOF tracks, up to one third of these HPP could experience GLOF discharges well above local design floods, as hydropower development continues to seek higher sites closer to glacial lakes. We compute that this systematic push of HPP into headwaters effectively doubles the uncertainty about GLOF peak discharge in these locations. Peak discharges farther downstream, in contrast, are easier to predict because GLOF waves attenuate rapidly. Considering this systematic pattern of regional GLOF exposure might aid the site selection of future Himalayan HPP. Our method can augment, and help to regularly update, current hazard assessments, given that global warming is likely changing the number and size of Himalayan meltwater lakes.
Himalayan water resources attract a rapidly growing number of hydroelectric power projects (HPP) to satisfy Asia's soaring energy demands. Yet HPP operating or planned in steep, glacier-fed mountain rivers face hazards of glacial lake outburst floods (GLOFs) that can damage hydropower infrastructure, alter water and sediment yields, and compromise livelihoods downstream. Detailed appraisals of such GLOF hazards are limited to case studies, however, and a more comprehensive, systematic analysis remains elusive. To this end we estimate the regional exposure of 257 Himalayan HPP to GLOFs, using a flood-wave propagation model fed by Monte Carlo-derived outburst volumes of >2300 glacial lakes. We interpret the spread of thus modeled peak discharges as a predictive uncertainty that arises mainly from outburst volumes and dam-breach rates that are difficult to assess before dams fail. With 66% of sampled HPP are on potential GLOF tracks, up to one third of these HPP could experience GLOF discharges well above local design floods, as hydropower development continues to seek higher sites closer to glacial lakes. We compute that this systematic push of HPP into headwaters effectively doubles the uncertainty about GLOF peak discharge in these locations. Peak discharges farther downstream, in contrast, are easier to predict because GLOF waves attenuate rapidly. Considering this systematic pattern of regional GLOF exposure might aid the site selection of future Himalayan HPP. Our method can augment, and help to regularly update, current hazard assessments, given that global warming is likely changing the number and size of Himalayan meltwater lakes.
The Tibetan Plateau, the world's largest orogenic plateau, hosts thousands of lakes that play prominent roles as water resources, environmental archives, and sources of natural hazards such as glacier lake outburst floods. Previous studies have reported that the size of lakes on the Tibetan Plateau has changed rapidly in recent years, possibly because of atmospheric warming. Tracking these changes systematically with remote sensing data is challenging given the different spectral signatures of water, the potential for confusing lakes with glaciers, and difficulties in classifying frozen or partly frozen lakes. Object-based image analysis (OBIA) offers new opportunities for automated classification in this context, and we have explored this method for mapping lakes from LANDSAT images and Shuttle Radar Topography Mission (SRTM) elevation data. We tested our algorithm for most of the Tibetan Plateau, where lakes in tectonic depressions or blocked by glaciers and sediments have different surface colours and seasonal ice cover in images obtained in 1995 and 2015. We combined a modified normalised difference water index (MNDWI) with OBIA and local topographic slope data in order to classify lakes with an area > 10 km(2). Our method derived 323 water bodies, with a total area of 31,258 km(2), or 2.6% of the study area (in 2015). The same number of lakes had covered only 24,892 km(2) in 1995; lake area has increased by -26% in the past two decades. The classification had estimated producer's and user's accuracies of 0.98, with a Cohen's kappa and F-score of 0.98, and may thus be a useful approximation for quantifying regional hydrological budgets. We have shown that our method is flexible and transferable to detecting lakes in diverse physical settings on several continents with similar success rates.
Fjords and old-growth forests store large amounts of organic carbon. Yet the role of episodic disturbances, particularly volcanic eruptions, in mobilizing organic carbon in fjord landscapes covered by temperate rainforests remains poorly quantified. To this end, we estimated how much wood and soils were flushed to nearby fjords following the 2008 eruption of Chaiten volcano in south-central Chile, where pyroclastic sediments covered >12km(2) of pristine temperate rainforest. Field-based surveys of forest biomass, soil organic content, and dead wood transport reveal that the reworking of pyroclastic sediments delivered similar to 66,500+14,600/-14,500tC of large wood to two rivers entering the nearby Patagonian fjords in less than a decade. A similar volume of wood remains in dead tree stands and buried beneath pyroclastic deposits (similar to 79,900+21,100/-16,900tC) or stored in active river channels (5,900-10,600tC). We estimate that bank erosion mobilized similar to 132,300(+21,700)/(-30,600)tC of floodplain forest soil. Eroded and reworked forest soils have been accreting on coastal river deltas at >5mmyr(-1) since the eruption. While much of the large wood is transported out of the fjord by long-shore drift, the finer fraction from eroded forest soils is likely to be buried in the fjords. We conclude that the organic carbon fluxes boosted by rivers adjusting to high pyroclastic sediment loads may remain elevated for up to a decade and that Patagonian temperate rainforests disturbed by excessive loads of pyroclastic debris can be episodic short-lived carbon sources. Plain Language Summary Fjords and old-growth forests are important sinks of organic carbon. However, the role of volcanic eruptions in flushing organic carbon in fjord landscapes remains unexplored. Here we estimated how much forest vegetation and soils were lost to fjords following the 2008 eruption ofunknownChaiten volcano in south-central Chile. Pyroclastic sediments obliterated near-pristine temperateunknownrainforest, and the subsequent reworking of these sediments delivered in less than a decade similar to 66,000 tC of large wood to the mountain rivers, draining into the nearby Patagonian fjords. A similar volume of wood remains in dead tree stands and buried beneath pyroclastic deposits or stored in active riverunknownchannels. We estimate that similar to 130,000 tC of organic carbon-rich soil was lost to erosion, thus adding to the carbon loads. While much of the wood enters the long-shore drift in the fjord heads, the finerunknownfraction from eroded forest soils is likely to be buried in the fjords at rates that exceed regional estimates by an order of magnitude. We anticipate that these eruption-driven fluxes will remain elevated forunknownthe coming years and that Patagonian temperate rainforests episodically switch from carbon sinks to hitherto undocumented carbon sources if disturbed by explosive volcanic eruptions.
Regional snow-avalanche detection using object-based image analysis of near-infrared aerial imagery
(2017)
Snow avalanches are destructive mass movements in mountain regions that continue to claim lives and cause infrastructural damage and traffic detours. Given that avalanches often occur in remote and poorly accessible steep terrain, their detection and mapping is extensive and time consuming. Nonetheless, systematic avalanche detection over large areas could help to generate more complete and up-to-date inventories (cadastres) necessary for validating avalanche forecasting and hazard mapping. In this study, we focused on automatically detecting avalanches and classifying them into release zones, tracks, and run-out zones based on 0.25m near-infrared (NIR) ADS80-SH92 aerial imagery using an object-based image analysis (OBIA) approach. Our algorithm takes into account the brightness, the normalised difference vegetation index (NDVI), the normalised difference water index (NDWI), and its standard deviation (SDNDWI) to distinguish avalanches from other land-surface elements. Using normalised parameters allows applying this method across large areas. We trained the method by analysing the properties of snow avalanches at three 4km−2 areas near Davos, Switzerland. We compared the results with manually mapped avalanche polygons and obtained a user's accuracy of >0.9 and a Cohen's kappa of 0.79–0.85. Testing the method for a larger area of 226.3km−2, we estimated producer's and user's accuracies of 0.61 and 0.78, respectively, with a Cohen's kappa of 0.67. Detected avalanches that overlapped with reference data by >80% occurred randomly throughout the testing area, showing that our method avoids overfitting. Our method has potential for large-scale avalanche mapping, although further investigations into other regions are desirable to verify the robustness of our selected thresholds and the transferability of the method.
Moderate to large earthquakes can increase the amount of water feeding stream flows, mobilizing excess water from deep groundwater, shallow groundwater, or the vadose zone. Here we examine the regional pattern of streamflow response to the Maule M8.8 earthquake across Chile's diverse topographic and hydro-climatic gradients. We combine streamflow analyses with groundwater flow modeling and a random forest classifier, and find that, after the earthquake, at least 85 streams had a change in flow. Discharge mostly increased () shortly after the earthquake, liberating an excess water volume of >1.1 km3, which is the largest ever reported following an earthquake. Several catchments had increased discharge of >50 mm, locally exceeding seasonal streamflow discharge under undisturbed conditions. Our modeling results favor enhanced vertical permeability induced by dynamic strain as the most probable process explaining the observed changes at the regional scale. Supporting this interpretation, our random forest classification identifies peak ground velocity and elevation extremes as most important for predicting streamflow response. Given the mean recurrence interval of ∼25 yr for >M8.0 earthquakes along the Peru–Chile Trench, our observations highlight the role of earthquakes in the regional water cycle, especially in arid environments.
Am Abend des 29. Mai 2016 wurde der Ort Braunsbach im Landkreis Schwäbisch-Hall (Baden-Württemberg) von einer Sturzflut getroffen, bei der mehrere Häuser stark beschädigt oder zerstört wurden. Die Sturzflut war eine der Unwetterfolgen, die im Frühsommer 2016 vom Tiefdruckgebiet Elvira ausgelöst wurden. Der vorliegende Bericht ist der zweite Teil einer Doppelveröffentlichung, welche die Ergebnisse zur Untersuchung des Sturzflutereignisses im Rahmen des DFG-Graduiertenkollegs “Naturgefahren und Risiken in einer sich verändernden Welt” (NatRiskChange, GRK 2043/1) der Universität Potsdam präsentiert. Während Teil 1 die meteorologischen und hydrologischen Ereignisse analysiert, fokussiert Teil 2 auf die geomorphologischen Prozesse und die verursachten Gebäudeschäden. Dazu wurden Ursprung und Ausmaß des während des Sturzflutereignisses mobilisierten und in den Ort getragenen Materials untersucht. Des Weiteren wurden zu 96 betroffenen Gebäuden Daten zum Schadensgrad sowie Prozess- und Gebäudecharakteristika aufgenommen und ausgewertet. Die Untersuchungen zeigen, dass bei der Betrachtung von Hochwassergefährdung die Berücksichtigung von Sturzfluten und ihrer speziellen Charakteristika, wie hoher Feststofftransport und sprunghaftes Verhalten insbesondere in bebautem Gelände, wesentlich ist, um effektive Schutzmaßnahmen ergreifen zu können.
Regional snow-avalanche detection using object-based image analysis of near-infrared aerial imagery
(2017)
Snow avalanches are destructive mass movements in mountain regions that continue to claim lives and cause infrastructural damage and traffic detours. Given that avalanches often occur in remote and poorly accessible steep terrain, their detection and mapping is extensive and time consuming. Nonetheless, systematic avalanche detection over large areas could help to generate more complete and up-to-date inventories (cadastres) necessary for validating avalanche forecasting and hazard mapping. In this study, we focused on automatically detecting avalanches and classifying them into release zones, tracks, and run-out zones based on 0.25 m near-infrared (NIR) ADS80-SH92 aerial imagery using an object-based image analysis (OBIA) approach. Our algorithm takes into account the brightness, the normalised difference vegetation index (NDVI), the normalised difference water index (NDWI), and its standard deviation (SDNDWI) to distinguish avalanches from other land-surface elements. Using normalised parameters allows applying this method across large areas. We trained the method by analysing the properties of snow avalanches at three 4 km−2 areas near Davos, Switzerland. We compared the results with manually mapped avalanche polygons and obtained a user's accuracy of > 0.9 and a Cohen's kappa of 0.79–0.85. Testing the method for a larger area of 226.3 km−2, we estimated producer's and user's accuracies of 0.61 and 0.78, respectively, with a Cohen's kappa of 0.67. Detected avalanches that overlapped with reference data by > 80 % occurred randomly throughout the testing area, showing that our method avoids overfitting. Our method has potential for large-scale avalanche mapping, although further investigations into other regions are desirable to verify the robustness of our selected thresholds and the transferability of the method.
Regional snow-avalanche detection using object-based image analysis of near-infrared aerial imagery
(2017)
Snow avalanches are destructive mass movements in mountain regions that continue to claim lives and cause infrastructural damage and traffic detours. Given that avalanches often occur in remote and poorly accessible steep terrain, their detection and mapping is extensive and time consuming. Nonetheless, systematic avalanche detection over large areas could help to generate more complete and up-to-date inventories (cadastres) necessary for validating avalanche forecasting and hazard mapping. In this study, we focused on automatically detecting avalanches and classifying them into release zones, tracks, and run-out zones based on 0.25 m near-infrared (NIR) ADS80-SH92 aerial imagery using an object-based image analysis (OBIA) approach. Our algorithm takes into account the brightness, the normalised difference vegetation index (NDVI), the normalised difference water index (NDWI), and its standard deviation (SDNDWI) to distinguish avalanches from other land-surface elements. Using normalised parameters allows applying this method across large areas. We trained the method by analysing the properties of snow avalanches at three 4 km−2 areas near Davos, Switzerland. We compared the results with manually mapped avalanche polygons and obtained a user's accuracy of > 0.9 and a Cohen's kappa of 0.79–0.85. Testing the method for a larger area of 226.3 km−2, we estimated producer's and user's accuracies of 0.61 and 0.78, respectively, with a Cohen's kappa of 0.67. Detected avalanches that overlapped with reference data by > 80 % occurred randomly throughout the testing area, showing that our method avoids overfitting. Our method has potential for large-scale avalanche mapping, although further investigations into other regions are desirable to verify the robustness of our selected thresholds and the transferability of the method.
The northern edge of the western central Tien Shan range is bounded by the Issyk-Ata fault situated south of Bishkek, the capital of Kyrgyzstan. Contraction in this thick-skinned orogen occurs with low-strain accumulation and long earthquake recurrence intervals. In the nineteenth to twentieth centuries, a sequence of large earthquakes with magnitudes between 6.9 and 8 affected the northern Tien Shan but left nearly the entire extent of the Issyk-Ata fault unruptured. Here, the only known historic earthquake ruptured in A.D. 1885 (M6.9) along the western end of the Issyk-Ata fault. Because earthquakes in low-strain regions often tend to cluster in time and may promote failure along nearby structures, the earthquake history of the northern Tien Shan represents an exceptional structural setting for studying fault behavior affected by an intraplate earthquake sequence. We present a paleoseismological study from one site (Belek) along the Issyk-Ata fault located east of the A.D. 1885 epicentral area. Our analysis combines a range of tools, including photogrammetry, differential Global Positioning System, 3D visualization, and age modeling with different dating methods (infrared stimulated luminescence, radiocarbon, U-series) to improve the reliability of an event chronology for the trench stratigraphy and fault geometry. We were able to distinguish three different surfacerupturing paleoearthquakes; these affected the area before 10.5 +/- 1.1 cal ka B.P., at similar to 5.6 +/- 1.0 cal ka B.P., and at similar to 630 +/- 100 cal B.P., respectively. Associated paleomagnitudes for the last two earthquakes range between M6.7 and 7.4, with a cumulative slip rate of 0.7 +/- 0.32 mm/yr. We did not find evidence for the A.D. 1885 event at Belek. Our study yielded two main overall results: first, it extends the regional historic and paleoseismic record; second, the documented rupture events along the Issyk-Ata fault suggest that this fault was not affected in its entirety; instead, these events indicate segmented rupture behavior.
Rock glaciers in semiarid mountains contain large amounts of ice and might be important water stores aside from glaciers, lakes, and rivers. Yet whether and how rock glaciers interact with river channels in mountain valleys remains largely unresolved. We examine the potential for rock glaciers to block or disrupt river channels, using a new inventory of more than 2000 intact rock glaciers that we mapped from remotely sensed imagery in the Karakoram (KR), Tien Shan (TS), and Altai (ALT) mountains. We find that between 5% and 14% of the rock glaciers partly buried, blocked, diverted or constricted at least 95 km of mountain rivers in the entire study area. We use a Bayesian robust logistic regression with multiple topographic and climatic inputs to discern those rock glaciers disrupting mountain rivers from those with no obvious impacts. We identify elevation and potential incoming solar radiation (PISR), together with the size of feeder basins, as dominant predictors, so that lower-lying and larger rock glaciers from larger basins are more likely to disrupt river channels. Given that elevation and PISR are key inputs for modelling the regional distribution of mountain permafrost from the positions of rock-glacier toes, we infer that river-blocking rock glaciers may be diagnostic of non-equilibrated permafrost. Principal component analysis adds temperature evenness and wet-season precipitation to the controls that characterise rock glaciers impacting on rivers. Depending on the choice of predictors, the accuracy of our classification is moderate to good with median posterior area-under-the-curve values of 0.71-0.89. Clarifying whether rapidly advancing rock glaciers can physically impound rivers, or fortify existing dams instead, deserves future field investigation. We suspect that rock-glacier dams are conspicuous features that have a polygenetic history and encourage more research on the geomorphic coupling between permafrost lobes, river channels, and the sediment cascades of semiarid mountain belts. (c) 2018 John Wiley & Sons, Ltd.
Mountain rivers respond to strong earthquakes by rapidly aggrading to accommodate excess sediment delivered by co-seismic landslides. Detailed sediment budgets indicate that rivers need several years to decades to recover from seismic disturbances, depending on how recovery is defined. We examine three principal proxies of river recovery after earthquake-induced sediment pulses around Pokhara, Nepal's second largest city. Freshly exhumed cohorts of floodplain trees in growth position indicate rapid and pulsed sedimentation that formed a fan covering 150 km2 in a Lesser Himalayan basin with tens of metres of debris between the 11th and 15th centuries AD. Radiocarbon dates of buried trees are consistent with those of nearby valley deposits linked to major medieval earthquakes, such that we can estimate average rates of re-incision since. We combine high-resolution digital elevation data, geodetic field surveys, aerial photos, and dated tree trunks to reconstruct geomorphic marker surfaces. The volumes of sediment relative to these surfaces require average net sediment yields of up to 4200 t km–2 yr–1 for the 650 years since the last inferred earthquake-triggered sediment pulse. The lithological composition of channel bedload differs from that of local bedrock, confirming that rivers are still mostly evacuating medieval valley fills, locally incising at rates of up to 0.2 m yr–1. Pronounced knickpoints and epigenetic gorges at tributary junctions further illustrate the protracted fluvial response; only the distal portions of the earthquake-derived sediment wedges have been cut to near their base. Our results challenge the notion that mountain rivers recover speedily from earthquakes within years to decades. The valley fills around Pokhara show that even highly erosive Himalayan rivers may need more than several centuries to adjust to catastrophic perturbations. Our results motivate some rethinking of post-seismic hazard appraisals and infrastructural planning in active mountain regions.
Quaternary glaciations have repeatedly shaped large tracts of the Andean foreland. Its spectacular large glacial lakes, staircases of moraine ridges, and extensive outwash plains have inspired generations of scientists to reconstruct the processes, magnitude, and timing of ice build-up and decay at the mountain front. Surprisingly few of these studies noticed many dozens of giant (≥108 m3) mass-wasting deposits in the foreland. We report some of the world's largest terrestrial landslides in the eastern piedmont of the Patagonian Ice Sheet (PIS) along the traces of the former Lago Buenos Aires and Lago Puyerredón glacier lobes and lakes. More than 283 large rotational slides and lateral spreads followed by debris slides, earthflows, rotational and translational rockslides, complex slides and few large rock avalanches detached some 164 ± 56 km3 of material from the slopes of volcanic mesetas, lake-bounding moraines, and river-gorge walls. Many of these landslide deposits intersect with well-dated moraine ridges or former glacial-lake shorelines, and offer opportunities for relative dating of slope failure. We estimate that >60% of the landslide volume (∼96 km3) detached after the Last Glacial Maximum (LGM). Giant slope failures cross-cutting shorelines of a large Late Glacial to Early Holocene lake (“glacial lake PIS”) likely occurred during successive lake-level drop between ∼11.5 and 8 ka, and some of them are the largest hitherto documented landslides in moraines. We conclude that 1) large portions of terminal moraines can fail catastrophically several thousand years after emplacement; 2) slopes formed by weak bedrock or unconsolidated glacial deposits bordering glacial lakes can release extremely large landslides; and 3) landslides still occur in the piedmont, particularly along postglacial gorges cut in response to falling lake levels.