Journal: Geochemistry, Geophysics, Geosystems

Loading...

Abbreviation

Geochem. Geophys. Geosyst.

Publisher

American Geophysical Union

Journal Volumes

ISSN

1525-2027

Description

Search Results

Publications 1 - 10 of 89
  • Gülcher, Anna; Golabek, Gregor J.; Thielmann, Marcel; et al. (2022)
    Geochemistry, Geophysics, Geosystems
    The rheological properties of Earth's lower mantle are key for mantle dynamics and planetary evolution. The main rock-forming minerals in the lower mantle are bridgmanite (Br) and smaller amounts of ferropericlase (Fp). Previous work has suggested that the large differences in viscosity between these minerals greatly affect the bulk rock rheology. The resulting effective rheology becomes highly strain-dependent as weaker Fp minerals become elongated and eventually interconnected. This implies that strain localization may occur in Earth's lower mantle. So far, there have been no studies on global-scale mantle convection in the presence of such strain-weakening (SW) rheology. Here, we present 2D numerical models of thermo-chemical convection in spherical annulus geometry including a new strain-dependent rheology formulation for lower mantle materials, combining rheological weakening and healing terms. We find that SW rheology has several direct and indirect effects on mantle convection. The most notable direct effect is the changing dynamics of weakened plume channels as well as the formation of larger thermochemical piles at the base of the mantle. The weakened plume conduits act as lubrication channels in the mantle and exhibit a lower thermal anomaly. SW rheology also reduces the overall viscosity, notable in terms of increasing convective vigor and core-mantle boundary heat flux. Finally, we put our results into context with existing hypotheses on the style of mantle convection and mixing. Most importantly, we suggest that the new kind of plume dynamics may explain the discrepancy between expected and observed thermal anomalies of deep-seated mantle plumes on Earth.
  • Tackley, Paul (2000)
    Geochemistry, Geophysics, Geosystems
    Presented here are self-consistent, three-dimensional simulations of mantle convection, some of which display an approximation of plate tectonic behavior that is continuous in space and time. Plate behavior arises through a reasonable material description of silicate deformation, with a simple yield stress being sufficient to give first-order plate-like behavior; however, the required yield strength or fault frictional coefficient is much less than experimentally determined values. Toroidal:poloidal ratios are within geologically observed limits. The sensitivity of the system to yield strength and the form of strength envelope is systematically investigated. Optimum plate character is obtained in a narrow range of yield strength, below which diffuse boundaries, and above which episodic behavior, and eventually a rigid lid, are observed. Models with mobile lids develop very long wavelength horizontal structure, the longest wavelength possible in the domain. Two-dimensional models display much greater time dependence than three-dimensional models.
  • Nakagawa, Takashi; Tackley, Paul J.; Deschamps, Frédéric; et al. (2009)
    Geochemistry, Geophysics, Geosystems
    Phase assemblages of mantle rocks calculated from the ratios of five oxides (CaO-FeO-MgO-Al₂O₃-SiO₂) by free energy minimization were used to calculate the material properties density, thermal expansivity, specific heat capacity, and seismic velocity as a function of temperature, pressure, and composition, which were incorporated into a numerical thermochemical mantle convection model in a 3-D spherical shell. The advantage of using such an approach is that thermodynamic parameters are included implicitly and self-consistently, obviating the need for ad hoc parameterizations of phase transitions which can be complex in regions such as the transition zone particularly if compositional variations are taken into account. Convective planforms for isochemical and thermochemical cases are, however, not much different from those computed using our previous, simple parameterized reference state, which means that our previous results are robust in this respect. The spectrum and amplitude of seismic velocity anomalies obtained using the self-consistently calculated material properties are more “realistic” than those obtained when seismic velocity is linearly dependent on temperature and composition because elastic properties are dependent on phase relationship of mantle minerals, in other words, pressure and temperature. In all cases, the spectra are dominated by long wavelengths (spherical harmonic degree 1 to 2), similar or even longer wavelength than seismic tomographic models of Earth, which is probably due to self-consistent plate tectonics and depth-dependent viscosity. In conclusion, this combined approach of mantle convection and self-consistently calculated mineral physics is a powerful and useful technique for predicting thermal-chemical-phase structures in Earth's mantle. However, because of uncertainties in various parameters, there are still some shortcomings in the treatment of the postperovskite phase transition. Additionally, transport properties such as thermal conductivity and viscosity are not calculated by this treatment and are thus subject to the usual uncertainties.
  • Cloetingh, Sierd; Koptev, Alexander; Kovács, István; et al. (2021)
    Geochemistry, Geophysics, Geosystems
    Although many different mechanisms for subduction initiation have been proposed, only few of them are viable in terms of consistency with observations and reproducibility in numerical experiments. In particular, it has recently been demonstrated that intra‐oceanic subduction triggered by an upwelling mantle plume could greatly contribute to the onset and operation of plate tectonics in the early and, to a lesser degree, modern Earth. On the contrary, the initiation of intra‐continental subduction still remains underappreciated. Here we provide an overview of 1) observational evidence for upwelling of hot mantle material flanked by downgoing proto‐slabs of sinking continental mantle lithosphere, and 2) previously published and new numerical models of plume‐induced subduction initiation. Numerical modeling shows that under the condition of a sufficiently thick (>100 km) continental plate, incipient downthrusting at the level of the lowermost lithospheric mantle can be triggered by plume anomalies of moderate temperatures and without significant strain‐ and/or melt‐related weakening of overlying rocks. This finding is in contrast with the requirements for plume‐induced subduction initiation within oceanic or thinner continental lithosphere. As a result, plume‐lithosphere interactions within continental interiors of Paleozoic‐Proterozoic‐(Archean) platforms are the least demanding (and thus potentially very common) mechanism for initiation of subduction‐like foundering in the Phanerozoic Earth. Our findings are supported by a growing body of new geophysical data collected in various intra‐continental areas. A better understanding of the role of intra‐continental mantle downthrusting and foundering in global plate tectonics and, particularly, in the initiation of “classic” ocean‐continent subduction will benefit from more detailed follow‐up investigations.
  • Curzi, Manuel; Caracausi, Antonio; Rossetti, Federico; et al. (2022)
    Geochemistry, Geophysics, Geosystems
    Post-orogenic back-arc magmatism is accompanied by hydrothermal ore deposits and mineralizations derived from mantle and crustal sources. We investigate Zannone Island (ZI), back-arc Tyrrhenian basin, Italy, to define the source(s) of mineralizing hydrothermal fluids and their relationships with the regional petrological-tectonic setting. On ZI, early Miocene thrusting was overprinted by late Miocene post-orogenic extension and related hydrothermal alteration. Since active submarine hydrothermal outflow is reported close to the island, Zannone provides an ideal site to determine the P-T-X evolution of the long-lived hydrothermal system. We combined field work with microstructural analyses on syn-tectonic quartz veins and carbonate mineralizations, X-ray diffraction analysis, microthermometry and element mapping of fluid inclusions (FIs), C, O, and clumped isotopes, and analyses of noble gases (He-Ne-Ar) and CO2 content in FIs. Our results document the evolution of a fluid system of magmatic origin with increasing mixing of meteoric fluids. Magmatic fluids were responsible for quartz veins precipitation at similar to 125 to 150 MPa and similar to 300 degrees C-350 degrees C. With the onset of extensional faulting, magmatic fluids progressively interacted with carbonate rocks and mixed with meteoric fluids, leading to (a) host rock alteration with associated carbonate and minor ore mineral precipitation, (b) progressive fluid neutralization, (c) cooling of the hydrothermal system (from similar to 320 degrees C to similar to 86 degrees C), and (d) embrittlement and fracturing of the host rocks. Both quartz and carbonate mineralizations show noble gases values lower than those from the adjacent active volcanic areas and submarine hydrothermal systems, indicating that the fossil-to-active hydrothermal history is associated with the emplacement of multiple magmatic intrusions.
  • Global variations in H2O/Ce
    Item type: Journal Article
    Cooper, Lauren B.; Ruscitto, Daniel M.; Plank, Terry; et al. (2012)
    Geochemistry, Geophysics, Geosystems
    We have calculated slab fluid temperatures for 51 volcanoes in 10 subduction zones using the newly developed H₂O/Ce thermometer. The slab fluid compositions were calculated from arc eruptives, using melt inclusion-based H₂O contents, and were corrected for background mantle contributions. The temperatures, adjusted to h, the vertical depth to the slab beneath the volcanic arc, range from ∼730 to 900°C and agree well (within 30°C on average for each arc) with sub-arc slab surface temperatures predicted by recent thermal models. The coherence between slab model and surface observation implies predominantly vertical transport of fluids within the mantle wedge. Slab surface temperatures are well reconciled with the thermal parameter (the product of slab age and vertical descent rate) andh. Arcs with shallow h (∼80 to 100 km) yield a larger range in slab surface temperature (up to ∼200°C between volcanoes) and more variable magma compositions than arcs with greater h (∼120 to 180 km). This diversity is consistent with coupling of the subducting slab and mantle wedge, and subsequent rapid slab heating, at ∼80 km. Slab surface temperatures at or warmer than the H₂O-saturated solidus suggest that melting at the slab surface is common beneath volcanic arcs. Our results imply that hydrous melts or solute-rich supercritical fluids, and not H₂O-rich aqueous fluids, are thus the agents of mass transport to the mantle wedge.
  • Di Giuseppe, Erika; van Hunen, Jeroen; Funiciello, Francesca; et al. (2007)
    Geochemistry, Geophysics, Geosystems
    Subduction zones are not static features, but trenches retreat (roll back) or advance. Here, we investigate the dominant dynamic controls on trench migration by means of two- and three-dimensional numerical modeling of subduction. This investigation has been carried out by systematically varying the geometrical and rheological model parameters. Our viscoplastic models illustrate that advancing style subduction is promoted by a thick plate, a large viscosity ratio between plate and mantle, and a small density contrast between plate and mantle or an intermediate width (w ∼ 1300 km). Advancing slabs dissipate ∼45% to ∼50% of the energy in the system. Thin plates with relatively low viscosity or relatively high density, or wide slabs (w ∼ 2300 km), on the other hand, promote subduction in the retreating style (i.e., slab roll-back). The energy dissipated by a retreating slab is ∼35% to ∼40% of the total dissipated energy. Most of the energy dissipation occurs in the mantle to accommodate the slab motion, whereas the lithosphere dissipates the remaining part to bend and “unbend.” With a simple scaling law we illustrate that this complex combination of model parameters influencing trench migration can be reduced to a single one: plate stiffness. Stiffer slabs cause the trench to advance, whereas more flexible slabs lead to trench retreat. The reason for this is that all slabs will bend into the subduction zone because of their low plastic strength near the surface, but stiff slabs have more difficulty “unbending” at depth, when arriving at the 660-km discontinuity. Those bent slabs tend to cause the trench to advance. In a similar way, variation of the viscoplasticity parameters in the plate may change the style of subduction: a low value of friction coefficient weakens the plate and results in a retreating style, while higher values strengthen the plate and promote the advancing subduction style. Given the fact that also on Earth the oldest (and therefore probably stiffest) plates have the fastest advancing trenches, we hypothesize that the ability of slabs to unbend after subduction forms the dominant control on trench migration.
  • Tewksbury-Christle, Carolyn M.; Behr, Whitney M.; Helper, Mark A. (2021)
    Geochemistry, Geophysics, Geosystems
    The architecture and mechanical properties of the subduction interface impact large‐scale subduction processes, including mass and volatile recycling, upper‐plate orogenesis, and seismic behavior. The nature of the deep subduction interface, where a dominantly frictional megathrust likely transitions to a distributed ductile shear zone, is poorly understood, due to a lack of constraints on rock types, strain distribution, and interface thickness in this depth range. We characterized these factors in the Condrey Mountain Schist, a Late Jurassic to Early Cretaceous subduction complex in northern California that consists of an upper and lower unit. The Lower Condrey unit is predominantly pelagic and hemipelagic metasediment with m‐to km‐scale metamafic and metaserpentinitic ultramafic lenses all deformed at epidote blueschist facies (0.7–1.1 GPa, 450°C). Major and trace element geochemistry suggest tectonic erosion of the overriding plate sourced all ultramafic and some mafic lenses. We identified two major ductile thrust zones responsible for Lower Condrey unit assembly, with earlier strain distributed across the structural thickness between the ductile thrusts. The Lower Condrey unit records distributed deformation across a sediment‐dominated, 2+ km thick shear zone, possibly consistent with low velocity zones observed in modern subduction zones, despite subducting along a sediment poor, tectonically erosive margin. Periodic strain localization occurred when rheological heterogeneities (i.e., km‐scale ultramafic lenses) entered the interface, facilitating underplating that preserved 10%–60% of the incoming sediment. Modern mass and volatile budgets do not account for erosive margin underplating, so improved quantification is crucial for predicting mass and volatile net flux to Earth′s interior.
  • Khan, Dean; Liebske, Christian; Connolly, James (2021)
    Geochemistry, Geophysics, Geosystems
    The compilation of thermodynamic models for geophysical applications is such a tedious and complex process that it is generally impractical for researchers to refit parameters in existing models in light of new constraints. To mitigate this difficulty, we develop a Bayesian algorithm that permits the modification of a thermodynamic model to account for additional observational constraints. This algorithm can be applied to any thermodynamic dataset and can utilize a wide variety of experimental constraints. To demonstrate the applicability of the algorithm it is used to revise the Stixrude and Lithgow-Bertelloni (2011, https://doi.org/10.1111/j.1365-246x.2010.04890.x), whole-mantle terrestrial thermodynamic model, using phase equilibrium constraints provided by Bertka and Fei (1997, https://doi.org/10.1029/96jb03270), for the more iron-rich compositions that are thought to be relevant to the Martian mantle. The revised thermodynamic model provides a more reliable prediction of phase equilibria in the Martian mantle. Seismic properties are calculated in an internally self-consistent manner along hot and cold areotherms to constrain the upper and lower bounds of these properties for different bulk silicate Mars compositional models.
  • Harrison, D.; Burnard, P.G.; Trieloff, M.; et al. (2003)
    Geochemistry, Geophysics, Geosystems
Publications 1 - 10 of 89