Sediment loading at the southern Chilean trench and its tectonic implications

Non erosive margins are characterized by heavily sedimented trenches which obscure the morphological expression of the outer rise; a forebulge formed by the bending of the subducting oceanic lithosphere seaward of the trench. Depending on the flexural rigidity (D) of the oceanic lithosphere and the thickness of the trench sedimentary fill, sediment loading can affect the lithospheric downward deflection in the vicinity of the trench and hence the amount of sediment subducted. We used seismic and bathymetric data acquired off south central Chile, from which representative flexural rigidities are estimated and the downward deflection of the oceanic Nazca plate is studied. By flexural modeling we found that efficient sediment subduction preferentially occurs in weak oceanic lithosphere (low D), whereas wide accretionary prisms are usually formed in rigid oceanic lithosphere (high D). In addition, well developed forebulges in strong oceanic plates behaves as barrier to seaward transportation of turbidites, whereas the absence of a forebulge in weak oceanic plates facilitates seaward turbidite transportation for distances >200 km.     

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.     

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.     

The thermal state and strength of the lithosphere in the Spanish Central System and Tajo Basin from crustal heat production and thermal isostasy

In this work we have modeled the thermal structure of the lithosphere of the Spanish Central System and the Tajo Basin, and their implications for lithospheric strength. For this, we have used refined heat-producing elements (HPE) values to obtain new estimates of heat production rates in the Spanish Central System and Tajo Basin areas, which have been used joined to the relation between topography and thermal structure of the lithosphere to calculate the best-fit surface heat flows in the study area. Moreover, we have implemented a temperature-dependent thermal conductivity (appropriate for olivine) for the lithospheric mantle to improve the calculations of temperature profiles in the mantle. The geotherms so obtained, together with the implementation of a new rheological law for the upper lithospheric mantle, have been used to calculate refined estimations of the strength and effective elastic thickness of the lithosphere. We have obtained surface heat flow values of 84 mW m⁻² and ∼82 mW m⁻² for the Spanish Central System and the Tajo Basin, respectively. The thermal state of the lithosphere affects mantle temperatures, and hence may be playing an important role in the uplift and maintenance of the Spanish Central System.     

Miocene to Late Quaternary Patagonian basalts (46–47°S): Geochronometric and geochemical evidence for slab tearing due to active spreading ridge subduction

Miocene to Quaternary large basaltic plateaus occur in the back-arc domain of the Andean chain in Patagonia. They are thought to result from the ascent of subslab asthenospheric magmas through slab windows generated from subducted segments of the South Chile Ridge (SCR). We have investigated three volcanic centres from the Lago General Carrera–Buenos Aires area (46–47°S) located above the inferred position of the slab window corresponding to a segment subducted 6 Ma ago. (1) The Quaternary Río Murta transitional basalts display major, trace elements, and Sr and Nd isotopic features similar to those of oceanic basalts from the SCR and from the Chile Triple Junction near Taitao Peninsula (e.g., (⁸⁷Sr/⁸⁶Sr)_o = 0.70396–0.70346 and εNd = + 5.5 − + 3.0). We consider them as derived from the melting of a Chile Ridge asthenospheric mantle source containing a weak subduction component. (2) The Plio-Quaternary (< 3.3 Ma) post-plateau basanites from Meseta del Lago Buenos Aires (MLBA), Argentina, likely derive from small degrees of melting of OIB-type mantle sources involving the subslab asthenosphere and the enriched subcontinental lithospheric mantle. (3) The main plateau basaltic volcanism in this region is represented by the 12.4–3.3-Ma-old MLBA basalts and the 8.2–4.4-Ma-old basalts from Meseta Chile Chico (MCC), Chile. Two groups can be distinguished among these main plateau basalts. The first group includes alkali basalts and trachybasalts displaying typical OIB signatures and thought to derive from predominantly asthenospheric mantle sources similar to those of the post-plateau MLBA basalts, but through slightly larger degrees of melting. The second one, although still dominantly alkalic, displays incompatible element signatures intermediate between those of OIB and arc magmas (e.g., La/Nb > 1 and TiO₂ < 2 wt.%). These intermediate basalts differ from their strictly alkalic equivalents by having lower High Field Strength Element (HFSE) and higher εNd (up to + 5.4). These features are consistent with their derivation from an enriched mantle source contaminated by ca. 10% rutile-bearing restite of altered oceanic crust. The petrogenesis of the studied Mio-Pliocene basalts from MLBA and MCC is consistent with contributions of the subslab asthenosphere, the South American subcontinental lithospheric mantle and the subducted Pacific oceanic crust to their sources. However, their chronology of emplacement is not consistent with an ascent through an asthenospheric window opened as a consequence of the subduction of segment SCR-1, which entered the trench at 6 Ma. Indeed, magmatic activity was already important between 12 and 8 Ma in MLBA and MCC as well as in southernmost plateaus, i.e., 6 Ma before the subduction of the SCR-1 segment. We propose a geodynamic model in which OIB and intermediate magmas derived from deep subslab asthenospheric mantle did uprise through a tear-in-the-slab, which formed when the southernmost segments of the SCR collided with the Chile Trench around 15 Ma. During their ascent, they interacted with the Patagonian supraslab mantle and, locally, with slivers of subducted Pacific oceanic crust that contributed to the geochemical signature of the intermediate basalts.     

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.     

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.     

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.     

Transition zone structure under a stationary hot spot: Cape Verde

We report on a two-year seismic deployment in the Cape Verde Islands, one goal of which was to study the upper mantle to determine its structure under a hot spot that is stationary in the hot spot reference frame. We find from analysis of P-to-S receiver functions estimated from broadband seismic recordings that, within uncertainty, the time separation between the 410 and 660 km discontinuities is normal compared to radial earth models. Thus, to exist, even stationary hot spots do not require vertical thermal anomalies from deep melting sources anchored in the lower mantle or at the core–mantle boundary or their anomalies are narrower than ～ 250 km in the upper mantle.     

Crustal intrusion beneath the Louisville hotspot track

We report here the first detailed 2D tomographic image of the crust and upper mantle structure of a Cretaceous seamount that formed during the interaction of the Pacific plate and the Louisville hotspot. Results show that at ～ 1.5 km beneath the seamount summit, the core of the volcanic edifice appears to be dominantly intrusive, with velocities faster than 6.5 km/s. The edifice overlies both high lower crustal (> 7.2–7.6 km/s) and upper mantle (> 8.3 km/s) velocities, suggesting that ultramafic rocks have been intruded as sills rather than underplated beneath the crust. The results suggest that the ratio between the volume of intra-crustal magmatic intrusion and extrusive volcanism is as high as ～ 4.5. In addition, the inversion of Moho reflections shows that the Pacific oceanic crust has been flexed downward by up to ～ 2.5 km beneath the seamount. The flexure can be explained by an elastic plate model in which the seamount emplaced upon oceanic lithosphere that was ～ 10 Myr at the time of loading. Intra-crustal magmatic intrusion may be a feature of hotspot volcanism at young, hot, oceanic lithosphere, whereas, magmatic underplating below a pre-existing Moho may be more likely to occur where a hotspot interacts with oceanic lithosphere that is several tens of millions of years old.     

Mantle transition zone beneath the Caribbean-South American plate boundary and its tectonic implications

We analyzed receiver-function data recorded by a temporary broadband array deployed as part of the BOLIVAR project and the permanent seismic network of Venezuela to study the mantle transition zone structure beneath the Caribbean-South American plate boundary and Venezuela. Significant topography on both the 410-km and the 660-km discontinuities was clearly imaged in the CCP (common-conversion-point) stacked images. Beneath the southeastern Caribbean, the 410-km is featured by a narrow (～ 200 km EW) ～ 25-km uplift extending in the NS direction around 63° west, while the 660-km is depressed by ～ 20 km in a narrow region slightly west to the uplift, a scenario that is more consistent with westward descent of the oceanic South American plate rather than a break-off of NNW dipping proto-Caribbean oceanic lithosphere along the El Pilar Fault. We also found a thick transition zone beneath the Falcon region in northwestern Venezuela, possibly associated with the subducted Nazca plate. A flat 410-km was observed beneath the Guayana shield, suggesting that the shield has a stable and moderately deep keel, which has little effect on the underlying transition zone structure.     

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.     

Geological,tomographic, kinematic and geodynamic constraints on the dynamics of sinking slabs

We use geodynamic models with imposed plate velocities to test the forward-modeled history of subduction based on a particular plate motion model against alternative seismic tomography models. We utilize three alternative published reference frames: a hybrid moving hotspot-palaeomagnetic, a hybrid moving hotspot-true polar wander corrected-palaeomagnetic, and a Subduction Reference Frame, a plate model including longitudinal shifts of subduction zones by matching subduction volumes imaged by P-wave tomography, to assess which model best predicts present day mantle structure compared with seismic tomography and volumetrically derived subduction history. Geodynamic modeling suggests paleo-longitudinal corrections applied to the Subduction Reference Frame result in lower mantle slab material beneath North America and East Asia accumulating up to 10–15° westward of that imaged by tomography, whereas the hybrid models develop material offset by 2–9°. However, the Subduction Reference Frame geodynamic model produces slab material beneath the Tethyan Domain coinciding with slab volumes imaged by tomography, whereas the hybrid reference frame models do not, suggesting regional paleo-longitudinal corrections are required to constrain slab locations. We use our models to test inferred slab sinking rates in the mantle focusing on well-constrained regions. We derive a globally averaged slab-sinking rate of 13 ± 3 mm/yr by combining the ages of onset and cessation of subduction from geological data and kinematic reconstructions with images of subducted slabs in the mantle. Our global average slab-sinking rate overlaps with the 15–20 mm/yr rate implied by mantle convection models using a lower mantle viscosity 100 times higher than the upper mantle.     

Performance of bridges with seismic isolation bearings during the Maule earthquake,Chile

On February 27, 2010 an earthquake of magnitude M_w=8.8, with epicenter in Cobquecura, Maule region, hit the central part of Chile. After the earthquake, a tsunami occurred that caused heavy casualties and damage to buildings and infrastructure. In particular, 4.5% of the overpasses located in the affected region suffered some type of damage and 25 bridges and several pedestrian bridges collapsed. At that time, there were about a dozen bridges with seismic isolation bearings in Chile, two of which were instrumented with accelerometer networks: the Marga Marga Bridge, located in Viña del Mar, and an elevated section of the Metro Line 5 in Santiago, at approximately 300 km and 400 km from the epicenter, respectively. This paper analyzes the acceleration records obtained at these instrumented structures and studies the effect of the seismic isolation on their dynamic response. The beneficial effect of the isolation system, especially in the longitudinal direction, is apparent. In addition, some flaws in the collapsed bridges are described.     

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.     

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.     

Estimates of the transition zone temperature in a mechanically mixed upper mantle

Partial melting of mantle peridotite generates a physically and chemically layered oceanic lithosphere that is cycled back into the mantle in subduction zones. Stirring times of the mantle are too long to allow for complete re-homogenization of subducted basalt and harzburgite, given the low chemical diffusivity of the solid mantle. This suggests that the Earth's mantle is a mechanical mixture of basaltic and harzburgitic components. Using a recently developed thermodynamic formulism we determine the phase equilibria and the seismic properties of a mantle comprised of a mechanical mixture of basalt and harzburgite (MM) and a homogeneous mantle (EA) with identical pyrolitic bulk chemistry. We use the theoretical shear velocity profiles as a new thermometer of the mantle below the magma-genetic zone by modeling the difference ΔT_410-660 between traveltimes of shear wave reflections off the 410-km and 660-km with the potential temperature T_P. ΔT_410-660 are measured from waveform stacks. They indicate that, over 1000+ km wave lengths, the temperature varies by about 200 K. Lowest and highest temperatures are resolved for the western Pacific subduction zones and the central Pacific, respectively. This variation is similar for the EA and MM and is in excellent agreement with estimates of transition zone thickness and shear velocity variations. The median value of T_P for the EA is 1720 K. It is about 1625 K for the MM, a value that is in better agreement with the Normal-MORB values of 1610 ± 40 K inferred from olivine-liquid equilibria given that our sampling region encompasses the Western Pacific subduction zones and the oldest parts of the Pacific Plate. We argue therefore that a mechanical mixed mantle, with generally higher velocities and steeper velocities gradients, represents a better physical reference model than a model based on a fully equilibrated assemblage.     

Crustal structure of the Reykjanes Ridge near 62°N,on the basis of seismic refraction and gravity data

Explosion deep seismic sounding data sections of high quality had been obtained with RV Meteor in the Reykjanes Iceland Seismic Project (RRISP77 [Angenheister, G., Gebrande, H., Miller, H., Goldflam, P., Weigel, W., Jacoby, W.R., Pálmason, G., Björnsson, S., Einarsson, P., Pavlenkova, N.I., Zverev, S., Litvinenko, I.V., Loncarecic, B., Solomon, S., 1980. Reykjanes Ridge Iceland Seismic Experiment (RRISP 77). J. Geophys. 47, 228–238]) which close an information gap near 62°N. Preliminary results were presented by Weigel [Weigel, W., 1980. Aufbau des Reykjanes Rückens nach refraktionsseismischen Messungen. In: Weigel, W. (Ed.), Reykjanes Rücken, Island, Norwegischer Kontinentalrand. Abschlusskolloquium, Hamburg zur Meteor-Expedition, vol. 45. DFG, Bonn, pp. 53–61], and here we report on the data and results of interpretation. Clear refracted phases to 90 km distance permit crustal and uppermost mantle structure to be modelled by ray tracing. The apparent P-wave velocities are around 4.5, 6–6.5, 7–7.6 and 8.2–8.7 km/s, but no wide-angle reflections have been clearly seen. Accompanying sparker reflection data reveal thin sediment ponds in the axial zone and up to 400 m thick sediments at 10 Ma crustal age. Ray tracing reveals the following model below the sediments: (1) a distinct, 1–2 km thick upper crust (layer 2A) with Vp increasing with age (to 10 Ma) from <3.4 to 4.9 km/s and with a vertical gradient of 0.1–0.2 km/s/km, (2) a lower crust or layer 3 beginning at depths of 2 (axis) to 4 km (10 Ma age) below sea level with 6.1–6.8 km/s and similar vertical gradients as above, (3) the lower crust bottoms at 5.2–9.5 km depth below sea level (0–10 Ma) with a marked discontinuity, underneath which (4) Vp rises from about 7.5–7.8 km/s (0–10 Ma) with a positive vertical gradient of, again, 0.1–0.2 km/s/km such that 8 km/s would be reached at 12 km and deeper near the axis. Our preferred interpretation is that the mantle begins at the distinct discontinuity (“Moho”), but a deeper “Moho” of Vp ≈ 8 km/s cannot be excluded. From Iceland southward to 60°N several experiments show a decrease of crustal thickness from 14 to 8 km. Velocity trends with age across the ridge reflect cooling and filling of cracks, and thickness trends probably suggest volcanic productivity variations as previously suggested.Gravity inversion concentrates on a profile across the ridge with the above seismic a priori information; with 0.2–0.5 km depth uncertainty it leads to a good fit (±2.5 mGal where seismic data exist). Best fitting densities are (in kg/m³) for sediments, 2180; upper crust, 2450–2570; lower crust, 2850–2940; mantle lithosphere, 3215–3240 with a deficit for an asthenospheric wedge of no more than −100 kg/m³. The morphological ridges and troughs superimposed on the SE ridge flank are partly correlated, partly anti-correlated with the Bouguer anomaly and suggest that variable crustal density variations accompany the morphology variations.     

New evidence of delamination in the Western Alboran Sea: Geodynamic evolution of the Alboran domain and its margins

The presence of continuous upper crustal blocks between the Iberian Betics and Moroccan Rif in the western and middle Alboran Sea, detected with tomography, can add new information about the lithosphere structure and geodynamic evolution in this region. A large volume of seismic data (P and S wave arrival times) has been collected for the period between 1 December 1988 and 31 December 2008 by 57 stations located in northern Morocco (National Institute of Geophysics, CNRST, Rabat), southern Portugal (Instituto de Meteorologia, Lisbon) and Spain (Instituto Geografico National, Madrid) and used to investigate the lithosphere in the western Alboran Sea region. We use a linearized inversion procedure comprising two steps: (1) finding the minimal 1-D model and simultaneous relocation of hypocenters and (2) determination of local velocity structure using linearized inversion. The model parameterization in this method assumes a continuous velocity field. The resolution tests indicate that the calculated images give near true structure imaged at 5 km depth for the Tanger peninsula, the Alhoceima region and southern Spain. At 15, 30 and 45 km depth we observe a near true structure imaged in northern Morocco, and southern Spain. At 60 and 100 km, southern Spain and the SW region of the Alboran Sea give a near true structure. The resulting tomographic image shows the presence of two upper crustal bodies (velocity 6.5 km/s) at 5–10 km depth between the Betics, Rif, western and central Alboran Sea. Low velocities at the base of these two bodies favor the presence of melt. This new evidence proves that the Tethysian ocean upper crust was not totally collapsed or broken down during the late Oligocene–early Miocene. These two blocks of upper crust were initially one block. The geodynamic process in the eastern of the Mediterranean is driven by slab rollback. The delamination process of the lithospheric mantle terminates with the proposed slab rollback in the western part of the Mediterranean. This can be explained by the removal of the major part of the lithosphere beneath the area, except in the SW part of the Alboran Sea where a small part of the lithospheric mantle is still attached and is extends and dips to SE beneath the Rif, slowly peeled back to the west. A second detached lithospheric mantle is located and extends to eastern part of the Rif and dips to the SE. The removal of lithosphere mantle from the base of the crust was replaced and heated by extrusion of asthenospheric material coming from depth to replace the part of crust detached. A combination of isostatic surface/topographic uplift and erosion induced a rapid exhumation and cooling of deep crustal rocks.