Mantle beneath the Gibraltar Arc from receiver functions

P and S receiver functions (PRF and SRF) from 19 seismograph stations in the Gibraltar Arc and the Iberian Massif reveal new details of the regional deep structure. Within the high-velocity mantle body below southern Spain the 660-km discontinuity is depressed by at least 20 km. The Ps phase from the 410-km discontinuity is missing at most stations in the Gibraltar Arc. A thin (similar to 50 km) low-S-velocity layer atop the 410-km discontinuity is found under the Atlantic margin. At most stations the S410p phase in the SRFs arrives 1.0-2.5 s earlier than predicted by IASP91 model, but, for the propagation paths through the upper mantle below southern Spain, the arrivals of S410p are delayed by up to +1.5 s. The early arrivals can be explained by elevated Vp/Vs ratio in the upper mantle or by a depressed 410-km discontinuity. The positive residuals are indicative of a low (similar to 1.7 versus similar to 1.8 in IASP91) Vp/Vs ratio. Previously, the low ratio was found in depleted lithosphere of Precambrian cratons. From simultaneous inversion of the PRFs and SRFs we recognize two types of the mantle: 'continental' and 'oceanic'. In the 'continental' upper mantle the S-wave velocity in the high-velocity lid is 4.4-4.5 km s(-1), the S-velocity contrast between the lid and the underlying mantle is often near the limit of resolution (0.1 km s(-1)), and the bottom of the lid is at a depth reaching 90 100 km. In the 'oceanic' domain, the S-wave velocities in the lid and the underlying mantle are typically 4.2-4.3 and similar to 4.0 km s(-1), respectively. The bottom of the lid is at a shallow depth (around 50 km), and at some locations the lid is replaced by a low S-wave velocity layer. The narrow S-N-oriented band of earthquakes at depths from 70 to 120 km in the Alboran Sea is in the 'continental' domain, near the boundary between the 'continental' and 'oceanic' domains, and the intermediate seismicity may be an effect of ongoing destruction of the continental lithosphere.

Unravelling the interaction between tectonic and sedimentary processes during lithospheric thinning in the Alpine Tethys margins

The discovery of exhumed continental mantle and hyper-extended crust in present-day magma-poor rifted margins is at the origin of a paradigm shift within the research field of deep-water rifted margins. It opened new questions about the strain history of rifted margins and the nature and composition of sedimentary, crustal and mantle rocks in rifted margins. Thanks to the benefit of more than one century of work in the Alps and access to world-class outcrops preserving the primary relationships between sediments and crustal and mantle rocks from the fossil Alpine Tethys margins, it is possible to link the subsidence history and syn-rift sedimentary evolution with the strain distribution observed in the crust and mantle rocks exposed in the distal rifted margins. In this paper, we will focus on the transition from early to late rifting that is associated with considerable crustal thinning and a reorganization of the rift system. Crustal thinning is at the origin of a major change in the style of deformation from high-angle to low-angle normal faulting which controls basin-architecture, sedimentary sources and processes and the nature of basement rocks exhumed along the detachment faults in the distal margin. Stratigraphic and isotopic ages indicate that this major change occurred in late Sinemurian time, involving a shift of the syn-rift sedimentation toward the distal domain associated with a major reorganization of the crustal structure with exhumation of lower and middle crust. These changes may be triggered by mantle processes, as indicated by the infiltration of MOR-type magmas in the lithospheric mantle, and the uplift of the Brianconnais domain. Thinning and exhumation of the crust and lithosphere also resulted in the creation of new paleogeographic domains, the Proto Valais and Liguria-Piemonte domains. These basins show a complex, 3D temporal and spatial evolution that might have evolved, at least in the case of the Liguria-Piemonte basin, in the formation of an embryonic oceanic crust. The re-interpretation of the rift evolution and the architecture of the distal rifted margins in the Alps have important implications for the understanding of rifted margins worldwide, but also for the paleogeographic reconstruction of the Alpine domain and its subsequent Alpine compressional overprint.

Upper-mantle seismic structure beneath SE and Central Brazil from P- and S-wave regional traveltime tomography

We present models for the upper-mantle velocity structure beneath SE and Central Brazil using independent tomographic inversions of P- and S-wave relative arrival-time residuals (including core phases) from teleseismic earthquakes. The events were recorded by a total of 92 stations deployed through different projects, institutions and time periods during the years 1992-2004. Our results show correlations with the main tectonic structures and reveal new anomalies not yet observed in previous works. All interpretations are based on robust anomalies, which appear in the different inversions for P-and S-waves. The resolution is variable through our study volume and has been analyzed through different theoretical test inversions. High-velocity anomalies are observed in the western portion of the Sao Francisco Craton, supporting the hypothesis that this Craton was part of a major Neoproterozoic plate (San Franciscan Plate). Low-velocity anomalies beneath the Tocantins Province (mainly fold belts between the Amazon and Sao Francisco Cratons) are interpreted as due to lithospheric thinning, which is consistent with the good correlation between intraplate seismicity and low-velocity anomalies in this region. Our results show that the basement of the Parana Basin is formed by several blocks, separated by suture zones, according to model of Milani & Ramos. The slab of the Nazca Plate can be observed as a high-velocity anomaly beneath the Parana Basin, between the depths of 700 and 1200 km. Further, we confirm the low-velocity anomaly in the NE area of the Parana Basin which has been interpreted by VanDecar et al. as a fossil conduct of the Tristan da Cunha Plume related to the Parana flood basalt eruptions during the opening of the South Atlantic.

Investigation of Mantle Dynamics from Platinum Group Elements and Rhenium-Osmium Isotope Systematics of Mantle Xenoliths from Oahu, Hawaii

Intraplate volcanism that has created the Hawaiian-Emperor seamount chain is generally thought to be formed by a deep-seated mantle plume. While the idea of a Hawaiian plume has not met with substantial opposition, whether or not the Hawaiian plume shows any geochemical signal of receiving materials from the Earth’s Outer Core and how the plume may or may not be reacting with the overriding lithosphere remain debatable issues. In an effort to understand how the Hawaiian plume works I report on the first in-situ sulfides and bulk rock Platinum Group Element (PGE) concentrations, together with Os isotope ratios on well-characterized garnet pyroxenite xenoliths from the island of Oahu in Hawaii. The sulfides are Fe-Ni Monosulfide Solid Solution and show fractionated PGE patterns. Based on the major elements, Platinum Group Elements and experimental data I interpret the Hawaiian sulfides as an immiscible melt that separated from a melt similar to the Honolulu Volcanics (HV) alkali lavas at a pressure-temperature condition of 1530 ± 100OC and 3.1±0.6 GPa., i.e. near the base or slightly below the Pacific lithosphere. The 187Os/188Os ratios of the bulk rock vary from subchondritic to suprachondritic (0.123-0.164); and the 187Os/188Os ratio strongly correlates with major element, High Field Strength Element (HFSE), Rare Earth Element (REE) and PGE abundances. These correlations strongly suggest that PGE concentrations and Os isotope ratios reflect primary mantle processes. I interpret these correlations as the result of melt-mantle reaction at the base of the lithosphere: I suggest that the parental melt that crystallized the pyroxenites selectively picked up radiogenic Os from the grain boundary sulfides, while percolating through the Pacific lithosphere. Thus the sampled pyroxenites essentially represent crystallized melts from different stages of this melt-mantle reaction process at the base of the lithosphere. I further show that the relatively low Pt/Re ratios of the Hawaiian sulfides and the bulk rock pyroxenites suggest that, upon ageing, such pyroxenites plus their sulfides cannot generate the coupled 186Os- 187Os isotope enrichments observed in Hawaiian lavas. Therefore, recycling of mantle sulfides of pyroxenitic parentage is unlikely to explain the enriched Pt-Re-Os isotope systematics of plume-derived lavas.

Melt/peridotite interaction in the Southern Lanzo peridotite: Field, textural and geochemical evidence

This paper presents field, petrographic-structural and geochemical data on spinet and plagioclase peridotites from the southern domain of the Lanzo ophiolitic peridotite massif (Western Alps). Spinet lherzolites, harzburgites and dunites crop out at Mt. Arpone and Mt. Musine. Field evidence indicates that pristine porphyroclastic spinet lherzolites are transformed to coarse granular spinet harzburgites, which are in turn overprinted by plagioclase peridotites, while strongly depleted spinet harzburgite and dunite bands and bodies replace the plagioclase peridotites. On the northern flank of Mt. Arpone, deformed, porphyroclastic (lithospheric) lherzolites, with diffuse pyroxenite banding, represent the oldest spinel-facies rocks. They show microstructures of a composite subsolidus evolution, suggesting provenance from deeper (asthenospheric) mantle levels and accretion to the lithosphere. These protoliths are locally transformed to coarse granular (reactive) spinet harzburgites and dunites, which show textures reminiscent of melt/rock reaction and geochemical characteristics suggesting that they are products of peridotite interaction with reactively percolating melts. Geochemical data and modelling suggest that <1-5% fractional melting of spinel-facies DMM produced the injected melts. Plagioclase peridotites are hybrid rocks resulting from pre-existing spinet peridotites and variable enrichment of plagioclase and micro-gabbroic material by percolating melts. The impregnating melts attained silica-saturation, as testified by widespread orthopyroxene replacement of olivine, during open system migration in the lithosphere. At Mt. Musine, coarse granular spinet harzburgite and dunite bodies replace the plagioclase peridotites. Most of these replacive, refractory peridotites have interstitial magmatic clinopyroxene with trace element compositions in equilibrium with MORB, while some Cpx have REE-depleted patterns suggesting transient geochemical features of the migrating MORB-type melts, acquired by interaction with the ambient plagioclase peridotite. These replacive spinet harzburgite and dunite bodies are interpreted as channels exploited for focused and reactive migration of silica-undersaturated melts with aggregate MORB compositions. Such melts were unrelated to the silica-saturated melts that refertilized the pre-existing plagioclase peridotites. Finally, MORB melt migration occurred along open fractures, now recorded as gabbroic dikes. Our data document the complexity of rock-types and mantle processes in the South Lanzo peridotite massif and describe a composite tectonic and magmatic scenario that is not consistent with the ``asthenospheric scenario'' proposed by previous authors. We envisage a ``transitional scenario'' in which extending subcontinental lithospheric mantle was strongly modified (both depleted and refertilized) by early melts with MORB-affinity formed by decompression partial melting of the upwelling asthenosphere, during pre-oceanic rifting and lithospheric thinning in the Ligurian Tethys realm. (C) 2006 Elsevier B.V. All rights reserved.

Geochemistry of ODP Leg 125 peridotites

Ocean Drilling Program Leg 125 recovered serpentined harzburgites and dunites from a total of jive sites on the crests and flanks of two serpen finite seamounts, Conical Seamount in the Mariana forearc and Torishima Forearc Seamount in the Izu-Bonin forearc. These are some of the first extant forearc peridotites reported in the literature and they provide a window into oceanic, supra-subduction zone (SSZ) mantle processes. Harzbutrgites from both seamounts are very refractory with low modal clinopyroxene (<4%), chrome-rich spinels (cx-number = 0.40-0.80), very low incompatible element contents, and (with the exception of amphibole-bearing samples) U-shaped rare earth element (REE) profiles with positive Eu anomalies. Both sets of peridotites have olivine-spinel equilibration temperatures that are low compared with abyssal peridotites, possibly because of water-assisted diffusional equilibration in the SSZ environment However, other features indicate that the harzburgites from the two seamounts have very different origins. Harzburgites from Conical Seamount are characterized by calculated oxygen fugacities between FMQ (fayalite- magnetite- quartz) - 1.1 (log units) and FMQ + 0.4 which overlap those of mid-ocean ridge basalt (MORB) peridotites. Dunites from Conical Seamotmt contain small amounts of clinopyroxene, orthopyroxene and amphibole and are light REE (LREE) enriched. Moreover; they are considerably more oxidized than the harzburgites to which they are spatially related, with calculated oxygen fugacities of FMQ -0.2 toFMQ + 1.2. Using textural and geochemical evidence, we interpret these harzburgites as residual MORB mantle (from 15 to 20 % fractional melting) which has subsequently been modified by interaction with boninitic melt ivithin the mantle wedge, and these dunites as zones of focusing of this melt in which pyroxene has preferentially been dissolved from the harzbutgite protolith. In contrast, harzburgites from Torishima Forearc Seamount give calculated oxygen fugacities between FMQ + 0.8 and FMQ + l.6, similar to those calculated for other subduction-zone related peridotites and similar to those calculated for the dunites (FMQ + 1.2 to FMQ + 1.8) from the same seamount. In this case, we interpret both the harzburgites and dunites as linked to mantle melting (20-25 % fractional melting) in a supra-subduction zone environment The results thus indicate that the forearc is underlain by at least two types of mantle lithosphere, one being trapped or accreted oceanic lithosphere, the other being lithosphere formed by subduction-related melting. They also demonstrate that both types of mantle lithosphere may have undergone extensive interaction with subduction-derived magmas.

Geochemistry of Miocene Alborán Basin volcanism

We present new major and trace element and O-Sr-Nd-isotope data for igneous rocks from the western Mediterranean Alborán Sea, collected during the METEOR 51/1 cruise, and for high-grade schists and gneisses from the continental Alborán basement, drilled during the Ocean Drilling Programme (ODP Leg 161, Site 976). The geochemical data allow a detailed examination of crustal and mantle processes involved in the petrogenesis of the lavas and for the first time reveal a zonation of the Miocene Alborán Sea volcanism: (1) a keel-shaped area of LREE-depleted (mainly tholeiitic series) lavas in the central Alborán Sea, generated by high degrees of partial melting of a depleted mantle source and involving hydrous fluids from subducted marine sediments, that is surrounded by (2) a horseshoe-shaped zone with LREE-enriched (mainly calc-alkaline series) lavas subparallel to the arcuate Betic-Gibraltar-Rif mountain belt. We propose that the geochemical zonation of the Miocene Alborán Basin volcanism results from eastward subduction of Tethys oceanic lithosphere coupled with increasing lithospheric thickness between the central Alborán Sea and the continental margins of Iberia and Africa.

Geochemical and Sr, Nd, Pb isotope investigation of the New Caledonia ophiolite

The New Caledonia ophiolite hosts one of the largest obducted mantle section in the world, hence providing a unique insight for the study of upper mantle processes. These mantle rocks belong to an “atypical” ophiolitic sequence, which is dominated by refractory harzburgites but it also includes minor spinel and plagioclase lherzolites. Upper crust is notably absent in the ophiolite, with the exception of some mafic-ultramafic cumulates cropping out in the southern part of the island. Although the New Caledonia ophiolite has been under investigation for decades, its ultra-depleted nature has made its characterization an analytical challenge, so that few trace element data are available, while isotopic data are completely missing. In this thesis a comprehensive geochemical study (major, trace element and Sr-Nd-Pb isotopes) of the peridotites and the associated intrusive mafic rocks from the New Caledonia ophiolite has been carried out. The peridotites are low-strain tectonites showing porphyroclastic textures. Spinel lherzolites are undepleted lithotypes, as attested by the presence of 7-8 vol% of Na2O and Al2O3-rich clinopyroxene (up to 0.5 wt% Na2O; 6.5 wt% Al2O3), Fo content of olivine (88.5-90.0 mol%) and low Cr# of spinel (13-17). Conversely, harzburgites display a refractory nature, proven by the remarkable absence of primary clinopyroxene, very high Fo content in olivine (90.9-92.9 mol%), high Mg# in orthopyroxene (89.8-94.2) and Cr# in spinel (39-71). REE contents show abyssal-type patterns for spinel lherzolites, while harzburgites display U-shaped patterns, typical of fore-arc settings. Spinel lherzolites REE compositions are consistent with relatively low degree (8-9%) of fractional melting of a DMM source, starting in the garnet stability field. Conversely, REE models for harzburgites indicate high melting degrees (20-25%) of a DMM mantle source under spinel faies conditions, consistent with hydrous melting in forearc setting. Plagioclase lherzolites exhibit melt impregnation microtextures, Cr- and TiO2-enriched spinels and REE, Ti, Y, Zr progressive increase with respect to spinel lherzolites. Impregnation models indicate that plagioclase lherzolites may derive from spinel lherzolites by entrapment of highly depleted MORB melts in the shallow oceanic lithosphere. Mafic intrusives are olivine gabbronorites with a very refractory composition, as attested by high Fo content of olivine (87.3-88.9 mol.%), very high Mg# of clinopyroxene (87.7-92.2) and extreme anorthitic content of plagioclase (An = 90-96 mol%). The high Mg#, low TiO2 concentrations in pyroxenes and the anorthitic composition of plagioclase point out an origin from ultra-depleted primitive magmas in a convergent setting. Geochemical trace element models show that the parental melts of gabbronorites are primitive magmas with striking depleted compositions, bearing only in part similarities with the primitive boninitic melts of Bonin Islands. The first Sr, Nd and Pb isotope data obtained for the New Caledonia ophiolite highlight the presence of DM mantle source variably modified by different processes. Nd-Sr-Pb isotopic ratios for the lherzolites (+6.98≤epsilon Ndi≤+10.97) indicate a DM source that suffered low-temperature hydrothermal reactions. Harzburgites are characterized by a wide variation of Sr, Nd and Pb isotopic values, extending from DM-type to EM2 compositions (-0.82≤ epsilon Ndi≤+17.55), suggesting that harzburgite source was strongly affected by subduction-related processes. Conversely, combined trace element and Sr-Nd-Pb isotopic data for gabbronorites indicate a derivation from a source with composition similar to Indian-type mantle, but affected by fluid input in subduction environment. These geochemical features point out an evolution in a pre-Eocenic marginal basin setting, possibly in the proximity of a transform fault, for the lherzolites. Conversely, the harzburgites acquired their main geochemical and isotopic fingerprint in subduction zone setting.

Total sediment thickness grid of the Southern Pacific Ocean off West Antarctica, links to NetCDF files

Paleotopographic models of the West Antarctic margin, which are essential for robust simulations of paleoclimate scenarios, lack information on sediment thickness and geodynamic conditions, resulting in large uncertainties. A new total sediment thickness grid spanning the Ross Sea-Amundsen Sea-Bellingshausen Sea basins is presented and is based on all the available seismic reflection, borehole, and gravity modeling data offshore West Antarctica. This grid was combined with NGDC's global 5 arc minute grid of ocean sediment thickness (Whittaker et al., 2013, doi:10.1002/ggge.20181) and extends the NGDC grid further to the south. Sediment thickness along the West Antarctic margin tends to be 3-4 km larger than previously assumed. The sediment volume in the Bellingshausen, Amundsen, and Ross Sea basins amounts to 3.61, 3.58, and 2.78 million km³, respectively. The residual basement topography of the South Pacific has been revised and the new data show an asymmetric trend over the Pacific-Antarctic Ridge. Values are anomalously high south of the spreading ridge and in the Ross Sea area, where the topography seems to be affected by persistent mantle processes. In contrast, the basement topography offshore Marie Byrd Land cannot be attributed to dynamic topography, but rather to crustal thickening due to intraplate volcanism. Present-day dynamic topography models disagree with the presented revised basement topography of the South Pacific, rendering paleotopographic reconstructions with such a limited dataset still fairly uncertain.

Anisotropic convection model for the earth's mantle

The paper presents a theory for modeling flow in anisotropic, viscous rock. This theory has originally been developed for the simulation of large deformation processes including the folding and kinking of multi-layered visco-elastic rock (Muhlhaus et al. [1,2]). The orientation of slip planes in the context of crystallographic slip is determined by the normal vector - the director - of these surfaces. The model is applied to simulate anisotropic mantle convection. We compare the evolution of flow patterns, Nusselt number and director orientations for isotropic and anisotropic rheologies. In the simulations we utilize two different finite element methodologies: The Lagrangian Integration Point Method Moresi et al [8] and an Eulerian formulation, which we implemented into the finite element based pde solver Fastflo (www.cmis.csiro.au/Fastflo/). The reason for utilizing two different finite element codes was firstly to study the influence of an anisotropic power law rheology which currently is not implemented into the Lagrangian Integration point scheme [8] and secondly to study the numerical performance of Eulerian (Fastflo)- and Lagrangian integration schemes [8]. It turned out that whereas in the Lagrangian method the Nusselt number vs time plot reached only a quasi steady state where the Nusselt number oscillates around a steady state value the Eulerian scheme reaches exact steady states and produces a high degree of alignment (director orientation locally orthogonal to velocity vector almost everywhere in the computational domain). In the simulations emergent anisotropy was strongest in terms of modulus contrast in the up and down-welling plumes. Mechanisms for anisotropic material behavior in the mantle dynamics context are discussed by Christensen [3]. The dominant mineral phases in the mantle generally do not exhibit strong elastic anisotropy but they still may be oriented by the convective flow. Thus viscous anisotropy (the main focus of this paper) may or may not correlate with elastic or seismic anisotropy.

Mantle source heterogeneity, magma generation and magmatic evolution at Terceira Island (Azoresarchipelago): Constraints from elemental and isotopic (Sr, Nd, Hf, and Pb) data

This work addresses the present-day (<100 ka) mantle heterogeneity in the Azores region through the study of two active volcanic systems from Terceira Island. Our study shows that mantle heterogeneities are detectable even when "coeval" volcanic systems (Santa Barbara and Fissural) erupted less than 10 km away. These volcanic systems, respectively, reflect the influence of the Terceira and D. Joao de Castro Bank end-members defined by Beier et at (2008) for the Terceira Rift Santa Barbara magmas are interpreted to be the result of mixing between a HIMU-type component, carried to the upper mantle by the Azores plume, and the regional depleted MORB magmas/source. Fissural lavas are characterized by higher Ba/Nb and Nb/U ratios and less radiogenic Pb-206/Pb-204, Nd-143/Nd-144 and Hf-176/Hf-177, requiring the small contribution of delaminated sub-continental lithospheric mantle residing in the upper mantle. Published noble gas data on lavas from both volcanic systems also indicate the presence of a relatively undegassed component, which is interpreted as inherited from a lower mantle reservoir sampled by the ascending Azores plume. As inferred from trace and major elements, melting began in the garnet stability field, while magma extraction occurred within the spinel zone. The intra-volcanic system's chemical heterogeneity is mainly explained by variable proportions of the above-mentioned local end-members and by crystal fractionation processes. (C) 2011 Elsevier By. All rights reserved.