Refine
Document Type
- Article (13)
- Postprint (4)
- Doctoral Thesis (1)
Language
- English (18)
Is part of the Bibliography
- yes (18) (remove)
Keywords
- 2008 Wenchuan earthquake (2)
- erosion (2)
- model (2)
- mountain belt (2)
- 1988 Saguenay (1)
- 1994 Northridge (1)
- 2010 Haiti earthquake (1)
- Erdbeben (1)
- Erdrutsch (1)
- Erosion (1)
Institute
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.
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.
The topographic signature of a mountain belt depends on the interplay of tectonic, climatic and erosional processes, whose relative importance changes over times, while quantifying these processes and their rates at specific times remains a challenge. The eastern Andes of central Bolivia offer a natural laboratory in which such interplay has been debated. Here, we investigate the Rio Grande catchment which crosses orthogonally the eastern Andes orogen from the Eastern Cordillera into the Subandean Zone, exhibiting a catchment relief of up to 5000 m. Despite an enhanced tectonic activity in the Subandes, local relief, mean and modal slopes and channel steepness indices are largely similar compared to the Eastern Cordillera and the intervening Interandean Zone. Nevertheless, a dataset of 57 new cosmogenic 10Be and 26AI catchment wide denudation rates from the Rio Grande catchment reveals up to one order of magnitude higher denudation rates in the Subandean Zone (mean 0.8 mm/yr) compared to the upstream physiographic regions. We infer that tectonic activity in the thrusting dominated Subandean belt causes higher denudation rates based on cumulative rock uplift investigations and due to the absence of a pronounced climate gradient. Furthermore, the lower rock strength of the Subandean sedimentary units correlates with mean slopes similar to the ones of the Eastern Cordillera and Interandean Zone, highlighting the fact, that lithology and rock strength can control high denudation rates at low slopes.
Low denudation rates measured at the outlet of the Rio Grande catchment (Abapo) are interpreted to be a result of a biased cosmogenic nuclide mixing that is dominated by headwater signals from the Eastern Cordillera and the Interandean zone and limited catchment sediment connectivity in the lower river reaches. Therefore, comparisons of short- (i.e., sediment yield) and millennial denudation rates require caution when postulating tectonic and/or climatic forcing without detailed studies. (C) 2015 The Authors. Published by Elsevier B.V.
Shallow earthquakes frequently disturb the hydrological and mechanical state of the subsurface, with consequences for hazard and water management. Transient post-seismic hydrological behavior has been widely reported, suggesting that the recovery of material properties (relaxation) following ground shaking may impact groundwater fluctuations. However, the monitoring of seismic velocity variations associated with earthquake damage and hydrological variations are often done assuming that both effects are independent. In a field site prone to highly variable hydrological conditions, we disentangle the different forcing of the relative seismic velocity variations delta v retrieved from a small dense seismic array in Nepal in the aftermath of the 2015 M-w 7.8 Gorkha earthquake. We successfully model transient damage effects by introducing a universal relaxation function that contains a unique maximum relaxation timescale for the main shock and the aftershocks, independent of the ground shaking levels. Next, we remove the modeled velocity from the raw data and test whether the corresponding residuals agree with a background hydrological behavior we inferred from a previously calibrated groundwater model. The fitting of the delta v data with this model is improved when we introduce transient hydrological properties in the phase immediately following the main shock. This transient behavior, interpreted as an enhanced permeability in the shallow subsurface, lasts for similar to 6 months and is shorter than the damage relaxation (similar to 1 yr). Thus, we demonstrate the capability of seismic interferometry to deconvolve transient hydrological properties after earthquakes from non-linear mechanical recovery.
The large, shallow earthquakes at Northridge, California (1994), Chi-Chi, Taiwan (1999), and Wenchuan, China (2008), each triggered thousands of landslides. We have determined the position of these landslides along hillslopes, normalizing for statistical bias. The landslide patterns have a co-seismic signature, with clustering at ridge crests and slope toes. A cross-check against rainfall-induced landslide inventories seems to confirm that crest clustering is specific to seismic triggering as observed in previous studies. In our three study areas, the seismic ground motion parameters and lithologic and topographic features used do not seem to exert a primary control on the observed patterns of landslide clustering. However, we show that at the scale of the epicentral area, crest and toe clustering occur in areas with specific geological features. Toe clustering of seismically induced landslides tends to occur along regional major faults. Crest clustering is concentrated at sites where the lithology along hillslopes is approximately uniform, or made of alternating soft and hard strata, and without strong overprint of geological structures. Although earthquake-induced landslides locate higher on hillslopes in a statistically significant way, geological features strongly modulate the landslide position along the hillslopes. As a result the observation of landslide clustering on topographic ridges cannot be used as a definite indicator of the topographic amplification of ground shaking.
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.
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.
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.
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.