skip to main content


Title: A new approach to mapping landslide hazards: a probabilistic integration of empirical and physically based models in the North Cascades of Washington, USA
Abstract. We developed a new approach for mapping landslide hazards by combiningprobabilities of landslide impacts derived from a data-driven statisticalapproach and a physically based model of shallow landsliding. Ourstatistical approach integrates the influence of seven site attributes (SAs) onobserved landslides using a frequency ratio (FR) method. Influential attributesand resulting susceptibility maps depend on the observations of landslidesconsidered: all types of landslides, debris avalanches only, or source areasof debris avalanches. These observational datasets reflect the detection ofdifferent landslide processes or components, which relate to differentlandslide-inducing factors. For each landslide dataset, a stability index (SI) is calculated as a multiplicative result of the frequency ratios for all attributes and is mapped across our study domain in the North Cascades National Park Complex (NOCA), Washington, USA. A continuous function is developed to relate local SI values to landslide probability based on a ratio of landslide and non-landslide grid cells. The empirical model probability derived from the debris avalanche source area dataset is combined probabilistically with a previously developed physically based probabilistic model. A two-dimensional binning method employs empirical andphysically based probabilities as indices and calculates a joint probabilityof landsliding at the intersections of probability bins. A ratio of thejoint probability and the physically based model bin probability is used asa weight to adjust the original physically based probability at each gridcell given empirical evidence. The resulting integrated probability oflandslide initiation hazard includes mechanisms not captured by the infinite-slope stability model alone. Improvements in distinguishing potentiallyunstable areas with the proposed integrated model are statisticallyquantified. We provide multiple landslide hazard maps that land managers canuse for planning and decision-making, as well as for educating the publicabout hazards from landslides in this remote high-relief terrain.  more » « less
Award ID(s):
1663859
NSF-PAR ID:
10185372
Author(s) / Creator(s):
; ;
Date Published:
Journal Name:
Natural Hazards and Earth System Sciences
Volume:
19
Issue:
11
ISSN:
1684-9981
Page Range / eLocation ID:
2477 to 2495
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. null (Ed.)
    Abstract. Landslides are the main source of sediment in most mountain ranges. Rivers then act as conveyor belts, evacuating landslide-derived sediment. Sediment dynamics are known to influence landscape evolution through interactions among landslide sediment delivery, fluvial transport and river incision into bedrock. Sediment delivery and its interaction with river incision therefore control the pace of landscape evolution and mediate relationships among tectonics, climate and erosion. Numerical landscape evolution models (LEMs) are well suited to study the interactions among these surface processes. They enable evaluation of a range of hypotheses at varying temporal and spatial scales. While many models have been used to study the dynamic interplay between tectonics, erosion and climate, the role of interactions between landslide-derived sediment and river incision has received much less attention. Here, we present HyLands, a hybrid landscape evolution model integrated within the TopoToolbox Landscape Evolution Model (TTLEM) framework. The hybrid nature of the model lies in its capacity to simulate both erosion and deposition at any place in the landscape due to fluvial bedrock incision, sediment transport, and rapid, stochastic mass wasting through landsliding. Fluvial sediment transport and bedrock incision are calculated using the recently developed Stream Power with Alluvium Conservation and Entrainment (SPACE) model. Therefore, rivers can dynamically transition from detachment-limited to transport-limited and from bedrock to bedrock–alluvial to fully alluviated states. Erosion and sediment production by landsliding are calculated using a Mohr–Coulomb stability analysis, while landslide-derived sediment is routed and deposited using a multiple-flow-direction, nonlinear deposition method. We describe and evaluate the HyLands 1.0 model using analytical solutions and observations. We first illustrate the functionality of HyLands to capture river dynamics ranging from detachment-limited to transport-limited conditions. Second, we apply the model to a portion of the Namche Barwa massif in eastern Tibet and compare simulated and observed landslide magnitude–frequency and area–volume scaling relationships. Finally, we illustrate the relevance of explicitly simulating landsliding and sediment dynamics over longer timescales for landscape evolution in general and river dynamics in particular. With HyLands we provide a new tool to understand both the long- and short-term coupling between stochastic hillslope processes, river incision and source-to-sink sediment dynamics. 
    more » « less
  2. Abstract

    Bedrock landsliding, including the formation of landslide dams, is a predominant geomorphic process in steep landscapes. Clarifying the importance of hydrologic and seismic mechanisms for triggering deep‐seated landslides remains an ongoing effort, and formulation of geomorphic metrics that predict dam preservation is crucial for quantifying secondary landslide hazards. Here, we identify >200 landslide‐dammed lakes in western Oregon and utilize dendrochronology and enhanced14C dating (“wiggle matching”) of “ghost forests” to establish slope failure timing at 20 sites. Our dated landslide dataset reveals bedrock landsliding has been common since the last Cascadia Subduction Zone earthquake in January 1700 AD. Our study does not reveal landslides that date to 1700 AD. Rather, we observe temporal clustering ofat leastfour landslides in the winter of 1889/1890 AD, coincident with a series of atmospheric rivers that generated one of the largest regionally recorded floods. We use topographic and field analyses to assess the relation between dam preservation and topographic characteristics of the impounded valleys. In contrast to previous studies, we do not observe systematic scaling between dam size and upstream drainage area, though dam stability indices for our sites correspond with “stable” dams elsewhere. Notably, we observe that dams are preferentially preserved at drainage areas of ∼1.5 to 13 km2and valley widths of ∼25 to 80 m, which may reflect the reduced downstream influence of debris flows and the accumulation of mature conifer trees upstream from landslide‐dammed lake outlets. We suggest that wood accumulation upstream of landslide dams tempers large stream discharges, thus inhibiting dam incision.

     
    more » « less
  3. null (Ed.)
    Abstract Many existing models that predict landslide hazards utilize ground-based sources of precipitation data. In locations where ground-based precipitation observations are limited (i.e., a vast majority of the globe), or for landslide hazard models that assess regional or global domains, satellite multisensor precipitation products offer a promising near-real-time alternative to ground-based data. NASA’s global Landslide Hazard Assessment for Situational Awareness (LHASA) model uses the Integrated Multisatellite Retrievals for Global Precipitation Measurement (IMERG) product to issue hazard “nowcasts” in near–real time for areas that are currently at risk for landsliding. Satellite-based precipitation estimates, however, can contain considerable systematic bias and random error, especially over mountainous terrain and during extreme rainfall events. This study combines a precipitation error modeling framework with a probabilistic adaptation of LHASA. Compared with the routine version of LHASA, this probabilistic version correctly predicts more of the observed landslides in the study region with fewer false alarms by high hazard nowcasts. This study demonstrates that improvements in landslide hazard prediction can be achieved regardless of whether the IMERG error model is trained using abundant ground-based precipitation observations or using far fewer and more scattered observations, suggesting that the approach is viable in data-limited regions. Results emphasize the importance of accounting for both random error and systematic satellite precipitation bias. The approach provides an example of how environmental prediction models can incorporate satellite precipitation uncertainty. Other applications such as flood and drought monitoring and forecasting could likely benefit from consideration of precipitation uncertainty. 
    more » « less
  4. Abstract

    The effects of strong ground shaking on hillslope stability can persist for many years after a large earthquake, leading to an increase in the rates of post earthquake land sliding. The factors that control the rate of post‐earthquake land sliding are poorly constrained, hindering our ability to reliably forecast how landscapes and landslide hazards and risk evolve. To address this, we use a unique data set comprising high‐resolution terrestrial laser scans and airborne lidar captured during and after the 2010–2011 Canterbury Earthquake Sequence, New Zealand. This earthquake sequence triggered thousands of rock falls, and rock and debris avalanches (collectively referred to as “rockfall”), resulting in loss‐of‐life and damage to residential dwellings, commercial buildings and other infrastructure in the Port Hills of Christchurch, New Zealand. This unique data set spans 5 years and includes five significant earthquakes. We used these data to (a) quantify the regional‐scale “rockfall” rates in response to these earthquakes and the postearthquake decay in rockfall rates with time; and (b) investigate the site‐specific factors controlling the location of seismic and nonseismic rockfalls using frequency ratios and logistic regression techniques. We found that rockfall rates increased significantly in response to the initial earthquake that generated the strongest shaking in the sequence—The MW6.2 22 February 2011 event—Compared to the long‐term background rates derived from the dating of pre‐2010 talus piles at the toe of the slopes. Non seismic rockfall rates also increased immediately after the 22 February 2011 earthquake and decayed with time following a power‐law trend. About 50% of the decay back to the pre‐earthquake rockfall rates occurred within 1–5 years after the 22 February 2011 earthquake. Our results show that the short‐term decay in rockfall rates over time, after the initial earthquake, was attributed to the subsequent erosion of seismically damaged rock mass materials caused by environmental processes such as rain. For earthquake‐induced rockfall at the regional‐scale, the peak ground accelerations is the most significant variable in forecasting rockfall volume, followed by the relative height above the base of the slope. For both earthquake and non‐seismic conditions at the site‐specific scale, the probability of rockfall increases when the adjacent areas have failed previously, indicating that accrued damage preconditions localized areas of the slope for subsequent failure. Such preconditioning is a crucial factor driving subsequent rockfalls; that is, future rockfalls are likely to cluster near areas that failed in the past.

     
    more » « less
  5. Abstract

    Bedrock landslides shape topography and mobilize large volumes of sediment. Yet, interactions between landslide‐produced sediment and fluvial systems that together govern large‐scale landscape evolution are not well understood. To explain morphological patterns observed in steep, landslide‐prone terrain, we explicitly model stochastic landsliding and associated sediment dynamics. The model accounts for several common landscape features such as slope frequency distributions, which include values in excess of regional stability limits, quasi‐planar hillslopes decorated with straight, closely spaced channel‐like features, and accumulation of sediment in valley networks rather than on hillslopes. Stochastic landsliding strongly affects the magnitude and timing of sediment supply to the fluvial system. We show that intermittent sediment supply is ultimately reflected in topography. At dynamic equilibrium, landslide‐derived sediment pulses generate persistent landscape dynamism through the formation and breaching of landslide dams and epigenetic gorges as landslides force shifts in channel positions. Our work highlights the importance of interactions between landslides and sediment dynamics that ultimately control landscape‐scale response to environmental change.

     
    more » « less