Effects of volcano topography on seismic broad‐band waveforms

Volcano seismology often deals with rather shallow seismic sources and seismic stations deployed in their near field. The complex stratigraphy on volcanoes and near‐field source effects have a strong impact on the seismic wavefield, complicating the interpretation techniques that are usually employed in earthquake seismology. In addition, as most volcanoes have a pronounced topography, the interference of the seismic wavefield with the stress‐free surface results in severe waveform perturbations that affect seismic interpretation methods. In this study we deal predominantly with the surface effects, but take into account the impact of a typical volcano stratigraphy as well as near‐field source effects. We derive a correction term for plane seismic waves and a plane‐free surface such that for smooth topographies the effect of the free surface can be totally removed. Seismo‐volcanic sources radiate energy in a broad frequency range with a correspondingly wide range of different Fresnel zones. A 2‐D boundary element method is employed to study how the size of the Fresnel zone is dependent on source depth, dominant wavelength and topography in order to estimate the limits of the plane wave approximation. This approximation remains valid if the dominant wavelength does not exceed twice the source depth. Further aspects of this study concern particle motion analysis to locate point sources and the influence of the stratigraphy on particle motions. Furthermore, the deployment strategy of seismic instruments on volcanoes, as well as the direct interpretation of the broad‐band waveforms in terms of pressure fluctuations in the volcanic plumbing system, are discussed.     

One-way representations of seismic data

A general one-way representation of seismic data can be obtained by substituting a Green's one-way wavefield matrix into a reciprocity theorem of the convolution type for one-way wavefields. From this general one-way representation, several special cases can be derived.
By introducing a Green's one-way wavefield matrix for primaries , a generalized Bremmer series representation is obtained. Terminating this series after the first-order term yields a primary representation of seismic reflection data. According to this representation, primary seismic reflection data are proportional to a reflection operator, 'modified' by primary propagators for downgoing and upgoing waves. For seismic imaging, these propagators need to be inverted. Stable inverse primary propagators can easily be obtained from a one-way reciprocity theorem of the correlation type.
By introducing a Green's one-way wavefield matrix for generalized primaries , an alternative representation is obtained in which multiple scattering is organized quite differently (in comparison with the generalized Bremmer series representation). According to the generalized primary representation, full seismic reflection data are proportional to a reflection operator, 'modified' by generalized primary propagators for downgoing and upgoing waves. Internal multiple scattering is fully included in the generalized primary propagators {either via a series expansion or in a parametrized way). Stable inverse generalized primary propagators can be obtained from the one-way reciprocity theorem of the correlation type. These inverse propagators are the nucleus for seismic imaging techniques that take the angle-dependent dispersion effects due to fine-layering into account.     

Inverse scattering of surface waves: imaging of near-surface heterogeneities

An efficient inverse scattering method is developed for imaging near-surface heterogeneities using scattered surface waves. Three dimensional elastodynamic wave propagation and scattering in a laterally invariant embedding medium is considered. The Born Approximation is used and the scattered wavefield is expressed as a domain type integral representation. The computation time of Green's tensor elements is reduced by considering the radial symmetry of the medium. The method is validated by numerical tests. Ultrasonic laboratory data obtained from a scale model experiment are used for imaging the near-surface inhomogeneities caused by an epoxy-filled hole in the surface of an aluminum block. Both synthetic and the scale model tests show that the location, the actual density contrast and the depth of the inhomogeneities are reasonably well estimated.     

Multipathing and spectral modulation of the downgoing wavefield in a complex crust: an example from the KTB (Germany) borehole

VSP data collected in the KTB (Germany) borehole to a depth of 8.5 km in 1999 show a surprising spectral modulation of the downgoing wavefield. After filtering the data with the singular value decomposition technique it was found that below about 6.2 km there are two depth intervals where the modulation can be explained in terms of a basic wavelet plus two weighted and delayed copies of that wavelet, with the delay for each wavelet remaining almost constant in each interval. The boundary between the two intervals is at about 7.25 km depth and above and below this depth the delay for the second wavelet is almost the same, while the delay for the third wavelet is significantly different. Neither the modulation nor its depth variation are source related and cannot be explained in terms of multiple reflections in a subhorizontal low-velocity layer. On the other hand, finite difference synthetic data show that subvertical layering (which is prevalent in the borehole area) provides a mechanism that can explain the observations. This mechanism has analogies with the generation of the standard refracted (i.e. head) waves. When a plane wave front propagates perpendicular to the boundaries of a vertical low-velocity layer surrounded by two vertical high-velocity layers, refracted wave fronts are generated in the low-velocity layer, which in turn generate secondary wave fronts in the high-velocity layers. These wave fronts trail the primary wave fronts by a constant delay whose magnitude has a simple dependence on the thickness of the low-velocity layer and the velocities involved. This process creates multipath arrivals that in geological settings with steeply inclined and faulted layers may appear and disappear rather abruptly, which may contribute to a scattered appearance of the wavefield.     

Inferring the relative measure of principal stress components

The phase velocity and the attenuation coefficient of compressional seismic waves, propagating in poroelastic, fluid-saturated, laminated sediments, are computed analytically from first principles. The wavefield is found to be strongly affected by the medium heterogeneity. Impedance fluctuations lead to poroelastic scattering; variations of the layer compressibilities cause inter-layer flow (a 1-D macroscopic local flow). These effects result in significant attenuation and dispersion of the seismic wavefield, even in the surface seismic frequency range, 10–100 Hz. The various attenuation mechanisms are found to be approximately additive, dominated by inter-layer flow at very low frequencies. Elastic scattering is important over a broad frequency range from seismic to sonic frequencies. Biot's global flow (the relative displacement of solid frame and fluid) contributes mainly in the range of ultrasonic frequencies. From the seismic frequency range up to ultrasonic frequencies, attenuation due to heterogeneity is strongly enhanced compared to homogeneous Biot models. Simple analytical expressions for the P -wave phase velocity and attenuation coefficient are presented as functions of frequency and of statistical medium parameters (correlation lengths, variances). These results automatically include different asymptotic approximations, such as poroelastic Backus averaging in the quasi-static and the no-flow limits, geometrical optics, and intermediate frequency ranges.     

Phase-array decomposition of a seismic section

A seismic re fraction/wide-angle reflection profile is analysed for the presence of correlated events ('phases'). The correlation problem is formulated in terms of temporally, spatially and frequency-local complex covariances. For robustness, the method concentrates on phase rather than amplitude information. This allows a computationally efficient algorithm that can make allowance for signal correlation length and can model curved wavefronts. A statistical test based on residual phase misfit across the analysed subarray is used to assess the probability that a detected event represents a real correlated signal.
With our chosen analysis parameters and confidence level (over 99.9 per cent). 1222 events were detected in the data. Using simple techniques based on 1-D earth models, detected events are associated with a small number of particular wave types. In this way, we have succeeded in classifying almost 95 per cent of the detected events. Those that remain describe those components of the data that are inconsistent with our simple ray paths in the 1-D assumption and with our prescribed tolerance. These include reverberations, near-surface guided waves and reflected waves from strongly laterally inhomogeneous structures. According to our modelling, about 25 per cent of the detected events are consistent with simple P -wave reflected energy, and these are to a very large extent (over 85 per cent) distinct from all the other wave-type models we have used. A direct mapping of the detected events into the offset-depth domain reveals dear internal and external consistencies among the detections for the various wave types. Estimated earth structure is consistent with models from previous analyses based on much larger data sets.
We have thus succeeded in extracting correlated events from the data and decomposing these, approximately but meaningfully, into distinct classes (ray paths)     

Numerical simulation of the propagation of P waves in fractured media

We study the propagation of P waves through media containing open fractures by performing numerical simulations. The important parameter in such problems is the ratio between crack length and incident wavelength. When the wavelength of the incident wavefield is close to or shorter than the crack length, the scattered waves are efficiently excited and the attenuation of the primary waves can be observed on synthetic seismograms. On the other hand, when the incident wavelength is greater than the crack length, we can simulate the anisotropic behaviour of fractured media resulting from the scattering of seismic waves by the cracks through the time delay of the arrival of the transmitted wave. The method of calculation used is a boundary element method in which the Green's functions are computed by the discrete wavenumber method. For simplicity, the 2-D elastodynamic diffraction problem is considered. The rock matrix is supposed to be elastic, isotropic and homogeneous, while the cracks are all empty and have the same length and strike direction. An iterative method of calculation of the diffracted wavefield is developed in the case where a large number of cracks are present in order to reduce the computation time. The attenuation factor Q ⁻¹ of the direct waves passing through a fractured zone is measured in several frequency bands. We observe that the attenuation factor Q ⁻¹ of the direct P wave peaks around kd = 2, where k is the incident wavenumber and d the crack length, and decreases proportionally to ( kd ) ⁻¹ in the high-wavenumber range. In the long-wavelength domain, the velocity of the direct P wave measured for two different crack realizations is very close to the value predicted by Hudson's theory on the overall elastic properties of fractured materials.     

Backscattering from spherical elastic inclusions and accuracy of the Kirchhoff approximation for curved interfaces

The Kirchhoff (or tangent plane) approximation, derived from the theoretically complete Kirchhoff–Helmholtz integral representation for the seismic wavefield, has been used extensively for the analysis of seismic-wave scattering from irregular interfaces; however, the accuracy of this method for curved interfaces has not been rigorously established. This paper describes an efficient Kirchhoff algorithm to simulate scattered waves from an arbitrarily curved interface in an elastic medium. Synthetic seismograms computed using this algorithm are compared with exact synthetics computed using analytical formulae for scattering of plane P waves by a spherical elastic inclusion. A windowing technique is used to remove strong internal reverberations from the analytical solution. Although the Kirchhoff method tends to underestimate the total scattering intensity, the accuracy of the approximation improves with increasing value of the wavenumber-radius product, kR . The arrival times and pulse shapes of primary reflections from the sphere are well approximated using the Kirchhoff approach regardless of curvature of the scattering surface, but the amplitudes are significantly underestimated for kR ≤ 5. The results of this work provide some new guidelines to assess the accuracy of Kirchhoff-synthetic seismograms for curved interfaces.     

An alternative method for calculation of Rayleigh and Love wave phase velocities by using three-component records on a single circular array without a central station

A method for the computation of phase velocities of surface waves from microtremor waveforms is shown. The technique starts from simultaneous three-component records obtained in a circular array without a central station. Then, Fourier spectra of vertical, radial and tangential components of motion are calculated for each station and considered as complex-valuated functions of the azimuthal coordinate. A couple of intermediate real physical quantities, B and C , can be defined from the 0- and ±1-order coefficients of the Fourier series expansion of such functions. Finally, phase velocities of Rayleigh and Love waves can be retrieved from B and C by solving respective one-unknown equations. The basic assumption is the possibility of expanding the wavefield as a sum of plane surface waves with Rayleigh and Love wavenumbers being univocal functions of the circular frequency. The method is tested in synthetic ambient noise wavefields confirming its reliability and robustness for passive seismic surveying.     

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.     

Inversion of seismic refraction data in planar dipping structure

Summary. A new method is presented for the direct inversion of seismic refraction data in dipping planar structure. Three recording geometries, each consisting of two common-shot profiles, are considered: reversed, split, and roll-along profiles. Inversion is achieved via slant stacking the common-shot wavefield to obtain a delay time–slowness (tau– p ) wavefield. The tau– p curves from two shotpoints describing the critical raypath of refracted and post-critically reflected arrivals are automatically picked using coherency measurements and the two curves are jointly used to calculate velocity and dip of isovelocity lines iteratively, thereby obtaining the final two-dimensional velocity model.
This procedure has been successfully applied to synthetic seismograms calculated for a dipping structure and to field data from central California. The results indicate that direct inversion of closely-spaced refraction/wide-aperture reflection data can practically be achieved in laterally inhomogeneous structures.     

Effect of a simple mountain range on underground seismic motion

Summary. The influence of a simple mountain range on seismic ground motion is studied. A two-dimensional model of the medium and vertically incident plane SH -waves are considered. Attention is devoted not only to the wavefield along the Earth's surface, but also within the medium. The wavefield is computed in two steps: (1) the computation of the impulse response by the finite-difference method, (2) the computation of the response to a time variation of the incident wave. approximately corresponding to a 'hundred-year' local earthquake at the site. Numerical results (the impulse response, the transfer function, the accelerograms, and their spectra) indicate strong spatial variability of the wavefield due to the topographical anomaly. The differences as large as 100 per cent in the peak amplitudes of the accelerograms, and of the order of 5 Hz in their predominant frequencies, at closely (∼ 200 m) spaced internal points of the medium have been found. Attention is focused also on the effect of causal absorption. Even unrealistically strong absorption seems to be unable to reduce significantly the spatial variability of the ground motion, caused by the topography. A variability like this. implying the occurrence of underground differential strains, might be of engineering importance in the antiseismic design of underground structures (tunnels, for example) in mountainous seismic regions. The ground-motion variability along the surface of the mountain is considerably smaller than within the medium.     

Regional tomographic inversion of the amplitude and phase of Rayleigh waves with 2-D sensitivity kernels

In this study, we test the adequacy of 2-D sensitivity kernels for fundamental-mode Rayleigh waves based on the single-scattering (Born) approximation to account for the effects of heterogeneous structure on the wavefield in a regional surface wave study. The calculated phase and amplitude data using the 2-D sensitivity kernels are compared to phase and amplitude data obtained from seismic waveforms synthesized by the pseudo-spectral method for plane Rayleigh waves propagating through heterogeneous structure. We find that the kernels can accurately predict the perturbation of the wavefield even when the size of anomaly is larger than one wavelength. The only exception is a systematic bias in the amplitude within the anomaly itself due to a site response.
An inversion method of surface wave tomography based on the sensitivity kernels is developed and applied to synthesized data obtained from a numerical simulation modelling Rayleigh wave propagation over checkerboard structure. By comparing recovered images to input structure, we illustrate that the method can almost completely recover anomalies within an array of stations when the size of the anomalies is larger than or close to one wavelength of the surface waves. Surface wave amplitude contains important information about Earth structure and should be inverted together with phase data in surface wave tomography.     

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.     

Crustal scattering at the KTB from a combined microearthquake and receiver analysis

In December of 1994 a fluid injection experiment which triggered several hundreds of microearthquakes was conducted at the KTB main borehole (Oberpfalz, Germany). These events were recorded with a temporal seismic network at the surface. Out of the complete data set, a cluster of five events recorded at four mini-arrays consisting of eight or nine stations was used to investigate the crustal scattering properties in the vicinity of the KTB. For this purpose, the 'Double Beam Method' (DBM; Krüger et al . 1993 , 1995 , 1996) and the 'Double Beam Imaging Method' (Scherbaum, Krüger & Weber 1997) were extended to curved wave fronts to drop the restriction of plane-wave propagation. This technique is used for imaging the crustal scattering strength using earthquake clusters recorded at close-by mini-arrays. The results of the array analysis show that the composition of the P coda is mainly affected by the site location of the arrays. Near-surface and deeper crustal scattering contribute in a very complicated pattern. Furthermore, with the present data set it was possible to identify reflections from the top of the Erbendorf Body. This is a very pronounced arrival in most of the recorded traces. In one of the arrays its amplitudes are even greater than the direct P phases. Five to eight coherent phases could be identified by the mini-arrays. Using only these phases, synthetic P -coda traces were constructed, which only contain the coherent part of the observed wavefield. By subtracting the synthetic coherent wavefield from the original traces we achieve a variance reduction in the P coda of up to 37 per cent. This leads to the conclusion that a large amount of the P coda at the KTB can be modelled by a simple deterministic single-scattering model using a small number of individual scatterers.     

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.     

Air-Coupled Seismic Waves at Long Range from Apollo Launchings

Summary Microphones and seismographs were co-located in arrays on Skidaway Island, Georgia, for the launchings of Apollo 13 and 14, 374 km to the south. Simultaneous acoustic and seismic waves were recorded for both events at times appropriate to the arrival of the acoustic waves from the source. Significant comparisons of the true signals are (1) the acoustic signal is relatively broadband compared to the nearly monochromatic seismic signal; (2) the seismic signal is much more continuous than the more pulse-like acoustic signal; (3) ground loading from the pressure variations of the acoustic waves is shown to be too small to account for the seismic waves; (4) the measured phase velocities of both acoustic and seismic waves across the local instrument arrays differ by less than 6 per cent and possibly 3 per cent if experimental error is included. It is concluded that the seismic waves are generated by resonant coupling to the acoustic waves along some 10 km of path on Skidaway Island. The thickness of unconsolidated sediment on the island is appropriate to a resonant ground wave frequency of 3.5 to 4 Hz, as observed. Under appropriate conditions, ground wave observations may prove more effective means of detecting certain aspects of acoustic signals in view of the filtering of wind noise and amplification through resonance.     

Caustics and focusing produced by sedimentary basins: applications of catastrophe theory to earthquake seismology

Summary. As high-frequency elastic waves propagate through real media, it is common for caustics and focusing to occur. Typically, rays may envelop a caustic surface in space, or exceptionally they may all coalesce at a focal point. In strong motion seismology, the observed large fluctuations in peak acceleration and intensity of ground shaking may just be a consequence of focusing and caustics created by waves propagating through irregularly shaped sedimentary basins. These basins, acting as deformed optical lenses, are capable of producing a complex network of patches and seemingly isolated pockets of intensified damage or high intensity shaking where the caustic intersects the ground surface. We adapt methods from optics and catastrophe theory to study the properties of caustics induced by typical sedimentary basins. Several hypothetical examples are shown that reflect the fact that these properties are useful to assess quantitatively the degree of wavefield amplification to be expected. A good correlation is found between actual damage patterns and caustic locations computed for the Caracas, Venezuela earthquake of 1967.     

Receiver functions at the stations of the German Regional Seismic Network (GRSN)

Several years of broad-band teleseismic data from the GRSN stations have been analysed for crustal structure using P -to- S converted waves at the crustal discontinuities. An inversion technique was developed which applies the Thomson-Haskell formalism for plane waves without slowness integration. The main phases observed are Moho conversions, their multiples in the crust, and conversions at the base of the sediments. The crustal thickness derived from these data is in good agreement with results from other studies. For the Gräfenberg stations, we have made a more detailed comparison of our model with a previously published model obtained from refraction seismic experiments. The refraction seismic model contains boundaries with strong velocity contrasts and a significant low-velocity zone, resulting in teleseismic waveforms that are too complicated as compared to the observed simple waveforms. The comparison suggests that a significant low-velocity zone is not required and that internal crustal boundaries are rather smooth.