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.     

Enriched,HIMU-type peridotite and depleted recycled pyroxenite in the Canary plume: A mixed-up mantle

The Earth's mantle is chemically and isotopically heterogeneous, and a component of recycled oceanic crust is generally suspected in the convecting mantle [Hofmann and White, 1982. Mantle plumes from ancient oceanic crust. Earth Planet. Sci. Lett. 57, 421–436]. Indeed, the HIMU component (high µ = ²³⁸U/²⁰⁴Pb), one of four isotopically distinct end-members in the Earth's mantle, is generally attributed to relatively old (≥ 1–2 Ga) recycled oceanic crust in the form of eclogite/pyroxenite, e.g. [Zindler and Hart, 1986. Chemical geodynamics. Ann. Rev. Earth Planet. Sci. 14, 493–571]. Although the presence of the recycled component is generally supported by element and isotopic data, little is known about its physical state at mantle depths. Here we show that the concentrations of Ni, Mn and Ca in olivine from the Canarian shield stage lavas, which can be used to assess the physical nature of the source material (peridotite versus olivine-free pyroxenite) [Sobolev et al., 2007. The amount of recycled crust in sources of mantle-derived melts. Science 316, 412–417], correlate strongly with bulk rock Sr, Nd and Pb isotopic ratios. The most important result following from our data is that the enriched, HIMU-type (having higher ²⁰⁶Pb/²⁰⁴Pb than generally found in the other mantle end-members) signature of the Canarian hotspot magmas was not caused by a pyroxenite/eclogite constituent of the plume but appears to have been primarily hosted by peridotite. This implies that the old (older than ~ 1 Ga) ocean crust, which has more evolved radiogenic isotope compositions, was stirred into/reacted with the mantle so that there is not significant eclogite left, whereas younger recycled oceanic crust with depleted MORB isotopic signature (< 1 Ga) can be preserved as eclogite, which when melted can generate reaction pyroxenite.     

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.     

Critical phenomena in thermal conductivity: Implications for lower mantle dynamics

Microscopic mechanisms for heat transport in dense minerals (phonon scattering and photon attenuation) exhibit aspects of threshold behavior, discussed qualitatively here. For all minerals examined so-far using laser-flash analysis, the lattice component of the thermal conductivity of the mantle asymptotes to a constant above a critical temperature of ∼1500 K. Radiative transfer calculated from absorption spectra has thresholds in both grain-size and Fe content, and a rather complex dependence on temperature. These critical phenomena impact convection of the lower mantle, because the lattice contribution tends to destabilize the cold boundary layers, whereas radiative transfer mostly promotes stability in the lower mantle, unless the grains are large and Fe-rich, which makes convection chaotic and time-dependent. The specific behavior suggests that flow in the lower mantle is sluggish, whereas flow in the upper mantle-transition zone is time-dependent. The decrease in k_rad as Fe/(Fe + Mg) increases beyond ∼0.1 may be connected with formation of lower mantle, thermo-chemical plumes through positive feedback.     

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.     

Plume Generation Zones at the margins of Large Low Shear Velocity Provinces on the core–mantle boundary

Large Igneous Province (LIP) eruption sites of the past 300 My lie vertically above 1% slow shear wave velocity (V_s) contours bounding the African and Pacific Large Low Shear Velocity Provinces (LLSVPs) at the core–mantle boundary (CMB), or in the cases of the Siberian and Columbia River LIPs, bounding one or other of two smaller, Low Shear Velocity Provinces (LSVPs). Steep gradients in V_s at the CMB coincide with those 1% slow contours. The sites of 24 active hotspot volcanoes project down to the same narrowly defined borders of the LLSVPs at the CMB. Plumes that have generated LIPs and major hotspot volcanoes have risen only from the immediate neighbourhoods of the 1% slow V_s contours at the CMB which thus define Plume Generation Zones (PGZs). PGZs projected vertically upward approximately match the + 10 m elevation contour of the geoid showing that the LLSVPs are a dominant control on the positively elevated geoid. Minima in the frequency distribution of shear wave velocities in the lowermost mantle near V_s = ? 1% indicate that regions with more negative velocities, forming ～ 2% of total mantle mass, are likely to be of material compositionally different from the rest of the mantle. Because all LIP eruption sites with ages younger than 300 Ma lie above the borders of LLSVPs or LSVPs at the CMB, PGZ footprints are inferred to have remained in the same places for the past 300 My. Because no plumes have risen from the interior of the LLSVPs and because no lithospheric slabs have penetrated those bodies the volumes of the LLSVPs are inferred to have also remained unchanged for the past 300 My. Because the LLSVPs are the dominant control on the positively elevated areas of the geoid those too must have remained as they now are since 300 Ma. The LLSVPs are not rising buoyant objects but stable features of the deep mantle. LIPs have been erupted throughout the past 2.5 Gy indicating that PGZs comparable to those of the past 0.3 Gy and LLSVPs (of which PGZs mark the margins at the CMB) have also existed for at least that long. LLSVPs could thus form the isolated reservoir invoked by some to explain the distinctive isotopic compositions of terrestrial rocks. PGZs lie at places where the boundaries of: (i) The outer core, (ii) one of the LLSVPs or LSVPs, and (iii) the seismically faster part of the deep mantle meet. Horizontal temperature gradients across the steeply inclined margins to the LLSVPs, the interiors of which are hotter than the surrounding mantle, at the CMB are key controls for the generation of plumes. Near the CMB the association of the high temperature of the outer core with an inclined thermal boundary layer at the margins of LLSVPs facilitates the generation of mantle plumes in the PGZs.     

Primitive neon and helium isotopic compositions of high-MgO basalts from the Kerguelen Archipelago,Indian Ocean

The geochemical characteristics of mildly alkalic basalts (24–25 Ma) erupted in the southeastern Kerguelen Archipelago are considered to represent the best estimate for the composition of the enriched Kerguelen plume end-member. A recent study of picrites and high-MgO basalts from this part of the archipelago highlighted the Pb and Hf isotopic variations and suggested the presence of mantle heterogeneities within the Kerguelen plume itself. We present new helium and neon isotopic compositions for olivines from these picrites and high-MgO basalts (6–17 wt.% MgO) both to constrain the enriched composition of the Kerguelen plume and to determine the origin of isotopic heterogeneities involved in the genesis of Kerguelen plume-related basalts. The olivine phenocrysts have extremely variable ⁴He / ³He compositions between MORB and primitive values observed in OIB (∼90,000 to 40,000; i.e., R / R_a ∼8 to 18) and they show primitive neon isotopic ratios (average ²¹Ne / ²¹Ne_ext ∼0.044). The neon isotopic systematics and the ⁴He / ³He ratios that are lower than MORB values for the Kerguelen basalts clearly suggest that the Kerguelen hotspot belongs to the family of primitive hotspots, such as Iceland and Hawaii. The rare gas signature for the Kerguelen samples, intermediate between MORB and solar, is apparently inconsistent with mixing of a primitive component with a MORB-like source, but may result from sampling a heterogeneous part of the mantle with solar ³He / ²²Ne and with a higher (U, Th) / ³He ratio compared to typically high R / R_a hotspot basalts such as those from Iceland and Hawaii.     

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.     

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.     

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.     

On the effects of preferential or barrier flow features on solute plumes in permeable porous media

Despite that discrete flow features (DFFs, e.g. fractures and faults) are common features in the subsurface, few studies have explored the influence of DFFs on solute plumes in otherwise permeable rocks (e.g. sandstone, limestone), compared to low-permeability rock settings (e.g. granite and basalt). DFFs can provide preferential flow pathways (i.e. ‘preferential flow features’; PFFs), or can act to impede flow (i.e. ‘barrier flow features’; BFFs). This research uses a simple analytical expression and numerical modelling to explore how a single DFF influences the steady-state distributions of solute plumes in permeable aquifers. The analysis quantifies the displacement and widening (or narrowing) of a steady-state solute plume as it crosses a DFF in idealised, 1 × 1 m moderately permeable rock aquifers. Previous research is extended by accounting for DFFs as 2D flow features, and including BFF situations. A range of matrix-DFF permeability ratios (0.01 to 100) and DFF apertures (0.25 mm to 2 cm), typical of sedimentary aquifers containing medium-to-large fractures, are considered. The results indicate that for the conceptual models considered here, PFFs typically have a more significant influence on plume distributions than BFFs, and the impact of DFFs on solute plumes generally increases with increasing aperture. For example, displacement of peak solute concentration caused by DFFs exceeds 20 cm in some PFF cases, compared to a maximum of 0.64 cm in BFF cases. PFFs widen plumes up to 9.7 times, compared to a maximum plume widening of 2.0 times in BFF cases. Plumes crossing a PFF are less symmetrical, and peak solute concentrations beneath PFFs are up to two orders of magnitude lower than plumes in BFF cases. This study extends current knowledge of the attenuating influence of DFFs in otherwise permeable rocks on solute plume characteristics, through evaluation of 2D flow effects in DFFs for a variety of DFF apertures, and by considering BFF situations.     

Geochemical component relationships in MORB from the Mid-Atlantic Ridge, 22–35°N

New trace element and Hf, Nd, and Pb isotope data are reported for 22 basalts collected between 22°N and 35°N on the Mid-Atlantic Ridge. (La / Sm)_N ratios identify the presence of enriched (E)-MORB in the northernmost part of this area and normal (N)-MORB elsewhere. A negative correlation is observed when ¹⁴³Nd / ¹⁴⁴Nd is plotted against ²⁰⁶Pb / ²⁰⁴Pb, ²⁰⁷Pb / ²⁰⁴Pb, and ²⁰⁸Pb / ²⁰⁴Pb, whereas ¹⁷⁶Hf / ¹⁷⁷Hf appears not to correlate with any of the other isotopic ratios. The E-MORB samples are characterized by high ²⁰⁶Pb / ²⁰⁴Pb, ²⁰⁷Pb / ²⁰⁴Pb, ²⁰⁸Pb / ²⁰⁴Pb, and low ¹⁴³Nd / ¹⁴⁴Nd. Principal Component Analysis (PCA) of Pb isotopes alone identifies three, and only three, significant geochemical end-members (‘components’). Including Nd and Hf isotopic data in the PCA produces spurious components, partly because of curved mixing relationships, and partly because of fractionation during melting. Our preferred interpretation of why ¹⁷⁶Hf / ¹⁷⁷Hf is decoupled from the other isotopic ratios is, as inferred from recent experimental data, that the Hf isotopic compositions of the melt and the residue fail to equilibrate during melting. A strong correlation between (Sr / Nd)_N and (Eu / Eu*)_N indicates that plagioclase is a residual phase of N-MORB, but not of E-MORB melting. The three end-members identified in this study are the depleted mantle, a common-type component, and an enriched plume-type end-member. The common, or ‘C’-type, end-member is characteristic of E-MORB and may itself be a mixture containing recycled oceanic crust (the MORB suite, terrigenous sediments, and/or oceanic plateaus). The plume-type end-member is likely to represent the lower mantle and may involve some primordial material. It is shown that mantle isochrons in general and the Pb–Pb isochron in particular do not characterize a specific geodynamic process acting to create mantle heterogeneities.     

Elastic thickness structure of South America estimated using wavelets and satellite-derived gravity data

We used a wavelet formulation of the classical spectral isostatic analysis to invert satellite-derived gravity and topography/bathymetry for elastic thickness (T_e) over South America and its surrounding plates. To provide a homogeneous representation of the gravity field for this vast region, we corrected free-air anomalies derived from a combination of terrestrial/marine gravity data with data from the GRACE and CHAMP satellite missions (model EIGEN-CG03C) by a simple Bouguer slab using a smoothed representation of surface relief (wavelengths > 125 km). The resulting Bouguer anomaly compares well with terrestrial data acquired in the Central Andes and allows T_e to be confidently estimated for values greater than 10 km. The T_e map resolves regional-scale features that are well-correlated with known surface structures and shows maximum values of 100 ± 15 km over the Archean–Neoproterozoic core of the continent, decreasing to less than 30 km around continental margins. Several regions of the oceanic plates and continental margins have an elastic thickness less than 10 km. We performed a quantitative analysis by comparing the elastic thickness with the thermal structure predicted from the age of oceanic crust and igneous–metamorphic rocks. This demonstrates that oceanic plates have been weakened by thermal interaction with hotspots and locally by fracturing and hydration near the trench. We observe that only the nucleus of the continent has resisted the thermomechanical weakening induced by the rifting of Africa and South America along the passive margin and the Andean orogeny along the active margin. This latter region shows along-strike variations in T_e that correlate with the geotectonic segmentation of the margin and with the pattern of crustal seismicity. Our results reveal that the rigidity structure follows the segmentation of the seismogenic zone along the subduction fault, suggesting a causal relationship that should be investigated in order to improve the understanding and predictability of great earthquakes and tsunamis.     

Hotspot distribution, gravity, mantle tomography: evidence for plumes

Thermal convection is the motor of Earth dynamics and therefore is the link between plate motions, hotspots, seismic velocity variations in the mantle, and anomalies of the gravity field. Small scale mantle anomalies, such as plumes, do, however, generally escape detection by tomographic methods. It is attempted to approach the problem of detection in a somewhat statistical manner. Correlations are sought between spherical harmonic expansions of the fields under study: the hotspot distribution, mantle velocity variations, gravity, heat flow. Using spherical harmonic representations of global fields implies integration and averaging over the whole globe. Thus, although relationships may remain masked in the space domain by a multitude of effects, tendencies may become visible in the spectra or in appropriate averages.The main results are the following: There is a significant long wavelength (n=2,3) negative correlation between the hotspot density and the P-wave velocity variation in the lower mantle. Positive hotspot density of degree 2 to 9 generally correlates with low seismic velocity in all depths of the upper mantle and with positive gravity. This fits well with plume-type convection. These results are also confirmed regionally for a number of individual mid-ocean ridges and hotspots. The hotspot density and the free air anomalies are distinctly positive above regions of low velocity extending to great depth. The effect is not distinct at ridges with shallow velocity anomalies. In a general way, we suggest that the antipodal upwellings (Pacific, Africa) are divided by downwelling currents around the shrinking Pacific. Plate boundaries can easily move away from their past connections with the deeper mantle. Small scale plume currents seem to be depicted in the hotspot expansion. © 1999 Elsevier Science Ltd. All rights reserved.     

Early crust on top of the Earth's core

In an effort to resolve the current conflict between geochemical requirements for an apparently isolated mantle reservoir and recent geophysical evidence for whole-mantle convection, we investigate the possibility that the region above the core-mantle boundary, termed D″, serves as an early-isolated, rare-gas- and incompatible-element-bearing reservoir, and we propose a mechanism for its formation that is a likely outcome of Earth accretion models. In these models, the most cataclysmic event in Earth history, the moon-forming giant impact on the proto-Earth of a Mars-size planet (perhaps as early as 4540 Ma ago) was followed by accretion of smaller bodies long afterwards (until ∼3900 Ma). Some collisions probably triggered melting, metal segregation and degassing. However, the small bodies, fragments, particles, dust, including those of chondrite-like composition, existed on near-earth orbits, were irradiated by intense solar wind, and finally fell on an early-formed, incompatible-element-enriched basaltic crust without causing extensive melting. The respective regions of the crust, loaded with chondrite-like debris, were therefore significantly enriched in iron. When this mixed material was subducted, the bulk density of the subducted lithosphere exceeded that of the bulk silicate mantle, which had already lost its metallic iron to the core. Segregation of this denser material at the base of the mantle was facilitated by the high temperatures at the core-mantle boundary, which greatly reduce the viscosity, as was quantitatively modelled by Christensen and Hofmann (Christensen, U.R., Hofmann, A.W., 1994. Segregation of subducted oceanic-crust in the convecting mantle. J. Geophys. Res.-Solid Earth 99 (B10), 19867–19884). Assuming a basalt/chondrite mass ratio of about 4/1, we obtain a density contrast of ∼7%, which would stabilize the subducted material between the metal core and silicate mantle.Mass balance considerations and preliminary results of geochemical modelling of the above scenario (similar to that performed by Tolstikhin and Marty [Tolstikhin, I.N., Marty, B., 1998. The evolution of terrestrial volatiles, A view from helium, neon, argon and nitrogen isotope modeling. Chem. Geol. 147, 27–52]) show the potential geochemical importance of D″. (1) Modelling of Pu–U–I–Xe isotope systematics predicts formation of this reservoir early in Earth history, ∼100 Ma after formation of the Solar system. (2) The total amount of heat-generating U, Th, K (and other highly incompatible elements) in D″ exceeds 20% of the Earth inventory, and a similar portion of terrestrial heat is being transferred from the core + D″ into the base of the overlying convecting mantle. (3) D″ is enriched in solar implanted rare gases because the small (re)-accreting debris with high surface/mass ratios will have been subjected to intense radiation by the early sun. (4) Rare gases diffuse from D″ into the overlying mantle and are then transferred into upwelling plumes, which originate above D″. In addition, small amounts of D″ material may be entrained by the mantle convective flow as was recently discussed by Schott et al. [Schott, B., Yuen, D.A., Braun, A., 2002. The influences of composition and temperature-dependent rheology in thermal-chemical convection on entrainment of the D″ layer. Physics Earth Planet. Inter. 129, 43–65]. From the rare-gas modelling it follows that initially (∼4500 Ma ago) D″ could have been more massive by a factor of ∼1.2 than at present (about 2 × 10²⁶ g). The present-day mass flux from D″ into the convecting mantle is estimated to be ≤0.05 × 10¹⁶ g year⁻¹, a factor of ∼100 less than the rate of ridge magmatism. This small contribution of D″ material makes it difficult to trace fingerprints of D″ even using such sensitive tracers as Pb isotope ratios. (5) The density contrast that stabilizes D″ is maintained by its higher intrinsic density due to the iron-rich chondrite-like component. Subduction of this material, its entrainment by convective mantle flow and mixing could also account for the preservation of the chondritic relative abundances of siderophile elements in the mantle. If D″ is partially molten, the density contrast may be caused by a high-density melt fraction.     

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.     

New constraints on the HIMU mantle from neon and helium isotopic compositions of basalts from the Cook–Austral Islands

High ⁴He/³He ratios of 100 000 to 160 000 found at HIMU ocean islands (“high μ,” where μ is the U/Pb ratio) are usually attributed to the presence of recycled oceanic crust in the HIMU mantle source. However, significantly higher ⁴He/³He ratios are expected in recycled crust after residence in the mantle for periods greater than 1 Ga. In order to better understand the helium isotopic signatures in HIMU basalts, we have measured helium and neon isotopic compositions in a suite of geochemically well-characterized basalts from the Cook–Austral Islands. We observe ⁴He/³He ratios ranging from 56 000 to 141 000, suggesting the involvement of mantle reservoirs both more and less radiogenic than the mantle source for mid-ocean ridge basalts (MORBs). In addition, we find that the neon isotopic compositions of HIMU lavas extend from the MORB range to compositions less nucleogenic than MORBs. The Cook-Austral HIMU He–Ne isotopic compositions, along with Sr, Nd, Pb, and Os isotopic compositions, indicate that in addition to recycled crust, a relatively undegassed mantle end-member (e.g., FOZO) is involved in the genesis of these basalts. The association of relatively undegassed mantle material with recycled crust provides an explanation for the close geographical association between HIMU lavas and EM (enriched mantle)-type lavas from this island chain: EM-type signatures represent a higher mixing proportion of relatively undegassed mantle material. Mixing between recycled material and relatively undegassed mantle material may be a natural result of entrainment processes and convective stirring in deep mantle.     

Subduction erosion of forearc mantle wedge implicated in the genesis of the South Sandwich Island (SSI) arc: Evidence from boron isotope systematics

The South Sandwich volcanic arc is sited on a young oceanic crust, erupts low-K tholeiitic rocks, is characterized by unexotic pelagic and volcanogenic sediments on the down-going slab, and simple tectonic setting, and is ideal for assessing element transport through subduction zones. As a means of quantifying processes attending transfer of subduction-related fluids from the slab to the mantle wedge, boron concentrations and isotopic compositions were determined for representative lavas from along the arc. The samples show variable fluid-mobile/fluid-immobile element ratios and high enrichments of B/Nb (2.7 to 55) and B/Zr (0.12 to 0.57), similar to those observed in western Pacific arcs. δ¹¹B values are among the highest so far reported for mantle-derived lavas; these are highest in the central part of the arc (+ 15 to + 18‰) and decrease toward the southern and northern ends (+ 12 to + 14‰). δ¹¹B is roughly positively correlated with B concentrations and with ⁸⁷Sr/⁸⁶Sr ratios, but poorly coupled with other fluid-mobile elements such as Rb, Ba, Sr and U. Peridotites dredged from the forearc trench also have high δ¹¹B (ca. + 10‰) and elevated B contents (38–140 ppm). Incoming pelagic sediments sampled at ODP Site 701 display a wide range in δ¹¹B (+ 5 to ? 13‰; average = ? 4.1‰), with negative values most common. The unusually high δ¹¹B values inferred for the South Sandwich mantle wedge cannot easily be attributed to direct incorporation of subducting slab materials or fluids derived directly therefrom. Rather, the heavy B isotopic signature of the magma sources is more plausibly explained by ingress of fluids derived from subduction erosion of altered frontal arc mantle wedge materials similar to those in the Marianas forearc. We propose that multi-stage recycling of high-δ¹¹B and high-B serpentinite (possibly embellished by arc crust and volcaniclastic sediments) can produce extremely ¹¹B-rich fluids at slab depths beneath the volcanic arc. Infiltration of such fluids into the mantle wedge likely accounts for the unusual magma sources inferred for this arc.     

Geoneutrinos and the energy budget of the Earth

The total energy loss of the Earth is well constrained by heat flux measurements on land, the plate cooling model for the oceans, and the buoyancy flux of hotspots. It amounts to 46 ± 2 TW. The main sources that balance the total energy loss are the radioactivity of the Earth's crust and mantle, the secular cooling of the Earth's mantle, and the energy loss from the core. Only the crustal radioactivity is well constrained. The uncertainty on each of the other components is larger than the uncertainty of the total heat loss. The mantle energy budget cannot be balanced by adding the best estimates of mantle radioactivity, secular cooling of the mantle, and heat flux from the core. Neutrino observatories in deep underground mines can detect antineutrinos emitted by the radioactivity of U and Th. Provided that the crustal contribution to the geoneutrino flux can be very precisely calculated, it will be possible to put robust constraints on mantle radioactivity and its contribution to the Earth's energy budget. Equally strong constraints could be obtained from a deep ocean observatory without the need of crustal correction. In the future, it may become possible to obtain directional information on the geoneutrino flux and to resolve radial variations in concentration of heat producing elements in the mantle.