Siberian traps: Hypotheses and seismology data

Siberian traps are the result of huge basalt eruptions which took place about 250 Ma ago over a vast territory of Siberia. The genesis of Siberian traps is attributed to a mantle plume with a center in the region of Iceland or beneath the central Urals in terms of their present coordinates. The eruption mechanism is associated with delamination—replacement of the mantle lithosphere by the deep magma material. The receiver function analysis of the records from the Norilsk seismic station (NRIL) allows comparing these hypotheses with the factual data on the depth structure of the region of Siberian traps. The S-wave velocity section place the seismic lithosphere/asthenosphere boundary (LAB) at a depth of 155–190 km, commensurate with the data for the other cratons. The mantle lithosphere has a high S-wave velocity characteristic of cratons (4.6–4.8 km/s instead of the typical value 4.5 km/s). The seismic boundary, which is located at a depth around 410 km beneath the continents is depressed by ~10 km in the region of the NRIL station. The phase diagram of olivine/wadsleyite transformation accounts for this depression by a 50–100°С increase in temperature. At the depths of 350–400 km, the S-wave velocity drops due to partial melting. A new reduction in the S-wave velocities is observed at a depth of 460 km. The similar anomalies (deepening of the 410-km seismic boundary and low shear wave velocity at depths of 350–400 and 460–500 km, respectively) were previously revealed in the other regions of the Meso-Cenozoic volcanism. In the case of a differently directed drift of the Siberian lithosphere and underlying mantle at depths down to 500 km, these anomalies are barely accountable. In particular, if the mantle at a depth ranging from 200 to 500 km is fixed, the anomalies should be observed at the original locations where they emerged 250 Ma ago, i.e. thousands of km from the Siberian traps. Our seismic data suggest that despite the low viscosity of the asthenosphere, the mantle drift at depths ranging from 200 to 500 km is correlated with the drift of the Siberian lithospheric plate. Furthermore, the position of the mantle plume beneath the Urals is easier to reconcile with the seismic data than its position beneath Iceland because of the Siberian traps being less remote from the Urals.     

Mixed mantle provenance: diverse garnet compositions in polymict peridotites, Kaapvaal craton, South Africa

Garnet compositions are used to understand mantle petrogenesis and to reconstruct the lithostratigraphy of the shallow mantle (<200 km). However, garnets in polymict peridotites from the Kaapvaal craton (>2500 Ma) have a centimeter-scale elemental and stable isotopic variability suggestive of a mixed mantle provenance. The chemical heterogeneity of the garnets is similar to that reported from rocks sampled over a considerable depth and temperature range within the lower lithosphere. For example garnets found in polymict peridotites are similar to garnets found in sheared and granular peridotites, ‘cold’ and ‘hot’ lherzolites, peridotitic (P-type) diamond inclusions, and garnets from polybaric (50-200 km) peridotites (i.e. spinel, garnet and diamond facies). These data indicate that the Kaapvaal cratonic root has been disturbed by complex processes possibly associated with crack propagation and entrainment that juxtaposed garnet-bearing lithologies of diverse petrogenesis, provenance and depth. This has preserved chemical disequilibrium in the high pressure minerals in what is, in effect, a mantle breccia possibly associated with kimberlite precursors.     

Testing the reference Moon model in respect of the thermal regime and chemical composition of the mantle: Thermodynamics versus seismology

The VPREMOON seismic reference Moon model (Garcia et al., 2011) has been tested with respect to the thermal regime and chemical composition of the mantle. Based on a self-consistent thermodynamic approach and petrological models of the lunar mantle covering a wide range of concentrations of CaO, Al₂O₃, and FeO, we convert the P- and S-wave velocity profiles to the temperature–depth profiles. The solution procedure relies on the method of the Gibbs free energy minimization and the equations of state for the mantle material which take into account the effects of phase transformations, anharmonicity, and anelasticity. We find that regardless of the chemical composition, the positive P- and S-wave velocity gradient in the lunar mantle leads to a negative temperature gradient, which has no physical basis. For adequate mantle temperatures, the P- and S-wave velocities should remain almost constant or slightly decrease with depth (especially V_S) as a result of the effects of the temperature, which grows faster than pressure. These findings underscore the importance of the relationship of the thermodynamics and physics of minerals with seismology.     

Upper mantle structure beneath central Western Europe from data on phase velocities of surface waves

The deep structure of the upper mantle is determined from data on phase velocities of Love and Rayleigh waves measured by a differential method on traces between two stations in central Western Europe. One-dimensional velocity structures are first constructed from data of each pair of stations, after which two-dimensional distributions of SH and SV velocities are calculated by the method of two-dimensional tomography from S wave velocities at fixed depths. The results are presented in the form of 2-D vertical structures of the average S wave velocity (S = (SV + SH)/2) constructed along profiles crossing the region in directions of the best resolution. The main structural features are a higher velocity zone at depths of 60–80 km in the area (48°–50°N, 9°–11°E) and a lower velocity zone in the western part of the region at depths of 100–150 km, probably extending farther beyond the studied area.     

A petrogenetic study of the garnet pyroxenite enclaves in spinel peridotite,North Dabieshan,China

Recently, garnet pyroxenite enclaves within peridotites occurring near Raobazhai, Huoshan County, have been discovered. The garnet pyroxenite is small pods, decimeters in size, enclosed within intensively serpentinized peridotites. Major mineral components comprise: garnet (Prp_25–35), sodium augite (Jd_10–25) with a small amount of ilmenite. There are two stages of retrometamorphism: the retrogressive granulite facies mineral assemblage is superimposed by that of amphibolite facies. The host rocks of the garnet pyroxenite are spinel peridotites, including spinel harzburgite and lherzolite. Due to intensive serpentinitization, only 5%–40% of the relic olivine (Fo_92–93) are preserved. The orthopyroxenes are Mg-rich (En_87–93) with bending of cleavages and granulation at their margins showing intracrystalline plasticity. On the basis of garnet-clinopyroxene Fe?Mg exchange equilibrium geothermometry proposed by Ellis & Green (1979) and Krogh (1988)K _D=4.06–5.28;T=793–919°C,P=1.5 GPa are estimated for the garnet pyroxenite. It is inferred that the peridotites are mantle rocks about 60 km in depth. During the exhumation of the orogenic belt, it was tectonically emplaced into the lower crust in the solid state and then uplifted to the shallow depth. Obviously, this kind of garnet pyroxenite must be petrogenetically related to its host rock. The REE distribution pattern and the Ni?Co?Sc diagram reveal that they are chemically equivalent to the basaltic melt and ultramafic residua respectively derived from partial melting of mantle rocks.     

Thermal structure and thickness of the lithospheric mantle underlying the Siberian Craton from the kraton and kimberlit superlong seismic profiles

A self-consistent approach is proposed for the investigation of the thermal conditions, chemical composition, and internal structure of the upper mantle of the Earth. Using this approach, the thermal state of the lithospheric mantle beneath the Siberian Craton (SC) is reconstructed from P velocities, taking into account the phase transitions, anharmonicity, and the effects of anelasticity. The velocities of seismic waves are more sensitive to temperature than to the composition of the mantle rocks, which allows the velocity models to be effectively used for reconstruction of the thermal regime of the mantle. The temperature at depths 100–300 km is reconstructed by inversion of the Kraton and Kimberlit superlong seismic profiles for compositions of the garnet harzburgite, lherzolite, and intermediate composition of garnet peridotite. The averaged temperature in the normal continental mantle is reconstructed by inversion of the IASP91 reference model for depleted and fertile substance. One-dimensional models and two-dimensional thermal fields undergo a substantial fall in temperature (～300–600°C) beneath the Siberian Craton as compared to the temperatures of the continental mantle and paleotemperatures inferred from the thermobarometry of xenoliths. Temperature profiles of the Siberian Craton deduced from seismic data lie between the conductive geotherms of 32.5–40.0 mW/m² and below the P(H)-T values obtained for low- and high-temperature xenoliths from the Mir, Udachnaya, and Obnazhennaya kimberlite pipes. The thickness of the thermal lithosphere estimated from the intersection with the potential adiabat is 300–320 km, which is consistent with the data on heat flows and seismotomographic observations. This provides grounds for the assumption that the low-temperature anomalies (thermal roots of continents) penetrate down to a depth of 300 km. The analysis of the sensitivity of seismic velocity and density to the variations in temperature, pressure, and chemical and phase composition of petrological models shows that recognition of fine differences in chemical composition of the lithospheric rocks by seismic methods is impossible.     

Formation time of the big mantle wedge beneath eastern China and a new lithospheric thinning mechanism of the North China craton—Geodynamic effects of deep recycled carbon

High-resolution P wave tomography shows that the subducting Pacific slab is stagnant in the mantle transition zone and forms a big mantle wedge beneath eastern China. The Mg isotopic investigation of large numbers of mantle-derived volcanic rocks from eastern China has revealed that carbonates carried by the subducted slab have been recycled into the upper mantle and formed carbonated peridotite overlying the mantle transition zone, which becomes the sources of various basalts. These basalts display light Mg isotopic compositions(δ26 Mg = –0.60‰ to –0.30‰) and relatively low87 Sr/86 Sr ratios(0.70314–0.70564) with ages ranging from 106 Ma to Quaternary, suggesting that their mantle source had been hybridized by recycled magnesite with minor dolomite and their initial melting occurred at 300-360 km in depth. Therefore, the carbonate metasomatism of their mantle source should have occurred at the depth larger than 360 km, which means that the subducted slab should be stagnant in the mantle transition zone forming the big mantle wedge before 106 Ma. This timing supports the rollback model of subducting slab to form the big mantle wedge. Based on high P-T experiment results, when carbonated silicate melts produced by partial melting of carbonated peridotite was raising and reached the bottom(180–120 km in depth) of cratonic lithosphere in North China, the carbonated silicate melts should have 25–18 wt% CO2 contents, with lower Si O2 and Al2 O3 contents, and higher Ca O/Al2 O3 values, similar to those of nephelinites and basanites, and have higher εNdvalues(2 to 6). The carbonatited silicate melts migrated upward and metasomatized the overlying lithospheric mantle, resulting in carbonated peridotite in the bottom of continental lithosphere beneath eastern China. As the craton lithospheric geotherm intersects the solidus of carbonated peridotite at 130 km in depth, the carbonated peridotite in the bottom of cratonic lithosphere should be partially melted, thus its physical characters are similar to the asthenosphere and it could be easily replaced by convective mantle. The newly formed carbonated silicate melts will migrate upward and metasomatize the overlying lithospheric mantle. Similarly, such metasomatism and partial melting processes repeat, and as a result the cratonic lithosphere in North China would be thinning and the carbonated silicate partial melts will be transformed to high-Si O2 alkali basalts with lower εNdvalues(to-2). As the lithospheric thinning goes on,initial melting depth of carbonated peridotite must decrease from 130 km to close 70 km, because the craton geotherm changed to approach oceanic lithosphere geotherm along with lithospheric thinning of the North China craton. Consequently, the interaction between carbonated silicate melt and cratonic lithosphere is a possible mechanism for lithosphere thinning of the North China craton during the late Cretaceous and Cenozoic. Based on the age statistics of low δ26 Mg basalts in eastern China, the lithospheric thinning processes caused by carbonated metasomatism and partial melting in eastern China are limited in a timespan from 106 to25 Ma, but increased quickly after 25 Ma. Therefore, there are two peak times for the lithospheric thinning of the North China craton: the first peak in 135-115 Ma simultaneously with the cratonic destruction, and the second peak caused by interaction between carbonated silicate melt and lithosphere mainly after 25 Ma. The later decreased the lithospheric thickness to about70 km in the eastern part of North China craton.     

The transport of water in subduction zones

The transport of water from subducting crust into the mantle is mainly dictated by the stability of hydrous minerals in subduction zones. The thermal structure of subduction zones is a key to dehydration of the subducting crust at different depths. Oceanic subduction zones show a large variation in the geotherm, but seismicity and arc volcanism are only prominent in cold subduction zones where geothermal gradients are low. In contrast, continental subduction zones have low geothermal gradients, resulting in metamorphism in cold subduction zones and the absence of arc volcanism during subduction. In very cold subduction zone where the geothermal gradient is very low(?5?C/km), lawsonite may carry water into great depths of ?300 km. In the hot subduction zone where the geothermal gradient is high(25?C/km), the subducting crust dehydrates significantly at shallow depths and may partially melt at depths of 80 km to form felsic melts, into which water is highly dissolved. In this case, only a minor amount of water can be transported into great depths. A number of intermediate modes are present between these two end-member dehydration modes, making subduction-zone dehydration various. Low-T/low-P hydrous minerals are not stable in warm subduction zones with increasing subduction depths and thus break down at forearc depths of ?60–80 km to release large amounts of water. In contrast, the low-T/low-P hydrous minerals are replaced by low-T/high-P hydrous minerals in cold subduction zones with increasing subduction depths, allowing the water to be transported to subarc depths of 80–160 km. In either case, dehydration reactions not only trigger seismicity in the subducting crust but also cause hydration of the mantle wedge. Nevertheless, there are still minor amounts of water to be transported by ultrahigh-pressure hydrous minerals and nominally anhydrous minerals into the deeper mantle. The mantle wedge overlying the subducting slab does not partially melt upon water influx for volcanic arc magmatism, but it is hydrated at first with the lowest temperature at the slab-mantle interface, several hundreds of degree lower than the wet solidus of hydrated peridotites. The hydrated peridotites may undergo partial melting upon heating at a later time. Therefore, the water flux from the subducting crust into the overlying mantle wedge does not trigger the volcanic arc magmatism immediately.     

Geochemical constraints on the model of the composition and thermal conditions of the Moon according to seismic data

A new method of reconstruction of the temperature profile in the lunar mantle from the velocities of seismic P- and S-waves for different models of chemical composition is developed. The procedure of the solution of an inverse problem is realized with the help of the minimization of the Gibbs free energy and the equations of state of a mantle substance, taking into account phase transformations, anharmonicity, and the effects of inelasticity. The geophysical and geochemical constraints on composition and temperature distribution in Moon’s mantle are established. The upper mantle can be composed of olivine pyroxenite, depleted by low-volatile oxides (∼2 wt % of CaO and Al₂O₃). On the contrary, the lower mantle must be enriched by low-volatile oxides (∼4–6 wt % of CaO and Al₂O₃). Its composition can be represented by a mineral association of the olivine + clinopyroxene + garnet or olivine + orthopyroxene + clinopyroxene + garnet type, which is close in composition to pyrolite. The temperature distribution at depths 50–1000 km are approximated by the equation: T(°C) = 351 + 1718[1–exp (−0.00082H)]. The constraints inferred make it possible to conclude that the published values of the velocities of P- and S-waves for the lunar mantle, obtained by processing the data of seismic experiments of the Apollo lunar mission are inconsistent with each other at depths below 300 km. Otherwise, the variations in the velocities of P- and S-waves disturb the symmetry between the petrological model (composition), the temperature profile, and the seismic profile.     

Crustal and upper mantle velocity model along the DOBRE-4 profile from North Dobruja to the central region of the Ukrainian Shield: 1. seismic data

For studying the structure of the lithosphere in southern Ukraine, wide-angle seismic studies that recorded the reflected and refracted waves were carried out under the DOBRE-4 project. The field works were conducted in October 2009. Thirteen chemical shot points spaced 35–50 km apart from each other were implemented with a charge weight varying from 600 to 1000 kg. Overall 230 recording stations with an interval of 2.5 km between them were used. The high quality of the obtained data allowed us to model the velocity section along the profile for P- and S-waves. Seismic modeling was carried out by two methods. Initially, trial-and-error ray tracing using the arrival times of the main reflected and refracted P- and S-phases was conducted. Next, the amplitudes of the recorded phases were analyzed by the finite-difference full waveform method. The resulting velocity model demonstrates a fairly homogeneous structure from the middle to lower crust both in the vertical and horizontal directions. A drastically different situation is observed in the upper crust, where the V _p velocities decrease upwards along the section from 6.35 km/s at a depth of 15–20 km to 5.9–5.8 km/s on the surface of the crystalline basement; in the Neoproterozoic and Paleozoic deposits, it diminishes from 5.15 to 3.80 km/s, and in the Mesozoic layers, it decreases from 2.70 to 2.30 km/s. The subcrustal V _p gradually increases downwards from 6.50 to 6.7–6.8 km/s at the crustal base, which complicates the problem of separating the middle and lower crust. The V _p velocities above 6.80 km/s have not been revealed even in the lowermost part of the crust, in contrast to the similar profiles in the East European Platform. The Moho is clearly delineated by the velocity contrast of 1.3–1.7 km/s. The alternating pattern of the changes in the Moho depths corresponding to Moho undulations with a wavelength of about 150 km and the amplitude reaching 8 to 17 km is a peculiarity of the velocity model.     

Modeling of the thermal structure of continental lithosphere

The thermal structure of continental lithosphere (the temperature, heat flows, and heat generation in the crust and lithosphere) is reconstructed from geothermal, seismic, and petrologic data. The first step is the determination of the temperature profile from absolute P and S wave velocities (T _{P, S} ). The T _{P, S} profile is then adjusted to a thermophysical model of conductive transfer. In addition, the surface heat flow and the T _{P, S} profile are used to determine heat generation, thicknesses of crustal layers, and heat flow components in the crust and lithosphere. A feature inherent in the solution of the thermophysical inverse problem obtained in this paper is the use of constraints derived from the temperature reconstruction by seismic data inversion. As a result, the analytical dependence of the temperature on depth, the intensity of radiogenic heat sources in the crust, and heat flow components in the crust and lithosphere are determined.     

Characteristics of subcontinental lithospheric mantle beneath Baegryeong Island, Korea: Spinel peridotite xenoliths

Abstract   Geochemical data on Baegryeong Island spinel peridotites found in Miocene alkali basalt provide the information for lithosphere composition, chemical processes, equilibrium pressure and temperature conditions. Spinel peridotite xenoliths, showing transitional textures between protogranular and porpyroclastic textures, were accidentally trapped by the ascending alkali basalt magma. The xenoliths originate at depths from 50 to 70 km with a temperature range from 800 to 1100°C. The variations of modal and mineral compositions of the spinel peridotite xenoliths indicate that the xenoliths have undergone 1–10% fractional melting. The spinel peridotites from Baegryeong Island have undergone cryptic mantle metasomatism subsequent to melt extraction. Metasomatic agent of enriched spinel peridotite xenoliths was carbonatite melt.     

Homologous temperature of olivine: Implications for creep of the upper mantle and fabric transitions in olivine

The homologues temperature of a crystalline material is defined as T/T_m, where T is temperature and T_m is the melting(solidus) temperature in Kelvin. It has been widely used to compare the creep strength of crystalline materials. The melting temperature of olivine system,(Mg,Fe)_2SiO_4, decreases with increasing iron content and water content, and increases with confining pressure. At high pressure, phase transition will lead to a sharp change in the melting curve of olivine. After calibrating previous melting experiments on fayalite(Fe_2SiO_4), the triple point of fayalite-Fe_2SiO_4 spinel-liquid is determined to be at 6.4 GPa and 1793 K. Using the generalized means, the solidus and liquidus of dry olivine are described as a function of iron content and pressure up to 6.4 GPa. The change of T/T_m of olivine with depth allows us to compare the strength of the upper mantle with different thermal states and olivine composition. The transition from semi-brittle to ductile deformation in the upper mantle occurs at a depth where T/T_m of olivine equals 0.5. The lithospheric mantle beneath cratons shows much smaller T/T_m of olivine than orogens and extensional basins until the lithosphere-asthenosphere boundary where T/T_m 0.66, suggesting a stronger lithosphere beneath cratons. In addition, T/T_m is used to analyze deformation experiments on olivine. The results indicate that the effect of water on fabric transitions in olivine is closely related with pressure. The hydrogen-weakening effect and its relationship with T/T_m of olivine need further investigation. Below 6.4 GPa(200 km), T/T_m of olivine controls the transition of dislocation glide from [100] slip to [001] slip. Under the strain rate of 10~(-12)–10~(-15) s~(-1) and low stress in the upper mantle, the [100](010) slip system(A-type fabric) becomes dominant when T/T_m 0.55–0.60. When T/T_m 0.55–0.60, [001] slip is easier and low T/T_m favors the operation of [001](100) slip system(C-type fabric). This is consistent with the widely observed A-type olivine fabric in naturally deformed peridotites, and the C-type olivine fabric in peridotites that experienced deep subduction in ultrahigh-pressure metamorphic terranes. However, the B-type fabric will develop under high stress and relatively low T/T_m. Therefore, the homologues temperature of olivine established a bridge to extrapolate deformation experiments to rheology of the upper mantle. Seismic anisotropy of the upper mantle beneath cratons should be simulated using a four-layer model with the relic A-type fabric in the upper lithospheric mantle, the B-type fabric in the middle layer, the newly formed A- or B-type fabric near the lithosphere-asthenosphere boundary, and the asthenosphere dominated by diffusion creep below the Lehmann discontinuity. Knowledge about transition mechanisms of olivine fabrics is critical for tracing the water distribution and mantle flow from seismic anisotropy.     

Structure of the upper mantle in Asia from phase and group velocities of Rayleigh waves

Dispersion curves of phase velocities of Rayleigh waves are determined by the method of frequency-time analysis in a range of periods of 10–200 s from data of 43 interstation traces in Central Asia. Because the joint use of phase and group velocities significantly decreases the uncertainty in the determination of S wave velocity structures, the same traces were used for calculating group velocities from tomographic reconstructions obtained in [Yanovskaya and Kozhevnikov, 2003, 2006] and determining average velocity structures along these traces. The velocity structures were calculated by the Monte Carlo and linear inversion methods, which gave consistent results. Using velocity values obtained at fixed depths by the 2-D tomography method, lateral variations in velocities at these depths were estimated, which allowed us to construct smoothed vertical velocity structures at some points in the region. The resulting structures were used as initial approximations for constructing local velocity structures solely from previously obtained local dispersion curves of group velocities in the area (32°–56°N, 80°–120°E). Based on these structures, we mapped the lateral distribution of velocity variations at upper mantle depths of 75–400 km and along three vertical profiles. The inferred velocity variations are in good agreement with data on the tectonics of the region.     

Experimental evidence for the exsolution of cratonic peridotite from high-temperature harzburgite

Phase equilibrium experiments were performed on typical ‘oceanic’ and ‘cratonic’ peridotite compositions and a Ca, Al-rich orthopyroxene composition, to test the proposal that garnet lherzolites exsolved from high-temperature harzburgites, and to further our understanding of the origin of ancient cratonic lithospheres. ‘Oceanic’ peridotites crystallize a garnet harzburgite assemblage at pressures above 5 GPa in the temperature range 1450–1600°C, but at 5 GPa and temperatures less than 1450°C, crystallize clinopyroxene to become true lherzolites. ‘Cratonic’ peridotites crystallize a garnet harzburgite assemblage at pressures above 5 GPa in the temperature range 1300–1600°C. Garnet-free harzburgite crystallizes from both ‘cratonic’ and ‘oceanic’ peridotite at temperatures above 1450°C and pressures below 4.5–5 GPa. Phase relations for the high Ca, Al-rich orthopyroxene composition essentially mirror those for ‘oceanic’ peridotite.The complete solution of garnet and clinopyroxene into orthopyroxene observed in all three starting compositions at temperatures near or above the mantle solidus at pressures less than 6 GPa supports the hypothesis that garnet lherzolite could have exsolved from harzburgite. The inferred cooling path for the original high-temperature harzburgite protoliths of garnet lherzolites differs depending on bulk composition. The precursor harzburgite protoliths of garnet lherzolites and harzburgites with ‘cratonic’ bulk compositions apparently experienced simple isobaric cooling from formation temperatures near the peridotite solidus to those at which most of these peridotites were sampled in the mantle (< 1200°C). The cooling histories for harzburgite protoliths of sheared garnet lherzolites with ‘oceanic’ compositional affinity are speculated to have involved convective circulation of mantle material to depths deeper than those at which it was originally formed.Phase equilibria and compositional relationships for orthopyroxenes produced in phase equilibrium experiments on peridotite and komatiite are consistent with an origin for ‘cratonic’ peridotite as a residue of Archean komatiite extraction, which has since cooled and exsolved clinopyroxene and garnet to become the common low-temperature, coarse-grained peridotite thought to comprise the bulk of the mantle lithosphere beneath the Archean Kaapvaal craton.     

Velocity structure of the lithosphere on the 2003 Mongolian-Baikal transect from <Emphasis Type="Italic">SV</Emphasis> waves

The S wave velocity distribution in the Earth’s crust and the first two hundred kilometers of the upper mantle is inferred from data of a seismological linear network including 18 broadband stations installed in the framework of the international teleseismic experiment carried out in 2003 in the south of Siberia and in Mongolia. Models were constructed by using P-to-S received function inversion beneath each station. Vertical cross sections of S wave velocities from the surface to depths of 65 and 270 km covering the entire 1000-km profile are constructed by the linear spline interpolation of individual velocity models. The vertical sections are also represented as anomalies relative to the standard velocity model. The most intense low velocity anomalies (from ?3 to ?6%) in the crust and upper mantle are identified beneath the Sayan, Khamar-Daban, and Khangai highlands and the Djida fold zone and agree both with other geophysical data and with the distribution of Late Cenozoic volcanic fields. The results of this work suggest that the activation of Mongolian-Siberian highlands is largely connected with uplift of the asthenosphere to the base of the crust.     

Effects of chemical,mineralogical, and temperature variations on the density and bulk sound velocity of the bottom lower mantle

In this study, we have modeled the density(ρ) and bulk sound velocity(VΦ) profiles of the bottom lower mantle using the experimental thermal equation of state(EoS) parameters of lower-mantle minerals, including bridgmanite, ferropericlase,CaSiO3-perovskite, and post-perovskite. We re-evaluated the literature pressure-volume-temperature relationships of these minerals using a self-consistent pressure scale in order to avoid the long-standing pressure scale problem and to provide more reliable constraints on the thermal EoS parameters. With the obtained thermal EoS parameters, we have constructed the ρ and VΦ profiles of the bottom lower mantle in different composition, mineralogy, and temperature models. Our modelling results show that the variations of chemistry, mineralogy, and temperature have different seismic signatures from each other. The Fe and Al enrichment at the bottom lower mantle can cause an increase in ρ but greatly lower VΦ. A change in mineralogy needs to be considered with the lateral variation in temperature. The cold slabs will be shown as denser regions compared to the normal mantle because of the combined effect of a lower temperature and the presence of a denser post-perovskite at a shallower depth,whereas the hot regions will have a 1–2% lower ρ than the normal mantle. VΦ of both cold slabs and hot regions will be lower than the normal mantle when bridgmanite is the dominant phase in the normal mantle, yet they will be greater once bridgmanite transforms into post-perovskite in the normal mantle. Our modeling also shows that the presence of a(Fe, Al)-enriched bridgmanite thermal pile above the core-mantle boundary will exhibit a seismic signature of enhanced ρ and VΦ, but a reduced VS,which is consistent with the observed seismic anomalies in the large-low-shear-velocity-provinces(LLSVPs). The existence of such a(Fe, Al)-enriched bridgmanite thermal pile thus can help to understand the origin of the LLSVPs. These results provide new insights for the chemical and structure of the deepest lower mantle.     

Velocity properties of the lithosphere in the ocean-continent transition zone in the Kamchatka region from seismic tomography data

Arrival times of P and S waves from local earthquakes in the Kamchatka area of the Kurile-Kamchatka Island Arc are used for calculating a spatial model of the elastic wave velocity distribution to a depth of 200 km. The lithosphere is shown to be strongly stratified in its velocity properties and laterally heterogeneous within the mantle wedge and seismic focal zone. A lower velocity layer (an asthenospheric wedge) is identified at depths of 70–130 km beneath the Eastern Kamchatka volcanic belt. The morphology of the Moho interface and the velocity properties of the crust are studied. The main tectonic structures of the region are shown to be closely interrelated with deep velocity heterogeneities. Regular patterns in the statistics of the earthquakes are analyzed in relation to variations in the elastic wave velocities in the focal layer. A mechanism of lithospheric block displacements along weakened zones of the lower crust and upper mantle is proposed.     

Generation of komatiite magma and gravitational differentiation in the deep upper mantle

The densities of silicate liquids with basic, picritic, and ultrabasic compositions have been estimated from the melting curves of minerals at high pressures. Silicate liquids generated by partial melting of the upper mantle are denser than olivine and pyroxenes at pressures higher than 70 kbar, and garnet is the only phase which is denser than the liquid at pressures from 70 kbar to at least 170 kbar. In this pressure range, garnet and some fraction of liquid separate from ascending partially molten diapirs. It is therefore suggested that aluminium-depleted komatiite with a high Ca/OAl₂O₃ ratio may be derived from diapirs which originated in the deep upper mantle at pressures from 70 kbar to at least 140 kbar (200–400 km in depth), where selective separation of pyropic garnet occurs effectively. On the other hand, aluminium-undepleted komatiite is probably derived from diapirs originating at shallower depths (< 200 km). Enrichment of pyropic garnet is expected at depths greater than 200 km by selective separation of garnet from ascending diapirs. The 200-km discontinuity in the seismic wave velocity profile may be explained by a relatively high concentration of pyropic garnet at depths greater than 200 km.     

A petrogenetic study of the garnet pyroxenite enclaves in spinel peridotite,North Dabieshan,China

Recently, garnet pyroxenite enclaves within peridotites occurring near Raobazhai, Huoshan County, have been discovered. The garnet pyroxenite is small pods, decimeters in size, enclosed within intensively serpentinized peridotites. Major mineral components comprise: garnet (Prp_25–35), sodium augite (Jd_10–25) with a small amount of ilmenite. There are two stages of retrometamorphism: the retrogressive granulite facies mineral assemblage is superimposed by that of amphibolite facies. The host rocks of the garnet pyroxenite are spinel peridotites, including spinel harzburgite and lherzolite. Due to intensive serpentinitization, only 5%–40% of the relic olivine (Fo_92–93) are preserved. The orthopyroxenes are Mg-rich (En_87–93) with bending of cleavages and granulation at their margins showing intracrystalline plasticity. On the basis of garnet-clinopyroxene Fe−Mg exchange equilibrium geothermometry proposed by Ellis & Green (1979) and Krogh (1988)K _D=4.06–5.28;T=793–919°C,P=1.5 GPa are estimated for the garnet pyroxenite. It is inferred that the peridotites are mantle rocks about 60 km in depth. During the exhumation of the orogenic belt, it was tectonically emplaced into the lower crust in the solid state and then uplifted to the shallow depth. Obviously, this kind of garnet pyroxenite must be petrogenetically related to its host rock. The REE distribution pattern and the Ni−Co−Sc diagram reveal that they are chemically equivalent to the basaltic melt and ultramafic residua respectively derived from partial melting of mantle rocks.