Evaluation of Ground Motion Amplification Effects in Slope Topography Induced by the Arbitrary Directions of Seismic Waves

The incident directions of seismic waves can change the ground motions of slope topography. To elaborate on the influences of the directions of seismic waves, a dynamic analysis of the slope topography was performed. Seismic waves were input using an equivalent nodal force method combined with a viscous-spring artificial boundary. The amplification of ground motions in double-faced slope topographies was discussed by varying the angles of incidence. Meanwhile, the components of seismic waves (P waves and SV waves), slope materials and slope geometries were all investigated with various incident earthquake waves. The results indicated that the pattern of the amplification of SV waves was stronger than that of P waves in the slope topography, especially in the greater incident angels of the incident waves. Soft materials intensely aggravate the acceleration amplification, and more scattered waves are produced under oblique incident earthquake waves. The variations in the acceleration amplification ratios on the slope crest were much more complicated at oblique incident waves, and the ground motions were underestimated by considering only the vertical incident waves. Therefore, in the evaluation of ground motion amplification of the slope topography, it is extremely important to consider the direction of incident waves.


Introduction
Earthquake-induced landslides are among the most hazardous secondary effects of strong seismicity in mountainous regions. Several historical cases highlight the potentially devastating consequences in the last 30 years [1][2][3][4][5][6] because nearly one-third of deaths are caused by earthquake-induced landslides among all disasters [7][8][9][10]. Many researchers have indicated that topographic effects have possible interactions between landslide mechanisms and triggering conditions [11][12][13][14][15][16]. Thus, the topographic effects, represented by the ground motions on the slope model, should be better understood.
Generally, convex topographies such as mountains, slopes and individual ridges lead to intense aggravation of the seismic responses irregularly along the ground surface [17]. Many acceleration-time histories recorded by the seismometers near the epicenter all proved that the peak acceleration at the slope crest showed intense amplification based on instrumental data from strong earthquakes such as the 1987 California earthquake, the 1994 Northridge earthquake, the 1995 Egion earthquake and the 2008 Wenchuan earthquake [12,[18][19][20]. In addition, statistical analysis of the seismological stations in Iran and Israel showed that the acceleration amplification ratios with respect to the slope crest and ground could reach as much as four times [21,22]. These measured records confirmed that the slope topography modified the seismic ground motions. The seismic waves were analyzed by two components of the seismic waves (considering only the body waves), that is, P waves and SV waves. The influence of the wave components was analyzed using the same model size. The materials of the slope were assumed to be isotropic, and the medium was linearly elastic. To investigate the influence of materials in the evaluation of the topographic effects, the material of the foundation was fixed, whereas the materials of the slope varied according to the variations in the materials in the foundations. The specific parameters of the materials are presented in Table 1.

Establishment of the Artificial Boundary and the Input Method
In a dynamic analysis of the numerical simulation, the study area is usually truncated from an infinite region with an artificial boundary used to absorb the scattered waves. In this study, a viscous-spring artificial boundary was adopted [54,55]. The elastic springs and dampers were established to implement the viscous-spring artificial boundary on the boundary nodes, as shown in Figure 2. The coefficients of the elastic spring (K) and damper (C) are defined as follows. The seismic waves were analyzed by two components of the seismic waves (considering only the body waves), that is, P waves and SV waves. The influence of the wave components was analyzed using the same model size. The materials of the slope were assumed to be isotropic, and the medium was linearly elastic. To investigate the influence of materials in the evaluation of the topographic effects, the material of the foundation was fixed, whereas the materials of the slope varied according to the variations in the materials in the foundations. The specific parameters of the materials are presented in Table 1.

Establishment of the Artificial Boundary and the Input Method
In a dynamic analysis of the numerical simulation, the study area is usually truncated from an infinite region with an artificial boundary used to absorb the scattered waves. In this study, a viscous-spring artificial boundary was adopted [54,55]. The elastic springs and dampers were established to implement the viscous-spring artificial boundary on the boundary nodes, as shown in Figure 2. The coefficients of the elastic spring (K) and damper (C) are defined as follows. where A and B are the correction coefficients; the good values of the coefficients are 0.8 and 1.1, respectively [56]; ρ is the mass density; R is the distance from the wave source to the boundary; c p and c s are the velocities of the compression wave and shear wave in the medium, respectively; G is the shear modulus; and subscripts T and N indicate the tangential and normal directions, respectively. A l represents the influence area at each node, for example, at node l, A l = (A 1 + A 2 )/2, as depicted in Figure 2.  (2) where A and B are the correction coefficients; the good values of the coefficients are 0.8 and 1.1, respectively [56];  is the mass density; R is the distance from the wave source to the boundary; cp and cs are the velocities of the compression wave and shear wave in the medium, respectively; G is the shear modulus; and subscripts T and N indicate the tangential and normal directions, respectively. Al represents the influence area at each node, for example, at node l, Al = (A1 + A2)/2, as depicted in Figure 2.
In case of the distribution of the wavefield on the artificial boundary, the displacement and stress on the artificial boundary can be divided into free-field motions (denoted by superscripts f) and scattered-field motions denoted by superscripts s), which are illustrated as follows: fs =+ u u u fs =+ σ σ σ (4) The equation of stress (specifically for the scattered-field motions) on the viscous-spring artificial boundary is given as: Equation (6) can be considered seismic input and non-radiation, and the right side of the equation represents the equivalent nodal forces. When the two directions of the viscous-spring artificial boundary (normal direction and tangential direction) are considered, the equivalent nodal force can be expressed as follows: In case of the distribution of the wavefield on the artificial boundary, the displacement and stress on the artificial boundary can be divided into free-field motions (denoted by superscripts f ) and scattered-field motions denoted by superscripts s), which are illustrated as follows: The equation of stress (specifically for the scattered-field motions) on the viscousspring artificial boundary is given as: Substituting Equations (4) and (5) into Equation (3), the motion equation on the artificial boundary is expressed as follows: Equation (6) can be considered seismic input and non-radiation, and the right side of the equation represents the equivalent nodal forces. When the two directions of the viscousspring artificial boundary (normal direction and tangential direction) are considered, the equivalent nodal force can be expressed as follows: where subscript i denotes the direction and subscript l denotes the node. K li and C li are the coefficients of the elastic spring and damper, which can be calculated using Equations (1) and (2) Based on the artificial boundary, the other key point in the numerical simulation of a dynamic problem is the implementation of seismic waves in finite element models. An efficient and high-precision method is to convert the seismic records, that is, acceleration records, together with the coefficients of the artificial boundary into equivalent nodal forces. Figure 3 presents a plane study area that was subjected to oblique incident SV waves. As shown in Figure 3, each truncated boundary was subjected to three waves that involved one incident wave and two reflected waves. When node A, one of the joints on the left truncated boundary, were considered, the incident seismic waves of node A reached the boundary directly (incident SV waves) or through reflection after first reaching the ground surface (reflected SV waves and reflected P waves). The displacement and velocity of the incident SV waves are denoted as u s (t) and . u s (t), respectively.
where subscript i denotes the direction and subscript l denotes the node. Kli and Cli are the coefficients of the elastic spring and damper, which can be calculated using Equations (1) and (2), respectively; f li σ , f li u and u f li represent the tensors of the stress, displacement and velocity, respectively.

SV Waves
Based on the artificial boundary, the other key point in the numerical simulation of a dynamic problem is the implementation of seismic waves in finite element models. An efficient and high-precision method is to convert the seismic records, that is, acceleration records, together with the coefficients of the artificial boundary into equivalent nodal forces. Figure 3 presents a plane study area that was subjected to oblique incident SV waves. As shown in Figure 3, each truncated boundary was subjected to three waves that involved one incident wave and two reflected waves. When node A, one of the joints on the left truncated boundary, were considered, the incident seismic waves of node A reached the boundary directly (incident SV waves) or through reflection after first reaching the ground surface (reflected SV waves and reflected P waves). The displacement and velocity of the incident SV waves are denoted as  Consequently, the free-field motion and the associated stress (in terms of x and y components) can be induced at the given point A (x0, y0) as follows:  Consequently, the free-field motion and the associated stress (in terms of x and y components) can be induced at the given point A (x 0 , y 0 ) as follows: where A 1 and A 2 are the amplitude ratios (denoted as the reflected SV waves and reflected P waves to the incident SV waves), ∆t is the delay time of waves from the wavefront at t = 0 to the boundary node (node A), α and β are the angles between the vertical direction and reflected SV waves, and the reflected P waves, respectively. A 1 , A 2 and β can be expressed as follows: where c s and c p are the wave velocities of P waves and SV waves, respectively. The delay times ∆t 1 , ∆t 2 , ∆t 3 , and variables S 1 to S 6 are all boundary-dependent, which can be depicted as follows: On the bottom boundary, they became: On the right boundary, the stresses were the same but in opposite directions to those on the left boundary, and for the displacement, the delay times were ∆t 1 , ∆t 2 , and ∆t 3 , each of them increased by L x sin α/c s for the additional travel distance. Figure 4 presents a plane study area that was subjected to oblique incident P waves. Similar to the oblique incident SV waves, each truncated boundary was subjected to three waves that involved one incident wave and two reflected waves, as shown in Figure 4. When node B, one of the joints on the left truncated boundary, were considered, the incident seismic waves of node B reached the boundary directly (incident P waves) or through reflection after first reaching the ground surface (reflected P waves and reflected SV waves). The displacement and velocity of the incident P waves are denoted as u p (t) and . u p (t), respectively. On the right boundary, the stresses were the same but in opposite directions to those on the left boundary, and for the displacement, the delay times were  Figure 4 presents a plane study area that was subjected to oblique incident P waves. Similar to the oblique incident SV waves, each truncated boundary was subjected to three waves that involved one incident wave and two reflected waves, as shown in Figure 4. When node B, one of the joints on the left truncated boundary, were considered, the incident seismic waves of node B reached the boundary directly (incident P waves) or through reflection after first reaching the ground surface (reflected P waves and reflected SV waves). The displacement and velocity of the incident P waves are denoted as  Consequently, the free-field motion and associated stress (in terms of their x and y components) can be induced at the given point B (x1, y1) as follows:   Consequently, the free-field motion and associated stress (in terms of their x and y components) can be induced at the given point B (x 1 , y 1 ) as follows:

P Waves
where A 3 and A 4 are the amplitude ratios (denoted as the reflected P waves and reflected SV waves to the incident P waves). ∆t denotes the delay time of waves from the wavefront at t = 0 to the boundary node (node B). α and β are the angles between the vertical direction and reflected P waves together with the reflected SV waves, respectively. A 3 , A 4 and β can be expressed as follows: where c s and c p are the wave velocities of P waves and SV waves, respectively. The delay times ∆t 4 , ∆t 5 , ∆t 6 , and the variables S 7 to S 12 are all boundary-dependent, which can be depicted as follows: On the bottom boundary, they became: On the right boundary, the stresses were the same but in directions opposite to those on the left boundary, and for the displacement, the delay times were ∆t 4 , ∆t 5 , and ∆t 6 " each of which increased by L x sin α/c p for the additional travel distance.

Verification
The incident seismic waves were implemented using the ABAQUS software, together with a self-developed MATLAB program. The specific implementation is described in the flowchart illustrated in Figure 5. Two numerical test examples were used to verify the validity of the oblique incident seismic waves. One was the propagating process of oblique incident waves in a truncated region, which was considered to explore the accuracy of the input method in a semi-infinite field. The other was the response of a semicircular canyon input by oblique incident waves, which was aimed at determining the accuracy of the simulation in the topography. process of oblique incident waves in a truncated region, which was considered to explore the accuracy of the input method in a semi-infinite field. The other was the response of a semicircular canyon input by oblique incident waves, which was aimed at determining the accuracy of the simulation in the topography.

Test Example 1
A truncated rectangular domain was adopted to simulate the propagation process of oblique incident waves in a semi-infinite field, as depicted in Figure 6a. The sizes of the truncated regions were as follows: Lx = 800 m and Ly = 400 m. The region was sufficiently large to detect the propagation of all the incident and reflected waves. The incident angles of the P waves and SV waves were assumed to be 30° and 18°, respectively. are the observation points. An impulse was adopted as an incident seismic wave, as shown in Figure 6b. The impulse equation is defined as follows [57]:

Test Example 1
A truncated rectangular domain was adopted to simulate the propagation process of oblique incident waves in a semi-infinite field, as depicted in Figure 6a. The sizes of the truncated regions were as follows: L x = 800 m and L y = 400 m. The region was sufficiently large to detect the propagation of all the incident and reflected waves. The incident angles of the P waves and SV waves were assumed to be 30 • and 18 • , respectively.
Energies 2021, 14, x FOR PEER REVIEW 9 of 29 process of oblique incident waves in a truncated region, which was considered to explore the accuracy of the input method in a semi-infinite field. The other was the response of a semicircular canyon input by oblique incident waves, which was aimed at determining the accuracy of the simulation in the topography.

Test Example 1
A truncated rectangular domain was adopted to simulate the propagation process of oblique incident waves in a semi-infinite field, as depicted in Figure 6a. The sizes of the truncated regions were as follows: Lx = 800 m and Ly = 400 m. The region was sufficiently large to detect the propagation of all the incident and reflected waves. The incident angles of the P waves and SV waves were assumed to be 30° and 18°, respectively. are the observation points. An impulse was adopted as an incident seismic wave, as shown in Figure 6b. The impulse equation is defined as follows [57]: The material parameters of the domain were as follows: Poisson's ratio = 0.25, elastic modulus = 1.25 Gpa and mass density, ρ = 2000 kg/m 3 . Points A (400, 400) and B (400, 0) are the observation points. An impulse was adopted as an incident seismic wave, as shown in Figure 6b. The impulse equation is defined as follows [57]: where G = (t/T 0 ) 3 H(t/T 0 ) H(t) is the Heaviside function, P 0 is the peak value of an impulse, P 0 is set to 1.0 m; and T 0 is the acting time of the impulse, where T 0 = 0.25 s. Figure 7 shows the displacement contours with the time of incident SV waves. As illustrated in the figures, the propagation processes of the incident waves and reflected waves at the ground surface were simulated effectively in a semi-infinite field. In addition, Energies 2021, 14, 6744 9 of 27 the displacement time histories of the observation points (points A and B) are shown in Figure 8. As shown in the figures, the numerical results proposed by the input method were in good agreement with the analytical solutions, where the analytical solutions were calculated using the elastic wave propagation theory [58,59].
is the Heaviside function, P0 is the peak value of an impulse, P0 is set to 1.0 m; and T0 is the acting time of the impulse, where T0 = 0.25 s. Figure 7 shows the displacement contours with the time of incident SV waves. As illustrated in the figures, the propagation processes of the incident waves and reflected waves at the ground surface were simulated effectively in a semi-infinite field. In addition, the displacement time histories of the observation points (points A and B) are shown in Figure 8. As shown in the figures, the numerical results proposed by the input method were in good agreement with the analytical solutions, where the analytical solutions were calculated using the elastic wave propagation theory [58,59].
is the Heaviside function, P0 is the peak value of an impulse, P0 is set to 1.0 m; and T0 is the acting time of the impulse, where T0 = 0.25 s. Figure 7 shows the displacement contours with the time of incident SV waves. As illustrated in the figures, the propagation processes of the incident waves and reflected waves at the ground surface were simulated effectively in a semi-infinite field. In addition, the displacement time histories of the observation points (points A and B) are shown in Figure 8. As shown in the figures, the numerical results proposed by the input method were in good agreement with the analytical solutions, where the analytical solutions were calculated using the elastic wave propagation theory [58,59]. The displacement contours with the time of P waves in Figure 9 show the propagation processes of the incident waves and reflected waves at the ground surface. Meanwhile, the displacement contours and the displacement time histories of the observation points (points A and B) in Figure 10 both confirm that the input method was effective and precise in the semi-infinite field. Therefore, the oblique incident SV waves and P waves in a semi-infinite field could be simulated using the proposed input method, and the The displacement contours with the time of P waves in Figure 9 show the propagation processes of the incident waves and reflected waves at the ground surface. Meanwhile, the displacement contours and the displacement time histories of the observation points (points A and B) in Figure 10 both confirm that the input method was effective and precise in the semi-infinite field. Therefore, the oblique incident SV waves and P waves in a semi-infinite field could be simulated using the proposed input method, and the ability of the viscous-spring artificial boundary to absorb the scattered waves was also verified. The displacement contours with the time of P waves in Figure 9 show the propagation processes of the incident waves and reflected waves at the ground surface. Meanwhile, the displacement contours and the displacement time histories of the observation points (points A and B) in Figure 10 both confirm that the input method was effective and precise in the semi-infinite field. Therefore, the oblique incident SV waves and P waves in a semi-infinite field could be simulated using the proposed input method, and the ability of the viscous-spring artificial boundary to absorb the scattered waves was also verified.

Test Example 2
A truncated rectangular domain with a semicircular canyon is shown in Figure 11a. The truncated region was 0 ≤ x ≤ 200 m, 0 ≤ y ≤ 200 m and the radius of the canyon was 25 m. The material parameters of the ground were as follows: Poisson's ratio = 0.3, elastic modulus = 208 MPa and mass density, ρ = 2000 kg/m 3 . A Ricker wavelet was adopted as the incident wave. The acceleration-time history of the wave is illustrated in Figure 11b. The time history of the Ricker wavelet was defined as follows: with an excitation frequency f 0 of 4.0 Hz. The incident angles had values of 25 • under SV waves and 30 • under P waves. Figures 12 and 13 show the horizontal displacement amplification ratios U h and vertical displacement amplification ratios U v of the ground surface on the canyon with incident SV waves and incident P waves, respectively. U h is defined as the ratio of the horizontal displacement on the ground surface to the input horizontal displacement, and U v is the ratio of the vertical displacement. The numerical results presented in Figures 12 and 13 were in good agreement with the analytical solutions proposed by Wong [60]. Therefore, in the topographic analysis, the proposed input method could effectively simulate the propagation process of oblique incident seismic waves, and the viscous-spring artificial boundary could absorb the scattered waves perfectly because of the topography. The incident angles had values of 25° under SV waves and 30° under P waves. Figures 12 and 13 show the horizontal displacement amplification ratios Uh and vertical displacement amplification ratios Uv of the ground surface on the canyon with incident SV waves and incident P waves, respectively. Uh is defined as the ratio of the horizontal displacement on the ground surface to the input horizontal displacement, and Uv is the ratio of the vertical displacement. The numerical results presented in Figures 12 and 13 were in good agreement with the analytical solutions proposed by Wong [60]. Therefore, in the topographic analysis, the proposed input method could effectively simulate the propagation process of oblique incident seismic waves, and the viscous-spring artificial boundary could absorb the scattered waves perfectly because of the topography. show the horizontal displacement amplification ratios Uh and vertical displacement amplification ratios Uv of the ground surface on the canyon with incident SV waves and incident P waves, respectively. Uh is defined as the ratio of the horizontal displacement on the ground surface to the input horizontal displacement, and Uv is the ratio of the vertical displacement. The numerical results presented in Figures 12 and 13 were in good agreement with the analytical solutions proposed by Wong [60]. Therefore, in the topographic analysis, the proposed input method could effectively simulate the propagation process of oblique incident seismic waves, and the viscous-spring artificial boundary could absorb the scattered waves perfectly because of the topography.

Discussion and Results
A schematic of the numerical model is illustrated in Figure 1. The mesh type of this numerical model is quadrilateral element. A Ricker wavelet with a predominant frequency of 4.0 Hz was employed as the incident wave to investigate the wave patterns, materials and slope geometries based on the seismic responses. The acceleration-time history of the incident wave, which was input from the left bottom of the model, is illustrated in Figure 10b. To analyze the topographic effects, the acceleration amplification ratio was applied, which is the ratio of the peak acceleration of the reflected waves at

Discussion and Results
A schematic of the numerical model is illustrated in Figure 1. The mesh type of this numerical model is quadrilateral element. A Ricker wavelet with a predominant frequency of 4.0 Hz was employed as the incident wave to investigate the wave patterns, materials and slope geometries based on the seismic responses. The acceleration-time history of the incident wave, which was input from the left bottom of the model, is illustrated in Figure 10b. To analyze the topographic effects, the acceleration amplification ratio was applied, which is the ratio of the peak acceleration of the reflected waves at

Discussion and Results
A schematic of the numerical model is illustrated in Figure 1. The mesh type of this numerical model is quadrilateral element. A Ricker wavelet with a predominant frequency of 4.0 Hz was employed as the incident wave to investigate the wave patterns, materials and slope geometries based on the seismic responses. The acceleration-time history of the incident wave, which was input from the left bottom of the model, is illustrated in Figure 10b. To analyze the topographic effects, the acceleration amplification ratio was applied, which is the ratio of the peak acceleration of the reflected waves at each point on the surface to that of the incident seismic waves [36,45]. The horizontal and vertical acceleration amplification ratios r h and r v are defined in Equations (20) and (21), respectively. For a given wave, r h and r v are the values of the acceleration amplification ratio r for the horizontal and vertical components, respectively. r = r h 2 + r v 2 defines zones of the net amplification of ground acceleration with respect to the energy carried by the incoming waves.
where max(|a h |) and max(|a v |) are the peak horizontal (h) and peak vertical (v) acceleration at the observation points along the slope ridges and slope crests, respectively; and a h,input and a v,input are the horizontal (h) and vertical (v) acceleration, respectively. In this study, the wave patterns were discussed based on the P waves and the SV waves owing to the specific focus on the effects of body waves. In addition, the impact of the materials was investigated in terms of the relative hardness and softness between the slopes and foundations. The shear waves of different materials on the slope topographies are presented in Table 1. Furthermore, the geometries of slopes with varying slope heights, widths and inclinations are discussed. The slope height H is a normalized height that is related to the wavelength λ [29,36], and it varies among 0.2λ, 0.5λ, 1.

Effects of Wave Patterns
Wave patterns, that is, P waves and SV waves, were discussed based on three different slope topographies with varied angles of incident waves. Three models with different slope widths, W = 0 m, 100 m and 400 m, were built to evaluate the effect of wave patterns on ground motion amplification. The original model had a height of 100 m and an inclination of 1.0 (45 • ); meanwhile, the materials of the slope were regarded as homogeneous. The acceleration amplification ratios of the P waves and SV waves at different seismic wave inclinations are summarized in Figure 14.
The acceleration amplification ratios of SV waves were larger than those of P waves at a specific slope topography. The slope width became wider and the amplification ratios gradually decreased for each slope pattern. With an increase in the angle of incidence, the acceleration amplification effects moved away from the epicenter and the maximum acceleration amplification ratios were obtained when the incident angles fluctuated between 20 • and 30 • , which were especially clear under SV waves. Thus, the pattern of amplification of the SV waves was verified to be stronger in the slope topography, and the evolution was the same as that of the findings of Meunier [35] (depicted as the slope configurations of Figure 14a). The acceleration amplification ratios of SV waves were larger than those of P waves at a specific slope topography. The slope width became wider and the amplification ratios gradually decreased for each slope pattern. With an increase in the angle of incidence, the acceleration amplification effects moved away from the epicenter and the maximum acceleration amplification ratios were obtained when the incident angles fluctuated between 20° and 30°, which were especially clear under SV waves. Thus, the pattern of amplification of the SV waves was verified to be stronger in the slope topography, and the evolution was the same as that of the findings of Meunier [35] (depicted as the slope configurations of Figure 14a).

Effects of Materials
In this section, materials are discussed based on the specific slope topography with varying directions of incident waves. Three models, wherein the slope is soft whereas the foundation is hard, the slope and foundation are both the same, and the slope is hard whereas the foundation is soft, were built to investigate the effects of the materials on the ground motion amplification of slope topography. The investigated model was 100 m high, 100 m wide and 1.0 (45°) in inclination. The so-called softness and hardness of the materials in the slopes or foundations in Figure 15 is in contrast to that of the materials in the foundation (the foundation materials are fixed); that is, the materials of the slope changed from soft to hard (from left to right in Figure 15). The acceleration amplification ratios of the different materials at different incident angles are summarized in Figure 15.

Effects of Materials
In this section, materials are discussed based on the specific slope topography with varying directions of incident waves. Three models, wherein the slope is soft whereas the foundation is hard, the slope and foundation are both the same, and the slope is hard whereas the foundation is soft, were built to investigate the effects of the materials on the ground motion amplification of slope topography. The investigated model was 100 m high, 100 m wide and 1.0 (45 • ) in inclination. The so-called softness and hardness of the materials in the slopes or foundations in Figure 15 is in contrast to that of the materials in the foundation (the foundation materials are fixed); that is, the materials of the slope changed from soft to hard (from left to right in Figure 15). The acceleration amplification ratios of the different materials at different incident angles are summarized in Figure 15. In Figure 15, the acceleration amplification ratios of the SV waves were larger than those of the P waves for each material. The ratios gradually decreased as the materials in the slope varied from soft to hard; that is, the soft materials intensified the amplification whereas the hard materials alleviated the amplification in the slope topography. In addition, the maximum acceleration amplification ratios were acquired when the incident  In Figure 15, the acceleration amplification ratios of the SV waves were larger than those of the P waves for each material. The ratios gradually decreased as the materials in the slope varied from soft to hard; that is, the soft materials intensified the amplification whereas the hard materials alleviated the amplification in the slope topography. In addition, the maximum acceleration amplification ratios were acquired when the incident angles fluctuated between 10 • and 25 • under SV waves incidence and fell between 0 • and 10 • under P waves incidence. Consequently, the slope topographies that covered soft materials were more harmful during earthquake events, and the maximum ground motions of slope topographies were usually in a certain direction of oblique incident seismic waves.
The synthetic accelerograms of the horizontal and vertical components of SV waves and P waves with incident angles of 0 • and 30 • are presented in Figures 16 and 17, respectively. In the figures, "soft-hard" was the abbreviation of "the slope is soft and the foundation is hard," which corresponded to the situation of "a" in Figure 15; "same" referred to the situation of "b" in Figure 15, and "hard-soft" was marked as the situation of "c" in Figure 15.
As illustrated in the figures, the soft materials led to intense amplification of the ground motions, whereas the hard materials caused de-amplification. The scattered waves under the soft materials were abundant, whereas they were weakened in the hard materials. These evolutions were observed in the horizontal components of the SV waves and the vertical components of the P waves. In addition, compared with the scattered waves at 0 • incidence and 30 • incidence in Figures 16 and 17, the scattered waves were much more abundant at the oblique incident waves owing to the added refraction effects by the oblique directions of incident waves. In summary, the intense ground motions with soft materials in slope topography could be explained by the abundantly scattered waves, and this phenomenon would aggravate at oblique incident waves. In Figure 15, the acceleration amplification ratios of the SV waves were larger than those of the P waves for each material. The ratios gradually decreased as the materials in the slope varied from soft to hard; that is, the soft materials intensified the amplification whereas the hard materials alleviated the amplification in the slope topography. In addition, the maximum acceleration amplification ratios were acquired when the incident angles fluctuated between 10° and 25° under SV waves incidence and fell between 0° and 10° under P waves incidence. Consequently, the slope topographies that covered soft materials were more harmful during earthquake events, and the maximum ground motions of slope topographies were usually in a certain direction of oblique incident seismic waves.
The synthetic accelerograms of the horizontal and vertical components of SV waves and P waves with incident angles of 0° and 30°are presented in Figures 16 and 17, respectively. In the figures, "soft-hard" was the abbreviation of "the slope is soft and the foundation is hard," which corresponded to the situation of "a" in Figure 15; "same" referred to the situation of "b" in Figure 15, and "hard-soft" was marked as the situation of "c" in Figure 15. As illustrated in the figures, the soft materials led to intense amplification of the ground motions, whereas the hard materials caused de-amplification. The scattered waves under the soft materials were abundant, whereas they were weakened in the hard materials. These evolutions were observed in the horizontal components of the SV waves and the vertical components of the P waves. In addition, compared with the scattered waves at 0° incidence and 30° incidence in Figures 16 and 17, the scattered waves were much more abundant at the oblique incident waves owing to the added refraction effects by the oblique directions of incident waves. In summary, the intense ground motions with soft materials in slope topography could be explained by the abundantly scattered waves, and this phenomenon would aggravate at oblique incident waves.

Effects of Slope Geometries
According to the analysis of the impacts of wave patterns, the pattern of P waves was weak because of the amplification ratios of the ground motions, whereas that of SV waves was strong. Thus, the effects of the slope geometries were discussed only based on the SV waves rather than the P waves. In addition, the slope geometry changed the wave paths at various directions of wave incidence; thus, the materials of the slope topographies were assumed to be homogenous to eliminate the impacts of non-geometric factors. Consequently, this section, discusses the investigation of the effects of the slope geometries using homogeneous materials under SV incident waves.

Effects of Slope Geometries
According to the analysis of the impacts of wave patterns, the pattern of P waves was weak because of the amplification ratios of the ground motions, whereas that of SV waves was strong. Thus, the effects of the slope geometries were discussed only based on the SV waves rather than the P waves. In addition, the slope geometry changed the wave paths at various directions of wave incidence; thus, the materials of the slope topographies were assumed to be homogenous to eliminate the impacts of non-geometric factors. Consequently, this section, discusses the investigation of the effects of the slope geometries using homogeneous materials under SV incident waves.

Influence of Slope Height H
In this subsection, the slope height varied from 20 m, 50 m, 100 m, to 200 m, whereas the other geometrical parameters (slope widths and slope inclinations) remained the same. The acceleration amplification ratios along the slope surface are illustrated in Figure 18 with varying angles of incident waves. In Figure 18, different distributions of ground motions are presented with varying slope heights under SV incident waves. The acceleration amplification ratios rose significantly with an increase in the slope height from 0.2  to 2.0  . The area of the maximum acceleration amplification ratios at each height moved to the center of the slope topography with an increase in slope height. In addition, the maximum acceleration amplification ratios could be obtained when the incident angles fell between 25° and 30° under a slope height of 0.2  , whereas the maximum acceleration amplification ratios In Figure 18, different distributions of ground motions are presented with varying slope heights under SV incident waves. The acceleration amplification ratios rose signifi-cantly with an increase in the slope height from 0.2λ to 2.0λ. The area of the maximum acceleration amplification ratios at each height moved to the center of the slope topography with an increase in slope height. In addition, the maximum acceleration amplification ratios could be obtained when the incident angles fell between 25 • and 30 • under a slope height of 0.2λ, whereas the maximum acceleration amplification ratios could be obtained when the incident angles fluctuated between 5 • and 10 • under a slope height of 2.0λ. Thus, the higher the slope heights, the closer the maximum acceleration amplification ratios were to the smaller inclinations of wave incidence. In summary, different slope heights influenced the wave aggregation at the crest; the higher the model, the larger the acceleration amplification ratios. The maximum ground motions were obtained in the oblique direction of the incident waves.
Moreover, the observation points (slope toe and slope crest) are presented with horizontal and vertical components in Figures 19 and 20, respectively, to analyze the comprehensive relations between the incident directions and slope heights. In Figure 18, different distributions of ground motions are presented with varying slope heights under SV incident waves. The acceleration amplification ratios rose significantly with an increase in the slope height from 0.2  to 2.0  . The area of the maximum acceleration amplification ratios at each height moved to the center of the slope topography with an increase in slope height. In addition, the maximum acceleration amplification ratios could be obtained when the incident angles fell between 25° and 30° under a slope height of 0.2  , whereas the maximum acceleration amplification ratios could be obtained when the incident angles fluctuated between 5° and 10° under a slope height of 2.0  . Thus, the higher the slope heights, the closer the maximum acceleration amplification ratios were to the smaller inclinations of wave incidence. In summary, different slope heights influenced the wave aggregation at the crest; the higher the model, the larger the acceleration amplification ratios. The maximum ground motions were obtained in the oblique direction of the incident waves.
Moreover, the observation points (slope toe and slope crest) are presented with horizontal and vertical components in Figures 19 and 20, respectively, to analyze the comprehensive relations between the incident directions and slope heights.  As shown in Figures 19 and 20, the variations in the amplification ratios at the slope toe were gentle. The horizontal amplification ratios at the left toe decreased with an increase in the angle of incidence and increased with an increase in slope height. However, the vertical amplification ratios at the left toe increased with an increase in the angle of incidence, whereas their variations remained almost unchanged with varying slope heights. The horizontal amplification ratios at the right toe increased with an increase in  Figures 19 and 20, the variations in the amplification ratios at the slope toe were gentle. The horizontal amplification ratios at the left toe decreased with an increase in the angle of incidence and increased with an increase in slope height. However, the vertical amplification ratios at the left toe increased with an increase in the angle of incidence, whereas their variations remained almost unchanged with varying slope heights. The horizontal amplification ratios at the right toe increased with an increase in the angle of wave incidence (at a slope height of 0.2), whereas they first increased and then decreased as the incident angle increased (at slope heights of 0.5, 1.0 and 2.0). The vertical amplification ratios of the right toe grew with an increase in the angle of incidence. Nevertheless, these ratios increased significantly at angles of 10 • and 20 • (at a slope height of 1.0). The variations above were in good agreement with the seismic responses of the ground motions in Figure 18 (at a slope height of 1.0), which indicated that the variations in the horizontal and vertical amplification ratios were non-monotonous at the slope toe.

As shown in
The variations in the horizontal and vertical amplification ratios at the slope crest are complicated. The horizontal amplification ratios at the left crest monotonously rose or decreased with an increase in the angle of incidence (at slope heights of 0.2, 0.5 and 1.0). However, the horizontal ratios first increased and then decreased with an increase in the angles of incidence (at a slope height of 2.0). The maximum horizontal amplification ratios fell between 5 • and 15 • (at a slope height of 2.0). The horizontal amplification ratios at the right crest first increased and then decreased as the incident angle increased, and the maximum horizontal amplification ratios moved from the an incident angle of 30 • to those at 15 • as the slope height increased. The vertical amplification ratios of the left crest increased with an increase in the angle of incidence, and the maximum vertical amplification ratios of the left crest were obtained at a slope height of 1.0. The vertical amplification ratios of the right crest increased with an increase in the angle of incidence (at slope heights of 0.2 and 0.5). The variations in the vertical amplification ratios of the right crest were intensely complicated (at slope heights of 1.0 and 2.0), and the maximum ratios were obtained when the incident angles fell between 10 • and 20 • .
In summary, the variations in the amplification ratios at the slope toe were monotonous, that is, the ground motions changed regularly with varying incident angles and heights. Nonetheless, the variations in ground motions at the slope crest were complicated owing to the interactions between the input waves and the reflected waves in the topography.

Influence of the Width of Slope Crest W
The width of the slope crest varied from 50 m, 100 m and 200 m to 400 m, whereas the other geometrical parameters (slope heights and slope inclinations) remained the same. The acceleration amplification ratios along the slope surface are illustrated in Figure 21 with varying incident angles.  As illustrated in Figure 21, the acceleration amplification ratios at a slope width of 1.0 H were the greatest. This phenomenon can be indicated by the resonance of the model sizes and wavelengths. The ground motions were concentrated on the slope crest, especially when the slope width was 4.0 H. In addition, the maximum acceleration amplification ratios were obtained when the incident angles were between 5° and 10° at a width of 0.5 H, whereas the ratios reached a maximum when the incident angles fluctuated between 25° and 30° at a slope width of 4.0 H. Thus, the wider the slope crests, the closer the maximum acceleration amplification ratios were to the greater inclinations of incident As illustrated in Figure 21, the acceleration amplification ratios at a slope width of 1.0 H were the greatest. This phenomenon can be indicated by the resonance of the model sizes and wavelengths. The ground motions were concentrated on the slope crest, especially when the slope width was 4.0 H. In addition, the maximum acceleration amplification ratios were obtained when the incident angles were between 5 • and 10 • at a width of 0.5 H, whereas the ratios reached a maximum when the incident angles fluctuated between 25 • and 30 • at a slope width of 4.0 H. Thus, the wider the slope crests, the closer the maximum acceleration amplification ratios were to the greater inclinations of incident waves. In summary, the energies of earthquake waves were mainly concentrated at the crest of the slope topography. The maximum ground motions occurred as the width of the slope model was nearly one wavelength, especially in certain directions of seismic incident waves.
Moreover, the observation points (slope toe and slope crest) are analyzed in Figures 22 and 23, respectively. All the observation points are illustrated based on the horizontal and vertical amplification ratios, to analyze the comprehensive relations between the incident angles and slope widths. Figure 21. Variations of the seismic responses along the slope ridges and slope crests of the real amplification ratios (r) of PGA with varied angles of wave incidence at different slope widths, (a) slope width of 50 m, (b) slope width of 100 m, (c) slope width of 200 m (d) slope width of 400 m As illustrated in Figure 21, the acceleration amplification ratios at a slope width of 1.0 H were the greatest. This phenomenon can be indicated by the resonance of the model sizes and wavelengths. The ground motions were concentrated on the slope crest, especially when the slope width was 4.0 H. In addition, the maximum acceleration amplification ratios were obtained when the incident angles were between 5° and 10° at a width of 0.5 H, whereas the ratios reached a maximum when the incident angles fluctuated between 25° and 30° at a slope width of 4.0 H. Thus, the wider the slope crests, the closer the maximum acceleration amplification ratios were to the greater inclinations of incident waves. In summary, the energies of earthquake waves were mainly concentrated at the crest of the slope topography. The maximum ground motions occurred as the width of the slope model was nearly one wavelength, especially in certain directions of seismic incident waves.
Moreover, the observation points (slope toe and slope crest) are analyzed in Figures  22 and 23, respectively. All the observation points are illustrated based on the horizontal and vertical amplification ratios, to analyze the comprehensive relations between the incident angles and slope widths.  The horizontal amplification ratios at the slope toe decreased monotonously with an increase in the incident angle, whereas they increased as the slope width increased. The vertical amplification ratios at the left toe increased with an increase in the incident angle. The vertical amplification ratios at the right toe reached a maximum when the slope widths were between 1.0 and 2.0. With an increase in the incident angle, the vertical amplification ratios were almost the same for each slope width.
Nonetheless, variations in the acceleration amplification ratios at the slope crest were complicated. The horizontal amplification ratios at the left crest first decreased and The horizontal amplification ratios at the slope toe decreased monotonously with an increase in the incident angle, whereas they increased as the slope width increased. The vertical amplification ratios at the left toe increased with an increase in the incident angle. The vertical amplification ratios at the right toe reached a maximum when the slope widths were between 1.0 and 2.0. With an increase in the incident angle, the vertical amplification ratios were almost the same for each slope width.
Nonetheless, variations in the acceleration amplification ratios at the slope crest were complicated. The horizontal amplification ratios at the left crest first decreased and then increased with an increase in the incident angle (at slope widths of 2.0 and 4.0). The horizontal ratios first increased and then decreased as the incident angle increased (at slope widths of 0.5 and 1.0). In addition, the ratios increased from an incident angle of 5 • to an incident angle of 20 • at a slope width of 0.5, which could be interpreted by the superposition of the scattered waves in the narrow models. The vertical amplification ratios at the left crest first decreased and then increased with an increase in the angle of incidence. Furthermore, the horizontal amplification ratios at the right crest first increased and then decreased with an increase in the angle of incidence. The maximum horizontal ratios were obtained between 15 • and 25 • of the incident angle (at a slope width of 1.0). The vertical amplification ratios at the right crest were almost the same for various incident angles. The vertical ratios reached a maximum between 15 • and 25 • of the incident angles (at a slope width of 1.0). The above variations were in good agreement with the seismic responses of ground motions of PGA in Figure 21b, and these evolutions could be interpreted by the correlations between the wavelength and the slope widths.
In summary, the variations in the amplification ratios were mainly concentrated on the slope crest with varying angles of incidence, and these complexities were aggravated when the width was almost one wavelength.

Influence of Slope Inclination I
The inclination of the slope, which is defined as the width-depth ratio, varied from 0.5, 0.7, 1.0 and 1.5 to 2.0, whereas the other geometrical parameters (slope heights and slope widths) remained the same. The acceleration amplification ratios along the slope surface are illustrated in Figure 24 with varying angles of incidence. In Figure 24, the acceleration amplification ratios moved from the left slope crest to the right slope crest when the slope model became gentle. The amplification ratios moved away from the wave source with an increase in the angles of incidence, especially in the gentle slope model. These evolutions are indicated by the scattered waves induced by slope inclinations. The ground motions reached a maximum when the incident angles fluctuated between 5° and 10° in the steep slope model (i < 1.0). However, the ground motions reached a maximum when the incident angles were between 15° and 20° in the gentle slope model (i > 1.0). Thus, the greater the slope inclinations, the closer the maximum acceleration amplification ratios were to the greater inclinations of the incident waves. In summary, the ground motions were mainly concentrated on the slope ridges. The oblique incident waves and slope inclinations changed the propagation of the seismic wave paths, and the scattered waves easily concentrated as the slope became gentle. The maximum acceleration amplification ratios were obtained in the oblique direction of the incident waves.
Similar to the analyses of the slope height and slope widths, the slope toe and slope crest are also presented with horizontal amplification ratios and vertical amplification ratios in Figures 25 and 26, respectively. The comprehensive relations between the incident angles and slope inclinations are discussed. In Figure 24, the acceleration amplification ratios moved from the left slope crest to the right slope crest when the slope model became gentle. The amplification ratios moved away from the wave source with an increase in the angles of incidence, especially in the gentle slope model. These evolutions are indicated by the scattered waves induced by slope inclinations. The ground motions reached a maximum when the incident angles fluctuated between 5 • and 10 • in the steep slope model (i < 1.0). However, the ground motions reached a maximum when the incident angles were between 15 • and 20 • in the gentle slope model (i > 1.0). Thus, the greater the slope inclinations, the closer the maximum acceleration amplification ratios were to the greater inclinations of the incident waves. In summary, the ground motions were mainly concentrated on the slope ridges. The oblique incident waves and slope inclinations changed the propagation of the seismic wave paths, and the scattered waves easily concentrated as the slope became gentle. The maximum acceleration amplification ratios were obtained in the oblique direction of the incident waves.
Similar to the analyses of the slope height and slope widths, the slope toe and slope crest are also presented with horizontal amplification ratios and vertical amplification ratios in Figures 25 and 26, respectively. The comprehensive relations between the incident angles and slope inclinations are discussed.
motions reached a maximum when the incident angles were between 15° and 20° in the gentle slope model (i > 1.0). Thus, the greater the slope inclinations, the closer the maximum acceleration amplification ratios were to the greater inclinations of the incident waves. In summary, the ground motions were mainly concentrated on the slope ridges. The oblique incident waves and slope inclinations changed the propagation of the seismic wave paths, and the scattered waves easily concentrated as the slope became gentle. The maximum acceleration amplification ratios were obtained in the oblique direction of the incident waves.
Similar to the analyses of the slope height and slope widths, the slope toe and slope crest are also presented with horizontal amplification ratios and vertical amplification ratios in Figures 25 and 26, respectively. The comprehensive relations between the incident angles and slope inclinations are discussed.  The horizontal amplification ratios of the slope toe monotonously decreased with an increase in the angle of incidence, whereas the horizontal ratios increased as the slope inclination increased. That is, the ground motions at the slope toe increased as the slope became gentle. However, the vertical amplification ratios of the slope toe increased monotonously with an increase in the angle of incidence, and the ratios gradually decreased when the slope became gentle.
The variations in the amplification ratios at the slope crest were complicated at varying incident angles. The horizontal amplification ratios at the left crest decreased with an increase in the angle of incidence, and the reduced amplitude intensified as the slope became steeper (i < 1.0). The horizontal amplification ratios at the right crest remained the same with an increase in the angle of incidence, and the maximum ratios fell between 1.0 and 1.5 of i. The vertical amplification ratios at the left crest and right crest decreased with an increase in slope inclination. The ratios were intensified at the left crest whereas they were alleviated at a slope inclination of 1.0, at an incident angle of 15°on the right crest. The horizontal amplification ratios of the slope toe monotonously decreased with an increase in the angle of incidence, whereas the horizontal ratios increased as the slope inclination increased. That is, the ground motions at the slope toe increased as the slope became gentle. However, the vertical amplification ratios of the slope toe increased monotonously with an increase in the angle of incidence, and the ratios gradually decreased when the slope became gentle.
The variations in the amplification ratios at the slope crest were complicated at varying incident angles. The horizontal amplification ratios at the left crest decreased with an increase in the angle of incidence, and the reduced amplitude intensified as the slope became steeper (i < 1.0). The horizontal amplification ratios at the right crest remained the same with an increase in the angle of incidence, and the maximum ratios fell between 1.0 and 1.5 of i. The vertical amplification ratios at the left crest and right crest decreased with an increase in slope inclination. The ratios were intensified at the left crest whereas they were alleviated at a slope inclination of 1.0, at an incident angle of 15 • on the right crest.
In summary, the variations in the amplification ratios were concentrated on the slope crests at varying angles of incidence because the varying inclinations of slope ridges changed the propagation of wave paths in the topography.

Summary and Conclusions
In this study, ground motion amplification of the slope topography was analyzed using dynamic FEMs. The wave patterns, materials and slope sizes were discussed based on varying angles of incidence. A viscous-spring artificial boundary was borrowed, and an equivalent nodal force method was proposed to implement oblique incident waves in the FEM. Two numerical examples were adopted to verify the validity of the input method and the accuracy of the artificial boundary. Subsequently, the ground motions in slope topography were investigated with arbitrary directions of incidence, considering the main impact factors of site effects (wave patterns, materials and sizes of slope topography).
The main conclusions and findings are as follows: (1) the amplification effects are underestimated by the vertical incident waves, and the maximum ground motions are closely related to the scattered waves, which depend on the ridge inclinations and the incident directions. (2) Owing to the intense amplification effects under SV waves, the ground motions of slope topography are preferred for analysis using SV waves rather than P waves. (3) With an increase in the angles of incidence, the amplification effects are more complicated at the crest of the slope topography, and the amplification regions are focused on the ridges away from the epicenter. (4) The ground motion amplification effects are aggravated in soft materials owing to the much more scattered waves produced. The directions of the incident waves are a key factor in the analysis of ground motions in slope topography. The amplification effects of slope topography should be discussed separately owing to the complex transmission paths induced by oblique incident waves.
The results obtained here can be considered in the analysis of the ground motion amplification effects of slope topography subjected to earthquakes. The above findings were discussed in terms of three influencing factors that impact the topographic effects: the patterns of the incident waves, the materials of the slope and the slope sizes. The wave patterns that included P waves and SV waves were sufficient for the analysis of topographic effects. Even though the surface waves (e.g., Rayleigh and Love waves) have not been directly analyzed, they were still considered in this study because the conclusions were based on the incident body waves, that is, P waves or SV waves, and were actually a comprehensive result of body waves and surface waves. That is, the proposed input method can also express the impacts of Rayleigh waves (in two dimensions). Thus, the topographic effects calculated using P waves and SV waves in this study are meaningful and necessary. Even though the samples in the analysis of slope materials and slope sizes were relatively few, the conclusions could interpret several regularity phenomena, achieving the purpose of this study. Furthermore, the lack of experimental evidence for the comparison and calibration of the findings in the numerical studies is a serious obstacle [36,61]. Hence, a specific model should be established according to the study area while considering real materials based on the specific earthquake record for a specific study area. Data from field measurements should be collected simultaneously to support numerical modeling results.