Probabilistic Analysis of Weathered Soil Slope in South Korea

Rainfall is a major trigger of shallow slope failures, and it is necessary to consider the spatial correlation of soil properties for probabilistic analysis of slope stability in heterogeneous soil. In this study, a case study of a weathered soil slope in Korea was performed to identify the rainfall-induced landslides considering the spatial variability of the soil properties and the probabilistic rainfall intensity depending on the return period and the rainfall duration. Various laboratory tests were performed to determine the physical properties of the site, and an electrical resistivity survey was carried out to understand the soil strata. Cohesion, friction angle, and permeability were considered as random variables considering the spatial variability, and the probabilistic rainfall intensities for return period of 2, 5, 10, 50, 100, and 200 years were used to consider the effects of rainfall infiltration. -e results showed that a probabilistic framework can be used to efficiently consider the spatial variability of soil properties, and various slope failure patterns were identified according to the spatial variability of the soil properties and the probabilistic rainfall intensity.


Introduction
Shallow slope failure (typically 1-3 m deep) due to heavy rainfall during rainstorms and typhoons is common in mountain areas and take the form of translational slides, which form parallel to the original surface [1,2].Over the past decade, many studies on rainfall pattern changes due to climate change have been carried out, and climate change has resulted in changes in rainfall patterns that can cause less frequent, but more intense rainfall events [3,4].
One of the main triggering factors for landslides is heavy rainfall [5][6][7][8], and the abnormal climate and localized heavy rainfall caused by climate change may lead to more frequent landslides [4,[8][9][10][11][12].For this reason, the prediction of rainfall-induced landslides is becoming more important, and the infinite slope model is usually implemented for the stability analysis of natural slopes; this model is appropriate when the horizontal dimensions of the surface are relatively much larger than the vertical depth of the potential failure slope.In addition, physically-based models have frequently been used in the framework of early warning systems devoted to rainfall-induced landslide hazard monitoring [13], and infinite slope model is commonly used to assess the factor of safety.However, the traditional infinite slope equation assumes homogeneous or averaged soil properties, and the potential failure surface is always fixed at the base of the slope [14] or wetting front depth (WFD).Although the uncertainties of the soil can be considered through the probabilistic analysis, the potential failure surface is still fixed and the probability of failure can be underestimated if the soil properties are treated as random variables without considering the spatial variability.
Geomorphological processes can lead to soil regions characterized by a degree of spatial heterogeneity [15], and the spatial variation can occur at various scales depending on the soil forming factors such as parent material, climate, topography, time, and so on.Describing and understanding the complexity of the interacting processes of soil formation has been a challenge, and geostatistics assumes that the variation of a property such as soil is continuous, which is generally more realistic for soil [16].If the spatial variability of soil properties is considered in the probabilistic analysis, the soil strength can vary depending on the depth considering the spatial correlation structure, and the potential failure surface can occur at various depths.Considering these features, some recent studies focused on the probabilistic analysis of slope stability in heterogeneous soil considering the spatial variability.
For example, Cho [17] identified the effect of spatial variability of unit weight and shear strength parameters using a limit equilibrium method, also, Griffiths et al. [18] investigated the effect of the spatial variability of shear strength parameters using the random finite-element method (RFEM), and Griffiths et al. [14] performed infinite slope analysis considering spatial variability of soil and showed that the potential failure surface can occur at various depths depending on spatial variability.Jiang et al. [19] identified the effect of spatial variability of shear strength parameters using a nonintrusive stochastic finite element method.Cho [20] discussed the effect of spatial variability of permeability on infinite slope stability and the distribution of failure depth due to rainfall infiltration.However, the saturated permeability is often depth-dependent, and Dou et al. [21] investigated the effect of spatial variability of permeability considering nonstationary random field of the saturated permeability as an extension of Cho's [20].In a similar vein, Li et al. [22] conducted reliability analysis of an infinite slope considering the linear trend of shear strength parameters in a random field.Li et al. [23] also proposed a multiple response-surface method for efficient slope reliability analysis considering spatial variability, and Cai et al. [24] conducted a cross-correlation analysis to determine the impact of heterogeneity of permeability, soil cohesion, and soil friction angle on slope stability.
In previous studies, the random variables most frequently considered when analyzing spatial variability of the soil properties are shear strength parameters (cohesion, c and friction angle, ϕ), and the permeability has also been considered in slope stability analysis due to rainfall infiltration.In practice, other important factors for estimating the probability of landslides occurrence are the intensity of rainfall and duration time [25], which are climatic factors.Rainfall is a major cause of landslides, and it is known that spatial and temporal variability is very high [26].erefore, it is important to consider and evaluate the appropriate probabilistic rainfall intensity for predicting the actual probability of landslide occurrence.However, little is known about the role of spatial variability of soil properties in probabilistic analysis of slope stability considering the probability rainfall intensity.
e purpose of this study was to identify the probability of failure of rainfall-induced landslides considering the spatial variability of soil properties and probabilistic rainfall intensity.A case study of shallow slope failure of weathered residual soil slope in Jangheung, Korea, was performed to verify the probabilistic analysis framework.e soil strata of the slope were identified, and the site investigation point was selected through electrical resistivity survey.en, the soil physical properties and infiltration characteristics of unsaturated soil on natural slope were investigated, and two shear strength parameters (c and ϕ) and permeability were considered as random fields.A slope stability analysis was performed using an infinite slope model, and the probabilistic rainfall intensity for 2, 5, 10, 50, 100, and 200 years frequency for the study area was considered.A series of Latin hypercube sampling-(LHS-) [27] based Monte Carlo simulations was conducted to investigate the effect of the spatial variability of soil properties on the mechanism of slope failure during periods of rain infiltration.Consequently, the probability of slope failure according to the probabilistic rainfall intensity and spatial variability of soil properties was estimated.How soil properties and climatic factors affect slope stability were then discussed.

Random Field
Nearly all natural soils are highly variable in their properties, and their variability shows a spatial correlation.e spatial variability of soil properties can be effectively considered using random field theory, and probabilistic analyses that incorporate the spatial variability of soil properties as random fields are more appropriate to consider the uncertainty of soil than those considering soil properties as a single random variable.

Spatial Variability of Soil.
Because of complex geological and environmental processes, soil is inherently heterogeneous, and its properties can be highly variable and spatially correlated in the vertical and horizontal directions.e spatial correlation of soil properties is known to influence the geotechnical response of soil, and it brings unavoidable uncertainty in design, leading to unexpected soil responses [15,28].ese uncertain spatial properties can be characterized using random field theory [29,30].Vanmarcke [29] used a scale of fluctuation (SOF, δ) to describe the extent of how soil properties are spatially correlated.Various methods are available to estimate the SOF, and the simplest approach is to fit the theoretical autocorrelation function (ACF, ρ) to the empirical ACF [31][32][33][34].However, determining a theoretical ACF may not be easily implemented because a large amount of data is required.erefore, some theoretical ACFs are usually used to characterize the spatial correlation of soil properties, and the single exponential ACF has been widely used to model the inherent spatial variability of soil properties in probabilistic analysis of slope stability [23]: where l represents the autocorrelation distance and reflects the rate at which the correlation decays between two points (x 1 , x 2 ).e SOF implied by the single exponential autocorrelation function is equal to twice the value of the autocorrelation distance (δ � 2l).

Random Fields.
In this study, the Karhunen-Loève expansion (KLE) was adopted to generate random fields because it is an efficient method for random field 2 Advances in Civil Engineering discretization with a desired level of accuracy and provides the greatest accuracy when an exponential ACF is used [35].e KLE of a random field with a mean value (μ ω ) and a variance (σ 2 ω ) is given by Spanos and Ghanem [36]: where λ i and f i (x) are the eigenvalues and eigenfunctions of the covariance function, respectively, and ξ i (θ) represents the uncorrelated zero mean random variables.For practical implementation, the discretization of the random field ω(x, θ) is obtained by truncating the series expansion at the Mth term: e accuracy of the represented random field depends on the number of terms used in the KLE expansion, and the number of required terms is determined according to the ratio of the correlation length and the domain size [37].
Normal random fields are often used for modeling uncertainties with spatial variability for mathematical convenience and due to a lack of available data, but they are not applicable in many situations where the random variable is always non-negative.erefore, the assumption of a lognormal distribution is appropriate as the soil properties used in this study are always non-negative [38,39].In (3), if μ ω is zero and σ ω is one, standard normal random fields are generated by KLE, and it can be transformed into lognormal random fields using (4) as follows: where μ ln(x) and σ ln(x) are equal to the mean and standard deviation of ln(x) of the underlying normal distribution.

Rainfall-Induced Landslide
In this section, the rainfall infiltration model and slope stability analysis model for unsaturated soils were described to assess the vulnerability of rainfall-induced landslides caused by rainfall, and then the probabilistic analysis procedure considering the spatial variability of soil properties was presented.

Rainfall Infiltration.
In order to perform infiltration analysis considering spatial variability, the flow of water in multilayered soils should be considered, and this can be performed using differential equations that satisfy Darcy's law [20,21,24,40].Alternatively, semianalytical and arealaveraged infiltration models have been established to study field-scale infiltration over soils with the variability of permeability based on the Green-Ampt [41] model [42,43].Chu and Marino [44] present a modified Green-Ampt model considering infiltration in multilayered soils. is model is able to deal with unsteady and steady rainfall events as well as ponding and nonponding conditions.erefore, in this study, the infiltration model by Chu and Marino [44] was adopted to consider the spatial variability of permeability.Figure 1 shows infiltration of rainfall into an n-layered soil profile with permeability (k) and initial volumetric water content (θ 0 ) for each layer.
When the wetting front is over layer n at location z (z n−1 < z ≤ z n ), the infiltration rate (i) under an initially ponded condition can be expressed as and cumulative infiltration (I) is e travel time to depth z can be expressed as where θ s is the saturated volumetric water content (L 3 /L 3 ), Ψ m is the pressure head (suction) (L), θ f is (θ s − θ 0 ), and the subscript n means the n-layer.In the case of nonponding, where the infiltration capacity of soil is greater than the rainfall intensity, all rainfall was assumed to infiltrate into the soil producing a piston flow.A detailed description of rainfall infiltration can be found in Chu and Marino [44].

Slope Stability of Unsaturated Soils.
Rainfall-induced landslides usually have shallow surface failure depth, being the main cause of the reduction of negative pore water pressure in unsaturated soil due to rainfall infiltration [5,45].Fredlund et al. [46] extended the shear strength criteria to unsaturated soils considering the matric suction, being possible to calculate the factor of safety along the potential failure surface (Figure 2) as where c ′ is the effective cohesion, σ ′ is the effective normal stress, ϕ ′ is the effective friction angle, u a is the pore air pressure, u w is the pore water pressure, W is the weight of the soil slice per unit base area, β is the slope angle, and ϕ b is the angle indicating the rate of increase of shear strength relative to increased matric suction.e shear strength contribution due to matric suction of unsaturated soil can be expressed as described by Fredlund and Rahardjo [45]: where θ is the volumetric water content at potential failure surface and θ s and θ r are the saturated and residual volumetric water content, respectively.Θ is normalized volumetric water Advances in Civil Engineering content, and the relationship between Θ and the matric suction (Ψ) was given by van Genuchten [47]: where α, n, and m are van Genuchten's soil-water characteristic curve (SWCC) parameters.Considering ( 9) and ( 10), ( 8) can be expressed as

Probabilistic Analysis.
e schematic view of the probabilistic slope analysis considering the spatial variability of soil properties and probabilistic rainfall intensity is shown in Figure 3.In order to consider the spatial variability of soil, random elds are generated for three soil properties (c, ϕ, and k s ), which are mainly considered as random variables in the existing studies on slope stability, and these are used as input values in the rainfall in ltration and slope stability analysis.
e probability rainfall intensity depending on the return period and the rainfall duration is applied to consider rainfall characteristics for the site, and the WFD due to rainfall inltration is estimated by Chu and Marino [44] model.Subsequently, the WFD is used for in nite slope stability analysis.As a result, the factor of safety for each depth and potential failure surface is estimated.is process is repeated to derive probabilistic results, and a series of LHS-based Monte Carlo simulations is conducted for e cient random sampling.

Soil Physical Properties.
e study area is located in Jangheung, Gyeonggi-do, South Korea (37 °45′03″N, 4 Advances in Civil Engineering

°56′21″E
). is area is underlain by a Gyeonggi gneiss complex that belongs to the Gyeonggi massif, and the bedrock is intensely weathered.e surface is covered by Pinus koraiensis, commonly known as the Korean pine.As a preliminary investigation, the slope length, depth, and gradient were measured for 5 slopes on the study area, and the slope which seems to be vulnerable to landslides (with steep and long slope) was chosen as the study slope.A topographic survey was carried out to clarify the geometrical characteristics of the slope using electronic distance measurement (EDM), slope length and angle approximately being 50 m and 31 °, respectively.Figure 4 shows the location and topographic surface of the study slope.
Electrical resistivity survey was carried out in order to identify soil strata from the resistivity distribution and to select the site investigation point.Resistar RS-100M equipment with multielectrode system ME-100 (Geofyzika A.s.) was used for a multielectrode survey using a Wenner-Schlumberger electrode array with electrode spacing of 2.5 m.Primary data were processed into resistivity cross sections by a 2-D inverse method using the software Res2DInv [48]. Figure 5 shows the resistivity tomography of the slope of the study area.
As shown in Figure 5, the direction of the soil strata can be regarded as parallel to the slope surface, and horizontal autocorrelation is higher than vertical autocorrelation for soil resistivity.
erefore, it is reasonable to generate Advances in Civil Engineering a random field for the vertical direction on the slope surface.Blue (red) indicates low (high) electrical resistivity.e area with the highest electrical resistivity was in good agreement with the location of the exposed outcrop.e upper slope shows a relatively low resistance value, and it can be assumed that the water content is high [49] and the strength is relatively weak.erefore, the cross section of the upper slope was selected as the analysis section.Hand auger borings were carried out primarily to identify soil depth and stratum, and two test pits with a size of 1 m × 1 m were dug.e test pits and soil sampling locations are depicted in Figure 6. e wet unit weight of each sampling depth was obtained by the sand replacement method, and three and two disturbed soil samples were obtained from the upper and lower slope test pits, respectively.e triaxial compression test and the constant head permeability test were carried out to determine the soil physical properties.In addition, the undisturbed samples were also obtained from the upper slope using block sampling and used for the SWCC test.e basic physical properties of soil and slope conditions are summarized in Table 1.
In accordance with the United Soil Classification System (USCS), the soil is classified as SM or SC, and the specific gravity (G s ) of soil is 2.62.e shear strength test was carried out on soil with a depth of 30, 60, and 120 cm, and the average values for cohesion and friction angle were 14.7 kPa and 17.9 °, respectively.e permeability test was also performed on the same depths, and the average saturated permeability (k s ) was 1.64 × 10 −6 m/s.e random variables selected for the stability analysis are two shear strength parameters (c ′ and ϕ ′ ) and saturated permeability.However, the measured data were small, and there was a limit to identifying the statistical characteristics of soil properties by using measured data.erefore, only the average values of the measured data for c, ϕ, and k s were used, and their COV and PDF were assumed through the literatures [15,30].e soil properties vary based on the depositional and stress history of the deposit, and its distribution characteristics can be quantitatively reflected through the probability density function (PDF).
e PDF type was assumed to be a lognormal distribution for all three parameters, and the coefficient of variation (COV) of c ′ , ϕ ′ , and k s , was assumed equal to 30%, 20%, and 50%, respectively.Phoon and Kulhawy [30] summarized the SOFs and COVs of inherent spatial variability of various soil properties; however, information of soil strength SOF is insufficient.erefore, the normalized SOFs (δ/D) were assumed to be 0.1, 0.2, 0.4, 0.8, 1.5, and 3.0 for vertical depth to slope surface (D) in order to observe the effect of spatial variability.e statistical properties of the three parameters are summarized in Table 2.
Rainfall infiltration into the soil begins from unsaturated conditions, and it is necessary to identify the permeability characteristics under unsaturated conditions.erefore, a laboratory test for SWCC was conducted, and parameters of the van Genuchten [47] model were estimated as α � 0.0291, m � 0.3473, and n � 1.532.e conditions of infiltration analysis are summarized in Table 3, and SWCC is shown in Figure 7.

Probabilistic Rainfall Intensity.
Rainfall information on the study area is very important because rainfall has a high spatial and temporal variability, and its probabilistic characteristics are very different depending on the region.e information of probability rainfall intensity at the study area was obtained from the Korea precipitation frequency data server of the Ministry of Land, Infrastructure and Transport (MOLIT, http://www.k-idf.re.kr).Figure 8 shows intensityduration-frequency (IDF) curves.
Rainfall intensity for the IDF curve increased as the return period increased, and the smaller the duration, the greater the difference in rainfall intensity.In this study, the probabilistic rainfall intensity for a return period (T r ) of 2, 5, 10, 50, 100, and 200 years was considered.Rainfall events with durations of 1, 2, 3, 6, 9, 12, 15, 18, and 24 hours were chosen for slope analysis, and rainfall intensity is summarized in Table 4.

Wetting Front Depth and the Factor of Safety.
Figure 9 shows the change of wetting front depth (WFD) and the factor of safety depending on rainfall duration.Soil  6 Advances in Civil Engineering saturation occurs from the surface due to rainfall infiltration, and the WFD increased with the increase of rainfall duration.e deterministic approach using Chu and Marino [44] model was carried out to estimate the WFD and the factor of safety at the baseline and WFD according to the rainfall duration.e gray line indicates the WDF changes and the dotted lines the changes of FS at the baseline and at the WFD.WFD increases approximately 2 m when rainfall duration is equal to 24 h, and the factor of safety at WFD was drastically reduced due to the dissipation of negative pore water pressure in unsaturated soil as the rainfall began to penetrate into the ground.e factor of safety at the baseline showed a very small reduction rate considering only an increase in unit weight due to the rainfall infiltration.In the deterministic approach, the factor of safety until a rainfall of 24 h was greater than 1. e change in the WFD with the rainfall duration and normalized SOF is shown in Figure 10.
e average WFD showed a significant increase with increasing rainfall duration, but it increased very slightly with normalized SOF.When the return period was low, the    Wetting front depth (m) 0.0 0.5 WFD Base line F s WFD F s 0 300 600 900 1200 1500 Figure 9: Change in wetting front depth and factor of safety depending on rainfall duration using a deterministic approach.
difference in average WFD according to the return period was not significant even though there was a large difference in rainfall intensity, whereas the average WFD significantly changed with the return period for high rainfall durations.
Although there was no significant difference in average WFD with normalized SOF, there was a large variation in standard deviation.e WFD appeared at various depths depending on the spatial variability of the permeability, and the distribution exhibited a wider range as the rainfall duration increased (Figure 11).

Effect of Spatial Variability and Probabilistic Rainfall
Intensity.In this study, the potential failure surface was defined as the deepest depth with a factor of safety of less than 1, and if the factor of safety is less than 1 at any depth, the slope failure was considered to have occurred.e probability of failure is significantly affected by the return period and the rainfall duration as shown in Figure 12. e effect of the return period increases with rainfall duration because WFD shows a large difference as the rainfall duration increased.However, the effect of the return period decreased as the return period increased, and there was no significant difference in the probability of failure when the return period was over 100 years.is is because the surface runoff occurs when the rainfall intensity exceeds a certain level, and the infiltration is dominated by the infiltration capability of the soil after reaching the initial ponding condition.e spatial variability of soil properties was found to have a major influence.If spatial variability of soil is not considered, the normalized SOF is considered infinity (δ/D � ∞), and its results showed the lowest probability of failure. is is because the soil is considered to be homogeneous, and it does not take into account the potential failure surface at various depths as the failure surface is fixed to the WFD or the baseline.e impact of the return period and the rainfall duration on the probability of failure also varied depending on the spatial variability of soil properties, and the variation of the failure probability with the return period and the rainfall duration was larger as the normalized SOF decreased.e variation of the probability of failure according to normalized SOF for the return period of 2 and 100 years is shown in Figure 13.e probability of failure decreased with increasing normalized SOF, and it rapidly decreased when the normalized SOF was less than approximately 1.As identified in various studies [14,20,22], the location of the potential failure surface appeared at various depths in consideration of the spatial variability of soil properties, and the larger the normalized SOF, the lower the probability of failure because the failure occurs at a narrower range of depths.e difference in the probability of failure according to the rainfall duration for the return period of 2 years is smaller than that of the return period of 100 years.e reason is that although the probability of failure is almost the same regardless of the return period for low rainfall duration under the same normalized SOF, the probability of failure increased with the increase of rainfall duration, and its growth rate was larger as the return period is higher as shown in Figure 12. Figure 14 shows an example of histograms for the frequency of the critical depth where the minimum factor of safety less than 1 occurs and the cumulative failure probability along the depth for 100-year return period and the rainfall duration of 18 hours.e left cumulative distribution and gray histogram are the results for the normalized SOF of 0.1, and the right cumulative distribution and black histogram are the results for the normalized SOF of 1.5.When the normalized SOF was 0.1, the failure frequency near the baseline was about 19%, whereas that of the normalized SOF of 1.5 was about 40%.e normalized SOF of 0.1 shows higher failure frequency at various depths than that of the normalized SOF of 1.5, and this is due to the effects of the shear strength parameter as found in the results of Griffiths et al. [14] and the distribution of WFD as shown in Figure 11.
e probability of failure according to the normalized SOF and the rainfall duration for 2 and 100 years return periods were compared with the case where the spatial variability of soil properties was not considered (i.e., homogeneous soil) for each return period, and the differences in probability of failure (Δp f ) are shown in Figure 15.e difference of the failure probability is significantly influenced by the normalized SOF, and the return period and the rainfall duration also have an effect on the difference in probability of failure when the normalized SOF is small.If the normalized SOF is small (i.e., δ/D � 0.1), a probabilistic difference of approximately 24∼30% can occur, while the spatial variability of soil properties has little effect on the normalized SOF of three or more.
is difference in the probability of failure is influenced by the slope and rainfall conditions and the statistical characteristics of soil properties, and therefore it can be different depending on the target slope.However, it was shown that the Advances in Civil Engineering probabilistic results of the slope stability exhibit noticeable differences that vary with the spatial variability of soil properties, and it is a well-known fact that soil is heterogeneous.erefore, there is a strong need to consider the spatial variability of soil in practice.In addition, rainfall is a major cause of landslides in South Korea, and the probability of failure of slopes is very different depending on the rainfall condition.erefore, it is also important to consider the probability rainfall intensity reflecting local rainfall characteristics.

Conclusions
e purpose of this study was to identify the failure of rainfall-induced landslides considering the spatial variability of soil properties and the probabilistic rainfall intensity through a case study of a weathered soil slope in Korea.An electrical resistivity survey was performed to understand the stratum of the slope.Although the electrical resistivity survey has some inherent limitations as it cannot provide direct information such as shear strength, the distribution of soil strata can be confirmed with images, and information can be provided on soil survey location and direction of random field generation.Herein, three soil properties (i.e., cohesion, friction angle, and permeability) were considered as random fields, and the IDF curves with a 2-, 5-, 10-50-, 100-, and 200-year return period were applied to identify the effects of probabilistic rainfall intensity.An infinite slope equation for unsaturated soil and multilayer infiltration model by Chu and Marino [44] was used to analyze the slope stability affected by rainfall, and a series of LHS-based Monte Carlo simulations was conducted for probabilistic analysis.
Probability of failure was significantly affected by the spatial variability of soil properties, and the results are similar to those of previous studies.e return period and rainfall duration also have a significant impact on slope stability due to the wetting front development.Although the rainfall conditions are the same, their effects on the slope stability were different according to the spatial variability of soil properties, and their effects were more sensitive with decreasing normalized SOF. e probability of failure increased as the return period and the rainfall duration increased.However, the increase rate of the probability of failure decreased as the return period increased, and there was no significant difference in probability of failure when the return period is over 100 years because rainfall infiltration is dominated by the infiltration capability of the soil.
Compared with the results of traditional probabilistic analysis without consideration of the spatial variability of soil properties, the difference in probability of failure was    Advances in Civil Engineering found to be approximately 24∼30% when the normalized SOF was small (i.e., δ/D � 0.1), and there was little effect of the spatial variability of soil properties for the normalized SOF more than 3. Naturally, the soil is inhomogeneous, and vertical spatial variability of soil properties is typically large.In addition, the probability of failure according to the probability rainfall intensity is greatly influenced by spatial variability of soil properties.erefore, probabilistic analysis with consideration of spatial variability of soil properties is very important to obtain reasonable results.Although this study applied various SOFs and identified their effects, it is necessary to apply appropriate SOFs to the site in order to accurately estimate the probability of failure.It is well known that the spatial variability in the vertical direction is larger than that in the horizontal direction.is fact was confirmed herein indirectly through an electrical resistivity survey, and it was also confirmed that it is reasonable to generate a random field for the vertical direction on the slope surface.
Additionally, this study used the IDF curve generated based on past rainfall data, but the impact of climate change was not reflected.Recently, however, a modified IDF curve has been proposed considering the climate change, and the impacts of climate change can be reflected by the modified IDF curve.

Figure 1 :Figure 2 :Figure 3 :
Figure 1: In ltration into a layered soil pro le under an initially ponded condition.

Figure 4 :Figure 5 :
Figure 4: Location of the study area and surface plot.

Figure 8 :
Figure 8: IDF curves for study site.

Figure 10 :Figure 11 :
Figure 10: Change in wetting front depth depending on rainfall duration and normalized SOF.

Figure 12 :Figure 13 :
Figure 12: Change in the probability of failure depending on the return period and the rainfall duration.

Figure 14 :
Figure 14: for the frequency of the critical depth and the cumulative failure probability along with the depth (T r � 100 years and 18 hours rainfall duration).

Figure 15 :
Figure 15: Difference in probability of failure with the results for homogeneous soil (δ/D � ∞) for (a) T r � 2 years and (b) T r � 100 years.

Table 1 :
Summary of basic soil properties and slope conditions.

Table 2 :
Statistical properties of shear strength parameters.

Table 3 :
Conditions of infiltration analysis.
s : saturated volumetric water content.

Table 4 :
Probabilistic rainfall intensity (mm/h) of the study area.