Response of a source on top of a vertically stratified half-space

The solution of the response of a source on top of a horizontally stratified half-space is well-known. However, when the half-space is vertically stratified, the problem can only be solved with numerical methods like the finite element method. Here a semi-analytic approach to solve such a problem is presented. The three-dimensional variation of the problem is reduced to two-dimensional variation by Fourier transform in one coordinate variable. The remaining two-dimensional problem is solved by finding the eigensolution in each of the half-spaces. The eigensolutions of each region are found from the partial differential equation directly using the same basis set of expansion functions. This makes the calculation of the reflection and transmission operators very efficient. The reflection and transmission operators account for the mode-conversion, reflection, and transmission of the waves. With the reflection and transmission operators, the field everywhere can be calculated. The solution reduces to that of the Sommerfeld´s half-space problem when the two half-spaces are homogeneous.