Seismic Tomographic Images Research Articles

Flat versus normal subduction zones: a comparison based on 3-D regional traveltime tomography and petrological modelling of central Chile and western Argentina (29°–35°S)

Our study compares the seismic properties between the flat and normal subduction regions in central Chile, to better understand the links between the slab geometry, surface deformation and the deeper structures. In comparison with previous studies, we show the most complete 3-D regional seismic tomography images for this region, in which we use (1) a larger seismic data set compiled from several short-term seismic catalogues, (2) a denser seismic array allowing a better resolution of the subduction zone from the trench to the backarc and into the upper ∼30 km of the slab and (3) a starting 1-D background velocity model specifically calculated for this region and refined over the years. We assess and discuss our tomography results using regional seismic attenuation models and estimating rock types on the basis of pressure and temperature conditions computed from thermomechanical models. Our results show significant seismic differences between the flat and normal subduction zones. As expected, the faster seismic velocities and increased seismicity within the flat slab and overriding lithosphere are generally consistent with a cooler thermal state. Our results are also consistent with dehydration of the mantle above the subducted Juan Fernandez Ridge at the eastern tip of the flat slab segment, indicating that the latter retains some fluids during subduction. However, fluids in the upper portion of the flat slab segment are not seismically detected, since we report instead fast slab seismic velocities which contradict the argument of its buoyancy being the cause of horizontal subduction. The forearc region, above the flat slab, exhibits high Vs and very low Vp/Vs ratios, uncorrelated with typical rock compositions, increased density or reduced temperature; this feature is possibly linked with the aftershock effects of the Mw7.1 1997 Punitaqui earthquake, the flat slab geometry and/or seismic anisotropy. At the surface, the seismic variations correlate with the geological terranes. The Andean crust is strongly reduced in seismic velocities along the La Ramada–Aconcagua deformation belt, suggesting structural damage. Slow seismic velocities along the Andean Moho match non-eclogitized hydrated rocks, consistent with a previous delamination event or a felsic composition, which in turn supports the extent of the Chilenia terrane at these depths. We confirm previous studies that suggest that the Cuyania terrane in the backarc region is mafic and contains an eclogitized lower crust below 50-km depth. We also hypothesize major Andean basement detachment faults (or shear zones) to extend towards the plate interface and canalize slab-derived fluids into the continental crust.

Emplacement ages and sources of kimberlites and related rocks in southern Africa: U–Pb ages and Sr–Nd isotopes of groundmass perovskite

Groundmass perovskite has been dated by LA-ICPMS in 135 kimberlites and related rocks from 110 localities across southern Africa. Sr and/or Nd isotopes have been analysed by LA-MC-ICPMS in a subset of these and integrated with published data. The age distribution shows peaks at 1,600–1,800, 1,000–1,200, 500–800 and 50–130 Ma. The major “bloom” of Group I kimberlites at ca 90 ± 10 Ma was preceded by a slow build-up in magmatic activity from ca 180 Ma. The main pulse of Group II kimberlites at 120–130 Ma was a distinct episode within this build-up. Comparison of the isotopic data with seismic tomography images suggests that metasomatized subcontinental lithospheric mantle (SCLM) with very low e Nd and high 87Sr/86Sr, (the isotopic signature of Group II kimberlites) was focused in low-Vs zones along translithospheric structures. Such metasomatized zones existed as early as 1,800 Ma, but were only sporadically tapped until the magmatic build-up began at ca 180 Ma, and contributed little to the kimberlitic magmas after ca 110 Ma. We suggest that these metasomatized volumes resided in the deep SCLM and that their low-melting point components were “burned off” by rising temperatures, presumably during an asthenospheric upwelling that led to SCLM thinning and a rise in the ambient geotherm between 120 and 90 Ma. The younger Group I kimberlites therefore rarely interacted with such SCLM, but had improved access to shallower volumes of differently metasomatized, ancient SCLM with low 87Sr/86Sr and intermediate e Nd (0–5). The kimberlite compositions therefore reflect the evolution of the SCLM of southern Africa, with metasomatic-enrichment events from as early as 1.8 Ga, through a major thermal and compositional change at ca 110 Ma, and the major kimberlite “bloom” around 90 Ma.

Inside Mt. Vesuvius: a new method to look at the seismic (velocity and attenuation) tomographic imaging

<p>New velocity and attenuation images of the geological structures below Mt. Vesuvius have been obtained using the programming facilities as well as the enhanced graphical power of Mathematica<span><sup>8TM</sup></span>. The velocity and attenuation space distributions, already calculated inverting respectively P-wave travel times and amplitude spectra of local VT quakes, are first optimally interpolated and then graphically represented in a new Mathematica<span><sup>8TM</sup></span> code notebook (a powerful computational document with more facilities than a simple code) developed by the present authors. The notebook aims at interactively and friendly representing 3D volume distributions of velocity and attenuation parameters. The user can easily obtain vertical sections (N-S, E-W, NE-SW and NW-SE oriented) and define color scales to represent velocity or attenuation variations or prefer iso-surface plots to represent the pattern of peculiar geological structures. The use of dynamic graphical representation, allowing the sliding of any (horizontal and/or vertical) slice through the volume under study, gives an unusual and powerful vision of any small velocity or attenuation anomaly. The (open source) code, coupled with the friendly use of internal routines of Mathematica, allows to adapt the graphical representation to any user necessity. The method appears to be particularly adapt to represent attenuation images, where the space variations of the parameters are strong with respect to their average. The 3-D plots of the interpolated velocity and attenuation fields enhance the image of Mt. Vesuvius structure, evidencing low-velocity associated with high attenuation anomalies which appeared unfocused in the plots reported by Scarpa et al. [2002] and De Siena et al. [2009].</p>

New insights into Mt. Vesuvius hydrothermal system and its dynamic based on a critical review of seismic tomography and geochemical features

<p>The seismic velocity and attenuation tomography images, calculated inverting respectively P-wave travel times and amplitude spectra of local VT quakes at Mt. Vesuvius have been reviewed and graphically represented using a new software recently developed using Mathematica<span><sup>8TM</sup></span>. The 3-D plots of the interpolated velocity and attenuation fields obtained through this software evidence low-velocity volumes associated with high attenuation anomalies in the depth range from about 1 km to 3 km below the sea level. The heterogeneity in the distribution of the velocity and attenuation values increases in the volume centred around the crater axis and laterally extended about 4 km, where the geochemical interpretation of the data from fumarole emissions reveals the presence of a hydrothermal system with temperatures as high as 400-450°C roughly in the same depth range (1.5 km to 4 km). The zone where the hydrothermal system is space-confined possibly hosted the residual magma erupted by Mt. Vesuvius during the recent eruptions, and is the site where most of the seismic energy release has occurred since the last 1944 eruption.</p>

Upper mantle compositional variations and discontinuity topography imaged beneath Australia from Bayesian inversion of surface‐wave phase velocities and thermochemical modeling

Here we discuss the nature of velocity heterogeneities seen in seismic tomography images of Earth's mantle whose origins and relation to thermochemical variations are yet to be understood. We illustrate this by inverting fundamental‐mode and higher‐order surface‐wave phase velocities for radial models of the thermochemical and anisotropic structure of the mantle to 450 km depth. Dispersion data are linked to thermochemical parameters through a thermodynamic formalism for computing mantle mineral phase equilibria and physical properties. The inverse problem is solved using a probabilistic inference approach whereby robust uncertainty estimates are obtained. We find that both compositional and thermal anomalies are required if observations are to be satisfied. Mantle thermochemical variations extend to 250 km depth beneath western and central Australia and are characterized by increased Mg/Fe and Mg/Si values relative to surrounding mantle. Correlated herewith are thermal variations that closely follow surface tectonics. We also observe a strong contribution to lateral variations in structure and topography across the “410 km” seismic discontinuity from thermochemically induced phase transformations that appear much stronger than lateral variations immediately above and below. Inside the transition zone, a general decoupling of structure relative to that of the upper mantle occurs driven by a reversal in Mg/Si, while thermal anomalies are smoothed out. Comparison of presently derived shear‐wave tomography models with other regional models is encouraging. Radial anisotropy is strongest at 150/200 km depth beneath oceanic/continental areas, respectively, and appears weak and homogeneous below. Finally, geoid anomalies are computed for a subset of sampled model and compared to observations.

Plate reconstructions in the Arctic region based on joint analysis of gravity, magnetic, and seismic anomalies

Based on the analysis of various geophysical data, namely, free-air gravity anomalies, magnetic anomalies, upper mantle seismic tomography images, and topography/bathymetry maps, we single out the major structural elements in the Circum Arctic and present the reconstruction of their locations during the past 200 million years. The configuration of the magnetic field patterns allows revealing an isometric block, which covers the Alpha–Mendeleev Ridges and surrounding areas. This block of presumably continental origin is the remnant part of the Arctida Plate, which was the major tectonic element in the Arctic region in Mesozoic time. We believe that the subduction along the Anyui suture in the time period from 200 to 120 Ma caused rotation of the Arctida Plate, which, in turn, led to the simultaneous closure of the South Anyui Ocean and opening of the Canadian Basin. The rotation of this plate is responsible for extension processes in West Siberia and the northward displacement of Novaya Zemlya relative to the Urals–Taimyr orogenic belt. The cratonic-type North American, Greenland, and European Plates were united before 130 Ma. At the later stages, first Greenland was detached from North America, which resulted in the Baffin Sea, and then Greenland was separated from the European Plate, which led to the opening of the northern segment of the Atlantic Ocean. The Cenozoic stage of opening of the Eurasian Basin and North Atlantic Ocean is unambiguously reconstructed based on linear magnetic anomalies. The counter-clockwise rotation of North America by an angle of ~15° with respect to Eurasia and the right lateral displacement to 200–250 km ensure an almost perfect fit of the contours of the deep water basin in the North Atlantic and Arctic Oceans.

Double layering of a thermochemical plume in the upper mantle beneath Hawaii

According to classical plume theory, purely thermal upwellings rise through the mantle, pond in a thin layer beneath the lithosphere, and generate hotspot volcanism. Neglected by this theory, however, are the dynamical effects of compositional heterogeneity carried by mantle plumes even though this heterogeneity has been commonly identified in sources of hotspot magmas. Numerical models predict that a hot, compositionally heterogeneous mantle plume containing a denser eclogite component tends to pool at ∼300–410 km depth before rising to feed a shallower sublithospheric layer. This double-layered structure of a thermochemical plume is more consistent with seismic tomographic images at Hawaii than the classical plume model. The thermochemical structure as well as time dependence of plume material rising from the deeper into the shallower layer can further account for long-term fluctuations in volcanic activity and asymmetry in bathymetry, seismic structure, and magma chemistry across the hotspot track, as are observed.

Caveats on tomographic images

AbstractGeological and geodynamic models of the mantle often rely on joint interpretations of published seismic tomography images and petrological/geochemical data. This approach tends to neglect the fundamental limitations of, and uncertainties in, seismic tomography results. These limitations and uncertainties involve theory, correcting for the crust, the lack of rays throughout much of the mantle, the difficulty in obtaining the true strength of anomalies, choice of what background model to subtract to reveal anomalies, and what cross‐sections to select for publication. The aim of this review is to provide a relatively non‐technical summary of the most important of these problems, collected together in a single paper, and presented in a form accessible to non‐seismologists. Appreciation of these issues is essential if final geodynamic models are to be robust, and required by the scientific observations.

Crustal complexity in the Lachlan Orogen revealed from teleseismic receiver functions

There is an ongoing debate about the tectonic evolution of southeast Australia, particularly about the causes and nature of its accretion to a much older Precambrian core to the west. Seismic imaging of the crust can provide useful clues to address this issue. Seismic tomography imaging is a powerful tool often employed to map elastic properties of the Earth's lithosphere, but in most cases does not constrain well the depth of discontinuities such as the Mohorovičić (Moho). In this study, an alternative imaging technique known as receiver function (RF) has been employed for seismic stations near Canberra in the Lachlan Orogen to investigate: (i) the shear-wave-velocity profile in the crust and uppermost mantle, (ii) variations in the Moho depth beneath the Lachlan Orogen, and (iii) the nature of the transition between the crust and mantle. A number of styles of RF analyses were conducted: H-K stacking to obtain the best compressional–shear velocity (V P /V S) ratio and crustal thickness; nonlinear inversion for the shear-wave-velocity structure and inversion of the observed variations in RFs with back-azimuth to investigate potential dipping of the crustal layers and anisotropy. The thick crust (up to 48 km) and the mostly intermediate nature of the crust−mantle transition in the Lachlan Orogen could be due to the presence of underplating at the base of the crust, and possibly to the existing thick piles of Ordovician mafic rocks present in the mid and lower crust. Results from numerical modelling of RFs at three seismic stations (CAN, CNB and YNG) suggest that the observed variations with back-azimuth could be related to a complex structure beneath these stations with the likelihood of both a dipping Moho and crustal anisotropy. Our analysis reveals crustal thickening to the west beneath CAN station which could be due to slab convergence. The crustal thickening may also be related to the broad Macquarie volcanic arc, which is rooted to the Moho. The crustal anisotropy may arise from a strong N–S structural trend in the eastern Lachlan Orogen and to the preferred crystallographic orientation of seismically anisotropic minerals in the lower and middle crust related to the paleo-Pacific plate convergence.

Why has the Nazca plate slowed since the Neogene?

The classic example of the not-well-understood rapid change of tectonic plate motion is the increase and then decrease of the convergence rate between the Nazca and South America plates during the past 25–20 m.y. that coincided with the growth of the Andes Mountains. Currently, the decrease in convergence rate is explained either by the increasing load of the Andes or by the appearance of flat slab segments beneath South America. Here, we present an alternative view based on a thermomechanical self-consistent (gravity driven) model of Nazca plate subduction. We explain the changes in the convergence rate as a natural consequence of the Nazca plate penetration into the transition zone and lower mantle after long-term oblique subduction of the Farallon plate. The model is consistent with seismic tomographic images of the Nazca plate beneath South America. Our model also shows that the presence of the Andes does not significantly affect the convergence rate between the Nazca and South America plates.

Tectonostratigraphic models of the Alpine terranes and subduction history of the Hellenides

The tectonostratigraphy of the Hellenic terranes is related to their paleogeodynamic and paleogeographic evolution, which comprises three major stages: (i) Continental rifting in the northern margin of Gondwana, characterised by volcanosedimentary complexes of Late Paleozoic–Triassic. (ii) Continental drifting and contemporaneous oceanic opening of Tethyan basins in between the continental terranes during Triassic–Paleogene with shallow-water carbonate platforms on the continental terranes and ophiolite suites interlayered with pelagic sediments within the Tethyan basins. (iii) Accretion of the tectonostratigraphic terranes with docking along the active European margin, characterised by flysch/melange sedimentation along the trenches from ?early Jurassic to Neogene. Two tectonostratigraphic models can be distinguished: one for the continental terranes/carbonate platforms and another for the oceanic basins. The duration of each geodynamic stage for each terrane is obtained from the chronology of the tectonostratigraphic facies change. The general trend is younger ages observed in the southern terranes and older ages towards the northern terranes. The two alternative tectonostratigraphic models are applied in the two groups of terranes and the chronology of the geodynamic–paleogeographic stages is estimated. The Hellenic subducted slab as shown in the seismic tomographic images is correlated to the terrane paleogeography using the estimated widths and the average subduction rates. The upper 1700km of the slab correspond to the three southern more external terranes whereas the remainingpart of the slab may represent the other more internal terranes. A palinspastic model of the Hellenides is presented taking into account the chronology of each stage for every terrane. Obduction of ophiolites over the southern platforms, is observed in all four oceanic basins. Blueschist exhumation, formation of core complexes and tectonic windows through extensional detachments occur after micro-collision of continental terranes/platforms and their isostatic uplift within the upper plate after their detachment from the subducting slab.

Relative arrival-time upper-mantle tomography and the elusive background mean

SUMMARY The interpretation of seismic tomographic images of upper-mantle seismic wave speed structure is often a matter of considerable debate because the observations can usually be explained by a range of hypotheses, including variable temperature, composition, anisotropy, and the presence of partial melt. An additional problem, often overlooked in tomographic studies using relative as opposed to absolute arrival-times, is the issue of the resulting velocity model’s zero mean. In shield areas, for example, relative arrival-time analysis strips off a background mean velocity structure that is markedly fast compared to the global average. Conversely, in active areas, the background mean is often markedly slow compared to the global average. Appreciation of this issue is vital when interpreting seismic tomographic images: ‘high’ and ‘low’velocityanomaliesshouldnotnecessarilybeinterpreted,respectively,as‘fast’and‘slow’ compared to ‘normal mantle’. This issue has been discussed in the seismological literature in detail over the years, yet subsequent tomography studies have still fallen into the trap of mis-interpreting their velocity models. I highlight here some recent examples of this and provide a simple strategy to address the problem using constraints from a recent global tomographic model, and insights from catalogues of absolute traveltime anomalies. Consultation of such absolute measures of seismic wave speed should be routine during regional tomographic studies, if only for the benefit of the broader Earth Science community, who readily follow the red = hot and slow, blue = cold and fast rule of thumb when interpreting the images for themselves.

Time-dependent convection models of mantle thermal structure constrained by seismic tomography and geodynamics: implications for mantle plume dynamics and CMB heat flux

SUMMARY One of the outstanding problems in modern geodynamics is the development of thermal convection models that are consistent with the present-day flow dynamics in the Earth’s mantle, in accord with seismic tomographic images of 3-D Earth structure, and that are also capable of providing a time-dependent evolution of the mantle thermal structure that is as ‘realistic’ (Earth-like) as possible. A successful realization of this objective would provide a realistic model of 3-D mantle convection that has optimal consistency with a wide suite of seismic, geodynamic and mineral physical constraints on mantle structure and thermodynamic properties. To address this challenge, we have constructed a time-dependent, compressible convection model in 3-D spherical geometry that is consistent with tomography-based instantaneous flow dynamics, using an updated and revised pseudo-spectral numerical method. The novel feature of our numerical solutions is that the equations of conservation of mass and momentum are solved only once in terms of spectral Green’s functions. We initially focus on the theory and numerical methods employed to solve the equation of thermal energy conservation using the Green’s function solutions for the equation of motion, with special attention placed on the numerical accuracy and stability of the convection solutions. A particular concern is the verification of the global energy balance in the dissipative, compressible-mantle formulation we adopt. Such validation is essential because we then present geodynamically constrained convection solutions over billion-year timescales, starting from present-day seismically constrained thermal images of the mantle. The use of geodynamically constrained spectral Green’s functions facilitates the modelling of the dynamic impact on the mantle evolution of: (1) depth-dependent thermal conductivity profiles, (2) extreme variations of viscosity over depth and (3) different surface boundary conditions, in this case mobile surface plates and a rigid surface. The thermal interpretation of seismic tomography models does not provide a radial profile of the horizontally averaged temperature (i.e. the geotherm) in the mantle. One important goal of this study is to obtain a steady-state geotherm with boundary layers which satisfies energy balance of the system and provides the starting point for more realistic numerical simulations of the Earth’s evolution. We obtain surface heat flux in the range of Earth-like values : 37 TW for a rigid surface and 44 TW for a surface with tectonic plates coupled to the mantle flow. Also, our convection simulations deliver CMB heat flux that is on the high end of previously estimated values, namely 13 TW and 20 TW, for rigid and plate-like surface boundary conditions, respectively. We finally employ these two end-member surface boundary conditions to explore the very-long-time scale evolution of convection over billion-year time windows. These billion-year-scale simulations will allow us to determine the extent to which a ‘memory’ of the starting tomography-based thermal structure is preserved and hence to explore the longevity of the structures in the present-day mantle. The two surface boundary conditions, along with the geodynamically inferred radial viscosity profiles, yield steady-state convective flows that are dominated by long wavelengths throughout the lower mantle. The rigid-surface condition yields a spectrum of mantle heterogeneity dominated by spherical harmonic degree 3 and 4, and the plate-like surface condition yields a pattern dominated by degree 1. Our exploration of the time-dependence of the spatial heterogeneity shows that, for both types of surface boundary condition, deep-mantle hot upwellings resolved in the present-day tomography model are durable and stable features. These deeply rooted mantle plumes show remarkable longevity over very long geological time spans, mainly owing to the geodynamically inferred high viscosity in the lower mantle.

Juan de Fuca subduction zone from a mixture of tomography and waveform modeling

Seismic tomography images of the upper mantle structures beneath the Pacific Northwestern United States display a maze of high‐velocity anomalies, many of which produce distorted waveforms evident in the USArray observations indicative of the Juan de Fuca (JdF) slab. The inferred location of the slab agrees quite well with existing contour lines defining the slab's upper interface. Synthetic waveforms generated from a recent tomography image fit teleseismic travel times quite well and also some of the waveform distortions. Regional earthquake data, however, require substantial changes to the tomographic velocities. By modeling regional waveforms of the 2008 Nevada earthquake, we find that the uppermost mantle of the 1D reference model AK135, the reference velocity model used for most tomographic studies, is too fast for the western United States. Here, we replace AK135 with mT7, a modification of an older Basin‐and‐Range model T7. We present two hybrid velocity structures satisfying the waveform data based on modified tomographic images and conventional slab wisdom. We derive P and SH velocity structures down to 660 km along two cross sections through the JdF slab. Our results indicate that the JdF slab is subducted to a depth of 250 km beneath the Seattle region, and terminates at a shallower depth beneath Portland region of Oregon to the south. The slab is about 60 km thick and has a P velocity increase of 5% with respect to mT7. In order to fit waveform complexities of teleseismic Gulf of Mexico and South American events, a slab‐like high‐velocity anomaly with velocity increases of 3% for P and 7% for SH is inferred just above the 660 discontinuity beneath Nevada.

Long-term preservation of slab signatures in the mantle inferred from hydrogen isotopes

Subduction transports water into the mantle, but it is uncertain whether the water is preserved in the slab or is rapidly diffused. Analysis of hydrogen and boron isotopes in volcanic rocks sourced from an ancient subducted slab beneath the southwestern Pacific Ocean provides evidence for the long-term preservation of subducted water in the mantle. Seismic tomographic images indicate that subducted lithosphere is transported into the deep mantle1. Petrologic modelling shows that water contained in subducted slabs can be carried to depths of at least 200 km (ref. 2); however, whether the hydrated slab signature is preserved at greater depths depends on diffusion rates. Experimental studies give conflicting results on the question of hydrogen preservation. On a small scale, hydrogen equilibration with ambient mantle should be rapid3,4, implying that the slab hydrogen signature may not be preserved in the deep mantle5. However, on large scales the time required for diffusive equilibration is longer and hydrogen anomalies may persist6,7. Here we present hydrogen and boron data from submarine volcanic glasses erupted in the Manus back-arc basin, southwestern Pacific Ocean. We find that samples with low hydrogen-isotope values also exhibit the geochemical signature of dehydrated, subducted lithosphere. Combined with additional geochemical and geophysical data, we interpret this as direct evidence for the preservation of hydrogen anomalies in an ancient slab in the mantle. Our geochemical data are consistent with experimental estimates of diffusion for the upper mantle6 and transition zone7. We conclude that hydrogen anomalies can persist in the mantle without suffering complete diffusive equilibration over timescales of up to a billion years.

The Pliocene‐Quaternary wedge‐top basins of southern Italy: an expression of propagating lateral slab tear beneath the Apennines

AbstractIn recent years, contrasting seismic tomographic images have given rise to an extensive debate about the occurrence and implications of migrating slab detachment beneath southern Italy. One of the most pertinent aspects of this process is the concentration of the slab pull force, and particularly its surface expression in terms of vertical motions and related basin subsidence/uplift. In this study we focused on shallow‐water to continental, Pliocene‐Quaternary basins that formed on top of the Apennine allochthonous wedge after its emplacement onto a large foreland carbonate platform domain (Apulian Platform). Due to the thick‐skinned style of deformation controlling the Pliocene‐Pleistocene stages of continental shortening, a high degree of coupling with the downgoing plate appears to characterize the late tectonic evolution of the southern Apennines. Therefore, the wedge‐top basins analysed in this study, although occurring on the deformed edge of the overriding plate, are capable of recording deep geodynamic processes affecting the slab. Detailed stratigraphic work on these wedge‐top basins points to a progressive SE‐ward migration of basin subsidence from c. 4 to c. 2.8 Ma over a distance of about 140 km along the strike of the Apennine belt. Such a migration is consistent with a redistribution of slab‐pull forces associated with the progressive lateral migration at a mean rate in the range of 12–14 cm y–1 of a slab tear within the down‐going Adriatic lithosphere. These results yield fundamental information on the rates of first‐order geodynamic processes affecting the slab, and on related surface response.

Seismic tomographic imaging of the crust and upper mantle under the central and western Tien Shan orogenic belt

[1] Detailed 3-D tomographic images of P and S wave velocity (Vp, Vs) and Poisson's ratio (σ) under the central and western Tien Shan orogenic belt are determined by using a large number of high-quality P and S wave arrival times from local earthquakes. The results show that under the Tien Shan orogenic belt high-Vp, high-Vs, and low-σ anomalies are revealed in the upper and middle crust, possibly indicating the existence of the Paleozoic crystalline basement rocks, while low-Vp, low-Vs, and high-σ anomalies appear in the lower crust and upper mantle, perhaps suggesting that the hot and wet material is upwelling under the Tien Shan orogenic belt from the mantle. Some high-Vp, high-Vs, and low-σ anomalies are tilted toward the Tien Shan along with the seismicity. These are found in the collision zones between the Tien Shan and the Tarim basin in the south and the Kazakh shield in the north and suggest the underthrusting of the Tarim and Kazakh lithosphere beneath the Tien Shan. Meanwhile, some low-Vp, low-Vs, and high-σ anomalies are imaged in other parts of these collision zones, perhaps indicating the intrusion of the hot and wet material into the crust from the upper mantle. These results indicate that both the upwelling of the hot and wet material and the underthrusting of the Tarim and Kazakh lithosphere may have played an important role in the mountain building. Under the Tarim and Fergana basins, low-Vp, low-Vs, and high-σ anomalies are revealed in the upper crust, while high-Vp, high-Vs, and low-σ anomalies are visible in the lower crust and upper mantle. These may reflect the existence of less compacted sedimentary material in the shallow crust and more highly compacted craton-like structures in the deeper crust and upper mantle under the basins. The Talas-Fergana fault shows an obvious tectonic boundary between central and western Tien Shan. The central Tien Shan displays high-Vp, high-Vs, and low-σ anomalies in the upper and middle crust, while western Tien Shan exhibits low-Vp, low-Vs, and high-σ anomalies. However, the pattern of seismic structure between central and western Tien Shan reverses in the lower crust. Such a correlation may extend down to the upper mantle, suggesting that the Talas-Fergana fault may be a lithospheric-scale boundary. Additionally, a columnar low-Vp and low-Vs anomaly is clearly observed around the turning point of the Talas-Fergana fault from the NWN to NWW trending orientations and may indicate that the fault provides a channel for the hot and wet material upwelling from the mantle to the surface.

Crustal-scale cross-sections across the NW Zagros belt: implications for the Arabian margin reconstruction

AbstractQuantified balanced and restored crustal cross-sections across the NW Zagros Mountains are presented in this work integrating geological and geophysical local and global datasets. The balanced crustal cross-section reproduces the surficial folding and thrusting of the thick cover succession, including the near top of the Sarvak Formation (~90 Ma) that forms the top of the restored crustal cross-section. The base of the Arabian crust in the balanced cross-section is constrained by recently published seismic receiver function results showing a deepening of the Moho from 42 ± 2 km in the undeformed foreland basin to 56 ± 2 km beneath the High Zagros. The internal parts of the deformed crustal cross-section are constrained by new seismic tomographic sections imaging a ~50° NE-dipping sharp contact between the Arabian and Iranian crusts. These surfaces bound an area of 10800 km2that should be kept constant during the Zagros orogeny. The Arabian crustal cross-section is restored using six different tectonosedimentary domains according to their sedimentary facies and palaeobathymetries, and assuming Airy isostasy and area conservation. While the two southwestern domains were directly determined from well-constrained surface data, the reconstruction of the distal domains to the NE was made using the recent margin model of Wrobel-Daveauet al. (2010) and fitting the total area calculated in the balanced cross-section. The Arabian continental–oceanic boundary, at the time corresponding to the near top of the Sarvak Formation, is located 169 km to the NE of the trace of the Main Recent Fault. Shortening is estimated at ~180 km for the cover rocks and ~149 km for the Arabian basement, including all compressional events from Late Cretaceous to Recent time, with an average shortening rate of ~2 mm yr−1for the last 90 Ma.

Continuing Colorado plateau uplift by delamination-style convective lithospheric downwelling

The Colorado plateau is a large, tectonically intact, physiographic province in the southwestern North American Cordillera that stands at ∼1,800-2,000 m elevation and has long been thought to be in isostatic equilibrium. The origin of these high elevations is unclear because unlike the surrounding provinces, which have undergone significant Cretaceous-Palaeogene compressional deformation followed by Neogene extensional deformation, the Colorado plateau is largely internally undeformed. Here we combine new seismic tomography and receiver function images to resolve a vertical high-seismic-velocity anomaly beneath the west-central plateau that extends more than 200 km in depth. The upper surface of this anomaly is seismically defined by a dipping interface extending from the lower crust to depths of 70-90 km. The base of the continental crust above the anomaly has a similar shape, with an elevated Moho. We interpret these seismic structures as a continuing regional, delamination-style foundering of lower crust and continental lithosphere. This implies that Pliocene (2.6-5.3 Myr ago) uplift of the plateau and the magmatism on its margins are intimately tied to continuing deep lithospheric processes. Petrologic and geochemical observations indicate that late Cretaceous-Palaeogene (∼90-40 Myr ago) low-angle subduction hydrated and probably weakened much of the Proterozoic tectospheric mantle beneath the Colorado plateau. We suggest that mid-Cenozoic (∼35-25 Myr ago) to Recent magmatic infiltration subsequently imparted negative compositional buoyancy to the base and sides of the Colorado plateau upper mantle, triggering downwelling. The patterns of magmatic activity suggest that previous such events have progressively removed the Colorado plateau lithosphere inward from its margins, and have driven uplift. Using Grand Canyon incision rates and Pliocene basaltic volcanism patterns, we suggest that this particular event has been active over the past ∼6 Myr.

On predicting mantle mushroom plumes

This study investigates the mechanism of formation of convection plumes of mushroom shape in sub-solidus mantle and their prediction. The seismic-tomographic images of columnar structures of several hundreds kilometers in diameter have been reported by several researchers, while the much cherished mushroom-shaped plume heads could only be found in computational geodynamics (CGD) models and simple small-scale laboratory analogue simulations. Our theory of transient instability shows that the formation of convection plumes is preceded by the onset of convection caused by unsteady-state heat conduction at the boundaries, from which filamentous plumes first appear. The plumes generated at the Core Mantle Boundary (CMB) and lithosphere rising and falling through the mantle have been predicted simply with our theory for various heat fluxes and viscosities, which still remain uncertain amongst geoscientists. The sizes of mushroom plumes in the sub-solidus mantle caused by heat fluxes of 20 and 120 mW/m 2 at the CMB are found to be 1842 km and 1173 km with critical times over 825 Myr and 334 Myr respectively. They are comparable to some large continental flood basalt provinces, and they number between 17 and 41. The thickness of the thermal boundary layers at the CMB from which convection plumes evolved are found to be 652 km and 415 km for 20 and 120 mW/m 2 respectively. Top cooling may produce plunging plumes of diameter of 585 km and at least 195 Myr old. The number of cold plumes is estimated to be 569, which has not been observed by seismic tomography or as cold spots. The cold plunging plumes may overwhelm and entrap some of the hot rising plumes from CMB, so that together they may settle in the transition zone.