Numerical modeling of seismic waves for the radial-heterogeneous spherical Earth (Q2729328)

The article is devoted to the numerical simulation of seismic waves propagation in the heterogeneous spherical Earth. The authors consider a spherical layer with coordinates \(R\), \(\theta\), \(\varphi\), where \(r_0\leq R\leq R_0\), \(0 < \theta \leq \pi\), \(0\leq \varphi < 2\pi\). The displacement vector is defined by the formula NEWLINE\[NEWLINE \vec U = U_{\varphi}(r,\theta,t)\vec e_{\varphi}. NEWLINE\]NEWLINE The function \(U_{\varphi}\) satisfies the following equation which describes the SH wave propagation in the radial-heterogeneous spherical Earth: NEWLINE\[NEWLINE \begin{multlined} \frac{\mu}{r^2}\frac{\partial}{\partial\theta}\left(\frac{ \partial U_{\varphi}} {\partial \theta} - U_{\varphi}\cot\theta\right) + \frac{\partial}{\partial r}\left[\mu\left(\frac{\partial U_{\varphi}} {\partial r} - \frac{U_{\varphi}}{r}\right)\right] + \\ 3\frac{\mu}{r}\left(\frac{\partial U_{\varphi}} {\partial r} - \frac{U_{\varphi}}{r}\right) + \frac{2\mu}{r^2}\left(\frac{\partial U_{\varphi}} {\partial \theta} - U_{\varphi}\cot\theta\right)\cot\theta = \rho\frac{\partial^2 U_{\varphi}}{\partial t^2}, \end{multlined} NEWLINE\]NEWLINE the initial data NEWLINE\[NEWLINE U_{\varphi_{t=0}} = \left|\frac{\partial U_{\varphi}}{\partial t}\right|_{t=0} = 0, NEWLINE\]NEWLINE and the boundary conditions NEWLINE\[NEWLINE \mu\left|\left(\frac{\partial U_{\varphi}} {\partial r} - \frac{U_{\varphi}}{r}\right)\right|_{r = R_0} = F(\theta,t), \quad \mu\left|\left(\frac{\partial U_{\varphi}} {\partial r} - \frac{U_{\varphi}}{r}\right)\right|_{r = r_0} = 0. NEWLINE\]NEWLINE Here the function \(F(\theta,t)\) is given by the impulse SH-torque: NEWLINE\[NEWLINE F(\theta,t) = \frac{\delta(\theta)}{2\pi R_0^2\sin\theta}\,f(t). NEWLINE\]NEWLINE The Lame parameter \(\mu(r)\) and the density \(\rho(r)\) are piecewise continuous functions.NEWLINENEWLINE To solve the problem the authors expose an algorithm based on using the integral Laguerre transform with respect to the time variable instead of the Fourier transform. The algebraic equations obtained have matrix independent of the degree of the Laguerre polynomials. For solving the algebraic system the authors employ the sweep method. The calculation results are presented in the form of synthetic seismograms of the SH-wave for the so-called preliminary reference Earth model and for the above-indicated source pulse. The synthetic seismograms are computed up to the epicentral distance of \(70'\) and predominant period \(T_0 = 30\) s. The time is measured in terms of the predominant period.