Refine
Language
- English (36)
Is part of the Bibliography
- yes (36)
Keywords
- erosion (4)
- seismic monitoring (3)
- 2008 Wenchuan earthquake (2)
- HVSR (2)
- correlation (2)
- environmental seismology (2)
- fatigue (2)
- floodplains (2)
- fundamental frequency (2)
- mass (2)
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.
Earthquakes impart an impressive force on epicentral landscapes, with immediate catastrophic hillslope response. However, their legacy on geomorphic process rates remains poorly constrained. We have determined the evolution of landslide rates in the epicentral areas of four intermediate to large earthquakes (M-w, 6.6-7.6). In each area, landsliding correlates with the cumulative precipitation during a given interval. Normalizing for this meteorological forcing, landslide rates have been found to peak after an earthquake and decay to background values in 1-4 yr, with the decay time scale probably proportional to the earthquake magnitude. The transient pulse of landsliding is not related to external forcing such as rainfall or aftershocks, and we tentatively attribute it to the reduction and subsequent recovery of ground strength. Observed geomorphic trends are not linked with groundwater level changes or root system damage, both of which could affect substrate strength. We propose that they are caused by reversible damage of rock mass and/or loosening of regolith. Qualitative accounts of ground cracking due to strong ground motion abound, and our observations are circumstantial evidence of its potential importance in setting landscape sensitivity to meteorological forcing after large earthquakes.
Inventories of individually delineated landslides are a key to understanding landslide physics and mitigating their impact. They permit assessment of area-frequency distributions and landslide volumes, and testing of statistical correlations between landslides and physical parameters such as topographic gradient or seismic strong motion. Amalgamation, i.e. the mapping of several adjacent landslides as a single polygon, can lead to potentially severe distortion of the statistics of these inventories. This problem can be especially severe in data sets produced by automated mapping. We present five inventories of earthquake-induced landslides mapped with different materials and techniques and affected by varying degrees of amalgamation. Errors on the total landslide volume and power-law exponent of the area-frequency distribution, resulting from amalgamation, may be up to 200 and 50 %, respectively. We present an algorithm based on image and digital elevation model (DEM) analysis, for automatic identification of amalgamated polygons. On a set of about 2000 polygons larger than 1000 m(2), tracing landslides triggered by the 1994 Northridge earthquake, the algorithm performs well, with only 2.7-3.6% incorrectly amalgamated landslides missed and 3.9-4.8% correct polygons incorrectly identified as amalgams. This algorithm can be used broadly to check landslide inventories and allow faster correction by automating the identification of amalgamation.
Inventories of individually delineated landslides are a key to understanding landslide physics and mitigating their impact. They permit assessment of area–frequency distributions and landslide volumes, and testing of statistical correlations between landslides and physical parameters such as topographic gradient or seismic strong motion. Amalgamation, i.e. the mapping of several adjacent landslides as a single polygon, can lead to potentially severe distortion of the statistics of these inventories. This problem can be especially severe in data sets produced by automated mapping. We present five inventories of earthquake-induced landslides mapped with different materials and techniques and affected by varying degrees of amalgamation. Errors on the total landslide volume and power-law exponent of the area–frequency distribution, resulting from amalgamation, may be up to 200 and 50%, respectively. We present an algorithm based on image and digital elevation model (DEM) analysis, for automatic identification of amalgamated polygons. On a set of about 2000 polygons larger than 1000 m2, tracing landslides triggered by the 1994 Northridge earthquake, the algorithm performs well, with only 2.7–3.6% incorrectly amalgamated landslides missed and 3.9–4.8% correct polygons incorrectly identified as amalgams. This algorithm can be used broadly to check landslide inventories and allow faster correction by automating the identification of amalgamation.
Large, compressional earthquakes cause surface uplift aswell as widespread mass wasting. Knowledge of their trade-off is fragmentary. Combining a seismologically consistent model of earthquake-triggered landsliding and an analytical solution of coseismic surface displacement, we assess how the mass balance of single earthquakes and earthquake sequences depends on fault size and other geophysical parameters. We find that intermediate size earthquakes (M-w 6-7.3) may cause more erosion than uplift, controlled primarily by seismic source depth and landscape steepness, and less so by fault dip and rake. Such earthquakes can limit topographic growth, but our model indicates that both smaller and larger earthquakes (M-w < 6, M-w > 7.3) systematically cause mountain building. Earthquake sequences with a Gutenberg-Richter distribution have a greater tendency to lead to predominant erosion, than repeating earthquakes of the same magnitude, unless a fault can produce earthquakes with M-w > 8 or more.
We present a new, seismologically consistent expression for the total area and volume of populations of earthquake-triggered landslides. This model builds on a set of scaling relationships between key parameters, such as landslide spatial density, seismic ground acceleration, fault length, earthquake source depth, and seismic moment. To assess the model we have assembled and normalized a catalog of landslide inventories for 40 shallow, continental earthquakes. Low landscape steepness causes systematic overprediction of the total area and volume of landslides. When this effect is accounted for, the model predicts the total landslide volume of 63% of 40 cases to within a factor 2 of the volume estimated from observations (R-2 = 0.76). The prediction of total landslide area is also sensitive to the landscape steepness, but less so than the total volume, and it appears to be sensitive to controls on the landslide size-frequency distribution, and possibly the shaking duration. Some outliers are likely associated with exceptionally strong rock mass in the epicentral area, while others may be related to seismic source complexities ignored by the model. However, the close match between prediction and estimate for about two thirds of cases in our database suggests that rock mass strength is similar in many cases and that our simple seismic model is often adequate, despite the variety of lithologies and tectonic settings covered. This makes our expression suitable for integration into landscape evolution models and application to the anticipation or rapid assessment of secondary hazards associated with earthquakes.
Bedrock landslides, by excavating deep below saprolite-rock interfaces, create conditions for weathering in which all mineral phases in a lithology are initially unweathered within landslide deposits. As a result, the most labile phases dominate the weathering immediately after mobilisation and during a transient period of depletion. This mode of dissolution can strongly alter the overall output of solutes from catchments and their contribution to global chemical cycles if landslide-derived material is retained in catchments for extended periods after mass wasting.
Linking together the processes of rapid physical erosion and the resultant chemical dissolution of rock is a crucial step in building an overall deterministic understanding of weathering in mountain belts. Landslides, which are the most volumetrically important geomorphic process at these high rates of erosion, can generate extremely high rates of very localised weathering. To elucidate how this process works we have taken advantage of uniquely intense landsliding, resulting from Typhoon Morakot, in the T'aimali River and surrounds in southern Taiwan. Combining detailed analysis of landslide seepage chemistry with estimates of catchment-by-catchment landslide volumes, we demonstrate that in this setting the primary role of landslides is to introduce fresh, highly labile mineral phases into the surface weathering environment. There, rapid weathering is driven by the oxidation of pyrite and the resultant sulfuric-acid-driven dissolution of primarily carbonate rock. The total dissolved load correlates well with dissolved sulfate - the chief product of this style of weathering - in both landslides and streams draining the area (R-2 = 0.841 and 0.929 respectively; p < 0.001 in both cases), with solute chemistry in seepage from landslides and catchments affected by significant landsliding governed by the same weathering reactions. The predominance of coupled carbonate-sulfuric-acid-driven weathering is the key difference between these sites and previously studied landslides in New Zealand (Emberson et al., 2016), but in both settings increasing volumes of landslides drive greater overall solute concentrations in streams.
Bedrock landslides, by excavating deep below saprolite-rock interfaces, create conditions for weathering in which all mineral phases in a lithology are initially unweathered within landslide deposits. As a result, the most labile phases dominate the weathering immediately after mobilisation and during a transient period of depletion. This mode of dissolution can strongly alter the overall output of solutes from catchments and their contribution to global chemical cycles if landslide-derived material is retained in catchments for extended periods after mass wasting.
A link between chemical weathering and physical erosion exists at the catchment scale over a wide range of erosion rates(1,2). However, in mountain environments, where erosion rates are highest, weathering may be kinetically limited(3-5) and therefore decoupled from erosion. In active mountain belts, erosion is driven by bedrock landsliding(6) at rates that depend strongly on the occurrence of extreme rainfall or seismicity(7). Although landslides affect only a small proportion of the landscape, bedrock landsliding can promote the collection and slow percolation of surface runoff in highly fragmented rock debris and create favourable conditions for weathering. Here we show from analysis of surface water chemistry in the Southern Alps of New Zealand that weathering in bedrock landslides controls the variability in solute load of these mountain rivers. We find that systematic patterns in surface water chemistry are strongly associated with landslide occurrence at scales from a single hillslope to an entire mountain belt, and that landslides boost weathering rates and river solute loads over decades. We conclude that landslides couple erosion and weathering in fast-eroding uplands and, thus, mountain weathering is a stochastic process that is sensitive to climatic and tectonic controls on mass wasting processes.
Long-term estimates of the dissolution of silicate rock are generally derived from a range of isotopic proxies, such as the radiogenic strontium isotope ratio (Sr-87/Sr-86), which are preserved in sediment archives. For these systems to fairly represent silicate weathering, the changes in isotopic ratios in terrestrial surface waters should correspond to changes in the overall silicate dissolution. This assumes that the silicate mineral phases that act as sources of a given isotope dissolve at a rate that is proportional to the overall silicate weathering. Bedrock landsliding exhumes large quantities of fresh rock for weathering in transient storage, and rapid weathering in these deposits is controlled primarily by dissolution of the most reactive phases. In this study, we test the hypothesis that preferential weathering of these labile minerals can decouple the dissolution of strontium sources from the actual silicate weathering rates in the rapidly eroding Western Southern Alps (WSA) of New Zealand. We find that rapid dissolution of relatively radiogenic calcite and biotite in landslides leads to high local fluxes in strontium with isotopic ratios that offer no clear discrimination between sources. These higher fluxes of radiogenic strontium are in contrast to silicate weathering rates in landslides that are not systematically elevated. On a mountain belt scale, radiogenic strontium fluxes are not coupled to volumes of recent landslides in large (>100 km(2)) catchments, but silicate weathering fluxes are. Such decoupling is likely due first to the broad variability in the strontium content of carbonate minerals, and second to the combination of radiogenic strontium released from both biotite and carbonate in recent landslides. This study supports previous work suggesting the limited utility of strontium isotopes as a system to study silicate weathering in the WSA. Crucially however, in settings where bedrock landsliding is a dominant erosive process there is potential for both random and systematic bias in isotope proxies if the most reactive phases exposed for dissolution by landslides disproportionately contribute to the proxy of choice. This clearly suggests that the isotopic composition of marine Sr is a proxy for periods of rapid mountain uplift and erosion rather than for the associated enhanced silicate weathering. (C) 2017 Elsevier Ltd. All rights reserved.
Prediction of the area affected by earthquake-induced landsliding based on seismological parameters
(2017)
We present an analytical, seismologically consistent expression for the surface area of the region within which most landslides triggered by an earthquake are located (landslide distribution area). This expression is based on scaling laws relating seismic moment, source depth, and focal mechanism with ground shaking and fault rupture length and assumes a globally constant threshold of acceleration for onset of systematic mass wasting. The seismological assumptions are identical to those recently used to propose a seismologically consistent expression for the total volume and area of landslides triggered by an earthquake. To test the accuracy of the model we gathered geophysical information and estimates of the landslide distribution area for 83 earthquakes. To reduce uncertainties and inconsistencies in the estimation of the landslide distribution area, we propose an objective definition based on the shortest distance from the seismic wave emission line containing 95% of the total landslide area. Without any empirical calibration the model explains 56% of the variance in our dataset, and predicts 35 to 49 out of 83 cases within a factor of 2, depending on how we account for uncertainties on the seismic source depth. For most cases with comprehensive landslide inventories we show that our prediction compares well with the smallest region around the fault containing 95% of the total landslide area. Aspects ignored by the model that could explain the residuals include local variations of the threshold of acceleration and processes modulating the surface ground shaking, such as the distribution of seismic energy release on the fault plane, the dynamic stress drop, and rupture directivity. Nevertheless, its simplicity and first-order accuracy suggest that the model can yield plausible and useful estimates of the landslide distribution area in near-real time, with earthquake parameters issued by standard detection routines.
Presentation and Analysis of a Worldwide Database of Earthquake-Induced Landslide Inventories
(2017)
Earthquake-induced landslide (EQIL) inventories are essential tools to extend our knowledge of the relationship between earthquakes and the landslides they can trigger. Regrettably, such inventories are difficult to generate and therefore scarce, and the available ones differ in terms of their quality and level of completeness. Moreover, access to existing EQIL inventories is currently difficult because there is no centralized database. To address these issues, we compiled EQIL inventories from around the globe based on an extensive literature study. The database contains information on 363 landslide-triggering earthquakes and includes 66 digital landslide inventories. To make these data openly available, we created a repository to host the digital inventories that we have permission to redistribute through the U.S. Geological Survey ScienceBase platform. It can grow over time as more authors contribute their inventories. We analyze the distribution of EQIL events by time period and location, more specifically breaking down the distribution by continent, country, and mountain region. Additionally, we analyze frequency distributions of EQIL characteristics, such as the approximate area affected by landslides, total number of landslides, maximum distance from fault rupture zone, and distance from epicenter when the fault plane location is unknown. For the available digital EQIL inventories, we examine the underlying characteristics of landslide size, topographic slope, roughness, local relief, distance to streams, peak ground acceleration, peak ground velocity, and Modified Mercalli Intensity. Also, we present an evaluation system to help users assess the suitability of the available inventories for different types of EQIL studies and model development.
Prediction of the area affected by earthquake-induced landsliding based on seismological parameters
(2017)
We present an analytical, seismologically consistent expression for the surface area of the region within which most landslides triggered by an earthquake are located (landslide distribution area). This expression is based on scaling laws relating seismic moment, source depth, and focal mechanism with ground shaking and fault rupture length and assumes a globally constant threshold of acceleration for onset of systematic mass wasting. The seismological assumptions are identical to those recently used to propose a seismologically consistent expression for the total volume and area of landslides triggered by an earthquake. To test the accuracy of the model we gathered geophysical information and estimates of the landslide distribution area for 83 earthquakes. To reduce uncertainties and inconsistencies in the estimation of the landslide distribution area, we propose an objective definition based on the shortest distance from the seismic wave emission line containing 95% of the total landslide area. Without any empirical calibration the model explains 56% of the variance in our dataset, and predicts 35 to 49 out of 83 cases within a factor of 2, depending on how we account for uncertainties on the seismic source depth. For most cases with comprehensive landslide inventories we show that our prediction compares well with the smallest region around the fault containing 95% of the total landslide area. Aspects ignored by the model that could explain the residuals include local variations of the threshold of acceleration and processes modulating the surface ground shaking, such as the distribution of seismic energy release on the fault plane, the dynamic stress drop, and rupture directivity. Nevertheless, its simplicity and first-order accuracy suggest that the model can yield plausible and useful estimates of the landslide distribution area in near-real time, with earthquake parameters issued by standard detection routines.
Erosion and tectonic uplift are widely thought to be coupled through feedbacks involving orographic precipitation, relief development, and crustal weakening. In many orogenic systems, it can be difficult to distinguish whether true feedbacks exist, or whether observed features are a consequence of tectonic forcing. To help elucidate these interactions, we examine Gongga Shan, a 7556 m peak on the eastern margin of the Tibetan Plateau where cosmogenic Be-10 basin-wide erosion rates reach >5 mm/yr, defining a region of localized rapid erosion associated with a restraining bend in the left-lateral Xianshuihe Fault. Erosion rates are consistent with topography, thermochronometry, and geodetic data, suggesting a stable pattern of uplift and exhumation over at least the past 2-3 My. Transpression along the Xianshuihe Fault, orographically enhanced precipitation, thermally weakened crust, and substantial local relief all developed independently in the Gongga region and existed there prior to the uplift of Gongga Shan. However, only where all of these conditions are present do the observed topographic and erosional extremes exist, and their relative timing indicates that these conditions are not a consequence of rapid uplift. We conclude that their collocation at 3-4 Ma set into motion a series of feedbacks between erosion and uplift that has resulted in the exceptionally high topography and rapid erosion rates observed today.
Mountain channels can be strongly coupled with adjacent hillslopes, exchanging both mass and energy. However, hypotheses of the underlying cause and effect relations are based on indirect observations that do not resolve the mechanics of channel-hillslope coupling at the process scale. Here we present direct observational data of a coupled channel-hillslope system in the catchment area of the Erlenbach, a mountain stream in Switzerland. A slow-moving landslide flanking this alpine stream failed after a flood had eroded an alluvial step in the channel at its base, representing evidence for an upsystem link in channel-hillslope coupling. Progressive accumulation of landslide debris in the channel resulted in a renewed step, stabilizing the hillslope and restoring the channel step in a downsystem link. Thus, upsystem and downsystem coupling mechanisms are joined in a negative feedback cycle. In this cycle, debuttressing and rebuttressing due to channel bed erosion and alluviation are the dominant controls on hillslope stability. Based on an order of magnitude estimate it is plausible that the observed feedback mechanism is a relevant process in the production of coarse (>2 mm) sediment in the Erlenbach.
Bedrock landsliding in mountain belts can elevate overall chemical weathering rates through rapid dissolution of exhumed reactive mineral phases in transiently stored deposits. This link between a key process of erosion and the resultant weathering affects the sequestering of carbon dioxide through weathering of silicate minerals and broader links between erosion in active orogens and climate change. Here we address the effect on the carbon cycle of weathering induced by bedrock landsliding in Taiwan and the Western Southern Alps of New Zealand. Using solute chemistry data from samples of seepage from landslide deposits and river discharge from catchments with variable proportions of landsliding, we model the proportion of silicate and carbonate weathering and the balance of sulfuric and carbonic acids that act as weathering agents. We correct for secondary precipitation, geothermal, and cyclic input, to find a closer approximation of the weathering explicitly occurring within landslide deposits. We find highly variable proportions of sulfuric and carbonic acids driving weathering in landslides and stable hillslopes. Despite this variability, the predominance of rapid carbonate weathering within landslides and catchments where mass wasting is prevalent results at best in limited sequestration of carbon dioxide by this process of rapid erosion. In many cases where sulfuric acid is a key weathering agent, a net release of CO2 to the atmosphere occurs. This suggests that a causal link between erosion in mountain belts and climate change through the sequestration of CO2, if it exists, must operate through a process other than chemical weathering driven by landsliding. Plain Language Summary There is a long-standing debate surrounding the link between erosion and climate. It is often suggested that as temperatures increase, rainier and stormier weather could increase erosion of rock; as that rock is exposed, silicate minerals within could break down, which, on long time scales, can remove CO2 from the atmosphere, lowering global temperatures and acting as a negative feedback. Recent studies have shown that landslide deposits are key locations for the link between chemical weathering and physical erosion in some mountain belts. To test how landslides affect the erosion-climate link, we used samples of water seeping through landslides in Taiwan and New Zealand to calculate the amount of carbon dioxide that is either absorbed or released through this chemical reaction. We find that the large amount of freshly exposed rock in Taiwanese landslide deposits contains significant carbonate rock and sulfide minerals; the net result of the weathering of these minerals is a release of carbon dioxide, which inverts the traditional perspective on the role erosion plays in controlling carbon dioxide release. In some mountain belts, it seems that increased erosion and resulting landsliding may act to increase carbon dioxide in the air, opening further questions into the nature of erosional-climatic links.
Microbial oxidation of lithospheric organic carbon in rapidly eroding tropical mountain soils
(2018)
Lithospheric organic carbon ("petrogenic"; OCpetro) is oxidized during exhumation and subsequent erosion of mountain ranges. This process is a considerable source of carbon dioxide (CO2) to the atmosphere over geologic time scales, but the mechanisms that govern oxidation rates in mountain landscapes are poorly constrained. We demonstrate that, on average, 67 +/- 11% of the OCpetro initially present in bedrock exhumed from the tropical, rapidly eroding Central Range of Taiwan is oxidized in soils, leading to CO2 emissions of 6.1 to 18.6 metric tons of carbon per square kilometer per year. The molecular and isotopic evolution of bulk OC and lipid biomarkers during soil formation reveals that OCpetro remineralization is microbially mediated. Rapid oxidation in mountain soils drives CO2 emission fluxes that increase with erosion rate, thereby counteracting CO2 drawdown by silicate weathering and biospheric OC burial.
Landslide hazard motivates the need for a deeper understanding of the events that occur before, during, and after catastrophic slope failures. Due to the destructive nature of such events, in situ observation is often difficult or impossible. Here, we use data from a network of 58 seismic stations to characterise a large landslide at the Askja caldera, Iceland, on 21 July 2014. High data quality and extensive network coverage allow us to analyse both long- and short-period signals associated with the landslide, and thereby obtain information about its triggering, initiation, timing, and propagation. At long periods, a landslide force history inversion shows that the Askja landslide was a single, large event starting at the SE corner of the caldera lake at 23:24:05 UTC and propagating to the NW in the following 2 min The bulk sliding mass was 7-16 x 10(10) kg, equivalent to a collapsed volume of 35-80 x 10(6) m(3). The sliding mass was displaced downslope by 1260 +/- 250 m. At short periods, a seismic tremor was observed for 30 min before the landslide. The tremor is approximately harmonic with a fundamental frequency of 2.3 Hz and shows time-dependent changes of its frequency content. We attribute the seismic tremor to stick-slip motion along the landslide failure plane. Accelerating motion leading up to the catastrophic slope failure culminated in an aseismic quiescent period for 2 min before the landslide. We propose that precursory seismic signals may be useful in landslide early-warning systems. The 8 h after the main landslide failure are characterised by smaller slope failures originating from the destabilised caldera wall decaying in frequency and magnitude. We introduce the term "afterslides" for this subsequent, declining slope activity after a large landslide.
Himalayan rivers are frequently hit by catastrophic floods that are caused by the failure of glacial lake and landslide dams; however, the dynamics and long-term impacts of such floods remain poorly understood. We present a comprehensive set of observations that capture the July 2016 glacial lake outburst flood (GLOF) in the Bhotekoshi/Sunkoshi River of Nepal. Seismic records of the flood provide new insights into GLOF mechanics and their ability to mobilize large boulders that otherwise prevent channel erosion. Because of this boulder mobilization, GLOF impacts far exceed those of the annual summer monsoon, and GLOFs may dominate fluvial erosion and channel-hillslope coupling many tens of kilometers downstream of glaciated areas. Long-term valley evolution in these regions may therefore be driven by GLOF frequency and magnitude, rather than by precipitation.
Soil degradation is a severe and growing threat to ecosystem services globally. Soil loss is often nonlinear, involving a rapid deterioration from a stable eco-geomorphic state once a tipping point is reached. Soil loss thresholds have been studied at plot scale, but for landscapes, quantitative constraints on the necessary and sufficient conditions for tipping points are rare. Here, we document a landscape-wide eco-geomorphic tipping point at the edge of the Tibetan Plateau and quantify its drivers and erosional consequences. We show that in the upper Kali Gandaki valley, Nepal, soil formation prevailed under wetter conditions during much of the Holocene. Our data suggest that after a period of human pressure and declining vegetation cover, a 20% reduction of relative humidity and precipitation below 200 mm/year halted soil formation after 1.6 ka and promoted widespread gullying and rapid soil loss, with irreversible consequences for ecosystem services.