Dec 10, 1994 - of travel time information alone leads to a wide class of models compatible with the ... jor discontinuites near 400 and 660 km depth. They also found ..... rays like the secret acoustic submarine monitoring net- works (known by ...
JOURNAL OF GEOPHYSICAL RESEARCH, VOL. 99, NO. B12, PAGES 23,753-23,766,DECEMBER 10, 1994
Seismic heterogeneity in the upper mantle G. Nolet Department of Geologicaland GeopkysicalSciences,Prince.*,on University,Princeton, New Jersey
1919-1994
S. P. Grand Department of Geological Sciences,University of Texas at Ar,stin B. L. N. Kenneff Research School of Earth Science, Australian National University, Canberra
Abstract. This paper gives an overview of recent seismologicalevidence on heterogeneityin the upper mantle. While there is a growing consensusamong seismologists on the low-order sphericalharmonic expansionof upper mantle structure, short- and intermediate-wavelengthfeaturesremain the subjectof intense debate. Regional studiesof Earth structure yield evidencethat the upper mantle is heterogeousat all resolvablescalesand that gl•)l•al models underestimatethe amplitudes of velocity anomalies. Some importanz problems still outstanding
concern:(1) the fate of subductingslabs,(2) the depthextent of low-velocity zonesunderspreadingcentersand (3) the depthat whichhotspotsoriginate.The grossradial layering of the upper mantle has beccnmmuch more evident to us in recent years, mainly as a result of the analysisof waveformsobservedby the modern generationoœdigital, broadbandstationsand by denselyspaced,narrow band arrays of seisinomelets.Nevertheless,important questionson the averageradial structure
still persist.The low shearvelocityzonehasbeen'mappedin detailbeneaththe continents,and Jordan's tectospherehypothesisis largely confirmedby regional
one-dimensional (l-D) and globaltomographic studies.However,it is lessclearto what extent these low shear velocitiesare matched by low compressionalvelocities. Anisotropy,for which compellh•gevidenceexistsin the observationsof split shear waves, may play a role in shaping our images of low velocity zones. Below the asthenosphere,the existenceoœthe Lehmann discc.ntinuityat about 220 km depth as a worldwide feature is unlikely, but early indications for a phase transition at 520 km have found some confixmationin stackedwave sections. Anomalouslylow velocitiesat greater depth have been observednear current,subductionzonesand even beneath more ancient suture zones. Recent progressin analysis of seismic
waveshasfollowedmajor advances in seismicinstrumentation of two kinds- (1) the installationoœvery broadbandstationsaroundthe world,and (2) the increased density oœstations in seismicnetworks. Future progresswill most likely depend on the expansionoœdensenetworks,including temporary and portable stations;on our ability to increasethe station coverageoœoceanicregions;and on advancingthe forward problemof computingwaveformsfor heterogeneous E•rth modelsto a level of efficiencythat allowsfor nonlinear optimization. Introduction The existence of lateral
notably the precisionof timing and the calibration of the instrument response. Early surface wave observavariations
in Earth's
seismic
tions(seeKnopoff[19721 andKovach[1979]forreviews)
velocitieswas knownevento seismologists suchas Jef- establishedthe presenceof a low shear velocity zone, or freys in the 1930s. I-Iowever•the study of lateral vari- LVZ, and showedconsiderablelateral variationsin the ations could not begin in earnest until the quality of thicknessof the high-velocitylid above the LVZ and in seismicobservationsimproved sufficientlyin the 1960s the minimum shear velocity found within it. Attempts were made to regionalize Earth into a small number Copyright1994 by the AmericanGeophysicalUnion.
of geophysical provinces(ancientshields,Phanerozoic
Papernumber94JB01892.
nentsand a subdivision accordingto agefor the oceans),
0148-0227/94/94JB-01892505.00
but efforts to define a common depth-dependentveloc-
platforms, and tectonically active areas for the conti-
23,753
23,754
NOLET
ET AL.: UPPER
MANTLE
HETEROGENEITY
ity profile for equivalent provincesmet with little suc-
this way needsto have a high degreeof continuityover the full distancerange. The structure in the uppermost The excellent seismic data that became available with mantle is evident in observationsout to 1000 km, but the installation of the World Wide Standard Seismoresolutionof strugturein the transition zonedown to graph Network (WWSSN) and the growingpowerof 660 km typically requiresobservationsto 3000 km and computerstriggered the developmentof powerfulalgo- beyond. rithms for the computation of synthetic seismograms. Refractedwavesare sensitiveto velocitygradientsas Detailed velocity models became available for many re- well as to velocity jumps at discontinuities,so there can gionsin Earth through the careful matchingof observed be a significanttrade-off between different aspectsof and computedseismograms,generallyby trial and er- the model during the interpretation process. The use ror. On the other hand, the growing volume of simple of travel time information alone leads to a wide class travel time observationsfor seismic body waves from of modelscompatible with the observations,which can local and global networks inspired the new technique be refined by the application of synthetic seismogram of seismictomography, pioneeredin the 1970sby Aki methods to match the amplitude behavior. Becauseof and Lee [1976],Aki et al. [1977],and Dziewonskiet the long profile range requiredto examinevelocitywith depth, interpretations have commonly been based on Small-scalelength heterogeneitiesin the near-surface single-endedprofiles. Where constraintsare available have always received attention from seismologists, from multiple sourcesalonga singleprofilethere is clear though often as a nuisancewhose effectsmust be re- evidenceof lateral variationsin structure in the upper moved. The useof array observations,particularly with mantle [Mechieet al., 1993]. Other evidencefor the the larger LASA and NORSAR arrays, led to a recog- presenceof significanthorizontalvariability in structure nition that small-scale velocity perturbations are not comesfrom the correlation properties of arrivals from confinedto the upper crust but extend at least into the singlesourcesrecordedat multiple receivers[Kennett lower lithosphere[e.g., Kennett, 1987; W• and and Bodyman,1990]. 1990, Table 1]. The resultsof seismictomographynow There have been two main styles of approach to the indicate that velocity anomalies exist throughout the use of refracted wave information. The first has been to upper mantle, and these observationscontribute to our useshort-periodwaveswhichofferthe possibilityof high understandingof chemicaland thermal heterogeneity. resolution but are also influenced by lateral variations In keepingwith the spirit of other papersin this J(7R in structure. Generally, such studies have been based 75th anniversary series,we shall review the evidencefor on a high density of recordinglocations, either by a lateral heterogeneityin the upper mantle at a tectonic large array for a single sourceor, more commonly,the scale, discussits geodynamicsignificance,and attempt synthesisof a large aperture array by the superposition
cess[Okal,1977;Jordan,1981].
to look ahead
into the near future.
of observations from manysources.In somecases[e.g.,
As is common for a sciencethat is in rapid transition, Bowmanand Kennett,1990;De•l et al., 1993],stacking there is disagreementon the interpretation of various of different observationshas been used to emphasise seismologicalobservations,even among the authors of thosefeatures of the wavefield commonto many events. this paper. Rather than hide such problems,we have A representativeselectionof upper mantle P wavevechosento give more emphasisto unansweredquestions locity modelsderivedfrom short-periodobservationson than to acceptedsolutions. Unlike many other scien- different continentsis shownin Figure 1. Construction
tists, seismologists have little control over their experiments. Consequently,a mixture of techniquesis often used on a patchwork of data from various instruments. Many methods attempt to interpret the data in terms
of each of these models has benefited
from the assem-
bly of a dense record section of observationsto constrain travel time branches,and in most cases,theoretical seismogramcomparisonshave been usedto conof a one-dimensional (l-D, i.e. depth-dependent) model strain the amplitude behaviorwith epicentraldistance. which is valid for a particular region, rather than in The complexfine structure in somemodelsmay result terms of a three-dimensional structure. While the latter from the mapping of three-dimensionalstructure into approachis more logical,the first one givesresultsthat
are often morerobustin the usualcasethat not enough data are availableto allowfor an unambiguous solution in three dimensions.We shall first discusssomeimportant differencesbetween the two philosophies,before directing our attention more specificallyto the results obtained.
Inference
From
1-D
Studies
Most studies of upper mantle structure have been based on the use of refracted waves from natural or man-made sources with receivers at substantial distance
from t'he source. Any 1-D structureto be resolvedin
a one-dimensionai velocitydepth profile [Kenneffand Bowman,1990].The modelsspana widerangeof geographicenvironmentsand interpretation techniquesbut consistentlyshowthe presenceof the two major discontinuities near 410 and 660 km which mark the upper mantle
transition
zone.
The second approach has been to use long-period body waves which are not severely affected by smallscalelateral heterogeneitybut which have limited vertical resolution; in such work, detailed analysis of a limited number of seismogramshas been the usualprocedure. This approach advancedstrongly during the 1970s with the developmentof efficient algorithmsto computesyntheticseismograms.Rapid changesin ave-
NOLET
ET AL.'
UPPER
MANTLE
HETEROGENEITY
C;
I
I
!:
ua• q•,doCl
23,755
23,756
NOLET
ET AL.' UPPER MANTLE
HETEROGENEITY
locity gradient, discontinuities, and low-velocity zones within the upper mantle will produce multiple arrivals at upper mantle distances. Synthetic seismogramscan accurately predict the relative amplitudes and timing
resulting in vastly more useful images of Earth's interior. These early tomographic methods used damped least squarestechniquesto solve the resulting system of linear equations, but computer speed and, most imof sucharrivalswhichcan then be directlycompared portantly, the large memory requirementsof orthodox to seismic data and reduce the nonuniquenessof seis- matrix solversposedlimitations to the size of the probmic modeling. Important early studies matched travel lem that could be handled. Iterative algorithms soon times as well as waveforms of P wavesat upper mantle brought the necessaryimprovementsin both computadistances[Helmbergerand Wiggins,1971; Wigginsand tion speed and problem size. The two techniquesused Helmberger,1973]. Upper mantle shearvelocity(SH) most at present were introduced 10 years ago in seismic profileswereobtainedby HelmbergerandEngen[1974]. tomography: the simultaneousiterative reconstruction These studies using perfectly elastic models generally technique(SIRT [Claytonand Comer,1983])and the resulted in simpler, i.e., smoother, models than earlier conjugategradientalgorithmLSQR [Nolet,1983].The work, although they did confirm the existenceof ma- matrix solver usedin tomographic inversionsdoesinflujor discontinuites near 400 and 660 km depth. They encethe result, but sucheffectsare now well understood also found evidence for a small jump in velocity near [Van der Slutsand Van der Vorst,1987;Nolet, 1993a]. 500 km depth and a low-velocity zone near 150 km Numerous tomography studies have been published depth. Since that time, there have been many studies that usehigh-frequencybody waves,rangingfrom global of the upper mantle beneath various different tectonic studiesthat includethe whole mantle [Inoue et al., provinces using synthetic modeling of upper mantle P 1990; Vascoet al., 1993]throughregionalscales(genand S waves including allowance for attenuation. In erally in subductionareas [e.g., Epakman,1991; Zhou North America, examples include those by Grand and and Cla•tton, 1990; Van der Hilst et al., 1991; Zhao
Helmberger[1984a,b],who studieddifferencesbetween andHasegawa,1993]to very local studiessuchasunder et al., 1983; Ellsworthand Ko•tshieldand tectonicstructure; Walck [1984],who stud- hotspots[Tr•tggvason ied upper mantle P velocity beneath the Gulf of Calianagi, 1977; I•ler et al., 1981]. Grand [1987,1994]has fornia;LeFevreandHelmberger[1989],whostudiedthe used a combination of refracted and vertical incidence upper mantle beneath the Canadian shield; and Zhao shearwavesto imagethe mantle. The shallowmantle in andHelmberger[1993],whostudiedthe P structurebe- these studies showsconsiderableheterogeneityon scale neath
the northwest
The structures
Atlantic.
for both
P and $ waves derived
from
lengths of 500 km or so, but the deeper upper mantle appearslargely uniform. For an extensivebibliography
longer-period observations (Figures2 and 3) showvari-
of tomographicstudies,seeI•ter and Hirahara [1993].
ability in the upper 300 km associatedwith the tectonic character of the surface regions. Shallow low-velocity zonesoccur in tectonically active regionsand beneath oceans, whereas high velocities axe sustained to some depth beneath shields. Comparisonof Figures2 and 3 shows that the low-velocity zones of P and S are different in character. S wave LVZs are generally broad
Low-order spherical harmonic inversionscan still be performed with noniterative methods and now include both waveformsand delay times from long-periodbody
wavesand their surfacemultiples[e.g., Woodhouse and Dziewonski, 1984; Woodwardand Masters, 1991; Suet
al., 1994]. The long-wavelength featuresof Earth show remarkable agreementamong these low-order models, (exceptfor shr14, which was actually derivedfrom P although the question has been raised that the limdata) and muchmoreconsistent from profileto profile. ited samplingfrom availablepaths might lead to spatial Continental models, such as k8 for Eurasia or s8 for aliasingeffectswhichbiasthe interpretation[Sniederet the northeasternUnited States, havenegativegradients al., 1991]. from 150 to 200 km depth but no well-developedlowAs a model that is representative for the most revelocity zone for P. From a petrological point of view cent sphericalharmonic inversioneffortswe have chosen this is to be expected if the low velocities are caused SH12/WM13 [Suet al., 1994]for occasionalreference by partial melting or subsolidusweakeningof the shear and comparison. This model expands the lateral hetmodulus. However, surface wave dispersion measure- erogeneity in Earth's whole mantle into sphericalharments, instrumental in mapping out the low-velocity monics to degree 12 for the horizontal coordinatesand zonesfor S waves, are insensitiveto P wave velocities; 13 Chebychev polynomials for the depth dependence. P wave LVZs are based on body wave analysesonly, The total of 2548 model parameters was estimated from and a bias causedby different analysis methods cannot 27,000 long-period seismogramsand 14,000 travel time completely be ruled out at this point. observations. Three-Dimensional
Studies
It has been arguedthat heterogeneitiesin the mantle occurpredominantlyat length scalesgreater than about
6000 km [Su and Dziewonski,199l]. The strongestargumentsfor this point of view are basedon the analyAbout 20 yearsago, Aki and Lee [1976],Aki et sis of long-period waves traveling over large distances, [1977],andDziewonski et al., [1977],usingtheincreased but it should be kept in mind that these waves themcomputing power available to them, made a dramatic selveshave a filtering effect becauseof the large size improvementin the interpretation of seismicdelaytimes of the Fresnelzonesinvolved(see,for example,Figure
NOLET ET AL.- UPPER MANTLE HETEROGENEITY
4
5
6
7
c•
km/s
8
9
10
23,757
11 t7 s25
k8
200
s8
nj?b ipremc
400
200
600
600
800
800
1000
1000 4,
I
I
I
I
I
I
I
5
6
7
8
9
10
11
o•
km/s
12
Figure 2. Upper mantle modelsfor P velocityderivedfrom long-periodor broadbandobservations. Regionalcoverage:s25, North American Shield; k8, Eurasia; s8, northeasternUnited States;njpb, northernAustralia; ipremc,Isotropicversionof PREM modelwith continentalcrust.
13 km/s 3
4
5
6
7
I
200
200
4OO
I
shr14
2
tna
3
nwa
4
sna
5 njpb 6 ipremc
$,..:. ............ •.:>.
•,400
600
600
800
800
1000
1000 3
I
I
I
I
4
5
6
7
13 km/s Figure 3. Upper mantlemodelsfor S velocityderivedfrom long-periodor broadbandobservations. Regionalcoverage:shr14, westernUnited States;tna, tectonicNorth America;sna, shield North America; njpb, northern Australia; ipremc, isotropicversionof PREM model with continental
crust.
23,758
NOLET
0
ET AL ' UPPER
MANTLE
HETEROGENEITY
Mantle S wave at 25 deg
200
4oo
600
800
lO
20
3O
distance (deg)
O.Os
>5.5s
Figure 4. Fresnelzonesfor mantle S wavescan be very large and complicated. This plot shows the travel time delay for a wave that is scattered once but otherwise follows Snell's law. The travel time delay is plotted at the hypothetical scatteringpoint. The darkest regionsin the plot have a delay •1 s. Because of Fermat's principle, low-velocity zones with a delay •1 s can be evadedby wavestraveling in the darkest zones. The fact that the darkest regionis split into two zonesis a consequence of the existenceof velocity gradientsin the mantle and the triplication of the travel
time
curve.
4), and the conclusion that smallerscaleheterogene- anomaliesin SH12/WM13 and those determinedwith ity is unimportant seemsincompatiblewith more local
waveformstudies[e.g.,Rial at al., 1984;Paulssen,1987; ßDey et al., 1993;gielhuisandNolet, 1994]and with the many resultsfrom delay time tomographyusingshortperiodbody waves.However,the spotty coverageof cal studiesstandsin the way of computingglobalpower
1-D methods. Local studies generally seem to produce larger velocity anomalies than are found in the low-order global inversions,indicating that the filtering action of such models cannot be neglected. Figure 5 gives a histogram for the $ velocity variations at var-
ious depth levelsin model SH12/WM13. At a depth of 136 km, the average$ velocityis 4.452 km/s with a ambiguously.We note, however,that Gudmuhdsson et standarddeviationof 0.068 km/s and extremesof 4.305 al. [1990]use global International Seismological Cen- and 4.660; at 310 km, the averageis 4.707 4- 0.049 km/s tre (ISC) delaytimesto arguefor a whiteheterogeneity with extremesof 4.609 and 4.851 km/s. The range of spectrumfor upper mantle heterogeneitiesat least out velocityvaluesdecreases with depth (Figure 6). Thus, to the resolvingpowerof high-frequency P waves(100- for SH12/WM13, the maximum variation in Vs at 136 200 km), implyingthat the upper mantle is heteroge- km is 4-4% over the whole Earth, with about 95% of all neous at all resolvablescales,none of which is clearly valuesin a narrow range of 4-3%, within 4-2% at 310 dominant. km, and within 4-1% at 660 km. Suchvaluesare easily Additional support for the view that the global exceededin regional and local studies, even outside of models underestimate lateral heterogeneityin the up- the major subduction zones. The models presented by per mantle comes from a comparisonbetween velocity Grand and Helmberger[1984a]haveshearvelocitiesof at high angular order and resolving the question un-
NOLET ET AL.' UPPER MANTLE
Vs (km/s) at 80 km '
I
'
I
'
I
'
I
'
136 km '
I
,
4.5
5.0
5.5
6.0
4.5
5.0
'
,
4.5
433 km
HETEROGENEITY
310 km I
'
I
,
5.0
I
'
I
,
5.5
I
'
I
,
6.0
4.5
533 km
5.5
6.0
4.5
23,759
5.0
5.0
5.5
6.0
5.5
6.0
703 km
5.5
6.0
4.5
5.0
Figure 5, Histogramsfor S wavevelocityvariationsin the globalmodelSH12/WM13 [S• et 1994]at variousdepth levels.
Model SH 1 2
100
200
300
400
500
600
700
Depth (km) Figure 6. A sequenceof 1-D modelsgeneratedby e•aiuatingthe depth dependence of model SH12/WM13 at 411pointson Earth, spacedabout 1100km apart. The velocitydiscontinuities are imposedby the background model(PREM), the depthdependence of t;helaterallyheterogeneity is smooth with a resolvingpower in the parametrization of about 230 km.
23,760
NOLET ET AL.' UPPER MANTLE HETEROGENEITY
4.29 km/s and 4.775 km/s at 125 km depthbeneath ters[1991].Suchrootsforman essentialelementin Jorthe eastern Pacific and Canadian shields respectively. dan's[1975]tectosphere hypothesis.Structuraldetails A modelderivedby Humphre,Iset al. [1984]for south- within the tectosphereare likely. Grand and Helmberger ern Californiashowsa thin (about 50 km wide) tabu- [1984a]founda negativevelocitygradientnear 175km lar structurewith a 2-3% fast anomalyextendingfrom
depth beneaththe Canadian shield. œeFe•reand Helm-
near the surface to 250 km depth beneath the Trans-
berger[1989]founda similarstructurefor P waves,also
verse Ranges. A 2-4% slow anomaly is also present in this area to 125 km depth giving a total range in mantle heterogeneityof 4-7% within the upper 150 km of southernCalifornia. Humphreysa•d Duecker[1994] have extended their studies to arrays throughout the western United States, finding peak to peak variations of mantle velocity as large as 8%. Over the whole of
beneath
the western United States there appears to be a periodicity in the structural variations with a wavelength of about 500 km. Sipkin and Jordan [1980] have found a similar scale length variation in the western Pacific
using multiple bounce $c$ waves. It is also interesting to note that for the western U.S. arrays, Humphre,Is and Duecker[1994]only reducethe RMS of the travel time residualsby about 50%, even though they allow variations
to occur over scales as short as 20 km.
The
unexplained delays appear coherent and are unlikely
to be due to randomerrorsin time measurements (E. D. Humphreys,personalcommunication, 1993),and it may be that there are large-amplitudeheterogeneities with very short wavelengthsthroughoutthe upper few hundred
kilometers
of the mantle.
Evidence for pervasive smaller-,•;cale heterogeneityis alsofound in observationsof short-period arrivals across
the Canadian
Low-Velocity
shield.
Zones
The existence of a low-velocity region beneath the Phanerozoic platforms has been establishedin many early surface wave studies. A preciseinterpretation of the dispersion data is hampered by,recognition that the upper mantle under many continental regions ex-
hibits seismicanisotropy[Babuskaet al., 1984; Regan and Anderson,1984] which, by itself, also may show largelateral changes,as observedfrom shearwavesplitting measurements [Silver and Chan, 1988]. For example, recent broadband observationsbeneath northern Australia [Kennett et al., 1994; 1994]haverevealeda verylargecontrastbetweena 210km-thick lithospherewith very low intrinsicattenuation and the astenosphere beneathwith a very high lossand slightlyreducedshearvelocities. Shearwavesplitting observations for the refractedS waves[Tonget al., 1994] are compatiblewith about 1% higher horizontalshear wavespeedsin the asthenospheric zone,whichmay well representthe locationof greater mantle flow. SH12/WM13 showsa pattern of low shearwavevelocitiesextendingto at least 250 km beneath ocean
portable arrays with a range of aperturesand station ridges,confirming earlysurfacewavestudies[Wielandt
spacings[Ke•mettand Bodyman, •990]. To matchthe
andKnopoff,1982]. 2hanga•d Tanimoto[1992]have
relatively rapid changesin amplitude behaviorfor refracted wavesfrom singleeventsacrosssucharrays,heterogeneityon scalelengthsof 200-300km with at least 1% variability is required. Kenneff a•d Nolet [1990] have demonstrated that such small-wavelengthheterogeneitydoesnot have a significantinfluenceon longerperiod surfacewavesand so would not be detectedin typical waveform tomography.
recentlychallengedthe view that the low velocitiesunder the ridgesextendto a depthwellin excessof 100km. However,the poor fit of the predictionsfor the Zhang
Major Upper Mantle Features In this section we review various results on the seismic
expressionof tectonicfeatures,and will focusattention to someimportant problemsstill outstanding. Deep Continental
Roots
and Tanimoto
model to observed $$ wave travel times
putstheir conclusion in doubt[$uet al., 1992].This debate over the depth extent of mid-oceanicridgesseems to mirror that over the power in small-scale heterogeneities:in both cases,data setswith differentresolution lead to incompatibleresults. The imagingof smallscaledetail of upwellingof magmabeneathoceanridges is made very difficult by the combinationof very few seismicstations located in oceanic areas, the poor reso-
lution at depth usingfundamentalmode surfacewaves, and the fact that imagingof low velocitiesis hampered by diffractionof wavesaroundthe low-velocityobstacles [Wielandt,1987]. The same is even more true for hotspots. Imaging ex-
ModelSH12/WM13 [$uet al., 1994]has400-km-deep perimentsbeneath the Hawaii hogspot[Ellsworthand keelsbeneath Canada, Baltica, much of Eurasia, Africa, Ko•anagi, 1977],Yellowstone [I•er et al., 1981],and et al., 1983]all indicatea sublithoBrazil, and the western part of Australia. The vertical Iceland[Tr•lggvason sphericorigin for the mantle plume, but resolutionis quicklylost at depth, denyingus an answerto the important questionof upper or lower mantle origin for hotspots.Surprisingly, Nataf and VanDecar[1993]reand Helmberger[1984a,b]and œerner-œam and Jordan port a very small but clear delay time signalthat may [1987],who showhigh 5H velocitiesextendingto such be causedby a plume beneaththe Bowiehotspot,indidepth, and fits the correlation of large SS-S residuals catinga lowermantle originfor the upwellingmaterial. resolution is that of the radial parametrization used, or about 230 km, making a depth extent for roots beneath shield areas to at least 300 km very probable. This is in agreement with regional models of Grand
with tectonic province shown by Wood•vardand Mas-
This observation could be made thanks to the fortunate
NOLET
ET AL.: UPPER
MANTLE
location of a denseseismicarray, the WashingtonRegional SeismicNetwork, with digital recording,making very accurate delay time readingspossible. Low velocities, together with low Q, have also been observedin or just above the transition zone in active subductionzones.Sacksand Okada[1974]report very low Qp valuesabovethe deepparts of the slabsin Japan and Fiji and similarly low Qp under Peru. Zhou and
it, t
orro-
orated these observationsby imaging low velocitiesin the west Pacificcorrespondingto the low Qj> zonesand alsofoundvelocityanomaliesoftenextendingawayfrom the slab. This is alsoobservedby Revenauõhand Sipkin
HETEROGENEITY
23,761
ily visible in their data and do showvariationsin depth with peak to peak amplitude variationsof about 30 km,
possiblyanticorrelated[Revenaugh and Jordan,1991a] with the largest uncertainty in the depth contributed
by the upper mantle velocity [Stammleret al., 1992]. A negative correlation would be in agreement with a thermal origin for the variations since the sign of the Clapeyron slope is positive for the olivine to •3 spinel transition at 400 km but almost certainly negative at
660 km [Ito and Takahashi,1989]. Evidencefor longwavelength topography was also found from underside
reflections PiooP [NeelearidSnieder,1991]. Recent modelsshowa relatively sharp changein ve-
[1994],whoreporta zoneof low SH velocityjust above locity at the major discontinuities,althoughthe details the 410-kin discontinuityextendingfrom the Japanese
of the
subduction
lengths involved. The model shr14 for shear velocity
zones beneath
China.
Zielhuis
and Nolet
transition
cannot
be resolved
with
the wave-
[1994]observea low S wavevelocityanomalyat tran-
[Helmbergerand Enger•,1974]did not includeattenu-
sition zone depths beneath an ancient suture zone in
ation, but an allowance for attenuation was included
easternEurope,while Taylor and Toksoz[1979]found in later studies[e.g., Grand and Helmberger,1984b]. low P wavevelocitiesat the locationof the subducting Kennett [1975]pointed out that when attenuationwas proto-Atlantic in the northeastern United States. No- includedin upper mantle models,the pattern of amplilet and Zielhuis[1994]conjecturethat suchzonescan tudes associatedwith the presenceof the major velocbe explained by a mechanismof water injection which ity transitionscould be matchedwith a rapid change would permanentlyerodethe tectosphereat its edges. compatible with the constraints on transition thickness inferred by Richards[1972]from observationsof preUpper Mantle Discontinuities cursorsto P'P' (PKPPKP). A high gradient over a A commonfeature of both short- and long-periodre- depth interval of 10 km will produceessentiallythe same fraction studiesis a relatively high gradient within the synthetic seismogramat wide anglesof incidenceas a transitionzone(from410 km to 660 km depth). The sharpdiscontinuityevenat frequenciesas high as I Hr. overall increasein velocity from about 425 to 625 km is [Burdickand Orcutt, 1979]. At longerperiods,such well resolved;however,the detailsof the gradientwithin as those used to study shear structure, the situation is this depth interval are uncertain, and the existenceof a worse,and a sharpjump in velocityor a high gradient over a depth interval of about 40 km cannot be distinguished. The strongestevidencefor the sharpnessof The magnitude of the jump in P wave velocity across discontinuities comes therefore from vertically incident the 410-km discontinuityis highly variableamongdif- waves. Adams[1968, 1971],Engdahland Flinn [1969], ferent models;it variesfrom as little as 0.23 km/s for and Whircomband Anderson[1970]were the first to PREM [Dziewonski andAnderson,1981]to as muchas observeprecursorsto PKPPKP waves at frequenciesof 1.20 km/s in a modelfor the Fiji plateau [Frohlichet about 1 Hr. which are interpreted as near-vertical incial., 1977]. Sucha large disparitymost probablyarises dence reflections from discontinuities within the upper from differencesin the style of interpretation in combi- mantle. Nakanishi [1989]and Benz and }ridale[1993] nation with the limited resolvingpower of the available have performed more recent studies on the sharpness data, rather than from significanttectonicsignals.The of the transition. In general, the 660-kin discontinuvariability for the jump in $ wave velocity at the 410- ity appearsto be a good reflectorof near-vertical, 1-Hr. km discontinuity among different models is much less P waves indicating that the transition is completed in an interval of 4 km or less[Richards,1972]. Support[Nolet and Wortel, 1987]. The data modeled with synthetic seismogramsare ing evidencefor the sharp nature of the discontinuity generally wide-anglereflectionsfrom discontinuitiesor also comesfrom P to S convertedwaves[Vinnik, 1977; turning rays. The travel times of suchphasesare very Stareruleret al., 1991]. Local focusingand defocusing sensitive to velocities in a r.one about 100 km thick of these waves has been used to infer a topography on above the discontinuities,resulting in a trade-off be- the 660-kin discontinuity of more than 10 km amplitude tween velocity and discontinuity depth. The travel over wavelengthsas short as 200 km beneath western times of near-vertical reflections from discontinuities are Europe[Paulssen,1988; Van der Lee et al., 1994]. Such better suited to estimating depths of discontinuities, pronouncedtopography could be expected to disrupt but depth estimates are still susceptibleto the influ- the refracted wave field and so alter the inferred sharpence of the reference velocity used to convert time to ness of the discontinuity. The 400-kin discontinuity is generally a much poorer d•pth. ••• •d •o•d• [•OO•,b], S•,• [•], and ShearerandMasters[1992]haveusedstackedlong- reflector, indicating that it is probably spreadover more period data to searchfor near-verticalupper mantle re- than 5 km depth or that it exhibits topography that flections.The 410- and 660-km depth reflectorsare eas- influences the amplitudes of reflected waves. It should
sharp discontinuityat about 520 km depth is not well established.
23,762
NOLET
ET AL.: UPPER
MANTLE
be noted, however,that in someregionsthe 400-kin dis-
continuitydoesreflectshort-periodwaves[Bowmanand Kennett, 1990; Vidale and Benz, 1992]. The variability of observationsof reflected or conv,;rtedphasesfrom the major discontinuitiesis a longstandingproblem in seis-
mology. The variations may reflect lateral changesin the character of the discontinuity, but a definite answer cannot yet be given. Someregionalmodels[e.g.,Mechieet al., 1993]have a high gradient zone or a discontinuity near 520 km depth with lower gradientsabove and below this depth. Shearer[1990]recentlyrevived the discussion on the existenceof a phasetranformation in the transition zone by his observationof stackedphasesoriginatingat this depth,whileRevenaugh andJordan[1991a]alsoobserve
HETEROGENEITY
is still considerabledebate on the fate of subducting slabsonce they are old and denseenoughto sink to the bottom of the upper mantle. In an effort to avoid the problemsposedby poor res-
olution, Creagerand Jordan [1984]averagethe effects of wave propagationoutside the mantle by projecting all delay times on a residualsphere. From their results and forward modeling,they concludethat slabpenetration into the lower mantle occurs, implying a convec-
tive regimethat involvesthe wholemantle rather than a stratified system. For portions of some subduction zonese.g., the northern Kurilesand the Marianas,penetration of slabs is confirmed by seismictomography. In other areas, such as Izu Bonin and central Japan, the slab appearsto deflect[Van der Hilst et al., 1991] and may remain stagnant on top of the 660-kin discona minor reflector at a mean depth of 520 km. Other modelshave a more uniform gradient which is greater tinuity [Fukaoet al, 1992]. The nature of the fate of than that predicted for self-compressionof the spinel subducted slabs appears to have a strong correlation phaseof olivine [Bassand Anderson,1984]. However, with the tectonichistoryof a region[Van der Hilst and Joneset al. [1992]and Cummins.et al. [1992]find no $eno, 1993]. compellingevidencefor a jump in velocitynear 520 km from analysisof waveformsreturned from the transition What Next? zone.
Secondarydiscontinuitiesare no• as well constrained.
A numberof profilesfrom shieldareas (8-13 in Figure 1) havestructurenear 210 km with a zoneof reduced velocity beneath. These observationsdirectly
contradictan earlierview by Anderson[1979]that the Lehmann discontinuity, with a positive jump in velocity, is a worldwide feature with a common cause and
Much of the progressin unraveling small-scaleseismic structure of the upper mantle was generated by a formidable increasein both the quantity and quality of seismicdata sincethe 1960s. The availability of almost 30 years of arrival times from a worldwide network of seismicstations reported to the International Seismo-
logical Centre made possiblemuch of the imaging of commoncharacteristics.tlevenaugh. andJordan[1991b] subductionzonesthrough delay time tomography.The have reported $H wave discontinuitiesnear 300 km be- installation of the WWSSN network, now largely reneath the shield region of westernAustralia, and near placed with the new digital, broadband stations, trig200 km in northernAustralia whereother studies[e.g., gered the analysisof long-periodwavesso important Leven, 1985; Bowmanand Kennett, 1990] (Figure 1) to the global imaging of structure in the upper mantle
have proposeda velocityjump usingrefracted P waves. and below. Increased access to observations is now beRecentinterpretationsof Russiandeepseismicsounding ing providedthroughthe Internet directly from stations data [Mechieet al., 1993]introducea complexfeature and through internationally operatingdata centerssuch in velocity near 210 km; see profiles 12 and 13 in Fig- as the IRIS Data Management Center in Seattle or the
ure 1. Vidale and Benz [1992]alsoreport observinga
Orfeus
reflector near 200 km depth within subduction zones. These discontinuities,however,have not been proposed as global features, and their interpretation is still an open question. Subduction
Despite all these gains in data a major problem in mapping heterogeneitycontinuesto be the limited resolving power of the available data sets. Local tomographic investigationsare hampered, often more than global studies, by the patchinessof data coverage.The
Subduction zones cause the most pronouncedvelocity anomalies in Earth's upper mantle. Gubbins and
The density of seismic•a,tions at the surfaceof Earth generally falls far short of what would be needed to
Data
Center
in the Nether]ands.
lackof resolutionis causedby a numberof factors:(1)
inferred,(2) $nieder [1991] report arrival times for compressional match the scalelengthsof heterogeneity wavesthat indicate path-averagedslab velocity anoma-
the sources,earthquakesor explosions,are not spread
lies in excessof 5% for the Tonga-Kermadecslab, although the layer with the highest velocity may be thin
uniformly,(3) the precisionof routinedelaytime readingsis oftenlow, and (4) thereare theoreticallimita-
(lessthan 15 km) as evidencedby the dispersive char- tionsto the resolvingpowerevenif we had an ideal coverageof wavepathsthroughoutthe planet. Progressis
acter of the wave. Since the anomaly is positive, the fast rays arrive before any diffracted waves, making slabseasierto image than low-velocityzones. Moreover, there is abundant seismicity at depth in most conver-
gent areas, resulting in increasedresolvingpower, despite the common lack of seismic stations seaward of the zone. Despite these fortuitous circumstances,there
to be expected on all these points. The density of seismicstations can be improvedconsiderablywith moderncommunicationmethods. Cellular phone use is expectedto grow stronglyin the next decade,and the planned use of a multisatellitesystem for commercial
communication
will make the transfer
NOLET ET AL.: UPPER MANTLE
of digital data from a massivenetwork of simple seismometers with accurate time stamps possible. Digital analysisof wavesenablesvery accurate measurementof
delay timesusingcross-correlation techniques[VanDecar a•d Cro•on, 1990]. By improving measurement procedures [here is a good chance that the accuracy of travel time observations can be upgradedby an order of magnitude, with a consequentimprovement in resolvingpower, whereas to achievethe sameimprovementin resolutionthrough statistical averagingof more observationswould require an increase in the quantity of the data by 2 orders of magnitude. At sea,deploymentof large arrays of hydrophonesis a provenmilitary technology,and the use of data from arrays like the secret acousticsubmarinemonitoring net-
HETEROGENEITY
23,763
arrays would actually be quite modest and might give us the seismologicalequivalent of a Hubble telescopefor only a fraction of the cost. The recent rate of improvement in the quality of seismic data has outstripped the comparable development of interpretation tools. Investments in recording technologywill be of greatestbenefitif it can be matchedby improved methods of modeling the influenceof heterogenity on the seismicwavefield. Nonlinear optimization
hasbeensuccesfully appliedto waveformfitting [Nolet, 1990],extendingits applicabilityto frequencies as high as 50 mHz. For the interpretation of higher frequencies, it will be necessary to incorporate scattered wave energy into the inversion. Current techniquesfor the coupling of normal modes are too inefficient to be applied directly in nonlinear optimization schemes. Tarantola
works(knownby acronymssuchas SOSUS)in scientific [1987]and Geller andHara [1993]haverecentlydevelinvestigationsawaits a relaxing of security requirments, and a more agressivepolicy of data preservation until then. Since such arrays are commonly located at the seaward side of oceanic trenches,a large increasein resolving power would become available that might help decide the fate of subductingslabs. We have to wait passively for earthquakes to occur, and just now seismologyhas gone through a few
oped efficient techniquesto invert waveform data which have already found application in low-orderglobal stud-
very quiet decades.For example,the deep(600 km)
lies of surface waves contain information that is only
earthquake that occurredon June 9, 1994, in northern Bolivia was the first large deep earthquake since the Colombian event of July, 31, 1970. Since then the very broadband, digital networks have been installed worldwide, and it is expected that these will yield superior observationsof normal mode overtones. Among other results, such data could greatly improve our knowledge of lateral variations in the attenuation factor Q. Global variations in this important seismicparameter
ies [Hara et al., 1993]. Alternatively, Sniecler[1986] has shown how to use linear inversion techniques for scattered fundamental mode Rayleigh waves in a firstorder Born approximation, and this method is potentially very powerful when large, dense seismicarrays become
available.
Polarization
and azimuthal
anoma-
beginningto be modeled[Lerner-Lainand Park, 1989; Le•shin et al., 1994; Laske and Masters, 1993]. Attempts to regionali•.eseismicanisotropyare crude and
likely to improvein the nearfuture [e.g.,Nishimuraand Fors•lth,1989;Park and Yu, 1992].
are now only knownto the lowestangularorder [e.g., Romaaowicz,1990], and more local studiesstill suffer
In order to improve the quantification of small-scale heterogeneity in the upper mantle it will be necessary to be able to model the propagation of seismicwaves at higher frequencies. This will need the extension of interpretation methods to include higher order scattering, involving mode-mode conversions,for both body
from a variety of instrumental and interpretation prob-
and surface waves.
lems [Nakanishi,1993]. In the absenceof earthquakesourcesin many regions of Earth, there is a growing demand for safe but powerful artificial sources. The technologyfor large vibrators has recievedlittle attention so far, but efforts by Russian seismologistshave been successful,and such instruments could be deployed in inactive areas. Finally, the limits on resolvingpower in current theory can be improved considerablyif we find ways to efficiently incorporatescatteredenergyin the interpretation. This does not necessarilyinvolve only theoretical progress. The example of the oil exploration industry showsus that the availability of a massiveoverkill in in-
strumentationallowsusto applysimplified(firstorder)
Acknowledgments. G.N. acknowledges partial support from NSF grant EAR9204386. References Adams, R. D., Early reflections of P'P' as an indication of upper mantle structure, Bull. Seisrnol. Soc. Am., 58, 19331947, 1968. Adams, R. D., Reflections from discontinuities beneath Antarctica, Bull. Seisrnol. Soc. Am., 61, 1441-1451, 1971. Aki, K., and W.H.K. Lee, Determination of three-dimensional velocity anomalies under a seismic array using first P arrival time from local earthquakes, 1, A homogeneous initial mode, J. Geophys. Res., 81, 4381-4399, 1976. Aki, K., A. Christoffersen, and E.S. Husebye, Determination
scattering theory on filtered data. of the three-dimensional seismic structure of the lithoWe concludethat most of the progresson decyphersphere, J. Geophys. Res., 8•, 277-296, 1977. ing upper mantle structure is to be gained by increasing the density of seismicstations, either permanently Anderson, D.L., The deep structure of continents, J. Geophys. Res., 8•, 7555-7560, 1979. or temporarily deployed. With continued development Babuska, V., J. Plomerova, and J. Sileny, Large-scaleoriof accurate seismicsensors,and the prospect of mass entated structures in the subcrustal lithosphere of central production of digital seismometerswith a reasonably Europe, Ann. Geophys., •, 649-662, 1984.
broadbandresponse[Nolet,1993b],the costof massive Bass, J.D., and D.L. Anderson, Composition of the up-
23,764
NOLET ET AL.: UPPER MANTLE HETEROGENEITY
per mantle: Geophysical tests of two petrological models, Geophys. Res. œett., 11, 237-240, 1984. Benz, H. M., and J. E. Vidale, The sharpnessof upper mantle discontinuities determined from high frequency P'P' precursors, Nature, $65, 147-150, 1993. Bowman, J.R., and B.L.N. Kennett, An investigation of the upper mantle beneath NW Australia using a hybrid seismograph array, Geophys. J. Int., 101, 411-424, 1990. Burdick, L.J., and D.V. Helmberger, The upper mantle P velocity structure of the western United States, J. Geophys. Res., 85, 1699-1712, 1978. Burdick, L.J., and J.A. Orcutt, A comparisonof the generalized ray and reflectivity methods of wave form synthesis, Geophys. J. R. Astron. Soc., 58, 261-278, 1979. Clayton, R.W., and R.P. Comer, A tomographic analysis of
mantle heterogeneities from body wavetravel times (abstract), Eos Trans. A GU, 6j, 776, 1983. Creager, K.C., and T.J. Jordan, Slab penetration into the lower mantle, J. Geophys. Res., 89, 3031-3049, 1984. Cummins, P.B.., B.L.N. Kennett, J.R. Bowman, and M.G. Boarock, The 520 km discontinuity?, Bull. Seistool. Soc. Am., 8•, 323-336, 1992. Dey, S.C., B.L.N. Kennett, J.R. Bowman, and A. Goody, Variations in the upper mantle velocity structure under northern Australia, Geophys. J. Int., 11j, 304-310, 1993. Dziewonski, A.M., and D.L. Anderson, Preliminary reference earth model, Phys. Earth Planet. Inter., •5, 297-356, 1981.
Dziewonski, A.M., B.H. Hayer, and R.J. O'Connell, Largescaleheterogeneitiesin the lower-mantle, J. Geophys.Res., 82, 239-255, 1977.
Ellsworth, W.L., and R.Y. Koyanagi, Three-dimensional crust and mantle structure of Kilauea Volcano, Hawaii, J. Geophys. Res., 8•, 5379-5394, 1977. Engdahl, E.R., and E.A. Flinn, Seismicwavesreflectedfrom discontinuities within the Earth's upper mantle, Science, 163, 177-179, 1969. Faber, S., and G. Muller, Sp phases from the transition zone between the upper and lower mantle, Bull. Seistool. Soc. Am., 70, 487-508, 1980.
Frohlich, C., M. Barazangi, and B. L. Isacks, Upper mantle structure beneath the Fiji Plateau: Seismic observations of second P arrivals from the olivine-spinel phase transition zone, Geophys. J. R. Astron. Soc., 50, 185-213, 1977. Fukao, Y., M. Obayashi, H. Inoue, and M. Nenbai, Subduction zones stagnant in the Earth's mantle, J. Geophys. Res., 97, 4809-4822, 1992. Geller, R.J., and T. Hats, Two efficient algorithms for iterative linearized inversion of seismic:waveform data, Geophys. J. Int., 115, 699-710, 1993. Grand, S. P., Tomographic inversion for shear velocity beneath the North American plate, J. Geophys. Res., 9•, 14,065-14,090, 1987.
Grand, S. P., Mantle shear structure beneath the Americas and surrounding oceans, J. Geophys.Res., 99, 11,591-11, 621, 1994.
Grand S., and D.V. Helmberger, Upper mantle shear structure of North America, Geophys. J. R. Astron. Soc., 76, 399-438, 1984a.
Grand, S. P., and D. V. Helmberger, Upper mantle shear structure beneath the northwest Atlantic Ocean, J. Geophys. Res., 89, 11,465-11,475, 1984b. Gubbins, D., and R. Snieder, Dispersion of P waves in subducted lithosphere: Evidence for an eclogitelayer, J. Geophys. Res., 96, 6321-6333, 1991. Gudmundsson, O., J.H. Davies, and R.W. Clayton, Stochastic analysisof global travel time data: Mantle heterogeneity and random errors in the ISC data, Geophys. Y. Int., 10•, 25-44, 1990.
Gudmundsson O., B.L.N. Kennett, and A. Goody, Broadband observationsof upper mantle seismicphasesin northern Australia and the attenuation structure of the mantle, Phys. Earth Planet. Inter., in press, 1994. Hats, T., S. Tsuboi, and R.J. Geller, Inversionfor laterally heterogeneous upper mantle S-velocitystructureusingiterative waveform inversion, Geophys. Y. Int.,115, 667-698, 1993.
Helmberger, D.V., and G.R. Engen, Upper mantle shear structure, Y. Geophys. Res., 79, 4017-4028, 1974.
Helmberger,D.V., and R. A. Wiggins, Upper mantle structure of the midwestern United States, Y. Geophys.Res., 76, 3229-3245, 1971. Humphreys, E. D., and K. G. Duecker, Western U.S. upper mantle structure, Y. Geophys. Res., 99, 9635-9650, 1994. Humphreys, E. D., R. W. Clayton, and B. H. HageL A
tomographicimage of mantle structurebeneathsouthern California, Geophys.Res. Left., 11, 625-627, 1984.
Inoue, H., Y. Fukao, K. Tanabe, and Y. Ogata, Whole mantle P-wavetravel time tomography,Phys. Earth Planet. Inter., 59, 294-328, 1990.
Ito, E., and E. Takahashi, Postspineltransformationsin the system Mg•.SiO4- Fe•.SiO4 and some geophysicalimplications, Y. Geophys.Res., 9j, 10,637-10,646, 1989.
Iyer, H.M., and K. Hirahara (Eds.), SeismicTomography, Theory and Practice,842 pp. Chapman and Hall, London, 1993.
Iyer, H.M., et al., A deep low velocity body under the Yellowstonecaldera, Wyoming: Delineation usingteleseismic P wave residuals and tectonic interpretation, Geol. Soc. Am. Bull., 9•, 792-798, 1981.
Jones, L. E., J. Mori, and D. V. Helmberger, Short-period constraints on the proposed transition zone discontinuity, J. Geophys. Res., 97, 8765-8774, 1992.
Jordan, T.H., The continental tectosphere,Rev. Geophys., 13, 1-12, 1975.
Jordan, T.H., Global tectonic regionalizationfor seismological data analysis, Bull. Seismol. Soc. Am., 71, 1131-1142, 1981.
Kennett, B.L.N., The effectsof attenuation on seismograms, Bull. Seismol. Soc. Am., 65, 1643-1651, 1975.
Kennett, B.L.N., Observational and theoretical constraints on crustal and upper-mantle heterogeneity, Phys. Earth Planet. Inter., j7, 319-332, 1987.
Kennett, B.L.N., and J.R. Bowman, The velocity structure and heterogeneityof the mantle, Phys. Earth Planet. Inter., 59, 429-465, 1990.
Kennett, B.L.N. and E.R. Engdahl, Travel times for global earthquake location and phaseidentification, Geophys.J. Int., 105, 429-465, 1991. Kennett, B.L.N., and G. Nolet, The interaction of the S wavefield with upper mantle heterogeneity, Geophys. J. Int., 101, 751-762, 1990. Kennett, B.L.N, O. Gudmundsson,and C. Tong, The upper mantle S and P velocity structure beneath northern Australia from broad-band observations Phys. Earth Planet. Inter., in press, 1994. Knopoff, L., Observation and inversion of surface wave dispersion, Tectonophysics,13, 497-519, 1972. Kovach, R.L., Seismicsurface waves and crustal and upper mantle structure, Rev. Geophys., 16, 1-13, 1979. Laske, G., and G. Masters, Polarization measurements of long-periodsurfacewavesand globalphase-velocitymaps,
Eos Trans. AGU, 7j(j3), Fall Meetingsuppl.,438, 1993. LeFevre, L.V., and D.V. Helmberger, Upper mantle P velocity structure of the Canadian shield, J. Geophys.Res., 9q{, 17,749-17,765, 1989.
Lerner-Lam, A.L., and T.H. Jordan, How thick are the continents?, J. Geophys.Res., 9•, 14,007-14,026, 1987.
NOLET
ET AL.: UPPER MANTLE
Lerner-Lam, A.L., and J.J. Park, Frequency dependent refraction and multipathing of 10-100 secondsurface waves in the western Pacific, Geophys. Res. œett., 16, 527-530, 1989.
Leven, J.H., The application of synthetic seismogramsin the interpretation of the upper mantle P-wave velocity structure in northern Australia, Phys. Earth Planet. Inter., 38, 9-27, 1985.
Levshin, A.L., M.H. Ritzwoller, and L.I. Ratnikova, The nature and cause of polarization anomalies of surface waves crossing northern and central Eurasia, Geophys. 3. lnt., ll7,
577-590, 1994.
Mechie, J., A.V. Egorkin, K. Fuchs,T. Ryberg, L. Solodilov, and F. Wenzel, P-wave mantle velocity structure beneath northern Eurasia from long-range recordings along the pro•e Quartz, Phys. Earth Planet. Inter., 79, 269-286, 1993.
Nakanishi, I., Reflections of P'P' from upper mantle discontinuities beneath the Mid-Atlantic Pddge, Geophys. 3. R. Astron. Soc., 93, 335-346, 1988. Nakanishi, I., A search for topography of the mantle discontinuities from precursorsto P'P', 3. Phys. Earth, 37, 297-301, 1989.
Nakanishi, I., Surface wave tomography: Velocity and Q, in Seismic Tomography, Theory and Practice, edited by H.M. Iyer and K. Hirahara, Chapman and Hall, London, pp. 92-132, 1993.
HETEROGENEITY
23,765
Revenaugh, J.S., and T.H. Jordan, Mantle layering from ScS reverberations, 2, The transition zone, 3. Geophys. Res., 96, 19,763-19,780, 1991a.
Revenaugh,J.S., and T.H. Jordan, Mantle layeringfrom ScS reverberations, 3, The upper mantle, 3. Geephys.Res., 96, 19,781-19,811, 1991b.
Revenaugh, J., and S.A. Sipkin, Seismicevidencefor silicate melt atop the 410-km mantle discontinuity, Nature, 369, 474-476, 1994.
Rial, J.A., S. Grand, and D.V. Helmberger, A note on lateral variations in upper mantle shear velocity acrossthe Alpine front, Geophys. 3. R. Astron. Soc., 77, 639-654, 1984. Palchards,P.G., Seismic waves reflected from velocity gradient anomalies with the Earth's upper mantle, Z. Geophys., •8, 517-527, 1972.
Romanowicz, B., The upper mantle degree 2: Constraints and inferences from global mantle wave attenuation measurements, 3. Geophys. Res., 95, 11,051-11,071, 1990. Sacks, I.S., and H. Okada, A comparison of anelasticity structure beneath western South America and Japan, Phys. Earth Placket.Inter., 9, 211-219, 1974. Shearer, P., Seismic imaging of upper-mantle structure with new evidence for a 520 km discontinuity, Nature, 344, 121126, 1990. Shearer, P.M., Constraints on upper mantle discontinuities from observationsof long-period reflected and converted phases, 3. Geophys. Res., 96, 18,147-18,182, 1991.
Natal, H.C., and J. VanDecaL Seismologicaldetection of a Shearer, P.M., and T. G. Masters, Global mapping of topography on the 660-km discontinuity, Nature, 355, 791-796, mantle plume?, Nature, 364, 115-120, 1993. 1992. Neele, F., and R. Snieder, Topography of the 400 km disSilver, P., and W.W. Chan, Implications for continental continuity from observationsof long period P4ooPphases, structure and evolution from seismic anisotropy, Nature, Geophys.3. R. Astron. Soc., •09, 670-682, 1991. 335, 34-39, 1988. Nishimura, C.E., and D.W. Forsyth, The anisotropic strucSipkin, S. A., and T. H. Jordan, Multiple ScS travel times in ture of the upper mantle in the Pacific, Geophys. 3. Int., the western Pacific: Implications for mantle heterogeneity, 96, 203-230, 1989. 3. Geophys. Res., 85, 853-861, 1980. Nolet, G., Inversion and resolution of linear tomographic Snieder, R., 3D linearized scattering of surface waves and a systems(abstract), Eos Trans. AGU, 6J, 775-776,1983. formalism for surfacewave holography, Geophys.3. R. AsNolet, G., Partitioned wave-form inversion and two-dimen-
sional structure under the NARS array, 3. Geophys.Res., 95, 8513-8526, 1990.
Nolet, G., Solving large linearized tomographic problems, in Seismic Tomography, Theory and Practice, edited by H.M. Iyer and K. Hirahara, Chapman and Hall, London, pp. 227-247, 1993a.
Nolet, G., A Volksseismometer?,IRIS Ne•vsl., •,
10-11,
1993b.
Nolet, G., and M.J.R. Wortel, Mantle, upper structure, in Encyclopediao• Solid Earth Geophysics,edited by D.E. James, pp. 775-788, Van Nostrand Reinhold, New York, 1989.
Nolet, G., and A. Zielhuis, Low S velocities under the Tornquist-Teisseyrezone: Evidence for water injection
into the transition zone by subduction,3. Geophys.Res., 99, 15, 813-15, 820, 1994.
Okal, E.A., The effect of intrinsic oceanicupper-mantle het-
erogeneityon regionalizationof long-periodrayleighwave phasevelocities,Geophys.3. R. Astron. Soc.,49, 357-370, 1977.
tron. Soc., 8J, 581-605, 1986. Snieder, R., J. Beckets, and F. Neele, The effect of smallscale structure on normal mode frequencies and global inversions, 3. Geophys. Res., 96, 501-515, 1991. Spakman, W., Delay-time tomography of the upper mantle below Europe, the Mediterranean, and Asia Minor, Geophys. 3. Int., •07, 309-332, 1991. Stammler, K., R. Kind, G.L. Kosarev, A. Piesinger, J. Horalek, L. Qiyuan, and L.P. Vinnik, Broadband observations of PS conversionsfrom the upper mantle in Eurasia, Geophys. 3. Int., •05, 801-804, 1991. Stammler, K., R. Kind, N. Petersen, G. Kosarev, L. Vinnik, and L. Qiyuan, The upper mantle discontinuities: Correlated or anticorrelated?, Geophys. Res. œett., •9, 15631566, 1992.
Su, W.-J., and A.M. Dziewonski, Predominance of longwavelengthheterogeneityin the mantle, Nature, 35•, 121126, 1991.
Su, W.-J., R.L. Woodward, and A.M. Dziewonski, Deep origin of mid-ocean-ridge seismicvelocity anomalies, Nature,
359, 149-152, 1992. Park, J., and Y. Yu, Anisotropy and coupled free oscillations: Simplifiedmodels and surfacewave observations, Su, W.-J., R.L. Woodward, and A.M. Dziewonski,Degree12 model of shear velocity heterogeneity in the mantle, Geophys.3. Int., •0, 401-420, 1992. 3. Geophys. Res., 99, 6945-6981, 1994. Paulssen,H., Lateral heterogeneityof Europe'supper mantle as inferredfrom modellingof broad-bandbody waves, Tarantola, A., Inverse Problem Theory: Methods •or Data Fitting and Model Parameter Estimation, Elsevier, New Geophys.3. R. Astron. Soc., 9•, 171-199, 1987. York, 1987. Paulssen,H., Evidencefor a sharp 670-km discontinuityas inferred from P-to-S convertedwaves,3. Geophys.Res., Taylor, S.R., and M.N. Toksoz, Three-dimensional crust and
93, 10,489-10,500, 1988.
upper mantle structure of the northeastern United States,
3. Geophys. Res., 84, 7627-7644, 1979. Regan, J., and D.L. Anderson,Anisotropicmodelsof the upper mantle, Phys. Earth Planet. Inter., 35, 227-263, Tong C., O. Gudmundsso,and B.L.N. Kennett, Shear-waves 1984.
splitting in refracted wavesreturned from the upper man-
23,766
NOLET ET AL.: UPPER MANTLE HETEROGENEITY
tle transition zone beneath northern Australia, Y. OeoRayleigh waves along the East Pacific Rise: Evidence for S wave velocity anomaliesto 450 km depth, J. Geophys. Res., 99, 15, 783-15, 799, 1994. phys. Res., 87, 8631-8641, 1982. Tryggvason, K., E.S. Husebye, and R. Stefansson,Seismic imaging of the hypothesized Icelandic hotspot, Tectono- Wiggins,R. A., and D. V. Helmberger,Upper mantlestructure of the western United States, J. Geophys.Res., 78, physics, 100, 97-118, 1983. 1870-1880, 1973. VanDecar, J.C., and R.S. Crosson, Determination of teleseismic relative phase arrival times using multi-channel Woodhouse,J.H., and A.M. Dziewonski,Mapping the upper mantle: Three-dimensional modeling of the Earth struccross-correlationand least squares, Bull. Seismol. Soc. Am., 80, 150-169, 1990.
Van der Hi]st, R., and T. Seno, Effects of relative plate motion on the deep structure and penetration depth of slabs below the Izu-Bonin and Marian island arcs, Earth Planet. Sci. œett., 1•0, 395-407, 1993. Van der Hi]st, R., R. Engdahl, W. •pakman, and G. Nolet, Tomographic imaging of subducted lithosphere below northwest Pacific island arcs, Nature, gSg, 37-43, 1991. Van der Lee, S., H. Paulssen, and G. Nolet, Variability of P660s phases as a consequenceof topography?, Phys. Earth Planet. Inter., in press, 1994. Van der Sluis, A., and H.A. Van der Vorst, Numerical solution of large, sparselinear algebraic systemsarising from tomographic problems, in: Seismic Tomography,edited by G. Nolet, pp. 49-84, D. Reidel, Norwell, Mass., 1987. Vasco, D.W., R.J. Pullism, and L.R. Johnson,Formal inversion of ISC arrival times for mantle P velocity structure, Geophys. J. Int., 113, 586-606, 1993. Vidale, J. E., and H. M. Benz, Upper-mantle seismicdiscontinuities and the thermal structure of subduction zones, Nature, 356, 678-683, 1992. Vinnik, L.P., Detection of waves converted from P to S Vin the mantle, Phys. Earth Planet. Ir.ter., 15, 39-45, 1977. Walck, M. C., The P-wave upper mantle structure beneath an active spreading center: The Gulf of California, Geophys. J. R. Astron. Soc., 76, 697-723, 1984. Whircomb, J., and D. L. Anderson, Reflection of P'P' seismic waves from discontinuities in the mantle, J. Geophys. Res., 75, 5713-5728, 1970. Wielandt, E., On the validity of the ray approximation for interpreting delay times, in Seismic Tomography,edited by G. Nolet, pp. 85-98, D. Reidel, Norwell, Mass., 1987. Wielandt, E., and L. Knopoff, Dispersionof very long-period
ture by inversionof seismicwaveforms,•. Geophys.Res., 89, 5953-5986, 1984.
Woodward, R.L., and G. Masters, Global upper mantle structurefrom long-perioddifferentialtraveltimes, J. Geophys. Res., 96, 6351-6378, 1991.
Wu, R.-S., and S.M. Flatte, Transmissionfluctuationsacross an array and heterogeneitiesin the crust and upper mantle, Pure Appl. Geophys.,13•, 175-196, 1990. Zhang, Y.-S., and T. Tanimoto, Ridges,hotspotsand their interpretation as observedin seismicvelocity maps, Nature, 355, 45-49, 1992.
Zhao, D., and A. Hasegawa, P wave tomographicimaging of the crust and upper mantle beneath the Japan islands, J. Geophys. Res., 98, 4333-4353, 1993.
Zhao, L. S., and D. V. Helmberger,Upper mantle compressional velocity structure beneath the northwest Atlantic Ocean, J. Geophys.Res., 98, 14,185-14,196, 1993. Zhou, H., and R.W. Clayton, P and S wave travel time inversionsfor subductingslab under the island arcsof the northwest Pacific, J. Geophys.Res., 95, 6829-6851, 1990.
Zielhuis, A., and G. Nolet, Shear wave velocity variations in the upper mantle beneath central Europe, Geophys.J. Int., 117, 695-715, 1994.
S.P. Grand, Department of Geological Sciences,University of Texas at Austin, Austin TX 78713. B.L.N. Kennett, ResearchSchoolof Earth Sciences,Australian National University, Canberra ACT 0200, Australia. G. Nolet, Department of Geologicaland GeophysicalSciences,Princeton University, Princeton NJ 08544.
(ReceivedMarch 7, 1994;revisedJuly 15, 1994; acceptedJuly 19, 1994.)