SH Wave Number Green’s Function for a Layered, Elastic Half-Space. Part I: Theory and Dynamic Canyon Response by the Discrete Wave Number Boundary Element Method

We present a closed-form frequency-wave number (ω – k) Green’s function for a layered, elastic half-space under SH wave propagation. It is shown that for every (ω – k) pair, the fundamental solution exhibits two distinctive features: (1) the original layered system can be reduced to a system composed by the uppermost superficial layer over an equivalent half-space; (2) the fundamental solution can be partitioned into three different fundamental solutions, each one carrying out a different physical interpretation, i.e., an equivalent half-space, source image impact, and dispersive wave effect, respectively. Such an interpretation allows the proper use of analytical and numerical integration schemes, and ensures the correct assessment of Cauchy principal value integrals. Our method is based upon a stiffness-matrix scheme, and as a first approach we assume that observation points and the impulsive SH line-source are spatially located within the uppermost superficial layer. We use a discrete wave number boundary element strategy to test the benefits of our fundamental solution. We benchmark our results against reported solutions for an infinitely long circular canyon subjected to oblique incident SH waves within a homogeneous half-space. Our results show an almost exact agreement with previous studies. We further shed light on the impact of horizontal strata by examining the dynamic response of the circular canyon to oblique incident SH waves under different layered half-space configurations and incident angles. Our results show that modifications in the layering structure manifest by larger peak ground responses, and stronger spatial variability due to interactions of the canyon geometry with trapped Love waves in combination with impedance contrast effects.