A Novel Approach To Probabilistic Seismic Landslide Hazard Mapping Using
A Novel Approach To Probabilistic Seismic Landslide Hazard Mapping Using
A Novel Approach To Probabilistic Seismic Landslide Hazard Mapping Using
Engineering Geology
journal homepage: www.elsevier.com/locate/enggeo
A R T I C L E I N F O A B S T R A C T
Keywords: As a potential destructive secondary hazard initiated by earthquakes, earthquake-triggered landslides can cause
Probabilistic analysis significant losses. Seismic landslide hazard maps are important for land use planning and landslide hazard
Seismic landslide reduction in tectonically active areas. A comprehensive probabilistic seismic landslide hazard mapping approach
Hazard map
incorporating the uncertainties of slope properties and displacement prediction is required to promote the
Monte Carlo simulation
Critical slope angle
application of hazard maps. Based on the definition of the critical slope angle, this paper proposes an alternative
seismic landslide hazard category considering the uncertainty of geotechnical properties of soil using Monte
Carlo simulations and incorporating the ground motion variability based on the displacement hazard curve. The
proposed new hazard category is applied to the region of Anchorage, Alaska, for the seismic landslide hazard
mapping process. Subsequently, parameter analyses are performed to explore the sensitivity among the factors
that contributed to the seismic landslide hazard map. The results of seismic landslide hazard maps show that the
area of hazardous zone of most geologic units in present study is smaller than that in the deterministic analysis.
Comparing with the probabilistic method using logic tree analysis, the proposed approach using Monte Carlo
simulations predicts different landslide hazardous zone due to the probabilistic distribution of geotechnical
properties of soil. The predicted hazardous zone is also sensitive to the resolution of the digital elevation model
and different coefficient of variation combinations.
* Corresponding author.
E-mail addresses: [email protected] (C. Li), [email protected] (G. Wang), [email protected] (J. He), [email protected] (Y. Wang).
https://doi.org/10.1016/j.enggeo.2022.106616
Received 27 September 2021; Received in revised form 17 February 2022; Accepted 8 March 2022
Available online 15 March 2022
0013-7952/© 2022 Elsevier B.V. All rights reserved.
C. Li et al. Engineering Geology 301 (2022) 106616
2017). However, these mapping methods of seismic landslide hazards analytical approach should be sufficient to consider other causative
are still considered deterministic approaches because a single empirical factors.
displacement model and a set of specific input parameters associated This study proposes an integral probabilistic method to promote the
with a slope are adopted in the calculation. That is, the epistemic un application of seismic landslide hazard maps that considers the
certainties of the slope properties and the aleatory variability of seismic epistemic uncertainties associated with slope properties using Monte
motions were not considered sufficiently. The influence of different Carlo simulation. The definition of the critical slope angle is described
displacement models and input parameters on the comparison between for a convenient application of region mapping incorporating the un
observed and predicted landslides is significant (Dreyfus et al., 2013). A certainties of different empirical displacement models and the aleatory
robust method incorporating the uncertainty in strength parameters variability of seismic motions. Instead of using sliding displacement
should be developed. threshold, a seismic landslide hazard category defined by critical slope
Capolongo et al. (2002) and Refice and Capolongo (2002) applied angle was proposed in this study. Subsequently, such categories were
physically based models to probabilistic seismic landslide hazard applied to the regional mapping process of Anchorage, Alaska, for
assessment, considering the uncertainties of input parameters. Wang validation and comparison with previous works. The effects of proba
et al. (2008) introduced a method to estimate the probabilities of bility distribution and digital elevation model (DEM) resolution are
displacement exceeding 10 cm, which were shown as a map of seismic discussed. Finally, a parametric analysis was conducted for the coeffi
landslide hazards using Monte Carlo simulation. Saygili and Rathje cient of variation (COV) and unit weight of the soil.
(2009) proposed a probabilistic approach using the concept of
displacement hazard curves to consider both the variability of ground 2. Probability-based seismic landslide hazard mapping
motion and the uncertainty in the prediction of sliding displacement. procedures
Wang and Rathje (2015) extended the works of Saygili and Rathje
(2009) to consider the uncertainties within slope properties and In this section, the construction of a seismic landslide hazard map
empirical displacement models using logic-tree analysis. The application ping procedure is elaborated by comprehensively expressing the un
of probability-based seismic landslide hazard mapping procedures has certainties of relevant factors. In regional analysis, the rigid sliding
received increasing attention (Liu et al., 2016; Martino et al., 2018; model based on the assumption of an infinite slope model has been
Wang et al., 2020). As indicated by Wasowski et al. (2011), the widely used in earthquake-triggered landslide hazard assessment and
improvement of the regional-scale assessment of seismic landslide sus seismic landslide hazard maps. This method was implemented in the
ceptibility and hazard represents another top research hotspot owing to present analysis. Following the general mapping logic of existing liter
their potential in guiding land use planning and the development of ature, this study includes three main steps for a region. A seismic
major lifelines. Hence, a thorough probability-based method is neces landslide hazard category was established after storing the DEM map
sary for seismic landslide hazard mapping. Meanwhile, the computation along with the topographic map in the ArcGIS platform, and then the
cost is best to be limited to a suitable range, and the extensibility of the mapping procedure was performed by combining the slope map and the
Fig. 1. Flowchart for seismic landslide hazard mapping procedures in the present study.
2
C. Li et al. Engineering Geology 301 (2022) 106616
hazard category. As shown in Fig. 1, the flowchart of the application of hypercube sampling (LHS) technique produced a realistic distribution
the proposed probability-based procedures for seismic landslide hazard from a relatively limited number of samples and allowed less loss of
zonation illustrates the main routine of the three steps. Among them, the accuracy, which led to efficient computation in the Monte Carlo simu
seismic landslide hazard category based on the critical slope angle is the lation process (Refice and Capolongo, 2002). In this study, the slope
most important work in this study. It consists of three sections including angle α is taken as a controlling variable, ranging from 1◦ to 90◦ with 1◦
the yield acceleration statistical samples, weighted mean values of the increments, and then each input parameter with assumed probability
annual rate of exceedance, and critical slope angle distributions. Each distribution is sampled using the Monte Carlo simulation together with
part of this section is detailed as follows. the LHS technique in the same geologic unit. These statistical samples,
along with a controlled slope angle, are substituted into eqs. (1) and (2),
2.1. Monte Carlo simulation of yield acceleration and the majority of the yield acceleration values that form the yield
acceleration samples can be calculated simultaneously. Many groups of
This is the first step in the seismic landslide hazard category. The ky samples exist, each of which is related to a specific geologic unit and a
epistemic uncertainty of input parameters (i.e., geotechnical properties selected slope angle. Thus, the repeat procedure leads to many ky values
of soil) is associated with the process of yield acceleration calculation. for all slope angle ranges and for the entire region composed of some
Based on the concept of permanent displacement, which embodies the geologic units.
characterizations of earth slopes subjected to seismic environments; the
susceptibility of slopes under earthquake loading can be characterized 2.2. Weighted mean annual rate of exceedance
by an index of yield acceleration, as stated by Newmark (Newmark,
1965). The yield acceleration (ky) can be explicitly described as The displacement hazard curve has been demonstrated as an effec
tive method to combine the displacement threshold of the described
ky = (FSstatic − 1)⋅sinα⋅g (1)
seismic slope performance with different seismic hazard levels, based on
the yield acceleration. It describes the seismic hazard level (mean annual
where g is the acceleration owing to gravity, α is the slope angle, and
rate of exceedance λD) for different displacement values x and was
FSstatic represents the static factor of safety of the slope deduced from an
considered as a fully probabilistic approach in permanent sliding
infinite slope model. The static factor of safety can be expressed as
displacement quantification with a specific hazard level. Wang and
′ ( )
Rathje (2015) detailed a complete description of the displacement
′
c tanφ γ
FSstatic = + 1− m w (2)
γtsinα tanα γ hazard curve (Section 2.4 associated with Eqs. 3–4 of Wang and Rathje
(2015)), and the relevant concept was applied in this study. As shown in
where c’ is the effective cohesion (kPa), φ’ is the effective friction angle, Fig. 2, one ky value corresponds to one displacement hazard curve,
γ is the unit weight of soil (kN/m3), γw is the unit weight of water (kN/ where λD(x) represents the mean annual rate of exceedance for a
m3), t is the sliding block thickness normal to the failure surface or displacement level of x cm. λD(x) increases with an increase in the yield
sliding surface in meters, and m is the proportion of the block thickness acceleration value given a selected displacement value (e.g., 5 cm in
saturated (saturation ratio). Eq. (2) implies that three elements Fig. 2 and λD(5) increases as ky increased from 0.1 g to 0.2 g). Un
contributed to the overall stability analysis, which contains the effects of doubtedly, numerous yield acceleration samples from the Monte Carlo
cohesive strength (first item of the right of the equation), friction angle simulation above deduce lots of λD(x) values with a selected displace
(second item of the right of the equation), and pore water pressure (the ment x under consideration within the same empirical displacement
last item of the right of the equation). model (e.g., RS09 model in Fig. 2).
A specific value of the yield acceleration of a slope can be calculated Owing to the importance of the displacement threshold (e.g., 1 cm, 5
from Eqs. 1–2 under a deterministic group of relevant parameters. cm, and 15 cm) for seismic landslide hazard zonation, a convenient
Therefore, the index of yield acceleration indicates the intrinsic char method that compromises the accuracy of dealing with numerous yield
acteristics of the slope that resist the effect of an earthquake. However, acceleration values and the annual rate of exceedance accompanied is to
uncertainties in slope properties play an important role in the quanti directly establish a regression relationship between ky and λD(x) for a
fication of seismic landslide hazard assessment. Ignoring it and selecting given threshold x. The regression relationship was detailed by Wang
a group of specific parameters would lead to a conservative or non- (2014) and Wang and Rathje (2015). In general, as shown in Fig. 2, after
conservative result for yield acceleration and landslide hazard assess
ment. To consider the epistemic uncertainties associated with these
parameters (c’, φ’, γ, t, m), Wang and Rathje (2015) applied a logic-tree
analysis method in which the branches represent possible selections of
input parameters, and each node is a specific parameter with an assigned
weight value. Unlike logic-tree treatment, which selects discrete limit
values, Monte Carlo simulation is an efficient approach to address un
certain distributed problems (Wang et al., 2008).
The proposed method adopts the Monte Carlo simulation method,
considering that each particular parameter is assumed to obey a
reasonable probability distribution. For example, Fenton and Griffiths
(2003 and 2004) promote a lognormal distribution suitable for geologic
parameters and cohesion while assuming an independent normal dis
tribution for the friction angle with a truncated range. Cherubini (2000)
recommended a lognormal distribution suitable for simulation and
advise that paid more attention to the connection between c’ and φ’. The
unit weight of soil (γ) can be considered to follow a normal distribution.
A similar probability distribution is suggested by assigning the geo
metric parameters of the slope (t and m), which adopt a uniform dis
tribution. Any suitable distribution with relevant statistical parameters
can be used with no restriction imposed on the selection of the assumed Fig. 2. Displacement hazard curves for ky values between 0.1 and 0.2 g
probability distribution for the input parameters. In addition, the Latin (adapted from Wang, 2014).
3
C. Li et al. Engineering Geology 301 (2022) 106616
As expressed in Eq. 4, within the same geologic unit (k) along with a λD (x) equaled the target hazard level λ*. Alternatively, the critical slope
specific displacement threshold (e.g., x = 5 cm), the value λD (x; α, k) angle is taken as the closest abscissa if only the target hazard level de
represents a type of correspondence between slope angle(α) ranging viates from the discrete value of λD (x). Taking the geologic unit G-1 in
from 1◦ to 90◦ and the weighted mean annual rate of exceedance. That Fig. 3 as an example, α*5 cm equals 9◦ and 21◦ according to two types of
is, one slope angle (α) corresponded to one λD (x) provided the selected hazard levels in the threshold x = 5 cm, whereas the critical
slope angle α*15 cm approaches 22◦ and 30◦ in threshold x = 15 cm,
displacement threshold x is determined. Once the value λD (x; α, k)
respectively. Therefore, the critical slope angle (α*x), corresponding to
equals a target hazard level λ*(e.g., λ* = 0.000404 1/year or λ* =
different displacement thresholds (x), for all geological units of a region
0.0021 1/year represent 2% probability and 10% probability of ex
can be obtained under the selected hazard level λ* with the following
ceedance in 50 years, respectively), the corresponding slope angle is
analysis procedures. The magnitude of the critical slope angle represents
defined as the critical slope angle (α*). Figs. 3 (a) and (b) show that
the levels of the seismic landslide hazard and works the same function
λD (x; α, k) varied with an increase in slope angle (α) for three illustrative
equivalent to the displacement threshold under a conditional hazard
geologic units (G-1, G-2, and G-3) under two different displacement
level. This suggests that the aforementioned critical slope angle can be
thresholds (i.e., x = 5 cm and x = 15 cm). The critical slope angle is
used for seismic landslide hazard mapping applications, similar to the
regarded as the abscissa value, which corresponds to the intersection
performance of the displacement threshold.
point of the horizontal target hazard level line and geologic unit λD
curve, considering the situation in which the slope angle varies with
discrete increments (increment set as 1◦ ). There are two situations 2.4. Seismic landslide hazard category
regarding the determination of the critical slope angle from the inter
section point. First, the critical slope angle is the abscissa value when In regional analysis or seismic landslide hazard mapping, instead of
4
C. Li et al. Engineering Geology 301 (2022) 106616
exact predictions of sliding displacement, landslide hazard categories proposed procedures is the derivation of equivalent seismic landslide
are more applicable. When introducing critical slope angles corre hazard categories and how the hazard category works in this region.
sponding to different displacement thresholds and to a selected geologic
unit, a modified assessment criterion for seismic landslide hazards can 3.2. Yield acceleration simulation
be obtained. The proposed seismic landslide hazard categories can be
applied to create seismic landslide hazard maps in the process of land As stated in Section 2, the first preparation is the samples of yield
planning and preliminary evaluation of potential seismic slope failures. acceleration that require a series of input parameters. Seventeen
As shown in Table 1, the middle column represented a common seismic geologic units of this study area along with the assigned shear strength
landslide hazard category used by USGS which is illustrated by parameters are listed in Table 2, whose values are considered as the
displacement thresholds (x = 1 cm, 5 cm, and 15 cm) and the right expected values. The four geologic units (bc, c-bl, l, and s) only provide
column is the equivalent seismic landslide hazard evaluation criterion the value of cohesion; therefore, these parameters are regarded as un
which is established from the concept of critical slope angle (α*1 cm, α*5 drained shear strength, whereas other groups of values with cohesion
cm, and α*15 cm) of the present study in the light of previous sections. and friction angles are considered as effective strength parameters. The
Each location can be assigned to a seismic landslide hazard zone by selection of the COV is consistent with the previous work of Wang and
combining the critical slope angle with the discrete values of the slope Rathje (2015). For example, 30% is selected for undrained shear
map developed from a DEM. For instance, a realistic slope angle of a strength (COV-undrained is 0.3 for bc, c-bl, l, and s), 10% was chosen as the
specified site located between α*1 cm to α*5 cm is defined as a moderate effective friction angle (COV-friction was 0.1), and 20% was assigned as
hazard, meaning the co-seismic sliding displacement of this grid cell is effective cohesion (COV-cohesion was 0.2). The mean value of the strength
located between 1 and 5 cm. The entire area landslide hazards under a parameters along with the selected COV above comprises a basic group
specific earthquake hazard level can be determined by comparing the of statistical parameters, which is known in advance in the present
real slope angle of each grid cell of geologic units with the corresponding study. Subsequently, universal probability distribution forms are
critical slope angle. assumed for strength parameters based on a consistent group of ex
pected values and corresponding COVs. In this study, strength parame
3. Application of the proposed approach in a study region of ters (c’, φ’) are considered to obey lognormal distributions (Fenton and
Anchorage, Alaska Griffiths, 2003 and 2004). Meanwhile, an independent normal distri
bution along with a joint normal distribution with an assumed correla
3.1. Relevant works tion coefficient (ρ = − 0.5) (Alamanis and Dakoulas, 2019) are assigned
to the strength parameters for a better sensitivity comparison. A more
Anchorage, Alaska, was selected as the case study for applicability rigorous probability distribution should be used in the simulation of
verification based on the availability of the required data and the geologic parameters, which is worthy of further study and is beyond the
accessibility of databases of soil properties in such areas. As a repre scope of this study. Based on the features of geologic units consisting of
sentative application, Jibson and Michael (2009) produced seismic different types of soils and on the relevant statistical results of shallow
landslide hazard maps for this region based on a deterministic approach soils of the region by Wang (2014), the sliding block thickness (t) is
that used a series of typical parameter values including shear strength,
sliding block thickness, and saturation ratio to compute yield accelera
Table 2
tion and introduced a displacement prediction model for the calculation
Geologic units and strength parameters (adapted from Wang and Rathje, 2015).
of Newmark displacement. Thus, the corresponding seismic landslide
Unit Cohesion Friction Unit description
hazard of a site was determined by comparing the calculated displace
angle
ment with the displacement thresholds. Wang and Rathje (2015) pro
posed a probabilistic procedure for the mapping of co-seismic landslide af 24 36 Deposits in alluvial fans, alluvial cones, and
emerged deltas
hazards in this region, improving the limitations of the deterministic
al 19 36 Alluvium in abandoned stream channels and in
method above. The aleatory variability in the sliding displacement terraces along modern streams
prediction and the epistemic uncertainty of input parameters and an 21 36 Coarse-grained surficial deposits
different displacement empirical models were considered in their map b 192 40 Bedrock
bc 120 0 Bootlegger Cove clay
ping process with the help of logic tree analysis, and a mean λD threshold
c-br 38 38 Colluvium derived from bedrock on slopes of the
approach was applied for seismic landslide hazard zonation in the study Chugach Mountains
area. c-bl 38 0 Colluvium derived from glacial materials along
Considered as an improvement to previous studies, this study is coastal bluffs
applied to the same study area for a better comparison. The basic f 48 34 Human made fills
ga 38 32 Glacial outwash in irregularly shaped hills
datasets that contributed to seismic landslide hazard mapping were
(including kames, eskers, and kame terraces)
mainly topographic maps, DEM of the study area, slope map derived gm 48 38 Glacial and (or) marine deposits, typically in
from DEM, and geologic unit map connected with input parameters to elongated hills
calculate yield acceleration. A relatively detailed description of such l 144 0 Lake and pond deposits
ls 24 30 Landslide deposits, similar to a unit
basic datasets was accessed from previous studies. The central work
m 43 38 Morainal deposits, generally in long ridges
contributing to the seismic landslide hazard mapping application of the marking the margins of former glaciers
mg 38 37 Marine, glacial, and (or) lacustrine deposits
s 72 0 Silt
Table 1 sh 24 34 Sand deposits in broad, low hills, and windblown
sand deposits in cliff-head dunes near Point
Seismic landslide hazard categories.
Campbell
Hazard category Sliding displacement (D) Critical slope angle (α*) sl 19 34 Sand deposits in a wide low-lying belt around
(used by USGS) (present approach) Connors Lake
Low 0–1 cm 0-α*1cm (a) The first column represents the abbreviation of the geologic unit description.
Moderate 1–5 cm α*1cm –α*5cm (b) The specific value is assumed to be the expectation in probability distribu
High 5–15 cm α*5cm –α*15cm
tions.
Very high >15 cm >α15cm
(c) Cohesion in the unit of Kilopascal (kPa).
Critical slope angle in the unit of degree. (d) Friction angle in the unit of degree.
5
C. Li et al. Engineering Geology 301 (2022) 106616
6
C. Li et al. Engineering Geology 301 (2022) 106616
Table 5
Regression coefficient and corresponding weights for different models (adapted from Wang, 2014).
Displacement threshold Empirical model a1 a2 a3 a4 a5 wi
Table 6
Critical slopes for all geologic units with selected hazard level (x = 5 cm).
Distribution pattern Hazard level af al an b bc c-br c-bl f ga gm l ls m mg s sh sl
Normal 2% 19 18 19 63 9 25 4 25 21 28 13 16 27 25 3 18 16
10% 24 22 24 76 20 30 12 30 26 33 26 20 31 30 11 23 21
Lognormal 2% 19 18 20 69 20 25 12 26 22 29 25 16 27 25 10 18 16
10% 24 22 24 80 27 30 17 31 27 33 33 20 32 29 14 23 21
Joint Normal 2% 20 19 20 66 9 27 4 27 23 30 13 16 28 26 3 19 17
10% 25 23 25 78 21 31 12 32 27 35 26 21 33 30 11 24 22
The symbol 2% represents a 2% probability of exceedance in 50 years and 10% represents a 10% probability of exceedance in 50 years.
7
C. Li et al. Engineering Geology 301 (2022) 106616
8
C. Li et al. Engineering Geology 301 (2022) 106616
Table 8
Distribution and comparison of high and very high seismic landslide hazard zone over geologic units of Anchorage, Alaska under 2% probability of exceedance in 50
years.
Geo units Area This study (Wang and Rathje, 2015) (Jibson and Michael, 2009)
proportion
Unit Contribution Unit Contribution Unit Contribution
proportion proportion proportion proportion proportion proportion
(a) Area proportion represents the percentage of the region covered by each geologic unit.
(b) Unit proportion represents the percentage of grid cells considered as high and very high hazards within each geologic unit.
(c) Contribution proportion represents the contribution of each geologic unit to all the number of D > 5 cm grid cells in this region.
9
C. Li et al. Engineering Geology 301 (2022) 106616
Fig. 7. Effect of the variability of the unit weight of soil on critical slope angle.
φ’. Other input parameters are the same as those mentioned above. The
reference value for the displacement threshold x is 5 cm and that for the
hazard level is 2% probability of exceedance in 50 years, and the
lognormal distribution is default satisfied for c’ and φ’.
Fig. 8 shows the critical slope angle comparisons of different COV
combinations and geologic units denoted by the continued number. The
variability of the COV combination has a significant influence on the
critical slope angle, especially for geologic units with friction of zero,
which are denoted as 5, 7, 11, and 15. For example, the result of the
geologic unit bc varies from 20◦ to 32◦ as the COV combination de
creases by 8◦ when incorporating a higher COV combination. However,
COV combinations have less effect on other geologic units, and the
critical slope angle decreases slightly with an increase in COV. For
example, geologic Unit 4 (geologic unit b) shows the best capacity to
resist instability during shaking and expresses only a one-degree in
crease between the lowest and highest combinations. This phenomenon
is attributed to the different selection of COV combinations in this sec
tion. The large step of COV-undrained from 0.1 to 0.5 led to a considerable
decrease in the results for these geologic units containing only
10
C. Li et al. Engineering Geology 301 (2022) 106616
undrained parameters (denoted as 5, 7, 11, and 15 in Fig. 8). The cor 6. Conclusions
responding outcomes of other geologic units, which have both strength
parameters and show minor variations, are caused by the slight decrease This paper introduces a comprehensive probabilistic method for
in COV- friction, from 0.15 to 0.05, and an unchanged value of 0.2, which seismic landslide hazard mapping procedures that include the prepara
is assigned to COV- cohesion. For example, if the friction angle is set to 30◦ , tion of a slope map from DEM and equivalent landslide hazard cate
then +σ leads to 34.5◦ , 33, ◦ and 31.5, ◦ for three COVs (corresponding gories. Instead of using displacement threshold, the critical slope angle,
to 0.15, 0.1, and 0.05, respectively). The difference between them is determined by selecting the seismic hazard level, leads to a modified
almost close to 1◦ , which coincides with the disparities of the critical seismic landslide hazard category, which is more convenient for seismic
slope angle in Fig. 8. Generally, the lower the COV value, the higher the landslide hazard mapping of a region. The proposed hazard category
slope angle, which indicates a low risk of co-seismic landslides, and less was applied to the region of Anchorage for the seismic landslide hazard
area will be regarded as a high landslide hazard. This suggests that the mapping process. The results in hazard maps show that the area of
accurate determination of soil strength parameters is of great signifi hazardous zone of most geologic units in present study is smaller than
cance for seismic landslide hazard mapping. The comparison results that in the deterministic analysis, except for the geologic units that were
further indicate that ignoring the uncertainty of slope properties leads to assigned a set of fixed values and ignored the variability of relevant
a considerable difference; therefore, a deterministic group value of input parameters. As far as these geologic units are concerned, logic-tree
parameters might result in an overestimation of resistance to landslides analysis shows an underestimated result, whereas the deterministic
or an underestimated mapping result. method overestimates the hazard levels of these units. Hence, an
appropriate probability distribution, such as a lognormal distribution, is
5. Discussion sufficient for the assessment of the occurrence of a seismic landslide.
However, a suitable DEM resolution needs to be considered in the
The proposed method provides a new type of probabilistic seismic mapping process for a considerable difference with a factor of 2 between
landslide hazard mapping procedure that incorporates epistemic un the comparisons of higher and lower resolutions. The seismic landslide
certainty associated with the variability of input parameters (i.e., c’, φ’, hazard mapping results or critical slope angles are more sensitive to
γ, t, and m) and includes aleatory variability of seismic motions and variations in the COV combination than the unit weight of soil. A single
different empirical displacement models. The influence of the input representative value of the unit weight of soil was sufficient for the
parameters and probability distributions on the final mapping results are seismic landslide hazard mapping process of a region.
contained in the middle process of the yield acceleration calculation. In This study describes a feasible method to explore the effects of DEM
addition, the effects of seismic motions and different empirical resolution. However, the proposed method is applicable to shallow
displacement models on the mapping procedure are expressed by landslides because of the hypothesis of the infinite slope model and rigid
displacement hazard curves along with regression empirical formula sliding block modeling. Hence, caution should be exercised when the
tions. An abnormal situation existed in the simulation of the strength proposed method is used to analyze deep landslides or soft stratum
parameters by the assumed normal distributions. A negative c’ and φ’ slopes. The development of more advanced technology together with the
might be obtained randomly, which leads to negative yield acceleration proposed method, which considers complex factors such as spatial cor
based on eqs. (1) and (2). However, this situation has been completely relations and topographic amplification, would considerably promote
avoided by the truncation operation for the ky distribution, which is the development of seismic landslide hazard maps on a regional scale.
discussed in the section above. This abnormal situation for negative c’
and φ’ in the Monte Carlo simulation could barely exist because of a CRediT authorship contribution statement
small COV and a large expectation value in this study, even for high COV
combinations. Chao Li: Conceptualization, Methodology, Software, Investigation,
As shown in Fig. 1, the entire process for seismic landslide hazard Formal analysis, Writing – original draft. Gongmao Wang: Visualiza
mapping is divided into two separate parts: the preparation of the slope tion, Investigation. Jianjian He: Resources. Yubing Wang: Conceptu
map and the critical slope angle distributions over all geologic units. alization, Funding acquisition, Resources, Supervision, Writing – review
Because the calculation of the critical slope angle and the accompanying & editing.
seismic landslide hazard category is independent of the mapping process
using a computation platform such as ArcGIS, it is convenient to adopt
the Monte Carlo method and to consider the effect of DEM resolutions Declaration of Competing Interest
without being limited to the computation cost. This differs significantly
from previous studies that were thought to be grid cell-dependent The authors declare that they have no known competing financial
analytical methods and were limited to the number of grid cells in a interests or personal relationships that could have appeared to influence
region. For instance, single and low resolution (such as 6 m, 10 m or 30 the work reported in this paper.
m) were usually adopted by previous studies (Jibson and Michael, 2009;
Wang et al., 2017; Wang and Rathje, 2015) to accelerate the computa Acknowledgments
tion speed and for the limited data storage capacity, considering that the
number of grid cells increases significantly as the resolution changes. The authors acknowledge the National Natural Science Foundation
However, the entire process with a resolution of 1.5 m, compared to a of China (No. 51808490) and the Chinese Program of Introducing Tal
resolution of 6 m adds few computational cost in this study. Therefore, it ents of Discipline to University (the 111 Project, B18047) for their
shows the superiority of the proposed approach in regional mapping, financial support. Technical discussions with Dr. Ellen M. Rathje of the
compared with common methods of grid cell dependency. The proposed University of Texas at Austin were very beneficial and gratefully
mapping procedure in this study employs Newmark’s method and an acknowledged.
infinite slope model. Hence, caution should be exercised when applied
to the analysis of deep landslides or soft stratum slopes. The variability References
of the input parameters was considered to be consistent within the same
geologic units in this study. Therefore, a better improvement in the Alamanis, N., Dakoulas, P., 2019. Effect of spatial variability of soil properties on the
simulation of regional application necessitates studies on the spatial stability and permanent seismic displacements of highway slopes. In: Proceedings of
the XVII ECSMGE-2019.
correlation between grid cells and a rigorous probability distribution. Bojadjieva, J., Sheshov, V., Bonnard, C., 2017. Hazard and risk assessment of
earthquake-induced landslides—case study. Landslides 15, 161–171.
11
C. Li et al. Engineering Geology 301 (2022) 106616
Bray, J.D., Travasarou, T., 2007. Simplified procedure for estimating earthquake-induced Peng, W.F., Wang, C.L., Chen, S.T., Lee, S.T., 2009. Incorporating the effects of
deviatoric slope displacements. J. Geotech. Geoenviron. 133, 381–392. topographic amplification and sliding areas in the modeling of earthquake-induced
Capolongo, D., Refice, A., Mankelow, J., 2002. Evaluating earthquake-triggered landslide landslide hazards, using the cumulative displacement method. Comput. Geosci. 35,
hazard at the basin scale through GIS in the Upper Sele river Valley. Surv. Geophys. 946–966.
23, 595–625. Phoon, K.K., Kulhawy, F.H., 1999a. Characterization of geotechnical variability. Can.
Cherubini, C., 2000. Reliability evaluation of shallow foundation bearing capacity on c’, Geotech. J. 36, 612–624.
φ’ soils. Can. Geotech. J. 37, 264–269. Phoon, K.K., Kulhawy, F.H., 1999b. Evaluation of geotechnical property variability. Can.
Dreyfus, D., Rathje, E.M., Jibson, R.W., 2013. The influence of different simplified Geotech. J. 36, 625–639.
sliding-block models and input parameters on regional predictions of seismic Rathje, E.M., Saygili, G., 2009. Probabilistic assessment of earthquake-induced sliding
landslides triggered by the Northridge earthquake. Eng. Geol. 163, 41–54. displacements of natural slopes. Bull. N. Z. Soc. Earthq. Eng. 42, 18–27.
Fenton, G.A., Griffiths, D.V., 2003. Bearing-capacity prediction of spatially random c - ϕ Refice, A., Capolongo, D., 2002. Probabilistic modeling of uncertainties in earthquake-
soils. Can. Geotech. J. 40, 54–65. induced landslide hazard assessment. Comput. Geosci. 28, 735–749.
Fenton, G.A., Griffiths, D.V., 2004. Reply to the discussion by R. Popescu on “Bearing Saade, A., Abou-Jaoude, G., Wartman, J., 2016. Regional-scale co-seismic landslide
capacity prediction of spatially random c - ϕ soils”. Can. Geotech. J. 41, 368–369. assessment using limit equilibrium analysis. Eng. Geol. 204, 53–64.
Gallen, S.F., Clark, M.K., Godt, J.W., Roback, K., Niemi, N.A., 2017. Application and Saygili, G., Rathje, E.M., 2008. Empirical predictive models for earthquake-induced
evaluation of a rapid response earthquake-triggered landslide model to the 25 April sliding displacements of slopes. J. Geotech. Geoenviron. 134, 790–803.
2015 Mw 7.8 Gorkha earthquake, Nepal. Tectonophysics 714-715, 173–187. Saygili, G., Rathje, E.M., 2009. Probabilistically based seismic landslide hazard maps: an
Jibson, R.W., 2007. Regression models for estimating coseismic landslide displacement. application in Southern California. Eng. Geol. 109, 183–194.
Eng. Geol. 91, 209–218. Shinoda, M., Miyata, Y., 2017. Regional landslide susceptibility following the Mid
Jibson, R.W., 2011. Methods for assessing the stability of slopes during earthquakes—a NIIGATA prefecture earthquake in 2004 with NEWMARK’S sliding block analysis.
retrospective. Eng. Geol. 122, 43–50. Landslides 14, 1887–1899.
Jibson, R.W., Michael, J.A., 2009. Maps showing seismic landslide hazards in Anchorage, Tao, W., Shuren, W., Jusong, S., Peng, X., 2013. Application and validation of seismic
Alaska. In: US Geological Survey Scientific Investigations Map 3077, scale 1:25,000, landslide displacement analysis based on Newmark model: a case study in Wenchuan
11-p. pamphlet. Available at URL. https://pubs.usgs.gov/sim/3077. earthquake. Acta Geologica Sinica (English Edition) 87, 393–397.
Jibson, R.W., Harp, E.L., Michael, J.A., 1998. A Method for Producing Digital Tsai, H.Y., Tsai, C.C., Chang, W.C., 2019. Slope unit-based approach for assessing
Probabilistic Seismic Landslide Hazard Maps: An Example from the Los Angeles, regional seismic landslide displacement for deep and shallow failure. Eng. Geol. 248,
California, Area. US Department of the Interior, US Geological Survey. 124–139.
Jibson, R.W., Harp, E.L., Michael, J.A., 2000. A method for producing digital Wang, Y., 2014. Probabilistic Assessments of the Seismic Stability of Slopes:
probabilistic seismic landslide hazard maps. Eng. Geol. 58, 271–289. Improvements to Site-Specific and Regional Analyses(Dotoral Dissertation).
Lee, S.T., Yu, T.T., Peng, W.F., Wang, C.L., 2010. Incorporating the effects of topographic University of Texas at Austin.
amplification in the analysis of earthquake-induced landslide hazards using logistic Wang, K.L., Lin, M.L., 2009. Development of shallow seismic landslide potential map
regression. Natural Hazards and Earth System Science 10, 2475–2488. based on Newmark’s displacement: the case study of Chi-Chi earthquake, Taiwan.
Liu, J.M., Gao, M.T., Wu, S.R., Wang, T., Wu, J., 2016. A Hazard Assessment Method for Environmental Earth Sciences 60, 775–785.
potential Earthquake-Induced Landslides - a Case Study in Huaxian County, Shaanxi Wang, Y., Rathje, E.M., 2015. Probabilistic seismic landslide hazard maps including
Province. Acta Geologica Sinica-English Edition 90, 590–603. epistemic uncertainty. Eng. Geol. 196, 313–324.
Liu, J., Shi, J., Wang, T., Wu, S., 2017. Seismic landslide hazard assessment in the Wang, H., Wang, G., Wang, F., Sassa, K., Chen, Y., 2008. Probabilistic modeling of
Tianshui area, China, based on scenario earthquakes. Bull. Eng. Geol. Environ. 77, seismically triggered landslides using Monte Carlo simulations. Landslides 5,
1263–1272. 387–395.
Martino, S., Battaglia, S., Delgado, J., Esposito, C., Martini, G., Missori, C., 2018. Wang, Y., Song, C., Lin, Q., Li, J., 2016. Occurrence probability assessment of
Probabilistic Approach to provide scenarios of Earthquake-Induced Slope failures earthquake-triggered landslides with Newmark displacement values and logistic
(PARSIFAL) Applied to the Alcoy Basin (South Spain). Geosciences 8. regression: the Wenchuan earthquake, China. Geomorphology 258, 108–119.
McCrink, T.P., 2001. Regional earthquake-induced landslide mapping using Newmark Wang, T., Liu, J., Shi, J., Wu, S., 2017. The influence of DEM resolution on seismic
displacement criteria, Santa Cruz County, California. In: Ferriz, H., Anderson, R. landslide hazard assessment based upon the Newmark displacement method: a case
(Eds.), Engineering Geology Practice in Northern California, California Division of study in the loess area of Tianshui, China. Environmental Earth Sciences 76.
Mines and Geology Bulletin 210/Association of Engineering Geologists Special Wang, T., Liu, J.M., Shi, J.S., Gao, M.T., Wu, S.R., 2020. Probabilistic seismic landslide
Publication, 12, pp. 77–92. hazard assessment: a case study in Tianshui, Northwest China. J. Mt. Sci. 17,
Meehan, C.L., Vahedifard, F., 2013. Evaluation of simplified methods for predicting 173–190.
earthquake-induced slope displacements in earth dams and embankments. Eng. Wasowski, J., Keefer, D.K., Lee, C.T., 2011. Toward the next generation of research on
Geol. 152, 180–193. earthquake-induced landslides: current issues and future challenges. Eng. Geol. 122,
Miles, S.B., Keefer, D.K., 2009. Evaluation of CAMEL — comprehensive areal model of 1–8.
earthquake-induced landslides. Eng. Geol. 104, 1–15. Zhang, Y.B., Xiang, C.L., Chen, Y.L., Cheng, Q.G., Xiao, L., Yu, P.C., Chang, Z.W., 2019.
Newmark, N.M., 1965. Effects of earthquakes on dams and embankments. Géotechnique Permanent displacement models of earthquake-induced landslides considering near-
15, 139–160. fault pulse-like ground motions. J. Mt. Sci. 16, 1244–1257.
Niño, M., Jaimes, M.A., Reinoso, E., 2014. Seismic-event-based methodology to obtain Zhao, H., Song, E.X., 2012. A method for predicting co-seismic displacements of slopes
earthquake-induced translational landslide regional hazard maps. Nat. Hazards 73, for landslide hazard zonation. Soil Dyn. Earthq. Eng. 40, 62–77.
1697–1713.
12