Thermal conditions of the central Sinai Microplate inferred from new surface heat-flow values and continuous borehole temperature logging in central and southern Israel

This paper reports ten new surface heat-flow density (qs) values for central and southern Israel (central Sinai Microplate), whose crystalline crust and lithosphere formed as part of the Neoproterozoic Arabian-Nubian Shield. Heat flow was calculated in Mesozoic sediments using the classical approach of heat-flow determination by implementing in the analysis high-precision continuous temperature logs obtained in air- and/or water-filled boreholes. Thermal conductivity (TC) measured for a large suite of rock samples of lithotypes making up the sequence was assigned to temperature gradients in intervals for which the lithology was known. The heat-flow values obtained for different depth intervals in a borehole as well as the average values for the individual borehole locations cover a narrow range, attesting heat-conduction conditions. A steady-state thermal model along an E–W crustal cross section through the area shows that the observed systematic spatial distribution of the qs values, which range between 50 and 62 mW m⁻², can primarily be explained by variations in the thickness of the upper crust and in the ratio between sedimentary and crystalline rocks therein. Given the time lapse of thermal heat transfer through the lithosphere, the qs data monitor the crustal thermal conditions prior to rift- and plume-related lithospheric thermal perturbations that have started in the larger area ca. 30 Ma ago. Observed and modeled qs display the best fit for a pre-Oligocene lithosphere–asthenosphere boundary (LAB) at ∼150 km, which would be at the upper end of LAB depths determined from stable areas of the Arabian Shield (150–120 km) not affected by the young, deep-seated thermal processes that have caused a further uprise of the LAB. Our data imply or predict that the surface heat flow of the Sinai Microplate generally tends to increase along N–S and W–E traverses, from ∼45–50 mW m⁻² to ∼55–60 mW m⁻². Surface heat flows on the order of 55–60 mW m⁻² may be common in the northern Arabian Shield, where it exhibits typical lithosphere structure and composition and is unaffected by young heating processes, compared to values of ≤45 mW m⁻² recently determined in the southern Arabian Plate for the Arabian Platform.     

Ancient Os isotope signatures from the Ontong Java Plateau lithosphere: Tracing lithospheric accretion history

In order to better understand the nature and formation of oceanic lithosphere beneath the Early Cretaceous Ontong Java Plateau, Re–Os isotopes have been analysed in a suite of peridotite xenoliths from Malaita, Solomon Islands. Geological, thermobarometric and petrological evidence from previous studies reveal that the xenoliths represent virtually the entire thickness of the southern part of subplateau lithospheric mantle (< 120 km). This study demonstrates that vertical Os isotopic variations correlate with compositional variations in a stratified lithosphere. The shallowest plateau lithosphere (< 85 km) is dominated by fertile lherzolites showing a restricted range of ¹⁸⁷Os/¹⁸⁸Os (0.1222 to 0.1288), consistent with an origin from ~ 160 Ma Pacific lithosphere. In contrast, the basal section of subplateau lithospheric mantle (~ 95–120 km) is enriched in refractory harzburgites with highly unradiogenic ¹⁸⁷Os/¹⁸⁸Os ratios ranging from 0.1152 to 0.1196, which yield Proterozoic model ages of 0.9–1.7 Ga. Although the whole range of Os isotope compositions of Malaita peridotites is within the variations seen in modern abyssal peridotites, the contrasting isotopic compositions of shallow and deep plateau lithosphere suggest their derivation from different mantle reservoirs. We propose that the subplateau lithosphere forms a genetically unrelated two-layered structure, comprising shallower, typical oceanic lithosphere underpinned by deeper impinged material, which included a component of recycled Proterozoic lithosphere. The impingement of residual but chemically heterogeneous mantle, mechanically coupled to the recently formed, thin lithosphere, may have a bearing on the anomalous initial uplift and late subsidence history of the seismically anomalous plateau root.     

Effective elastic thicknesses of the lithosphere in the Central Iberian Peninsula from heat flow: Implications for the rheology of the continental lithospheric mantle

The traditional view of the rheology of the continental lithosphere, sometimes known as the “jelly sandwich model”, consists of a strong upper crust, a weak lower crust, and a strong upper lithospheric mantle. Some authors argue, however, that the lithospheric mantle is weak and contributes little to the total strength and the effective elastic thickness of the lithosphere; this weakness is claimed to be due to the mantle being wet or subjected to temperatures higher than usually believed. This paper uses the relationship between rheology of the lithosphere and heat flow to calculate theoretical effective elastic thicknesses for three regions of the central Iberian Peninsula (the Duero Basin, the Spanish Central System and the Tajo Basin), taking into account the contribution of the crust and the lithospheric mantle, for dry and wet rheologies. We found that a wet peridotite rheology for the lithospheric mantle is generally consistent with independent (based on Bouguer coherence or flexural modeling) estimates of the effective elastic thickness for the study area, whereas a dry peridotite rheology cannot be reconciled with them. Moreover, the contribution of the mantle to the bending moment of the lithosphere, and therefore to both the effective elastic thickness and the total strength of the lithosphere, is important, and it may even be the dominant contribution. Therefore, the jelly sandwich model may be considered valid for the central Iberian Peninsula.     

Three-dimensional numerical thermal and rheological modelling in the central Fennoscandian Shield

A three-dimensional model for the central Fennoscandian Shield was constructed for analysing the thermal, the rheological and the structural conditions in the lithosphere. The mesh covers a rectangular area in the southern Finland with horizontal dimensions of 500 km × 400 km and a depth extent of 100 km. Structural boundaries are derived from the several deep seismic soundings carried out in the area. Constructed model is first used in the calculation of the thermal and the rheological models and secondly in analysing the stress and the deformational conditions with the obtained rheology. Thermal and structural models are solved with the finite element method. The calculated surface HFD is between 40 and 48 mW m⁻² in the Proterozoic southern part and below 40 mW m⁻² in the older and northern Archaean part of the model. The calculated rheological strength shows a layered structure with two individual rheologically weak layers in the crust and strong layer in the upper part of the lower crust. The minimum brittle–ductile transition (BDT) depth is around 10 km in the southern part of the model while in the north and north-eastern parts the BDT depth is around 45–50 km. Comparison with the focal depth data shows that as most of the earthquakes occur no deeper than the depth of 10 km are they located in the brittle regime. Resulting stress conditions and possible regions of deformation after the model is subjected to pressure of 50 MPa reveals that the stress field is quite uniformly distributed in different crustal layers and that the elastic parameters control more the state of the stress than the applied rheological structure. In the upper crust, the stress intensity has values between 42 and 45 MPa whereas in the middle crust the values are around 50 MPa. Comparison of the 3-D model with earlier 2-D models shows that some differences in the results are to be expected.     

Lithospheric contributions to high-MgO basanites from the Cumbre Vieja Volcano,La Palma,Canary Islands and evidence for temporal variation in plume influence

New geochemical and isotopic data are presented from the oldest part of the Cumbre Vieja volcano, La Palma (Canary Islands), located near the assumed emergence of the Canary mantle plume. The volcanics comprise a suite dominated by basanite flows with subordinate amounts of phono-tephrite, tephri-phonolite and phonolite flows and intrusives. Two compositionally different basanite groups have been identified, both with HIMU (high-μ)-type incompatible trace element characteristics: Primitive high-MgO basanites (10.7–12.1% MgO), found only at the base of a stratigraphic profile near Fuencaliente on the south coast, and intermediate-MgO basanites (6.0–7.3% MgO), exposed in the upper part of the profile and widespread on the east coast of La Palma. The high-MgO basanites are interpreted as near-primary mantle melts (primary composition 14–15% MgO) derived by progressive melting (2.9% to 4.5%) of a common lithospheric mantle source. Model calculations indicate that it is not possible to generate the intermediate-MgO basanites from the high-MgO group by crystal fractionation of observed phenocrysts. Relative to intermediate-MgO basanites, the high-MgO flows have lower concentrations of LIL and HFS elements, except for Ti, which is markedly enriched in the primitive rocks (3.7–4.7% TiO₂ vs 3.4–3.9% TiO₂). Fuencaliente volcanics display limited temporal isotopic variations suggested to be a result of mixing of melts originating from the rising plume and the metazomatized lithospheric mantle. ⁸⁷Sr / ⁸⁶Sr and ¹⁴³Nd / ¹⁴⁴Nd ratios range 0.70305–0.70311 and 0.51285–0.51291, respectively, while the corresponding ranges in Pb-isotope ratios are ²⁰⁶Pb / ²⁰⁴Pb = 19.46–19.64, ²⁰⁷Pb / ²⁰⁴Pb = 15.55–15.61, and ²⁰⁸Pb / ²⁰⁴Pb = 39.16–39.53. The overall variation of the Cumbre Vieja isotopic data can be accounted for by mixtures of three mantle components in the proportions 72–79% plume source (LVC = low velocity component), 9–16% depleted mantle (DM) and up to 12% enriched mantle (EMI). Negative Δ7 / 4 Pb (− 0.6 to − 5.4) in the Cumbre Vieja volcanics suggest derivation from a young HIMU mantle source. The relative abundance of plume source material increase in younger rocks in the Fuencaliente section, suggesting waning plume–lithosphere interaction during the emplacement of this part of the Cumbre Vieja volcano. The high-MgO volcanics define regular and systematic geochemical trends, interpreted as partial melting trends, when plotted against abundances of highly incompatible elements (P, Ce). Evaluation of minor and trace element variation in consecutive melts suggests control by residual amphibole, phlogopite, garnet and a Ti-bearing phase, possibly ilmenite. The melting mode changed gradually, allowing increasing input from residual phlogopite during partial melting. The residual mineralogy constrains the source region of the high-MgO basanites to the lowermost oceanic lithospheric mantle, presumably around 100 km depths.     

Deccan plume,lithosphere rifting,and volcanism in Kutch,India

Kutch (northwest India) experienced lithospheric thinning due to rifting and tholeiitic and alkalic volcanism related to the Deccan Traps K/T boundary event. Alkalic lavas, containing mantle xenoliths, form plug-like bodies that are aligned along broadly east–west rift faults. The mantle xenoliths are dominantly spinel wehrlite with fewer spinel lherzolite. Wehrlites are inferred to have formed by reaction between transient carbonatite melts and lherzolite forming the lithosphere. The alkalic lavas are primitive (Mg# = 64–72) relative to the tholeiites (Mg# = 38–54), and are enriched in incompatible trace elements. Isotope and trace element compositions of the tholeiites are similar to what are believed to be the crustally contaminated Deccan tholeiites from elsewhere in India. In terms of Hf, Nd, Sr, and Pb isotope ratios, all except two alkalic basalts plot in a tight cluster that largely overlap the Indian Ridge basalts and only slightly overlap the field of Reunion lavas. This suggests that the alkalic magmas came largely from the asthenosphere mixed with Reunion-like source that welled up beneath the rifted lithosphere. The two alkalic outliers have an affinity toward Group I kimberlites and may have come from an old enriched (metasomatized) asthenosphere. We present a new model for the metasomatism and rifting of the Kutch lithosphere, and magma generation from a CO₂-rich lherzolite mantle. In this model the earliest melts are carbonatite, which locally metasomatized the lithosphere. Further partial melting of CO₂-rich lherzolite at about 2–2.5 GPa from a mixed source of asthenosphere and Reunion-like plume material produced the alkalic melts. Such melts ascended along deep lithospheric rift faults, while devolatilizing and exploding their way up through the lithosphere. Tholeiites may have been generated from the main plume head further south of Kutch.     

Mantle dynamics of the Paleoproterozoic North China Craton: A perspective based on seismic tomography

We investigate the mantle dynamics beneath the North China Craton (NCC) and surrounding regions based on a synthesis of recent P-wave mantle tomographic data down to depths of 600–800 km and their correlation with the surface geological features, with particular reference to the Paleoproterozoic tectonic events associated with the incorporation of the NCC within the Columbia supercontinent amalgam. From the tomographic images, we identify a hot corridor in the mantle transition zone beneath the central region of the Western Block of the NCC sandwiched between two cold corridors. This scenario is similar to the donut-shaped high-velocity anomaly surrounding a region of low-velocity anomaly in the lowermost mantle under the Pacific and suggests that the cold regions might represent slab graveyards which provide the fuel for the plumes rising from the center. A tomographic transect along the collisional suture of the NCC with the Columbia supercontinent, covering the Yinshan-Ordos Blocks in the Western Block through the Central Orogenic Belt and into the Eastern Block of the NCC reveals a ca. 250 km thick lithospheric keel below the Ordos Block defined by a prominent high-velocity anomaly. We identify slab break-off and asthenospheric upwelling in this region and suggest that this process probably initiated the thermal and material erosion of the tectosphere beneath the Eastern Block from the Paleoproterozoic, which was further intensified during the Mesozoic when a substantial part of the sub-continental mantle lithosphere was lost. We visualize heat input from asthenosphere and interaction between asthenosphere and overlying carbonated tectosphere releasing CO₂-rich fluids for the preservation of ultra-high temperature (ca. 1000 °C) metamorphic rocks enriched in CO₂ as well as high-pressure mafic granulites as a paired suite in this region. We also identify a hot swell of the asthenosphere rooted to more than 200 km depth and reaching up to the shallow mantle in the tomographic section along 35°N latitude at a depth of 800 km. This zone represents a cross-section through the southern part of the NCC. The surface distribution of Paleoproterozoic Xiong’er lavas and mafic dykes in this region would indicate that this region might have evidenced similar upwellings in the past. Our study has important implications in understanding the evolution of the NCC and suggests that the extensive modification of the mantle architecture and lithospheric structure beneath one of the fundamental Precambrian nuclei of Asia had a prolonged history probably dating from the Paleoproterozoic suturing of the NCC within the Columbia supercontinent amalgam.     

Depth constraints on azimuthal anisotropy in the Great Basin from Rayleigh-wave phase velocity maps

We present fundamental-mode Rayleigh-wave azimuthally anisotropic phase velocity maps obtained for the Great Basin region at periods between 16 s and 102 s. These maps offer the first depth constraints on the origin of the semi-circular shear-wave splitting pattern observed in central Nevada, around a weak azimuthal anisotropy zone. A variety of explanations have been proposed to explain this signal, including an upwelling, toroidal mantle flow around a slab, lithospheric drip, and a megadetachment, but no consensus has been reached. Our phase velocity study helps constrain the three-dimensional anisotropic structure of the upper mantle in this region and contributes to a better understanding of the deformation mechanisms taking place beneath the western United States. The dispersion measurements were made using data from the USArray Transportable Array. At periods of 16 s and 18 s, which mostly sample the crust, we find a region of low anisotropy in central Nevada coinciding with locally reduced phase velocities, and surrounded by a semi-circular pattern of fast seismic directions. Away from central Nevada the fast directions are ～ N–S in the eastern Great Basin, NW–SE in the Walker Lane region, and they transition from E–W to N–S in the northwestern Great Basin. Our short-period phase velocity maps, combined with recent crustal receiver function results, are consistent with the presence of a semi-circular anisotropy signal in the lithosphere in the vicinity of a locally thick crust. At longer periods (28–102 s), which sample the uppermost mantle, isotropic phase velocities are significantly reduced across the study region, and fast directions are more uniform with an ～ E–W fast axis. The transition in phase velocities and anisotropy can be attributed to the lithosphere–asthenosphere boundary at depths of ～ 60 km. We interpret the fast seismic directions observed at longer periods in terms of present-day asthenospheric flow-driven deformation, possibly related to a combination of Juan de Fuca slab rollback and eastward-driven mantle flow from the Pacific asthenosphere. Our results also provide context to regional SKS splitting observations. We find that our short-period phase velocity anisotropy can only explain ～ 30% of the SKS splitting times, despite similar patterns in fast directions. This implies that the origin of the regional shear-wave splitting signal is complex and must also have a significant sublithospheric component.     

Bending-related faulting of the incoming oceanic plate and its effect on lithospheric hydration and seismicity: A passive and active seismological study offshore Maule,Chile

We have studied the dependency between incoming plate structure, bending-related faulting, lithospheric hydration, and outer rise seismic activity offshore Maule, Chile. We derived a 2D Poisson's ratio distribution from P- and S-wave seismic wide angle data collected in the trench-outer rise. High values of Poisson's ratio in the uppermost mantle suggest that the oceanic lithosphere is highly hydrated due to the water infiltration through bending-related normal faults outcropping at the seafloor. This process is presumably facilitated by the presence of a seamount in the area. We conclude that water infiltrates deep into the lithosphere, when it approaches the Chile trench, producing a reduction of crustal and upper mantle velocities, supporting serpentinization of the upper mantle. Further, we observed a mantle Vp anisotropy of 8%, with the fast velocity axis running normal to the abyssal hill fabric and hence in spreading direction, indicating that outer rise processes have yet not affected anisotropy.The first weeks following the megatrust Mw = 8.8 Maule earthquake in 2010 were characterized by a sudden increase of the outer rise seismic activity, located between 34° S and 35°30′ S. We concluded that this phenomenon is a result of an intensification of the water infiltration process in the outer rise, presumably triggered by the main shock, whose epicenter was located some 100 km to the south east of the cluster.     

Intraplate seismicity in the western Bohemian Massif (central Europe): A possible correlation with a paleoplate junction

Locations of the Eger Rift, Cheb Basin, Quaternary volcanoes, crustal earthquake swarms and exhalation centers of CO₂ and ³He of mantle origin correlate with the tectonic fabric of the mantle lithosphere modelled from seismic anisotropy. We suggest that positions of the seismic and volcanic phenomena, as well as of the Cenozoic sedimentary basins, correlate with a “triple junction” of three mantle lithospheres distinguished by different orientations of their tectonic fabric consistent within each unit. The three mantle domains most probably belong to the originally separated microcontinents – the Saxothuringian, Teplá-Barrandian and Moldanubian – assembled during the Variscan orogeny. Cenozoic extension reactivated the junction and locally thinned the crust and mantle lithosphere. The rigid part of the crust, characterized by the presence of earthquake foci, decoupled near the junction from the mantle probably during the Variscan. The boundaries (transitions) of three mantle domains provided open pathways for Quaternary volcanism and the ascent of ³He- and CO₂-rich fluids released from the asthenosphere. The deepest earthquakes, interpreted as an upper limit of the brittle–ductile transition in the crust, are shallower above the junction of the mantle blocks (at about 12 km) than above the more stable Saxothuringian mantle lithosphere (at about 20 km), probably due to a higher heat flow and presence of fluids.     

On the Vp/Vs–Mg# correlation in mantle peridotites: Implications for the identification of thermal and compositional anomalies in the upper mantle

We use thermodynamically self-consistent and hybrid methods to analyze the correlation of important physical parameters (e.g. bulk density, elastic moduli) with bulk Mg# and modal composition in mantle peridotites at upper mantle conditions. Temperature (anharmonic and anelastic), pressure and compositional derivatives for all these parameters are evaluated. The results show that the widely used correlations between Vp/Vs and Mg# in peridotites are strictly valid only for garnet-bearing assemblages at temperatures < 900 °C. The correlation breaks down when: i) spinel is the stable Al-rich phase in the assemblage and ii) when anelastic attenuation of seismic velocities becomes important (T ? 900 °C). This implies that the range of applicability of published Vp/Vs–Mg# correlations for the upper mantle is limited to a depth interval between the spinel–garnet phase transition and the 900 °C isotherm. We use numerical simulations to show that this depth interval is virtually nonexistent in lithospheres thinner than ～ 140 km and can comprise up to ～ 50% of the lithospheric mantle in thick (> 220 km) lithospheric domains. In addition, we show that for most of the upper mantle the expected Δ(Vp/Vs) values associated with compositional variations are smaller than the resolution limit of current seismological methods. All these considerations suggest that the Vp/Vs ratio is not a reliable measure of compositional variations and that for large parts of the upper mantle compositional anomalies cannot be separated from thermal anomalies on the basis of seismological studies only. We further confirm that the only reliable indicator of compositional anomalies in a peridotitic mantle is the ratio of density to shear wave velocities (ρ/Vs). Our results demonstrate that geophysical–petrological models (forward or inverse) that model these two fields (i.e. density and Vs) self-consistently within a robust thermodynamic framework are necessary for characterizing the small-scale thermal and compositional structure of the lithosphere and sublithospheric upper mantle.     

Atmospheric noble gas signatures in deep Michigan Basin brines as indicators of a past thermal event

Atmospheric noble gases (e.g., ²²Ne, ³⁶Ar, ⁸⁴Kr, ¹³⁰Xe) in crustal fluids are only sensitive to subsurface physical processes. In particular, depletion of atmospheric noble gases in groundwater due to boiling and steam separation is indicative of the occurrence of a thermal event and can thus be used to trace the thermal history of stable tectonic regions. We present noble gas concentrations of 38 deep brines (~ 0.5–3.6 km) from the Michigan Basin. The atmospheric noble gas component shows a strong depletion pattern with respect to air saturated water. Depletion of lighter gases (²²Ne and ³⁶Ar) is stronger compared to the heavier ones (⁸⁴Kr and ¹³⁰Xe). To understand the mechanisms responsible for this overall atmospheric noble gas depletion, phase interaction models were tested. We show that this atmospheric noble gas depletion pattern is best explained by a model involving subsurface boiling and steam separation, and thus, consistent with the occurrence of a past thermal event of mantle origin as previously indicated by both high ⁴He/heat flux ratios and the presence of primordial mantle He and Ne signatures in the basin. Such a conceptual model is also consistent with the presence of past elevated temperatures in the Michigan Basin (e.g., ~ 80–260 °C) at shallow depths as suggested by previous thermal studies in the basin. We suggest that recent reactivation of the ancient mid-continent rift system underneath the Michigan Basin is likely responsible for the release of both heat and mantle noble gases into the basin via deep-seated faults and fracture zones. Relative enrichment of atmospheric Kr and Xe with respect to Ar is also observed, and is interpreted as reflecting the addition of sedimentary Kr and Xe from associated hydrocarbons, following the hydrothermal event. This study pioneers the use of atmospheric noble gases in subsurface fluids to trace the thermal history of stable tectonic regions.     

Geochemistry of a new enriched mantle type locality in the northern hemisphere: Implications for the origin of the EM-I source

Late Cretaceous (66.2 ± 0.5 Ma amphibole and 66.7 ± 0.2 Ma phlogopite ⁴⁰Ar/³⁹Ar ages) nephelinitic volcanic rocks from Godzilla Seamount in the eastern North Atlantic (34°N latitude) have trace element and Sr–Nd–Pb–Hf-isotope compositions similar to the Enriched Mantle I (EM-I) endmember, except for their low ²⁰⁷Pb/²⁰⁴Pb relative to ²⁰⁶Pb/²⁰⁴Pb ratios (²⁰⁶Pb/²⁰⁴Pb_in = 17.7, ²⁰⁷Pb/²⁰⁴Pb_in = 15.34) plotting below the Northern Hemisphere Reference Line on the uranogenic Pb isotope diagram. O isotope data on amphibole separates are mantle-like (δ¹⁸O = 5.6–5.8‰). Age and location of the isolated Godzilla Seamount, however, preclude it from being derived from the Madeira or Canary hotspots, making a lower-mantle origin unlikely. Therefore we propose derivation from a shallow (lithospheric/asthenospheric) melting anomaly. As observed in mid-ocean-ridge and ocean-island basalts, there is a systematic decrease of ²⁰⁷Pb/²⁰⁴Pb ratios (and Δ7/4) in the individual EM-I endmember type localities towards northern latitudes with Godzilla lying on the extension of this trend. This trend is mirrored in ultra-potassic volcanic rocks such as lamproites and kimberlites, which reflect the composition of enriched subcontinental lithospheric mantle. Therefore, a global pattern in ²⁰⁷Pb/²⁰⁴Pb ratios and Δ7/4 is suggested. The geochemical composition of EM-I endmember type localities, including Godzilla lavas, and the enriched (DUPAL) anomaly in the southern hemisphere could reflect derivation from ancient, metasomatized subcontinental lithospheric mantle. We propose a two-stage model to explain the trace element and isotopic composition of the EM-I mantle endmember localities worldwide: 1) during the early history of the Earth, subcontinental lithosphere was metasomatized by melts from subducted slabs along convergent margins generating high μ (²³⁸U/²⁰⁴Pb) sources, and 2) as the Earth cooled, hydrous fluids replaced hydrous melts as the main slab component metasomatizing the subcontinental lithospheric mantle (generating EM-I sources with lower μ). In accordance with this model, the global variations in ²⁰⁷Pb/²⁰⁴Pb ratios and Δ7/4 could reflect geographic differences in μ and/or the age at which the transition from stages 1 to 2 took place in the Archaean lithosphere. The model would require a re-definition of the EM-I endmember to low ²⁰⁶Pb/²⁰⁴Pb, high ²⁰⁸Pb/²⁰⁴Pb (positive Δ8/4) but variable ²⁰⁷Pb/²⁰⁴Pb (positive and negative Δ7/4).     

Radiogenic heat production in the lithosphere of Sulu ultrahigh-pressure metamorphic belt

The Chinese Continental Scientific Drilling (CCSD) project is located at the Sulu ultrahigh-pressure metamorphic (UHPM) belt. It offers a unique opportunity for studying the radiogenic heat production of both shallower and deeper rocks. Based on the concentrations of radiogenic elements U, Th and K on 349 samples from main hole of CCSD (CCSD MH), pilot holes and exposures, we determined radiogenic heat productions of all major rock types in the Sulu UHPM belt. Results show the mean values of orthogneiss and paragneiss are respectively 1.65 ± 0.81 and 1.24 ± 0.61 µW m^? 3. Due to different composition and grade of retrogressive metamorphism, the eclogites display significant scatter in radiogenic heat production, ranging from 0.01 to 2.85 µW m^? 3, with a mean of 0.44 ± 0.55 µW m^? 3. The radiogenic heat production in ultramafic rocks also varies within a large range of 0.02 to 1.76 µW m^? 3, and the average turns out to be 0.18 ± 0.31 µW m^? 3. Based on the measurements and crustal petrologic model, the vertical distribution model of heat production in Sulu crust is established. The resulting mean heat production (0.76 µW m^? 3) contributes 24 mW m^? 2 to the surface heat flow. 1-D thermal model indicates that the temperature at the Moho reaches above 750 °C, and the thermal thickness of the lithosphere is ~ 75 km, in good agreement with the geophysical results. The high teat flow (~ 75 mW m^? 2) together with thin lithosphere presents strong support for the extension events during the late Cretaceous and Cenozoic.     

The formation of volcanic centers at the Colorado Plateau as a result of the passage of aqueous fluid through the oceanic lithosphere and the subcontinental mantle: New implications for the planetary water cycle in the western United States

We provide new petrological evidence for the strong influence of water on the formation of the oceanic lithospheric mantle, the subcontinental mantle above, and the continental lithosphere. Our analysis throws new light on the hypothesis that new continental lithosphere was formed by the passage of silicate-rich aqueous fluid through the sub-continental mantle. In order to investigate this hypothesis, we analyzed a representative collection of lherzolite and harzburgite xenoliths from the sample volcano known as “The Thumb”, located in the center of the Colorado Plateau, western United States. The studied sample collection exhibits multi-stage water enrichment processes along point, line and planar defect structures in nominally anhydrous minerals and the subsequent formation of the serpentine polymorph antigorite along grain boundaries and in totally embedded annealed cracks. Planar defect structures act like monomineralic and interphase grain boundaries in the oceanic lithosphere and the subcontinental mantle beneath the North American plate, which was hydrated by the ancient oceanic Farallon plate during the Cenozoic and Mesozoic eras. We used microspectroscopical, petrological, and seismological techniques to confirm multi-stage hydration from a depth of ∼150 km to just below the Moho depth. High-resolution mapping of the water distribution over homogeneous areas and fully embedded point, line and planar defects in olivine crystals of lherzolitic and harzburgitic origin by synchrotron infrared microspectroscopy enabled us to resolve local wet spots and thus reconstruct the hydration process occurring at a depth of ∼150 km (T ≈ 1225 °C). These lherzolites originated from the middle part of the Farallon mantle slab; they were released during the break up of the Farallon mantle slab, caused by the instability of the dipping slab. The background hydration levels in homogeneous olivines reached ∼138 ppm wt H₂O, and the water concentration at the planar defects could reach up to ∼1000 ppm wt H₂O. However, the formation of antigorite in grain boundaries was found to be the primary hydration mechanism for harzburgitic samples originating from the subcontinental mantle (for hydration, T ≈ 600 °C). Additionally, the formation of antigorite in lherzolites could be found in annealed cracks. From these observations, we conclude that hydration induces multi-stage water enrichment of the mantle wedge by a process that is dominated by the growth and movement of ubiquitous cracks, which acts as planar defects. Cracks in the mantle seem to be the an important feature in both the water cycle of the subduction zone and the formation of the continental lithosphere.     

Pliocene climate change of the Southwest Pacific and the impact of ocean gateways

The transition from the early Pliocene “Warmhouse” towards the present “Icehouse” climate and the role of Gateway dynamics are intensively debated. Both, the constrictions of the Central American Seaway and the Indonesian Gateway affected ocean circulation and climate during the Pliocene epoch. Here, we use combined δ¹⁸O and Mg/Ca ratios of planktonic foraminifera (marine protozoa) from surface and subsurface levels to reconstruct the thermal structure and changes in salinities from the Southwest Pacific Deep Sea Drilling Project (DSDP) Site 590B from 6.5 to 2.5 Ma. Our data suggest a gradual cooling of ~ 2 °C and freshening of the sea surface during ~ 4.6–4 Ma with an increased meridional temperature gradient between the West Pacific Warm Pool and the Southwest Pacific when the closing of the Central American Seaway reached a critical threshold. After ~ 3.5 Ma, the restricted Indonesian Gateway might have amplified the East Australian Current, allowing enhanced heat transport towards the Southwest Pacific with reduced meridional temperature gradients when the global climate gradually cooled. At the same time our data suggest a cooling and freshening of Subantarctic Mode Water (SAMW) or/and an increased northward flow of SAMW towards Site 590B, possibly a first step towards the present Antarctic Frontal System.     

Scattering coefficients in the mantle revealed from the seismogram envelope analysis based on the multiple isotropic scattering model

At longer periods, scattered ScS waves sometimes dominate over coda waves at large lapse times. Examining recordings of seismic envelopes at 9 IRIS seismic stations of regional earthquakes with focal depths deeper than 150 km in periods from 1 to 20 s for a wide lapse time range up to 2000 s, we found significant frequency dependence. The coda decay gradient at short periods is steeper than that at longer periods; however, the change of coda gradient associated with the ScS arrival becomes distinct as the period becomes longer. In particular, a clear offset of coda amplitude appears in central Asia for 10 and 15 s period bands. The multiple isotropic scattering process of S-waves in the heterogeneous mantle can be simply simulated by using the Monte Carlo simulation method based on the radiative transfer theory in scattering media. Assuming a two-plane-layer attenuation structure and smoothed velocity model of the PREM, we estimated the average total scattering coefficients of S-waves such as 7.52 × 10^− 4∼1.32 × 10^− 3 km^− 1 and 2.08 × 10^− 4∼6.23 × 10^− 4 km^− 1 at 4 s, and 4.51 × 10^− 4∼7.37 × 10^− 4 km^− 1 and 2.80 × 10^− 5∼2.71 × 10^− 4 km^− 1 at 10 s, for the lithosphere and the upper mantle and for the lower mantle, respectively. Our results indicate that scattering occurs mostly in the lithosphere and the upper mantle and support that medium heterogeneity spreads over the whole mantle though its scattering power is small. Strong scattering occurs beneath central Asia and Papua New Guinea, whereas the scattering beneath Italy and regions of east Russia is much weaker. The numerical calculation enables us to confirm that much stronger scattering than intrinsic attenuation causes the offset behavior with coda decay gradient change after the ScS arrival for 4 and 10 s period bands in some regions.     

Volatile and trace elements in basaltic glasses from Samoa: Implications for water distribution in the mantle

We report volatile (H₂O, CO₂, F, S, Cl) and trace element data for submarine alkalic basalt glasses from the three youngest Samoan volcanoes, Ta'u, Malumalu and Vailulu'u. Most samples are visibly sulfide saturated, so have likely lost some S during fractionation. Cl / K ratios (0.04–0.15) extend to higher values than pristine MORBs, but are suspected to be partly due to source differences since Cl / K roughly varies as a function of ⁸⁷Sr / ⁸⁶Sr. There are no resolvable differences in the relative enrichment of F among sources, and compatibility of F during mantle melting is established to be nearly identical to Nd. Shallow degassing has affected CO₂ in all samples, and H₂O only in the most shallowly erupted samples from Vailulu'u. Absolute water contents are high for Samoa (0.63–1.50 wt.%), but relative enrichment of water compared to equally incompatible trace elements (Ce, La) is low and falls entirely below normal MORB values. H₂O / Ce (58–157) and H₂O / La (120–350) correlate inversely with ⁸⁷Sr / ⁸⁶Sr compositions (0.7045–0.7089). This leads us to believe that, because of very fast diffusion of hydrogen in mantle minerals, recycled lithospheric material with high initial water and trace element content will lose water to the drier ambient mantle during storage within the inner Earth. The net result is the counter-intuitive appearance of greater dehydration with greater mantle enrichment. We expect that subducted slabs will experience a two-stage dehydration history, first within subduction zones and then in the ambient mantle during long-term convective mixing.     

Osmium isotope and highly siderophile element systematics of the lunar crust

Coupled ¹⁸⁷Os/¹⁸⁸Os and highly siderophile element (HSE: Os, Ir, Ru, Pt, Pd, and Re) abundance data are reported for pristine lunar crustal rocks 60025, 62255, 65315 (ferroan anorthosites, FAN) and 76535, 78235, 77215 and a norite clast in 15455 (magnesian-suite rocks, MGS). Osmium isotopes permit more refined discrimination than previously possible of samples that have been contaminated by meteoritic additions and the new results show that some rocks, previously identified as pristine, contain meteorite-derived HSE. Low HSE abundances in FAN and MGS rocks are consistent with derivation from a strongly HSE-depleted lunar mantle. At the time of formation, the lunar floatation crust, represented by FAN, had 1.4 ± 0.3 pg g^? 1 Os, 1.5 ± 0.6 pg g^? 1 Ir, 6.8 ± 2.7 pg g^? 1 Ru, 16 ± 15 pg g^? 1 Pt, 33 ± 30 pg g^? 1 Pd and 0.29 ± 0.10 pg g^? 1 Re (～ 0.00002 × CI) and Re/Os ratios that were modestly elevated (¹⁸⁷Re/¹⁸⁸Os = 0.6 to 1.7) relative to CI chondrites. MGS samples are, on average, characterised by more elevated HSE abundances (～ 0.00007 × CI) compared with FAN. This either reflects contrasting mantle-source HSE characteristics of FAN and MGS rocks, or different mantle–crust HSE fractionation behaviour during production of these lithologies. Previous studies of lunar impact-melt rocks have identified possible elevated Ru and Pd in lunar crustal target rocks. The new results provide no supporting evidence for such enrichments.If maximum estimates for HSE in the lunar mantle are compared with FAN and MGS averages, crust–mantle concentration ratios (D-values) must be ≤ 0.3. Such D-values are broadly similar to those estimated for partitioning between the terrestrial crust and upper mantle, with the notable exception of Re. Given the presumably completely different mode of origin for the primary lunar floatation crust and tertiary terrestrial continental crust, the potential similarities in crust–mantle HSE partitioning for the Earth and Moon are somewhat surprising. Low HSE abundances in the lunar crust, coupled with estimates of HSE concentrations in the lunar mantle implies there may be a ‘missing component’ of late-accreted materials (as much as 95%) to the Moon if the Earth/Moon mass-flux estimates are correct and terrestrial mantle HSE abundances were established by late accretion.     

Thermal history and origin of the Tanzanian Craton from Pb isotope thermochronology of feldspars from lower crustal xenoliths

Common and radiogenic Pb isotopic compositions of plagioclase and anti-perthitic feldspars from granulite-facies lower crustal xenoliths from the Labait Volcano on the eastern margin of the Tanzanian Craton have been measured via laser ablation MC-ICP-MS. Common Pb in plagioclase and a single stage Pb evolution model indicate that the lower crust of the Tanzanian Craton was extracted from mantle having a ²³⁸U/²⁰⁴Pb of 8.1 ± 0.3 and a ²³²Th/²³⁸U of 4.3 ± 0.1 at 2.71 ± 0.09 Ga (all uncertainties are 2σ). Since 2.4 Ga, some orthoclase domains within anti-perthites have evolved with a maximum ²³⁸U/²⁰⁴Pb of 6 and ²³²Th/²³⁸U of 4.3. The spread in Pb isotopic composition in the anti-perthitic feldspars yields single crystal Pb–Pb isochrons of ～ 2.4 Ga, within uncertainty of U–Pb zircon ages from the same sample suite. The Pb isotopic heterogeneities imply that these granulites resided at temperatures < 600 °C in the lower crust of the Tanzanian Craton from ca. 2.4 Ga to the present. In concert with the chemistry of surface samples, mantle xenoliths, and lower crustal xenoliths, our data imply that the cratonic lithosphere in Tanzania formed ca. ～ 2.7 Ga, in a convergent margin setting, and has remained undisturbed since 2.7 Ga.