Displacement, strain and stress fields due to shear and tensile dislocations in a viscoelastic half-space

Okada (1992) provided expressions for the displacement and strain fields due to a finite rectangular source in an elastic, homogeneous and isotropic half-space. Starting with these results, we applied the correspondence principle of linear viscoelasticity to derive the quasi-static displacement, strain and stress fields in a viscoelastic, homogeneous and isotropic half-space. We assume that the medium deforms viscoelastically with respect to both the shear and the normal stresses but keeps a constant bulk modulus; in particular, the shear modulus relaxes as Maxwell fluid. We presented the viscoelastic effect on displacement, displacement gradient and stress fields, for a choice of parameter values. The viscoelastic effect due to the sudden dislocation reaches a limit value after about 10 times the Maxwell time. The expressions obtained here provide tools for the study of viscoelastic relaxation of lithosphere associated with seismic and volcanic phenomena.     

Surface-wave propagation in a cracked poroelastic half-space lying under a uniform layer of fluid

The effect of cracks on the elastic properties of an isotropic elastic solid is studied when the cracks are saturated with a soft fluid. A polynomial equation in effective Poisson's ratio is obtained, whose coefficients are functions of Poisson's ratio of the uncracked solid, crack density and saturating fluid parameter. Elastic and dynamical constants used in Blot's theory of wave propagation in poroelastic solids are modified for the introduction of cracks. The effects of cracks on the velocities of three types of waves are observed numerically. The frequency equation is derived for the propagation of Rayleigh-type surface waves in a saturated poroelastic half-space lying under a uniform layer of liquid. Dispersion curves for a particular model of oceanic crust containing cracks are plotted. The effects of variations in crack density and saturation on the phase and group velocity are also analysed.     

Transient electromagnetic fields about an infinitesimally long grounded horizontal electric dipole on the surface of a uniform half-space

Summary. Laplace and Bessel Transforms are used to solve for the transient behaviour of the electromagnetic fields after switching off a steady current in a grounded infinitesimal horizontal dipole on the surface of a uniformly conducting half-space. Simple analytic expressions, which are valid for times sufficiently long after the switch that displacement terms can be ignored, are obtained on the surface of the half-space for the electric field and the time derivative of the magnetic field. At the instant of switching an infinitesimally long image becomes established directly under the source dipole. It is the diffusion of this image which gives the vertical magnetic field and horizontal electric fields their transient behaviour. During the transient, there is also a decaying charge distribution on the surface.     

Exact synthetic seismograms for an inhomogeneity in a layered elastic half-space

Summary. The propagation of a pulsed elastic wave in the following geometry is considered. An elastic half-space has a surface layer of a different material and the layer furthermore contains a bounded 3-D inhomogeneity. The exciting source is an explosion, modelled as an isotropic pressure point source with Gaussian behaviour in time.
The time-harmonic problem is solved using the null field approach (the T matrix method), and a frequency integral then gives the time-domain response. The main tools of the null field approach are integral representations containing the free space Green's dyadic, expansions in plane and spherical vector wave functions, and transformations between plane and spherical vector wave functions. It should be noted that the null field approach gives the solution to the full elastodynamic equations with, in principle, an arbitrarily high accuracy. Thus no ray approximations or the like are used. The main numerical limitation is that only low and intermediate frequencies, in the sense that the diameter of the inhomogeneity can only be a few wavelengths, can be considered.
The numerical examples show synthetic seismograms consisting of data from 15 observation points at increasing distances from the source. The normal component of the velocity field is computed and the anomalous field due to the inhomogeneity is sometimes shown separately. The shape of the inhomogeneity, the location and depth of the source, and the material parameters are all varied to illustrate the relative importance of the various parameters. Several specific wave types can be identified in the seismograms: Rayleigh waves, direct and reflected P -waves, and head waves.     

Ground deformation in a viscoelastic medium composed of a layer overlying a half-space: a comparison between point and extended sources

We obtain and compare analytical and numerical solutions for ground displacement caused by an overpressurized magma chamber placed in a linear viscoelastic medium composed of a layer over a half-space. Different parameters such as size, depth and shape of the chamber, crustal rheology and topography are considered and discussed. Numerical solutions for an axisymmetric extended source are computed using a finite element method (FEM). Analytical solutions for a point source are obtained using the dislocation theory and the propagator matrix technique. In both cases, the elastic solutions are used together with the correspondence principle of linear viscoelasticity to obtain the solution in the Laplace transform domain. Viscoelastic solutions in the time domain are derived inverting the Laplace transform using the Prony series method. The differences between the results allow us to constrain the applicability of the point source and the flat surface hypothesis, which are usually implicitly assumed when analytical solutions are derived. The effect of the topography is also considered. The results obtained show that neglecting the topographic effects may, in some cases, introduce an error greater than that implicit in the point-source hypothesis. Therefore, for an adequate modelling and interpretation of the time-dependent displacements, topography must be considered.     

Dynamic stress variations due to shear faults in a plane-layered medium

A complete set of expressions is presented for the computation of elastic dynamic stress in plane-layered media. We use a discrete-wavenumber reflectivity method to compute the stress field radiated by arbitrary moment-tensor sources. The expressions derived here represent an interesting tool for both-the observational and theoretical analysis of dynamic stress changes associated with earthquake phenomena. Dynamic stress changes associated with a strike-slip fault having unilateral rupture are shown. This modelling, which is similar to the 1992 Landers California earthquake, illustrates the effects of distance, directivity and depth on transient stress changes.     

Transient electromagnetic field above a permeable and conducting half-space

Summary. The transient fields resulting from an abrupt current switch-off in a vertically oriented finite loop and a magnetic dipole above a permeable and conducting half-space have been investigated by Fourier methods, utilizing an appropriate Green's function. Two partial fields are identified, one 'radiative' and evanescent and one 'diffusive'. Asymptotic formulae for the diffusive field above the interface, applicable for late times, are developed, and from these the effects of height of the transmitting source and the permeability contrast are calculated. It is shown that for 'late enough' times and at sufficient distances the dipole formulae provide adequate approximations to the field due to a finite loop.     

Quasi-static deformation of a poroelastic half-space with anisotropic permeability by two-dimensional surface loads

An analytical solution is obtained of the fully coupled diffusion–deformation system of equations governing the quasi-static plane strain deformation of a poroelastic half-space with anisotropic permeability and compressible constituents. The stresses and the pore pressure are taken as the basic state variables. Displacements are obtained by integrating the coupled constitutive relations. The problem of surface loads is discussed in detail. Explicit analytical solutions are derived for normal line loading, shear line loading and normal strip loading. The permeability anisotropy is found to have a significant effect on the quasi-static deformation of the half-space. However, in the drained and undrained limits, the anisotropy has no effect. The stresses in the drained and undrained states are independent of the poroelastic parameters. Numerical computation of the pore pressure indicates that ignoring permeability anisotropy may lead to an overestimation of the pore pressure at points vertically below the point of normal loading. Further, anisotropy in permeability may lead to a dilution in the theoretical prediction of the Mandel–Cryer Effect.     

Implementation of perfectly matched layers in an arbitrary geometrical boundary for elastic wave modelling

The perfectly matched layer (PML) absorbing boundary condition is incorporated into an irregular-grid elastic-wave modelling scheme, thus resulting in an irregular-grid PML method. We develop the irregular-grid PML method using the local coordinate system based PML splitting equations and integral formulation of the PML equations. The irregular-grid PML method is implemented under a discretization of triangular grid cells, which has the ability to absorb incident waves in arbitrary directions. This allows the PML absorbing layer to be imposed along arbitrary geometrical boundaries. As a result, the computational domain can be constructed with smaller nodes, for instance, to represent the 2-D half-space by a semi-circle rather than a rectangle. By using a smooth artificial boundary, the irregular-grid PML method can also avoid the special treatments to the corners, which lead to complex computer implementations in the conventional PML method. We implement the irregular-grid PML method in both 2-D elastic isotropic and anisotropic media. The numerical simulations of a VTI lamb's problem, wave propagation in an isotropic elastic medium with curved surface and in a TTI medium demonstrate the good behaviour of the irregular-grid PML method.     

Transient SH-waves in a wedge-shaped layer

Summary. A generalized ray theory for transient SH -waves in a wedge-shaped layer over an elastic half-space is developed in this paper. The ray integrals for multiply reflected waves in the layer are derived in terms of two systems of coordinates and two sets of local wavenumbers, one along the free surface and the other along the sloped interface. All local wavenumbers are then transformed to a common wavenumber in all ray integrals which are evaluated by the Cagniard method. Results for the first motion approximation are in agreement with previous investigations.     

Seismic body waves in anisotropic media: reflection and refraction at a plane interface

Summary. A formulation is derived for calculating the energy division among waves generated by plane waves incident on a boundary between generally anisotropic media. A comprehensive account is presented for P, SV and SH waves incident from an isotropic half-space on an orthorhombic olivine half-space, where the interface is parallel to a plane of elastic symmetry. For comparison, a less anisotropic medium having transverse isotropy with a horizontal axis of symmetry is also considered. The particle motion polarizations of waves in anisotropic medium differ greatly from the polarizations in isotropic media, and are an important diagnostic of the presence of anisotropy. Incident P and SV waves generate quasi- SH waves, and incident SH waves generate quasi- P and quasi- SV waves, often of considerable relative magnitude. The direction of energy transport diverges from the propagation direction.     

Coseismic well-level changes due to the 1992 Roermond earthquake compared to static deformation of half-space solutions

The M _w 5.4 Roermond earthquake of 1992 April 13 was one of the strongest events during the last 500 years in Central Europe. For the period March–May 1992, we collected records of 194 continuously operating well-level sensors, mostly located within 120  km of the epicentre. Nearly all wells penetrate unconfined or poorly confined Quaternary deposits with high hydraulic conductivities. 81 out of 194 raw data sets show a significant dynamic or step-like response of centimetre amplitude to the passage of seismic waves. Precursory anomalies are not obvious in these records. Coseismic well-level fluctuations could reflect a redistribution of stress and pore pressure in the brittle crust. Systematic analyses of such fluctuations may improve our knowledge of the role of pore fluids in crustal rheology and earthquake mechanics. The rather high number of individual observational records for a single event allows a regional correlation of the signs and amplitudes of the coseismic steps to changes in volume strain caused by the earthquake. The coseismic strain field at the surface was calculated for a homogeneous and a layered half-space. The results show reasonable agreement in the sign of the well-level steps but the amplitudes predicted from the wells' volumetric strain responses are much smaller than those that were recorded. Clearly, the coseismic well-level steps cannot be explained by volume strain changes, as derived from linear elastic models.     

The correspondence between gravitational attraction and surface displacement due to volume expansion

Summary. It is shown that the formulas for the calculation of surface displacement due to a volume source in a half-space may be derived from a potential of the same form as the gravitational potential. The similarity between these expressions allows the use of methods developed for gravity in the study of surface displacement. As an example of the applicability of the correspondence, an efficient algorithm is developed for the calculation of uplift due to volume expansion within a homogeneous half-space. The method may be generalized to three-dimensional volumes undergoing expansion in a half-space.     

Synthetic SH seismograms in a laterally varying medium by the discrete wavenumber method

Summary. We present a new method to calculate the SH wavefield produced by a seismic source in a half-space with an irregular buried interface. The diffracting interface is represented by a distribution of body forces. The Green's functions needed to solve the boundary conditions are evaluated using the discrete wavenumber method. Our approach relies on the introduction of a periodicity in the source-medium configuration and on the discretization of the interface at regular spacing. The technique developed is applicable to boundaries of arbitrary shapes and is valid at all frequencies. Some examples of calculation in simple configurations are presented showing the capabilities of the method.     

Spectral aspects of the Gaussian beam method: reflection from a homogeneous half-space

Summary. Use of paraxially approximated Gaussian beams continues to be actively pursued for construction of synthetic seismograms in complicated environments. How to select the beams in the stack remains a source of difficulty which has primarily been addressed by semi-heuristic considerations. In this paper, the classical example of line-source field reflection from a homogeneous half-space that can sustain a head wave is examined from a plane-wave spectral point of view. The individual beam fields are modelled exactly by the complex source point technique, which emphasizes the complex spectral content of these wave objects. The quality of the paraxial approximation of a typical reflected (Gaussian) beam characterized by different parameters is examined from this perspective, and is compared with uniform and non-uniform asymptotics generated from the exact beam field spectral integral. With this information as background, the reflected field for a real line-source is synthesized by beam superposition. Except for the immediate vicinity of the critical reflection angle, the well-known failure of narrow paraxial beams, no matter how densely stacked, to reproduce the head wave effects is shown to be due to the inadequate spectral content of these beams and not to the failure of beam stacking per se. When the rigorous solutions are used for the narrow-waist beams, even relatively few suffice to yield agreement with the exact solution. This circumstance emphasizes the importance of fully understanding the spectral implications of various beam stacking schemes.     

Love waves excited by discontinuous propagation of a rupture front

Summary. Analytical results are presented for Love waves generated by sudden changes of the rate of advance of a curved rupture front in an inclined fault plane that is embedded in an elastic half-space. The boundary condition at the surface of the half-space approximates the presence of an overlying layer. The calculation consists of two parts. First, ray theory is used to calculate far-field approximations to the horizontally polarized wavefields which are emitted when the speed of the rupture front suddenly changes. These fields can be expressed as products of emission coefficients (which govern the angular dependence) and propagation terms. Secondly, a representation integral for the Love wave over a surface enclosing the rupture front is constructed, using the emitted signal and an appropriate Green's function. This integral is evaluated asymptotically. The resulting approximate Love-wave spectrum shows an explicit dependence on the nature of the rupture process, on the rupture-front and fault-plane geometry, and on the magnitude of a sudden change in the rate of advance of the rupture front.     

Scattering of horizontally-polarized shear waves by surface irregularities

Summary. The problem of the scattering of harmonic SH waves by an arbitrary surface irregularity in an otherwise semi-infinite elastic, homogeneous, isotropic two-dimensional half-space is examined in this study in order to ascertain the effect of topography on this type of seismic ground motion and to develop a useful scheme which can realistically handle arbitrary two-dimensional topography. Three geometric models are considered: a semicircular hill which is of academic interest; a mountain with a Gaussian shape which utilizes realistic dimensions and the combination of a ridge and a depression that models a region in Sylmar, California.
A singular Fredholm integral equation of the second kind for the displacement at the free surface is developed and solved numerically. In the case of the semicircular hill, horizontal ground motion can be more than twice that occurring in the case of smooth topography. The mountain simulated by a Gaussian profile experiences at its crest amplifications for certain angles of incidence and de-amplifications for other angles of incidence, as well as displacements whose amplitudes vary slowly with frequency on the side of the mountain which is in the same direction as the incident waves. The ridge-depression combination which is approximated by a sixth-order polynomial actually experienced shattered earth at its ridge crest during the San Fernando, California earthquake of 1971. This amplification is also exhibited by the results of the analysis which, predicts amplifications of over 75 per cent at the top of the ridge for waves arriving on the same side as the ridge.     

B-polarization induction in two generalized thin sheets at the surface of a conducting half-space

Summary. A new closed-form solution is obtained analytically for a B- polarization induction problem of geophysical interest, in which a local region of the Earth is represented by a generalized thin sheet at the surface of and in electrical contact with a uniformly conducting half-space. The generalized sheet, first introduced by Ranganayaki & Madden, is a mathematical idealization of a double layer which consists, in this problem, of two adjacent half-planes with distinct conductances representing a surface conductivity discontinuity such as an ocean—coast boundary, underlain by a uniform sheet of finite integrated resistivity representing the lower crust. The resistive sheet exerts a considerable mathematical influence on the solution causing, under certain conditions, an additional pole to appear in one of the forms of contour integral by which the solution can be expressed; it also weakens or eliminates field singularities that would otherwise occur at the conductance discontinuity. A numerical calculation is made for model parameters typifying an ocean—coast boundary underlain by a highly resistive crust. It is found that the residue of the pole associated with the resistive sheet dominates the solution for this example, the main consequence of which is a huge increase in the horizontal range over which the induced currents adjust themselves between the different 'skin-effect' distributions at infinity on either side of the model. Moreover the solution shows that this 'adjustment distance' has a more complicated dependence on the conductance and integrated resistivity of the sheet than that given simply by the square root of their product which was the length parameter proposed by Ranganayaki & Madden.     

A perturbation theory for Love waves in anelastic media

Summary. We give a systematic formulation and a rigorous justification of a perturbation technique for the computation of the eigenvalues and eigenfunctions of Love waves (and toroidal oscillations by an appropriate change for variables) in an anelastic medium with a constitutive law modelling geophysical media of current interest such as the Kelvin—Voigt Solid, the Maxwell Solid, the Standard Linear Solid, and the Standard Linear Solid with a continuous spectrum of relaxation times. We develop expressions relating the eigenvalues of eigenfunctions for Love waves in a continuously varying vertically stratified anelastic half-space to the corresponding elastic eigenvalues and eigenfunctions. Analytically, our correspondence principle has the form of a regular perturbation expansion in terms of a parameter for both the eigenvalues and eigenfunctions. The identification of ɛ is motivated by the dissipativity principle of viscoelasticity theory. Moreover, we show that our correspondence principle applies respectively only in the high and low frequency range for the Maxwell and Kelvin—Voigt Solids. Outside of the applicable range of frequencies, our correspondence principle yields no useful information. For the family of Standard Linear Solids it is uniformly applicable for all non-zero frequencies.
We also derive an explicit formula to estimate the radius of convergence of our perturbation expansions. This estimate of the radius of convergence for each eigenvalue and eigenfunction is functionally defined by the constitutive model for the anelastic medium. The estimate is frequency dependent and depends on the separation distance between the eigenvalue and the remainder of the spectrum of the corresponding elastic problem.