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.     

Long-period corrections to body waves: theory

Summary. The usual asymptotic methods used to correct the high-frequency solutions of the wave equation are unsatisfactory as they do not give the low-frequency, partial reflections expected from a region of high velocity gradient. A new iterative solution is obtained which uses the first term of the Langer asymptotic expansion as the zeroth iterate. This satisfactorily gives the partial reflections from a region of high velocity gradient, even when they are generated near the turning point of the ray. Although the results are somewhat complicated in the frequency domain, in the time domain all types of wave interaction are described by six universal time functions. For any problem, these functions are scaled in time according to the depth of the interaction, and in strength according to the magnitude of the coupling parameter. Numerical results and approximations are given for these functions. Coupling parameters are investigated for acoustic and elastic waves in a plane model, and acoustic and elastic-gravitational waves in a spherical model. The same universal time functions allow the excitation of elastic waves to be studied when the source is in a region of high velocity gradient or is near the wave's turning point. Results are given for a moment tensor, point source in plane and spherical models.     

Diffraction of P, S and Rayleigh waves by three-dimensional topographies

The diffraction of P, S and Rayleigh waves by 3-D topographies in an elastic half-space is studied using a simplified indirect boundary element method (IBEM). This technique is based on the integral representation of the diffracted elastic fields in terms of single-layer boundary sources. It can be seen as a numerical realization of Huygens principle because diffracted waves are constructed at the boundaries from where they are radiated by means of boundary sources. A Fredholm integral equation of the second kind for such sources is obtained from the stress-free boundary conditions. A simplified discretization scheme for the numerical and analytical integration of the exact Green's functions, which employs circles of various sizes to cover most of the boundary surface, is used.
The incidence of elastic waves on 3-D topographical profiles is studied. We analyse the displacement amplitudes in the frequency, space and time domains. The results show that the vertical walls of a cylindrical cavity are strong diffractors producing emission of energy in all directions. In the case of a mountain and incident P, SV and SH waves the results show a great variability of the surface ground motion. These spatial variations are due to the interference between locally generated diffracted waves. A polarization analysis of the surface displacement at different locations shows that the diffracted waves are mostly surface and creeping waves.     

A near-regional explosion source model for tuff

Summary. A variety of near-regional (300 km) data, including spectral amplitudes of Pg , surface-wave forms, and close-in (5–10 km) accelerograms have been used to build an elastic seismic source model for a 1-Mton explosion in tuff at near-regional distances. The model consists of: (1) a pressure pulse which injects 3 × 10¹² cm³ of volume into the medium, (2) a vertical, upward force impulse that imparts 10¹⁸ dyn-s of momentum to the medium, each source component having a time duration of 0.6 s and a depth of 1.3 km. The force impulse appears to be required by two considerations: (a) the striking similarity, apart from sign, of explosion surface waves with those of their cavity collapses, (b) the observation of considerable SV energy leaving the source of the 1-Mton explosions JORUM and HANDLEY . Scaling curves have been constructed which fit the proposed source model. These scaling curves employ: very slow decrease, as (yield)^−0.10 of the primary corner frequency; decay as (frequency)⁴ or (frequency)³ to high frequency. While these scaling curves are unconventional, they appear to be the only ones which can satisfy the near-regional data. The slow scaling with yield of the spectral carner frequency suggests that it is caused by something other than the equivalent elastic radius, e.g. the time duration of motion at the source. The results, at odds with similar studies at teleseismic distances, suggest that significantly different equivalent elastic sources are required at near-regional (as compared with teleseismic) distances; therefore, the effect of the upward impulse might not be seen at teleseismic distances. Consequently, these results probably do not pertain to the seismic discrimination problem at teleseismic distances.     

On the decoupling of P and S waves in inhomogeneous elastic media

Summary. A representation derived by Richards (1974) for P-SV wave displacements in spherically symmetric elastic media is extended to general displacements in a general inhomogeneous isotropic elastic medium, with suitably differentiable elastic constants. This representation in terms of appropriate potentials gives rise to a partial decoupling of the P and S (weighted) displacements, which satisfy simple second order wave equations with lower order coupling terms. The highest order P and S components satisfy homogeneous wave equations that depend only on the P and S velocities α and β respectively and are unaffected by the density other than at the source and observer positions.     

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.     

Elastic wave propagation in model sediments —I

Summary. A fluid-saturated packing of like elastic spheres is used as a model of an oceanic sediment and a method is presented for calculating the effective velocities of elastic waves in such a medium. In particular the method is applied to low-frequency waves travelling vertically down a cubic packing, saturated with an inviscid fluid and initially at rest under a uniform compressive force. It is found that two waves propagate and moreover, that their velocities are not related through the usual equations of classical elasticity to the effective elastic moduli for static deformation of the packing. For a dry packing, there is found to exist a 'cut-off' frequency above which the wave decays with depth. An extension of the method to slightly viscous fluids is also given.     

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.     

Formulation of the multiple anisotropic scattering process in two dimensions for anisotropic source radiation

The radiative transfer theory (RTT) describes the energy transport through a random heterogeneous medium, neglecting phase information. It provides an adequate framework for modelling high-frequency seismogram envelopes. For isotropic scattering and sources, the radiative transfer equation (RTE) has been formulated analytically and numerically simulated using Monte Carlo methods for acoustic and elastic media. Here, we derive an exact analytical solution of the RTE in 2-D space for the acoustic case, including anisotropic scattering for a anisotropic point-like impulsive source. For this purpose, we generalize the path integral method, which has been used before in the isotropic case, to take into account the anisotropy of both the source radiation pattern and scattering processes, simultaneously. Then we obtain a general solution, which is written in a closed form in the Fourier space. To illustrate the theoretical results, we compute the full space and time evolution of the specific intensity for an arbitrary case. We also compare the time traces computed from our general solution with cases in which the source and/or the scattering process are isotropic. The importance of taking into account both anisotropies simultaneously becomes obvious in our examples. We also show that at long lapse time, our example approaches the solution of the diffusion equation.     

The PWaves From the Amchitka Island Explosions

Summary. There is evidence that the equivalent seismic sources of the Amchitka Island explosions — Longshot, Milrow and Cannikin — depart significantly from the simple model of a point compressional-source in a layered elastic-medium. Consequently modelling the observed seismograms using standard source-models may not be the most efficient method of determining source properties. Here an alternative to modelling is used to obtain information on the seismic sources due to the explosions. Broad-band (BB) estimates of the P signals are obtained from the short-period (SP) seismograms, corrected for attenuation, and interpreted in terms of P, pP and radiation from secondary sources. the main conclusions are:
(i) BB estimates of the radiated displacement from the explosions can be obtained with only a small reduction in the signal-to-noise ratio seen on SP seismograms;
(ii) observations of differences in pulse amplitudes and spectra are not necessarily due to differences in anelastic attenuation;
(iii) P and pP at a given station may differ in shape so that notches in the signal spectrum may not be related to source depth;
(iv) there is evidence of arrivals that others have identified as due to slap-down but which could be interpreted as an overshoot to pP;
(v) direct interpretation of the estimated ground displacement is a better procedure for determining the seismic source properties of explosions than modelling SP seismograms using idealised models as a starting point.     

Windowed f-k spectra of a three-dimensional wavefield excited by a point source in a two-dimensional multilayered elastic medium

When full 3-D modelling is too costly or cumbersome, computations of 3-D elastic wave propagation in laterally heterogeneous, multilayered 2-D geological structures may enhance considerably our ability to predict strong ground motion for seismological and engineering purposes. Towards this goal, we extend the method based on the combination of the thin-layer finite-element and boundary-element methods (TLFE-BEM) and calculate windowend f - k spectra of the 3-D wavefield. The windowed f - k spectra are spatially localized spectra from which the local properties of the wavefield can be extracted. The TLFE-BEM is particularly suited for calculating the complete wavefield where surface waves are dominant in multilayered media. The computations are performed in the frequency domain, providing the f - k spectra directly. From the results for the 3-D wavefield excited by a point source in a 2-D multilayered, sloped structure, it can be said that the phase velocity of the fundamental-mode Rayleigh wave in a laterally heterogeneous multilayered medium, estimated from the windowed f - k spectra, varies with the location of the point source. For the model calculated in this article, the phase velocity varies between the value for the flat layered structure of the thick-layer side and that for the structure just under the centre of the window. The exact subsurface structure just under the centre of an array in a laterally heterogeneous medium cannot be obtained if we use the f - k spectral analysis assuming a flat layered structure.     

Cagniard de-Hoop method for a Haskell type vertical fault

Summary. There are two ways to apply the Cagniard de-Hoop method when generating synthetic seismograms due to a source in a three-dimensional medium. One is the Hankel transform method (hereafter called 'the cylindrical wave representation method'), which utilizes the property of the modified Bessel function. The other is 'the plane wave representation method', which replaces the Bessel function by a superposition of plane waves. In extending a point source to a finite dimensional source, the latter method is extremely useful, because it enables the integrations on the fault surface to be performed analytically.
Using the latter method, expressions for displacements due to a Haskell type vertical fault in a uniform half-space are obtained. A solution is given as a sum of four quadrantal source contributions, which is similar to Madariaga's solution for a source in the whole space. Each contribution consists of a single finite range integration or a single integration plus a pole contribution, depending on the location of the observation point with respect to the source. The procedure can be extended to other fault types.     

Point Sources of Elastic Deformation: Elementary Sources, Static Displacements

Summary. We construct a catalogue of all possible elementary point sources for static deformation in an elastic solid. the familiar double-couples, CLVD's, centres of compression and dilatation, etc., are all members of the complete catalogue. the sources are classified according to the rank of the seismic moment tensor, and according to the weight (or order) of the irreducible tensor representation of the 3-D rotation group. These sources can be classified as belonging to one of three general classes. the static excitation functions are calculated for an infinite, homogeneous, isotropic medium for all these sources. We show that, except for sources belonging to these three general classes, all other sources — which are numerous for the tensors of high rank — are null static sources. That is, sources that do not produce any static displacement outside of the source region. Due to the presence of null sources, an inversion of the static deformation data is non-unique. the expansion of the equivalent-force tensors and the stress glut tensors (or seismic moment tensors) into a set of the symmetric trace-free source tensors is proposed. the sources corresponding to seismic moment tensors of the second, third and fourth ranks are considered in more detail. We identify the third-rank sources with rotational dislocations or disclinations.     

Time dependent piezomagnetic fields in viscoelastic medium

We investigated time dependent piezomagnetic fields due to volcanic sources embedded in a viscoelastic, homogeneous half-space. Especially in volcanic areas, the presence of inhomogeneous materials and high temperatures produce a lower effective viscosity of the Earth's crust that calls for considering anelastic properties of the medium. Piezomagnetic properties are carried by grains of titano-magnetite, which occupy only a small fraction of ordinary rock volume and are supposed to be elastic, while the non-magnetic surrounding matrix is assumed to be viscoelastic. From all the possible rheological models, we investigated two cases in which the bulk modulus is purely elastic and the shear modulus relaxes as: (i) a Maxwell solid and (ii) a standard linear solid (SLS). We applied the Correspondence Principle to the analytical elastic solutions for pressurized spherical sources and dislocation sources in order to determine the time dependent piezomagnetic fields in a viscoelastic medium. The piezomagnetic field completely vanishes after the relaxation process for a Maxwell rheology, whereas it is found to decrease over time and reach some finite offset value for a SLS rheology. These different behaviours provide helpful hints in understanding the temporal evolution of piezomagnetic anomalies in volcanic regions.     

Synthetic logs of multipole sources in boreholes based on the Kelvin–Voigt stress–strain relation

We design a numerical algorithm for wave simulation in a borehole due to multipole sources. The stress–strain relation of the formation is based on the Kelvin–Voigt mechanical model to describe the attenuation. The modelling, which requires two anelastic parameters and twice the spatial derivatives of the lossless case, simulates 3-D waves in an axisymmetric medium by using the Fourier and Chebyshev methods to compute the spatial derivatives along the vertical and horizontal directions, respectively. Instabilities of the Chebyshev differential operator due to the implementation of the fluid–solid boundary conditions are solved with a characteristic approach, where the characteristic variables are evaluated at the source central frequency. The algorithm uses two meshes to model the fluid and the solid. The presence of the logging tool is modelled by imposing rigid boundary conditions at the inner surface of the fluid mesh. Examples illustrating the propagation of waves are presented, namely, by using monopoles, dipoles and a quadrupoles as sources in hard and soft formations. Moreover, the presence of casing and layers is considered. The modelling correctly simulates the features—traveltime and attenuation—of the wave modes observed in sonic logs, namely, the P and S body waves, the Stoneley wave, and the dispersive S waves in the case of multipole sources.     

Kirchhoff–Helmholtz reflection seismograms in a laterally inhomogeneous multi-layered elastic medium – II. Computations

Summary. High-frequency reflection and refraction seismograms for laterally variable multi-layered elastic media are computed by using the frequency domain elastic Kirchhoff–Helmholtz (KH) theory of Frazer and Sen. Both source and receiver wavefields are expanded in series of generalized rays and then elastic (KH) theory is applied to determine the coupling between each source ray and each receiver ray at each interface. The motion at the receiver is given as a series of integrals, one for each generalized ray. We use geometrical optics and plane wave reflection and transmission coefficients for rapid evaluation of the integrand. When the source or the receiver ray field has caustics on the surface of integration geometrical ray theory breaks down and this gives rise to singularities in the KH integrand. We repair this using methods suggested by Frazer and Sen.
Examples of reflection seismograms for 2-D structures computed by elastic KH theory are shown. Those for a vertical fault scarp structure are compared with the seismograms obtained by physical modelling. Then OBS data obtained from the mid-America trench offshore Guatemala area are analysed by computing KH synthetics for a velocity model that has been proposed for that area. Our analysis indicates the existence of a small low-velocity zone off the trench axis.
No head wave arrivals are obtained in our KH synthetics since we do not consider multiple interactions of a ray with an interface. The nearly discontinuous behaviour of elastic R/T coefficients near the critical angle causes small spurious phases which arrive later than the correct arrivals.     

Comparison of global synthetic seismograms calculated using the spherical 2.5-D finite-difference method with observed long-period waveforms including data from the intra-Antarctic region

We have been developing an accurate and efficient numerical scheme, which uses the finite-difference method (FDM) in spherical coordinates, for the computation of global seismic wave propagation through laterally heterogeneous realistic Earth models. In the field of global seismology, traditional axisymmetric modeling has been used widely as an efficient approach since it can solve the 3-D elastodynamic equation in spherical coordinates on a 2-D cross-section of the Earth, assuming structures to be invariant with respect to the axis through the seismic source. However, it has the severe disadvantages that asymmetric structures about the axis cannot be incorporated and the source mechanisms with arbitrary shear dislocation have not been attempted for a long time. Our scheme is based on the framework of axisymmetric modeling but has been extended to treat asymmetric structures, arbitrary moment-tensor point sources, anelastic attenuation, and the Earth center which is a singularity of wave equations in spherical coordinates. All these types of schemes which solve 3-D wavefields on a 2-D model cross-section are classified as 2.5-D modeling, so we have named our scheme the spherical 2.5-D FDM. In this study, we compare synthetic seismograms calculated using our FDM scheme with three-component observed long-period seismograms including data from stations newly installed in Antarctica in conjunction with the International Polar Year (IPY) 2007–2008. Seismic data from inland Antarctica are expected to reveal images of the Earth's deep interior with enhanced resolution because of the high signal-to-noise ratio and wide extent of this region, in addition to the rarity of sampling paths along the rotation axis of the Earth. We calculate synthetic seismograms through the preliminary reference earth model (PREM) including attenuation using a moment-tensor point source for the November 9, 2009 Fiji earthquake. Our results show quite good agreement between synthetic and observed seismograms, which indicates the accuracy of observations in the Antarctica, as well as the feasibility of the spherical 2.5-D modeling scheme.     

Modelling of ground deformation and gravity fields using finite element method: an application to Etna volcano

Elastic finite element models are applied to investigate the effects of topography and medium heterogeneities on the surface deformation and the gravity field produced by volcanic pressure sources. Changes in the gravity field cannot be interpreted only in terms of gain of mass disregarding the ground deformation of the rocks surrounding the source. Contributions to gravity changes depend also on surface and subsurface mass redistribution driven by dilation of the volcanic source. Both ground deformation and gravity changes were firstly evaluated by solving a coupled axisymmetric problem to estimate the effects of topography and medium heterogeneities. Numerical results show significant discrepancies in the ground deformation and gravity field compared to those predicted by analytical solutions, which disregard topography, elastic heterogeneities and density subsurface structures. With this in mind, we reviewed the expected gravity changes accompanying the 1993–1997 inflation phase on Mt Etna by setting up a fully 3-D finite element model in which we used the real topography, to include the geometry, and seismic tomography, to infer the crustal heterogeneities. The inflation phase was clearly detected by different geodetic techniques (EDM, GPS, SAR and levelling data) that showed a uniform expansion of the overall volcano edifice. When the gravity data are integrated with ground deformation data and a coupled FEM modelling was solved, a mass intrusion could have occurred at depth to justify both ground deformation and gravity observations.