1
|
Vezinet A, Chugunov AV, Sobolev AV, Jain C, Sobolev SV, Batanova VG, Asafov EV, Koshlyakova AN, Arndt NT, Danyushevsky LV, Valley JW. Growth of continental crust and lithosphere subduction in the Hadean revealed by geochemistry and geodynamics. Nat Commun 2025; 16:3850. [PMID: 40280901 PMCID: PMC12032353 DOI: 10.1038/s41467-025-59024-6] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/31/2024] [Accepted: 04/09/2025] [Indexed: 04/29/2025] Open
Abstract
The rates of continental crust growth and recycling on early Earth remain unclear due to the lack of information resulting from the extensive alteration of ancient rocks. Melt inclusions trapped and shielded from alteration in Archean high-Mg olivine crystals offer a solution to this problem. We report an unprecedented unradiogenic Sr mantle source component (87Sr/86Sr = 0.69932 ± 0.00024, 95% confidence interval) of melts included in olivine from 3.27 Ga komatiitic lava flows in the Barberton Greenstone Belt, South Africa. This component indicates a model age of 4.31 ± 0.19 Ga and significant chemical fractionation (Nb/U = 36.9 ± 1.5, Ce/Pb=16.7 ± 1.1), suggesting up to 80% ± 16% of the present-day continental crust's mass was extracted by the late Hadean from the whole mantle. Geodynamic models support this finding, explaining geochemical data by producing 40% to 70% of the present-day continental crust mass during the Hadean in a variable tectonic regime with tens of millions of years-long periods of massive impulsive subduction induced by mantle plumes.
Collapse
Affiliation(s)
- Adrien Vezinet
- Univ. Grenoble Alpes, Univ. Savoie Mont Blanc, CNRS, IRD, Univ. Gustave Eiffel, ISTerre, Grenoble, France.
| | - Aleksandr V Chugunov
- Univ. Grenoble Alpes, Univ. Savoie Mont Blanc, CNRS, IRD, Univ. Gustave Eiffel, ISTerre, Grenoble, France.
| | - Alexander V Sobolev
- Univ. Grenoble Alpes, Univ. Savoie Mont Blanc, CNRS, IRD, Univ. Gustave Eiffel, ISTerre, Grenoble, France.
| | - Charitra Jain
- GFZ Helmholtz Centre for Geosciences, Geodynamic Modelling Section, Potsdam, Germany
| | - Stephan V Sobolev
- GFZ Helmholtz Centre for Geosciences, Geodynamic Modelling Section, Potsdam, Germany.
- University of Potsdam, Institute of Geosciences, Potsdam, Germany.
| | - Valentina G Batanova
- Univ. Grenoble Alpes, Univ. Savoie Mont Blanc, CNRS, IRD, Univ. Gustave Eiffel, ISTerre, Grenoble, France
| | - Evgeny V Asafov
- Vernadsky Institute of Geochemistry and Analytical Chemistry, Russian Academy of Sciences, Moscow, Russia
| | - Alina N Koshlyakova
- Vernadsky Institute of Geochemistry and Analytical Chemistry, Russian Academy of Sciences, Moscow, Russia
| | - Nicholas T Arndt
- Univ. Grenoble Alpes, Univ. Savoie Mont Blanc, CNRS, IRD, Univ. Gustave Eiffel, ISTerre, Grenoble, France
| | | | - John W Valley
- WiscSIMS Lab, Dept. of Geoscience, Univ. of Wisconsin, Madison, WI, USA
| |
Collapse
|
2
|
Erdős Z, Huismans RS, Wolf SG, Faccenna C. Terrane accretion explains thin and hot ocean-continent back-arcs. SCIENCE ADVANCES 2025; 11:eadq8444. [PMID: 40279427 PMCID: PMC12024691 DOI: 10.1126/sciadv.adq8444] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 06/03/2024] [Accepted: 03/21/2025] [Indexed: 04/27/2025]
Abstract
The origin of hot ocean-continent back-arc regions with very thin mantle lithosphere and very high surface heat flow in both extensional and contractional ocean-continent subduction systems is highly enigmatic and unresolved. These first-order characteristics have often been explained with either convective mantle lithosphere removal or by back-arc extension. However, it is unclear what may cause the proposed convective thinning and/or delamination of eclogitic lower crust over very wide regions, whereas back-arc extension is either not observed or insufficient to explain the observed very thin mantle lithosphere. Notably, many of these ocean-continent systems have a long history of terrane accretion. Here, we show, using thermomechanical model experiments, that terrane accretion provides a consistent explanation for the observed key characteristics and naturally leads to rheologically weak back-arcs with continental crust directly on top of the hot sublithospheric mantle.
Collapse
Affiliation(s)
- Zoltán Erdős
- HUN-REN Institute of Earth Physics and Space Science (HUN-REN EPSS), Sopron, Hungary
- Department of Geophysics and Space Science, Eötvös Loránd University, Budapest, Hungary
- GFZ Helmholtz Centre for Geosciences, Potsdam, Germany
| | | | | | - Claudio Faccenna
- GFZ Helmholtz Centre for Geosciences, Potsdam, Germany
- Department of Science, University Roma Tre, Roma, Italy
| |
Collapse
|
3
|
Tobie G, Auclair-Desrotour P, Běhounková M, Kervazo M, Souček O, Kalousová K. Tidal Deformation and Dissipation Processes in Icy Worlds. SPACE SCIENCE REVIEWS 2025; 221:6. [PMID: 39830012 PMCID: PMC11739232 DOI: 10.1007/s11214-025-01136-y] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 12/28/2023] [Accepted: 12/30/2024] [Indexed: 01/22/2025]
Abstract
Tidal interactions play a key role in the dynamics and evolution of icy worlds. The intense tectonic activity of Europa and the eruption activity on Enceladus are clear examples of the manifestation of tidal deformation and associated dissipation. While tidal heating has long been recognized as a major driver in the activity of these icy worlds, the mechanism controlling how tidal forces deform the different internal layers and produce heat by tidal friction still remains poorly constrained. As tidal forcing varies with orbital characteristics (distance to the central planet, eccentricity, obliquity), the contribution of tidal heating to the internal heat budget can strongly change over geological timescales. In some circumstances, the tidally-produced heat can result in internal melting and surface activity taking various forms. Even in the absence of significant heat production, tidal deformation can be used to probe the interior structure, the tidal response of icy moons being strongly sensitive to their hydrosphere structure. In the present paper, we review the methods to compute tidal deformation and dissipation in the different layers composing icy worlds. After summarizing the main principle of tidal deformation and the different rheological models used to model visco-elastic tidal response, we describe the dissipation processes expected in rock-dominated cores, subsurface oceans and icy shells and highlight the potential effects of tidal heating in terms of thermal evolution and activity. We finally anticipate how data collected by future missions to Jupiter's and Saturn's moons could be used to constrain their tidal response and the consequences for past and present activities.
Collapse
Affiliation(s)
- G. Tobie
- Laboratoire de Planétologie et Géosciences, UMR 6112, CNRS, Nantes Université, Université d’Angers, Le Mans Université, Nantes, France
| | - P. Auclair-Desrotour
- IMCCE, CNRS, Observatoire de Paris, PSL University, Sorbonne Université, Paris, France
| | - M. Běhounková
- Faculty of Mathematics and Physics, Department of Geophysics, Charles University, V Holesšovičkách 2, Praha, Praha 8 180 00 Czech Republic
| | - M. Kervazo
- Laboratoire de Planétologie et Géosciences, UMR 6112, CNRS, Nantes Université, Université d’Angers, Le Mans Université, Nantes, France
| | - O. Souček
- Faculty of Mathematics and Physics, Mathematical Institute, Charles University, Sokolovská 83, Praha, Praha 8 186 75 Czech Republic
| | - K. Kalousová
- Faculty of Mathematics and Physics, Department of Geophysics, Charles University, V Holesšovičkách 2, Praha, Praha 8 180 00 Czech Republic
| |
Collapse
|
4
|
Paul J, Golabek GJ, Rozel AB, Tackley PJ, Katsura T, Fei H. Effect of bridgmanite-ferropericlase grain size evolution on Earth's average mantle viscosity: implications for mantle convection in early and present-day Earth. PROGRESS IN EARTH AND PLANETARY SCIENCE 2024; 11:64. [PMID: 39640495 PMCID: PMC11615032 DOI: 10.1186/s40645-024-00658-3] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 06/19/2024] [Accepted: 10/10/2024] [Indexed: 12/07/2024]
Abstract
Recent experimental investigations of grain size evolution in bridgmanite-ferropericlase assemblages have suggested very slow growth for these bimodal phases. Despite numerous speculations on grain size-dependent viscosity, a comprehensive test with realistic grain size evolution parameters compatible with the lower mantle has been lacking. In this study, we develop self-consistent 2-D spherical half-annulus geodynamic models of Earth's evolution using the finite volume code StagYY to assess the role of grain size on lower mantle viscosity. We explore several models with and without grain size evolution to compare their effects on mantle viscosity. In models with grain size evolution, we consider three scenarios: (1) uniform grain growth throughout the entire mantle with a composite rheology, (2) different grain growth in the upper and lower mantle with a composite rheology, and (3) different grain growth in the upper and lower mantle with purely diffusion creep rheology. In the case of different grain size evolution, the upper mantle's grain size evolution law is controlled by forsterite-enstatite grain growth, while the lower mantle's grain size evolution law is controlled by bridgmanite-ferropericlase grain growth. Our results suggest that mantle viscosity is primarily controlled by temperature, whereas grain size has a minor effect compared to the effect of temperature. We attribute two primary reasons for this: First, the bridgmanite-ferropericlase growth is very slow in the lower mantle and the grain size variation is too small to significantly alter the mantle viscosity. Secondly, if grains grow too fast, thus the mantle deforms in the dislocation creep regime, making viscosity grain size-independent. To establish the robustness of this finding we vary several other model parameters, such as surface yield strength, phase transition grain size reset, different transitional stresses for creep mechanisms, pressure dependence on grain growth, and different grain damage parameters. For all our models, we consistently find that grain size has a very limited effect on controlling lower mantle viscosity in the present-day Earth. However, large grain size may have affected the lower mantle viscosity in the early Earth as larger grains of single phase bridgmanite could increase the viscosity of the early mantle delaying the onset of global convection.
Collapse
Affiliation(s)
- Jyotirmoy Paul
- Bayerisches Geoinstitut, University of Bayreuth, Universitätsstraße 30, 95447 Bayreuth, Germany
- Centre for Planetary Habitability, University of Oslo, 0316 Oslo, Norway
| | - Gregor J. Golabek
- Bayerisches Geoinstitut, University of Bayreuth, Universitätsstraße 30, 95447 Bayreuth, Germany
| | - Antoine B. Rozel
- Institute of Geophysics, Department of Earth and Planetary Sciences, ETH Zürich, Sonneggstrasse 5, 8092 Zürich, Switzerland
| | - Paul J. Tackley
- Institute of Geophysics, Department of Earth and Planetary Sciences, ETH Zürich, Sonneggstrasse 5, 8092 Zürich, Switzerland
| | - Tomoo Katsura
- Bayerisches Geoinstitut, University of Bayreuth, Universitätsstraße 30, 95447 Bayreuth, Germany
| | - Hongzhan Fei
- School of Earth Sciences, Zhejiang University, Hangzhou, China
| |
Collapse
|
5
|
Yuan X, Li Y, Brune S, Li K, Pons M, Wolf SG. Coordination between deformation, precipitation, and erosion during orogenic growth. Nat Commun 2024; 15:10362. [PMID: 39609430 PMCID: PMC11604927 DOI: 10.1038/s41467-024-54690-4] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/19/2024] [Accepted: 11/18/2024] [Indexed: 11/30/2024] Open
Abstract
Crustal thickening associated with orogenic growth elevates topography, causing orographic enhancement of precipitation, which in turn facilitates local erosion and possibly intensifies localization of deformation. How these three processes-deformation, precipitation, and erosion-coordinate during orogenic growth remains unknown. Here, we present a numerical model where tectonics, surface processes, and orographic precipitation are tightly coupled, and explore the impact on low, intermediate, and high erodibility orogens. We show that, for intermediate erodibility models, rock uplift rates and precipitation rates correlate well with erosion rates during the formation of orogenic plateaus with high correlation coefficients of ~0.9 between rock uplift and erosion rates, and ~0.8 between precipitation and erosion rates. We demonstrate a cyclicity of correlation evolution among uplift, precipitation, and erosion rates through the development of new faults propagating outward. These results shed insights into the relative tectonic or climatic control on erosion in active orogens (e.g., Himalayas, Central Andes, and Southern Alps of New Zealand), and provide a plausible explanation for several conflicting data and interpretations in the Himalayas, which depend on the stage of maturity of the newest fault and the relative locations to old faults.
Collapse
Affiliation(s)
- Xiaoping Yuan
- State Key Laboratory of Biogeology and Environmental Geology, Hubei Key Laboratory of Critical Zone Evolution, School of Earth Sciences, China University of Geosciences, Wuhan, China.
| | - Yuqiang Li
- State Key Laboratory of Biogeology and Environmental Geology, Hubei Key Laboratory of Critical Zone Evolution, School of Earth Sciences, China University of Geosciences, Wuhan, China
| | - Sascha Brune
- Helmholtz Centre Potsdam - GFZ German Research Centre for Geosciences, Potsdam, Germany
- Institute of Geosciences, University of Potsdam, Potsdam, Germany
| | - Kai Li
- Helmholtz Centre Potsdam - GFZ German Research Centre for Geosciences, Potsdam, Germany
- Institute of Geological Sciences, Freie Universität Berlin, Berlin, Germany
| | - Michaël Pons
- Helmholtz Centre Potsdam - GFZ German Research Centre for Geosciences, Potsdam, Germany
| | - Sebastian G Wolf
- Helmholtz Centre Potsdam - GFZ German Research Centre for Geosciences, Potsdam, Germany
- Department of Earth Science, University of Bergen, Bergen, Norway
| |
Collapse
|
6
|
Li Y, Gurnis M. Rapid shear zone weakening during subduction initiation. Proc Natl Acad Sci U S A 2024; 121:e2404939121. [PMID: 39453749 PMCID: PMC11536125 DOI: 10.1073/pnas.2404939121] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/09/2024] [Accepted: 09/12/2024] [Indexed: 10/27/2024] Open
Abstract
Subduction zones play a pivotal role in the mechanics of plate tectonics by providing the driving force through slab pull and weak megathrusts that facilitate the relative motion between tectonic plates. The initiation of subduction zones is intricately linked to the accumulation of slab pull and development of weakness at plate boundaries and, by consequence, the largest changes in the energetics of mantle convection. However, the transient nature of subduction initiation accompanied by intense subsequent tectonic activity, leaves critical evidence poorly preserved and making subduction initiation difficult to constrain. We overcome these limitations through a comprehensive analysis focused on Puysegur, a well-constrained extant example of subduction initiation offshore South Island, New Zealand. Through time-dependent, three-dimensional thermo-mechanical computations and quantitative comparison to new geophysical and geological observations, including topography, stratigraphy, and seismicity, we demonstrate that subduction initiation develops with a fast strain weakening described with a small characteristic displacement ([Formula: see text] 4 to 8 km). Potential physical mechanisms contributing to the strain weakening are explored and we find that the observed fast weakening may arise through a combination of grain-size reduction within the lower lithosphere and fluid pressurization at shallower depths. With the shared commonality in the underlying physics of tectonic processes, the rapid strain weakening constrained at Puysegur offers insights into the formation of the first subduction during early Earth and the onset of plate tectonics.
Collapse
Affiliation(s)
- Yida Li
- Division of Geological and Planetary Sciences, Seismological Laboratory, California Institute of Technology, Pasadena, CA91125
| | - Michael Gurnis
- Division of Geological and Planetary Sciences, Seismological Laboratory, California Institute of Technology, Pasadena, CA91125
| |
Collapse
|
7
|
Gernon TM, Hincks TK, Brune S, Braun J, Jones SM, Keir D, Cunningham A, Glerum A. Coevolution of craton margins and interiors during continental break-up. Nature 2024; 632:327-335. [PMID: 39112622 PMCID: PMC11306106 DOI: 10.1038/s41586-024-07717-1] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/25/2022] [Accepted: 06/13/2024] [Indexed: 08/10/2024]
Abstract
Many cratonic continental fragments dispersed during the rifting and break-up of Gondwana are bound by steep topographic landforms known as 'great escarpments'1-4, which rim elevated plateaus in the craton interior5,6. In terms of formation, escarpments and plateaus are traditionally considered distinct owing to their spatial separation, occasionally spanning more than a thousand kilometres. Here we integrate geological observations, statistical analysis, geodynamic simulations and landscape-evolution models to develop a physical model that mechanistically links both phenomena to continental rifting. Escarpments primarily initiate at rift-border faults and slowly retreat at about 1 km Myr-1 through headward erosion. Simultaneously, rifting generates convective instabilities in the mantle7-10 that migrate cratonward at a faster rate of about 15-20 km Myr-1 along the lithospheric root, progressively removing cratonic keels11, driving isostatic uplift of craton interiors and forming a stable, elevated plateau. This process forces a synchronized wave of denudation, documented in thermochronology studies, which persists for tens of millions of years and migrates across the craton at a comparable or slower pace. We interpret the observed sequence of rifting, escarpment formation and exhumation of craton interiors as an evolving record of geodynamic mantle processes tied to continental break-up, upending the prevailing notion of cratons as geologically stable terrains.
Collapse
Affiliation(s)
- Thomas M Gernon
- School of Ocean & Earth Science, University of Southampton, Southampton, UK.
| | - Thea K Hincks
- School of Ocean & Earth Science, University of Southampton, Southampton, UK
| | - Sascha Brune
- Helmholtz Centre Potsdam - GFZ German Research Centre for Geosciences, Potsdam, Germany
- University of Potsdam, Potsdam, Germany
| | - Jean Braun
- Helmholtz Centre Potsdam - GFZ German Research Centre for Geosciences, Potsdam, Germany
- University of Potsdam, Potsdam, Germany
| | - Stephen M Jones
- School of Geography, Earth and Environmental Sciences, University of Birmingham, Birmingham, UK
| | - Derek Keir
- School of Ocean & Earth Science, University of Southampton, Southampton, UK
- Dipartimento di Scienze della Terra, Universita degli Studi di Firenze, Florence, Italy
| | - Alice Cunningham
- School of Ocean & Earth Science, University of Southampton, Southampton, UK
| | - Anne Glerum
- Helmholtz Centre Potsdam - GFZ German Research Centre for Geosciences, Potsdam, Germany
| |
Collapse
|
8
|
Rhoden AR, Ferguson SN, Bottke W, Castillo-Rogez JC, Martin E, Bland M, Kirchoff M, Zannoni M, Rambaux N, Salmon J. Geologic Constraints on the Formation and Evolution of Saturn's Mid-Sized Moons. SPACE SCIENCE REVIEWS 2024; 220:55. [PMID: 39036784 PMCID: PMC11255024 DOI: 10.1007/s11214-024-01084-z] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 04/05/2023] [Accepted: 06/13/2024] [Indexed: 07/23/2024]
Abstract
Saturn's mid-sized icy moons have complex relationships with Saturn's interior, the rings, and with each other, which can be expressed in their shapes, interiors, and geology. Observations of their physical states can, thus, provide important constraints on the ages and formation mechanism(s) of the moons, which in turn informs our understanding of the formation and evolution of Saturn and its rings. Here, we describe the cratering records of the mid-sized moons and the value and limitations of their use for constraining the histories of the moons. We also discuss observational constraints on the interior structures of the moons and geologically-derived inferences on their thermal budgets through time. Overall, the geologic records of the moons (with the exception of Mimas) include evidence of epochs of high heat flows, short- and long-lived subsurface oceans, extensional tectonics, and considerable cratering. Curiously, Mimas presents no clear evidence of an ocean within its surface geology, but its rotation and orbit indicate a present-day ocean. While the moons need not be primordial to produce the observed levels of interior evolution and geologic activity, there is likely a minimum age associated with their development that has yet to be determined. Uncertainties in the populations impacting the moons makes it challenging to further constrain their formation timeframes using craters, whereas the characteristics of their cores and other geologic inferences of their thermal evolutions may help narrow down their potential histories. Disruptive collisions may have also played an important role in the formation and evolution of Saturn's mid-sized moons, and even the rings of Saturn, although more sophisticated modeling is needed to determine the collision conditions that produce rings and moons that fit the observational constraints. Overall, the existence and physical characteristics of Saturn's mid-sized moons provide critical benchmarks for the development of formation theories.
Collapse
Affiliation(s)
| | | | - William Bottke
- Southwest Research Institute, 1050 Walnut St, Boulder, CO 80302 USA
| | | | - Emily Martin
- Center for Earth and Planetary Studies, National Air and Space Museum, Smithsonian Institution, Washington, DC USA
| | - Michael Bland
- U.S. Geological Survey, Astrogeology Science Center, Flagstaff, AZ USA
| | | | - Marco Zannoni
- Dipartimento di Ingegneria Industriale, Alma Mater Studiorum – Università di Bologna, Forlì, Italy
| | - Nicolas Rambaux
- IMCCE, CNRS, Observatoire de Paris, PSL Université, Sorbonne Université, Université de Lille 1, UMR 8028 du CNRS, 77 Denfert-Rochereau, 75014 Paris, France
| | - Julien Salmon
- Southwest Research Institute, 1050 Walnut St, Boulder, CO 80302 USA
| |
Collapse
|
9
|
Hu J, Rudi J, Gurnis M, Stadler G. Constraining Earth's nonlinear mantle viscosity using plate-boundary resolving global inversions. Proc Natl Acad Sci U S A 2024; 121:e2318706121. [PMID: 38968110 PMCID: PMC11252765 DOI: 10.1073/pnas.2318706121] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/10/2023] [Accepted: 05/20/2024] [Indexed: 07/07/2024] Open
Abstract
Variable viscosity in Earth's mantle exerts a fundamental control on mantle convection and plate tectonics, yet rigorously constraining the underlying parameters has remained a challenge. Inverse methods have not been sufficiently robust to handle the severe viscosity gradients and nonlinearities (arising from dislocation creep and plastic failure) while simultaneously resolving the megathrust and bending slabs globally. Using global plate motions as constraints, we overcome these challenges by combining a scalable nonlinear Stokes solver that resolves the key tectonic features with an adjoint-based Bayesian approach. Assuming plate cooling, variations in the thickness of continental lithosphere, slabs, and broad scale lower mantle structure as well as a constant grain size through the bulk of the upper mantle, a good fit to global plate motions is found with a nonlinear upper mantle stress exponent of 2.43 [Formula: see text] 0.25 (mean [Formula: see text] SD). A relatively low yield stress of 151 [Formula: see text] 19 MPa is required for slabs to bend during subduction and transmit a slab pull that generates asymmetrical subduction. The recovered long-term strength of megathrusts (plate interfaces) varies between different subduction zones, with South America having a larger strength and Vanuatu and Central America having lower values with important implications for the stresses driving megathrust earthquakes.
Collapse
Affiliation(s)
- Jiashun Hu
- Department of Earth and Space Sciences, Southern University of Science and Technology, Shenzhen518055, China
| | - Johann Rudi
- Department of Mathematics, Virginia Tech, Blacksburg, VA24061
| | - Michael Gurnis
- Seismological Laboratory, Division of Geological and Planetary Sciences, California Institute of Technology, Pasadena, CA91125
| | - Georg Stadler
- Courant Institute of Mathematical Sciences, New York University, New York, NY10012
| |
Collapse
|
10
|
Lee C, Kim Y. Understanding subduction infancy to mature subduction in Southwest Japan via the self-consistent formation of a weak slab interface. Sci Rep 2023; 13:21425. [PMID: 38052949 DOI: 10.1038/s41598-023-48746-6] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/05/2023] [Accepted: 11/29/2023] [Indexed: 12/07/2023] Open
Abstract
The weak slab interface controls long-term subduction dynamics. A weak hydrous layer at the slab interface promotes mechanical decoupling between the forearc mantle and the subducting slab and converts a hot forearc mantle to a cold mantle. Often referred to as a cold nose, the cold forearc mantle, plays a key role in the transition from subduction infancy to mature subduction. This study was the first to numerically demonstrate the self-consistent formation of a weak hydrous layer with permeability anisotropy based on the Southwest Japan subduction zone case, where transition-related geological features were present. Our models showed that mechanical decoupling by spontaneous downdip growth of the weak hydrous layer created a cold nose by converting a hot forearc mantle to a cold mantle. The emergence of a cold nose explained the migration of the forearc-to-arc volcanic front, expressed as the formation of mid-Miocene forearc high-magnesium andesite and Quaternary arc adakite. Furthermore, the weak hydrous layer providing a pathway for free-water transport toward the mantle wedge tip elucidates slab/mantle-derived geochemical components in deep groundwater as well as large S-wave delay times and non-volcanic seismic tremors in the forearc.
Collapse
Affiliation(s)
- Changyeol Lee
- Department of Earth System Sciences, Yonsei University, 50 Yonsei-ro Seodaemun-gu, Seoul, 03722, Republic of Korea
| | - YoungHee Kim
- School of Earth and Environmental Sciences, Seoul National University, 1 Gwanak-ro, Gwanak-gu, Seoul, 08826, Republic of Korea.
| |
Collapse
|
11
|
Wu X, Hu J, Chen L, Liu L, Liu L. Paleogene India-Eurasia collision constrained by observed plate rotation. Nat Commun 2023; 14:7272. [PMID: 37949864 PMCID: PMC10638303 DOI: 10.1038/s41467-023-42920-0] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/29/2023] [Accepted: 10/24/2023] [Indexed: 11/12/2023] Open
Abstract
The Cenozoic India-Eurasia collision has had profound impacts on shaping the Tibetan plateau, but its early history remains controversial due to uneven availability of constraints. Recent plate reconstructions reveal two prominent counterclockwise rotation (azimuthal change) rate peaks of the Indian plate at 52-44 and 33-20 Ma, respectively, which could bear key information about this collision history. Using fully dynamic three-dimensional numerical modeling, we show that the first rotation rate peak reflected the initial diachronous collision from the western-central to eastern Indian front, and the second peak reflected the full collision leading to strong coupling between India and Eurasia. Further comparison with observation suggests that the initial and complete India-Eurasia collision likely occurred at 55 ± 5 and 40 ± 5 Ma, respectively, an inference consistent with key geological observations. We suggest that this collision history is instructive for studying the tectonic history of the Tibetan plateau and its surrounding areas.
Collapse
Affiliation(s)
- Xiaoyue Wu
- State Key Laboratory of Lithospheric Evolution, Institute of Geology and Geophysics, Chinese Academy of Sciences, 100029, Beijing, China
- Department of Earth and Space Sciences, Southern University of Science and Technology, 518055, Shenzhen, China
- College of Earth and Planetary Sciences, University of Chinese Academy of Sciences, 100049, Beijing, China
| | - Jiashun Hu
- Department of Earth and Space Sciences, Southern University of Science and Technology, 518055, Shenzhen, China.
| | - Ling Chen
- State Key Laboratory of Lithospheric Evolution, Institute of Geology and Geophysics, Chinese Academy of Sciences, 100029, Beijing, China
- College of Earth and Planetary Sciences, University of Chinese Academy of Sciences, 100049, Beijing, China
| | - Liang Liu
- State Key Laboratory of Isotope Geochemistry, Guangzhou Institute of Geochemistry, Chinese Academy of Sciences, 510640, Guangzhou, China
| | - Lijun Liu
- State Key Laboratory of Lithospheric Evolution, Institute of Geology and Geophysics, Chinese Academy of Sciences, 100029, Beijing, China.
- Department of Earth Science & Environmental Change, University of Illinois at Urbana-Champaign, Champaign, IL, 61820, USA.
| |
Collapse
|
12
|
Samuel H, Drilleau M, Rivoldini A, Xu Z, Huang Q, Garcia RF, Lekić V, Irving JCE, Badro J, Lognonné PH, Connolly JAD, Kawamura T, Gudkova T, Banerdt WB. Geophysical evidence for an enriched molten silicate layer above Mars's core. Nature 2023; 622:712-717. [PMID: 37880437 PMCID: PMC10600000 DOI: 10.1038/s41586-023-06601-8] [Citation(s) in RCA: 4] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/10/2022] [Accepted: 08/31/2023] [Indexed: 10/27/2023]
Abstract
The detection of deep reflected S waves on Mars inferred a core size of 1,830 ± 40 km (ref. 1), requiring light-element contents that are incompatible with experimental petrological constraints. This estimate assumes a compositionally homogeneous Martian mantle, at odds with recent measurements of anomalously slow propagating P waves diffracted along the core-mantle boundary2. An alternative hypothesis is that Mars's mantle is heterogeneous as a consequence of an early magma ocean that solidified to form a basal layer enriched in iron and heat-producing elements. Such enrichment results in the formation of a molten silicate layer above the core, overlain by a partially molten layer3. Here we show that this structure is compatible with all geophysical data, notably (1) deep reflected and diffracted mantle seismic phases, (2) weak shear attenuation at seismic frequency and (3) Mars's dissipative nature at Phobos tides. The core size in this scenario is 1,650 ± 20 km, implying a density of 6.5 g cm-3, 5-8% larger than previous seismic estimates, and can be explained by fewer, and less abundant, alloying light elements than previously required, in amounts compatible with experimental and cosmochemical constraints. Finally, the layered mantle structure requires external sources to generate the magnetic signatures recorded in Mars's crust.
Collapse
Affiliation(s)
- Henri Samuel
- Université Paris Cité, Institut de physique du globe de Paris, CNRS, Paris, France.
| | - Mélanie Drilleau
- Institut Supérieur de l'Aéronautique et de l'Espace ISAE-SUPAERO, Toulouse, France
| | | | - Zongbo Xu
- Université Paris Cité, Institut de physique du globe de Paris, CNRS, Paris, France
| | - Quancheng Huang
- Department of Geophysics, Colorado School of Mines, Golden, CO, USA
- University of Maryland, College Park, MD, USA
| | - Raphaël F Garcia
- Institut Supérieur de l'Aéronautique et de l'Espace ISAE-SUPAERO, Toulouse, France
| | | | | | - James Badro
- Université Paris Cité, Institut de physique du globe de Paris, CNRS, Paris, France
| | - Philippe H Lognonné
- Université Paris Cité, Institut de physique du globe de Paris, CNRS, Paris, France
| | | | - Taichi Kawamura
- Université Paris Cité, Institut de physique du globe de Paris, CNRS, Paris, France
| | - Tamara Gudkova
- Schmidt Institute of Physics of the Earth, Russian Academy of Sciences, Moscow, Russia
| | - William B Banerdt
- Jet Propulsion Laboratory, California Institute of Technology, Pasadena, CA, USA
| |
Collapse
|
13
|
Cathles L, Fjeldskar W, Lenardic A, Romanowicz B, Seales J, Richards M. Influence of the asthenosphere on earth dynamics and evolution. Sci Rep 2023; 13:13367. [PMID: 37591899 PMCID: PMC10435468 DOI: 10.1038/s41598-023-39973-y] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/31/2023] [Accepted: 08/02/2023] [Indexed: 08/19/2023] Open
Abstract
The existence of a thin, weak asthenospheric layer beneath Earth's lithospheric plates is consistent with existing geological and geophysical constraints, including Pleistocene glacio-isostatic adjustment, modeling of gravity anomalies, studies of seismic anisotropy, and post-seismic rebound. Mantle convection models suggest that a pronounced weak zone beneath the upper thermal boundary layer (lithosphere) may be essential to the plate tectonic style of convection found on Earth. The asthenosphere is likely related to partial melting and the presence of water in the sub-lithospheric mantle, further implying that the long-term evolution of the Earth may be controlled by thermal regulation and volatile recycling that maintain a geotherm that approaches the wet mantle solidus at asthenospheric depths.
Collapse
Affiliation(s)
- Lawrence Cathles
- Department of Earth and Atmospheric Sciences, Cornell University, Ithaca, USA
| | | | | | - Barbara Romanowicz
- Department of Earth and Planetary Science, University of California, Berkeley, USA
| | - Johnny Seales
- Department of Earth Science, Rice University, Houston, USA
| | - Mark Richards
- Department of Earth and Space Sciences, University of Washington, Seattle, USA
| |
Collapse
|
14
|
Briaud A, Ganino C, Fienga A, Mémin A, Rambaux N. The lunar solid inner core and the mantle overturn. Nature 2023; 617:743-746. [PMID: 37138080 DOI: 10.1038/s41586-023-05935-7] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/19/2022] [Accepted: 03/08/2023] [Indexed: 05/05/2023]
Abstract
Seismological models from Apollo missions provided the first records of the Moon inner structure with a decrease in seismic wave velocities at the core-mantle boundary1-3. The resolution of these records prevents a strict detection of a putative lunar solid inner core and the impact of the lunar mantle overturn in the lowest part of the Moon is still discussed4-7. Here we combine geophysical and geodesic constraints from Monte Carlo exploration and thermodynamical simulations for different Moon internal structures to show that only models with a low viscosity zone enriched in ilmenite and an inner core present densities deduced from thermodynamic constraints compatible with densities deduced from tidal deformations. We thus obtain strong indications in favour of the lunar mantle overturn scenario and, in this context, demonstrate the existence of the lunar inner core with a radius of 258 ± 40 km and density 7,822 ± 1,615 kg m-3. Our results question the evolution of the Moon magnetic field thanks to its demonstration of the existence of the inner core and support a global mantle overturn scenario that brings substantial insights on the timeline of the lunar bombardment in the first billion years of the Solar System8.
Collapse
Affiliation(s)
- Arthur Briaud
- Université Côte d'Azur, Observatoire de la Côte d'Azur, CNRS, Géoazur, Valbonne, France.
| | - Clément Ganino
- Université Côte d'Azur, Observatoire de la Côte d'Azur, CNRS, Géoazur, Valbonne, France
| | - Agnès Fienga
- Université Côte d'Azur, Observatoire de la Côte d'Azur, CNRS, Géoazur, Valbonne, France.
- IMCCE, Observatoire de Paris, Sorbonne Université, PSL University, CNRS, Paris, France.
| | - Anthony Mémin
- Université Côte d'Azur, Observatoire de la Côte d'Azur, CNRS, Géoazur, Valbonne, France
| | - Nicolas Rambaux
- IMCCE, Observatoire de Paris, Sorbonne Université, PSL University, CNRS, Paris, France
| |
Collapse
|
15
|
Gianni GM, Likerman J, Navarrete CR, Gianni CR, Zlotnik S. Ghost-arc geochemical anomaly at a spreading ridge caused by supersized flat subduction. Nat Commun 2023; 14:2083. [PMID: 37045842 PMCID: PMC10097660 DOI: 10.1038/s41467-023-37799-w] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/07/2022] [Accepted: 03/27/2023] [Indexed: 04/14/2023] Open
Abstract
The Southern Atlantic-Southwest Indian ridges (SASWIR) host mid-ocean ridge basalts with a residual subduction-related geochemical fingerprint (i.e., a ghost-arc signature) of unclear origin. Here, we show through an analysis of plate kinematic reconstructions and seismic tomography models that the SASWIR subduction-modified mantle source formed in the Jurassic close to the Georgia Islands slab (GI) and remained near-stationary in the mantle reference frame. In this analysis, the GI lies far inboard the Jurassic Patagonian-Antarctic Peninsula active margin. This was formerly attributed to a large-scale flat subduction event in the Late Triassic-Early Jurassic. We propose that during this flat slab stage, the subduction-modified mantle areas beneath the Mesozoic active margin and surrounding sutures zones may have been bulldozed inland by >2280 km. After the demise of the flat slab, this mantle anomaly remained near-stationary and was sampled by the Karoo mantle plume 183 Million years (Myr) ago and again since 55 Myr ago by the SASWIR. We refer to this process as asthenospheric anomaly telescoping. This study provides a hitherto unrecognized geodynamic effect of flat subduction, the viability of which we support through numerical modeling.
Collapse
Affiliation(s)
- Guido M Gianni
- Instituto Geofísico Sismológico Ing. Fernando Volponi (IGSV), Universidad Nacional de San Juan, San Juan, Argentina
- National Scientific and Technical Research Council (CONICET), Capital Federal, Argentina
| | - Jeremías Likerman
- National Scientific and Technical Research Council (CONICET), Capital Federal, Argentina
- Instituto de Estudios Andinos Don Pablo Groeber, Universidad de Buenos Aires, Capital Federal, Argentina
| | - César R Navarrete
- National Scientific and Technical Research Council (CONICET), Capital Federal, Argentina
- Laboratorio Patagónico de Petro-Tectónica, Universidad Nacional de la Patagonia "San Juan Bosco", Comodoro Rivadavia, Chubut, Argentina
| | - Conrado R Gianni
- Instituto Geofísico Sismológico Ing. Fernando Volponi (IGSV), Universidad Nacional de San Juan, San Juan, Argentina
| | - Sergio Zlotnik
- Laboratori de Cálcul Numéric, Escola Técnica Superior d'Enginyers de Camins, Canals i Ports, Universitat Politécnica de Catalunya, Barcelona, Spain.
- Centre Internacional de Métodes Numérics a l'Enginyeria (CIMNE), Barcelona, Spain.
| |
Collapse
|
16
|
Bonnet Gibet V, Michaut C, Wieczorek M, Lognonné P. A Positive Feedback Between Crustal Thickness and Melt Extraction for the Origin of the Martian Dichotomy. JOURNAL OF GEOPHYSICAL RESEARCH. PLANETS 2022; 127:e2022JE007472. [PMID: 37033153 PMCID: PMC10078261 DOI: 10.1029/2022je007472] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 07/25/2022] [Revised: 10/20/2022] [Accepted: 11/17/2022] [Indexed: 06/19/2023]
Abstract
A North/South difference in crustal thickness is likely at the origin of the Martian dichotomy in topography. Recent crustal thickness maps were obtained by inversion of topography and gravity data seismically anchored at the InSight station. On average, the Martian crust is 51-71 km thick with a southern crust thicker by 18-28 km than the northern one. The origin of this crustal dichotomy is still debated although the hypothesis of a large impact is at present very popular. Here, we propose a new mechanism for the formation of this dichotomy that involves a positive feedback between crustal growth and mantle melting. As the crust is enriched in heat-producing elements, the lid of a one-plate planet is hotter and thinner where the crust is thicker, inducing a larger amount of partial melt below the lid and hence a larger rate of melt extraction and crustal growth. We first demonstrate analytically that larger wavelength perturbations, that is, hemispherical perturbations, grow faster because smaller wavelengths are more attenuated by thermal diffusion. We then use a parameterized thermal evolution model with a well-mixed mantle topped by two different lids characterized by their thermal structures and thicknesses to study the growth of the crust in the two hemispheres. Our results demonstrate that this positive feedback can generate a significant crustal dichotomy.
Collapse
Affiliation(s)
- Valentin Bonnet Gibet
- Laboratoire de Géologie de LyonTerre, Planètes, EnvironnementEcole Normale Supérieure de LyonCNRSUniversité de LyonUniversité Claude Bernard Lyon 1Université Jean MonetLyonFrance
| | - Chloé Michaut
- Laboratoire de Géologie de LyonTerre, Planètes, EnvironnementEcole Normale Supérieure de LyonCNRSUniversité de LyonUniversité Claude Bernard Lyon 1Université Jean MonetLyonFrance
- Institut Universitaire de FranceParisFrance
| | - Mark Wieczorek
- Laboratoire LagrangeObservatoire de la Côte d’AzurCNRSUniversité Côte d’AzurNiceFrance
| | - Philippe Lognonné
- CNRSInstitut de physique du globe de ParisUniversité de ParisParisFrance
| |
Collapse
|
17
|
Bahadori A, Holt WE, Austermann J, Campbell L, Rasbury ET, Davis DM, Calvelage CM, Flesch LM. The role of gravitational body forces in the development of metamorphic core complexes. Nat Commun 2022; 13:5646. [PMID: 36163354 PMCID: PMC9513114 DOI: 10.1038/s41467-022-33361-2] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/01/2022] [Accepted: 09/14/2022] [Indexed: 12/01/2022] Open
Abstract
Within extreme continental extension areas, ductile middle crust is exhumed at the surface as metamorphic core complexes. Sophisticated quantitative models of extreme extension predicted upward transport of ductile middle-lower crust through time. Here we develop a general model for metamorphic core complexes formation and demonstrate that they result from the collapse of a mountain belt supported by a thickened crustal root. We show that gravitational body forces generated by topography and crustal root cause an upward flow pattern of the ductile lower-middle crust, facilitated by a detachment surface evolving into low-angle normal fault. This detachment surface acquires large amounts of finite strain, consistent with thick mylonite zones found in metamorphic core complexes. Isostatic rebound exposes the detachment in a domed upwarp, while the final Moho discontinuity across the extended region relaxes to a flat geometry. This work suggests that belts of metamorphic core complexes are a fossil signature of collapsed highlands.
Collapse
Affiliation(s)
- Alireza Bahadori
- Lamont-Doherty Earth Observatory, Columbia University in the City of New York, Palisades, NY, USA.
- Department of Geosciences, Stony Brook University, Stony Brook, NY, USA.
| | - William E Holt
- Department of Geosciences, Stony Brook University, Stony Brook, NY, USA
| | - Jacqueline Austermann
- Lamont-Doherty Earth Observatory, Columbia University in the City of New York, Palisades, NY, USA
| | - Lajhon Campbell
- Department of Geosciences, Stony Brook University, Stony Brook, NY, USA
| | - E Troy Rasbury
- Department of Geosciences, Stony Brook University, Stony Brook, NY, USA
| | - Daniel M Davis
- Department of Geosciences, Stony Brook University, Stony Brook, NY, USA
| | | | - Lucy M Flesch
- Department of Earth, Atmospheric, and Planetary Sciences, Purdue University, West Lafayette, IN, USA
| |
Collapse
|
18
|
Peltier WR, Wu PPC, Argus DF, Li T, Velay-Vitow J. Glacial isostatic adjustment: physical models and observational constraints. REPORTS ON PROGRESS IN PHYSICS. PHYSICAL SOCIETY (GREAT BRITAIN) 2022; 85:096801. [PMID: 35820343 DOI: 10.1088/1361-6633/ac805b] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 11/18/2021] [Accepted: 07/12/2022] [Indexed: 06/15/2023]
Abstract
By far the most prescient insights into the interior structure of the planet have been provided on the basis of elastic wave seismology. Analysis of the travel times of shear or compression wave phases excited by individual earthquakes, or through analysis of the elastic gravitational free oscillations that individual earthquakes of sufficiently large magnitude may excite, has been the central focus of Earth physics research for more than a century. Unfortunately, data provide no information that is directly relevant to understanding the solid state 'flow' of the polycrystalline outer 'mantle' shell of the planet that is involved in the thermally driven convective circulation that is responsible for powering the 'drift' of the continents and which controls the rate of planetary cooling on long timescales. For this reason, there has been an increasing focus on the understanding of physical phenomenology that is unambiguously associated with mantle flow processes that are distinct from those directly associated with the convective circulation itself. This paper reviews the past many decades of work that has been invested in understanding the most important of such processes, namely that which has come to be referred to as 'glacial isostatic adjustment' (GIA). This process concerns the response of the planet to the loading and unloading of the high latitude continents by the massive accumulations of glacial ice that have occurred with almost metronomic regularity over the most recent million years of Earth history. Forced by the impact of gravitationaln-body effects on the geometry of Earth's orbit around the Sun through the impact upon the terrestrial regime of received solar insolation, these surface mass loads on the continents have left indelible records of their occurrence in the 'Earth system' consisting of the oceans, continents, and the great polar ice sheets on Greenland and Antarctica themselves. Although this ice-age phenomenology has been clearly recognized since early in the last century, it was for over 50 years considered to be no more than an interesting curiosity, the understanding of which remained on the periphery of the theoretical physics of the Earth. This was the case in part because no globally applicable theory was available that could be applied to rigorously interpret the observations. Equally important to understanding the scientific lethargy that held back the understanding of this phenomenon involving mantle flow processes was the lack of appreciation of the wide range of observations that were in fact related to GIA physics. This paper is devoted to a review of the global theories of the GIA process that have since been developed as a means of interpreting the extensive variety of observations that are now recognized as being involved in the response of the planet to the loading and unloading of its surface by glacial ice. The paper will also provide examples of the further analyses of Earth physics and climate related processes that applications of the modern theoretical structures have enabled.
Collapse
Affiliation(s)
| | | | - Donald F Argus
- Jet Propulsion Laboratory, Caltech, United States of America
| | - Tanghua Li
- Earth Observatory of Singapore, Nanyang Technological University, Singapore
| | | |
Collapse
|
19
|
Kendall E, Faccenda M, Ferreira AMG, Chang S. On the Relationship Between Oceanic Plate Speed, Tectonic Stress, and Seismic Anisotropy. GEOPHYSICAL RESEARCH LETTERS 2022; 49:e2022GL097795. [PMID: 36247518 PMCID: PMC9539886 DOI: 10.1029/2022gl097795] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 01/20/2022] [Revised: 06/17/2022] [Accepted: 06/25/2022] [Indexed: 06/16/2023]
Abstract
Seismic radial anisotropy (the squared ratio between the speeds of horizontally and vertically polarized shear waves,ξ = V S H 2 V S V 2 ) is a powerful tool to probe the direction of mantle flow and accumulated strain. While previous studies have confirmed the dependence of azimuthal anisotropy on plate speed, the first order control on radial anisotropy is unclear. In this study, we develop 2D ridge flow models combined with mantle fabric calculations to report that faster plates generate higher tectonics stresses and strain rates which lower the dislocation creep viscosity and lead to deeper anisotropy than beneath slower plates. We apply the SGLOBE-rani tomographic filter, resulting in a flat depth-age trend and stronger anisotropy beneath faster plates, which correlates well with 3D global anisotropic mantle models. Our predictions and observations suggest that as plate speed increases from 2 to 8 cm/yr, radial anisotropy increases by ∼0.01-0.025 in the upper 100-200 km of the mantle between 10 and 60 Ma.
Collapse
Affiliation(s)
- E. Kendall
- Department of Earth SciencesUniversity College LondonLondonUK
- GFZ German Research Centre for GeosciencesPotsdamGermany
| | - M. Faccenda
- Dipartimento di GeoscienzeUniversità di PadovaPaduaItaly
| | - A. M. G. Ferreira
- Department of Earth SciencesUniversity College LondonLondonUK
- CERISInstituto Superior TécnicoUniversidade de LisboaLisbonPortugal
| | - S.‐J. Chang
- Department of GeophysicsKangwon National UniversityChuncheonSouth Korea
| |
Collapse
|
20
|
Coupled influence of tectonics, climate, and surface processes on landscape evolution in southwestern North America. Nat Commun 2022; 13:4437. [PMID: 35915061 PMCID: PMC9343645 DOI: 10.1038/s41467-022-31903-2] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/22/2021] [Accepted: 07/07/2022] [Indexed: 11/08/2022] Open
Abstract
The Cenozoic landscape evolution in southwestern North America is ascribed to crustal isostasy, dynamic topography, or lithosphere tectonics, but their relative contributions remain controversial. Here we reconstruct landscape history since the late Eocene by investigating the interplay between mantle convection, lithosphere dynamics, climate, and surface processes using fully coupled four-dimensional numerical models. Our quantified depth-dependent strain rate and stress history within the lithosphere, under the influence of gravitational collapse and sub-lithospheric mantle flow, show that high gravitational potential energy of a mountain chain relative to a lower Colorado Plateau can explain extension directions and stress magnitudes in the belt of metamorphic core complexes during topographic collapse. Profound lithospheric weakening through heating and partial melting, following slab rollback, promoted this extensional collapse. Landscape evolution guided northeast drainage onto the Colorado Plateau during the late Eocene-late Oligocene, south-southwest drainage reversal during the late Oligocene-middle Miocene, and southwest drainage following the late Miocene.
Collapse
|
21
|
Vilacís B, Hayek JN, Stotz IL, Bunge HP, Friedrich AM, Carena S, Clark S. Evidence for active upper mantle flow in the Atlantic and Indo-Australian realms since the Upper Jurassic from hiatus maps and spreading rate changes. Proc Math Phys Eng Sci 2022; 478:20210764. [PMID: 35756875 PMCID: PMC9199074 DOI: 10.1098/rspa.2021.0764] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/01/2021] [Accepted: 05/12/2022] [Indexed: 11/12/2022] Open
Abstract
Histories of large-scale horizontal and vertical lithosphere motion hold important information on mantle convection. Here, we compare continent-scale hiatus maps as a proxy for mantle flow induced dynamic topography and plate motion variations in the Atlantic and Indo-Australian realms since the Upper Jurassic, finding they frequently correlate, except when plate boundary forces may play a significant role. This correlation agrees with descriptions of asthenosphere flow beneath tectonic plates in terms of Poiseuille/Couette flow, as it explicitly relates plate motion changes, induced by evolving basal shear forces, to non-isostatic vertical motion of the lithosphere. Our analysis reveals a timescale, on the order of a geological series, between the occurrence of continent-scale hiatus and plate motion changes. This is consistent with the presence of a weak upper mantle. It also shows a spatial scale for interregional hiatus, on the order of 2000-3000 km in diameter, which can be linked by fluid dynamic analysis to active upper mantle flow regions. Our results suggest future studies should pursue large-scale horizontal and vertical lithosphere motion in combination, to track the expressions of past mantle flow. Such studies would provide powerful constraints for adjoint-based geodynamic inverse models of past mantle convection.
Collapse
Affiliation(s)
- Berta Vilacís
- Department of Earth and Environmental Sciences, Ludwig-Maximilians-Universität München, Theresienstraße 41 and Luisenstraße 37, Munich 80333 Germany
| | - Jorge N Hayek
- Department of Earth and Environmental Sciences, Ludwig-Maximilians-Universität München, Theresienstraße 41 and Luisenstraße 37, Munich 80333 Germany
| | - Ingo L Stotz
- Department of Earth and Environmental Sciences, Ludwig-Maximilians-Universität München, Theresienstraße 41 and Luisenstraße 37, Munich 80333 Germany
| | - Hans-Peter Bunge
- Department of Earth and Environmental Sciences, Ludwig-Maximilians-Universität München, Theresienstraße 41 and Luisenstraße 37, Munich 80333 Germany
| | - Anke M Friedrich
- Department of Earth and Environmental Sciences, Ludwig-Maximilians-Universität München, Theresienstraße 41 and Luisenstraße 37, Munich 80333 Germany
| | - Sara Carena
- Department of Earth and Environmental Sciences, Ludwig-Maximilians-Universität München, Theresienstraße 41 and Luisenstraße 37, Munich 80333 Germany
| | - Stuart Clark
- University of New South Wales Sydney, Minerals and Energy Res. Eng., Kensington, New South Wales 2052, Australia
| |
Collapse
|
22
|
Liu CZ, Dick HJ, Mitchell RN, Wei W, Zhang ZY, Hofmann AW, Yang JF, Li Y. Archean cratonic mantle recycled at a mid-ocean ridge. SCIENCE ADVANCES 2022; 8:eabn6749. [PMID: 35648865 PMCID: PMC9159695 DOI: 10.1126/sciadv.abn6749] [Citation(s) in RCA: 8] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 12/12/2021] [Accepted: 04/15/2022] [Indexed: 05/26/2023]
Abstract
Basalts and mantle peridotites of mid-ocean ridges are thought to sample Earth's upper mantle. Osmium isotopes of abyssal peridotites uniquely preserve melt extraction events throughout Earth history, but existing records only indicate ages up to ~2 billion years (Ga) ago. Thus, the memory of the suspected large volumes of mantle lithosphere that existed in Archean time (>2.5 Ga) has apparently been lost somehow. We report abyssal peridotites with melt-depletion ages up to 2.8 Ga, documented by extremely unradiogenic 187Os/188Os ratios (to as low as 0.1095) and refractory major elements that compositionally resemble the deep keels of Archean cratons. These oceanic rocks were thus derived from the once-extensive Archean continental keels that have been dislodged and recycled back into the mantle, the feasibility of which we confirm with numerical modeling. This unexpected connection between young oceanic and ancient continental lithosphere indicates an underappreciated degree of compositional recycling over time.
Collapse
Affiliation(s)
- Chuan-Zhou Liu
- State Key Laboratory of Lithospheric Evolution, Institute of Geology and Geophysics, Chinese Academy of Sciences, Beijing 100029, China
- Laboratory for Marine Geology, Qingdao National Laboratory for Marine Science and Technology, Qingdao 266061, China
- CAS Center for Excellence in Deep Earth Sciences, Guangzhou 510640, China
- University of Chinese Academy of Sciences, Beijing 100049, China
| | - Henry J.B. Dick
- Department of Geology and Geophysics, Woods Hole Oceanographic Institution, Woods Hole, MA 02543, USA
| | - Ross N. Mitchell
- State Key Laboratory of Lithospheric Evolution, Institute of Geology and Geophysics, Chinese Academy of Sciences, Beijing 100029, China
| | - Wu Wei
- State Key Laboratory of Lithospheric Evolution, Institute of Geology and Geophysics, Chinese Academy of Sciences, Beijing 100029, China
- University of Chinese Academy of Sciences, Beijing 100049, China
| | - Zhen-Yu Zhang
- State Key Laboratory of Lithospheric Evolution, Institute of Geology and Geophysics, Chinese Academy of Sciences, Beijing 100029, China
- University of Chinese Academy of Sciences, Beijing 100049, China
| | | | - Jian-Feng Yang
- State Key Laboratory of Lithospheric Evolution, Institute of Geology and Geophysics, Chinese Academy of Sciences, Beijing 100029, China
| | - Yang Li
- State Key Laboratory of Lithospheric Evolution, Institute of Geology and Geophysics, Chinese Academy of Sciences, Beijing 100029, China
| |
Collapse
|
23
|
Wolf SG, Huismans RS, Braun J, Yuan X. Topography of mountain belts controlled by rheology and surface processes. Nature 2022; 606:516-521. [PMID: 35650431 DOI: 10.1038/s41586-022-04700-6] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/04/2020] [Accepted: 03/29/2022] [Indexed: 11/09/2022]
Abstract
It is widely recognized that collisional mountain belt topography is generated by crustal thickening and lowered by river bedrock erosion, linking climate and tectonics1-4. However, whether surface processes or lithospheric strength control mountain belt height, shape and longevity remains uncertain. Additionally, how to reconcile high erosion rates in some active orogens with long-term survival of mountain belts for hundreds of millions of years remains enigmatic. Here we investigate mountain belt growth and decay using a new coupled surface process5,6 and mantle-scale tectonic model7. End-member models and the new non-dimensional Beaumont number, Bm, quantify how surface processes and tectonics control the topographic evolution of mountain belts, and enable the definition of three end-member types of growing orogens: type 1, non-steady state, strength controlled (Bm > 0.5); type 2, flux steady state8, strength controlled (Bm ≈ 0.4-0.5); and type 3, flux steady state, erosion controlled (Bm < 0.4). Our results indicate that tectonics dominate in Himalaya-Tibet and the Central Andes (both type 1), efficient surface processes balance high convergence rates in Taiwan (probably type 2) and surface processes dominate in the Southern Alps of New Zealand (type 3). Orogenic decay is determined by erosional efficiency and can be subdivided into two phases with variable isostatic rebound characteristics and associated timescales. The results presented here provide a unified framework explaining how surface processes and lithospheric strength control the height, shape, and longevity of mountain belts.
Collapse
Affiliation(s)
- Sebastian G Wolf
- Department of Earth Science, University of Bergen, Bergen, Norway.
| | | | - Jean Braun
- Helmholtz Centre Potsdam, German Research Centre for Geosciences, Potsdam, Germany.,Institute of Geosciences, University of Potsdam, Potsdam, Germany
| | - Xiaoping Yuan
- Helmholtz Centre Potsdam, German Research Centre for Geosciences, Potsdam, Germany.,Hubei Key Laboratory of Critical Zone Evolution, School of Earth Sciences, China University of Geosciences, Wuhan, China
| |
Collapse
|
24
|
Theunissen T, Huismans RS. Mantle exhumation at magma-poor rifted margins controlled by frictional shear zones. Nat Commun 2022; 13:1634. [PMID: 35347122 PMCID: PMC8960832 DOI: 10.1038/s41467-022-29058-1] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/31/2021] [Accepted: 02/17/2022] [Indexed: 12/02/2022] Open
Abstract
The transition zone from continental crust to the mature mid-ocean ridge spreading center of the Iberia-Newfoundland magma-poor rifted margins is mostly composed of exhumed mantle characterized by highs and domes with varying elevation, spacing and shape. The mechanism controlling strain localization and fault migration explaining the geometry of these peridotite ridges is poorly understood. Here we show using forward geodynamic models that multiple out-of-sequence detachments with recurring dip reversal form during magma-poor rifting and mantle exhumation as a consequence of the strength competition between weak frictional-plastic shear zones and the thermally weakened necking domain beneath the exhuming footwall explaining geometry of these peridotite ridges. Model behaviour also shows that fault types and detachment styles vary with spreading rate and fault strength and confirm that these results can be compared to other magma poor passive margins such as along Antarctica-Australia and to ultra-slow mid-ocean spreading systems as the South-West Indian Ridge.
Collapse
Affiliation(s)
- Thomas Theunissen
- Department of Earth Science, University of Bergen, Allegaten 41, Postboks 7803, N-5020, Bergen, Norway.
| | - Ritske S Huismans
- Department of Earth Science, University of Bergen, Allegaten 41, Postboks 7803, N-5020, Bergen, Norway
| |
Collapse
|
25
|
Michaut C, Neufeld JA. Formation of the Lunar Primary Crust From a Long-Lived Slushy Magma Ocean. GEOPHYSICAL RESEARCH LETTERS 2022; 49:e2021GL095408. [PMID: 35865331 PMCID: PMC9286579 DOI: 10.1029/2021gl095408] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 09/06/2021] [Revised: 12/07/2021] [Accepted: 12/10/2021] [Indexed: 06/15/2023]
Abstract
Classical fractional crystallization scenarios of early lunar evolution suggest crustal formation by the flotation of light anorthite minerals from a liquid magma ocean. However, this model is challenged by the> 200 Myr age range of primitive ferroan anorthosites, their concordance with Mg-suite magmatism and by the compositional diversity observed in lunar anorthosites. Here, we propose a new model of slushy magma ocean crystallization in which crystals remain suspended in the lunar interior and crust formation only begins once a critical crystal content is reached. Thereafter crustal formation occurs by buoyant melt extraction and magmatism. The mixture viscosity strongly depends on temperature and solid fraction driving the development of a surface stagnant lid where enhanced solidification and buoyant ascent of melt lead to an anorthite-enriched crust. This model explains lunar anorthosites heterogeneity and suggests a crustal formation timescale of 100s Ma, reconciling anorthosite ages with an early age of the Moon.
Collapse
Affiliation(s)
- Chloé Michaut
- Ecole Normale Supérieure de LyonUniversité de LyonUniversité Claude Bernard Lyon 1Laboratoire de Géologie de Lyon, Terre, Planètes, EnvironnementLyonFrance
- Institut Universitaire de FranceParisFrance
| | - Jerome A. Neufeld
- Centre for Environmental and Industrial FlowsUniversity of CambridgeCambridgeUK
- Department of Earth SciencesUniversity of CambridgeCambridgeUK
- Department of Applied Mathematics and Theoretical PhysicsUniversity of CambridgeCambridgeUK
| |
Collapse
|
26
|
Gerya TV, Bercovici D, Becker TW. Dynamic slab segmentation due to brittle-ductile damage in the outer rise. Nature 2021; 599:245-250. [PMID: 34759365 DOI: 10.1038/s41586-021-03937-x] [Citation(s) in RCA: 12] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/30/2020] [Accepted: 08/20/2021] [Indexed: 11/09/2022]
Abstract
Subduction is the major plate driving force, and the strength of the subducting plate controls many aspects of the thermochemical evolution of Earth. Each subducting plate experiences intense normal faulting1-9 during bending that accommodates the transition from horizontal to downwards motion at the outer rise at trenches. Here we investigate the consequences of this bending-induced plate damage using numerical subduction models in which both brittle and ductile deformation, including grain damage, are tracked and coupled self-consistently. Pervasive slab weakening and pronounced segmentation can occur at the outer-rise region owing to the strong feedback between brittle and ductile damage localization. This slab-damage phenomenon explains the subduction dichotomy of strong plates and weak slabs10, the development of large-offset normal faults6,7 near trenches, the occurrence of segmented seismic velocity anomalies11 and distinct interfaces imaged within subducted slabs12,13, and the appearance of deep, localized intraplate areas of reduced effective viscosity14 observed at trenches. Furthermore, brittle-viscously damaged slabs show a tendency for detachment at elevated mantle temperatures. Given Earth's planetary cooling history15, this implies that intermittent subduction with frequent slab break-off episodes16 may have been characteristic for Earth until more recent times than previously suggested17.
Collapse
Affiliation(s)
- T V Gerya
- Department of Earth Sciences, Swiss Federal Institute of Technology Zurich, Zurich, Switzerland.
| | - D Bercovici
- Department of Earth and Planetary Sciences, Yale University, New Haven, CT, USA
| | - T W Becker
- Institute for Geophysics, Jackson School of Geosciences, The University of Texas at Austin, Austin, TX, USA.,Department of Geological Sciences, Jackson School of Geosciences, The University of Texas at Austin, Austin, TX, USA.,Oden Institute for Computational Engineering & Sciences, The University of Texas at Austin, Austin, TX, USA
| |
Collapse
|
27
|
Was There Land on the Early Earth? Life (Basel) 2021; 11:life11111142. [PMID: 34833018 PMCID: PMC8623345 DOI: 10.3390/life11111142] [Citation(s) in RCA: 12] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/29/2021] [Revised: 10/22/2021] [Accepted: 10/25/2021] [Indexed: 11/17/2022] Open
Abstract
The presence of exposed land on the early Earth is a prerequisite for a certain type of prebiotic chemical evolution in which the oscillating activity of water, driven by short-term, day–night, and seasonal cycles, facilitates the synthesis of proto-biopolymers. Exposed land is, however, not guaranteed to exist on the early Earth, which is likely to have been drastically different from the modern Earth. This mini-review attempts to provide an up-to-date account on the possibility of exposed land on the early Earth by integrating recent geological and geophysical findings. Owing to the competing effects of the growing ocean and continents in the Hadean, a substantial expanse of the Earth’s surface (∼20% or more) could have been covered by exposed continents in the mid-Hadean. In contrast, exposed land may have been limited to isolated ocean islands in the late Hadean and early Archean. The importance of exposed land during the origins of life remains an open question.
Collapse
|
28
|
Lo Bue R, Faccenda M, Yang J. The Role of Adria Plate Lithospheric Structures on the Recent Dynamics of the Central Mediterranean Region. JOURNAL OF GEOPHYSICAL RESEARCH. SOLID EARTH 2021; 126:e2021JB022377. [PMID: 35845546 PMCID: PMC9285053 DOI: 10.1029/2021jb022377] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 05/09/2021] [Revised: 08/25/2021] [Accepted: 09/20/2021] [Indexed: 06/15/2023]
Abstract
The Tertiary tectonic evolution of the Central Mediterranean has been relatively well constrained by abundant geological data. Yet, several uncertainties persist about the mechanisms that led to the present-day surface morphology and deep slab geometry. Here, we combine geodynamic and seismological numerical modeling techniques to reproduce the recent large-scale evolution of the Central Mediterranean and the associated strain-induced upper mantle fabrics and seismic anisotropy. 3D thermo-mechanical subduction models were designed and calibrated according to paleogeographic-tectonic reconstructions and seismological observations available in the literature. It is found that, although the opening of back-arc extensional basins in response to the retreat of the Ionian slab is a common feature in all models, structural heterogeneities within the Adria plate and/or the geometry of its Tyrrhenian passive margin profoundly impact on the segmentation of the subducting slab and the amount of Ionian trench retreat. More, in general, this study highlights the importance of coupling geodynamic and seismological modeling to better constrain the tectonic evolution of complex convergent margins such as the Central Mediterranean.
Collapse
Affiliation(s)
- Rosalia Lo Bue
- Dipartimento di GeoscienzeUniversità di PadovaPadovaItaly
| | | | - Jianfeng Yang
- Dipartimento di GeoscienzeUniversità di PadovaPadovaItaly
| |
Collapse
|
29
|
Differentiating induced versus spontaneous subduction initiation using thermomechanical models and metamorphic soles. Nat Commun 2021; 12:4632. [PMID: 34330929 PMCID: PMC8324860 DOI: 10.1038/s41467-021-24896-x] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/30/2020] [Accepted: 07/08/2021] [Indexed: 02/07/2023] Open
Abstract
Despite the critical role of subduction in plate tectonics, the dynamics of its initiation remains unclear. High-temperature low-pressure metamorphic soles are vestiges of subduction initiation, providing records of the pressure and temperature conditions along the subducting slab surface during subduction initiation that can possibly differentiate the two end-member subduction initiation modes: spontaneous and induced. Here, using numerical models, we show that the slab surface temperature reaches 800-900 °C at ~1 GPa over a wide range of parameter values for spontaneous subduction initiation whereas for induced subduction initiation, such conditions can be reached only if the age of the overriding plate is <5 Ma. These modeling results indicate that spontaneous subduction initiation would be more favorable for creating high-temperature conditions. However, the synthesis of our modeling results and geological observations indicate that the majority of the metamorphic soles likely formed during induced subduction initiation that involved a young overriding plate.
Collapse
|
30
|
Kobayashi Y, Gomyo H, Arai N. Molecular Insight into the Possible Mechanism of Drag Reduction of Surfactant Aqueous Solution in Pipe Flow. Int J Mol Sci 2021; 22:ijms22147573. [PMID: 34299196 PMCID: PMC8307477 DOI: 10.3390/ijms22147573] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/07/2021] [Revised: 07/03/2021] [Accepted: 07/12/2021] [Indexed: 11/16/2022] Open
Abstract
The phenomenon of drag reduction (known as the “Toms effect”) has many industrial and engineering applications, but a definitive molecular-level theory has not yet been constructed. This is due both to the multiscale nature of complex fluids and to the difficulty of directly observing self-assembled structures in nonequilibrium states. On the basis of a large-scale coarse-grained molecular simulation that we conducted, we propose a possible mechanism of turbulence suppression in surfactant aqueous solution. We demonstrate that maintaining sufficiently large micellar structures and a homogeneous radial distribution of surfactant molecules is necessary to obtain the drag-reduction effect. This is the first molecular-simulation evidence that a micellar structure is responsible for drag reduction in pipe flow, and should help in understanding the mechanisms underlying drag reduction by surfactant molecules under nonequilibrium conditions.
Collapse
|
31
|
Lee C, Kim Y. Role of warm subduction in the seismological properties of the forearc mantle: An example from southwest Japan. SCIENCE ADVANCES 2021; 7:7/28/eabf8934. [PMID: 34244142 PMCID: PMC8270487 DOI: 10.1126/sciadv.abf8934] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 11/26/2020] [Accepted: 05/27/2021] [Indexed: 06/13/2023]
Abstract
A warm slab thermal structure plays an important role in controlling seismic properties of the slab and mantle wedge. Among warm subduction zones, most notably in southwest Japan, the spatial distribution of large S-wave delay times and deep nonvolcanic tremors in the forearc mantle indicate the presence of a serpentinite layer along the slab interface. However, the conditions under which such a layer is generated remains unclear. Using numerical models, we here show that a serpentinite layer begins to develop by the slab-derived fluids below the deeper end of the slab-mantle decoupling interface and grows toward the corner of the mantle wedge along the interface under warm subduction conditions only, explaining the large S-wave delay times in the forearc mantle. The serpentinite layer then allows continuous free-fluid flow toward the corner of the mantle wedge, presenting possible mechanisms for the deep nonvolcanic tremors in the forearc mantle.
Collapse
Affiliation(s)
- Changyeol Lee
- Department of Earth System Sciences, Yonsei University, 50 Yonsei-ro Seodaemun-gu, Seoul 03722, Republic of Korea
| | - YoungHee Kim
- School of Earth and Environmental Sciences, Seoul National University, 1 Gwanak-ro, Gwanak-gu, Seoul 08826, Republic of Korea.
| |
Collapse
|
32
|
Lu G, Huismans RS. Melt volume at Atlantic volcanic rifted margins controlled by depth-dependent extension and mantle temperature. Nat Commun 2021; 12:3894. [PMID: 34162843 PMCID: PMC8222230 DOI: 10.1038/s41467-021-23981-5] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/29/2020] [Accepted: 05/21/2021] [Indexed: 02/05/2023] Open
Abstract
Breakup volcanism along rifted passive margins is highly variable in time and space. The factors controlling magmatic activity during continental rifting and breakup are not resolved and controversial. Here we use numerical models to investigate melt generation at rifted margins with contrasting rifting styles corresponding to those observed in natural systems. Our results demonstrate a surprising correlation of enhanced magmatism with margin width. This relationship is explained by depth-dependent extension, during which the lithospheric mantle ruptures earlier than the crust, and is confirmed by a semi-analytical prediction of melt volume over margin width. The results presented here show that the effect of increased mantle temperature at wide volcanic margins is likely over-estimated, and demonstrate that the large volumes of magmatism at volcanic rifted margin can be explained by depth-dependent extension and very moderate excess mantle potential temperature in the order of 50-80 °C, significantly smaller than previously suggested.
Collapse
Affiliation(s)
- Gang Lu
- grid.7914.b0000 0004 1936 7443Department of Earth Science, Bergen University, Bergen, Norway
| | - Ritske S. Huismans
- grid.7914.b0000 0004 1936 7443Department of Earth Science, Bergen University, Bergen, Norway
| |
Collapse
|
33
|
Ojha L, Karunatillake S, Karimi S, Buffo J. Amagmatic hydrothermal systems on Mars from radiogenic heat. Nat Commun 2021; 12:1754. [PMID: 33741920 PMCID: PMC7979869 DOI: 10.1038/s41467-021-21762-8] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/22/2020] [Accepted: 02/04/2021] [Indexed: 11/30/2022] Open
Abstract
Long-lived hydrothermal systems are prime targets for astrobiological exploration on Mars. Unlike magmatic or impact settings, radiogenic hydrothermal systems can survive for >100 million years because of the Ga half-lives of key radioactive elements (e.g., U, Th, and K), but remain unknown on Mars. Here, we use geochemistry, gravity, topography data, and numerical models to find potential radiogenic hydrothermal systems on Mars. We show that the Eridania region, which once contained a vast inland sea, possibly exceeding the combined volume of all other Martian surface water, could have readily hosted a radiogenic hydrothermal system. Thus, radiogenic hydrothermalism in Eridania could have sustained clement conditions for life far longer than most other habitable sites on Mars. Water radiolysis by radiogenic heat could have produced H2, a key electron donor for microbial life. Furthermore, hydrothermal circulation may help explain the region's high crustal magnetic field and gravity anomaly.
Collapse
Affiliation(s)
- Lujendra Ojha
- Department of Earth and Planetary Sciences. Rutgers, The State University of New Jersey, Piscataway, NJ, USA.
| | - Suniti Karunatillake
- Department of Geology and Geophysics, Louisiana State University, Baton Rouge, LA, USA
| | - Saman Karimi
- Department of Earth and Planetary Sciences, Johns Hopkins University, Baltimore, MD, USA
| | - Jacob Buffo
- Thayer School of Engineering, Dartmouth College, Hanover, NH, USA
| |
Collapse
|
34
|
Thermochemical lithosphere differentiation and the origin of cratonic mantle. Nature 2020; 588:89-94. [PMID: 33268867 DOI: 10.1038/s41586-020-2976-3] [Citation(s) in RCA: 18] [Impact Index Per Article: 3.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/25/2020] [Accepted: 09/17/2020] [Indexed: 11/09/2022]
Abstract
Cratons record the early history of continental lithosphere formation, yet how they became the most enduring part of the lithosphere on Earth remains unknown1. Here we propose a mechanism for the formation of large volumes of melt-depleted cratonic lithospheric mantle (CLM) and its evolution to stable cratons. Numerical models show large decompression melting of a hot, early Earth mantle beneath a stretching lithosphere, where melt extraction leaves large volumes of depleted mantle at depth. The dehydrated, stiffer mantle resists further deformation, forcing strain migration and cooling, thereby assimilating depleted mantle into the lithosphere. The negative feedback between strain localization and stiffening sustains long-term diffused extension and emplacement of large amounts of depleted CLM. The formation of CLM at low pressure and its deeper re-equilibration reproduces the evolution of Archaean lithosphere constrained by depth-temperature conditions1,2, whereas large degrees of depletion3,4 and melt volumes5 in Archaean cratons are best matched by models with lower lithospheric strength. Under these conditions, which are otherwise viable for plate tectonics6,7, thermochemical differentiation effectively prevents yielding and formation of margins: rifting and lithosphere subduction are short lived and embedded in the cooling CLM as relict structures, reproducing the recycling and reworking environments that are found in Archaean cratons8,9. Although they undergo major melting and extensive recycling during an early stage lasting approximately 500 million years, the modelled lithospheres progressively differentiate and stabilize, and then recycling and reworking become episodic. Early major melting and recycling events explain the production and loss of primordial Hadean lithosphere and crust10, whereas later stabilization and episodic reworking provides a context for the creation of continental cratons in the Archaean era4,8.
Collapse
|
35
|
Glerum A, Brune S, Stamps DS, Strecker MR. Victoria continental microplate dynamics controlled by the lithospheric strength distribution of the East African Rift. Nat Commun 2020; 11:2881. [PMID: 32513970 PMCID: PMC7280495 DOI: 10.1038/s41467-020-16176-x] [Citation(s) in RCA: 18] [Impact Index Per Article: 3.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/31/2019] [Accepted: 04/09/2020] [Indexed: 11/24/2022] Open
Abstract
The Victoria microplate between the Eastern and Western Branches of the East African Rift System is one of the largest continental microplates on Earth. In striking contrast to its neighboring plates, Victoria rotates counterclockwise with respect to Nubia. The underlying cause of this distinctive rotation has remained elusive so far. Using 3D numerical models, we investigate the role of pre-existing lithospheric heterogeneities in continental microplate rotation. We find that Victoria’s rotation is primarily controlled by the distribution of rheologically stronger zones that transmit the drag of the major plates to the microplate and of the mechanically weaker mobile belts surrounding Victoria that facilitate rotation. Our models reproduce Victoria’s GPS-derived counterclockwise rotation as well as key complexities of the regional tectonic stress field. These results reconcile competing ideas on the opening of the rift system by highlighting differences in orientation of the far-field divergence, local extension, and the minimum horizontal stress. One of the largest continental microplates on Earth is situated in the center of the East African Rift System, and oddly, the Victoria microplate rotates counterclockwise with respect to the neighboring African tectonic plate. Here, the authors' modelling results suggest that Victoria microplate rotation is caused by edge-driven lithospheric processes related to the specific geometry of rheologically weak and strong regions.
Collapse
Affiliation(s)
- Anne Glerum
- Helmholtz Centre Potsdam - GFZ German Research Centre for Geosciences, Potsdam, Germany.
| | - Sascha Brune
- Helmholtz Centre Potsdam - GFZ German Research Centre for Geosciences, Potsdam, Germany.,University of Potsdam, Potsdam-Golm, Germany
| | | | | |
Collapse
|
36
|
Uncertainty Quantification in Planetary Thermal History Models: Implications for Hypotheses Discrimination and Habitability Modeling. ACTA ACUST UNITED AC 2020. [DOI: 10.3847/1538-4357/ab822b] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/18/2022]
|
37
|
Maunder B, Prytulak J, Goes S, Reagan M. Rapid subduction initiation and magmatism in the Western Pacific driven by internal vertical forces. Nat Commun 2020; 11:1874. [PMID: 32312969 PMCID: PMC7170853 DOI: 10.1038/s41467-020-15737-4] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/30/2019] [Accepted: 03/12/2020] [Indexed: 11/25/2022] Open
Abstract
Plate tectonics requires the formation of plate boundaries. Particularly important is the enigmatic initiation of subduction: the sliding of one plate below the other, and the primary driver of plate tectonics. A continuous, in situ record of subduction initiation was recovered by the International Ocean Discovery Program Expedition 352, which drilled a segment of the fore-arc of the Izu-Bonin-Mariana subduction system, revealing a distinct magmatic progression with a rapid timescale (approximately 1 million years). Here, using numerical models, we demonstrate that these observations cannot be produced by previously proposed horizontal external forcing. Instead a geodynamic evolution that is dominated by internal, vertical forces produces both the temporal and spatial distribution of magmatic products, and progresses to self-sustained subduction. Such a primarily internally driven initiation event is necessarily whole-plate scale and the rock sequence generated (also found along the Tethyan margin) may be considered as a smoking gun for this type of event. The magmatic progression produced during the initiation of the Izu-Bonin-Marianas subduction zone took place rapidly over 1 million years, but it has been unclear why. Here, using numerical models, the authors show that subduction initiation was dominated by vertical forces, internal to the system itself, progressing to self-sustained subduction.
Collapse
Affiliation(s)
- B Maunder
- Imperial College London, Department of Earth Science and Engineering, Royal School of Mines, Prince Consort Road, South Kensington, London, SW7 2BP, UK.
| | - J Prytulak
- Durham University, Department of Earth Sciences, Science Labs, Lower Mountjoy, South Road, Durham, DH1 3LE, UK
| | - S Goes
- Imperial College London, Department of Earth Science and Engineering, Royal School of Mines, Prince Consort Road, South Kensington, London, SW7 2BP, UK
| | - M Reagan
- University of Iowa, Department of Earth and Envitonmental Sciences, 115 Trowbridge Hall, Iowa City, IA, 52240, USA
| |
Collapse
|
38
|
Immoor J, Marquardt H, Miyagi L, Speziale S, Merkel S, Schwark I, Ehnes A, Liermann HP. An improved setup for radial diffraction experiments at high pressures and high temperatures in a resistive graphite-heated diamond anvil cell. THE REVIEW OF SCIENTIFIC INSTRUMENTS 2020; 91:045121. [PMID: 32357741 DOI: 10.1063/1.5143293] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 12/22/2019] [Accepted: 04/07/2020] [Indexed: 06/11/2023]
Abstract
We present an improved setup for the experimental study of deformation of solids at simultaneous high pressures and temperatures by radial x-ray diffraction. This technique employs a graphite resistive heated Mao-Bell type diamond anvil cell for radial x-ray diffraction in combination with a water-cooled vacuum chamber. The new chamber has been developed by the sample environment group at PETRA III and implemented at the Extreme Conditions Beamline P02.2 at PETRA III, DESY (Hamburg, Germany). We discuss applications of the new setup to study deformation of a variety of materials, including ferropericlase, calcium perovskite, bridgmanite, and tantalum carbide, at high-pressure/temperature.
Collapse
Affiliation(s)
- J Immoor
- Bayerisches Geoinstitut BGI, University of Bayreuth, 95440 Bayreuth, Germany
| | - H Marquardt
- Department of Earth Sciences, University of Oxford, Oxford OX1 3AN, United Kingdom
| | - L Miyagi
- University of Utah, 115 So. 1460 E., Salt Lake City, Utah UT84112-0111, USA
| | - S Speziale
- German Research Center for Geosciences GFZ, 14473 Potsdam, Germany
| | - S Merkel
- Univ. Lille, CNRS, INRAE, Centrale Lille, UMR 8207 - UMET - Unité Matériaux et Transformations, F-59000 Lille, France
| | - I Schwark
- Photon Sciences, Deutsches Elektronen-Synchrotron (DESY), 22607 Hamburg, Germany
| | - A Ehnes
- Photon Sciences, Deutsches Elektronen-Synchrotron (DESY), 22607 Hamburg, Germany
| | - H-P Liermann
- Photon Sciences, Deutsches Elektronen-Synchrotron (DESY), 22607 Hamburg, Germany
| |
Collapse
|
39
|
Intraplate volcanism originating from upwelling hydrous mantle transition zone. Nature 2020; 579:88-91. [PMID: 32103183 DOI: 10.1038/s41586-020-2045-y] [Citation(s) in RCA: 59] [Impact Index Per Article: 11.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/19/2019] [Accepted: 12/02/2019] [Indexed: 11/08/2022]
Abstract
Most magmatism occurring on Earth is conventionally attributed to passive mantle upwelling at mid-ocean ridges, to slab devolatilization at subduction zones, or to mantle plumes. However, the widespread Cenozoic intraplate volcanism in northeast China1-3 and the young petit-spot volcanoes4-7 offshore of the Japan Trench cannot readily be associated with any of these mechanisms. In addition, the mantle beneath these types of volcanism is characterized by zones of anomalously low seismic velocity above and below the transition zone8-12 (a mantle level located at depths between 410 and 660 kilometres). A comprehensive interpretation of these phenomena is lacking. Here we show that most (or possibly all) of the intraplate and petit-spot volcanism and low-velocity zones around the Japanese subduction zone can be explained by the Cenozoic interaction of the subducting Pacific slab with a hydrous mantle transition zone. Numerical modelling indicates that 0.2 to 0.3 weight per cent of water dissolved in mantle minerals that are driven out from the transition zone in response to subduction and retreat of a tectonic plate is sufficient to reproduce the observations. This suggests that a critical amount of water may have accumulated in the transition zone around this subduction zone, as well as in others of the Tethyan tectonic belt13 that are characterized by intraplate or petit-spot volcanism and low-velocity zones in the underlying mantle.
Collapse
|
40
|
Hertgen S, Yamato P, Guillaume B, Magni V, Schliffke N, van Hunen J. Influence of the Thickness of the Overriding Plate on Convergence Zone Dynamics. GEOCHEMISTRY, GEOPHYSICS, GEOSYSTEMS : G(3) 2020; 21:e2019GC008678. [PMID: 32714097 PMCID: PMC7375164 DOI: 10.1029/2019gc008678] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 09/20/2019] [Revised: 01/23/2020] [Accepted: 01/23/2020] [Indexed: 06/11/2023]
Abstract
The important role played by the upper plate in convergence zones dynamics has long been underestimated but is now more and more emphasized. However, the influence of its thickness and/or strength on orogenic systems evolution remains largely unknown. Here we present results from 3D thermo-mechanical numerical simulations of convergence zones (including oceanic subduction followed by continental subduction/collision), in which we vary the rheological profile of the overriding plate (OP). For this, we systematically modify the crustal thickness of the overriding lithosphere and the temperature at the Moho to obtain a thermal thickness of the overriding lithosphere ranging from 80 to 180 km. While all models share a common global evolution (i.e., slab sinking, interaction between slab and the 660 km discontinuity, continental subduction/collision, and slab breakoff), they also highlight first-order differences arising from the variations in the OP strength (thermal thickness). With a thin/weak OP, slab rollback is favored, the slab dip is low, the mantle flow above the slab is vigorous, and the trench migrates at a high rate compared to a thick/strong OP. In addition, slab breakoff and back-arc basin formation events occur significantly earlier than in models involving a thick OP. Our models therefore highlight the major role played by the thickness/strength of the OP on convergence zone dynamics and illustrate its influence in a quantitative way.
Collapse
Affiliation(s)
- Solenn Hertgen
- Univ Rennes, CNRS, Géosciences Rennes ‐ UMR 6118RennesFrance
| | - Philippe Yamato
- Univ Rennes, CNRS, Géosciences Rennes ‐ UMR 6118RennesFrance
- Institut Universitaire de France (IUF), ParisFrance
| | | | - Valentina Magni
- The Centre for Earth Evolution and Dynamics, Department of GeosciencesUniversity of OsloOsloNorway
| | | | | |
Collapse
|
41
|
Chen NC, Yang TF, Hong WL, Yu TL, Lin IT, Wang PL, Lin S, Su CC, Shen CC, Wang Y, Lin LH. Discharge of deeply rooted fluids from submarine mud volcanism in the Taiwan accretionary prism. Sci Rep 2020; 10:381. [PMID: 31942015 PMCID: PMC6962449 DOI: 10.1038/s41598-019-57250-9] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/24/2019] [Accepted: 12/24/2019] [Indexed: 11/14/2022] Open
Abstract
Qualitative and quantitative assessments of fluid cycling are essential to address the role and transport of deeply sourced fluids in subduction systems. In this study, sediment cores distributed across a submarine mud volcano (SMV) offshore southwestern Taiwan were investigated to determine the characteristics of fluids generated through the convergence between the Eurasian and Phillippine Sea Plates. The low dissolved chloride concentration combined with the enrichment of 18O, and depletion of 2H of pore fluids suggest the discharge of deep freshwater formed by smectite dehydration at an equilibrium temperature of 100 to 150 °C. The upward fluid velocities, decreasing from 2.0 to 5.0 cm yr−1 at the center to a negligible value at margin sites, varied with the rate and efficiency of anaerobic methanotrophy, demonstrating the impact of fluid migration on biogeochemical processes and carbon cycling. By extrapolating the velocity pattern, the flux of fluids exported from 13 SMVs into seawater amounted up to 1.3–2.5 × 107 kg yr−1, a quantity accounting for 1.1–28.6% of the smectite-bound water originally stored in the incoming sediments. Our results imply that SMVs could act as a conduit to channel the fluids produced from great depth/temperature into seafloor environments in a subduction system of the western Pacific Ocean.
Collapse
Affiliation(s)
- Nai-Chen Chen
- Department of Geoscience, National Taiwan University, Taipei, Taiwan
| | | | - Wei-Li Hong
- Geological Survey of Norway, Trondheim, Norway
| | - Tsai-Luen Yu
- Department of Geoscience, National Taiwan University, Taipei, Taiwan.,Research Center for Future Earth, National Taiwan University, Taipei, Taiwan
| | - In-Tian Lin
- Exploration and Development Research Institute, CPC, Taiwan
| | - Pei-Ling Wang
- Research Center for Future Earth, National Taiwan University, Taipei, Taiwan.,Institute of Oceanography, National Taiwan University, Taipei, Taiwan
| | - Saulwood Lin
- Institute of Oceanography, National Taiwan University, Taipei, Taiwan
| | - Chih-Chieh Su
- Research Center for Future Earth, National Taiwan University, Taipei, Taiwan.,Institute of Oceanography, National Taiwan University, Taipei, Taiwan
| | - Chuan-Chou Shen
- Department of Geoscience, National Taiwan University, Taipei, Taiwan.,Research Center for Future Earth, National Taiwan University, Taipei, Taiwan
| | | | - Li-Hung Lin
- Department of Geoscience, National Taiwan University, Taipei, Taiwan. .,Research Center for Future Earth, National Taiwan University, Taipei, Taiwan.
| |
Collapse
|
42
|
Whitehouse PL, Gomez N, King MA, Wiens DA. Solid Earth change and the evolution of the Antarctic Ice Sheet. Nat Commun 2019; 10:503. [PMID: 30700704 PMCID: PMC6353952 DOI: 10.1038/s41467-018-08068-y] [Citation(s) in RCA: 14] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/01/2017] [Accepted: 12/15/2018] [Indexed: 11/17/2022] Open
Abstract
Recent studies suggest that Antarctica has the potential to contribute up to ~15 m of sea-level rise over the next few centuries. The evolution of the Antarctic Ice Sheet is driven by a combination of climate forcing and non-climatic feedbacks. In this review we focus on feedbacks between the Antarctic Ice Sheet and the solid Earth, and the role of these feedbacks in shaping the response of the ice sheet to past and future climate changes. The growth and decay of the Antarctic Ice Sheet reshapes the solid Earth via isostasy and erosion. In turn, the shape of the bed exerts a fundamental control on ice dynamics as well as the position of the grounding line—the location where ice starts to float. A complicating issue is the fact that Antarctica is situated on a region of the Earth that displays large spatial variations in rheological properties. These properties affect the timescale and strength of feedbacks between ice-sheet change and solid Earth deformation, and hence must be accounted for when considering the future evolution of the ice sheet. The evolution of the Antarctic Ice Sheet is driven by a combination of climate forcing and non-climatic feedbacks. In this review, the authors focus on feedbacks between the Antarctic Ice Sheet and the solid Earth, and the role of these feedbacks in shaping the response of the ice sheet to past and future climate changes.
Collapse
Affiliation(s)
| | - Natalya Gomez
- Department of Earth and Planetary Sciences, McGill University, Montreal, H3A 0E8, Canada
| | - Matt A King
- School of Technology, Environments and Design, University of Tasmania, Hobart, TAS, 7001, Australia
| | - Douglas A Wiens
- Department of Earth and Planetary Sciences, Washington University, St Louis, MO, 63130, USA
| |
Collapse
|
43
|
LIP formation and protracted lower mantle upwelling induced by rifting and delamination. Sci Rep 2018; 8:16578. [PMID: 30410091 PMCID: PMC6224380 DOI: 10.1038/s41598-018-34194-0] [Citation(s) in RCA: 19] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/13/2018] [Accepted: 10/12/2018] [Indexed: 11/08/2022] Open
Abstract
Large Igneous Provinces (LIPs) are commonly attributed to mantle plumes, hot upwellings from the deep lower mantle, apparently unrelated to plate motions. However, LIPs often form in association with rifting and breakup. Using numerical modelling, we introduce a novel idea that explains plume-like mantle upwelling by plate tectonic processes. Our model indicates that rifting-induced delamination of orogenic lithosphere can perturb the thermochemical mantle stratification and induce lower mantle upwelling which causes syn-rift LIP formation followed by protracted and enhanced mid ocean ridge basalt (MORB) generation. Our model provides an explanation for the geographical correlation between the Caledonian suture, the North Atlantic Igneous Province (NAIP) and present-day Icelandic magmatism.
Collapse
|
44
|
Korenaga J. Crustal evolution and mantle dynamics through Earth history. PHILOSOPHICAL TRANSACTIONS. SERIES A, MATHEMATICAL, PHYSICAL, AND ENGINEERING SCIENCES 2018; 376:rsta.2017.0408. [PMID: 30275159 PMCID: PMC6189559 DOI: 10.1098/rsta.2017.0408] [Citation(s) in RCA: 20] [Impact Index Per Article: 2.9] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Accepted: 05/31/2018] [Indexed: 05/13/2023]
Abstract
Resolving the modes of mantle convection through Earth history, i.e. when plate tectonics started and what kind of mantle dynamics reigned before, is essential to the understanding of the evolution of the whole Earth system, because plate tectonics influences almost all aspects of modern geological processes. This is a challenging problem because plate tectonics continuously rejuvenates Earth's surface on a time scale of about 100 Myr, destroying evidence for its past operation. It thus becomes essential to exploit indirect evidence preserved in the buoyant continental crust, part of which has survived over billions of years. This contribution starts with an in-depth review of existing models for continental growth. Growth models proposed so far can be categorized into three types: crust-based, mantle-based and other less direct inferences, and the first two types are particularly important as their difference reflects the extent of crustal recycling, which can be related to subduction. Then, a theoretical basis for a change in the mode of mantle convection in the Precambrian is reviewed, along with a critical appraisal of some popular notions for early Earth dynamics. By combining available geological and geochemical observations with geodynamical considerations, a tentative hypothesis is presented for the evolution of mantle dynamics and its relation to surface environment; the early onset of plate tectonics and gradual mantle hydration are responsible not only for the formation of continental crust but also for its preservation as well as its emergence above sea level. Our current understanding of various material properties and elementary processes is still too premature to build a testable, quantitative model for this hypothesis, but such modelling efforts could potentially transform the nature of the data-starved early Earth research by quantifying the extent of preservation bias.This article is part of a discussion meeting issue 'Earth dynamics and the development of plate tectonics'.
Collapse
Affiliation(s)
- Jun Korenaga
- Department of Geology and Geophysics, Yale University, New Haven, CT 06520, USA
| |
Collapse
|
45
|
Nebel O, Capitanio FA, Moyen JF, Weinberg RF, Clos F, Nebel-Jacobsen YJ, Cawood PA. When crust comes of age: on the chemical evolution of Archaean, felsic continental crust by crustal drip tectonics. PHILOSOPHICAL TRANSACTIONS. SERIES A, MATHEMATICAL, PHYSICAL, AND ENGINEERING SCIENCES 2018; 376:rsta.2018.0103. [PMID: 30275165 PMCID: PMC6189554 DOI: 10.1098/rsta.2018.0103] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Accepted: 08/10/2018] [Indexed: 05/18/2023]
Abstract
The secular evolution of the Earth's crust is marked by a profound change in average crustal chemistry between 3.2 and 2.5 Ga. A key marker for this change is the transition from Archaean sodic granitoid intrusions of the tonalite-trondhjemite-granodiorite (TTG) series to potassic (K) granitic suites, akin (but not identical) to I-type granites that today are associated with subduction zones. It remains poorly constrained as to how and why this change was initiated and if it holds clues about the geodynamic transition from a pre-plate tectonic mode, often referred to as stagnant lid, to mobile plate tectonics. Here, we combine a series of proposed mechanisms for Archaean crustal geodynamics in a single model to explain the observed change in granitoid chemistry. Numeric modelling indicates that upper mantle convection drives crustal flow and subsidence, leading to profound diversity in lithospheric thickness with thin versus thick proto-plates. When convecting asthenospheric mantle interacts with lower lithosphere, scattered crustal drips are created. Under increasing P-T conditions, partial melting of hydrated meta-basalt within these drips produces felsic melts that intrude the overlying crust to form TTG. Dome structures, in which these melts can be preserved, are a positive diapiric expression of these negative drips. Transitional TTG with elevated K mark a second evolutionary stage, and are blends of subsided and remelted older TTG forming K-rich melts and new TTG melts. Ascending TTG-derived melts from asymmetric drips interact with the asthenospheric mantle to form hot, high-Mg sanukitoid. These melts are small in volume, predominantly underplated, and their heat triggered melting of lower crustal successions to form higher-K granites. Importantly, this evolution operates as a disseminated process in space and time over hundreds of millions of years (greater than 200 Ma) in all cratons. This focused ageing of the crust implies that compiled geochemical data can only broadly reflect geodynamic changes on a global or even craton-wide scale. The observed change in crustal chemistry does mark the lead up to but not the initiation of modern-style subduction.This article is part of a discussion meeting issue 'Earth dynamics and the development of plate tectonics'.
Collapse
Affiliation(s)
- O Nebel
- School of Earth, Atmosphere and Environment, Monash University, Clayton, 3800 Victoria, Australia
| | - F A Capitanio
- School of Earth, Atmosphere and Environment, Monash University, Clayton, 3800 Victoria, Australia
| | - J-F Moyen
- Laboratoire Magmas et Volcans, Université de Lyon, UJM-UCA-CNRS-IRD, 23 rue Dr. Paul Michelon, 42023 Saint Etienne, France
| | - R F Weinberg
- School of Earth, Atmosphere and Environment, Monash University, Clayton, 3800 Victoria, Australia
| | - F Clos
- School of Earth, Atmosphere and Environment, Monash University, Clayton, 3800 Victoria, Australia
| | | | - P A Cawood
- School of Earth, Atmosphere and Environment, Monash University, Clayton, 3800 Victoria, Australia
| |
Collapse
|
46
|
Foley BJ. The dependence of planetary tectonics on mantle thermal state: applications to early Earth evolution. PHILOSOPHICAL TRANSACTIONS. SERIES A, MATHEMATICAL, PHYSICAL, AND ENGINEERING SCIENCES 2018; 376:20170409. [PMID: 30275160 PMCID: PMC6189558 DOI: 10.1098/rsta.2017.0409] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.1] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Accepted: 07/24/2018] [Indexed: 05/23/2023]
Abstract
For plate tectonics to operate on a planet, mantle convective forces must be capable of forming weak, localized shear zones in the lithosphere that act as plate boundaries. Otherwise, a planet's mantle will convect in a stagnant lid regime, where subduction and plate motions are absent. Thus, when and how plate tectonics initiated on the Earth is intrinsically tied to the ability of mantle convection to form plate boundaries; however, the physics behind this process are still uncertain. Most mantle convection models have employed a simple pseudoplastic model of the lithosphere, where the lithosphere 'fails' and develops a mobile lid when stresses in the lithosphere reach the prescribed yield stress. With pseudoplasticity high mantle temperatures and high rates of internal heating, conditions relevant for the early Earth, impede plate boundary formation by decreasing lithospheric stresses, and hence favour a stagnant lid for the early Earth. However, when a model for shear zone formation based on grain size reduction is used, early Earth thermal conditions do not favour a stagnant lid. While lithosphere stress drops with increasing mantle temperature or heat production rate, the deformational work, which drives grain size reduction, increases. Thus, the ability of convection to form weak plate boundaries is not impeded by early Earth thermal conditions. However, mantle thermal state does change the style of subduction and lithosphere mobility; high mantle temperatures lead to a more sluggish, drip-like style of subduction. This 'sluggish lid' convection may be able to explain many of the key observations of early Earth crust formation processes preserved in the geologic record. Moreover, this work highlights the importance of understanding the microphysics of plate boundary formation for assessing early Earth tectonics, as different plate boundary formation mechanisms are influenced by mantle thermal state in fundamentally different ways.This article is part of a discussion meeting issue 'Earth dynamics and the development of plate tectonics'.
Collapse
Affiliation(s)
- Bradford J Foley
- Department of Geosciences, Pennsylvania State University, University Park, PA 16802, USA
| |
Collapse
|
47
|
Abstract
In the present-day Earth, some subducting plates (slabs) are flattening above the upper–lower mantle boundary at ~670 km depth, whereas others go through, indicating a mode between layered and whole-mantle convection. Previous models predicted that in a few hundred degree hotter early Earth, convection was likely more layered due to dominant slab stagnation. In self-consistent numerical models where slabs have a plate-like rheology, strong slabs and mobile plate boundaries favour stagnation for old and penetration for young slabs, as observed today. Here we show that such models predict slabs would have penetrated into the lower mantle more easily in a hotter Earth, when a weaker asthenosphere and decreased plate density and strength resulted in subduction almost without trench retreat. Thus, heat and material transport in the Earth’s mantle was more (rather than less) efficient in the past, which better matches the thermal evolution of the Earth. The subducting plates can either penetrate straight into the lower mantle or flatten in the mantle transition zone, yet slab dynamics in the past remains unclear. Here, using subduction models, the authors predict that a hotter early Earth was probably more favourable to lower mantle slab penetration.
Collapse
|
48
|
Foley BJ, Smye AJ. Carbon Cycling and Habitability of Earth-Sized Stagnant Lid Planets. ASTROBIOLOGY 2018; 18:873-896. [PMID: 30035642 DOI: 10.1089/ast.2017.1695] [Citation(s) in RCA: 12] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/08/2023]
Abstract
Models of thermal evolution, crustal production, and CO2 cycling are used to constrain the prospects for habitability of rocky planets, with Earth-like size and composition, in the stagnant lid regime. Specifically, we determine the conditions under which such planets can maintain rates of CO2 degassing large enough to prevent global surface glaciation but small enough so as not to exceed the upper limit on weathering rates provided by the supply of fresh rock, a situation which would lead to runaway atmospheric CO2 accumulation and an inhospitably hot climate. The models show that stagnant lid planets with initial radiogenic heating rates of 100-250 TW, and with total CO2 budgets ranging from ∼10-2 to 1 times Earth's estimated CO2 budget, can maintain volcanic outgassing rates suitable for habitability for ≈1-5 Gyr; larger CO2 budgets result in uninhabitably hot climates, while smaller budgets result in global glaciation. High radiogenic heat production rates favor habitability by sustaining volcanism and CO2 outgassing longer. Thus, the results suggest that plate tectonics may not be required for establishing a long-term carbon cycle and maintaining a stable, habitable climate. The model is necessarily highly simplified, as the uncertainties with exoplanet thermal evolution and outgassing are large. Nevertheless, the results provide some first-order guidance for future exoplanet missions, by predicting the age at which habitability becomes unlikely for a stagnant lid planet as a function of initial radiogenic heat budget. This prediction is powerful because both planet heat budget and age can potentially be constrained from stellar observations. Key Words: Exoplanets-Habitability-Stagnant lid tectonics-Carbon cycle-Volcanism. Astrobiology 18, 873-896.
Collapse
Affiliation(s)
- Bradford J Foley
- Department of Geosciences, Pennsylvania State University, University Park , Pennsylvania
| | - Andrew J Smye
- Department of Geosciences, Pennsylvania State University, University Park , Pennsylvania
| |
Collapse
|
49
|
Qiu Q, Moore JDP, Barbot S, Feng L, Hill EM. Transient rheology of the Sumatran mantle wedge revealed by a decade of great earthquakes. Nat Commun 2018. [PMID: 29520095 PMCID: PMC5843651 DOI: 10.1038/s41467-018-03298-6] [Citation(s) in RCA: 54] [Impact Index Per Article: 7.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/17/2022] Open
Abstract
Understanding the rheological properties of the upper mantle is essential to develop a consistent model of mantle dynamics and plate tectonics. However, the spatial distribution and temporal evolution of these properties remain unclear. Here, we infer the rheological properties of the asthenosphere across multiple great megathrust earthquakes between 2004 and 2014 along the Sumatran subduction zone, taking advantage of decade-long continuous GPS and tide-gauge measurements. We observe transient mantle wedge flow following these earthquakes, and infer the temporal evolution of the effective viscosity. We show that the evolution of stress and strain rate following these earthquakes is better matched by a bi-viscous than by a power-law rheology model, and we estimate laterally heterogeneous transient and background viscosities on the order of ~1017 and ~1019 Pa s, respectively. Our results constitute a preliminary rheological model to explain stress evolution within earthquake cycles and the development of seismic hazard in the region. The rheology of the upper mantle is key to understanding how plate tectonics may evolve. Here, using GPS and tide-gauge measurements along the Sumatran subduction zone, the authors’ show that a bi-viscous rheology model is needed to explain the stress and strain evolution of the upper mantle following earthquakes.
Collapse
Affiliation(s)
- Qiang Qiu
- Earth Observatory of Singapore, Nanyang Technological University, Singapore, 639798, Singapore. .,Asian School of the Environment, Nanyang Technological University, Singapore, 639798, Singapore.
| | - James D P Moore
- Earth Observatory of Singapore, Nanyang Technological University, Singapore, 639798, Singapore
| | - Sylvain Barbot
- Earth Observatory of Singapore, Nanyang Technological University, Singapore, 639798, Singapore. .,Asian School of the Environment, Nanyang Technological University, Singapore, 639798, Singapore.
| | - Lujia Feng
- Earth Observatory of Singapore, Nanyang Technological University, Singapore, 639798, Singapore
| | - Emma M Hill
- Earth Observatory of Singapore, Nanyang Technological University, Singapore, 639798, Singapore.,Asian School of the Environment, Nanyang Technological University, Singapore, 639798, Singapore
| |
Collapse
|
50
|
Nzogang BC, Bouquerel J, Cordier P, Mussi A, Girard J, Karato S. Characterization by Scanning Precession Electron Diffraction of an Aggregate of Bridgmanite and Ferropericlase Deformed at HP-HT. GEOCHEMISTRY, GEOPHYSICS, GEOSYSTEMS : G(3) 2018; 19:582-594. [PMID: 29937698 PMCID: PMC5993221 DOI: 10.1002/2017gc007244] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 09/18/2017] [Accepted: 01/27/2018] [Indexed: 06/08/2023]
Abstract
Scanning precession electron diffraction is an emerging promising technique for mapping phases and crystal orientations with short acquisition times (10-20 ms/pixel) in a transmission electron microscope similarly to the Electron Backscattered Diffraction (EBSD) or Transmission Kikuchi Diffraction (TKD) techniques in a scanning electron microscope. In this study, we apply this technique to the characterization of deformation microstructures in an aggregate of bridgmanite and ferropericlase deformed at 27 GPa and 2,130 K. Such a sample is challenging for microstructural characterization for two reasons: (i) the bridgmanite is very unstable under electron irradiation, (ii) under high stress conditions, the dislocation density is so large that standard characterization by diffraction contrast are limited, or impossible. Here we show that detailed analysis of intracrystalline misorientations sheds some light on the deformation mechanisms of both phases. In bridgmanite, deformation is accommodated by localized, amorphous, shear deformation lamellae whereas ferropericlase undergoes large strains leading to grain elongation in response to intense dislocation activity with no evidence for recrystallization. Plastic strain in ferropericlase can be semiquantitatively assessed by following kernel average misorientation distributions.
Collapse
Affiliation(s)
- B. C. Nzogang
- Univ. Lille, CNRS, INRA, ENSCL, UMR 8207‐UMET‐Unité Matériaux et TransformationsLilleFrance
| | - J. Bouquerel
- Univ. Lille, CNRS, INRA, ENSCL, UMR 8207‐UMET‐Unité Matériaux et TransformationsLilleFrance
| | - P. Cordier
- Univ. Lille, CNRS, INRA, ENSCL, UMR 8207‐UMET‐Unité Matériaux et TransformationsLilleFrance
| | - A. Mussi
- Univ. Lille, CNRS, INRA, ENSCL, UMR 8207‐UMET‐Unité Matériaux et TransformationsLilleFrance
| | - J. Girard
- Department of Geology and GeophysicsYale UniversityNew HavenConnecticutUSA
| | - S. Karato
- Department of Geology and GeophysicsYale UniversityNew HavenConnecticutUSA
| |
Collapse
|