A Novel Radial Basis Function Approach for Inﬁltration-Induced Landslides in Unsaturated Soils

: In this article, the modeling of inﬁltration–induced landslides, in unsaturated soils using the radial basis function (RBF) method, is presented. A novel approach based on the RBF method is proposed to deal with the nonlinear hydrological process in the unsaturated zone. The RBF is ﬁrst adopted for curve ﬁtting to build the representation of the soil water characteristic curve (SWCC) that corresponds to the best estimate of the relationship between volumetric water content and matric suction. The meshless method with the RBF is then applied to solve the nonlinear Richards equation with the inﬁltration boundary conditions. Additionally, the ﬁctitious time integration method is adopted in the meshless method with the RBF for tackling the nonlinearity. To model the stability of the landslide, the stability analysis of inﬁnite slope coupled with the nonlinear Richards equation considering the ﬂuctuation of transient pore water pressure is developed. The validation of the proposed approach is accomplished by comparing with exact solutions. The comparative analysis of the factor of safety using the Gardner model, the van Genuchten model and the proposed RBF model is provided. Results illustrate that the RBF is advantageous for reconstructing the SWCC with better estimation of the relationship than conventional parametric Gardner and van Genuchten models. We also found that the computed safety factors signiﬁcantly depend on the representation of the SWCC. Finally, the stability of landslides is highly affected by matric potential in unsaturated soils during the inﬁltration process.


Introduction
As severe rainfall events become common due to climate change, rainfall-induced landslides are found more frequent around the world, causing damages and fatalities exposed to landslides dramatically increased [1,2]. The infiltration-induced landslides in unsaturated soils is usually very close to the ground surface in which the shallow landslide may occur in the period of heavy rainfall subjected to the infiltration process. Shallow landslides mainly occurred in the foothill areas involving sliding surface which is commonly occurred within the vadose zone in unsaturated soils [3,4]. The triggering of landslides from rainfall requires the infiltration process in which the slope failure is determined by the fluctuation of matric potential in unsaturated soils [5].
The conventional slope stability analyses mainly focus on the mechanical behavior of the slope materials to calculate factors of safety [6]. The infinite slope model is often used to analyze the shallow landslide in which the sliding surface along the failure plane. A number of studies, presuming a steady groundwater level associated with the stability analysis of infinite slope, have been conducted to recognize the rainfall-induced landslides [7][8][9]. Iverson (2000) proposed the exact solution of the simplified form of Richards equation for one-dimensional linear diffusion equation [10]. Later, Baum et al. (2002) developed the transient rainfall infiltration and grid-based regional slope-stability (TRIGRS) model for transient rainfall infiltration and grid-based regional slope-stability analysis [11]. The TRIGRS model incorporated with the geographical information system has been successfully applied for assessing shallow landslides [12]. The linear diffusion equation may not comprehensively consider the influence of the soil water characteristic curve (SWCC) and the matric potential [13,14]. Accordingly, the importance of infiltration and pore water pressure fluctuation through unsaturated soil is neglected.
Studies showed that the Richards equation is highly nonlinear due to the high nonlinearity of physical behavior of unsaturated soils [15,16]. The SWCC is the reason for the nonlinear relationship of the hydraulic conductivity in unsaturated soils. Several empirical models, such as the Gardner model, the van Genuchten model, et al., have been proposed to describe the SWCC [17][18][19]. Due to the complexity, the Richards equation considering the SWCC in unsaturated soils is usually solved using the numerical methods. Mesh-based numerical techniques such as the finite difference method and the finite element method are well documented and typically used to solve the unsaturated flow equation in the past [20][21][22].
In this study, we propose an innovative approach for rainfall-induced landslides, emphasizing on the importance of infiltration processes to predict landslide hazard. With the superior capability of dealing with different kinds of partial differential equations, the radial basis function (RBF) method is proposed for solving the Richards equation considering the SWCC. Proposed by Hardy in 1971, the RBF was first used for scattered data interpolation [23]. Instead of using the empirical Gardner model or the van Genuchten models, we first adopt the RBF for curve fitting to build the representation of the SWCC. The RBF method is then applied to solve the nonlinear Richards equation with the infiltration using designated boundary conditions. Additionally, the fictitious time integration method [24] is integrated into the meshless method with the RBF for tackling the nonlinearity. To model the stability of the landslide, the infinite slope stability analysis using the nonlinear Richards equation with the fluctuation of transient pore water pressure is developed. The infinite slope model is one of the most common models for calculating the factor of safety based on limit equilibrium analysis [7,25]. The infinite slope model is to determine the balance between shear stress and shear strength of the sliding failure plane [26]. The elastoplastic constitutive model for the sliding failure plane is based on the Mohr-Coulomb model which is originally used in the infinite slope model. The relationship between the SWCC and the shear strength has been developed in which the shear strength of unsaturated soils at any specified value of suction is described [27,28]. Accordingly, this study emphasizes shallow translational landslides, considering the infinite slope model for the infiltration in unsaturated soils. The validity of the proposed approach is accomplished by comparing with exact solutions. Application examples are also conducted. The methodology is described as following section.

The Richards Equation
The unsaturated flow in soils can be described by the following variably saturated flow equation [21] where K is hydraulic conductivity, H is groundwater head, described as H = h e + h, h denotes suction head (or matric potential), h e denotes elevation head, S s is specific storage, S w is effective saturation, C(h) is specific capacity, described as C(h) = ∂θ /∂h, θ is volumetric water content, t is time, and ∇ is gradient. By considering the soils to be anisotropic and heterogeneous in Equation (1), the following equation is acquired in which K x (h), K y (h), and K z (h) are functions of hydraulic conductivity in lateral and vertical directions, respectively. To model the variably saturated flow in a slope, the following definition of the elevation head by Iverson (2000) is considered [10] h e = x sin(β) + z cos(β), where β denotes slope angle, x direction is along the slope and z direction is normal to the slope. By substituting Equation (3) into Equation (2), we obtain the following equation For one-dimensional flows in z direction through an unsaturated inclined slope, Equation (4) is simplified to The above equation can be rewritten as follows.
where K s denotes saturated hydraulic conductivity, and K r denotes relative hydraulic conductivity, described as K r = K z (h)/K s . Since relative hydraulic conductivity and volumetric water content are both functions of pressure head for unsaturated soils, the Richards equation becomes a highly nonlinear partial differential equation. The high nonlinearity of physical behavior of unsaturated soils can be described using SWCC.
To tackle the Richards equation as depicted in Equation (6), the empirical characteristic functions including the relative hydraulic conductivity function and SWCC are required.

The SWCC Model
The SWCC describes the relationship between matric potential and volumetric water content or effective saturation. The effective saturation is defined as in which S e denotes effective saturation, θ s and θ r are saturated and residual water contents, respectively. Since the SWCC contains the primary information for describing physical behavior of unsaturated soils, several SWCC models have been developed such as the Gardner model, the van Genuchten model, and the Fredlund and Xing model [17][18][19].
In these SWCC models, empirical fitting constants or fitted parameters are required to describe the SWCC. Proposed by Gardner in 1958, the Gardner exponential model was one of the mathematical models used to describe the SWCC [17]. The Gardner exponential model can be described as follows.
S G e = e α g h , where α g denotes pore size distribution parameter, and S G e denotes effective saturation proposed by Gardner. The Gardner model has a particularly simple form with only one fitting parameter. Another commonly used SWCC model was proposed by van Genuchten, described as where S v e denotes effective saturation proposed by van Genuchten [18], α v is related to the pore size, n v is a function of the spread of pore size distribution, and m v = 1 − 1/n v . The relative hydraulic conductivity from the Gardner model and van Genuchten model [29] can then be described as follows where K G r and K v r denote relative hydraulic conductivity from the Gardner and van Genuchten models, respectively. The Gardner model depicted in Equation (10) can be used to linearize the Richards equation. It is therefore widely utilized to formulate exact solution for the linearized Richards equation and in a wide range of engineering applications. However, the Gardner model may fit experimental data only within a limited range of matric suction [30]. For the van Genuchten model as depicted in Equation (11), the application of the quasi-linear analyses may often be limited to relatively simple initial and boundary conditions as well as unsaturated flow problems with simple pore geometry.
Instead of using the empirical Gardner model or the van Genuchten models, we propose the RBF to build the representation of the SWCC for curve fitting in this study. The RBF approach was originally presented by Hardy [23]. Numerical modeling with RBFs emphasizes on the reformation of unknown functions from data. The effective saturation to be solved is approximated as a linear combination of the RBF expressed as follows.
where S RBF e (r j ) denotes predicted effective saturation approximated by the RBF, M denotes the number of centers, κ j denotes weighting coefficients, r j denotes radial distance to the jth center point, r j = x − x s j , x s j denotes the jth center point, and φ(r j ) denotes RBF. Several RBFs including the multiquadric, polyharmonic spline, Gaussian, inverse multiquadric and radial polynomials (RP) have been developed [31], as listed in Table 1. In this study, we considered the RP RBF described as where k is the order of the RP, and N is the number of terms for the RP. By substituting Equation (13) in Equation (12), we obtain the following equation where β denotes a p × q coefficient matrix, κ denotes a q × 1 vector of unknown coefficients, p denotes the number of experimental data, and q = M × N, b denotes a p × 1 vector described EXP_p e denote experimental effective saturation. To investigate the accuracy of the proposed approach, the root mean square error (RMSE) and maximum relative error (MRE) between the experimental data of the effective saturation and the predicted effective saturation approximated by the RBF fitting model are defined as where S EXP_j e refers to jth experimental effective saturation data. Table 1. Types of RBF.

Type of RBF RBF
Notation: c denotes shape parameter.
To illustrate the proposed RBF model, we compare the proposed RBF fitting model to several SWCC models, including the Gardner model, van Genuchten model, and experimental data by Lu and Likos [32]. The soil type is glacial till. Best fitted parameters of the above SWCC models were obtained by minimizing the RMSE between experimental and predicted SWCC. The predicted SWCC, along with the testing data, is demonstrated in Figure 1. The fitted parameters and the RMSE are presented in Table 2. We then compare the RMSE of the proposed RBF model and those of the Gardner model and van Genuchten model. From Table 2, the RMSE of the Gardner model and van Genuchten model are only on the order of 10 −2 . The RMSE of the proposed RBF model, however, can reach the order of 10 −10 . It is significant that excellent agreement and highly accurate SWCC achieved by the proposed RBF model can be obtained. From the curve fitting results, we found that the RBF achieves the best estimate for volumetric water content and matric suction to build the representation of the SWCC of glacial till among those SWCC models, as shown in Figure 1. Table 2. The fitted parameters and RMSE in the example.

SWCC Model
Fitted Parameters RMSE

The RBF Methods for the Richards Equation
With the superior capability of dealing with partial differential equation, the RBF collocation method is one of the prominent methods for solving partial differential equations [33]. We assume that the pressure head can be estimated using the RBF approximation as where i λ denotes expansion coefficients. The RP RBF depicted in Equation (13) is utilized. The above equation can also be written as denotes the interpolation matrix, and T N ] ,..., . Proposed by Sarra [34], the n th derivative of Equation (17) is then computed as

The RBF Methods for the Richards Equation
With the superior capability of dealing with partial differential equation, the RBF collocation method is one of the prominent methods for solving partial differential equations [33]. We assume that the pressure head can be estimated using the RBF approximation as where λ i denotes expansion coefficients. The RP RBF depicted in Equation (13) is utilized. The above equation can also be written as where A = [φ(r 1 ), . . . , φ(r N )] denotes the interpolation matrix, and λ = [λ 1 , . . . , λ N ] T . Proposed by Sarra [34], the nth derivative of Equation (17) is then computed as By the same manner, we have the matrix form as where . It is often more efficient to use the derivative matrix as where D (n) denotes derivative matrix, and A −1 denotes inverse of the A matrix. The spatial derivative of the Equation (17) is then approximated as The nonlinear Richards equation, as depicted in Equation (6), is then discretized in space and time domains utilizing RBF to acquire the following equation as Since the above equation is nonlinear, the iteration method is required to obtain the solution of the nonlinear Richards equation. The solution process of the iteration method is described as follows.

The Fictitious Time Integration Method
For dealing with the nonlinear algebraic equations, as depicted in Equation (23), the iteration methods are mostly applied. An iteration method named the fictitious time integration method (FTIM) was developed by Liu and Atluri [24]. The FTIM tackles the nonlinear Richards equation by introducing a fictitious time. From Equation (7), we may obtain K r = K z (h)/K s and θ = S e (θ s − θ r ) + θ r . It is obvious that K r and θ are both functions of the matric potential. In the Gardner model K r = S G e where S G e is used for the fitting from the SWCC. However, K r can also be obtained directly from the experimental data using the RBF fitting. Accordingly, we may rewrite Equation (23) as follows.
where i = 1, . . . , Q is the collocation number. The iterative scheme of the FTIM for solving the Richards equation can then be expressed as follows.
where h k i denotes h at the kth discrete time, m denotes a factor ranging from zero to one, δ denotes fictitious time step size, ν denotes a fictitious damping coefficient, and τ k denotes fictitious time, defined as τ k = kδ. Equation (25) is adopted to solve the system of nonlinear Richards equation resulted from the RBF spatial discretization and Euler temporal discretization.
Beginning by giving an initial guess of the pressure head, we then utilize Equation (25) from t = 0 to the final elapsed time. The iterative procedure of the FTIM ends while one of the following convergence criteria is achieved.

Landslide Stability for Unsaturated Soils
To model the infiltration-induced landslides in unsaturated soils, the landslide stability analysis linked to the RBF approach with the fluctuation of transient pressure head is developed. The infiltration-induced landslides in unsaturated soils is usually very close to the ground surface in which the shallow landslide may occur in the period of heavy rainfall. Accordingly, this study emphasizes on shallow translational landslide considering the infinite slope model for the infiltration in unsaturated soils. The infinite slope model is used to analyze the shallow landslide in which the sliding surface along the failure plane [6]. The safety factor is defined as [7] FS = c + zγ t cos 2 (β) − hγ w cos 2 (β)] tan φ zγ t cos(β) sin(β) , where FS denotes safety factor, c denotes effective cohesiveness, φ denotes effective friction angle, z denotes soil thickness, γ w and γ t denote unit weight of water and soil, respectively. Fredlund et al. proposed a linear shear strength equation for the unsaturated soil as follows [26]: where τ f denotes the shear strength of an unsaturated soil, (σ n − u a ) denotes the net normal stress on the plane of failure at failure, (u a − u w ) denotes the matric suction of the soil on the plane of failure, and φ b denotes the angle of shearing resistance with respect to matric suction. According to previous studies [27], results show nonlinear shear strength behavior when experimental test are performed over a wide range of suction. Soils that are resistant to desaturation essentially exhibits a linear shear strength behavior over a relatively large range of soil suction. The above equation can be adopted to describe the shear strength behavior with respect to specified soil suction range. The angle of shearing resistance with respect to soil suction exhibit nonlinear shear strength behavior. Therefore, the relationship between the SWCC and the shear strength has been developed. The shear strength of unsaturated soils at any specified value of suction can be described as [28] where σ represents the total stress. The above revised shear strength equation is then substituted into Equation (28) for analyzing infiltration-induced landslides in unsaturated slope [27,35] FS = c zγ t cos(β) sin(β) Figure 2 illustrates the landslide stability for a variably saturated infinite slope. As demonstrated in Figure 2, Equations (28) and (31) are equations of the safety factor for the saturated and unsaturated slope, respectively. The schematic diagram for modeling infiltration-induced landslides in unsaturated soils is depicted in Figure 3.  Figure 2 illustrates the landslide stability for a variably saturated infinite slope. As demonstrated in Figure 2, Equations (28) and (31) are equations of the safety factor for the saturated and unsaturated slope, respectively. The schematic diagram for modeling infiltration-induced landslides in unsaturated soils is depicted in Figure 3.

One-Dimensional Steady-State Infiltration Problem
The one-dimensional steady-state flow in the unsaturated soil is described by the Richards equation as ∂ ∂z K r ∂h ∂z By maintaining the pressure head to be zero, the infiltration is remained at the ground surface, described as where L is 100 m which is the soil thickness. The boundary condition at the bottom side of the unsaturated soils is in dry condition, defined as follows.
where h d represents dry pressure head and h d is −10 5 m. To describe the nonlinearity of relative hydraulic conductivity, the scattered data produced by the Gardner model are utilized [17]. The type of the unsaturated soils is silty loam, where α g is 3.2 × 10 −4 . The schematic illustration of one-dimensional steady-state infiltration problem is illustrated in Figure 4.
( By maintaining the pressure head to be zero, the infiltration is remained at ground surface, described as where L is 100 m which is the soil thickness. The boundary condition at the bottom s of the unsaturated soils is in dry condition, defined as follows.
where d h represents dry pressure head and d h is 5 10 − m. To describe the nonlinea of relative hydraulic conductivity, the scattered data produced by the Gardner model utilized [17]. The type of the unsaturated soils is silty loam, where g α is 4 10 . T schematic illustration of one-dimensional steady-state infiltration problem is illustra in Figure 4. Figure 5 shows the predicted SWCC based on the RBF along with the scattered d by the Gardner model. It is clear that the predicted SWCC based on the RBF closely agr with the scattered data by the Gardner model. The relative error between the scatte data by the Gardner and RBF model is depicted in Figure 6. The maximum relative er of the predicted SWCC based on the RBF is in the order of 11 10 − . It is significant t excellent agreement is achieved and highly accurate SWCC can be obtained by the RB   Figure 5 shows the predicted SWCC based on the RBF along with the scattered data by the Gardner model. It is clear that the predicted SWCC based on the RBF closely agrees with the scattered data by the Gardner model. The relative error between the scattered data by the Gardner and RBF model is depicted in Figure 6. The maximum relative error of the predicted SWCC based on the RBF is in the order of 10 −11 . It is significant that excellent agreement is achieved and highly accurate SWCC can be obtained by the RBF.     We then compare the accuracy of the results with the following exact solution [36].
The numerical parameters including N = 4, δ = 1.5, v = −1, m = 0.01, and ε = 5 × 10 −6 are considered in this numerical implementation. We considered two boundary points, 34 interior points, and 34 center points. In this example, the SWCC using the RBF, Gardner and van Genuchten model are adopted. The pore size distribution parameter of the Gardner model is 3.2 × 10 −4 . The fitted parameters of the van Genuchten model are α v = 0.05 and n v = 1.5. Figure 7 illustrates the comparison of the computed pressure using the proposed RBF, Gardner, and van Genuchten model with the exact solution. From Figure 7, it is clear that the computed pressure head using the proposed RBF and the Gardner model agree well with the exact solution. However, the computed pressure head using the van Genuchten model may not exactly fit with the exact solution. The main reason is that the above exact solution adopted is basically developed by the Gardner model. This is the reason why excellent agreement is achieved and accurate numerical solutions can be obtained by using the Gardner model. Results obtained shows that the validity of the proposed RBF approach is achieved in the one-dimensional steady-state infiltration problem. We then compare the accuracy of the results with the following exact solution [36].
The numerical parameters including , and 6 10 5 − × = ε are considered in this numerical implementation. We considered two boundary points, 34 interior points, and 34 center points. In this example, the SWCC using the RBF, Gardner and van Genuchten model are adopted. The pore size distribution parameter of the Gardner model is . Figure 7 illustrates the comparison of the computed pressure using the proposed RBF, Gardner, and van Genuchten model with the exact solution. From Figure 7, it is clear that the computed pressure head using the proposed RBF and the Gardner model agree well with the exact solution. However, the computed pressure head using the van Genuchten model may not exactly fit with the exact solution. The main reason is that the above exact solution adopted is basically developed by the Gardner model. This is the reason why excellent agreement is achieved and accurate numerical solutions can be obtained by using the Gardner model. Results obtained shows that the validity of the proposed RBF approach is achieved in the one-dimensional steady-state infiltration problem.

The Green-Ampt Problem
The second example is the analyzing of the one-dimensional Green-Ampt problem [37], as defined in Equation (6). The initial condition is as

The Green-Ampt Problem
The second example is the analyzing of the one-dimensional Green-Ampt problem [37], as defined in Equation (6). The initial condition is as The boundary data are given as follows.
As depicted in Equation (38), the ponding on the ground surface is assigned. The unsaturated soil is in dry condition until the infiltration starts to gradually flow through the unsaturated soil.
The exact solution of the one-dimensional Green-Ampt problem is as [36] h(z, t) = 1 where µ k = (α 2 g /4 + λ 2 k )/c, c = α g (θ s − θ r )/K s , and λ k = πk/L. The soil thickness is 100 m. The type of the unsaturated soil is sand. To define the nonlinearity of relative hydraulic conductivity, the scattered data produced by the Gardner model are utilized, where α g = 1.6 × 10 −2 , θ s = 0.5, θ r = 0.1, K s = 1 m/h, and total time is 0.25 h. The schematic illustration of one-dimensional Green-Ampt problem is illustrated in Figure 8.   The number of the interior points, center points, and boundary points are considered to be 48, 48, and two. The numerical parameters including N = 4, δ = 1.5, v = −2.5 × 10 −3 , m = 0.01, ε = 5 × 10 −6 , and ∆t = 10 −4 are adopted. Figure 9 shows the predicted SWCC based on the RBF along with the scattered data by the Gardner model. It is clear that the predicted SWCC based on the RBF closely agrees with the scattered data by the Gardner model. The maximum relative errors of the predicted SWCC based on the RBF is in the order of 10 −11 . It is significant that excellent agreement is achieved and highly accurate SWCC can be obtained by the RBF. Figure 8. A schematic of one-dimensional Green-Ampt problem. Figure 9. Comparing effective saturation using RBF with scattered data by the Gardner model To clearly demonstrate the predicted pressure head, the computed results using the RBF, Gardner and van Genuchten model on different simulation time were compare with the exact solution. The fitted parameters of the van Genuchten model are α v = 0.05 and n v = 1.8. Figure 10 demonstrates the comparison of one-dimensional transient flow in unsaturated soil with the exact solution. It is obvious that the computed pressure head using the proposed RBF and the Gardner model agree well with the exact solution, as depicted in Figure 10. We compare the relative error of the RBF with that of the exact solution, as depicted in Figure 11. The maximum relative error is in the order of 10 −3 to 10 −6 .

Infiltration-Induced Landslide Problem
The modeling of infiltration-induced landslide in unsaturated soils is presented. The illustration of the unsaturated inclined slope is depicted in Figure 12. We consider the thickness of the unsaturated inclined slope and the slope angle to be 10 m and 35 degrees, respectively. Equation (6) is adopted for this example.    . Relative error of the computed results using the RBF. Figure 11. Relative error of the computed results using the RBF. [32]. The soil type is sandy soil. The parameter of pore size distribution is 10 6 . 1 × , the saturated water content is 0.50, the residual water content is 0.11, the saturated hydraulic conductivity is 1 m/h, the effective cohesiveness is 4.6 kPa, the effective friction angle is 30 degrees, the soil unit weight is 21.5 kN/m 3 , and total simulation time is 0.   The unsaturated soils are in dry condition until the infiltration starts to gradually flow through the unsaturated soil. The initial and boundary conditions are assigned as follows.
h(z, t = 0) = −10. (40) To depict the nonlinearity of unsaturated soils, the experimental data are utilized [32]. The soil type is sandy soil. The parameter of pore size distribution is 1.6 × 10 −2 , the saturated water content is 0.50, the residual water content is 0.11, the saturated hydraulic conductivity is 1 m/h, the effective cohesiveness is 4.6 kPa, the effective friction angle is 30 degrees, the soil unit weight is 21.5 kN/m 3 , and total simulation time is 0.1 h. The fitted parameters of the van Genuchten model are α v = 0.05 and n v = 2.12. There are 48 interior points, 48 center points, and two boundary points. The numerical parameters including N = 4, δ = 1.5, v = −2.5 × 10 −3 , m = 10 −2 , ε = 5 × 10 −6 , and ∆t = 10 −4 are utilized. Figure 13 shows the predicted SWCCs using the RBF fitting model, the Gardner model and the van Genuchten model comparing with the experimental data. It is significant that the predicted SWCC using the proposed RBF fitting model presents the best fitting with the MAE in the order of 10 −10 among those empirical Gardner model and van Genuchten model. We then evaluate the computed matric potential using the Gardner model, van Genuchten model, and the proposed RBF. Figure 14 demonstrates a plot of matric potential versus depth at different simulation time. From Figure 14, it is found that the matric potential obtained by the Gardner model, van Genuchten model, and the proposed RBF are different. The conventional Gardner and van Genuchten models have simple forms considering only one and three fitting parameters, respectively. Accordingly, these SWCC models may be only applicable for fitting the SWCC of certain soil textures. However, the advantage of the proposed RBF model lies in its applicability in fitting almost any soil types due to the fact that many weighting coefficients can be utilized in the RBF model for approximating the high nonlinearity of the SWCC. It appears that the matric potential obtained by the proposed RBF may depict the physical property of unsaturated soils.
To analyze the infiltration-induced landslides in unsaturated soils, the hydrological model considering the stability analysis of infinite slope is then carried out. To evaluate the safety factor for the unsaturated slope, Equations (28) and (31) are utilized. Figure 15 depicts the computed results of safety factor. As illustrated in Figure 15, we also found that the computed FS obtained by different methods are very different because Gardner or van Genuchten models may have limitations to fit the SWCC. Since excellent agreement and highly accurate SWCC using the proposed RBF model can be yielded (as depicted in Figure 13), the FS obtained by the proposed RBF may comprehensively consider the physical property of unsaturated soils. Furthermore, the FS obtained by the proposed RBF demonstrates that landslide would be triggered rapidly when FS is less than 1.0 and the lowest FS may occur at the z = 8 to 10 m during the infiltration, as shown in Figure 15. From Figures 14 and 15, it is found that the slope stability is highly related to the fluctuation of matric potential in unsaturated soils during the infiltration process.
Water 2022, 14, x FOR PEER REVIEW 18 of 22 cant that the predicted SWCC using the proposed RBF fitting model presents the best fitting with the MAE in the order of 10 10 − among those empirical Gardner model and van Genuchten model. We then evaluate the computed matric potential using the Gardner model, van Genuchten model, and the proposed RBF. Figure 14 demonstrates a plot of matric potential versus depth at different simulation time. From Figure 14, it is found that the matric potential obtained by the Gardner model, van Genuchten model, and the proposed RBF are different. The conventional Gardner and van Genuchten models have simple forms considering only one and three fitting parameters, respectively. Accordingly, these SWCC models may be only applicable for fitting the SWCC of certain soil textures. However, the advantage of the proposed RBF model lies in its applicability in fitting almost any soil types due to the fact that many weighting coefficients can be utilized in the RBF model for approximating the high nonlinearity of the SWCC. It appears that the matric potential obtained by the proposed RBF may depict the physical property of unsaturated soils.  To analyze the infiltration-induced landslides in unsaturated soils, the hydrological model considering the stability analysis of infinite slope is then carried out. To evaluate the safety factor for the unsaturated slope, Equations (28) and (31) are utilized. Figure 15 depicts the computed results of safety factor. As illustrated in Figure 15, we also found that the computed FS obtained by different methods are very different because Gardner or van Genuchten models may have limitations to fit the SWCC. Since excellent agreement and highly accurate SWCC using the proposed RBF model can be yielded (as depicted in Figure 13), the FS obtained by the proposed RBF may comprehensively consider the physical property of unsaturated soils. Furthermore, the FS obtained by the proposed RBF demonstrates that landslide would be triggered rapidly when FS is less than 1.0 and the lowest FS may occur at the z = 8 to 10 m during the infiltration, as shown in Figure 15. From Figures 14 and 15, it is found that the slope stability is highly related to the fluctuation of matric potential in unsaturated soils during the infiltration process.

Conclusions
A novel RBF approach for infiltration-induced landslides in unsaturated soils is presented. The study emphasizes on the importance of infiltration processes to predict landslide hazard. Significant findings are concluded as follows: (1) This article is a leading attempt to model infiltration-induced landslides in unsaturated soils utilizing the RBF method. Since the soil water characteristic curve is a primary factor defining the nonlinearity of unsaturated soils, we first utilize the RBF for Figure 15. Computed FS.

Conclusions
A novel RBF approach for infiltration-induced landslides in unsaturated soils is presented. The study emphasizes on the importance of infiltration processes to predict landslide hazard. Significant findings are concluded as follows: (1) This article is a leading attempt to model infiltration-induced landslides in unsaturated soils utilizing the RBF method. Since the soil water characteristic curve is a primary factor defining the nonlinearity of unsaturated soils, we first utilize the RBF for curve fitting to build the representation of the SWCC. Results demonstrate that the RBF seems to be advantageous for reconstructing the SWCC of glacial till and sand soil with better estimation of the relationship than conventional parametric Gardner and van Genuchten models. (2) For solving the nonlinear Richards equation, we developed a novel work using the RBF approach with the infiltration boundary conditions. The fictitious time integration method is adopted in the RBF approach for tackling the nonlinearity of the Richards equation. We propose the novel RBF approach for reconstructing the SWCC to directly solving the nonlinear Richards equation such that the conventional parametric Gardner model and van Genuchten model are not required in the proposed approach. (3) The hydrological model considering the stability analysis of infinite slope is conducted to analyze the infiltration-induced landslides with emphasis on the behavior of unsaturated soils. Additionally, the proposed RBF method for the nonlinear Richards equation in an inclined slope can be widely applied for disparate types of unsaturated soils. We also found that the stability of landslides is highly related to the variation of matric potential in unsaturated soils during the infiltration process.