Propagation of 2-D SH waves in random media: Insights from ab initio numerical simulations and transport theory

Increasing the frequency range of physics-based predictions of earthquake ground motions requires to account for small-scale heterogeneities, which can only be described in a stochastic way. Although many studies have addressed the impact of random heterogeneities on ground motion intensity parameters obtained by numerical simulation, very few have verified the accuracy of their numerical solutions or controlled the scattering regime they were simulating. Here we present a comprehensive analysis of SH wave propagation in 2D random media which covers a broad range of propagation regimes from ballistic to diffusive. The coherent and incoherent parts of the wavefield are examined independently. The random media consist in correlated density and velocity fluctuations described by a von Kármán autocorrelation function with a Hurst coefficient of 0.25 and a correlation length a = 500 m. The Birch correlation coefficient which relates density to velocity fluctuations takes 4 possible values between 0.5 and 1, and the standard deviation of the perturbations is either 5% or 10%. Spectral element simulations of SH wave propagation excited by a plane wave are performed for normalized wavenumbers (ka) up to 5. Analysis of the decay of the coherent wave amplitude, obtained through different averaging procedures, allows for a direct measure of the scattering attenuation, which we successfully compare with the predictions of the Dyson mean field theory. We also present the comparison between the energy envelopes measured from the synthetics and their theoretical counterpart provided by the radiative transfer theory and the diffusion approximation. Excellent agreement is found between numerical simulations and theoretical predictions of radiative transfer theory for the mean intensity. The numerical study highlights the difference of attenuation length between the mean field and the mean intensity. In the forward scattering regime, the peak intensity appears to decay exponentially over a length scale known as the transport mean free path. Furthermore, the fluctuations of intensity in the ballistic peak exhibit a transition from Log-normal to Exponential statistics. This transition occurs for a propagation distance of the order of the mean free path, which offers an alternative method of estimating this parameter.