Age-dependent seismic thickness and mechanical strength of the Australian lithosphere

We present constraints on the regional variations of the seismic and mechanical thickness of the Australian lithosphere. We infer the seismic thickness from a waveform tomographic model of S-wave speed, and as a proxy for the elastic thickness we use the wavelength at which the coherence of surface topography and Bouguer gravity drops below half of its long-wavelength maximum. Our results show that on scales <1000 km the relationship between the age of the crust and the thickness of the lithosphere is more complicated than longer-wavelength or global averages suggest. Recent geochemical and geodynamical evidence for small-scale secular variations of the composition and stability of continental cratons further illustrates the complexity of the age dependence of seismo-mechanical lithospheric properties on regional scales

Mesozoic subducted slabs under Siberia

Recent results from seismic tomography demonstrate that subducted oceanic lithosphere can be observed globally as slabs of relatively high seismic velocity in the upper as well as lower mantle(1,2). The Asian mantle is no exception, with high-velocity slabs being observed downwards from the west Pacific subduction zones under the Kurile Islands, Japan and farther south(3-5), as well as under Asia's ancient Tethyan margin. Here we present evidence for the presence of slab remnants of Jurassic age that were subducted when the Mongol-Okhotsk and Kular-Nera oceans closed between Siberia, the combined Mongolia-North China blocks and the Omolon block(6-8). We identify these proposed slab remnants in the lower mantle west of Lake Baikal down to depths of at least 2,500 km, where they join what has been interpreted as a 'graveyard'(9) of subducted lithosphere at the bottom of the mantle. Our interpretation implies that slab remnants in the mantle can still be recognized some 150 million years or more after they have been subducted and that such structures may be useful in associating geodynamic to surface-tectonic processes.Peer Reviewedhttp://deepblue.lib.umich.edu/bitstream/2027.42/62524/1/397246a0.pd

Multiple volcanic episodes of flood basalts caused by thermochemical mantle plumes

The hypothesis that a single mushroom-like mantle plume head can generate a large igneous province within a few million years has been widely accepted(1). The Siberian Traps at the Permian Triassic boundary(2) and the Deccan Traps at the Cretaceous Tertiary boundary(3) were probably erupted within one million years. These large eruptions have been linked to mass extinctions. But recent geochronological data(4-11) reveal more than one pulse of major eruptions with diverse magma flux within several flood basalts extending over tens of million years. This observation indicates that the processes leading to large igneous provinces are more complicated than the purely thermal, single-stage plume model suggests. Here we present numerical experiments to demonstrate that the entrainment of a dense eclogite-derived material at the base of the mantle by thermal plumes can develop secondary instabilities due to the interaction between thermal and compositional buoyancy forces. The characteristic timescales of the development of the secondary instabilities and the variation of the plume strength are compatible with the observations. Such a process may contribute to multiple episodes of large igneous provinces.Peer Reviewedhttp://deepblue.lib.umich.edu/bitstream/2027.42/62705/1/nature03697.pd

Extremely high He isotope ratios in MORB-source mantle from the proto-Iceland plume

The high &lt;sup&gt;3&lt;/sup&gt;He/&lt;sup&gt;4&lt;/sup&gt;He ratio of volcanic rocks thought to be derived from mantle plumes is taken as evidence for the existence of a mantle reservoir that has remained largely undegassed since the Earth's accretion. The helium isotope composition of this reservoir places constraints on the origin of volatiles within the Earth and on the evolution and structure of the Earth's mantle. Here we show that olivine phenocrysts in picritic basalts presumably derived from the proto-Iceland plume at Baffin Island, Canada, have the highest magmatic &lt;sup&gt;3&lt;/sup&gt;He/&lt;sup&gt;4&lt;/sup&gt;He ratios yet recorded. A strong correlation between &lt;sup&gt;3&lt;/sup&gt;He/&lt;sup&gt;4&lt;/sup&gt;He and &lt;sup&gt;87&lt;/sup&gt;Sr/&lt;sup&gt;86&lt;/sup&gt;Sr, &lt;sup&gt;143&lt;/sup&gt;Nd/&lt;sup&gt;144&lt;/sup&gt;Nd and trace element ratios demonstrate that the &lt;sup&gt;3&lt;/sup&gt;He-rich end-member is present in basalts that are derived from large-volume melts of depleted upper-mantle rocks. This reservoir is consistent with the recharging of depleted upper-mantle rocks by small volumes of primordial volatile-rich lower-mantle material at a thermal boundary layer between convectively isolated reservoirs. The highest &lt;sup&gt;3&lt;/sup&gt;He/&lt;sup&gt;4&lt;/sup&gt;He basalts from Hawaii and Iceland plot on the observed mixing trend. This indicates that a &lt;sup&gt;3&lt;/sup&gt;He-recharged depleted mantle (HRDM) reservoir may be the principal source of high &lt;sup&gt;3&lt;/sup&gt;He/&lt;sup&gt;4&lt;/sup&gt;He in mantle plumes, and may explain why the helium concentration of the 'plume' component in ocean island basalts is lower than that predicted for a two-layer, steady-state model of mantle structure

Kimberlites reveal 2.5-billion-year evolution of a deep, isolated mantle reservoir

The widely accepted paradigm of Earth's geochemical evolution states that the successive extraction of melts from the mantle over the past 4.5 billion years formed the continental crust, and produced at least one complementary melt-depleted reservoir that is now recognized as the upper-mantle source of mid-ocean-ridge basalts1. However, geochemical modelling and the occurrence of high 3He/4He (that is, primordial) signatures in some volcanic rocks suggest that volumes of relatively undifferentiated mantle may reside in deeper, isolated regions2. Some basalts from large igneous provinces may provide temporally restricted glimpses of the most primitive parts of the mantle3,4, but key questions regarding the longevity of such sources on planetary timescales—and whether any survive today—remain unresolved. Kimberlites, small-volume volcanic rocks that are the source of most diamonds, offer rare insights into aspects of the composition of the Earth’s deep mantle. The radiogenic isotope ratios of kimberlites of different ages enable us to map the evolution of this domain through time. Here we show that globally distributed kimberlites originate from a single homogeneous reservoir with an isotopic composition that is indicative of a uniform and pristine mantle source, which evolved in isolation over at least 2.5 billion years of Earth history—to our knowledge, the only such reservoir that has been identified to date. Around 200 million years ago, extensive volumes of the same source were perturbed, probably as a result of contamination by exogenic material. The distribution of affected kimberlites suggests that this event may be related to subduction along the margin of the Pangaea supercontinent. These results reveal a long-lived and globally extensive mantle reservoir that underwent subsequent disruption, possibly heralding a marked change to large-scale mantle-mixing regimes. These processes may explain why uncontaminated primordial mantle is so difficult to identify in recent mantle-derived melts

Seismic evidence for a lower-mantle origin of the Iceland plume

Iceland, one of the most thoroughly investigated hotspots1,2,3, is generally accepted to be the manifestation of an upwelling mantle plume4. Yet whether the plume originates from the lower mantle or from a convective instability at a thermal boundary layer between the upper and lower mantle near 660 km depth5,6 remains unconstrained. Tomographic inversions of body-wave delay times show that low seismic velocities extend to at least 400 km depth beneath central Iceland7,8, but cannot resolve structure at greater depth. Here we report lateral variations in the depths of compressional-to-shear wave conversions at the two seismic discontinuities marking the top and bottom of the mantle transition zone beneath Iceland. We find that the transition zone is 20 km thinner than in the average Earth9 beneath central and southern Iceland, but is of normal thickness beneath surrounding areas, a result indicative of a hot and narrow plume originating from the lower mantle.This work was supported by the US NSFPeer Reviewe

Surface wave array tomography in SE Tibet from ambient seismic noise and two-station analysis - II. Crustal and upper-mantle structure

We determine the 3-D shear wave speed variations in the crust and upper mantle in the southeastern borderland of the Tibetan Plateau, SW China, with data from 25 temporary broad-band stations and one permanent station. Interstation Rayleigh wave (phase velocity) dispersion curves were obtained at periods from 10 to 50 s from empirical Green's function (EGF) derived from (ambient noise) interferometry and from 20 to 150 s from traditional two-station (TS) analysis. Here, we use these measurements to construct phase velocity maps (from 10 to 150 s, using the average interstation dispersion from the EGF and TS methods between 20 and 50 s) and estimate from them (with the Neighbourhood Algorithm) the 3-D wave speed variations and their uncertainty. The crust structure, parametrized in three layers, can be well resolved with a horizontal resolution about of 100 km or less. Because of the possible effect of mechanically weak layers on regional deformation, of particular interest is the existence and geometry of low (shear) velocity layers (LVLs). In some regions prominent LVLs occur in the middle crust, in others they may appear in the lower crust. In some cases the lateral transition of shear wave speed coincides with major fault zones. The spatial variation in strength and depth of crustal LVLs suggests that the 3-D geometry of weak layers is complex and that unhindered crustal flow over large regions may not occur. Consideration of such complexity may be the key to a better understanding of relative block motion and patterns of seismicity. © Journal compilation © 2008 RAS

The signal of mantle anisotropy in the coupling of normal modes

We investigate whether the coupling of normal mode (NM) multiplets can help us constrain mantle anisotropy. We first derive explicit expressions of the generalized structure coefficients of coupled modes in terms of elastic coefficients, including the Love parameters describing radial anisotropy and the parameters describing azimuthal anisotropy (Jc, Js, Kc, Ks, Mc, Ms, Bc, Bs, Gc, Gs, Ec, Es, Hc, Hs, Dc and Ds). We detail the selection rules that describe which modes can couple together and which elastic parameters govern their coupling. We then focus on modes of type 0Sl-0Tl+1 and determine whether they can be used to constrain mantle anisotropy. We show that they are sensitive to six elastic parameters describing azimuthal anisotropy, in addition to the two shear-wave elastic parameters L and N (i.e. VSV and VSH). We find that neither isotropic nor radially anisotropic mantle models can fully explain the observed degree two signal. We show that the NM signal that remains after correction for the effect of the crust and mantle radial anisotropy can be explained by the presence of azimuthal anisotropy in the upper mantle. Although the data favour locating azimuthal anisotropy below 400km, its depth extent and distribution is still not well constrained by the data. Consideration of NM coupling can thus help constrain azimuthal anisotropy in the mantle, but joint analyses with surface-wave phase velocities is needed to reduce the parameter trade-offs and improve our constraints on the individual elastic parameters and the depth location of the azimuthal anisotropy. © 2008 The Authors Journal compilation © 2008 RAS

The signal of mantle anisotropy in the coupling of normal modes

We investigate whether the coupling of normal mode (NM) multiplets can help us constrain mantle anisotropy. We first derive explicit expressions of the generalized structure coefficients of coupled modes in terms of elastic coefficients, including the Love parameters describing radial anisotropy and the parameters describing azimuthal anisotropy (J , J , K , K , M , M , B , B , G , G , E , E , H , H , D and D ). We detail the selection rules that describe which modes can couple together and which elastic parameters govern their coupling. We then focus on modes of type S T and determine whether they can be used to constrain mantle anisotropy. We show that they are sensitive to six elastic parameters describing azimuthal anisotropy, in addition to the two shear-wave elastic parameters L and N (i.e. V and V ). We find that neither isotropic nor radially anisotropic mantle models can fully explain the observed degree two signal. We show that the NM signal that remains after correction for the effect of the crust and mantle radial anisotropy can be explained by the presence of azimuthal anisotropy in the upper mantle. Although the data favour locating azimuthal anisotropy below 400km, its depth extent and distribution is still not well constrained by the data. Consideration of NM coupling can thus help constrain azimuthal anisotropy in the mantle, but joint analyses with surface-wave phase velocities is needed to reduce the parameter trade-offs and improve our constraints on the individual elastic parameters and the depth location of the azimuthal anisotropy. © 2008 The Authors Journal compilation © 2008 RAS. c s c s c s c s c s c s c s c s 0 l-0 l+1 SV S