Weighted Residual Method for Diffraction of Plane P-Waves in a 2 D Elastic Half-Space Revisited : On an Almost Circular Arbitrary-Shaped Canyon

Scattering and diffraction of elastic in-plane Pand SV-waves by a surface topography such as an elastic canyon at the surface of a half-space is a classical problem which has been studied by earthquake engineers and strong-motion seismologists for over forty years. The case of out-of-plane SH-waves on the same elastic canyon that is semicircular in shape on the half-space surface is the first such problem that was solved by analytic closed-form solutions over forty years ago by Trifunac. The corresponding case of in-plane Pand SV-waves on the same circular canyon is a much more complicated problem because the in-plane Pand SVscattered-waves have different wave speeds and together they must have zero normal and shear stresses at the half-space surface. It is not until recently in 2014 that analytic solution for such problem is found by the author in the work of Lee and Liu. This paper uses the technique of Lee and Liu of defining these stress-free scattered waves to solve the problem of the scattering and diffraction of these in-plane waves on an on an almost-circular surface canyon that is arbitrary in shape.


Introduction
This paper studies the subject on the diffraction of in-plane P-waves in an elastic half-space by arbitrary-shaped canyons using the weighted residual method.It presents a solution for any arbitrary-shaped canyons where the depth of the canyon is approximately half the width of the canyon, such as a semicircle, ellipse, or trapezoid.
Researchers continue to study the effects of scattering and diffraction of waves on two-dimensional canyons in an elastic, isotropic, and homogeneous medium.These studies, which assist researchers to understand earthquake ground motions in and around topographic features, initially, addressed incident SH-waves [1][2][3].In solving SH-waves, the method of images, which assumes equal and opposite scattered waves upon reflection, has been used [4].Because of the mode conversion, an incident P-wave, which produces both a reflected P-and SV-wave, is more complex.Thus, the P-and SV-waves diffraction problems cannot be solved using the method of images.Much of the latest and most recent work on diffraction and soil-structure interaction for topographies involving analytic methods is so far thus limited mostly to SH-wave problems [5].
Recently in 2014, Lee and Liu analyzed the harmonic motion induced by an incident P-wave for a two-dimensional diffraction around a semicircular canyon in an elastic halfspace using an analytic solution to satisfy the zero-stress boundary conditions [6].In past approaches, numerical approximations of geometry and/or wave functions were made to satisfy the half-space boundary condition using wave functions that are a function of both sine and cosine [3,7].Lee and Liu's new method redefines the cylindrical-wave functions for both the longitudinal P-and shear SV-wave so that they are now a function of sine or cosine, but not a combination of both.Using the Fourier half-range expansion-and because each function itself is orthogonal in the half-space-the functions satisfy the zero-stress boundary condition along the half-space.The normal-and shear-stress boundary condition along the half-space in Figure 1 will be zero since the stress function contains sin(), which equals zero when  = 0, .This paper expands on Lee and Liu's new approach for the semicircular canyon replacing it with an arbitrary-shaped canyon with harmonic motion induced by an incident Pwave.In this study, the weighted residual method is applied for the solution of the wave function for the arbitrary-shaped canyon [8,9].Lee and Wu's method defined the scattered wave potentials as a combination of both the sine and the cosine functions, and the origin was defined above the halfspace for the case of the shallow canyon problem.This study uses the same principles from Lee and Liu's paper for which the P-and SV-cylindrical-wave functions are defined by the sine function only.In addition, Lee and Wu [9] define two sets of scattered P-and SV-waves, (a total of four sets of waves) so as to satisfy four sets of boundary conditions, two on the half-space and two on the canyon.The coordinate system in this study is located at the half-space surface, which allows for arbitrary-shaped canyons to be solved.The new method uses only one set of scattered P-and SV-waves, using the method of Lee and Liu [6] to automatically satisfy the freestress boundary conditions at the half-space surface.Using this improved weighted residual method, the results for Lee and Liu's semicircle were verified and new results for an ellipse, trapezoid, and rectangle are presented.

Model
The model for the canyon has no restrictions on shape other than the surface of the canyon which must be continuous and defined by a sequential number of points whose polar coordinates have an increasing value of .In Figure 1, the two-dimensional, arbitrary-shaped canyon in an elastic halfspace ( > 0) is defined.Each point will have an (, ) coordinate and once transformed into polar coordinates, each point will have an (, ) location.The geometry of the canyon is transformed from the rectangular coordinate system into a cylindrical coordinate system with the same origin at , shown as follows: ) . ( The figure encompasses the following: an incident plane Pwave defined by the potential   , with incidence angle   .
The incident angle is measured with respect to the horizontal -axis.The waves have a circular frequency  = 2, a longitudinal-wave velocity , and a transverse-wave velocity .
The half-space is elastic, isotropic, and homogeneous, with the following material properties: Lame constants  and  and mass density .The longitudinal-wave velocity  and transverse-wave velocity  are Therefore, the constants   = / (P-wave number) and   = / (SV-wave number) can be determined.

Harmonic Motion Induced by Incident P-Wave
Plane longitudinal-(P-) waves enter the half-space at angle   , resulting in displacements and a propagation in the - plane.The incident plane P-wave is characterized by a potential   and is defined in the - plane as The time factor is exp(−), where  = √ −1 and  is the time coordinate.The time factor is removed from the equation because we are only interested in studying the amplitude of the waves as a function of the canyon geometry and not the time aspect , which describes the time-dependent oscillation of the wave due to the mode conversion along the half-space surface at  = 0; the incident P-wave produces a reflected plane P-wave with P-wave potential   , reflection angle   , reflected plane SV-wave with a SV-wave potential   , and reflection angle   .The mode conversion takes place in order to satisfy the stress-free boundary conditions   = 0 and   = 0.The reflected P-and SV-wave potentials are derived in Achenbach's paper [10].The reflection angles,   and   , are calculated using Snell's law, for  >  implies   >   : On the flat half-space surface, the stress-free boundary conditions can be satisfied for the incident P and the reflected Pand SV-wave by utilizing the boundary conditions: This results in the following set of equations for the reflected P-wave and SV-wave potentials defined in the - plane as Reflected Plane P-Waves Reflected Plane SV-Waves The reflection coefficients  1 and  2 derived by Cao and Lee and are defined as [3] The free-field waves are unaffected by the canyon and become a combination of the input P-wave and reflected P-and SVwaves.The free-field P-wave potential  ff and the free-field SV-wave potential  ff are as follows: Both of these potentials must satisfy the zero-stress boundary conditions along the half-space surface; namely, for all  ≥  (same as (5)), The solution to the wave problems is to start by deriving the Helmholtz wave equation Prior research has shown that the solution Θ() to the Helmholtz wave equation for the scattered P-wave potential  = (, ) − in 2D cylindrical (polar) coordinates in the half-space region  ≥ 0, 0 ≤  ≤ , and in the full space − ≤  ≤ , is a combination of both the sin  and cos  terms included as two independent solutions with  and  as constants and shown as Lee and Liu [6] observed that in the half-space region, where  ≥ 0, 0 ≤  ≤ , the {sin } and {cos } terms are no longer two independent solutions.The authors concluded that the solution will either be a sine series or a cosine series, but not a combination of both.Therefore the solution becomes To solve the semicircular canyon, it is imperative that Sommerfeld's radiation condition (waves that are outgoing towards infinity) is satisfied.The Hankel function of the first kind,  (1)   (  ), is used since it represents an outgoing wave that decreases in amplitude, unlike the Hankel function of the second kind,  (2)   (  ), which represents an incoming wave.Therefore, using the theory asserted by Lee and Liu [6] that, in the half-space, the sine functions are complete and orthogonal, the scattered wave potentials will take only the sine terms: (1)   (  ) sin , with  0 =  0 = 0.Here  (1)   (⋅) are Hankel functions of the first kind corresponding to outgoing waves to infinity satisfying Sommerfeld's radiation condition.They can be defined in terms of Bessel functions of the 1st and 2nd kind which in turn are computed by algorithms defined and found in Abramowitz and Stegun [11].

Boundary Conditions for the Canyon Surface
Boundary conditions on the canyon surface must be satisfied for the free-field and scattered wave potentials.To satisfy these boundary conditions for no force on the canyon surface, the traction components-radial   , and angular   -on the surface are computed (see Figure 2(a)).Both traction components must satisfy the condition that there is zero stress as derived in the papers by Lee and Wu [8,9].The equations for the traction components (14) are derived from the stress equation (15): In Figure 2(b), the angle  is the angle between the radial vector r and the normal vector n at each point of the surface that is measured positive in the counterclockwise direction from the radial vector.The equations for the stresses are defined as [10,12]   = ∇ The stresses  +  ,  +  ,  +  are the in-plane stresses induced by the incident P-wave and reflected P-wave and SV-wave potentials.These stresses are calculated directly from (15) and by substituting the free-field potentials in (8).
On the canyon, the boundary condition of zero stress that must be satisfied on the surface of the arbitrary shape is [9] as follows: (18)

Application of Weighted Residual Method
The functions that define the stresses are an infinite summation.Therefore, the traction equations are an infinite summation.An approximate solution would use a finite summation with  terms and thus 2 unknowns,   , and   ,  = 1 to .The procedure for solving these two equations, traction equation (18), is a special case of the method of moments defined in Roger Harrington's paper, "Matrix methods for field problems" [13].A set of weighting functions,  1 ,  2 ,  3 , . . .,   , in the range of 0 to , is defined and applied to the traction terms.This results in 2 equations that require integration from 0 to .Applying the method of weighted residuals, the weighting function is chosen as For  = 1, 2, . . ., , the function is Since sin() equals zero when  = 0, (20) represent 2 pairs of equations with 2 unknowns,   and   , for  = 1, 2, . . ., .

Numerical Solutions
The equations for traction, (20), form a set of complex simultaneous equations with unknowns   and   , and coefficient matrix with terms   ,   and constant terms   , are shown as . . .
. . . ) . ( The terms of the coefficient matrix are determined by substituting the stress equations for the scattered waves (16) into the traction equations,   and   , (17) and multiplying each equation by the weighting function (19).Each term requires integration along the boundary of the canyon from 0 to .
The equations are integrated using the Gaussian quadrature method.The surface of integration-the canyon surfaceis divided into 400 segments defined by 401 points.Each segment is subdivided again by the 10-point Gauss-Legendre integration.This integration technique is highly accurate when the integrand is very smooth, which is the case for the sine and cosine functions, as the method converges much more quickly than in other integration schemes [14].
The constant terms  1 to  2 are calculated from tractions computed from the free-field waves.These free-field stresses are substituted into the equations for traction (17) and then multiplied by the weighting function (19) and integrated along the boundary from 0 to .
The set of 2 complex equations are solved for the unknowns  1 to   , and  1 to   .

Comparison of Results to Previous and Existing Studies
From the above analysis, which determined the coefficient   and   of the P-and SV-wave potentials, the displacement amplitudes can now be determined.These displacement amplitudes of interest occur on the surface of the half-space and within the canyon surface.The amplitudes are important in studying the variability of ground motions in the vicinity of a canyon.
The free-field displacements are calculated by substituting the incident-and reflected-wave potentials into The scatter-wave displacements are calculated by substituting the P-and SV-wave potentials into The resulting displacements for the scatter waves are in the following form [12]: with the displacement functions [12] given by  (3)  71 =  (3)  71 (,   ) =    (1)  −1 (  ) −  (1)   (  ) ,  (3)  72 =  (3)  72 (,   ) =  (1)   (  ) ,  (3)  81 =  (3)  81 (,   ) =  (1)   (  ) , In order to transform the scatter-wave displacements   ,   in polar coordinates to displacements   and   in rectangular coordinates, a transformation is used.The transformation is as follows: The total displacement amplitudes are a linear combination of the scattered and the free-field waves, and the result of the real (Re) and imaginary (Im) parts The phase angle of the points on the canyon and the surface of the half-space can be calculated by using phase These two-dimensional displacement amplitudes, horizontal (|  |) and vertical (|  |), and the phase angle are plotted versus the dimensionless horizontal distance / for a specific dimensionless frequency  and an angle of incidence   .The following plots use the dimensionless frequency parameters : where  is the radius of the canyon,   is the wave length of the shear wave at frequency ,  is the shear wave speed,   is the wave number, and  is the cyclic frequency defined as  = 2, with a Poisson ratio of ] = 0.25.To verify the validity of this numerical method, each plot was generated using an increasing value of  max -the total number of equationsuntil convergence was achieved.The figures have a depth-tohalf-width ratio of ℎ/ = 1 and their displacement amplitudes on the half-space surface and the surface of the canyon ( 1 = ) are plotted along the horizontal -axis in the interval −4 ≤ / ≤ 4. The point / = −1 corresponds to the left rim of the canyon, / = 0 to the bottom, and / = 1 to the right rim.The incident P-waves are assumed to arrive from the left (/ < 0) in all cases.The displacement amplitudes were computed and compared with the results obtained by the closed-form analytic solutions presented in Lee and Liu [6] for the model of a semicircular canyon in Figure 3.The amplitudes from Lee and Liu [6] are shown in Figure 4, and the calculated results here are shown in Figure 5, to be described in the next paragraph.
In Figure 5, the results are shown for the circular canyon with a coordinate system at the half-space (Figure 3) for   = 60,  = 10, and  max = 112 using the weighted residual method verified by the results of Lee and Liu [6].The results shown in Figure 5 can be compared with the results in Figure 4 for the incident P-wave on a 2D semicircular canyon, which was computed using Lee and Liu's analytic solution with the matrix equations of order from  = 104-112.The results matched Lee and Liu's results using  = 112 terms.

The Case of the Semicircular Canyon
The results for the case of the semicircular canyon can also be compared theoretically by matching the boundary condition equations derived from the weighted residual method with those derived from Lee and Liu's exact closed-form solution [6].For the incident P-wave, the total waves are a combination of the free-field waves and scattered waves as follows: (1)   (  ) sin , This assumes that the constants   ,   ,   , and   are defined as follows: and taking into consideration the stress-free boundary condition at the surface of the semicircular canyon at  = , for the case of a semicircular canyon,  = 0.The orthogonality of the sine function gives Therefore, the two boundary condition equations reduce down to ( (1)  11 ()   +  (1)  12 = −  2 (− (1)  41 ()   +  (1)  42 ( (1)  41 ()   +  (1)  42 ()   )   .
Equations ( 34) are two boundary conditions that are identical to those derived in Lee and Liu's [6] exact solution for a semicircular canyon.Therefore, with the choice of the sine function as the weighting function, the weighted residual methods resulted in the exact closed-form solution for the case of a semicircular canyon in an elastic half-space.

Application to Other Arbitrary Shapes
This new methodology for solving the arbitrary-shaped canyon is applicable to any arbitrary-shaped canyon.This section of the paper will look at an elliptical-shaped canyon and a trapezoidal-shaped canyon.For each shape, the number of equations,  max , increases until convergence is reached or the problem becomes numerically unstable.The procedure that has been applied in this paper, the method of moments using weighting functions, is an approximate technique and its accuracy is a function of convergence as the infinite series are truncated to a finite number of unknown terms:  1 , . . .,   and  1 , . . .,   .Three factors that affect the results are the number of terms, the stability of the Bessel functions, and integration.The cases of shallow canyons will be described and presented in the next paper.Figure 7 shows the results for an elliptical canyon of depth / = 1.25, a dimensionless frequency of  = 2, and an incident angle of  = 5 ∘ , 30 ∘ , 60 ∘ , and 90 ∘ .For the graphs with incident angle  = 5 ∘ , the -component displacement amplitudes slightly oscillate about the free-field amplitude on the front side of the canyon and produce a shadowy behavior along the backside of the canyon.Within the canyon surface, both the and -displacement amplitudes are oscillatory, producing a spike on the backside rim of the canyon at / = +1.The -component displacement amplitudes slightly oscillate about the free-field amplitude on the front side of the canyon.On the backside of the canyon, the oscillations are very steady, showing a shadowy behavior.Within the canyon surface, the displacement amplitudes are oscillatory with a spike at the rims of the canyon at / = +1.For the graphs with incident angle  = 30 ∘ , the component displacement amplitudes tend to oscillate about the free-field amplitude on the front side of the canyon but gradually oscillate and produce a shadowy behavior along the backside of the canyon.Within the canyon surface both the and -displacement amplitudes are highly oscillatory and produce a spike on the front side and backside rim of the canyons.The -component displacement amplitudes tend to oscillate about the free-field amplitude on the front side of the canyon and produce a shadowy behavior along the backside of the canyon.Within the canyon surface, the displacement amplitudes are highly oscillatory with spikes at the rims of the canyon / = +1 and / = −1.For the graphs with incident angle  = 60 ∘ , the and -component displacement amplitudes exhibit the same behavior as the graphs for incident angle 30 ∘ , with slightly smaller amplitudes and small corner spikes.For the graphs with incident angle of  = 90 ∘ , the and -component displacement amplitudes are symmetric about the origin 0. The free-field amplitudes exhibit decaying oscillatory amplification on the front side of the canyon and shadowy behavior on the backside of the canyon.Within the canyon surface, the displacement x/a amplitudes are slightly oscillatory with spikes at both rims of the canyon for displacement of / = +1 and / = −1 for the -component displacement amplitude.
Figure 8 shows the results for an elliptical canyon of depth / = 1.25, dimensionless frequency  = 8, and incidence angles of  = 5 ∘ , 30 ∘ , 60 ∘ , and 90 ∘ .For the graphs with incident angle  = 5 ∘ , the -component displacement amplitudes oscillate about the free-field amplitude on the front side of the canyon and produce a shadowy behavior along the backside.Within the canyon surface, both the and -displacement amplitudes are oscillatory, producing spikes on the front side of the canyon rim at / = −1 and the backside rim at / = +1.The -component displacement amplitudes oscillate about the free-field amplitude on the  x/a x/a Angle of incidence 14 Journal of Earthquakes front side and the backside of the canyon.The oscillations are very steady, showing a shadowy behavior.Within the canyon surface, the displacement amplitudes are oscillatory with spikes at the rims of the canyon at / = +1 and / = −1.For the graphs with incident angle  = 30 ∘ , the -component displacement amplitudes tend to rapidly oscillate about the free-field amplitude on the front side of the canyon but gradually oscillate on the backside of the canyon.Within the canyon surface, both the and -displacement amplitudes are highly oscillatory and only produce a spike on the backside rim of the canyon at / = +1.The component displacement amplitudes tend to oscillate about the free-field amplitude on the front side of the canyon.On the backside of the canyon, the oscillations are very steady, showing a shadowy behavior.Within the canyon surface, the displacement amplitudes are highly oscillatory with spikes at the canyon rims at / = +1 and / = −1.These spikes are due to the sharp corners, which result in a change of curvature.For the graphs with incident angle  = 60 ∘ , the and -component displacement amplitudes exhibit the same behavior as the graphs for the incident angle of 30 ∘ with slightly smaller amplitudes and no corner spike.For the graphs with incident angle  = 90 ∘ , the and -component displacement amplitudes are symmetric about the origin 0. The free-field amplitudes exhibit decaying oscillatory amplification on the front side of the canyon and a shadowy behavior on the backside of the canyon.Within the canyon surface, the displacement amplitudes are slightly oscillatory with spikes at both rims of the canyon for displacement of / = +1 and / = −1.
Figures 9(a) and 9(b) show the results for an elliptical canyon of depth / = 1.5, a dimensionless frequency of  = 2, and an incidence angle of  = 5 ∘ , 30 ∘ , 60 ∘ , and 90 ∘ .The displacement amplitudes tend to look similar along the free-field surface of the front side of the canyon and backside of the canyon as seen in Figure 7 for the ellipse with  = 1.25.Figure 9 shows significant spikes at the rims of the canyons for an angle of incidence of  = 60 ∘ .
In Figure 9, the phase diagrams are shown alongside the corresponding displacement amplitudes.All phase diagrams have been scaled by  and shifted arbitrarily to have a zero phase angle at / = 0.
Figure 10 shows the results for elliptical canyon depth / = 1.5, dimensionless frequency  = 8, and angles of incidence  = 5 ∘ , 30 ∘ , 60 ∘ , and 90 ∘ .The displacement amplitudes tend to look similar along the free-field surface of the front side of the canyon and backside of the canyon as those seen in Figure 8 for the ellipse with ℎ = 1.25.Figure 8 shows significant spikes at the rims of the canyons and much less oscillatory behavior within the surface of the canyon, producing shadow-like behavior.

Trapezoidal Canyons.
The results for the trapezoidal canyon, as shown in Figure 11, are presented in this section.The trapezoidal canyon is defined by ℎ/ = 1, with sloping sides of 60 ∘ and 45 ∘ measured from the horizontal axis.The incident P-waves are assumed to arrive from the front side of the canyon for all cases.Each figure shows the and -displacement amplitudes for angles of incidence  = 30 ∘ , 60 ∘ , and 90 ∘ .The case for  = 30 ∘ is for an oblique incident P-wave, and the case of  = 90 ∘ is for a vertical incident P-wave.The and -displacement amplitudes are plotted in the figures on the horizontal axis from / = −4 to / = +4.For convenience, each graph has three areas of focus: the front side, from / = −4 to / − 1, the backside, from / + 1 to / = +4, and the surface of the canyon, from / = −1 to / + 1.
Figure 12 shows the results for a trapezoidal canyon with depth of ℎ/ = 1, dimensionless frequency  = 2, a slope of 60 ∘ , and incidence angles of  = 30 ∘ , 60 ∘ , and 90 ∘ .The and -displacement amplitudes are plotted in the figures on the horizontal axis from / = −4 to / = +4.For the graphs with incident angle  = 5 ∘ , the -component displacement amplitudes gradually oscillate about the free-field amplitude on the front side of the canyon and produce a shadowy behavior along the backside of the canyon.Within the canyon surface, both the and -displacement amplitudes are oscillatory, producing clustered spikes at the canyon rims.The component displacement amplitudes oscillate about the freefield amplitude on the front side of the canyon.On the backside of the canyon, oscillations are very steady and show shadowy behavior.Within the canyon surface, the displacement amplitudes are oscillatory with clustered spikes at the canyon rims.For the graphs with incident angles  = 30 ∘ and  = 60 ∘ , similar trends are observed for the -component displacement amplitudes and -component displacement amplitude with larger amplitude-clustered spikes on the front-sloping side of the canyon (/ = −1 to / = −0.5)and on the backsloping side of the canyon (/ = +0.5 to / = +1).For the graphs with incident angle  = 90 ∘ , the -component displacement amplitudes are symmetric about the origin.They exhibit a shadowy behavior about the free-field amplitudes on the front side and backside of the canyon.Within the canyon surface, the displacement amplitudes are highly oscillatory with very large clustered spikes on the sloping sides.The -component displacement amplitudes are symmetric about the origin and exhibit a shadowy behavior along the freefield amplitudes on the front side and backside of the canyon.Within the canyon surface, displacement amplitudes are highly oscillatory with clustered spikes on the sloping sides.

Journal of Earthquakes
x/a Angle of incidence For the graphs with an incident angle of  = 5 ∘ , the component displacement amplitudes oscillate about the freefield amplitude on the front side of the canyon, and produce a shadowy behavior along the backside of the canyon.Within the canyon surface, both the and -displacement amplitudes are oscillatory, producing clustered spikes at the canyon rims.The -component displacement amplitudes oscillate about the free-field amplitude on the front side of the canyon.On the backside of the canyon, the oscillations are very steady, showing a shadowy behavior.Within the canyon surface, the displacement amplitudes are oscillatory with clustered spikes at the rims of the canyon.For the graphs with incident angle  = 30 ∘ , the -component displacement amplitudes tend to rapidly oscillate about the free-field amplitude on the front side of the canyon but gradually oscillate on the backside of the canyon.Within the canyon surface, the displacement amplitudes are highly oscillatory with clustered spikes on the front-sloping side (/ = −1 to / = −0.5)and on the back-sloping side of the canyon (/ = +0.5 to / = +1).The oscillations are much more oscillatory than the ones observed for deep elliptical canyons.The -component displacement amplitudes tend to oscillate about the free-field amplitude on the front side of the canyon.On the backside, the oscillations are very steady, showing a shadowy behavior.Within the canyon surface, the displacement amplitudes are highly oscillatory with small clustered spikes at the frontsloping side of the canyon and larger clustered spikes at the back-sloping side.The oscillations are at a higher frequency than the ones observed for the deep elliptical canyons.The graphs with incident angle  = 60 ∘ display the same trends as the ones for incidence angle  = 30 ∘ except without the spikes at the sloping sides of the trapezoidal canyon.For the graphs with incident angle  = 90 ∘ , the -component displacement amplitudes are symmetric about the origin and exhibit no amplitude change about the free-field amplitudes on the front side and backside of the canyon.Within the canyon surface, the displacement amplitudes are highly oscillatory with very large clustered spikes on the sloping sides.The -component displacement amplitudes are symmetric about the origin and exhibit a shadowy behavior along the free-field amplitudes on the front side and backside of the canyon.Within the canyon surface, displacement amplitudes are highly oscillatory with clustered spikes on the sloping sides.
Figure 14 shows the results for a trapezoidal canyon of depth ℎ/ = 1, dimensionless frequency  = 2, a slope of 45 ∘ , and angles of incidence  = 5 ∘ , 30 ∘ , 60 ∘ , and 90 ∘ .For the graphs with incident angle  = 5 ∘ , the -component displacement amplitudes gradually oscillate about the freefield amplitude on the front side of the canyon and produce a shadowy behavior along the backside of the canyon.Within the canyon surface both the and -displacement amplitudes are oscillatory.The -component displacement amplitudes exhibit a very steady shadowy behavior on the front side and backside of the canyon.Within the canyon surface, the displacement amplitudes are oscillatory.For the graphs with the incident angles of  = 30 ∘ and  = 60 ∘ , the and component displacement amplitudes exhibit the same trends as those seen for the incident angle  = 5 ∘ .For the graphs with incident angle  = 90 ∘ , the -component displacement amplitudes are symmetric about the origin and tend to exhibit a shadowy behavior along the free-field amplitudes on the front side and backside of the canyon.Within the canyon surface, the displacement amplitudes are highly oscillatory with very large clustered spikes at the rims of the canyon at / = +1 and / = −1.The -component displacement amplitudes are symmetric about the origin and they also exhibit a shadowy behavior about the free-field amplitudes on the front side and backside of the canyon.Within the canyon surface, the displacement amplitudes are highly oscillatory with clustered spikes at the center of the canyon.
Figure 15 shows the results for a trapezoidal canyon of depth ℎ/ = 1, dimensionless frequency  = 6, slope of 45 ∘ , and angles of incidence  = 5 ∘ , 30 ∘ , 60 ∘ , and 90 ∘ .For the graphs with an incident angle of  = 5 ∘ , the component displacement amplitudes oscillate about the freefield amplitude on the front side of the canyon and produce a shadowy behavior along the backside of the canyon.Within the canyon surface, both the and -displacement amplitudes are oscillatory, producing clustered spikes on the rims of the canyon rim.The -component displacement amplitudes oscillate about the free-field amplitude on the front side of the canyon.On the backside of the canyon, the oscillations are very steady, showing a shadowy behavior.Within the canyon surface, the displacement amplitudes are oscillatory with clustered spikes that peak at the center of the canyon.For the graphs with incident angle  = 30 ∘ , the component displacement amplitudes slightly oscillate about the free-field amplitudes on the front side of the canyon and exhibit a shadowy behavior along the backside of the canyon.Within the canyon surface, the displacement amplitudes are highly oscillatory, oscillating much more rapidly than the ones observed for the trapezoidal canyons with a slope of 60 ∘ .The -component displacement amplitudes exhibit a shadowy behavior about the free-field amplitudes on the front side and along the backside of the canyon.Within the canyon surface, the displacement amplitudes are highly oscillatory.For the graphs with incident angle  = 60 ∘ , the component displacement amplitudes rapidly oscillate about the free-field amplitudes on the front side of the canyon, but gradually oscillate on the backside of the canyon.Within the canyon surface, the displacement amplitudes are highly oscillatory with large clustered spikes at the rims of the canyon at / = −1 and / = +1.The -component displacement amplitudes gradually oscillate about the freefield amplitudes on the front side of the canyon and along the backside of the canyon.Within the canyon surface, the displacement amplitudes are highly oscillatory.For the graphs with the incident angle  = 90 ∘ , the -component displacement amplitudes are symmetric about the origin, and they tend to exhibit a shadowy behavior along the free-field amplitudes on the front side and the backside of the canyon.Within the canyon surface, the displacement amplitudes are highly oscillatory with very large clustered spikes at the canyon rims of / = +1 and / = −1.The -component displacement amplitudes are symmetric about the origin, and they also exhibit a shadowy behavior about the freefield amplitudes on the front side and around the backside of the canyon.Within the canyon surface, the displacement  amplitudes are highly oscillatory with clustered spikes at the center of the canyon.

Conclusions
Results of the current research have been presented for various shaped canyons: semicircular, elliptical, and trapezoidal.These results have been presented in Figures 4-15 and the following general conclusions are proposed: (1) This is the first time that Lee and Liu's redefined cylindrical-wave function has been used to satisfy the zero-stress boundary condition along the half-space for arbitrary-shaped canyons.The application of the weighted residual method for the solution of the wave function allows for a solution to noncircular arbitrary-shaped canyons.
(2) This analytical approach was applied to a semicircular canyon with  = 10,  max = 112, and  = 60 ∘ .The results were compared to those of Lee and Liu [6] and shown to match exactly.The same amplification of surface displacements, as seen in previous studies, is also demonstrated in these results.
(3) This analytical approach was applied to ellipticalshaped canyons with varying / (depth/width) ratios of 1.25 and 1.5.The resulting plots for the displacement amplitudes are for elliptical canyons with a higher frequency range than those for semicircular canyons in previous studies [14].
(4) The elliptical canyon with / = 1.25 produced total displacement amplitudes on both the and components for  = 8, which were highly oscillatory around the free-field surface on the front side of the canyon, and produced a more shadowy behavior along the backside of the canyon.Within the canyon, there were oscillations with spikes at the rims of the canyon at / = −1 and / = +1.The oscillations along the surface became more rapid as  increased.
Increasing the elliptical depth to  = 1.5 decreased the oscillations within the canyon surface and produced larger spikes with a greater magnitude at the frontside canyon rim at / = −1.
(5) This analytical approach was then applied to the trapezoidal canyon.The total displacement amplitudes of both the and -components on the half-space showed a trend of being oscillatory on the front side of the canyon around the free-field amplitudes.
On the backside of the canyon, the -component amplitudes were oscillatory about the free field, but the -component amplitudes produced a shadowy behavior around the free field.Both the and components are highly oscillatory about the canyon surface.
(6) At the trapezoidal canyon rims, spikes in the displacement amplitudes were observed.The spikes are consistent with the spikes that were observed for the semicircular and elliptical canyons.These spikes tend to increase in amplitude at higher dimensionless frequencies.
(7) Changing the slope of the trapezoidal canyon walls from 60 ∘ to 45 ∘ demonstrated that the amplitudes of both the component and -component appeared sensitive to the angle of incidence and the slope of the canyon sides.The amplification of the -component amplitudes and -component amplitudes changed from the canyon sides to the canyon bottom.The degree of the spikes of the canyon rims also changed.
(8) For the trapezoidal canyon, larger values of  (greater than 6) caused convergence problems.Further study is necessary.
(9) This method provides good results for arbitraryshaped canyons in which the radius to the canyon surface is approximately equal to .The solutions became numerically unstable as the radius to the canyon boundary became small, because the Bessel functions caused numerical problems.(10) The method used in this paper, which is an expansion of Lee and Liu's [6] redefined cylindrical-wave function with the application of the weighted residual method, can be extended to other types of incident waves such as SH and SV, as well as surface Rayleigh waves in an elastic homogeneous half-space.(11) The method used in this paper can be applied to canyons of any shape, but convergence of the solutions may be difficult.This is in particular true when the irregular canyon has nonsmooth sharp corners, and when parts of the canyon are too shallow being close to the origin.The sequel to this paper will involve arbitrary-shaped canyons that are much shallower in depth from the half-space surface.
(12) This method can next be applied to arbitrary-shaped cylindrical valleys, where the valley medium is softer or harder than that of the half-space.
(13) This method can also be extended for problems involving three-dimensional elastic-wave propagations using the same concept and methodology but (different) three-dimensional spherical Hankel functions instead of the two-dimensional cylindrical Hankel functions here.
(14) For underground arbitrary-shaped cavities, the above concept of weighted residues can be applied, but the stress-free boundary conditions at the half-space surface will have to be addressed and taken care of, in addition to that at the cavity's surface.

Summary
The two-dimensional diffraction of incident P-waves around an arbitrary-shaped canyon in an elastic half-space is presented in this paper.The scattered wave potentials for the resulting P-waves and S-waves are defined by an infinite series of terms with Hankel functions and only sine terms.Using the zero-stress boundary conditions, a solution is made using

Figure 1 :
Figure 1: Arbitrary-shaped with coordinates at the half-space.