Seismic Travel Time Evidence for Lateral Inhomogeneity in the Deep Mantle

B. R. Julian and M. K. Sengupta

Editor’s Note

The Earth’s lower mantle is inaccessible for direct study. The best method for probing it is through the measurement of seismic waves, which travel at different speeds through materials of different density or composition. Although there were some earlier signs that the lower mantle is not uniform but somehow “lumpy”, Bruce Julian and Mrinal Sengupta at the Massachusetts Institute of Technology here provide clear evidence from seismic travel times that the lowest part of the mantle is inhomogeneous on scales of around 1,000 km. They could not offer definitive reasons, but suggest that at least some of the irregularities might be caused by convection structures: specifically, by plumes of rising, hot material, which others had proposed earlier. 中文

EVIDENCE from seismic body wave and surface wave data has long indicated that the Earth’s upper mantle (depth <700 km) is strongly laterally heterogeneous. Lateral variations of the compressional wave velocity as large as 10% have been reported for the upper 200 km, and the shear velocity probably varies even more. For depths greater than 700 km, however, the existence of lateral variations has been more difficult to establish, although such variations have been invoked by various workers ^1,2 to explain the scatter of some seismological data. Greenfield and Sheppard ³ , in a study of d T /d∆ measurements made at the Large Aperture Seismic Array in Montana, found a pronounced difference between data from events to the northwest and the southeast for epicentral distances greater than 60°, which could not be attributed to the structure beneath the array, and seems explicable only in terms of heterogeneities in the lower mantle. Davies and Sheppard ⁴ have presented a more extensive collection of data of this type in the form of an “array diagram”, on which d T /d∆ and azimuth anomalies are represented as vectors in slowness space. Many anomalies are found which are too large to be effects of upper mantle heterogeneities in the source regions; on the other hand, the anomalies often vary rapidly with the direction of approach of the waves, implying that structure directly beneath the array is not responsible. Further evidence has come from a study of the diffraction of compressional waves by the Earth’s core, in which Alexander and Phinney ⁵ found that the region of the core-mantle boundary beneath the Pacific Basin is distinctly different from the region beneath the North Atlantic and Africa. 中文

Travel Time Anomalies

We have found evidence that significant lateral variations occur in the lowest few hundred km of the mantle, this region being much more heterogeneous than that which lies above it. The methods of this travel time study are described in detail elsewhere ⁶ . 中文

We restricted the study to data from deep focus earthquakes in an attempt to avoid systematic errors caused by near source velocity variations in the upper mantle, which are particularly severe in seismically active regions. About 3,300 arrival time data from 47 earthquakes with depths between 450 and 650 km were used, all events being located in the deepest parts of their respective seismic zones. For most of the 18 seismic regions involved, two or more events were available, and for each station a consistency check was made to eliminate data contaminated by gross errors. An iterative procedure, similar to the one described by Herrin et al. ⁷ , was then used to determine the earthquake locations, the travel time curve for a 550 km focal depth, and a set of “station corrections”(each represented by a constant) to account for the effect of lateral variations in the upper mantle beneath the stations. 中文

The travel time curve thus determined is shown in Fig. 1 (in terms of deviation from the standard Jeffreys-Bullen tables), together with the data means and their standard errors for 2°distance intervals. This curve is similar in shape to those found in other recent studies ^7,8 except beyond 85°, where most other curves remain approximately parallel to the Jeffreys-Bullen curve, but ours becomes progressively earlier by about 0.9 s. Fig. 2 shows some of the data which have contributed to the determination of the travel time curve; the observed times show a striking dependence upon the location of the earthquakes. The difference between the curves in Fig. 1 results from differences in the geographic regions sampled. Possible causes of a regional dependence of this nature are velocity variations in the upper mantle in the source or receiver regions or in the lower mantle, and mislocation of the events (caused by uneven station distribution, and so on). Event mislocation and structure in the receiver regions can be ruled out because they would be expected to produce similar effects at all epicentral distances. Although the observed variations are most striking beyond 85°, we shall show that they are much smaller at distances less than 70°. It is conceivable that structure in the source regions could produce a distance dependent regional variation of this kind, if the velocity anomalies were systematically located relative to the earthquake hypocentres (as indeed they are beneath island arcs). In that case the variations would have to be localized in a very small region beneath the hypocentres, because a 10° distance interval maps into about a 5° difference in angle at the focus. Even if the anomalous regions are as deep as 1,000 km, the velocity change must occur over a horizontal distance of only 50 km or so. This possibility may be ruled out because all the earthquakes in each source region yield a similar pattern of travel time residuals, even though the epicentre locations in each region are typically distributed over >200 km. Velocity variations near the focus are further ruled out because early arrivals beyond 85° are not restricted to observations of deep earthquakes; they also occur, for example, in data from nuclear explosions in the Marshall Islands ⁹ . 中文

中文

Deep Mantle Structure

It seems, then, that lateral variations of compressional velocity in the middle or lower mantle are required to explain the travel time anomalies. But because of the uneven distribution to seismological observatories and deep earthquakes, the sampling of the mantle provided by available data is uneven, and it is impossible to determine uniquely the complete three dimensional velocity structure of the mantle. What can be determined is the average travel time residual for each of a number of “bundles” of rays following nearly identical paths from a seismic region to a group of stations, and from this information we infer the most probable cause of the variations. Table 1 summarizes the travel time data for all paths for which 9 or more observations are available. For each path a Student’s t test has been used to evaluate the hypothesis that the mean travel time (after station corrections have been applied) is the value given by the curve in Fig. 1 and that deviations from this curve can be attributed to random measuring errors. Those ray paths for which the hypothesis could be rejected at the 99.5% confidence level are indicated in Table 1. Fig. 3 shows histograms of the residual distribution for these anomalous paths. For observations at distances beyond 70°, 16 paths (out of 34 tested) showed significant variations from the average curve, whereas for smaller distances, only 3 anomalous paths (out of 22) were found. This strongly suggests that most of the scatter originates in the deep mantle (depth >2,000 km). The possibility of the variations occurring at a shallower depth cannot be absolutely disproved, but the velocity distribution in the Earth would have to be such that, given the distribution of earthquakes and seismic stations, all observed waves bottoming at the depth of the heterogeneities happen to be unaffected by them (even though P waves spend about 25% of their travel time traversing the lowest 10% of the ray path), while waves penetrating beneath the heterogeneities are affected. It seems unreasonable to assume such a conspiratorial behaviour on the part of the velocity variations when a much simpler hypothesis is available. 中文

中文

N =Number of observations.

=Mean travel time residual after station correction (s).

s =Standard deviation of residuals (s).

t =

t _99.5 =99.5% confidence limit for | t | if true mean is zero.

* Indicates paths with mean significantly different from zero.

中文

中文

The actual details of the velocity distribution cannot, however, be determined precisely, because rays travel a great distance at approximately the same depth near their turning points. Fig. 4 shows the regions of the lower mantle sampled by the various ray bundles (defined arbitrarily as the central 30°of each path) and indicates which paths correspond to early and late arrivals. These are in most cases probably the regions where the actual velocity anomalies occur. Where regions overlap on the figure, they are generally consistent with each other (for example, regions 1-6 and 1-7, regions 2-2 and 2-3, and regions 4-1 and 4-6). This consistency is encouraging, in that it supports our argument that the travel time anomalies originate in the deep mantle and are not the result of some other type of systematic error. An apparent inconsistency exists between regions 1-3 and 1-16, but this is not surprising because of the uncertainty as to exactly where the travel time anomalies actually originate. The rays following path 1-3 also pass through regions 2-2 and 2-3 further to the north, and it is likely that the travel time anomaly actually originates there. Another interesting feature of Fig. 4 is a correlation between the anomalies in the two greatest distance ranges (for example, regions 1-7 and 2-7, regions 1-4 and 2-4, and regions 1-5 and 2-6), suggesting that the structure in the deep mantle has a spatial “coherence” of at least a few hundred km vertically. 中文

中文

The mean travel times in Table 1 show a variation of about 1.5 s for rays bottoming below 2,600 km, and about 0.6 s for rays bottoming between 2,000 and 2,600 km. These numbers are somewhat uncertain, but it is likely that the true travel times vary by at least 1 s. The amount by which the actual velocity in the mantle varies depends on the size of the regions within which the variations occur. The data of Fig. 4 suggest that the size of some of the anomalies, at least, is about 1,000 km or less, in which case the velocity must vary by at least 1%. This is a lower bound both because we have probably overestimated the scale of the inhomogeneities and because we are measuring averages of the velocity in rather large regions and some cancellation of the effects of positive and negative velocity anomalies is likely. Combined interpretation of travel time and d T /d∆ measurements can probably improve the resolution of structural details. 中文

Might the deep mantle variations be related in some way to the convection plumes hypothesized by Wilson ¹⁰ and Morgan ¹¹ to exist in the deep mantle? To answer this the region of the Hawaiian Islands provides the best data, and here they indeed indicate a pronounced lateral variation, the velocities being high to the northwest of Hawaii (regions 2-2, 2-3, and probably 1-3) and low in the vicinity of the islands (regions 1-4, 2-4, and 1-6). Interestingly, the d T /d∆ data presented by Davies and Sheppard ⁴ also indicate a horizontal velocity contrast of this sort in the vicinity of Hawaii. Unfortunately, no other proposed plumes are well sampled by our data. Region 1-13 includes the Mascarene Islands, but the data here are highly scattered, and no conclusion can be drawn. Regions 1-5 and 2-6, both with apparently high velocities, are located slightly to the northeast of the Galapagos Islands, so it is not clear what relation, if any, this velocity anomaly may bear to a possible plume. If travel time data for rays passing through more proposed plume regions can be obtained, they may provide valuable evidence relevant to the Wilson-Morgan hypothesis. 中文

Except, perhaps, for the Hawaiian Islands, no geological or tectonic features show an obvious correlation with the inferred deep mantle velocity anomalies. At shallower depths, however, this is not the case; regions 3-1, 4-1, and 4-6, all seeming to have low velocities, lie beneath the Kurile and Aleutian Island arcs. The only other island arc adequately sampled at these depths lies beneath Middle America (region 4-2) and is associated with early arrivals (though they are not significant at the 99.5% confidence level). The data thus suggest that low velocities may be characteristic of island arcs at depths greater than 1,000 km. 中文

The data considered here do not support any correlation between velocity variations below 2,000 km and global gravity anomalies or geoid heights. At shallower depths such a correlation does exist, but it is merely another manifestation of the low velocities beneath the Kurile and Aleutian Islands, because the concave sides of island arcs are generally the sites of prominent positive free air gravity anomalies. 中文

It is not possible to make a direct comparison between these results and those of Alexander and Phinney ⁵ . The region of the North Atlantic found by them to be anomalous is further east than the corresponding region sampled by our data. Further, travel time measurements such as ours provide a measure of the average velocity in a region, whereas the behaviour of core diffracted waves depends on features such as the velocity gradient in the lower mantle. Further studies of variations in the “visibility” within the core shadow would be a useful complement to travel time and d T /d∆ studies of the lower mantle. 中文

We thank Dr. David Davies and Dr. M. Nafi Toksöz for helpful suggestions. This work was sponsored by the Advanced Research Projects Agency of the Department of Defense. 中文

( 242 , 443-447; 1973)

Bruce R. Julian and Mrinal K. Sengupta

Lincoln Laboratory, Massachusetts Institute of Technology, Lexington, Massachusetts 02173

Received January 2, 1973.

---

References: U1IF9DiDOT8e+bTtqtqlbySXtNIN/L85SoniywazujchU5Mjf11WsdpTJ6mMrQLj

1. Toksöz, M. Nafi, Chinnery, M. A., and Anderson, D. L., Geophys. J. Roy. Astron. Soc ., 13 , 31(1967).
2. Hales, A. L., and Roberts, J. L., Bull, Seismol. Soc. Amer. , 60 , 1427(1970).
3. Greenfield, R. J., and Sheppard, R. M., Bull. Seismol. Soc. Amer. , 59 (1969).
4. Davies, D., and Sheppard, R. M., Nature , 239 , 318(1972).
5. Alexander, Shelton S., and Phinney, R. A., J . Geophys. Res. , 7 , 5943(1966).
6. Sengupta, Mrinal K., thesis, MIT(1972).
7. Herrin, E., Tucker, W., Taggart, J., Gordon, D. W., and Lobdell, J. L., Bull. Seismol. Soc. Amer. , 58 , 1273(1968).
8. Lilwal, R. C., and Douglas, A., Geophys. J. Roy. Astron. Soc. , 19 , 165(1970).
9. Carder, D. S., Gordon, D. W., and Jordan, J. N., Bull. Seismol. Soc. Amer. , 56 , 815(1966).
10. Wilson, J. T., Phil. Trans. Roy. Soc. , A, 258 , 145(1965).
11. Morgan, W. J., Nature , 230 , 42(1971).