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.     

Earthquake processes of the Himalayan collision zone in eastern Nepal and the southern Tibetan Plateau

Focal mechanisms determined from moment tensor inversion and first motion polarities of the Himalayan Nepal Tibet Seismic Experiment (HIMNT) coupled with previously published solutions show the Himalayan continental collision zone near eastern Nepal is deforming by a variety of styles of deformation. These styles include strike-slip, thrust and normal faulting in the upper and lower crust, but mostly strike-slip faulting near or below the crust–mantle boundary (Moho). One normal faulting earthquake from this experiment accommodates east–west extension beneath the Main Himalayan Thrust of the Lesser Himalaya while three upper crustal normal events on the southern Tibetan Plateau are consistent with east–west extension of the Tibetan crust. Strike-slip earthquakes near the Himalayan Moho at depths >60 km also absorb this continental collision. Shallow plunging P -axes and shallow plunging EW trending T -axes, proxies for the predominant strain orientations, show active shearing at focal depths ∼60–90 km beneath the High Himalaya and southern Tibetan Plateau. Beneath the southern Tibetan Plateau the plunge of the P -axes shift from vertical in the upper crust to mostly horizontal near the crust–mantle boundary, indicating that body forces may play larger role at shallower depths than at deeper depths where plate boundary forces may dominate.     

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

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

P-Wave Velocities in the Upper Mantle Beneath Fennoscandia and Western Russia

A detailed and extensive record section constructed from recordings at the NORSAR array of presumed explosions in continental Russia exhibits two distinct ( T , Δ) triplications. The reliable identification of these upper mantle travel-time branches is possible because of the dense areal sampling of the NORSAR configuration. A simple upper mantle P- velocity model which can account adequately for the data involves velocity discontinuities at depths of 420 km and 690 km, and fairly uniform velocity gradients elsewhere. For this model, the first arrival branch for Δ≤ 21° extends as a second arrival to a distance of about 33°, at which distance it is terminated by the 420-km discontinuity. Rays bottoming between depths of 420 and 690 km span the distance range 16°≤Δ≤ 28°, and give first arrivals in the range 21°≤, Δ, 24°. Rays which penetrate the 690-km discontinuity give rise to secondary arrivals in the range 19°≤Δ≤ 25°, and first arrivals for distances Δ≤ 25°.     

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 .     

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

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

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.     

Radial variation of compressional and shear velocities in the Earth's lower mantle

Summary. This paper extends an earlier study (Sengupta & Julian) of travel times of P waves of deep-focus earthquakes to include shear waves. Primary advantage of deep-focus earthquakes is the reduction of anomalies caused by complex structures near the source. The standard deviations of travel times and station anomalies of this study are about half as large as those determined from the data of shallow-focus earthquakes (e.g. Herrin et al.; Hales & Roberts). Spherically-symmetric velocity models derived from the travel times by a linearized inverse technique have resolving lengths of about 70 km for standard errors in velocity of about 0.02 km/s. No pronounced reversal of either compressional or shear velocity was required at the base of the mantle to satisfy the data, though a small velocity decrease could not be entirely ruled out. Some anomalous rapid changes in compressional velocity gradient were, however, found centred around the depths of 2400 and 2600 km. The models derived in this study agree most closely with that of Herrin et al . for compressional velocity and the model 1066B of Gilbert & Dziewonski for shear velocity.     

Upper mantle shear structure of North America

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

Upper mantle heterogeneities beneath Central Europe

Summary. An inversion of ISC travel-time data from selected earthquakes in the distance range 30°-90° to 53 stations in Central Europe has been used to model velocity down to 600 km depth. The model explains 0.1–0.2s of the residuals, as for other array studies, leaving 0.5 s unexplained as noise. The uppermost 100 km of the mantle and crust contains inhomogeneities that correlate remarkably well with the geology. This may be due to deep-seated thermal anomalies or, in some areas, to delays introduced by passage of the rays through sedimentary cover. The deeper anomalies are smaller and unrelated to those in the lithosphere, which suggests that the asthenosphere is decoupled from the rigid lithosphere. The structure at 600 km depth is again quite inhomogeneous and might be due to undulations of the 650 km discontinuity. The models show some suggestion of a high velocity slab trending from east to west beneath the Alps.     

The deep seismicity of the Tyrrhenian Sea

The study reappraises the deep seismicity of the Tyrrhenian Sea. Careful examination of the quality of reported hypocentres shows that the earthquakes define a zone dipping NW, about 200 km along strike, 50 km thick, and reaching a depth of about 500 km. The zone is slightly concave to the NW at a depth of 300 km, but, contrary to many previous reports, is not tightly concave, nor are there significant spatial gaps in the seismicity, which is effectively continuous with depth. Seismicity is, however, concentrated in the depth interval 250–300 km, where the dip of the seismic zone changes from 70° (above 250 km) to a more gentle dip of 45° at greater depths. Seven fault-plane solutions are available for the largest earthquakes in this depth interval, all of them consistent with a P -axis down the dip of the seismic zone, and all of them requiring movement on faults out of the plane of the subducting slab.
Two deep earthquakes near Naples lie well outside the main zone of activity; for one of which a fault-plane solution is available that has a P -axis not aligned with the dip of the seismic zone. The tightly concave slab-geometry favoured by other reports is supported mainly by the location of these events near Naples, which we think may represent deformation in a separate, probably shallower dipping, piece of subducted lithosphere.
The lack of shallow seismicity, and particularly of thrust faulting earthquakes, at the surface projection of the Benioff zone suggests that active subduction has ceased. Estimates of the convergence rate responsible for subduction in the last 10 Myr far exceed the present convergence rate of Africa and Eurasia, suggesting that the subduction was related instead to the stretching and thinning of the crust in the Tyrrhenian Sea.     

Geometry and brittle deformation of the subducting Nazca Plate, Central Chile and Argentina

We use data from the Chile Argentina Geophysical Experiment (CHARGE) broad-band seismic deployment to refine past observations of the geometry and deformation within the subducting slab in the South American subduction zone between 30°S and 36°S. This region contains a zone of flat slab subduction where the subducting Nazca Plate flattens at a depth of ∼100 km and extends ∼300 km eastward before continuing its descent into the mantle. We use a grid-search multiple-event earthquake relocation technique to relocate 1098 events within the subducting slab and generate contours of the Wadati-Benioff zone. These contours reflect slab geometries from previous studies of intermediate-depth seismicity in this region with some small but important deviations. Our hypocentres indicate that the shallowest portion of the flat slab is associated with the inferred location of the subducting Juan Fernández Ridge at 31°S and that the slab deepens both to the south and the north of this region. We have also determined first motion focal mechanisms for ∼180 of the slab earthquakes. The subhorizontal T -axis solutions for these events are almost entirely consistent with a slab pull interpretation, especially when compared to our newly inferred slab geometry. Deviations of T -axes from the direction of slab dip may be explained with a gap within the subducting slab below 150 km in the vicinity of the transition from flat to normal subducting geometry around 33°S.     

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.     

Seismic structure of the upper mantle beneath the United States by three-dimensional inversion of body wave arrival times

Summary. Teleseismic P and S arrival times to North American stations are obtained from the ISC bulletins for the 10-yr period 1964–73, and relative travel-time delays are calculated with respect to standard tables. Station anomalies as well as variations of the delays with azimuth and epicentral distance from station are analysed, and the location of the velocity anomalies responsible for them is discussed. Inversion of the P delays to infer upper mantle velocity structure down to a depth of 700 km is obtained using three-dimensional blocks, as proposed by Aki, Christofferson & Husebye. Three layers can be resolved in this depth range. It is found that the heterogeneities responsible for the travel-time delays are primarily located in the first 250 km of the upper mantle, and that they correlate with surface features. Significant heterogeneities subsist to depths of at least 700 km and their broad scale pattern also correlates with the surface features: in the third layer (500 to 700 km depth) there is an increase of velocity from the West to the East of the United States, while the second layer (250 to 450 km depth) exhibits a reversed pattern. A tentative interpretation of these deeper anomalies is made, as being due mainly to topography of the major upper mantle discontinuities, near 400 and 650 km depth.     

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.     

State of stress in the Southern Tyrrhenian subduction zone from fault-plane solutions

In this paper we present revised locations and original focal mechanisms computed for intermediate and deep earthquakes that occurred within the Southern Tyrrhenian subduction zone between 1988 and 1994, in order to improve our knowledge of the state of stress for this compressional margin. In particular, we define the stress distribution within a large portion of the descending slab, between 40 and about 450 km depth. The seismicity distribution reveals a continuous 40–50 km thick slab that abruptly increases its dip from subhorizontal in the Ionian Sea to a constant 70° dip in the Tyrrhenian. We computed focal mechanisms for events with magnitudes ranging from 2.7 and 5.7, obtaining the distribution of P - and T -axes for many events for which centroid moment tensor (CMT) solutions are not available, thus enabling the sampling of a larger depth range compared to previous studies. We define three portions of the slab characterized by different distributions of P - and T -axes. A general down-dip compression is found between 165 and 370 km depth, whereas in the upper part of the slab (40–165 km depth) the fault-plane solutions are strongly heterogeneous. Below 370 km the P -axes of the few deep events located further to the north have a shallower dip and are not aligned with the 70° dipping slab, possibly suggesting that they belong to a separated piece of subducted lithosphere. There is a good correspondence between the depth range in which the P -axes plunge closer to the slab dip (∼ 70°) and the interval characterized by the highest seismic energy release (190–370 km).     

Constraints on a plume in the mid-mantle beneath the Iceland region from seismic array data

Teleseismic P waves passing through low-wave-speed bodies in the mantle are refracted, causing anomalies in their propagation directions that can be measured by seismometer arrays. Waves from earthquakes in the eastern Pacific and western North America arriving at the NORSAR array in Norway and at seismic stations in Scotland pass beneath the Iceland region at depths of ∼ 1000–2000 km. Waves arriving at NORSAR have anomalous arrival azimuths consistent with a low-wave-speed body at a depth of ∼ 1500 km beneath the Iceland–Faeroe ridge with a maximum diameter of ∼250 km and a maximum wave-speed contrast of ∼ 1.5 per cent. This agrees well with whole-mantle tomography results, which image a low-wave-speed body at this location with a diameter of ∼ 500 km and a wave-speed anomaly of ∼ 0.5 per cent, bearing in mind that whole-mantle tomography, because of its limited resolution, broadens and weakens small anomalies. The observations cannot resolve the location of the body, and the anomaly could be caused in whole or in part by larger bodies farther away, for example by a body imaged beneath Greenland by whole-mantle tomography.     

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.     

Ivrea seismic array: a study of continental crust and upper mantle

A seismic-array study of the continental crust and upper mantle in the Ivrea-Yerbano and Strona-Ceneri zones (northwestern Italy) is presented. A short-period network is used to define crustal P - and S -wave velocity models from earthquakes. The analysis of the seismic-refraction profile LOND of the CROP-ECORS project provided independent information and control on the array-data interpretation.
Apparent-velocity measurements from both local and regional earthquakes, and time-term analysis are used to estimate the velocity in the lower crust and in the upper mantle. The geometry of the upper-lower crust and Moho boundaries is determined from the station delay times.
We have obtained a three-layer crustal seismic model. The P -wave velocity in the upper crust, lower crust and upper mantle is 6.1±0.2 km s⁻¹, 6.5±0.3 km s⁻¹ and 7.8±0.3 km s⁻¹ respectively. Pronounced low-velocity zones in the upper and lower crust are not observed. A clear change in the velocity structure between the upper and lower crust is documented, constraining the petrological interpretation of the Ivrea-type reflective lower continental crust derived from small-scale petrophysical data. Moreover, we found a V _P/ V _S ratio of 1.69±0.04 for the upper crust and 1.82±0.08 for the lower crust and upper mantle. This is consistent with the structural and petrophysical differences between a compositionally uniform and seismically transparent upper crust and a layered and reflective lower crust. The thickness of the lower crust ranges from about 8 km in front of the Ivrea body (ARVO, Arvonio station) in the northern part of the array to a maximum of about 15 km in the southern part of the array. The lower crust reaches a minimum depth of 5 km below the PROV (Provola) station.     

Correlation between seismic anisotropy and Bouguer gravity anomalies in Tibet and its implications for lithospheric structures

As a baseline measurement for understanding the Himalayan–Tibetan orogen, a product of continent–continent collision between India and Eurasia, we analyse digital seismic data in order to constrain the seismic anisotropy of the Indian shield. Based on spatially sparse data that are currently available in the public domain, there is little shear-wave birefringence for SKS phases under the Indian shield, even though it is part of a fast-moving plate in the hotspot frame of reference. If most of the northern Indian mantle has little transverse anisotropy, the onset of significant anisotropy under Tibet marks the northern terminus of intact Indian lithosphere that is thrusting under the Himalayan–Tibetan orogen. Beyond this terminus, tectonic fabric such as that associated with the deforming lithospheric mantle of Eurasia must be present in the upper mantle. Along the profile from Yadong to Golmud, the only profile in Tibet where a number of shear-wave birefringence data are available, the amount of birefringence shows two marked increases, near 30° and 33°N, between which a local high in Bouguer gravity anomaly is observed. Such a correlation between patterns of shear-wave birefringence and gravity anomalies is explained by the juxtaposition of Indian lithosphere against the overlying Eurasian lithosphere: while the Eurasian lithospheric mantle appears only to the north of 30°N, the Indian lithospheric mantle extends northwards to near 33°N.