Earth Motion Caused by Local Atmospheric Pressure Changes

Summary Observations have been made of the local atmospheric pressure field and the long-period seismic noise fields both on the surface of the Earth and in a mine at a depth of 183 metres. The observations show that during windy intervals and in the period range 20–100 s there is a strong correlation between local atmospheric pressure changes and the noise recorded by a vertical seismograph located on the surface. In contrast, over the same range of periods there is no correlation between the seismic noise recorded in the mine and local atmospheric pressure changes except during the passage of acoustic waves. It is shown that the noise in this pass band is not due to the buoyant response of the seismograph, but is caused by the motion of the Earth responding to atmospheric pressure changes.     

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.     

On reduction of long-period horizontal seismic noise using local barometric pressure

Tilt from atmospheric loading has long been known to be the major source of long-period horizontal seismic noise. We try to quantify these effects for seismic data from the Black Forest Observatory (BFO), which is known to be a very quiet station. Experimental transfer functions between local barometric pressure and horizontal seismic noise are estimated for two long time-series by standard methods. Two simple analytical physical models are developed: the local deformation model (LDM) and the acoustic-gravity wave model (TWM). Subsequently these models, with only two free parameters are fit using least squares to the observed seismic noise for time-series of widely differing lengths. The results are variable, sometimes rather dramatic variance reductions are obtained and sometimes the reduction is hardly significant. The method produces the best results when barometrically induced noise is high. The resulting admittances for the LDM are compared to finite element calculations. Since the methods are simple and can result in conspicuous reductions in noise we provide one more reason for installing barometers at even the best broad-band seismic stations.     

The utility of regional Chinese seismograms for source and path studies in central Asia

Summary. Broadband seismograms from the National Seismic Network of the People's Republic of China (PRC) have recently become available through a data exchange programme between NOAA and the State Seismological Bureau of the PRC. In this study, regional surface waves recorded at the Urumchi station located about 700 km north of the Tibetan Plateau in the Sinkiang Province are used to study East Kazakh explosions and wave propagation in central Asia. The data consist of broadband (flat to displacement between 0.1 and 10 Hz), photographic records from an SK Kirnos galvanometric system. Simultaneous inversion of Rayleigh wave phase and group velocities for the path from East Kazakh through the Dzhungarian Basin yields a crustal model dominated by the presence of very low velocities and a strong positive velocity gradient above 15 km depth. Velocities below 15 km depth are not significantly different from other continental structures underlain by Palaeozoic or Precambrian basement. Seismic moments were estimated for seven East Kazakh explosions using models of explosion sources with associated tectonic strain release. The largest explosion studied occurred on 1980 September 14 and had an m_b of 6.2 and a seismic moment of 2.7 × 10²³ dyn cm. The observed amplitude spectra of Rayleigh waves are richer in high frequencies than spectra calculated from our models. This could be caused by a path effect involving seismic wave focusing by the Dzhungarian Basin, although source medium effects cannot be ruled out.     

Flexure of the ocean lithosphere from island uplift, bathymetry and geoid height observations: the Society Islands

Summary. The response of a stratified elastic medium can be conveniently characterized by the Green's tensor for the medium. For coupled seismic wave propagation ( P—SV or fully anisotropic), the Green's tensor may be constructed directly from two matrices of linearly independent displacement solutions. Rather simple forms for the Green's tensor can be found if each displacement matrix satisfies one of the boundary conditions on the seismic field. This approach relates directly to 'reflection matrix' representations of the seismic field.
For a stratified elastic half space the Green's tensor is used to give a spectral representation for coupled seismic waves. By means of a contour integration a general completeness relation is obtained for the 'body wave' and 'surface wave' parts of the seismic field. This relation is appropriate for SH and P–SV waves in an isotropic medium and also for full anisotropy.     

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.     

Migration - why doesn't it work for deep continental data?

Summary. Conventional migration of deep seismic reflection data often produces disappointingly poor results even when the original unmigrated data are of high quality and are relatively noise free. Surprisingly, deep data often subjectively appear to be best migrated at velocities which are up to 50% less than appropriate interval velocities derived from crustal refraction experiments or directly from stacking velocities. The explanation for this behaviour is that near surface features distort and attenuate the seismic wave field and produce apparent discontinuities in deep reflections. Since these discontinuities are spurious they are not associated with the appropriate large diffractions which real discontinuities at depth would produce. During the process of migration reflections are invented in order to cancel out the missing diffractions thereby producing a "smiley" section which appears overmigrated. Since the lateral extent of individual smiles increases with increasing migration velocity, increasing two-way-time, and increasing seismic wavelength, the effect" is almost unnoticeable on conventional shallow seismic data but is overwhelming for deep crustal data.     

Double-planed deep seismic zone and upper-mantle structure in the Northeastern Japan Arc

Summary. The ScSp wave converted from the ScS wave at the boundary between the descending lithospheric slab and the mantle above it was clearly observed from a nearby deep earthquake with magnitude 7.7 at some stations of the seismic network of Tohoku University which covers the Tohoku District, the northeastern part of Honshu, Japan. By applying the three-dimensional seismic-ray tracing method, the location of this boundary was determined from the difference in arrival time between the ScS and ScSp waves. The result shows that the upper boundary of the descending slab lies exactly on the upper plane of the double-planed deep seismic zone found in the Northeastern Japan Arc.
There is an additional evidence that the boundary is located on the upper plane of the double-planed deep seismic zone. The hypocentre distribution of intermediate-depth earthquakes located by the small-scale seismic-array observation is extremely different from that obtained by the relatively large-scale seismic network. The discrepancy in the distribution of hypocentres of the same earthquake independently located is well explained by the inclined lithospheric slab model derived from the difference in arrival time between the ScS and ScSp waves.
The earthquakes with reverse faulting or with down-dip compressional stresses occur at the upper boundary of the descending slab. Within the descending slab, the earthquakes with down-dip extensional stresses also occur in a very narrow zone from 30 to 40 km below the dipping boundary in the depth range from 50 to about 200 km, and these shocks form the lower plane of the double-planed deep seismic zone.     

Anisotropy of the mantle inferred from observations of P to S converted waves

Summary. Seismic anisotropy has been previously studied at depths usually not exceeding 100 or 150 km. In this paper we present a method of analysis of seismic records which is very sensitive to azimuthal anisotropy and is applicable at almost any depth range. The idea of the method is to detect and analyse the SH -component of the waves, converted from P to S in the mantle. The procedure of record processing includes frequency filtering, axis rotation, transformation of the record to a standard form, stacking the standardized SH -component records of many seismic events, and the harmonic analysis of amplitude as a function of the direction of wave propagation. When applied to the long-period records of NORSAR the procedure detected a converted wave with the properties implying the possibility of its propagation in a transversely isotropic medium with a horizontal axis of symmetry . Our preferred model postulates anisotropy of ∼ 1 per cent in a layer 50 km thick at the base of the upper mantle.     

Array analysis using circular-wave-front geometry:an application to locate the nearby seismo-volcanic source

The zero-lag cross-correlation technique, used for array analysis in the hypothesis of plane waves, has been modified to allow the wave front to be circular. Synthetic tests have been performed to check the capability of the method, which returns the input test data when the source–array distances are not greater than two or three times the array aperture. For this distance range the method furnishes an estimate of the apparent velocity and the epicentral coordinates of the source. For more distant sources the method becomes equivalent to that based on the planar-wave approximation, which gives an estimate of the backazimuth to the source and the apparent velocity. The method has been applied to seismic data recorded at the active volcano located at Deception Island, Antarctica. 35 volcanic long-period events occurring in a small swarm were selected. Results show that the epicentres are close to the array (between 0.4 and 2 km) and aligned in a SW direction, in agreement with one of the main directions of the fracture system of Deception volcano.     

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.     

Traveltime measurements from noise correlation: stability and detection of instrumental time-shifts

We test the feasibility of using Green's functions extracted from records of ambient seismic noise to monitor temporal changes in the Earth crust properties by repeated measurements at regional distances. We use about 11 yr of continuous recordings to extract surface waves between three pairs of stations in California. The correlations are computed in a moving 1-month window and we analyse the temporal evolution of measured interstation traveltimes. The comparison of the arrival times in the positive and negative correlation time of Rayleigh and Love waves allows us to separate time-shifts associated with any form of physical change in the medium, those resulting from clock drift or other instrumental errors, and those due to change in the localization of the noise sources. This separation is based on the principle of time symmetry. When possible, we perform our analysis in two different period bands: 5–10 and 10–20 s. The results indicate that significant instrumental time errors (0.5 s) are present in the data. These time-shifts can be measured and tested by closure relation and finally corrected independently of any velocity model. The traveltime series show a periodic oscillation that we interpret as the signature of the seasonal variation of the region of origin of the seismic noise. Between 1999 and 2005, the final arrival time fluctuations have a variance of the order of 0.01 s. This allows us to measure interstation traveltimes with errors smaller than 0.3 per cent of the interstation traveltime and smaller than 1 per cent of the used wave period. This level of accuracy was not sufficient to detect clear physical variation of crustal velocity during the considered 11 yr between the three stations in California. Such changes may be more easily detectable when considering pairs of stations more closely located to each other and in the vicinity of tectonically active faults or volcanoes.     

Surface-wave propagation across the Mexican Volcanic Belt and the origin of the long-period seismic-wave amplification in the Valley of Mexico

A network of nine broad-band seismographs was operated from March to May 1994 to study the propagation of seismic waves across the Mexican Volcanic Belt (MVB) in the region of the Valley of Mexico. Analysis of the data from the network reveals an amplification of seismic waves in a wide period band al the stations situated in the southern part of the MVB.
The group velocities of the fundamental mode of the Rayleigh wave in the period range 2–13 s are found to be lower in the southern part of the MVB than in its northern part and in the region south of the MVB. The inversion of dispersion curves shows that the difference in group velocities is due to the presence of a superficial low-velocity layer (with an average S -wave velocity of 1.7 km s^-1 and an average thickness of 2 km) beneath the southern part of the MVB. This low-velocity zone is associated with the region of active volcanism.
Numerical simulations show that this superficial low-velocity layer causes a regional amplification of 8–10 s period signals, which is of the same order as the amplification measured from the data. This layer also increases the signal duration significantly because of the dispersion of the surface waves. These results confirm the hypothesis of Singh et al. (1995), who suggested that the regional amplification observed in the Valley of Mexico is due to the anomalously low shear-wave velocity of the shallow volcanic rocks in the southern MVB     

Transitional continental–oceanic structure beneath the Norwegian Sea from inversion of surface wave group velocity data

We have analysed the fundamental mode of Love and Rayleigh waves generated by 12 earthquakes located in the mid-Atlantic ridge and Jan Mayen fracture zone. Using the multiple filter analysis technique, we isolated the Rayleigh and Love wave group velocities for periods between 10 and 50  s. The surface wave propagation paths were divided into five groups, and average group velocities calculated for each group. The average group velocities were inverted and produced shear wave velocity models that correspond to a quasi-continental oceanic structure in the Greenland–Norwegian Sea region. Although resolution is poor at shallow depth, we obtained crustal thickness values of about 18  km in the Norwegian Sea area and 9  km in the region between Svalbard and Iceland. The abnormally thick crust in the Norwegian Sea area is ascribed to magmatic underplating and the thermal blanketing effect of sedimentary layers. Maximum crustal shear velocities vary between 3.5 and 3.9  km  s⁻¹ for most paths. An average lithospheric thickness of 60  km was observed, which is lower than expected for oceanic-type structure of similar age. We also observed low shear wave velocities in the lower crust and upper mantle. We suggest that high heat flow extending to depths of about 30  km beneath the surface can account for the thin lithosphere and observed low velocities. Anisotropy coefficients of 1–5 per cent in the shallow layers and >7 per cent in the upper mantle point to the existence of polarization anisotropy in the region.     

Resolution potential of surface wave phase velocity measurements at small arrays

The deployment of temporary arrays of broadband seismological stations over dedicated targets is common practice. Measurement of surface wave phase velocity across a small array and its depth-inversion gives us information about the structure below the array which is complementary to the information obtained from body-wave analysis. The question is however: what do we actually measure when the array is much smaller than the wave length, and how does the measured phase velocity relates to the real structure below the array? We quantify this relationship by performing a series of numerical simulations of surface wave propagation in 3-D structures and by measuring the apparent phase velocity across the array on the synthetics. A principal conclusion is that heterogeneities located outside the array can map in a complex way onto the phase velocities measured by the array. In order to minimize this effect, it is necessary to have a large number of events and to average measurements from events well-distributed in backazimuth. A second observation is that the period of the wave has a remarkably small influence on the lateral resolution of the measurement, which is dominantly controlled by the size of the array. We analyse if the artefacts created by heterogeneities can be mistaken for azimuthal variations caused by anisotropy. We also show that if the amplitude of the surface waves can be measured precisely enough, phase velocities can be corrected and the artefacts which occur due to reflections and diffractions in 3-D structures greatly reduced.     

Seismic coda due to non-linear elasticity

Non-linear elastic response of rocks has been widely observed in laboratory, but very few seismic studies are reported in the literature, even though it is the most natural environment where this feature could be observed. Analytic solutions to the non-linear wave propagation phenomena are not readily available, and there is a need to use approximated techniques. It is clear that when a seismic wave propagates through a homogeneous non-linear elastic media, it will be perturbed by the non-linearity. This perturbation can be treated as a source of scattering, spreading the energy of the primary wave in space and time, contributing to the seismic coda. This is in some sense similar to the effect of heterogeneities. The properties of the coda due to the non-linearity depend on the amount of non-linearity and the seismic moment. Using a perturbation approach we calculate the amplitude of the scattered waves, and show that it can describe reasonably well the main features of real seismic codas.     

Regional variations of higher Rayleigh-mode phase velocities: a spatial-filtering method

Summary. Application of both spatial filtering and multiple-frequency filtering techniques allows one to isolate a higher-mode surface wave and to find the regional variation of its dispersion. The method is applied to four sets of long-period records across the United States for intermediate earthquakes located in the New Hebrides. The first and second higher Rayleigh modes show lower phase velocities in the western part of the States than in the eastern part. Furthermore, an oceanic phase-velocity curve is determined between 50 and 130 s for the first higher Rayleigh mode on the full path across the Pacific Ocean.     

Observations of Infrasound and Subsonic Disturbances Related to Severe Weather

Summary During the past 10 years the Geoacoustics Group of NOAA's Wave Propagation Laboratory studied travelling low-frequency pressure variations related to thunderstorms and severe weather. Two general categories of waves were associated with severe weather conditions: 'subsonic' pressure disturbances and infrasonic waves with acoustic velocities. The low-frequency pressure variations were measured at the Earth's surface using microphone arrays located at times thousands of kilometres from the severe-weather disturbance. The radiated infra- sound was related to thunderstorms penetrating the tropopause and spectral analyses were performed on several signals. Possible practical applications to storm warning and classification are discussed for both infrasound and 'subsonic' pressure disturbances. Past measurements of these signals are reviewed.     

Horizontal layered structure with heterogeneity beneath continents and island arcs from particle orbits of long-period P waves

We investigate the particle orbits of long-period (about 20 s) P waves observed with the global seismic network. By analysing 84 three-component seismograms recorded at 25 stations from 60 earthquakes occurring beneath 300 km, we quantitatively evaluate the orbits by three sets of eigenvalues and eigenvectors, using a covariance matrix method. The eigenvalues for P waves recorded at stations located on continents are explained by the standard horizontal layered structure model (iasp91). On the other hand, the orbits observed at stations close to island arcs are affected not only by the horizontal layered structure but also by heterogeneity due to subducting plates, mantle diapirs and so on. On the basis of a single-scattering model for a plane P wave, we quantify the heterogeneities by an isotropic scattering coefficient g0. Fitting the theoretical eigenvalues to the observed ones, we estimate g0 for the crust and upper mantle beneath continents to be less than 0.0005 km^-1, and the mean g₀ for the structure beneath island arcs to be about 0.0015 to 0.003 km^-1.     

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.