An anisotropic viscous representation of Mohr-Coulomb failure for use is modelling coupled mantle-continent dynamics

In mantle convection models it has become common to make use of a modified (pressure sensitive, Boussinesq) von Mises yield criterion to limit the maximum stress the lithosphere can support. This approach allows the viscous, cool thermal boundary layer to deform in a relatively plate-like mode even in a fully Eulerian representation. In large-scale models with embedded continental crust where the mobile boundary layer represents the oceanic lithosphere, the von Mises yield criterion for the oceans ensures that the continents experience a realistic broad-scale stress regime. In detailed models of crustal deformation it is, however, more appropriate to choose a Mohr-Coulomb yield criterion based upon the idea that frictional slip occurs on whichever one of many randomly oriented planes happens to be favorably oriented with respect to the stress field. As coupled crust/mantle models become more sophisticated it is important to be able to use whichever failure model is appropriate to a given part of the system. We have therefore developed a way to represent Mohr-Coulomb failure within a code which is suited to mantle convection problems coupled to large-scale crustal deformation. Our approach uses an orthotropic viscous rheology (a different viscosity for pure shear to that for simple shear) to define a prefered plane for slip to occur given the local stress field. The simple-shear viscosity and the deformation can then be iterated to ensure that the yield criterion is always satisfied. We again assume the Boussinesq approximation - neglecting any effect of dilatancy on the stress field. An additional criterion is required to ensure that deformation occurs along the plane aligned with maximum shear strain-rate rather than the perpendicular plane which is formally equivalent in any symmetric formulation. It is also important to allow strain-weakening of the material. The material should remember both the accumulated failure history and the direction of failure. We have included this capacity in a Lagrangian-Integration-point finite element code and will show a number of examples of extension and compression of a crustal block with a Mohr-Coulomb failure criterion, and comparisons between mantle convection models using the von Mises versus the Mohr-Coulomb yield criteria. The formulation itself is general and applies to 2D and 3D problems, although it is somewhat more complicated to identify the slip plane in 3D.

Anisotropy model for mantle convection

The paper presents a new theory for modeling flow in anisotropic, viscous rock. This theory has originally been developed for the simulation of large deformation processes including folding and kinking in multi-layered visco-elastic rock. The orientation of slip planes in the context of crystallographic slip is determined by the normal vector, the so-called director of these surfaces. The model is applied to simulate anisotropic natural mantle convection. We compare the evolution of the director and approximately steady states of isotropic and anisotropic convection. The isotropic case has a simple steady state solution, whereas the orthotropic convection model produces a continuously evolving patterning in tile core of the convection cell which makes only a near-steady condition possible, in which the thermal boundary layer appears to be well aligned with the flow and hence as observed in seismic tomomgraphy strong anistropic.

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.

Arc crust-magma interaction in the Andean Southern Volcanic Zone from thermobarometry, mineral composition, radiogenic isotope and rare earth element systematics of the Azufre-Planchon-Peteroa volcanic complex, Chile

The Andean Southern Volcanic Zone (SVZ) is a vast and complex continental arc that has been studied extensively to provide an understanding of arc-magma genesis, the origin and chemical evolution of the continental crust, and geochemical compositions of volcanic products. The present study focuses on distinguishing the magma/sub-arc crustal interaction of eruptive products from the Azufre-Planchon-Peteroa (APP 35°15'S) volcanic center and other major centers in the Central SVZ (CSVZ 37°S–42°S), Transitional SVZ (TSVZ 34.3–37.0°S), and Northern SVZ (NSVZ 33°S–34°30'S). New Hf and Nd isotopic and trace element data for SVZ centers are consistent with former studies that these magmas experienced variable depths of crystal fractionation, and that crustal assimilation is restricted to the lower crustal depths with an apparent role of garnet. Thermobarometric calculations applied to magma compositions constrain the depth of magma separation from mantle sources in all segments of the SVZ to(70-90 km). Magmatic separation at the APP complex occurs at an average depth of ~50 km which is confined to the mantle lithosphere and the base of the crust suggesting localized thermal abrasion both reservoirs. Thermobarometric calculations indicate that CSVZ primary magmas arise from a similar average depth of (~54 km) which confines magma separation to the asthenospheric mantle. The northwards along-arc Sr-Nd-Hf isotopic data and LREE enrichment accompanied with HREE depletion of SVZ mafic magmas correlates well with northward increasing crustal thickness and decreasing primary melt separation from mantle source regions indicating an increased involvement of lower crustal components in SVZ magma petrogenesis. ^ The study concludes that the development of mature subduction zones over millions of years of continuous magmatism requires that mafic arc derived melts stagnate at lower crustal levels due to density similarities and emplace at lower crustal depths. Basaltic underplating creates localized hot zone environments below major magmatic centers. These regions of high temperature/partial melting, and equilibration with underplated mafic rocks provides the mechanism that controls trace element and isotopic variability of primary magmas of the TSVZ and NSVZ from their baseline CSVZ-like precursors.^

Rare earth element concentrations and geochemistry of carbonate minerals from the Logatchev Hydrothermal Field and the Gakkel Ridge

The ultramafic-hosted Logatchev Hydrothermal Field (LHF) at 15°N on the Mid-Atlantic Ridge and the Arctic Gakkel Ridge (GR) feature carbonate precipitates (aragonite, calcite, and dolomite) in voids and fractures within different types of host rocks. We present chemical and Sr isotopic compositions of these different carbonates to examine the conditions that led to their formation. Our data reveal that different processes have led to the precipitation of carbonates in the various settings. Seawater-like 87Sr/86Sr ratios for aragonite in serpentinites (0.70909 to 0.70917) from the LHF are similar to those of aragonite from the GR (0.70912 to 0.70917) and indicate aragonite precipitation from seawater at ambient conditions at both sites. Aragonite veins in sulfide breccias from LHF also have seawater-like Sr isotope compositions (0.70909 to 0.70915), however, their rare earth element (REE) patterns show a clear positive europium (Eu) anomaly indicative of a small (< 1%) hydrothermal contribution. In contrast to aragonite, dolomite from the LHF has precipitated at much higher temperatures (~100 °C), and yet its 87Sr/86Sr ratios (0.70896 to 0.70907) are only slightly lower than those of aragonite. Even higher temperatures are calculated for the precipitation of deformed calcite veins in serpentine-talc fault schists form north of the LHF. These calcites show unradiogenic 87Sr/86Sr ratios (0.70460 to 0.70499) indicative of precipitation from evolved hydrothermal fluids. A simple mixing model based on Sr mass balance and enthalpy conservation indicates strongly variable conditions of fluid mixing and heat transfers involved in carbonate formation. Dolomite precipitated from a mixture of 97% seawater and 3% hydrothermal fluid that should have had a temperature of approximately 14 °C assuming that no heat was transferred. The much higher apparent precipitation temperatures based on oxygen isotopes (~ 100 °C) may be indicative of conductive heating, probably of seawater prior to mixing. The hydrothermal calcite in the fault schist has precipitated from a mixture of 67% hydrothermal fluid and 33% seawater, which should have had an isenthalpic mixing temperature of ~ 250 °C. The significantly lower temperatures calculated from oxygen isotopes are likely due to conductive cooling of hydrothermal fluid discharging along faults. Rare earth element patterns corroborate the results of the mixing model, since the hydrothermal calcite, which formed from waters with the greatest hydrothermal contribution, has REE patterns that closely resemble those of vent fluids from the LHF. Our results demonstrate, for the first time, that (1) precipitation from pure seawater, (2) conductive heating of seawater, and (3) conductive cooling of hydrothermal fluids in the sub-seafloor all can lead to carbonate precipitation within a single ultramafic-hosted hydrothermal system.

Seismic waveform inversion for core-mantle boundary topography

International audience

Savani: A variable resolution whole-mantle model of anisotropic shear velocity variations based on multiple data sets

International audience

Electrochemistry of Heteroleptic Tris(phthalocyaninato) Rare Earth(III) Complexes

The electrochemical characteristics of a series of heteroleptic tris(phthalocyaninato) complexes with identical rare earths or mixed rare earths (Pc)M(OOPc)M(OOPc) [M = Eu...Lu, Y; H2Pc = unsubstituted phthalocyanine, H2(OOPc) = 3,4,12,13,21,22,30,31-octakis(octyloxy)phthalocyanine] and (Pc)Eu(OOPc)Er(OOPc) have been recorded and studied comparatively by cyclic voltammetry (CV) and differential pulse voltammetry (DPV) in CH2Cl2 containing 0.1 M tetrabutylammonium perchlorate (TBAP). Up to five quasi-reversible one-electron oxidations and four one-electron reductions have been revealed. The half-wave potentials of the first, second and fifth oxidations depend on the size of the metal center, but the fifth changes in the opposite direction to that of the first two. Moreover, the difference in redox potentials of the first oxidation and first reduction for (Pc)M(OOPc)M(OOPc), 0.85−0.98 V, also decreases linearly along with decreasing rare earth ion radius, clearly showing the rare earth ion size effect and indicating enhanced π−π interactions in the triple-deckers connected by smaller lanthanides. This order follows the red-shift seen in the lowest energy band of triple-decker compounds. The electronic differences between the lanthanides and yttrium are more apparent for triple-decker sandwich complexes than for the analogous double-deckers. By comparing triple-decker, double-decker and mononuclear [ZnII] complexes containing the OOPc ligand, the HOMO−LUMO gap has been shown to contract approximately linearly with the number of stacked phthalocyanine ligands.

Sandwich complexes of naphthalocyanine with the rare earth metals

Over the past two decades and in particular the past five years, numerous sandwich-type rare earth complexes containing naphthalocyanine ligands have been synthesized. The more extended delocalized π-electron system of naphthalocyanine in comparison with phthalocyanine generates unique physical, spectroscopic, electrochemical and photoelectrochemical properties which have aroused significant research interest in these compounds. This review summarizes recent progress in research on this important class of molecular materials and overviews the current status of the field.

Comparative Electrochemical Study of Unsubstituted and Substituted Bis(phthalocyaninato) Rare Earth(III) Complexes

The electrochemistry of homoleptic substituted phthalocyaninato rare earth double-decker complexes M(TBPc)2 and M(OOPc)2 [M = Y, La...Lu except Pm; H2TBPc = 3(4),12(13),21(22),30(31)-tetra-tert-butylphthalocyanine, H2OOPc = 3,4,12,13,21,22,30,31-octakis(octyloxy)phthalocyanine] has been comparatively studied by cyclic voltammetry (CV) and differential pulse voltammetry (DPV) in CH2Cl2 containing 0.1 M tetra-n-butylammonium perchlorate (TBAP). Two quasi-reversible one-electron oxidations and three or four quasi-reversible one-electron reductions have been revealed for these neutral double-deckers of two series of substituted complexes, respectively. For comparison, unsubstituted bis(phthalocyaninato) rare earth analogues M(Pc)2 (M = Y, La...Lu except Pm; H2Pc = phthalocyanine) have also been electrochemically investigated. Two quasi-reversible one-electron oxidations and up to five quasi-reversible one-electron reductions have been revealed for these neutral double-decker compounds. The three bis(phthalocyaninato)cerium compounds display one cerium-centered redox wave between the first ligand-based oxidation and reduction. The half-wave potentials of the first and second oxidations and first reduction for double-deckers of the tervalent rare earths depend on the size of the metal center. The difference between the redox potentials of the second and third reductions for MIII(Pc)2, which represents the potential difference between the first oxidation and first reduction of [MIII(Pc)2]−, lies in the range 1.08−1.37 V and also gradually diminishes along with the lanthanide contraction, indicating enhanced π−π interactions in the double-deckers connected by the smaller, lanthanides. This corresponds well with the red-shift of the lowest energy band observed in the electronic absorption spectra of reduced double-decker [MIII(Pc′)2]− (Pc′ = Pc, TBPc, OOPc).