Relative Arrival Time Residuals Research Articles

Post-collisional lithospheric delamination in eastern Iran, revealed by non-linear teleseismic tomography and residual topography

The Eastern Iranian Mountain Ranges (EIR) emerged as a consequence of the Late Cretaceous collision between the Afghan and Lut blocks. However, the response of the uppermost mantle to this collision remains enigmatic. Additionally, although petrological evidence suggests that post-collisional delamination is possible, it has not been conclusively identified in prior regional seismic imagery. This observation leads us to further explore this possibility using a dense seismic network. To gain insight into the geodynamic implications for eastern Iran and address knowledge gaps, we extensively investigated the seismic structure of the uppermost mantle beneath the EIR using a dense seismic network of 34 temporary stations, complemented by data from nine additional local permanent stations. By meticulously analyzing 6589 relative arrival time residuals from teleseismic records with favorable signal-to-noise ratios, we applied a non-linear tomography method to map P-wave velocity perturbations in a relative sense. Our tomographic images unveiled distinct instances of rapid high-velocity anomalies beneath low-velocity regions in the shallow mantle, suggesting the potential occurrence of lithospheric dripping, followed by subsequent asthenospheric upwelling. This observation offers a plausible explanation for the observed post-collisional magmatism over the Lut Block. Furthermore, to maintain the approximately 1.5-km positive residual topography across the EIR, beyond the influence of crustal properties, additional support from the hot and buoyant asthenosphere becomes crucial, particularly in the absence of a substantial lithospheric mantle.

Upper Mantle Velocity Structure Beneath the Yarlung–Tsangpo Suture Revealed by Teleseismic P-Wave Tomography

We applied teleseismic tomography to investigate the 3D P-wave velocity (Vp) structure of the crust and upper mantle at depths of 50–400 km beneath the Yarlung–Tsangpo suture (YTS), by using 6164 P-wave relative travel-time residuals collected from 495 teleseismic events recorded at 20 three-component broadband seismograms. A modified multi-channel cross-correlation method was adopted to automatically calculate the relative arrival-time residuals of all teleseismic events, which significantly improved the efficiency and precision of the arrival-time data collection. Our results show that alternating low- and high-Vp anomalies are visible beneath the Himalayan and Lhasa blocks across the YTS, indicating that strong lateral heterogeneities exist beneath the study region. A significant high-Vp zone is visible beneath the southern edge of the Lhasa block at 50–100 km depths close to the YTS, which might indicate the rigid Tibetan lithosphere basement. There exists a prominent low-Vp zone beneath the Himalayan block to the south of the YTS extending to ~150 km depth, which might be associated with the fragmentation of the underthrusting Indian continental lithosphere (ICL) and induce localized upwelling of asthenospheric materials from the upper mantle. In addition, significant low-Vp anomalies were observed beneath the Yadong–Gulu rift and the Cona–Sangri rift extending to ~300 km depth, indicating that the tearing of the subducted ICL might provide pathways for the localized asthenospheric materials upwelling, which contributes to the widespread distribution of north–south trending rifts and geothermal activities in southern Tibet.

Upper mantle structure of the northeastern Arabian Platform from teleseismic body-wave tomography

A thick Infra-Cambrian to Phanerozoic sedimentary sequence in eastern Arabia conceals a poorly known Neoproterozoic basement that once separated the currently exposed Arabian Shield (western Arabia) from NW India in Gondwana. In this study, the Arabian Platform (eastern Arabia) is targeted by relatively high-resolution passive-source seismic tomography for the first time, which aims to illuminate crust and upper mantle structure in order to make inferences on their nature and origin. This is facilitated by a new broadband seismic network of 55 stations in the United Arab Emirates (UAE), which has recorded sufficient teleseismic body wave data to allow for the determination of three-dimensional seismic structure. Using a grid-based eikonal solver and a subspace inversion technique implemented in FMTOMO, relative arrival-time residuals from teleseismic earthquakes are mapped as 3-D P-wave velocity perturbations in the upper mantle beneath the seismic network. The resultant tomographic images reveal a marked transition from higher velocities in the Gulf of Oman to lower velocities in the west that we interpret to indicate the lithospheric boundary between Cretaceous Tethyan oceanic lithosphere and the Arabian passive continental margin. The northern extent of Tethyan oceanic lithosphere appears to terminate against a possible offshore continuation of the Dibba fault, until now only identified onshore. Further inland, the tomographic model exhibits a number of low- and high- velocity anomalies with a predominant NW orientation. Significantly, one well-defined high velocity anomaly is situated along strike of the Neoproterozoic exposure of volcanic-arc granodiorites found in northern Oman. We relate this anomaly to an intra-oceanic arc created during the Tonian subduction of the Mozambique Ocean, which strongly support the idea of continental growth in the Arabian Platform via magmatic arc accretion.

Precambrian Plate Tectonics in Northern Hudson Bay: Evidence From P and S Wave Seismic Tomography and Analysis of Source Side Effects in Relative Arrival‐Time Data Sets

AbstractThe geology of northern Hudson Bay, Canada, documents more than 2 billion years of history including the assembly of Precambrian and Archean terranes during several Paleoproterozoic orogenies, culminating in the Trans‐Hudson Orogen (THO) ∼1.8 Ga. The THO has been hypothesized to be similar in scale and nature to the ongoing Himalaya‐Karakoram‐Tibetan orogen, but the nature of lithospheric terrane boundaries, including potential plate‐scale underthrusting, is poorly understood. To address this problem, we present new P and S wave tomographic models of the mantle seismic structure using data from recent seismograph networks stretching from northern Ontario to Nunavut (60–100∘W and 50–80∘N). The large size of our network requires careful mitigation of the influence of source side structure that contaminates our relative arrival time residuals. Our tomographic models reveal a complicated internal structure in the Archean Churchill plate. However, no seismic wave speed distinction is observed across the Snowbird Tectonic Zone, which bisects the Churchill. The mantle lithosphere in the central region of Hudson Bay is distinct from the THO, indicating potential boundaries of microcontinents and lithospheric blocks between the principal colliders. Slow wave speeds underlie southern Baffin Island, the leading edge of the generally high wave speed Churchill plate. This is interpreted to be Paleoproterozoic material underthrust beneath Baffin Island in a modern‐style subduction zone setting.

P- and S-wave tomography of the Hainan and surrounding regions: Insight into the Hainan plume

High-resolution 3-D P- and S-wave velocity structures beneath Hainan and surrounding regions down to a depth of 700km are determined using teleseismic tomographic inversion of 10,873 P and 9147 S relative arrival time residuals from 165 teleseismic events recorded by 88 permanent seismic stations belonging to three regional seismic networks of Hainan, Guangdong and Guangxi. The present tomographic models show that a low-velocity (low-V) anomaly extends below the transition zone (TZ) in and around the Hainan hotspot. The low-V zone does not exhibit a simple vertical columnar shape but a complex and deflected image toward the northeast. At depths of 13 to 170km, the main low-V anomaly exists right beneath the Hainan hotspot, but in the depth range of 170 to 700km, the low-V anomaly gradually deflects to the coastal area of Guangdong, which is to the northeast of the Hainan hotspot, an area with a diameter of approximately 200km. I think that this area may be the extent of the Hainan plume from the deep upper mantle to the bottom of the TZ. On the basis of the present tomographic models, combined with the receiver function and geochemical results, I believe that the Hainan plume may originate from the low mantle, and is characterized by a high-temperature anomaly.

The P and S wave velocity structure of the mantle beneath eastern Africa and the African superplume anomaly

P and S relative arrival time residuals from teleseismic earthquakes recorded on over 60 temporary AfricaArray broadband seismic stations deployed in Uganda, Tanzania, and Zambia between 2007 and 2011 have been inverted, together with relative arrival time residuals from earthquakes recorded by previous deployments, for a tomographic image of mantle wave speed variations extending to a depth of 1200 km beneath eastern Africa. The image shows a low‐wave speed anomaly (LWA) well developed at shallow depths (100–200 km) beneath the Eastern and Western branches of the Cenozoic East African rift system and northwestern Zambia, and a fast wave speed anomaly at depths ≤ 350 km beneath the central and northern parts of the East African Plateau and the eastern and central parts of Zambia. At depths ≥350 km the LWA is most prominent under the central and southern parts of the East African Plateau and dips to the southwest beneath northern Zambia, extending to a depth of at least 900 km. The amplitude of the LWA is consistent with a ∼150–300 K thermal perturbation, and its depth extent indicates that the African superplume, originally identified as a lower mantle anomaly, is likely a whole mantle structure. A superplume extending from the core‐mantle boundary to the surface implies an origin for the Cenozoic extension, volcanism, and plateau uplift in eastern Africa rooted in the dynamics of the lower mantle.

Transportable seismic array tomography in southeast Australia: Illuminating the transition from Proterozoic to Phanerozoic lithosphere

The Phanerozoic Tasmanides of eastern Australia is comprised of a series of orogenic belts that developed along the east margin of Gondwana following the breakup of the supercontinent Rodinia and subsequent formation of the Pacific Ocean. The tectonic complexities of this region have been well studied, but most work has been confined to evidence collected from the near surface, where extensive Mesozoic and Cenozoic basin cover masks large tracts of Palaeozoic basement. We apply teleseismic tomography to distant earthquake data recorded by WOMBAT – the largest transportable seismic array experiment in the southern hemisphere – to image P-wavespeed variations in the mantle lithosphere beneath the southern portion of the Tasmanides in detail. In order to seamlessly suture together the teleseismic datasets from each of the 14 sub-arrays of WOMBAT, we use P-wavespeeds from the AuSREM model to construct a laterally heterogeneous starting model that captures the long wavelength structural variations that would otherwise be lost through the use of relative arrival time residuals. Synthetic resolution tests indicate good horizontal resolution of ~50km within most of the array between depths of 50–350km. A key feature of the 3-D P-wave model is a pronounced easterly high velocity salient in the mantle lithosphere beneath the northern limit of the New England Orocline, which may indicate the presence of underpinning Proterozoic lithosphere that was instrumental in its formation. Another pronounced high velocity anomaly underlies the Curnamona Province, a large crustal block with a strong Archean provenance, which is clearly separated from the Gawler Craton to the west at upper mantle depths by a low velocity zone beneath the Adelaide Fold Belt. We also estimate the location of the eastern boundary of Precambrian Australia at depth, and show that it extends eastward further than previously thought.

Lithospheric structure beneath the Zagros collision zone resolved by non-linear teleseismic tomography

The upper-mantle structure across the Zagros collision zone, in southwest Iran, is investigated using a non-linear weighted damped least-squares teleseismic tomography approach. The resolution of the structures/transitions in the upper mantle is enhanced significantly by correcting the teleseismic relative arrival time residuals for an a priori crustal velocity model and then performing the inversion with fixed crustal blocks. To investigate whether or not the lithospheric blocks and major transitions in the resulting model are required by the data or are artefacts of the inversion, the data were inverted using two different inverse methods (singular value decomposition and a quadratic programming method). New high-quality seismic velocity models show apparent correlation between surface geological features and seismic velocity structures at lithospheric depth across the Zagros collision zone. The image shows a sharp lithospheric boundary at the Main Zagros Thrust between 100 km and 250 km depth with P-wave velocity about 3 per cent faster within the Arabian Shield to the south. A step-like increase in lithospheric thickness across the Zagros collision zone is assumed to separate two different mantle structures namely the Arabian (to the south) and the Eurasian (to the north) domains. The most striking feature resolved is a north-dipping slab-like positive velocity anomaly.

Teleseismic tomography of the mantle in the Carpathian-Pannonian region of central Europe

SUMMARY Subducted slab roll-back, lithospheric instability and asthenospheric extrusion have all been proposed as mechanisms that explain the evolution of the extensional Pannonian Basin, within the convergent arc of the Alpine–Carpathian mountain system in central Europe. We determine the P- and S-wave velocity structure of the mantle to depths of 850 km beneath this region using tomographic inversion of relative arrival-time residuals from 225 (P waves) and 124 (S waves) teleseismic earthquakes recorded by 56 stations of the Carpathian Basins Project (CBP) temporary seismic network (16-month duration) and 44 permanent seismic stations. The observed median P-wave relative arrival-time residuals vary between −1.13 s (early) in the Alps and 1.12 s (late) at the western end of the Carpathians; S-wave relative arrival-time residuals are about twice as large (−2.13 s and 3.39 s). We tested the effect of deterministic corrections on our relative arrival-time residuals using crustal velocity models from controlled source experiments, but show that the use of station terms in the inversion provides a robust method of correcting for near-surface crustal variation. Our tomographic models reduce the P-wave rms residual by 71 per cent to 0.130 s and our S-wave rms residual by 59 per cent to 0.624 s. At shallow sublithospheric depths we image several localized lower velocity regions, correlated with higher heat flow and interpreted as upwelling asthenosphere. We image a high velocity structure down to depths of about 350 km beneath the Eastern Alps. Further east, beneath the Pannonian Basin, a deeper continuation of the Eastern Alps fast anomaly is imaged trending E–W from ∼300 km depth and extending into the mantle transition zone (MTZ). In the MTZ we image a fast anomaly extending outwards as far as the Carpathians, the Dinarides and the Eastern Alps. This higher velocity mantle material is interpreted as being produced by a mantle downwelling, whose detachment from the lithosphere above may have triggered the extension of the Pannonian Basin.

Lithosphere-asthenosphere interaction beneath Ireland from joint inversion of teleseismic P-wave delay times and GRACE gravity

The nature and extent of the regional lithosphere-asthenosphere interaction beneath Ireland and Britain remains unclear. Although it has been established that ancient Caledonian signatures pervade the lithosphere, tertiary structure related to the Iceland plume has been inferred to dominate the asthenosphere. To address this apparent contradiction in the literature, we image the 3-D lithospheric and deeper upper-mantle structure beneath Ireland via non-linear, iterative joint teleseismic-gravity inversion using data from the ISLE (Irish Seismic Lithospheric Experiment), ISUME (Irish Seismic Upper Mantle Experiment) and GRACE (Gravity Recovery and Climate Experiment) experiments. The inversion combines teleseismic relative arrival time residuals with the GRACE long wavelength satellite derived gravity anomaly by assuming a depth-dependent quasilinear velocity-density relationship. We argue that anomalies imaged at lithospheric depths probably reflect compositional contrasts, either due to terrane accretion associated with Iapetus Ocean closure, frozen decompressional melt that was generated by plate stretching during the opening of the north Atlantic Ocean, frozen Iceland plume related magmatic intrusions, or a combination thereof. The continuation of the anomalous structure across the lithosphere-asthenosphere boundary is interpreted as possibly reflecting sub-lithospheric small-scale convection initiated by the lithospheric compositional contrasts. Our hypothesis thus reconciles the disparity which exists between lithospheric and asthenospheric structure beneath this region of the north Atlantic rifted margin.

Upper-mantle seismic structure beneath SE and Central Brazil from P- and S-wave regional traveltime tomography

SUMMARY We present models for the upper-mantle velocity structure beneath SE and Central Brazil using independent tomographic inversions of P- and S-wave relative arrival-time residuals (including core phases) from teleseismic earthquakes. The events were recorded by a total of 92 stations deployed through different projects, institutions and time periods during the years 1992–2004. Our results show correlations with the main tectonic structures and reveal new anomalies not yet observed in previous works. All interpretations are based on robust anomalies, which appear in the different inversions for P- and S-waves. The resolution is variable through our study volume and has been analyzed through different theoretical test inversions. High-velocity anomalies are observed in the western portion of the Sao Francisco Craton, supporting the hypothesis that this Craton was part of a major Neoproterozoic plate (San Franciscan Plate). Low-velocity anomalies beneath the Tocantins Province (mainly fold belts between the Amazon and Sao Francisco Cratons) are interpreted as due to lithospheric thinning, which is consistent with the good correlation between intraplate seismicity and low-velocity anomalies in this region. Our results show that the basement of the Parana Basin is formed by several blocks, separated by suture zones, according to model of Milani & Ramos. The slab of the Nazca Plate can be observed as a high-velocity anomaly beneath the Parana Basin, between the depths of 700 and 1200 km. Further, we confirm the low-velocity anomaly in the NE area of the Parana Basin which has been interpreted by VanDecar et al. as a fossil conduct of the Tristan da Cunha Plume related to the Parana flood basalt eruptions during the opening of the South Atlantic.

S and P velocity heterogeneities within the upper mantle below the Baltic Shield

Upper mantle structure beneath the Baltic (Fennoscandian) Shield is investigated using non-linear tomographic inversion of relative arrival-time residuals. 52 selected teleseismic earthquakes recorded by 45 broadband stations of the Swedish National Seismological Network (SNSN) provide 1532 good quality S-wave relative arrival times. SV and SH arrival-time residuals were initially analyzed independently, providing two separate models. These reveal several consistent major features, many of which are also consistent with P-wave results. Lateral velocity variations of ± 3–4% are observed to depths of at least 470 km. The correlation between the SH and SV models is investigated and shows a pattern of minor but significant differences down to around 150–200 km depth, below which the models are essentially similar. Direct cell by cell comparison of the model velocities reveals a similar pattern, with velocity differences between the models of up to 4%. Numerical tests show that differences in the two S-wave models can only be partially attributed to noise and limited resolution, and some features are attributed to the effect of large scale anisotropy. One of the significant and sharp discrepancies between the S models coincides with a presumed boundary between Archean and Proterozic domains, suggesting different anisotropic characteristics in the two regions.

Teleseismic tomography of the upper mantle beneath the southern Lachlan Orogen, Australia

With the goal of better understanding the deep structure and tectonic setting of the Lachlan Orogen, 50 short period seismic stations were deployed across the southern end of the Great Dividing Range in Victoria (southeast Australia) between 2005 and 2006 to record distant earthquakes. A total of 7452 relative P-wave arrival time residuals from 169 teleseismic events have been extracted from the continuous records using an adaptive stacking technique, which exploits the coherency of global phases across the array. These residuals are mapped as three-dimensional perturbations in P-wavespeed in the upper mantle beneath the array using a recently developed iterative non-linear tomographic procedure, which combines a grid based eikonal solver and a subspace inversion technique. The capability of the new scheme to include interface geometry is utilised in order to investigate the effects of a priori Moho topography on the resolution of upper mantle structure. The resultant images show a pattern of P-wavespeed anomalies that lacks a predominant orientation, and therefore does not favour a purely W–E subduction–accretion model for the formation of the Lachlan Orogen. One of the main features of the three-dimensional model is a zone of elevated wavespeed beneath the northern end of the array, which extends to a depth of approximately 150 km, and contrasts with significantly lower wavespeeds to the south. This anomaly, which does not appear to be an artifact of relative arrival time residual contributions from the adjoining mountainous terrane, may reflect the presence of a substantial piece of Proterozoic lithosphere incorporated within the Phanerozoic subduction–accretion setting of the Lachlan Orogen. Another key feature of the solution model is a zone of relatively low velocity beneath the newer volcanic province northwest of Melbourne, which extends from the crust to a depth of approximately 200 km. This is likely to represent the signature of elevated temperatures associated with a diffuse mantle source for the Quaternary volcanism in Victoria.

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

SUMMARY 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.

Insights into the structure of the upper mantle beneath the Murray basin from 3D teleseismic tomography

Distant earthquake records from a passive array of 20 short-period seismometers are used to image the 3D P-wave velocity structure of the upper mantle beneath the Murray Basin in southern New South Wales and northern Victoria. During the five-month deployment period of the array, 158 teleseisms with good signal-to-noise ratios were recorded, allowing 3085 relative arrival time residuals to be picked with high accuracy. These arrival time residuals are mapped as 3D perturbations in P-wave velocity with respect to the ak135 global reference model using teleseismic tomography. The resulting images exhibit lateral variations in wave speed in the upper mantle to depths of nearly 300 km, and help to reveal the deep structure beneath a region almost devoid of Palaeozoic outcrop. Of particular significance is an east – west high – low – high relative wave speed variation that is present throughout the upper mantle between 70 and 250 km depth. The transition from faster to slower wave speeds in the west is indicative of a change from Proterozoic to Phanerozoic lithosphere that has also been observed further south in a previous teleseismic study. The transition from slower to faster wave speeds in the east has a less obvious association, but may well signal an underlying change in lithospheric structure related to the convergence of the western and central Lachlan subprovinces. At shallower depths (<70 km), the structural pattern is dominated by a fast perturbation (>3.5%) which overlies the slower region across the northern margin of the Bendigo Zone, and may be a consequence of localised thrusting and emplacement of denser lithosphere during the Early – Middle Devonian Tabberabberan Orogeny.

Lithospheric structure of the Tornquist Zone resolved by nonlinear P and S teleseismic tomography along the TOR array

The main aim of the TOR project is to study the lithospheric–asthenospheric boundary structure under the Sorgenfrei–Tornquist Zone, across northern Germany, Denmark and southern Sweden. Relative arrival-time residuals of teleseismic P and S phases from 51 earthquakes, recorded by 150 seismic stations along the TOR array, were used to delineate the transition zone in the studied area. The effects of crustal structures were investigated by correcting the teleseismic residuals for travel-time variations in the crust based on a 3D crustal model derived from other data. The inversion was carried out for S phases. The results were then compared with the corresponding P-wave models. As expected, the derived models show that the relatively old and cold Baltic Shield has higher velocity at depth than the younger lithosphere farther South. The models show two sharp and distinct increases in depth to velocities which are low compared to our reference model, as we move from South to North. The location and sharpness of these boundaries suggests that the features resolved are, at least partially, compositional in origin, presumably related to mantle depletion. A sharp and steep subcrustal boundary is found roughly coincident with the southern edge of Sweden. This is below where the edge of the Baltic Shield is usually placed, based on surface geological evidence (the Sorgenfrei–Tornquist Zone). Another less significant transition is recognised more or less beneath the Elbe-lineament. Relatively high d( V p / V s) ratios under the central part of the profile (Denmark) indicate relatively low S-velocity in an area where a gravity high supports the hypothesis of extensive mafic intrusions.

Upper-mantle seismic structure in a region of incipient continental breakup: northern Ethiopian rift

SUMMARY The northern Ethiopian rift forms the third arm of the Red Sea, Gulf of Aden triple junction, and marks the transition from continental rifting in the East African rift to incipient oceanic spreading in Afar. We determine the P- and S-wave velocity structure beneath the northern Ethiopian rift using independent tomographic inversion of P- and S-wave relative arrival-time residuals from teleseismic earthquakes recorded by the Ethiopia Afar Geoscientific Lithospheric Experiment (EAGLE) passive experiment using the regularised non-linear least-squares inversion method of VanDecar. Our 79 broad-band instruments covered an area 250 × 350 km centred on the Boset magmatic segment ∼70 km SE of Addis Ababa in the centre of the northern Ethiopian rift. The study area encompasses several rift segments showing increasing degrees of extension and magmatic intrusion moving from south to north into the Afar depression. Analysis of relative arrival-time residuals shows that the rift flanks are asymmetric with arrivals associated with the southeastern Somalian Plate faster (∼0.65 s for the P waves; ∼ 2s for the S waves) than the northwestern Nubian Plate. Our tomographic inversions image a 75 km wide tabular low-velocity zone (δVP ≈− 1.5 per cent, δVS ≈− 4 per cent) beneath the less-evolved southern part of the rift in the uppermost 200‐250 km of the mantle. At depths of >100 km, north of 8.5 ◦ N, this low-velocity anomaly broadens laterally and appears to be connected to deeper low-velocity structures under the Afar depression. An off-rift low-velocity structure extending perpendicular to the rift axis correlates with the eastern limit of the E‐W trending reactivated Precambrian Ambo‐Guder fault zone that is delineated by Quaternary eruptive centres. Along axis, the low-velocity upwelling beneath the rift is segmented, with low-velocity material in the uppermost 100 km often offset to the side of the rift with the highest rift flank topography. Our observations from this magmatic rift zone, which is transitional between continental and oceanic rifting, do not support detachment fault models of lithospheric extension but instead point to strain accommodation via magma assisted rifting.

Regional teleseismic tomography of the western Lachlan Orogen and the Newer Volcanic Province, southeast Australia

SUMMARY From 1998 May to September a portable array of 40 short-period digital seismograph stations was operated in western Victoria, southeast Australia, across the western end of the midPaleozoic Lachlan Foldbelt and the Newer Volcanic Province. Consisting of four parallel, almost W‐E-oriented receiver lines, the array covered an area of about 270 × 150 km 2 . The major aim of the LF98 (Lachlan Foldbelt survey 1998) project is to map lateral variations in P-wave speeds (Vp) in the crust and upper mantle using teleseismic arrival time tomography, primarily in order to investigate whether the major surface structural zones are associated with seismic velocity signatures at depth. Little a priori information from seismic profiling is available. We invert 4067 relative arrival time residuals for a minimum structure Vp model in the upper few hundred km using non-linear iteration and 3-D ray tracing. The most prominent negative anomaly (∼3.8 per cent) in Vp is found at a depth of about 45 km underneath the eastern part of the Newer Volcanic Province. It correlates spatially with the highest density of Pliocene and Pleistocene eruption centres northwest of Melbourne, and is therefore interpreted as a hotspot-related high-temperature anomaly causing reduced mantle velocities. The related coherent volume of significantly lower than average velocities extends down to depths greater than 100 km in the east, and extends west underneath the Newer Volcanic Province. A strong velocity contrast, with average velocities ∼2 per cent greater in the west, is found down to about 100 km across the Moyston Fault Zone, which forms the major structural boundary between the early-Paleozoic Delamerian Orogen in the west and the Lachlan Orogen in the east. This result suggests that the Moyston Fault Zone should be seen as a major lithospheric boundary. In the south this boundary is also expressed by a distinct discontinuity in Sr-isotopic ratios of xenoliths (the so-called Mortlake discontinuity) and a change in the geochemistry of plutons of similar age. However, if the east to west velocity contrast originally existed in this southern zone, it is now overprinted by the thermally reduced mantle velocities beneath the Newer Volcanic Province.