Excitation mechanism of atmospheric pressure waves from the 1980 Mount St Helens eruption

Summary. Atmospheric pressure waves from the Mount St Helens eruption 1980 May 18 have been clearly recorded by a sensitive microbarograph at Berkeley, California. The record shows three types of waves with different group velocities. The pressure waves can be interpreted in terms of direct waves A1, antipodean travelling waves A2 and circumnavigating waves A3, all of which are composed of several acoustic-gravity modes propagated in the lower atmosphere. Synthetic barograms appropriate to the Berkeley station have been calculated on the basis of the dynamic response of the lower atmospheric structure, together with various assumptions of source properties. Comparisons between synthetic and observed barograms provide estimates for ranges of the time history of upward particle velocity at the source, source dimensions and the velocity of the source spreading over the blast zone, as well as for the average dissipation effects over the circumferential path. The results suggest that two major compression pulses on the A1 record correlate with the arrival of pressure waves from the first (lateral) blast and second (vertical) blast, although the inferred interblast time interval is not consistent with that estimated from seismic observations.     

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.     

Guided wave propagation in laterally varying media - I. Theoretical development

Summary. The propagation of surface waves in a laterally varying medium can be described by representing the wavetrain as a superposition of modal contributions for a reference structure. As the guided waves propagate through a heterogeneous zone the modal coefficients needed to describe the wavetrain vary with position, leading to interconversions between modes and reflection into backward travelling modes. The evolution of the modal terms may be described by a set of first-order differential equations which allow for coupling to both forward and backward travelling waves; the coefficients in these equations depend on the differences between the actual structure and the reference structure. This system is established using the orthogonality properties of the modal eigenfunctions and is valid for SH -waves, P - SV -waves and full anisotropy.
The reflected and transmitted wavefields for a region of heterogeneity can be related to the incident wave by introducing reflection and transmission matrices which connect the modal coefficients in these fields to those in the incident wavetrain. By considering a sequence of models with increasing width of heterogeneity we are able to derive a set of Ricatti equations for the reflection and transmission matrices which may be solved by initial value techniques. This avoids an awkward two-point boundary value problem for a large number of coupled equations. The method is demonstrated for 1 Hz Lg - and Sn -waves in a multilayered model for which there are 19 coupled modes.
The method is applicable to three-dimensional heterogeneity, and we are able to show that the interconversion between Love and Rayleigh waves, in the presence of gradients in seismic properties transverse to the propagation path, leads to a net rate of increase of the transverse components of the seismogram at the expense of the other components.     

Experimental research on acoustic wave velocity of frozen soils during the uniaxial loading process

Ultrasonic P-wave tests of frozen silt and frozen sand were conducted during uniaxial loading by using an RSM^®-SY5(T) nonmetal ultrasonic test meter to study the velocity characteristics of P-waves. The experimental results indicate that the P-wave velocity is affected by soil materials, temperature, and external loads, so the P-wave velocity is different in frozen silt and frozen sand, but all decrease with an increase of temperature and increase at first and then decrease with strain during the loading process. There is an exponential relationship between uniaxial compressive strength and P-wave velocity, and the correlation between them is very good. The characteristic parameters of acoustic waves can, to some extent, reflect the development of internal cracks in frozen soils during loading.     

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.     

Seismic wave properties in time-dependent porosity homogeneous media

It is quantified the properties of seismic waves in fully saturated homogeneous porous media within the framework of Sahay's modified and reformulated poroelastic theory. The computational results comprise amplitude attenuation, velocity dispersion and seismic waveforms. They show that the behaviour of all four waves modelled as a function of offset, frequency, porosity, fluid viscosity and source bandwidth depicts realistic dissipation within the sonic–ultrasonic band. Therefore, it appears that there is no need to include material heterogeneity to model attenuation. By inference it is concluded that the fluid viscosity effects may be enhanced by dynamic porosity.     

Acoustoelasticity in cracked solids

A new model that accounts for the stress dependence of the phase velocity of elastodynamic waves propagating in a cracked solid under compression is presented. The phase velocities of longitudinal and shear waves are derived from the effective elastic properties of a cracked solid, which are evaluated within the framework of Kachanov's approach. Following Kachanov, the extra-compliance tensor of the cracked solid is related to the crack compliances, which display a marked non-linear behaviour when subjected to a compressive load. Such non-linear behaviour is shown to be derived from the elastic interaction between the contacting crack faces under compression. This work does not address the effect of mutual interaction among cracks and the generation of higher harmonics due to the medium non-linearity. Numerical examples are presented that illustrate the phase velocity changes occurring in a solid with a random distribution of parallel cracks as a function of an external compressive load. A distinctive feature of the acoustoelastic effect in solids with large parallel fractures and in solids with distributions of aligned microcracks is also illustrated.     

Stress drops, wave spectra and recurrence intervals of great earthquakes — implications of the Etorofu earthquake of 1958 November 6

Summary . The great Etorofu earthquake of 1958 November 6 is characterized by a relatively small aftershock area (70 × 150 km²) and an extremely large felt area. The felt area is more extensive than those of any other large earthquakes which have occurred in the southern Kurile to northern Japan arc since the beginning of this century. The mechanism is a pure thrust fault typical of most great earthquakes in island arcs. A body wave magnitude of m _b = 8.2 is obtained at periods around 6 s using more than 40 observations, although an m _b value of only 7.6–7.7 would be expected empirically from the observed surface wave magnitude of M _s= 8.1–8.2. Both an unusually large felt area and a high m _b indicate a dominance of high-frequency components in the seismic waves. A seismic moment of M _o= 4.4 × 10²⁸ dyne cm is determined from long-period surface waves from which a high stress drop of Δσ = 78 bar is obtained using a relatively small aftershock area. Historic data indicate an anomalously long time interval between the 1958 event and any earlier great earthquake from the same source region. The observed high stress drop can be interpreted as a consequence of this long intervening period through which strain built up. The dominance of the high-frequency seismic waves can then be interpreted as a result of this high stress drop. Stress drops, seismic wave spectra and recurrence intervals of great earthquakes are in this way closely related to each other. The 1958 event may represent a high strength extreme of stochastic fluctuation of fracture strength relevant to great earthquakes.     

Geometrical structure of shear wave surfaces near singularity directions in anisotropic media

There are three types of surfaces which are used for studying wave propagation in anisotropic media: normal surfaces, slowness surfaces and wave surfaces. Normal surfaces and slowness surfaces have been researched in detail. Wave surfaces are the most complicated and comparatively poorly known compared with the other two. Areas of complicated geometrical structure of the wave surfaces are located in the vicinity of conical acoustic axes. There is an elliptical hole on the quick shear wave surface and complicated folds and cusps on the slow shear wave surface. Decomposition of the slow shear wave surface into smooth sheets is used for the study of its geometrical structure. Complexity of shear wave surfaces can be expressed by the number of waves corresponding to a fixed ray. An original approach to the calculation of wave normals depending on ray direction is presented.     

Propagation of harmonic plane waves in a general anisotropic porous solid

Wave propagation is studied in a general anisotropic poroelastic solid. The presence of dissipation due to fluid-viscosity as well as hydraulic anisotropy of pore permeability are also considered. Biot's theory is used to derive a system of modified Christoffel equations for the propagation of plane harmonic waves in porous media. A non-trivial solution of this system is ensured by a determinantal equation. This equation is separated into two different polynomial equations. One is the quartic equation whose roots represent the complex velocities of four attenuating waves in the medium. The other is a eighth-degree polynomial whose roots represent the vertical slowness values for the four waves propagating upward and downward in a finite porous medium. Procedure is explained to associate the numerically obtained roots with the waves propagating in the medium. The slowness surfaces of waves reflected at the boundary of the medium are computed for a realistic numerical model. The behaviours of phase velocity surfaces are analysed with the help of numerical examples.     

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.     

Experimental Noise Reducers for an Active Microbarograph Array

Summary Brief outlines are given of the problem of identifying low frequency acoustic waves from atmospheric background noise and of the methods used by UKAEA to improve detection of these waves.
An important aspect of the experimental work has been to develop devices to reduce atmospheric pressure noise, suitable for use at existing microbarograph recording sites, and to evaluate their performance in field tests. Some interim results from these experiments are given.     

Seismo-acoustic analysis of the Buncefield oil depot explosion in the UK, 2005 December 11

A massive vapour cloud explosion occurred at the Buncefield fuel depot near Hemel Hempstead, UK, in the morning of 2005 December 11. The explosion was the result of an overflow from one of the storage tanks with the release of over 300 tons of petrol and generating a vapour cloud that spread over an area of 80 000 m², before being ignited. Considerable damage was caused in the vicinity of the explosion and a total of 43 people were injured. The explosion was detected by seismograph stations in the UK and the Netherlands and by infrasound arrays in the Netherlands. We analysed the seismic recordings to determine the origin time of 06:01:31.45 ±0.5 s (UTC) from P -wave arrival times. Uncertainties in determination of origin time from acoustic arrival times alone were less than 10 s. Amplitudes of P -, Lg and primary acoustic waves were measured to derive decay relationships as function of distance. From the seismic amplitudes we estimated a yield of 2–10 tons equivalent to a buried explosion. Most seismic stations recorded primary and secondary acoustic waves. We used atmospheric ray tracing to identify the various travel paths, which depend on temperature and wind speed as function of altitude, leading to directional variation. Refracted waves were observed from the troposphere, stratosphere and thermosphere with a good match between observed and calculated traveltimes. The various wave types were also identified through array processing, which provides backazimuth and slowness, of recordings from an infrasound array in the Netherlands. The amplitude of stratospheric refracted acoustic waves recorded by the array microbarometers was used to estimate a yield of 21.6 (±5) tons TNT equivalent. We have demonstrated through joint seismo-acoustic analysis of the explosion that both the seismic velocity model and the atmospheric model are sufficient to explain the observed traveltimes.     

Scattering of waves in elastic media containing passive and active cracks

Scattering of wavefields in a 3-D medium that includes passive and/or active structures, is numerically solved by using the boundary integral equation method (BIEM). The passive structures are velocity anomalies that generate scattered waves upon incidence, and the active structures contain endogenous fracture sources, which are dynamically triggered by the dynamic load due to the incident waves. Simple models are adopted to represent these structures: passive cracks act as scatterers and active cracks as fracture sources. We form cracks using circular boundaries, which consist of many boundary elements. Scattering of elastic waves by the boundaries of passive cracks is treated as an exterior problem in BIEM. In the case of active cracks, both the exterior and interior problems need to be solved, because elastic waves are generated by fracturing with stress drop, and the growing crack boundaries scatter the incident waves from the outside of the cracks. The passive cracks and/or active cracks are randomly distributed in an infinite homogeneous elastic medium. Calculations of the complete waveform considering a single scatter show that the active crack has weak influence on the attenuation of first arrivals but strong influence on the amplitudes of coda waves, as compared with those due to the passive crack. In the active structures, multiple scattering between cracks and the waves triggered by fracturing strongly affect the amplitudes of first arrivals and coda waves. Compared to the case of the passive structures, the attenuation of initial phase is weak and the coda amplitudes decrease slowly.     

Impact of storms on mixed carbonate and siliciclastic shelves: insights from combined diffusive and fluid-flow transport stratigraphic forward model

A quantitative stratigraphic model of mixed carbonate/siliciclastic continental shelves is presented to investigate the relationships between depositional processes and stratigraphic responses at long‐term, large spatial scales. A diffusion model is combined with a fluid‐flow approach to simulate both long‐term factors, i.e. the processes controlling large‐scale architecture, and short‐term processes, i.e. sediment redistribution by storms. Any net sediment accumulation is the result of the succession of a storm and a fair‐weather period. Sediments are mobilized by waves and advected by low‐frequency currents during storm events. Sediments are then reworked and redistributed downslope by diffusive processes during fair‐weather period. The results are successful in capturing several major characteristics of both modern and ancient depositional systems (geometry, differential preservation, net accumulation rates). The study highlights the importance of waves and unidirectional currents. Depositional geometry and shelf morphology depend on the balance between available sediment supply (generated in situ or detrital) and the transport energy, which is related to the style of sediment transport (diffusive or advective), and to the magnitude and frequency of storms.     

On the relative scattering of P- and S-waves

Summary. Using a single scattering approximation, we derive equations for the scattering attenuation coefficients of P- and S -body waves. We discuss our results in the light of some recent energy renormalization approaches to seismic wave scattering. Practical methods for calculating the scattering attenuation coefficients for various earth models are emphasized. The conversions of P - to S -waves and S- to P -waves are included in the theory. The earth models are assumed to be randomly inhomogeneous, with their properties known only through their average wavenumber power spectra. We approximate the power spectra with piecewise constant functions, each segment of which contributes to the net, frequency-dependent, scattering attenuation coefficient. The smallest and largest wavenumbers of a segment can be plotted along with the wavevectors of the incident and scattered waves on a wavenumber diagram. This diagram gives a geometric interpretation for the frequency behaviour associated with each spectral segment, including a 'transition' peak that is due entirely to the wavenumber limits of the segment. For regions of the earth where the inhomogeneity spectra are concentrated in a band of wavenumbers, it should be possible to observed such a peak in the apparent attenuation of seismic waves. We give both the frequency and distance limits on the accuracy of the theoretical results.     

The Acoustics of Pipe Arrays

Summary A theoretical analysis is given for the acoustical behaviour of the pipe-microbarograph systems used to detect acoustic gravity waves and other modes of infrasound. It is shown how to compute the response of the microbarograph to a fluctuating pressure at any one inlet port of the pipe and how the results of such computations may be used to calculate the response to a plane sound wave traversing the system.
The analysis is illustrated by numerical examples obtained by means of a computer program. These examples confirm that the tapered tube modelled after Daniels' line microphone has very good characteristics, but that good results may also be obtained using pipes of uniform bore. The work leans heavily on Benade's calculations of sound propagation in a circular conduit.     

Observations and origin of Rayleigh-wave amplitude anomalies

This is a report of observations of amplitude anomalies of fundamental-mode Rayleigh waves ( R 1) between periods of 17 and 100  s. The anomalies are with respect to amplitudes predicted by Rayleigh-wave excitation for a reference earth model and catalogued centroid earthquake source parameters, such as are used in large-scale waveform inversions. The observations indicate that the amplitude anomalies are consistent for nearby recordings of the same event, while there is no obvious relation between the observed anomalies and the paths travelled by the waves. This is in contrast to Rayleigh-wave phase anomalies, which are consistent for similar propagation paths, and hence form the input in many inversions for along-path structure. The observations in this paper show that a similar inversion of intermediate-period amplitude anomalies for along- and near-path structure is not warranted without eliminating source effects, since the amplitude anomalies are dominated by scattering off near-source earth structure and by possible uncertainties in the source parameters. Sensitivity kernels that take the coupling between the moment tensor and displacement field into account demonstrate that Rayleigh-wave amplitude sensitivity is largest near the source. This report argues that the interaction between source-radiated Rayleigh waves and near-source earth structure may not be ignored in amplitude inversion procedures.     

塔克拉玛干沙漠腹地湍流能量耗散率和结构函数参数特征

利用塔中气象站垂直探测系统的湍流观测资料,计算并分析了塔克拉玛干沙漠腹地不同稳定层结条件的湍涡特征长度尺度、能量耗散率及湍流结构函数参数特征。结果表明：湍涡特征长度尺度在弱不稳定或近中性条件时最大,随不稳定程度的增强有明显减小的趋势,随着稳定程度的增加有先迅速减小后又缓慢增加的趋势,且热量特征长度尺度比动量特征长度尺度总体要大一个数量级;无因次湍流热量耗散率的拟合函数形式与其他试验一致,但经验系数略有差异;利用Kaimal公式和湍流谱方法间接计算的稳定条件湍流能量耗散率及湍流结构函数参数具有较好的一致性。     

Finite-frequency sensitivity kernels for head waves

Head waves are extremely important in determining the structure of the predominantly layered Earth. While several recent studies have shown the diffractive nature and the 3-D Fréchet kernels of finite-frequency turning waves, analogues of head waves in a continuous velocity structure, the finite-frequency effects and sensitivity kernels of head waves are yet to be carefully examined. We present the results of a numerical study focusing on the finite-frequency effects of head waves. Our model has a low-velocity layer over a high-velocity half-space and a cylindrical-shaped velocity perturbation placed beneath the interface at different locations. A 3-D finite-difference method is used to calculate synthetic waveforms. Traveltime and amplitude anomalies are measured by the cross-correlation of synthetic seismograms from models with and without the velocity perturbation and are compared to the 3-D sensitivity kernels constructed from full waveform simulations. The results show that the head wave arrival-time and amplitude are influenced by the velocity structure surrounding the ray path in a pattern that is consistent with the Fresnel zones. Unlike the 'banana–doughnut' traveltime sensitivity kernels of turning waves, the traveltime sensitivity of the head wave along the ray path below the interface is weak, but non-zero. Below the ray path, the traveltime sensitivity reaches the maximum (absolute value) at a depth that depends on the wavelength and propagation distance. The sensitivity kernels vary with the vertical velocity gradient in the lower layer, but the variation is relatively small at short propagation distances when the vertical velocity gradient is within the range of the commonly accepted values. Finally, the depression or shoaling of the interface results in increased or decreased sensitivities, respectively, beneath the interface topography.