Theoretical and numerical investigations of global and regional seismic wave propagation in weakly anisotropic earth models

Smith and Dahlen demonstrated that in a weakly anisotropic earth model the relative surface wave phase‐speed perturbation δc/c may be written in the form , where An and Bn are frequency‐dependent depth integrals and ζ denotes the ray azimuth. In this approximation, surface wave anisotropy is governed by 13 elastic parameters and the azimuthal dependence of the phase speed is represented by an even Fourier series in ζ involving degrees zero (five elastic parameters), two (six elastic parameters), and four (two elastic parameters). Jech and Ps̆enc̆ík demonstrated that in such a weakly anisotropic earth model the relative compressional‐wave phase‐speed perturbation may be expressed as δc/c= (2c2)−1B33, whereas the relative shear wave phase‐speed perturbations are given by δc/c= (4c2)−1{B11+B22±[(B11−B22)2+ 4B212]1/2}. We demonstrate that the coefficients B33, B11, B22, and B12 may be written in the generic form , where ζ denotes the local azimuth and i the local angle of incidence. For B11, B22 and B33 the coefficients an(i) and bn(i) are an even Fourier series of degrees zero, two and four in i, but for B12 they are an odd Fourier series of degrees one and three. Like surface waves, the azimuthal (ζ) dependence of body waves involves even degrees zero (five elastic parameters), two (six elastic parameters), and four (two elastic parameters), but, unlike surface waves, it also involves the odd degrees one (six elastic parameters) and three (two elastic parameters). Thus, weakly anisotropic body‐wave propagation involves all 21 independent elastic parameters. We use spectral‐element simulations of global and regional seismic wave propagation to assess the validity of these asymptotic body and surface wave results. The numerical simulations and asymptotic predictions are in good agreement for anisotropy at the 5 per cent level.