Time-averaged and time-dependent energy-related quantities of harmonic waves in inhomogeneous viscoelastic anisotropic media

The energy–flux vector and other energy-related quantities play an important role in various wave propagation problems. In acoustics and seismology, the main attention has been devoted to the time-averaged energy flux of time-harmonic wavefields propagating in non-dissipative, isotropic and anisotropic media. In this paper, we investigate the energy–flux vector and other energy-related quantities of wavefields propagating in inhomogeneous anisotropic viscoelastic media. These quantities satisfy energy-balance equations, which have, as we show, formally different forms for real-valued wavefields with arbitrary time dependence and for time-harmonic wavefields. In case of time-harmonic wavefields, we study both time-averaged and time-dependent constituents of the energy-related quantities. We show that the energy-balance equations for time-harmonic wavefields can be obtained in two different ways. First, using real-valued wavefields satisfying the real-valued equation of motion and stress–strain relation. Second, using complex-valued wavefields satisfying the complex-valued equation of motion and stress–strain relation. The former approach yields simple results only for particularly simple viscoelastic models, such as the Kelvin–Voigt model. The latter approach is considerably more general and can be applied to viscoelastic models of unrestricted anisotropy and viscoelasticity. Both approaches, when applied to the Kelvin–Voigt viscoelastic model, yield the same expressions for the time-averaged and time-dependent constituents of all energy-related quantities and the same energy-balance equations. This indicates that the approach based on complex-valued representation of the wavefield may be used for time harmonic waves quite universally. This study also shows importance of joint consideration of time-averaged and time-dependent constituents of the energy-related quantities in some applications.     

A modified Lax-Wendroff correction for wave propagation in media described by Zener elements

A modified Lax-Wendroff correction for wave propagation in attenuating and dispersive media described by Zener elements is presented. As opposed to the full correction, this new technique is explicit and offers large computational savings. The technique may be applied to a wide variety of hyperbolic problems. Here, the concept is illustrated for wave propagation in visco-acoustic media.     

Stress transfer relations among the earthquakes that occurred in Kerman province, southern Iran since 1981

We explore the possible stress triggering relationship of the   M ≥ 6.4  earthquakes that occurred in Kerman Province, southern Iran since 1981. We calculated stress changes due to both coseismic sudden movement in the upper crust and the time-dependent viscous relaxation of the lower crust and/or upper mantle following the event. Four events of   M ≥ 6.4  occurred between 1981 and 2005, on and close to the Gowk fault, show a clear Coulomb stress load to failure relationship. The  2003 M = 6.5  Bam earthquake, however, which occurred approximately 95 km SW of the closest Gowk event, shows a very weak stress relation to preceding earthquakes. The coseismic static stress change at the hypocentre of the Bam earthquake is quite small (∼0.006 bars). The time-dependent post-seismic stress change could be 26 times larger or 7 times lower than that of coseismic static stress alone depending on the choice of viscoelastic crustal model and the effective coefficient of friction. Given the uncertainties in the viscoelastic earth models and the effective coefficient of friction, we cannot confidently conclude that the 2003 Bam event was brought closer to failure through coseismic or post-seismic stress loading. Interestingly, the southern Gowk segment with a similar strike to that of the Bam fault, experienced a stress load of up to 8.3 bars between 1981 and 2003, and is yet to have a damaging earthquake.     

Effect of the viscoelastic lithosphere on polar wander speed caused by the Late Pleistocene glacial cycles

Previous studies of the wander of the rotation pole associated with the Late Pleistocene glacial cycles indicate that the predicted polar wander speed is sensitive to the density jump at the 670 km discontinuity, the thickness of the elastic lithosphere, and the lower mantle viscosity. In particular, the M1 mode related to the density jump at 670 km depth has been shown to contribute a dominant portion of predicted polar wander speed for sufficiently small lower mantle viscosities. In this study, we examine the sensitivity of polar wander to variations in the viscosity of the viscoelastic lithosphere using simplified compressible Maxwell viscoelastic earth models. Model calculations for earth models with a viscoelastic lithosphere of finite viscosity indicate that the contribution of the M1 mode is similar to those associated with the density discontinuity at the core–mantle boundary (C0 mode) and the lithosphere (L0 mode). We speculate that this is due to the interaction between the M1 mode and the transient mode associated with the viscoelastic lithosphere, which reduces the magnitude of polar wander rates. Therefore, the M1 mode does not contribute a dominant portion of the predicted polar wander speed for earth models with a viscoelastic lithosphere of finite viscosity. In this case, predictions of polar wander speed as a function of lower mantle viscosity exhibit the qualitative form of an 'inverted parabola', as predicted for the J ˙₂ curve. We caution, however, that these results are obtained for simplified earth models, and the results for seismological earth models such as PREM may be complicated by the interaction between the M1 mode and the large set of transient modes.