Kinematic constraints on the Rodinia to Gondwana transition

Earth's plate tectonic history during the breakup of the supercontinent Pangea is well constrained from the seafloor spreading record, but evolving plate configurations during older supercontinent cycles are much less well understood. A relative paucity of available palaeomagnetic and geological data for deep time reconstructions necessitates innovative approaches to help discriminate between competing plate configurations. More difficult is tracing the journeys of individual continents during the amalgamation and breakup of supercontinents. Typically, deep-time reconstructions are built using absolute motions defined by palaeomagnetic data, and do not consider the kinematics of relative motions between plates, even for occasions where they are thought to be ‘plate-pairs’, either rifting apart leading to the formation of conjugate passive margins separated by a new ocean basin, or brought together by collision and orogenesis. Here, we use open-source software tools (GPlates/pyGPlates) to assess quantitative plate kinematics inherent within alternative reconstructions, such as rates of relative plate motion. We analyse the Rodinia-Gondwana transition during the Neoproterozoic, investigating the proposed Australia-Laurentia configurations during Rodinia, and the motion of India colliding with Gondwana. We find that earlier rifting times provide more optimal kinematic results. The AUSWUS and AUSMEX configurations with rifting at 800 Ma are the most kinematically supported configurations for Australia and Laurentia (average rates of 57 and 64 mm/a respectively), and angular rotation of ∼1.4°/Ma, compared to a SWEAT configuration (average spreading rate ∼76 mm/a) and Missing-Link configuration (∼90 mm/a). Later rifting, at, or after, 725 Ma necessitates unreasonably high spreading rates of >130 mm/a for AUSWUS and AUSMEX and ∼150 mm/a for SWEAT and Missing-Link. Using motion paths and convergence rates, we create a kinematically reasonable (convergence below 70 mm/a) tectonic model that is built upon a front-on collision of India with Gondwana, while also incorporating sinistral strike-slip motion against Australia and East Antarctica. We use this simple tectonic model to refine a global model for the breakup of western Rodinia and the transition to eastern Gondwana. © 2017 Elsevier B.V.This manuscript is a contribution to IGCP projects 628 (Gondwana Map) and 648 (Supercontinent Cycles and Global Geodynamics). This research was supported by the Science Industry Endowment Fund (RP 04-174) Big Data Knowledge Discovery Project, Australian Research Council grant DP130101946 (RDM) and the AuScope NCRIS project. ASM is supported by a CSIRO-Data61 Postgraduate Scholarship. ASC's contribution forms TRaX Record #379 and was funded by an Australian Research Council Future Fellowship FT120100340

Kinematic constraints on the Rodinia to Gondwana transition

Earth's plate tectonic history during the breakup of the supercontinent Pangea is well constrained from the seafloor spreading record, but evolving plate configurations during older supercontinent cycles are much less well understood. A relative paucity of available palaeomagnetic and geological data for deep time reconstructions necessitates innovative approaches to help discriminate between competing plate configurations. More difficult is tracing the journeys of individual continents during the amalgamation and breakup of supercontinents. Typically, deep-time reconstructions are built using absolute motions defined by palaeomagnetic data, and do not consider the kinematics of relative motions between plates, even for occasions where they are thought to be ‘plate-pairs’, either rifting apart leading to the formation of conjugate passive margins separated by a new ocean basin, or brought together by collision and orogenesis. Here, we use open-source software tools (GPlates/pyGPlates) to assess quantitative plate kinematics inherent within alternative reconstructions, such as rates of relative plate motion. We analyse the Rodinia-Gondwana transition during the Neoproterozoic, investigating the proposed Australia-Laurentia configurations during Rodinia, and the motion of India colliding with Gondwana. We find that earlier rifting times provide more optimal kinematic results. The AUSWUS and AUSMEX configurations with rifting at 800 Ma are the most kinematically supported configurations for Australia and Laurentia (average rates of 57 and 64 mm/a respectively), and angular rotation of ∼1.4°/Ma, compared to a SWEAT configuration (average spreading rate ∼76 mm/a) and Missing-Link configuration (∼90 mm/a). Later rifting, at, or after, 725 Ma necessitates unreasonably high spreading rates of >130 mm/a for AUSWUS and AUSMEX and ∼150 mm/a for SWEAT and Missing-Link. Using motion paths and convergence rates, we create a kinematically reasonable (convergence below 70 mm/a) tectonic model that is built upon a front-on collision of India with Gondwana, while also incorporating sinistral strike-slip motion against Australia and East Antarctica. We use this simple tectonic model to refine a global model for the breakup of western Rodinia and the transition to eastern Gondwana. © 2017 Elsevier B.V.This manuscript is a contribution to IGCP projects 628 (Gondwana Map) and 648 (Supercontinent Cycles and Global Geodynamics). This research was supported by the Science Industry Endowment Fund (RP 04-174) Big Data Knowledge Discovery Project, Australian Research Council grant DP130101946 (RDM) and the AuScope NCRIS project. ASM is supported by a CSIRO-Data61 Postgraduate Scholarship. ASC's contribution forms TRaX Record #379 and was funded by an Australian Research Council Future Fellowship FT120100340

Dynamic redox and nutrient cycling response to climate forcing in the Mesoproterozoic ocean

Controls on Mesoproterozoic ocean redox heterogeneity, and links to nutrient cycling and oxygenation feedbacks, remain poorly resolved. Here, we report ocean redox and phosphorus cycling across two high-resolution sections from the ~1.4 Ga Xiamaling Formation, North China Craton. In the lower section, fluctuations in trade wind intensity regulated the spatial extent of a ferruginous oxygen minimum zone, promoting phosphorus drawdown and persistent oligotrophic conditions. In the upper section, high but variable continental chemical weathering rates led to periodic fluctuations between highly and weakly euxinic conditions, promoting phosphorus recycling and persistent eutrophication. Biogeochemical modeling demonstrates how changes in geographical location relative to global atmospheric circulation cells could have driven these temporal changes in regional ocean biogeochemistry. Our approach suggests that much of the ocean redox heterogeneity apparent in the Mesoproterozoic record can be explained by climate forcing at individual locations, rather than specific events or step-changes in global oceanic redox conditions

Transient mobilization of subcrustal carbon coincident with Palaeocene–Eocene Thermal Maximum

Plume magmatism and continental breakup led to the opening of the northeast Atlantic Ocean during the globally warm early Cenozoic. This warmth culminated in a transient (170 thousand year, kyr) hyperthermal event associated with a large, if poorly constrained, emission of carbon called the Palaeocene–Eocene Thermal Maximum (PETM) 56 million years ago (Ma). Methane from hydrothermal vents in the coeval North Atlantic Igneous Province (NAIP) has been proposed as the trigger, though isotopic constraints from deep sea sediments have instead implicated direct volcanic carbon dioxide (CO2) emissions. Here we calculate that background levels of volcanic outgassing from mid-ocean ridges and large igneous provinces yield only one-fifth of the carbon required to trigger the hyperthermal. However, geochemical analyses of volcanic sequences spanning the rift-to-drift phase of the NAIP indicate a sudden ~220 kyr-long intensification of magmatic activity coincident with the PETM. This was likely driven by thinning and enhanced decompression melting of the sub-continental lithospheric mantle, which critically contained a high proportion of carbon-rich metasomatic carbonates. Melting models and coupled tectonic–geochemical simulations indicate that >104 gigatons of subcrustal carbon was mobilized into the ocean and atmosphere sufficiently rapidly to explain the scale and pace of the PETM

Closure of the Proterozoic Mozambique Ocean was instigated by a late Tonian plate reorganization event

International audiencePlate reorganization events involve fundamental changes in lithospheric plate-motions and can influence the lithosphere-mantle system as well as both ocean and atmospheric circulation through bathymetric and topographic changes. Here, we compile published data to interpret the geological record of the Neoproterozoic Arabian-Nubian Shield and integrate this with a full-plate tectonic reconstruction. Our model reveals a plate reorganization event in the late Tonian period about 720 million years ago that changed plate-movement directions in the Mozambique Ocean. After the reorganization, Neoproterozoic India moved towards both the African cratons and Australia-Mawson and instigated the future amalgamation of central Gondwana about 200 million years later. This plate kinematic change is coeval with the breakup of the core of Rodinia between Australia-Mawson and Laurentia and Kalahari and Congo. We suggest the plate reorganization event caused the long-term shift of continents to the southern hemisphere and created a pan-northern hemisphere ocean in the Ediacaran

Assembly of the basal mantle structure beneath Africa

International audiencePlate tectonics shapes Earth's surface, and is linked to motions within its deep interior1,2. Cold oceanic lithosphere sinks into the mantle, and hot mantle plumes rise from the deep Earth, leading to volcanism3,4. Volcanic eruptions over the past 320 million years have been linked to two large structures at the base of the mantle presently under Africa and the Pacific Ocean5,6. This has led to the hypothesis that these basal mantle structures have been stationary over geological time7,8, in contrast to observations and models suggesting that tectonic plates9,10, subduction zones11-14 and mantle plumes15,16 have been mobile, and that basal mantle structures are presently deforming17,18. Here we reconstruct mantle flow from one billion years ago to the present day to show that the history of volcanism is statistically as consistent with mobile basal mantle structures as with fixed ones. In our reconstructions, cold lithosphere sank deep into the African hemisphere between 740 and 500 million years ago, and from 400 million years ago the structure beneath Africa progressively assembled, pushed by peri-Gondwana slabs, to become a coherent structure as recently as 60 million years ago. Our mantle flow models suggest that basal mantle structures are mobile, and aggregate and disperse over time, similarly to continents at Earth's surface9. Our models also predict the presence of continental material in the mantle beneath Africa, consistent with geochemical data19,20

Author Correction: Assembly of the basal mantle structure beneath Africa (Nature, (2022), 603, 7903, (846-851), 10.1038/s41586-022-04538-y)

In the version of this article initially published, there were labeling errors in the x-axis tick mark labels for Fig. 3b–d. The Fig. 3b labels now reading “0.70, 0.80, 0.90” appeared initially as “0.65, 0.70, 0.75,” the Fig. 3c labels now reading “–5, 0” originally read “–6, –4”, and the Fig. 3d labels now reading “0, 1.0” originally read “0, 0.5”. The x-axis labels have been corrected in the HTML and PDF versions of the article

Pulsated Global Hydrogen and Methane Flux at Mid-Ocean Ridges Driven by Pangea Breakup

International audienceMolecular hydrogen production occurs through the serpentinization of mantle peridotite exhumed at mid-ocean ridges. Hydrogen is considered essential to sustain microbial life in the subsurface; however, estimates of hydrogen flux through geological time are unknown. Here we present a model of the primary, abiotic production of molecular hydrogen from the serpentinization of oceanic lithosphere using full-plate tectonic reconstructions for the last 200 Ma. We find significant variability in hydrogen fluxes (1-70 • 1016 mol/Ma or 0.2-14.1 • 105 Mt/a), which are a function of the sensitivity of evolving ocean basins to spreading rates and can be correlated with the opening of key ocean basins during the breakup of Pangea. We suggest that the primary driver of this hydrogen flux is the continental reconfiguration during Pangea breakup, as this produces ocean basins more conducive to exhuming and exposing mantle peridotite at slow and ultraslow spreading ridges. Consequently, present-day flux estimates are ~7 • 1017 mol/Ma (1.4 • 106 Mt/a), driven primarily by the slow and ultraslow spreading ridges in the Atlantic, Indian, and Arctic oceans. As methane has also been sampled alongside hydrogen at hydrothermal vents, we estimate the methane flux using methane-to-hydrogen ratios from present-day hydrothermal vent fluids. These ratios suggest that methane flux ranges between 10 and 100% of the total hydrogen flux, although as the release of methane from these systems is still poorly understood, we suggest a lower estimate, equivalent to around 7-12 • 1016 mol/Ma (1.1-1.9 • 107 Mt/Ma) of methane

Evolving Marginal Terranes During Neoproterozoic Supercontinent Reorganization: Constraints From the Bemarivo Domain in Northern Madagascar

International audienceMadagascar is a key area for unraveling the geodynamic evolution of the transition between the Rodinia and Gondwana supercontinents as it contains several suites of c. 850-700 Ma magmatic rocks that have been postulated to correlate with other Rodinian terranes. The Bemarivo Domain of northern Madagascar contains the youngest of these units that date to c. 750-700 Ma. We present zircon Hf and O isotope data to understand northern Madagascar's place in the Neoproterozoic plate tectonic reconfiguration. We demonstrate that the northern component of the Bemarivo Domain is distinct from the southern part of the Bemarivo Domain and have therefore assigned new names—the Bobakindro Terrane and Marojejy Terrane, respectively. Magmatic rocks of the Marojejy Terrane and Anaboriana Belt are characterized by evolved ɛHf(t) signatures and a range of δ18O values, similar to the Imorona-Itsindro Suite of central Madagascar. These magmatic suites likely formed together in the same long-lived volcanic arc. In contrast, the Bobakindro Terrane contains juvenile ɛHf(t) and mantle-like δ18O values, with no probable link to the rest of Madagascar. We propose that the Bobakindro Terrane formed in a juvenile arc system that included the Seychelles, the Malani Igneous Suite of northwest India, Oman, and the Yangtze Belt of south China, which at the time were all outboard from continental India and south China. The final assembly of northern Madagascar and amalgamation of the Bobakindro Terrane and Marojejy Terrane occurred along the Antsaba subduction zone, with collision occurring at c. 540 Ma

Extending full-plate tectonic models into deep time: Linking the Neoproterozoic and the Phanerozoic

Recent progress in plate tectonic reconstructions has seen models move beyond the classical idea of continental drift by attempting to reconstruct the full evolving configuration of tectonic plates and plate boundaries. A particular problem for the Neoproterozoic and Cambrian is that many existing interpretations of geological and palaeomagnetic data have remained disconnected from younger, better-constrained periods in Earth history. An important test of deep time reconstructions is therefore to demonstrate the continuous kinematic viability of tectonic motions across multiple supercontinent cycles. We present, for the first time, a continuous full-plate model spanning 1 Ga to the present-day, that includes a revised and improved model for the Neoproterozoic–Cambrian (1000–520 Ma) that connects with models of the Phanerozoic, thereby opening up pre-Gondwana times for quantitative analysis and further regional refinements. In this contribution, we first summarise methodological approaches to full-plate modelling and review the existing full-plate models in order to select appropriate models that produce a single continuous model. Our model is presented in a palaeomagnetic reference frame, with a newly-derived apparent polar wander path for Gondwana from 540 to 320 Ma, and a global apparent polar wander path from 320 to 0 Ma. We stress, though while we have used palaeomagnetic data when available, the model is also geologically constrained, based on preserved data from past-plate boundaries. This study is intended as a first step in the direction of a detailed and self-consistent tectonic reconstruction for the last billion years of Earth history, and our model files are released to facilitate community development