Seismic tomographic imaging of P- and S-waves velocity perturbations in the upper mantle beneath Iran

The inverse tomography method has been used to study the P - and S -waves velocity structure of the crust and upper mantle underneath Iran. The method, based on the principle of source–receiver reciprocity, allows for tomographic studies of regions with sparse distribution of seismic stations if the region has sufficient seismicity. The arrival times of body waves from earthquakes in the study area as reported in the ISC catalogue (1964–1996) at all available epicentral distances are used for calculation of residual arrival times. Prior to inversion we have relocated hypocentres based on a 1-D spherical earth's model taking into account variable crustal thickness and surface topography. During the inversion seismic sources are further relocated simultaneously with the calculation of velocity perturbations. With a series of synthetic tests we demonstrate the power of the algorithm and the data to reconstruct introduced anomalies using the ray paths of the real data set and taking into account the measurement errors and outliers. The velocity anomalies show that the crust and upper mantle beneath the Iranian Plateau comprises a low velocity domain between the Arabian Plate and the Caspian Block. This is in agreement with global tomographic models, and also tectonic models, in which active Iranian plateau is trapped between the stable Turan plate in the north and the Arabian shield in the south. Our results show clear evidence of the mainly aseismic subduction of the oceanic crust of the Oman Sea underneath the Iranian Plateau. However, along the Zagros suture zone, the subduction pattern is more complex than at Makran where the collision of the two plates is highly seismic.     

Insights into the lithospheric structure and tectonic setting of the Barents Sea region from isostatic considerations

We study the tectonic setting and lithospheric structure of the greater Barents Sea region by investigating its isostatic state and its gravity field. 3-D forward density modelling utilizing available information from seismic data and boreholes shows an apparent shift between the level of observed and modelled gravity anomalies. This difference cannot be solely explained by changes in crustal density. Furthermore, isostatic calculations show that the present crustal thickness of 35–37 km in the Eastern Barents Sea is greater than required to isostatically balance the deep basins of the area (>19 km). To isostatically compensate the missing masses from the thick crust and deep basins and to adequately explain the gravity field, high-density material (3300–3350 kg m⁻³) in the lithospheric mantle below the Eastern Barents Sea is needed. The distribution of mantle densities shows a regional division between the Western and Eastern Barents and Kara Seas. In addition, a band of high-densities is observed in the lower crust along the transition zone from the Eastern to Western Barents Sea. The distribution of high-density material in the crust and mantle suggests a connection to the Neoproterozoic Timanide orogen and argues against the presence of a Caledonian suture in the Eastern Barents Sea. Furthermore, the results indicate that the basins of the Western Barents Sea are mainly affected by rifting, while the Eastern Barents Sea basins are located on a stable continental platform.     

Regional three-dimensional gravity modelling of the NE Atlantic margin

A series of three‐dimensional models has been constructed for the structure of the crust and upper mantle over a large region spanning the NE Atlantic passive margin. These incorporate isostatic and flexural principles, together with gravity modelling and integration with seismic interpretations. An initial isostatic model was based on known bathymetric/topographic variations, an estimate of the thickness and density of the sedimentary cover, and upper mantle densities based on thermal modelling. The thickness of the crystalline crust in this model was adjusted to equalise the load at a compensation depth lying below the zone of lateral mantle density variations. Flexural backstripping was used to derive alternative models which tested the effect of varying the strength of the lithosphere during sediment loading. The models were analysed by comparing calculated and observed gravity fields and by calibrating the predicted geometries against independent (primarily seismic) evidence. Further models were generated in which the thickness of the sedimentary layer and the crystalline crust were modified in order to improve the fit to observed gravity anomalies. The potential effects of igneous underplating and variable upper mantle depletion were explored by a series of sensitivity trials. The results provide a new regional lithospheric framework for the margin and a means of setting more detailed, local investigations in their regional context. The flexural modelling suggests lateral variations in the strength of the lithosphere, with much of the margin being relatively weak but areas such as the Porcupine Basin and parts of the Rockall Basin having greater strength. Observed differences between the model Moho and seismic Moho along the continental margin can be interpreted in terms of underplating. A Moho discrepancy to the northwest of Scotland is ascribed to uplift caused by a region of upper mantle with anomalously low density, which may be associated with depletion or with a temperature anomaly.     

Crustal and upper-mantle structure in the Eastern Mediterranean from the analysis of surface wave dispersion curves

The dispersive properties of surface waves are used to infer earth structure in the Eastern Mediterranean region. Using group velocity maps for Rayleigh and Love waves from 7 to 100 s, we invert for the best 1-D crust and upper-mantle structure at a regular series of points. Assembling the results produces a 3-D lithospheric model, along with corresponding maps of sediment and crustal thickness. A comparison of our results to other studies finds the uncertainties of the Moho estimates to be about 5 km. We find thick sediments beneath most of the Eastern Mediterranean basin, in the Hellenic subduction zone and the Cyprus arc. The Ionian Sea is more characteristic of oceanic crust than the rest of the Eastern Mediterranean region as demonstrated, in particular, by the crustal thickness. We also find significant crustal thinning in the Aegean Sea portion of the backarc, particularly towards the south. Notably slower S -wave velocities are found in the upper mantle, especially in the northern Red Sea and Dead Sea Rift, central Turkey, and along the subduction zone. The low velocities in the upper mantle that span from North Africa to Crete, in the Libyan Sea, might be an indication of serpentinized mantle from the subducting African lithosphere. We also find evidence of a strong reverse correlation between sediment and crustal thickness which, while previously demonstrated for extensional regions, also seems applicable for this convergence zone.     

A three-dimensional image of shallow subduction: crustal structure of the Raukumara Peninsula, New Zealand

Earthquake arrival time data from a 36-station deployment of portable seismographs on the Raukumara Peninsula have been used to determine the 3-D Vp and Vp/Vs structure of this region of shallow subduction. A series of inversions have been performed, starting with an inversion for 1-D structure, then 2-D, and finally 3-D. This procedure ensures a smooth regional model in places of low resolution. The subducted plate is imaged as a northwest-dipping feature, with Vp consistently greater than 8.5  km  s⁻¹ in the uppermost mantle of the plate. Structure in the overlying plate changes significantly along strike. In the northeast, there is an extensive low-velocity zone in the lower crust underlying the most rapidly rising part of the Raukumara Range. It is bounded on its arcward side by an upwarp of high velocity. A viable explanation for the low-velocity zone is that it represents an accumulation of underplated subducted sediment, while serpentinization of the uppermost mantle may be responsible for the adjacent high-velocity region. The low-velocity zone decreases and the adjacent high-velocity region is less extensive in the southwest. This change is interpreted to be related to a change in the thickness of the crust of the overlying plate. In the northeast the crust is thinner, and subducted sediment ponds against relatively strong uppermost mantle, while in the southwest the crust is thicker, and the relatively weak lower crust allows sediment subduction to greater depths. A narrow zone of high Vp/Vs parallels the shallow part of the plate interface. This suggests elevated fluid pressures, with the distribution of earthquakes about this zone further suggesting that these pressures may be close to lithostatic. The plate interface at 20  km depth beneath the Raukumara Peninsula may thus be a closed system for fluid flow, similar to that seen at much shallower depths in other subduction décollements.     

Small-scale mantle flows and induced lithospheric stress near island arcs

Summary. This paper explores the middle ground between complex thermally-coupled viscous flow models and simple corner flow models of island arc environments. The calculation retains the density-driven nature of convection and relaxes the geometrical constraints of corner flow, yet still provides semianalytical solutions for velocity and stress. A novel aspect of the procedure is its allowance for a coupled elastic lithosphere on top of a Newtonian viscous mantle. Initially, simple box-like density drivers illustrate how vertical and horizontal forces are transmitted through the mantle and how the lithosphere responds by trench formation. The flexural strength of the lithosphere spatially broadens the surface topography and gravity anomalies relative to the functional form of the vertical flow stresses applied to the plate base. I find that drivers in the form of inclined subducting slabs cannot induce self-driven parallel flow; however, the necessary flow can be provided by supplying a basal drag of 1–5 MPa to the mantle from the oceanic lithosphere. These basal drag forces create regional lithospheric stress and they should be quantifiable through seismic observations of the neutral surface. The existence of a shallow elevated phase transition is suggested in two slab models of 300 km length where a maximum excess density of 0.2 g cm⁻³ was needed to generate an acceptable mantle flow. A North New Hebrides subduction model which satisfies flow requirements and reproduces general features of topography and gravity contains a high shear stress zone (75 MPa) around the upper slab surface to a depth of 150 km and a deviatoric tensional stress in the back arc to a depth of 70 km. The lithospheric stress state of this model suggests that slab detachment is possible through whole plate fracture.     

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.     

Solution of the inverse problem of seismology for laterally inhomogeneous media

Summary. The Backus-Gilbert method has been extended to the estimation of the seismic wave velocity distribution in 2-D or 3-D inhomogeneous media from a finite set of travel-time data. The method may be applied to the inversion of body wave as well as surface wave data. The problem of determining a local average of the unknown velocity corrections may be reduced to a choice of a suitable δ-ness criterion for the averaging kernel. For 2-D and 3-D inhomogeneous media the simplest criterion is to minimize a sum of 'spreads' over all the coordinates. The use of this criterion requires the solution (the averaged velocity corrections) to be represented as a sum of functions, each of which depends only on one coordinate. This is a basic restriction of the method. In practice it is possible to achieve good agreement between the solution and a real velocity distribution by a reasonable choice of the coordinate system.
Numerical tests demonstrate the efficiency of the method. Some examples of the application of the method to the inversion of real seismological data for body and surface waves are given.     

Double-planed deep seismic zone and upper-mantle structure in the Northeastern Japan Arc

Summary. The ScSp wave converted from the ScS wave at the boundary between the descending lithospheric slab and the mantle above it was clearly observed from a nearby deep earthquake with magnitude 7.7 at some stations of the seismic network of Tohoku University which covers the Tohoku District, the northeastern part of Honshu, Japan. By applying the three-dimensional seismic-ray tracing method, the location of this boundary was determined from the difference in arrival time between the ScS and ScSp waves. The result shows that the upper boundary of the descending slab lies exactly on the upper plane of the double-planed deep seismic zone found in the Northeastern Japan Arc.
There is an additional evidence that the boundary is located on the upper plane of the double-planed deep seismic zone. The hypocentre distribution of intermediate-depth earthquakes located by the small-scale seismic-array observation is extremely different from that obtained by the relatively large-scale seismic network. The discrepancy in the distribution of hypocentres of the same earthquake independently located is well explained by the inclined lithospheric slab model derived from the difference in arrival time between the ScS and ScSp waves.
The earthquakes with reverse faulting or with down-dip compressional stresses occur at the upper boundary of the descending slab. Within the descending slab, the earthquakes with down-dip extensional stresses also occur in a very narrow zone from 30 to 40 km below the dipping boundary in the depth range from 50 to about 200 km, and these shocks form the lower plane of the double-planed deep seismic zone.     

Massif Central (France): new constraints on the geodynamical evolution from teleseismic tomography

The Massif Central, the most significant geomorphological unit of the Hercynian belt in France, is characterized by graben structures which are part of the European Cenozoic Rift System (ECRIS) and also by distinct volcanic episodes, the most recent dated at 20 Ma to 4000 years BP. In order to study the lithosphere-asthenosphere system beneath this volcanic area, we performed a teleseismic field experiment.
During a six-month period, a joint French-German team operated a network of 79 mobile short-period seismic stations in addition to the 14 permanent stations. Inversion of P -wave traveltime residuals of teleseismic events recorded by this dense array yielded a detailed image of the 3-D velocity structure beneath the Massif Central down to 180 km depth. The upper 60 km of the lithosphere displays strong lateral heterogeneities and shows a remarkable correlation between the volcanic provinces and the negative velocity perturbations. The 3-D model reveals two channels of low velocities, interpreted as the remaining thermal signature of magma ascent following large lithospheric fractures inherited from Hercynian time and reactivated during Oligocene times. The teleseismic inversion model yields no indication of a low-velocity zone in the mantle associated with the graben structures proper. The observation of smaller velocity perturbations and a change in the shape of the velocity pattern in the 60–100 km depth range indicates a smooth transition from the lithosphere to the asthenosphere, thus giving an idea of the lithosphere thickness. A broad volume of low velocities having a diameter of about 200 km from 100 km depth to the bottom of the model is present beneath the Massif Central. This body is likely to be the source responsible for the volcanism. It could be interpreted as the top of a plume-type structure which is now in its cooling phase.     

Three-dimensional image of the Moho undulations beneath the Gulf of Bothnia using wide-angle seismic data

The BABEL marine seismic experiment has been carried out to investigate the lithospheric structure and antecedent tectonic signatures of the Baltic Shield, including the Archaean-Proterozoic collisional structure in the northern part of the Gulf of Bothnia.
Lithospheric seismic-reflection streamer data and simultaneously recorded wide-angle reflection and refraction data collected in the Gulf of Bothnia as part of the BABEL project have been used for 3-D modelling. The distribution of land stations around the Gulf provides a good 3-D ray coverage of the PMP reflection data recorded at the eight stations in the area and allows an estimation of strikes and dips of the Moho boundary in the area. The traveltimes of reflected phases are calculated using a method that utilizes the finite-difference solution of the eikonal equation. The Moho wide-angle-reflection (PMP) traveltimes are modelled using an inversion method. A 2-D model from the Gulf of Bothnia extended into the third dimension is used as an initial model. During the inversion the velocity is kept constant and only the Moho boundary is allowed to vary. To estimate the strike of the Moho boundary and the stability of the inversion, two initial models with different strikes are examined.
The results indicate that the Moho depth in the Gulf of Bothnia undulates and has a maximum depth of 55 km in the south, rising to 42 km in the north. The Moho depth variations seem to be step-like. This change in the Moho depth coincides with the location of the presumed fossil subduction zone in the area. The crustal-thickness variations seem to be well approximated by a nearly 2-D structure striking parallel to a postulated subduction zone immediately to the south of the Skellefte area. The presence of the step at the crust/mantle boundary can be interpreted as a result of a plate-collision event at about 2 Ga.     

Axial Intrusion Zone beneath the Median Valley of the Mid-Atlantic Ridge at 37°N Detected by Explosion Seismology

Summary. Four seismic refraction lines, three of which had shots every 250 m, were shot across, along and parallel to the median valley of the Mid-Atlantic Ridge at 37° N. A method has been developed for calculating the effect on the travel times of the rough sea-floor relief beneath the profiles and has been used to correct all the travel times for this effect. Most arrivals were from a main refractor of apparent velocity 5·4 to 6·3 km s⁻¹; only beyond 35 km were faster arrivals observed from an 8·09 ± 36 km s⁻¹ refractor. The main refractor corresponds in depth, at least approximately, to the top of Layer 3 of the ocean basins but its velocity is significantly less than normal for Layer 3, perhaps due to dip. A study of time residuals along two profiles across the median valley indicates the presence of a 2 to 3 km wide low velocity zone (about 3·2 km s⁻¹) beneath the median valley floor. This zone extends over the upper 2·5 km of the crust and is believed to represent a zone of intrusion through which magma passes on its way to the sea floor.     

A joint hypocentre location and velocity inversion technique applied to the Central Aleutians

Summary. Travel times of seismic waves from teleseismically recorded events that occurred within the Adak Island region of the Central Aleutian Islands are used to constrain both the velocity structure of the region and the locations of events. It is found that the P-wave velocity within the descending lithosphere is approximately 7–8 per cent faster than that in the surrounding mantle for a slab length of 300–400 km.
The velocity structure so determined is then used as a model for the relocation of events recorded only locally. It is found that features of the Benioff zone seen in locations made with a spherically symmetric model can be seriously in error. In particular, an apparent increase in dip of the zone with depth disappears when the locations are made with the inhomogeneous model. The new locally determined hypocentres are now also more consistent with the locations based on teleseismic data.     

Upper Mantle Velocity Structure Estimated From Ps-Converted Wave Beneath the North-Eastern Japan Arc

Summary. The upper boundary of the descending oceanic plate is located by using PS -waves (converted from P to S at the boundary) in the Tohoku District, the north-eastern part of Honshu, Japan. the observed PS-P time data are well explained by a two-layered oceanic plate model composed of a thin low-velocity upper layer whose thickness is less than 10 km and a thick high-velocity lower layer; the upper and lower layers respectively have 6 per cent lower and 6 per cent higher velocity than the overriding mantle. the estimated location of the upper boundary is just above the upper seismic plane of the double-planed deep seismic zone. This result indicates that events in the upper seismic plane, at least in the depth range from 60 to 150 km, occur within the thin low-velocity layer on the surface of the oceanic plate.     

Two-dimensional modelling of subduction zone anisotropy with application to southwestern Japan

We present a series of 2-D numerical models of viscous flow in the mantle wedge induced by a subducting lithospheric plate. We use a kinematically defined slab geometry approximating the subduction of the Philippine Sea plate beneath Eurasia. Through finite element modelling we explore the effects of different rheological and thermal constraints (e.g. a low-viscosity region in the wedge corner, power law versus Newtonian rheology, the inclusion of thermal buoyancy forces and a temperature-dependent viscosity law) on the velocity and finite strain field in the mantle wedge. From the numerical flow models we construct models of anisotropy in the wedge by calculating the evolution of the finite strain ellipse and combining its geometry with appropriate elastic constants for effective transversely isotropic mantle material. We then predict shear wave splitting for stations located above the model domain using expressions derived from anisotropic perturbation theory, and compare the predictions to ∼500 previously published shear wave splitting measurements from seventeen stations of the broad-band F-net array located in southwestern Japan. Although the use of different model parameters can have a substantial effect on the character of the finite strain field, the effect on the average predicted splitting parameters is small. However, the variations with backazimuth and ray parameter of individual splitting intensity measurements at a given station for different models are often different, and rigorous analysis of details in the splitting patterns allows us to discriminate among different rheological models for flow in the mantle wedge. The splitting observed in southwestern Japan agrees well with the predictions of trench-perpendicular flow in the mantle wedge along with B-type olivine fabric dominating in a region from the wedge corner to about 125 km from the trench.     

The reflection seismic survey of project TIPTEQ—the inventory of the Chilean subduction zone at 38.2° S

We describe results of an active-source seismology experiment across the Chilean subduction zone at 38.2°S. The seismic sections clearly show the subducted Nazca plate with varying reflectivity. Below the coast the plate interface occurs at 25 km depth as the sharp lower boundary of a 2–5 km thick, highly reflective region, which we interpret as the subduction channel, that is, a zone of subducted material with a velocity gradient with respect to the upper and lower plate. Further downdip along the seismogenic coupling zone the reflectivity decreases in the area of the presumed 1960 Valdivia hypocentre. The plate interface itself can be traced further down to depths of 50–60 km below the Central Valley. We observe strong reflectivity at the plate interface as well as in the continental mantle wedge. The sections also show a segmented forearc crust in the overriding South American plate. Major features in the accretionary wedge, such as the Lanalhue fault zone, can be identified. At the eastern end of the profile a bright west-dipping reflector lies perpendicular to the plate interface and may be linked to the volcanic arc.     

LISPB – IV. Crustal structure of Northern Britain

Summary. This paper presents those results from the 1974 Lithospheric Seismic Profile in Britain (LISPB) which relate to the compressional velocity structure of the crust and uppermost mantle beneath Northern Britain. A combination of interpretation techniques suitable for modelling laterally inhomogeneous media, including two-dimensional ray-tracing and time-term analysis, has resulted in a detailed seismic cross-section across the Caledonian orogenic belt. The main features of this section are a possible horizontal discontinuity in the Pre-Caledonian'basement, a change in the relationship between the lower crust and the uppermost mantle from north to south and a considerable thickening of the crust beneath the Caledonian fold belt. These results place considerable constraints upon tectonic models for the evolution of the Caledonides in particular in their implication of differing crustal structures north and south of the Southern Uplands and their indication of the primary significance of the Southern Uplands Fault.     

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.     

Computation of high-frequency seismic wavefields in 3-D laterally inhomogeneous anisotropic media

Summary. An algorithm for the computation of travel times, ray amplitudes and ray synthetic seismograms in 3-D laterally inhomogeneous media composed of isotropic and anisotropic layers is described. All 21 independent elastic parameters may vary within the anisotropic layers. Rays and travel times are evaluated by numerical solution of the ray tracing equations. Ray amplitudes are determined by evaluating reflection/ transmission coefficients and the geometrical spreading along individual rays. The geometrical spreading is computed approximately by numerical measurement of the cross-sectional area of the ray tube formed by three neighbouring rays. A similar approximate procedure is used for the determination of the coefficients of the paraxial ray approximation. The ray paraxial approximation makes computation of synthetic seismograms on the surface of the model very efficient. Examples of ray synthetic seismograms computed with a program package based on the described algorithm are presented.     

Lithospheric structure of the Arabian and Eurasian collision zone in eastern Turkey from S-wave receiver functions

Crustal and upper-mantle seismic discontinuities beneath eastern Turkey are imaged using teleseismic S -to- P converted phases. Three crustal phases are observed: the Moho with depth ranging between 30 and 55 km, indicating variable tectonic regimes within this continental collision zone; an upper-crustal discontinuity at approximately 10 km depth; and various crustal low-velocity zones, possibly associated with recent Quaternary volcanism. Imaging of the upper mantle is complicated by the 3-D geometry of the region, in particular due to the Bitlis–Zagros suture zone. However, several upper-mantle S -to- P converted phase are identified as being the signature of the lithosphere–asthenosphere boundary (LAB). The inferred LAB for the Eastern Anatolian Accretionary Complex indicates that eastern Turkey has an anomalously thin (between ∼60 and 80 km) lithosphere which is consistent with an oceanic slab detachment model. The observed LAB phases for the Arabian shield and Iranian plateau indicate that lithospheric thickness for these stable regions is on the order of 100 to 125 km thick, which is typical of continental margins.