Crust and upper mantle shear velocities from controlled sources

b
A two ship refraction profile was undertaken on the Australian continental shelf during the Banda Sea geophysical program, carried out by the Woods Hole Oceanographic Institution, the Scripps Institution of Oceanography and the Geological Survey of Indonesia. S waves originating close to the sea bottom were observed to distances of up to 1150 km at an array of stations in northern Australia.
These observations are interpreted as implying S mantle velocities of 4.60 km s^-1 from a depth of 45 km to a depth of 76 km and 4.72 km s^-1 below a depth of 76 km.
Ratios of the P and S travel times (V_p/V_s) have been determined to be 1.74 in the crust rising to a value of greater than 1.79 below a velocity discontinuity at a depth of 200 km. It is inferred that this high value arises because the effect of temperature is greater for S than for P .
Using the data from this and other studies in the shield region of Northern Australia it has been found that the _S travel times are significantly less than predicted by the Jeffreys—Bullen tables.     

Travel times of Pn -waves in the Aegean and surrounding area

Summary. Based on accurately located 23 very shallow earthquakes ( h = 1–14 km) in northern and central Greece by portable networks of seismic stations and by the joint epicentre method, the travel times of the P_n -waves from the foci of these earthquakes to the sites of 54 permanent stations in the Balkan region have been determined. The travel times of P_n -waves in the central and eastern part of the area (eastern Greece, south-eastern Yugoslavia, the Aegean Sea, Bulgaria, southern Romania, western Turkey) fit a straight line very well with the P_n velocity equal to 7.9 ± 0.1 km s^-1. On the contrary, the travel times of P_n -waves to stations in the western part of the area (Albania, western Greece) do not fit this curve because the P_n -waves travelling to these stations are delayed by more than 1 s due to the thicker crust under the Dinarides–Hellenides mountain range. Time delays for P_n -waves have been calculated for each permanent station in the Balkan area with respect to the mean travel-time curve of these waves in the central and eastern part of the area. Corrections of the travel times for these delays contribute very much to the improvement of the accuracy in the location of the shallow earthquakes in the Aegean and surrounding area.     

Detailed S-wave structure in the Dublin Basin and its northern margin

Summary. The seismic structure has been measured to a depth of about 3 km along a 30 km seismic profile in east central Ireland. This profile is unusual in that it is the S -wave velocity—depth structure that has been measured to a degree of precision more normally associated with P -wave results. One reason for this is that the sources used were quarry blasts which generated strong S -waves and short-period surface waves but rather weak P -waves.
The results show a layer of Carboniferous limestone with shear velocity 2.65 km⁻¹ s overlying a layer with a velocity of 3.06 km s⁻¹. This second layer was interpreted as Lower Palaeozoic strata (Silurian/Ordovician) since this velocity was evident in an inlier seen at the surface at the northern end of the line. A third refraction horizon, shear velocity 3.45 km s⁻¹ and displaying a basinal structure, was also recognized. This may be Cambrian or Precambrian basement.     

The traveltime perturbations for seismic body waves in factorized anisotropic inhomogeneous media

The traveltime perturbation equations for the quasi-compressional and the two quasi-shear waves propagating in a factorized anisotropic inhomogeneous (FAI) media are derived. The concept of FAI media simplifies considerably these equations. In the FAI medium, the density normalized elastic parameters a _ijkl ( X _i ) can be described by the relation a _ijkl ( X _i) = f ²( x _i ) A _ijkl, where A _ijkl are constants, independent of coordinates x _i and f ²( x _i) is a continuous smooth function of x _i . The types of anisotropy ( A _ijkl ) and inhomogeneity [ f ( x _i)] are not restricted. The traveltime perturbations of individual seismic body waves ( q ^P , qS 1 and qS 2) propagating in the FAI medium depend, of course, both on the structural pertubations [δ f ²( x _i)] and on the anisotropy perturbations (δ A _ijkl ), but both these effects are fully separated. The perturbation equations for the time delay between the two qS -waves propagating in the FAI medium are simplified even more. If the unperturbed (background) medium is isotropic, the perturbation of the time delay does not depend on the structural perturbations (δ f ²( x _i) at all. This striking result, valid of course only in the framework of first-order perturbation theory, will simplify considerably the interpretation of the time delay between the two split qS -waves in inhomogeneous anisotropic media. Numerical examples are presented.     

Anisotropy in the oceanic lithosphere of the North-western Pacific Basin

Summary. A series of long-range explosion seismological experiments has been conducted by the use of specially designed ocean bottom seismographs (OBSs) in the Western Pacific. OBS studies of apparent velocity measurements by the use of natural earthquakes have also been made. The experiments have made clear that large-scale P -wave anisotropy exists in the entire thickness of the oceanic lithosphere. The existence of the large-scale anisotropy in the oceanic lithosphere has been demonstrated for the first time by seismic body-wave studies. Previously, anisotropy had been found only in the uppermost oceanic mantle in the Eastern Pacific.
The azimuth of the maximum velocity, 8.6 km s^-1, is about 155° clock-wise from north. The direction is perpendicular to the magnetic lineation of the region, however, the direction differs from the direction of the present plate motion by about 30°. So it appears that the anisotropy has been 'frozen' at least since the change of the plate motion that occurred 40 Myr ago. The frozen anisotropy should set important constraints on the mechanical properties of the lithosphere such as the viscosity and temperature of the lower lithosphere.     

Crustal structure of the Faeroe–Shetiand Channel

Summary. The deep structure of the Faeroe–Shetland Channel has been investigated as part of the North Atlantic Seismic Project. Shot lines were fired along and across the axis of the Channel, with recording stations both at sea and on adjacent land areas. At 61°N, 1.7 km of Tertiary sediments overlies a 3.9–4.5 km s^-1 basement interpreted as the top of early Tertiary volcanics. A main 6.0–6.6 km s^-1 crustal refractor interpreted as old oceanic crust occurs at about 9 km depth. The Moho (8.0 ° 0.2 km s^-1) is at about 15–17 km depth. There is evidence that P _n may be anisotropic beneath the Faeroe–Shetland Channel. Arrivals recorded at land stations show characteristics best explained by scattering at an intervening boundary which may be the continent–ocean crustal contact or the edge of the volcanics.
The Moho delay times at the shot points, determined by time-term analysis, show considerable variation along the axis of the Channel. They correlate with the basement topography, and the greatest delays occur over the buried extension of the Faeroe Ridge at about 60° 15'N, where they are nearly 1 s more than the delays at 61°N after correction for the sediments. The large delays are attributed to thickening of the early Tertiary volcanic layer with isostatic downsagging of the underlying crust and uppermost mantle in response to the load, rather than to thickening of the main crustal ayer.
The new evidence is consistent with deeply buried oceanic crust beneath the Faeroe–Shetland Channel, forming a northern extension of Rockall Trough. The seabed morphology has been grossly modified by the thick and laterally variable pile of early Tertiary volcanic rocks which swamped the region, accounting for the anomalous shallow bathymetry, the transverse ridges and the present narrowness of the Channel.     

Earthquake intensity attenuation pattern in the United States

Summary. Attenuation of earthquake intensities with epicentral distance was studied by analysing the intensity data for 39 earthquakes in the United States. Attenuation of MM intensity ( I ) with distance (Δ km) follows a simple relation of the type log I = log I ₀ - m Δ, where I ₀ is the intensity at the epicentre and m is a constant. Slope m is found to be inversely proportional to the square of the focal depth. Intensity attenuation pattern in the United States in general can be represented by a unified relation I/I ₀= exp [-(0.8999/ h ²+ 0.0014)Δ] where 16km ≤ h ≤ 60km. Intensities were calculated with the help of this equation and a good agreement with the observed intensities were found. A comparative study has also been made between the attenuation relations applicable to India and the United States.     

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     

Comparison of the Upper Mantle Shear Wave Velocity Structure of the Indian Shield and the Tibetan Plateau and Tectonic Implications

Summary. Travel times and waveforms of long-period SH -waves recorded at distances of 10–30° and some SS waveforms are used to constrain the upper mantle velocities down to a depth of 400km beneath both the Indian Shield and the Tibetan Plateau. the shear velocity in the uppermost mantle beneath both the Indian Shield and the Tibetan Plateau is high and close to 4.7 km s⁻¹. the Indian Shield has a fairly thick high velocity lid, and the mean velocity between 40 and 250 km is between 4.58 and 4.68 km s⁻¹. In contrast, S -wave travel times and waveforms of S -waves, as well as a few for SS , show that the mean velocity between 70 and 250km beneath the central and northern part of the Tibetan Plateau is slower by 4 per cent or more than that beneath the Indian Shield and probably is between 4.4 and 4.5km s⁻¹. No large differences in the structure of the two areas below 250 km are required to explain both the arrival times and the waveforms of SH phases crossing Tibet or the Indian Shield. These results show that the structure of Tibet is not that of a shield and imply that the Indian plate is not underthrusting the whole of the Tibetan Plateau at the present time.     

Seismic anisotropy within the uppermost mantle of southern Germany

This paper presents an updated interpretation of seismic anisotropy within the uppermost mantle of southern Germany. The dense network of reversed and crossing refraction profiles in this area made it possible to observe almost 900 traveltimes of the P_n phase that could be effectively used in a time-term analysis to determine horizontal velocity distribution immediately below the Moho. For 12 crossing profiles, amplitude ratios of the P_n phase compared to the dominant crustal phase were utilized to resolve azimuthally dependent velocity gradients with depth. A P -wave anisotropy of 3–4 per cent in a horizontal plane immediately below the Moho at a depth of 30 km, increasing to 11 per cent at a depth of 40 km, was determined. For the axis of the highest velocity of about 8.03 km s⁻¹ at a depth of 30 km a direction of N31°F was obtained. The azimuthal dependence of the observed P_n amplitude is explained by an azimuth-dependent sub-Moho velocity gradient decreasing from 0.06 s⁻¹ in the fast direction to 0 s⁻¹ in the slow direction of horizontal P -wave velocity. From the seismic results in this study a petrological model suggesting a change of modal composition and percentage of oriented olivine with depth was derived.     

Evidence of spatial variations of attenuation in the western Pyrenean range

Summary. Spectral attenuation of coda waves has been studied in the range 2–40 Hz from local events recorded in the western Pyrenean range from 1980 to 1982. Q _c was obtained using a single scattering model of S -waves for different segments of the coda. An increase of Q _c with lapse time was found and attributed to a rapid increase of Q _β with depth.
Three groups of events were selected from distinct focal areas. Two data sets are mainly composed of aftershocks of moderate earthquakes of magnitude 5.1 and 4.8, respectively. No moderate earthquake occurred in the third area in the few years preceding or following the selected events. Use of stations close to epicentres allowed sampling of the coda at very short lapse times and then study of small, distinct scattering volumes. Noticeable differences were found between the three studied areas and attributed to spatial rather than temporal variations.
The Q _c frequency dependence was studied according to Q _c= qf ^α. α is found to range from 0.7 to 1.1 and q from 30 to 140. These values are in agreement with those found in other tectonic areas. It is shown that scattering is the dominant attenuation process below 10Hz.     

Source history of the 1905 great Mongolian earthquakes (Tsetserleg, Bolnay)

Two great Mongolian earthquakes, Tsetserleg and Bolnay, occurred on 1905 July 9 and 23. We determined the source history of these events using body waveform inversion. The Tsetserleg rupture (azimuth N60°) correspond to a N60° oriented branch of the long EW oriented Bolnay fault.
Historical seismograms recorded by Wiechert instruments are digitized and corrected for the geometrical deformation due to the recording system. We use predictive filters to recover the signals lost at the minute marks.
The total rupture length for the Tsetserleg earthquake may reach up to 190 km, in order to explain the width of the recorded body waves. This implies adding 60 km to the previously mapped fault. The rupture propagation is mainly eastward. It starts at the southwest of the central subsegment, showing a left lateral strike-slip with a reverse component. The total duration of the modelled source function is 65 s. The seismic moment deduced from the inversion is 10²¹ N m, giving a magnitude   M _w = 8  .
The nucleation of the Bolnay earthquake was at the intersection between the main fault (375 km left lateral strike-slip) and the Teregtiin fault (N160°, 80 km long right lateral strike-slip with a vertical component near the main fault). The rupture was bilateral along the main fault: 100 km to the west and 275 km to east. It also propagated 80 km to the southeast along the Teregtiin fault. The source duration was 115 s. The moment magnitude Mw varies between 8.3 and 8.5.
The nucleation and rupture depths remain uncertain. We tested three cases: (1) nucleation and rupture depth limited to the seismogenic zone; (2) nucleation in the seismogenic zone and rupture propagation going to the base of the crust and (3) nucleation within the crust–upper mantle interface and rupture propagation within the upper mantle.     

Crust and upper mantle structure of the central Iberian Meseta (Spain)

Summary. Quarry blasts recorded along three lines on the central Iberian Meseta are used in an attempt to interpret the crustal structure. The results of the interpretation of the data, together with published surface wave and earthquake data, suggest a layered structure of the crust having the following features: the basement, in some areas covered by up to 4 km of sediments, has a P -velocity of 6.1 km s⁻¹; a low-velocity layer, between 7 and 11 km depth, seems to exist on the basis of both P and S interpretation of seismic data; a thick middle crust of 12 km has a P -velocity of 6.4 km s⁻¹ and overlies a lower crust with a mean P -velocity of 6.9 km s⁻¹ and a possible slight negative gradient; the mean v _p/ v _s ratio for the crust is about 1.75; the Moho is reached at about 31 km depth and consists of a transition zone at least 1.5 km thick. The P -velocity of the upper mantle is close to 8.1 km s⁻¹ and the S -velocity about 4.5 km s⁻¹, which gives a v _p /v _s ratio of 1.8 for the uppermost mantle. A tentative petrological interpretation of the velocities and composition of the layers is given.     

Ray-theoretical relations between intrinsic and apparent dissipation factors in the Earth

Summary. A ray-theoretical approach is successfully employed for obtaining the relations between the intrinsic dissipation factors of P - and S -waves, Q _α and Q _β, and the apparent dissipation factors of toroidal and spheroidal oscillations, Q _Tor and Q _SphThese are given in quite simple forms and are expressed, to zero-order approximatiom which neglects the effect of reflection of waves at an interface and/or at a free surface), as where α and β are the P - and S -wave velocities and the integrals are evaluated long the surface to surface P - and S -rays.
Some numerical comparisons support the validity of these formulae for higher modes of toroidal and radial oscillations of a realistic earth model.
A first-order approximation reveals that the existence of a discontinuity in an earth induces a systematic fluctuation in normal-mode Q values for any given phase velocity as a function of frequency (or radial mode number).     

Analysis of the Spatial Coherence of Short-Period Acoustic-Gravity Waves in the Atmosphere

Summary The coherence of atmospheric acoustic-gravity waves has been measured in the period range 10–100 s at the Large Aperture Microbarograph Array in south-eastern Montana. The acoustic-gravity waves observed were signals generated by presumed nuclear explosions. The decrease of coherence with increasing distance between pairs of microbarographs is less rapid in the direction of wave propagation than transverse to it. Variation of direction of arrival over a small range of azimuth (±5°) explains the spatial behaviour of coherence in the direction normal to the wave propagation; variation of phase velocity of ±10 ms^-1 explains the behaviour along the direction of wave propagation. Both effects may be due to inhomogeneities in the atmosphere; the velocity variation may be due to the presence in the signal of several normal modes of acoustic- gravity waves, each travelling at a slightly different phase velocity in the range 300–330 ms^-1.     

Explosive volcanic eruptions — V. Observations of plume dynamics during the 1979 Soufrière eruption, St Vincent

Summary During the 1979 eruption of the Soufrière of St Vincent, West Indies, a 14 min period of explosive activity occurred on April 22. The first three minutes of this activity were filmed. Five individual explosions formed distinct plumes which fed an eruption column which eventually ascended to a height of over 18 km. The ascent velocities of the fronts of these plumes ranged from 8.5 to 61.7 ms^-1. The half-angle of spreading of the plume fronts ranged from 21.5 to 24°. One of the plumes was observed to 8 km height and decelerated steadily from 51.5 to 23 ms^-1. The main column fed by individual explosions was, however, observed to accelerate from 8.5 to 58 m s^-1 in the same height interval.
A theoretical treatment of volcanic plume motion is presented. Measurements of dimensions, velocities and spreading rates from the film are used to estimate plume parameters such as temperature, particle content and volume discharge rate of magma from the theoretical relationships. These calculations show that the individual plumes became hotter and richer in juvenile ash with time. The acceleration of the main eruption column was the result of being fed by increasingly hotter and more ash-rich explosions. An average volume discharge rate of 12600m³ s^-1 is estimated from the analysis of the plume motions. This value agrees closely with an estimate of discharge rate based on the heat flux required to form an 18 km high column. This agreement suggests that theories of convective motions in plumes can be successfully adapted to the volcanic case, as long as the effect of hot entrained ash particles is considered.     

Upper mantle shear structure of North America

Summary. The waveforms and travel times of S and SS phases in the range 10°–60° have been used to derive upper mantle shear velocity structures for two distinct tectonic provinces in North America. Data from earthquakes on the East Pacific Rise recorded at stations in western North America were used to derive a tectonic upper mantle model. Events on the north-west coast of North America and earthquakes off the coast of Greenland provided the data to investigate the upper mantle under the Canadian shield. All branches from the triplications due to velocity jumps near 400 and 660 km were observed in both areas. Using synthetic seismograms to model these observations placed tight constraints on heterogeneity in the upper mantle and on the details of its structure. SS–S travel-time differences of 30 s along with consistent differences in waveforms between the two data sets require substantial heterogeneity to at least 350 km depth. Velocities in the upper 170 km of the shield are about 10 per cent higher than in the tectonic area. At 250 km depth the shield velocities are still greater by about 4.5 per cent and they gradually merge near 400 km. Below 400 km no evidence for heterogeneity was found. The two models both have first-order discontinuities of 4.5 per cent at 405 km and 7.5 per cent at 695 km. Both models also have lids with lower velocities beneath. In the western model the lid is very thin and of relatively low velocity. In the shield the lid is 170 km thick with very high elocity (4.78 km s^-1); below it the velocity decreases to about 4.65 km s^-1. Aside from these features the models are relatively smooth, the major difference between them being a larger gradient in the tectonic region from 200 to 400 km.     

Calibration of the Yellowknife Seismic Array with First Zone Explosions

Summary Recordings from a crustal seismic experiment, which was conducted in the Yellowknife area in 1966, were used for calibration of the Yellow-knife seismic array. In the immediate vicinity of the array the crust is found to be very uniform. A superficial layer with an intercept time of 0–172 ± 0–012s and unknown velocity is underlain by a crust with a P wave velocity of 6.04 ± 0–01 km s^-1 near the top: assuming this velocity constant throughout the second layer, the total thickness of the crust is about 34 ± 2 km. The Mohorovicic discontinuity is horizontal under the array within the resolution of this experiment and the apparent P_n velocity is 8.15 km s^-1. At a distance of a few tens of kilometres the crustal uniformity breaks down. The distances are such that, for most teleseismic signals, the effect of these in homogeneities should be negligible.     

Crustal structure of the Mid-Atlantic ridge crest at 37° N

Summary. A structural model of the Mid-Atlantic Ridge at 37° N is proposed on the basis of travel-time data and synthetic seismograms. At the ridge axis the crust is only 3 km thick and overlies material with an anomalously low'upper mantle'velocity of 7.2 km s⁻¹. Crustal thickening and the formation of layer 3 and a layer with velocity 7.2–7.3 km s⁻¹ takes place within a few kilometres of the axis, producing a 6–7 km thick crust by less than 10 km from the axis. A normal upper mantle velocity of 8.1 km s⁻¹ exists within 10 km of the axis. Shear waves propagate across the axis, thus precluding the existence of any sizeable magma chamber at shallow depth.