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.     

Localized velocity anomalies in the lower mantle

Summary. Two localized regions of velocity heterogeneity in the lower mantle with scale lengths of 1000–2000 km and 2 per cent velocity contrasts are detected and isolated through comparison of S, ScS, P and PcP travel times and amplitudes from deep earthquakes in Peru, Bolivia, Argentina and the Sea of Okhotsk. Comparison of the relative patterns of ScS-S differential travel times and S travel-time residuals across North American WWSSN and CSN stations for the different source regions provides baselines for interpreting which phases have anomalous times. A region of low S and P velocities is located beneath Northern Brazil and Venezuela at depths of 1700–2700 km. This region produces S -wave delays of up to 4 s for signals from deep Argentine events recorded at eastern North American stations. The localized nature of the anomaly is indicated by the narrow bounds in azimuth (15°) and take-off angle (13°) of the arrivals affected by it. The long period S -waves encountering this anomaly generally show 30–100 per cent amplitude enhancement, while the short-period amplitudes show no obvious effect. The second anomaly is a high-velocity region beneath the Caribbean originally detected by Jordan and Lynn, who used travel times from deep Peruvian events. The data from Argentine and Bolivian events presented here constrain the location of the anomaly quite well, and indicate a possible short- and long-period S -wave amplitude diminution associated with it. When the travel-time data are corrected for the estimated effects of these two anomalies, a systematic regional variation in ScS-S station residuals is apparent between stations east of and west of the Rocky Mountains. One possible explanation of this is a long wavelength lateral variation in the shear velocity structure of the lower mantle at depths greater than 2000 km beneath North America.     

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.     

Upper-mantle structure of the Baltic Shield below the Swedish National Seismological Network (SNSN) resolved by teleseismic tomography

Upper-mantle structure under the Baltic Shield is studied using non-linear high resolution teleseismic P -phase tomography. Observed relative arrival-time residuals from 52 teleseismic earthquakes recorded by the Swedish National Seismological Network (SNSN) are inverted to delineate the structure of the upper mantle. The network consists of 47 (currently working) three-component broad-band stations located in an area about 450 km wide and 1450 km long. In order to reduce complications due to possible significant three-dimensionality of Earth structure, events chosen for this study lay close to in-line with the long-axis of the array  (±30°)  . Results indicate P -wave velocity perturbations of ±3 per cent down to at least 470 km below the network. The size of the array allows inversion for structures even at greater depths, and lateral variations of velocity at depths of up to 680 km appear to be resolved. Below the central part of the array (60°–64° N), where ray coverage is best, the data reveals a large region of relatively low velocity at depths of over about 300 km. At depths less than about 250–300 km, the models include a number of features, including an apparent slab-like structure dipping gently towards the north.     

Three-dimensional image of the anisotropic bodies beneath central Honshu, Japan

Shear-wave splitting from local deep earthquakes is investigated to clarify the volume and the location of two anisotropic bodies in the mantle wedge beneath central Honshu, Japan. We observe a spatial variation in splitting parameters depending on the combination of sources and receivers, nearly N–S fast in the northern region, nearly E–W fast in the southern region and small time delays in the eastern region. Using forward modelling, two models with 30 and 10 per cent anisotropy are tested by means of a global search for the locations of anisotropic bodies with various volumes. The optimum model is obtained for 30 per cent anisotropy, which means a 5 per cent velocity difference between fast and slow polarized waves. The northern anisotropic body has a volume of 1.00° (longitude) × 0.5° (latitude) × 75 km (depth), with the orientation of the symmetry axis being N20°E. The southern anisotropic body has a volume of 1.25° × 1.25° × 100 km with the symmetry axis along N95°E. Our results show that the anisotropic bodies are located in low-velocity and low- Q regions of the mantle. This, together with petrological data and the location of volcanoes in the arc, suggests that the possible cause of the anisotropy is the preferred alignment of cracks filled with melt.     

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.     

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.     

The shear-wave velocity structure in the upper mantle beneath Eurasia

We develop an approach that allows us to invert for the mantle velocity structure within a finely parametrized region as a perturbation with respect to a low-resolution, global tomographic model. We implement this technique to investigate the upper-mantle structure beneath Eurasia and present a new model of shear wave velocity, parametrized laterally using spherical splines with ∼2.9° spacing in Eurasia and ∼11.5° spacing elsewhere. The model is obtained from a combined data set of surface wave phase velocities, long-period waveforms and body-wave traveltimes. We identify many features as narrow as few hundred kilometres in diameter, such as subducting slabs in eastern Eurasia and slow-velocity anomalies beneath tectonically active regions. In contrast to regional studies in which these features have been identified, our model encompasses the structure of the entire Eurasian continent. Furthermore, including mantle- and body-wave waveforms helped us constrain structures at depths larger than 250 km, which are poorly resolved in earlier models. We find that up to +9 per cent faster-than-average anomalies within the uppermost ∼200 km of the mantle beneath cratons and some orogenic regions are separated by a sharp gradient zone from deeper, +1 to +2 per cent anomalies. We speculate that this gradient zone may represent a boundary separating the lithosphere from the continental root, which might be compositionally distinct from the overlying lithosphere and remain stable either due to its compositional buoyancy or due to higher viscosity compared with the suboceanic mantle. Our regional model of anisotropy is not significantly different from the global one.     

Anelasticity of the Crust and Upper Mantle of South America From the Inversion of Observed Surface Wave Attenuation

Fundamental-mode Rayleigh and Love waves generated by several earthquakes situated along great-circle paths between pairs of seismograph stations have been analysed to obtain coefficients of attenuation, group velocities, phase velocities, and specific quality factors in the period range 18–80s in two regions of the South American continent. One set of paths crosses the shield region which lies on the eastern coast and another set traverses the mountainous region inland. the average attenuation coefficient values are clearly higher in the tectonically active western region throughout the entire period range than in the eastern or shield region.
Inversion of the attenuation data yielded shear wave internal friction ( Q ^-1_β) models as a function of depth in the crust and upper mantle in both regions. A low- Q zone below the lithosphere is prominent in both regions. the results show that substantial variations of Q _β occur in the two regions of South America. the Q_β values were found to be inversely related to the heat flow values or to the temperature.     

Seismotectonics and analysis of earthquakes from the Hindukush region

Summary. Over 80 earthquakes, exclusively from the Hindukush focal region, which were recorded at the Gauribidanur seismic array (GBA) have been used in this study. These events have similar epicentral distances and a narrow azimuthal range from GBA but varying focal depths from 10 to 240 km. A fault plane dipping steeply (75°) in the north-west direction and striking N 66° E has been investigated on the basis of the spatial distribution of earthquakes in two vertical planes through 68° E and 32° N. Short period P -wave recordings up to 30 s were processed using the adaptive cross-correlation filtering technique. Slowness and azimuthal anomalies were obtained for first arrivals. These anomalies show positive as well as negative bias and are attributed to a steep velocity gradient in the upper mantle between the 400–700 km depth range where the seismic rays have their maximum penetration. Relative time residuals between the stations of GBA owe their origin very near to the surface beneath the array. A search of the signals across the array revealed that most of the events occurring at shallower depths had complex signatures as compared to the deeper events. The structure near the source region, complicated source functions and the scattering confined to the crust—upper mantle near source are mainly responsible for the complexity of the Hindukush earthquakes as the transmission zone of the ray tubes from turning point to the recording station is practically the same.     

A lower mantle S-wave triplication and the shear velocity structure of D"

Summary. A lower mantle S-wave triplication detected with short- and long-period WWSSN and CSN recordings indicates a substantial shear velocity discontinuity near 280 km above the core–mantle boundary. The triplication can be observed in rotated SH seismograms from intermediate and deep focus events throughout the distance range from 70° to 95°. Three distinct source region–receiver array combinations that have been investigated in detail demonstrate consistent travel time and relative amplitude behaviour of the triplication, with slight systematic shifts in the triplication indicating up to 40 km variations in the depth of the discontinuity. Modelling of the observations with synthetic seismograms produced with the Cagniard de Hoop and reflectivity methods constrains the shear velocity increase to be 235 ± 0.25 per cent, comparable to upper mantle discontinuities. Short-period observations indicate that the velocity increase may be a sharp first-order discontinuity, or may extend over a transition zone no more than 50 km thick. The shear velocity gradient below the discontinuity, within the D" layer, is not well-constrained by the SH data, but slightly positive or near zero velocity gradients are consistent with the long-period amplitude ratios of ScSH/SH .     

Shear and compressional velocity models of the mantle from cluster analysis of long-period waveforms

We present a new technique for the efficient measurement of the traveltimes of long period body wave phases. The technique is based on the fact that all arrivals of a particular seismic phase are remarkably similar in shape for a single event. This allows the application of cross-correlation techniques that are usually used in a regional context to measure precise global differential times. The analysis is enhanced by the inclusion of a clustering algorithm that automatically clusters waveforms by their degree of similarity. This allows the algorithm to discriminate against unusual or distorted waveforms and makes for an extremely efficient measurement technique.
This technique can be applied to any seismic phase that is observed over a reasonably large distance range. Here, we present the results of applying the algorithm to the long-period channels of all data archived at the IRIS DMC from 1976 to 2005 for the seismic phases S and P (from 23° to 100°) and SS and PP (from 50° to 170°). The resulting large data sets are inverted along with existing surface wave and updated differential traveltime measurements for new mantle models of S and P velocity. The resolution of the new model is enhanced, particularly, in the mid-mantle where SS and PP turn. We find that slow anomalies in the central Pacific and Africa extend from the core–mantle boundary to the upper mantle, but their direct connection to surface hotspots is beyond our resolution. Furthermore, we find that fast anomalies that are likely associated with subducting slabs disappear between 1700 and 2500 km, and thus are not continuous features from the upper to lower mantle despite our extensive coverage and high resolution of the mid-mantle.     

The upper mantle structure of the North Sea region from Rayleigh wave dispersion

Summary. The phase velocity dispersion of fundamental mode Rayleigh waves (period range 13–127 s) is determined by the interstation method for three profiles that traverse the North Sea region of northwest Europe. The resulting observations have been combined to produce a regional phase velocity curve with 95 per cent confidence intervals, which belongs to the aseismic continental platform category of Knopoff.
Inversions of the regional phase velocity curve by the'Hedgehog'method indicate that the North Sea region is characterized by an upper mantle low-velocity zone of S -wave velocity 4.35–4.45 km/s between depths of approximately 85–200 km.     

A three-dimensional radially anisotropic model of shear velocity in the whole mantle

We present a 3-D radially anisotropic S velocity model of the whole mantle (SAW642AN), obtained using a large three component surface and body waveform data set and an iterative inversion for structure and source parameters based on Non-linear Asymptotic Coupling Theory (NACT). The model is parametrized in level 4 spherical splines, which have a spacing of ∼ 8°. The model shows a link between mantle flow and anisotropy in a variety of depth ranges. In the uppermost mantle, we confirm observations of regions with   V_SH > V_SV   starting at ∼80 km under oceanic regions and ∼200 km under stable continental lithosphere, suggesting horizontal flow beneath the lithosphere. We also observe a   V_SV > V_SH   signature at ∼150–300 km depth beneath major ridge systems with amplitude correlated with spreading rate for fast-spreading segments. In the transition zone (400–700 km depth), regions of subducted slab material are associated with   V_SV > V_SH   , while the ridge signal decreases. While the mid-mantle has lower amplitude anisotropy (<1 per cent), we also confirm the observation of radially symmetric   V_SH > V_SV   in the lowermost 300 km, which appears to be a robust conclusion, despite an error in our previous paper which has been corrected here. The 3-D deviations from this signature are associated with the large-scale low-velocity superplumes under the central Pacific and Africa, suggesting that   V_SH > V_SV   is generated in the predominant horizontal flow of a mechanical boundary layer, with a change in signature related to transition to upwelling at the superplumes.     

Travel-time anomalies in the mantle under the North Atlantic

Summary. Lateral heterogeneity exists in the Earth's mantle, and may result in seismic velocity anomalies up to several per cent. If convection cells and plumes extend down to the core, then these features may be associated with local inhomogeneities observed in the lower mantle.
Published data for direct and core-reflected P -wave residuals are used to delineate velocity anomalies in the middle—lower mantle under the North Atlantic. Differential ( PcP — P ) residuals indicate travel-time anomalies near the core—mantle transition, and may be due to core topography or lateral variations in velocity. It is assumed that the anomalies occur near the midpoints of the ray paths. The main source of error in the data set may arise from phases which have been identified incorrectly. Hence trend-surfaces are fitted to the residual data to show only the large-scale trends in anomaly values, with wavelengths of the order of 1000 km.
The Azores and Colorado hot spots occur in a region covered by the data. A possible interpretation of the trend maps is that an anomalous zone extends from a relatively fast region at the core boundary at 35° N, 50° W up to these hot spots, at about 30 degrees from the vertical. This may agree with the suggestion of Anderson that plumes are chemical rather than thermal in origin. If inclⁱned plumes do exist, the deviation from the ideal vertical plume or convection cell boundary may imply that lateral shear or other distortion effects exist in the mantle.     

Three-dimensional seismic structure beneath the Australasian region from refracted wave observations

The earthquakes in the seismicity belt extending through Indonesia, New Guinea, Vanuatu and Fiji to the Tonga–Kermadec subduction zone recorded at the 65 portable broad-band stations deployed during the Skippy experiment from 1993–1996 provide good coverage of the lithosphere and mantle under the Australian continent, Coral Sea and Tasman Sea.
The variation in structure in the upper part of the mantle is characterized by deter-mining a suite of 1-D structures from stacked record sections utilizing clear P and S arrivals, prepared for all propagation paths lying within a 10° azimuth band. The azimuth of these bands is rotated by 20° steps with four parallel corridors for each azimuth. This gives 26 separate azimuthal corridors for which 15 independent 1-D seismic velocity structures have been derived, which show significant variation in P and S structure.
The set of 1-D structures is combined to produce a 3-D representation by projecting the velocity values along the ray path using a turning point approximation and stacking into 3-D cells (5° by 50 km in depth). Even though this procedure will tend to underestimate wave-speed perturbations, S -velocity deviations from the ak135 reference model exceed 6 per cent in the lithosphere.
In the uppermost mantle the results display complex features and very high S -wave speeds beneath the Precambrian shields with a significant low-velocity zone beneath. High velocities are also found towards the base of the transition zone, with high S -wave speeds beneath the continent and high P -wave speeds beneath the ocean. The wave-speed patterns agree well with independent surface wave studies and delay time tomography studies in the zones of common coverage.     

Regional 3-D heterogeneities by waveform inversion – application to the Atlantic area

The waveform inversion method described in Woodhouse & Dziewonski (1984) was modified to retrieve regional scale 3-D heterogeneities by using the minor arc part of seismograms. The lateral heterogeneities are expanded horizontally into blocks (10°× 10°) and radially into Legendre polynomials up to order 3 (0–670 km), and thus the results show much fine details of lateral variation than previous global scale studies. We assumed that the heterogeneities produce the perturbation of eigenfrequencies which are the minor arc average of local eigenfrequency shift. We applied the method to the upper mantle beneath the Atlantic Ocean and its environments. Care was taken about the weighting of the data set. We found that the fit of each seismogram became better when the weighting of each seismogram is proportional to the inverse of initial data residuals. Resolution is good in the triangular region surrounded by South America, Europe, and North America. Resolution is not good in the South Atlantic because of the poor path coverage. Depth resolution is not clear, because of the use of Legendre polynomials, though the results suggest a broad half-width of the order of 200 km or more. We found some similarities between previous global studies and our results. For example, low velocities beneath the East Pacific Rise, Chile Rise and Azores triple junction and a high velocity Canadian shield are obtained. However, there are also differences; the high-velocity zone beneath the Brazilian shield at shallow depth is not a prominent feature in this study. Instead, we found a somewhat unexpected feature near the Romanche and Vema fracture zones where shallow positive anomalies exist. Smoothed results calculated by the spherical harmonic expansion are also shown for the purpose of comparison with global studies.     

Experimental study of the anisotropy of longitudinal and transverse waves from local earthquake records

Summary. The paper gives the results of a study of the anisotropy of seismic wave velocities within the Ashkhabad test field in Central Asia. The anisotropy was studied by analysing variations in the values of apparent velocities of first arrivals for epicentral distances ranging from 30 to 130 km and by analysing the delays (Δ t_s1-_s2 ) between the arrival times of shear waves with different polarizations.
The velocities of P -waves vary with azimuth from 5.3 to 6.27 km s^-1 and the velocities of S -waves vary from 3.15 to 3.5 km s^-1.
The delay times Δ t_S1 - _S2 depend on the direction of the propagation. The character of the variation of the propagation velocity of the longitudinal wave, the presence of two differently polarized shear waves S 1 and S 2 propagating at different velocities, and the character of the distribution of Δ t_S1 - _S2 on the stereogram suggest that the symmetry of the anisotropic medium is close to hexagonal with a nearly horizontal symmetry axis coinciding with the direction of maximal velocity. The azimuth of the symmetry axis of the medium is 140° and coincides with the direction of geological faults.     

The three-dimensional shear wave structure in the mantle by overtone waveform inversion - I. Radial seismogram inversion

Summary. The three-dimensional (3-D) shear wave structure of the mantle, down to the depth of about 900 km, is obtained by inverting waveforms of radial component seismograms. Radial component seismograms contain large amplitude overtone signals which circle the Earth as wave packets and are sometimes called X1, X2, X3, … We use data which contain R1, X1 and X2 and filtered between 2 and 10mHz. It is shown that, unless each seismogram is weighted, all seismograms are not fitted uniformly. Only data from large earthquakes are fitted and the final velocity anomalies are biased by the small number of large earthquake data. Resolution is good at shallow depths, becomes worse in the intermediate depth range between about 400 and 500 km and then becomes better at greater depth ranges (600–900km). Even though we use only spheroidal mode data, velocity anomalies in the shallow structure show excellent correlation with the age of the surface rocks of the Earth. In the deeper regions, between about 600 and 900km, South America shows a fast velocity anomaly which may indicate the slab penetration beyond 700 km there. Another region which shows a fast velocity anomaly is the Mariana trench, but other subduction regions do not show such features.     

Crustal and upper mantle structure of the northwestern North Island, New Zealand, from seismic refraction data

The crustal and upper mantle structure of the northwestern North Island of New Zealand is derived from the results of a seismic refraction experiment; shots were fired at the ends and middle of a 575 km-long line extending from Lake Taupo to Cape Reinga. The principal finding from the experiment is that the crust is 25 ± 2 km thick, and is underlain by what is interpreted to be an upper mantle of seismic velocity 7.6 ± 0.1 km s⁻¹, that increases to 7.9 km s⁻¹ at a depth of about 45 km. Crustal seismic velocities vary between 5.3 and 6.36 km s⁻¹ with an average value of 6.04 km s⁻¹. There are close geophysical and geological similarities between the north-western North Island of New Zealand and the Basin and Range province of the western United States. In particular, the conditions of low upper-mantle seismic velocities, thin crust with respect to surface elevation, and high heat-flow (70–100 mW m⁻²) observed in these two areas can be ascribed to their respective positions behind an active convergent margin for about the past 20 Myr.