Oceanic crustal carbon cycle drives 26-million-year atmospheric carbon dioxide periodicities

Atmospheric carbon dioxide (CO2) data for the last 420 million years (My) show long-term fluctuations related to supercontinent cycles as well as shorter cycles at 26 to 32 My whose origin is unknown. Periodicities of 26 to 30 My occur in diverse geological phenomena including mass extinctions, flood basalt volcanism, ocean anoxic events, deposition of massive evaporites, sequence boundaries, and orogenic events and have previously been linked to an extraterrestrial mechanism. The vast oceanic crustal carbon reservoir is an alternative potential driving force of climate fluctuations at these time scales, with hydrothermal crustal carbon uptake occurring mostly in young crust with a strong dependence on ocean bottom water temperature. We combine a global plate model and oceanic paleo-age grids with estimates of paleo-ocean bottom water temperatures to track the evolution of the oceanic crustal carbon reservoir over the past 230 My. We show that seafloor spreading rates as well as the storage, subduction, and emission of oceanic crustal and mantle CO2 fluctuate with a period of 26 My. A connection with seafloor spreading rates and equivalent cycles in subduction zone rollback suggests that these periodicities are driven by the dynamics of subduction zone migration. The oceanic crust-mantle carbon cycle is thus a previously overlooked mechanism that connects plate tectonic pulsing with fluctuations in atmospheric carbon and surface environments. Copyright © 2018 The Authors

Dynamic topography of passive continental margins and their hinterlands since the Cretaceous

Even though it is well accepted that the Earth\u27s surface topography has been affected by mantle-convection induced dynamic topography, its magnitude and time-dependence remain controversial. The dynamic influence to topographic change along continental margins is particularly difficult to unravel, because their stratigraphic record is dominated by tectonic subsidence caused by rifting. We follow a three-fold approach to estimate dynamic topographic change along passive margins based on a set of seven global mantle convection models. We first demonstrate that a geodynamic forward model that includes adiabatic and viscous heating in addition to internal heating from radiogenic sources, and a mantle viscosity profile with a gradual increase in viscosity below the mantle transition zone, provides a greatly improved match to the spectral range of residual topography end-members as compared with previous models at very long wavelengths (spherical degrees 2-3). We then combine global sea level estimates with predicted surface dynamic topography to evaluate the match between predicted continental flooding patterns and published paleo-coastlines by comparing predicted versus geologically reconstructed land fractions and spatial overlaps of flooded regions for individual continents since 140 Ma. Modelled versus geologically reconstructed land fractions match within 10% for most models, and the spatial overlaps of inundated regions are mostly between 85% and 100% for the Cenozoic, dropping to about 75-100% in the Cretaceous. Regions that have been strongly affected by mantle plumes are generally not captured well in our models, as plumes are suppressed in most of them, and our models with dynamically evolving plumes do not replicate the location and timing of observed plume products. We categorise the evolution of modelled dynamic topography in both continental interiors and along passive margins using cluster analysis to investigate how clusters of similar dynamic topography time series are distributed spatially. A subdivision of four clusters is found to best reveal end-members of dynamic topography evolution along passive margins and their hinterlands, differentiating topographic stability, long-term pronounced subsidence, initial stability over a dynamic high followed by moderate subsidence and regions that are relatively proximal to subduction zones with varied dynamic topography histories. Along passive continental margins the most commonly observed process is a gradual motion from dynamic highs towards lows during the fragmentation of Pangea, reflecting the location of many passive margins now over slabs sinking in the lower mantle. Our best-fit model results in up to 500 (± 150) m of total dynamic subsidence of continental interiors while along passive margins the maximum predicted dynamic topographic change over 140 million years is about 350 (± 150) m of subsidence. Models with plumes exhibit clusters of transient passive margin uplift of about 200 ± 200 m, but are mainly characterised by long-term subsidence of up to 400 m. The good overall match between predicted dynamic topography to geologically mapped paleo-coastlines makes a convincing case that mantle-driven topographic change is a critical component of relative sea level change, and indeed the main driving force for generating the observed geometries and timings of large-scale continental inundation through time

Influence of mantle flow on the drainage of eastern Australia since the Jurassic Period

Recent studies of the past eastern Australian landscape from present-day longitudinal river profiles and from mantle flow models suggest that the interaction of plate motion with mantle convection accounts for the two phases of large-scale uplift of the region since 120 Ma. We coupled the dynamic topography predicted from one of these mantle flow models to a surface process model to study the evolution of the eastern Australian landscape since the Jurassic Period. We varied the rainfall regime, erodibility, sea level variations, dynamic topography magnitude, and elastic thickness across a series of experiments. The approach accounts for erosion and sedimentation and simulates catchment dynamics. Despite the relative simplicity of our model, the results provide insights on the fundamental links between dynamic topography and continental-scale drainage evolution. Based on temporal and spatial changes in longitudinal river profiles as well as erosion and deposition maps, we show that the motion of the Australian plate over the convecting mantle has resulted in significant reorganization of the eastern Australian drainage. The model predicts that the Murray river drained eastward between 150 and ∼120 Ma, and switched to westward draining due to the tilting of the Australian plate from ∼120 Ma. First order comparisons of eight modeled river profiles and of the catchment shape of modeled Murray-Darling Basin are in agreement with present-day observations. The predicted denudation of the eastern highlands is compatible with thermochronology data and sedimentation rates along the southern Australian margin are consistent with cumulative sediment thickness. © 2017. American Geophysical Union. All Rights Reserved.We thank Robert Moucha, Nicky White and an anonymous reviewer for their constructive comments and suggestions which have greatly improve this paper. The authors were supported by ARC grants IH130200012 and DE160101020

Special Libraries, May 1916

Volume 7, Issue 5https://scholarworks.sjsu.edu/sla_sl_1916/1004/thumbnail.jp

Oblique rifting: The rule, not the exception

Movements of tectonic plates often induce oblique deformation at divergent plate boundaries. This is in striking contrast with traditional conceptual models of rifting and rifted margin formation, which often assume 2-D deformation where the rift velocity is oriented perpendicular to the plate boundary. Here we quantify the validity of this assumption by analysing the kinematics of major continent-scale rift systems in a global plate tectonic reconstruction from the onset of Pangea breakup until the present day. We evaluate rift obliquity by joint examination of relative extension velocity and local rift trend using the script-based plate reconstruction software pyGPlates. Our results show that the global mean rift obliquity since 230 Ma amounts to 34° with a standard deviation of 24°, using the convention that the angle of obliquity is spanned by extension direction and rift trend normal. We find that more than ∼ 70% of all rift segments exceeded an obliquity of 20° demonstrating that oblique rifting should be considered the rule, not the exception. In many cases, rift obliquity and extension velocity increase during rift evolution (e.g. Australia-Antarctica, Gulf of California, South Atlantic, India-Antarctica), which suggests an underlying geodynamic correlation via obliquity-dependent rift strength. Oblique rifting produces 3-D stress and strain fields that cannot be accounted for in simplified 2-D plane strain analysis. We therefore highlight the importance of 3-D approaches in modelling, surveying, and interpretation of most rift segments on Earth where oblique rifting is the dominant mode of deformation. © 2018 Author(s).Acknowledgements. This research has been funded by the German Academic Exchange Service (DAAD), project no. 57319603. Sascha Brune was supported through the Helmholtz Young Investigators Group CRYSTALS (VH-NG-1132). Simon E. Willliams and R. Dietmar Müller were supported by Australian Research Council grant IH130200012. We thank two anonymous reviewers and editor Federico Rossetti for their constructive and motivating comments that significantly helped to improve this manuscript

The role of deep Earth dynamics in driving the flooding and emergence of New Guinea since the Jurassic

The paleogeography of New Guinea indicates fluctuating periods of flooding and emergence since the Jurassic, which are inconsistent with estimates of global sea level change since the Eocene. The role of deep Earth dynamics in explaining these discrepancies has not been explored, despite the strongly time-dependent geodynamic setting within which New Guinea has evolved. We aim to investigate the role of subduction-driven mantle flow in controlling long-wavelength dynamic topography and its manifestation in the regional sedimentary record, within a tectonically complex region leading to orogeny. We couple regionally refined global plate reconstructions with forward geodynamic models to compare trends of dynamic topography with estimates of eustasy and regional paleogeography. Qualitative corroboration of modelled mantle structure with equivalent tomographic profiles allows us to ground-truth the models. We show that predicted dynamic topography correlates with the paleogeographic record of New Guinea from the Jurassic to the present. We find that subduction at the East Gondwana margin locally enhanced the high eustatic sea levels from the Early Cretaceous (∼145 Ma) to generate long-term regional flooding. During the Miocene, however, dynamic subsidence associated with subduction of the Maramuni Arc played a fundamental role in causing long-term inundation of New Guinea during a period of global sea level fall. © 2017 Elsevier B.VThis research was undertaken with the assistance of resources from the National Computational Infrastructure (NCI), which is supported by the Australian Government . SZ and RDM were supported by Australian Research Council grant IH130200012 and DP130101946 . NF was supported by Australian Research Council grant DE16010102

The role of deep Earth dynamics in driving the flooding and emergence of New Guinea since the Jurassic

The paleogeography of New Guinea indicates fluctuating periods of flooding and emergence since the Jurassic, which are inconsistent with estimates of global sea level change since the Eocene. The role of deep Earth dynamics in explaining these discrepancies has not been explored, despite the strongly time-dependent geodynamic setting within which New Guinea has evolved. We aim to investigate the role of subduction-driven mantle flow in controlling long-wavelength dynamic topography and its manifestation in the regional sedimentary record, within a tectonically complex region leading to orogeny. We couple regionally refined global plate reconstructions with forward geodynamic models to compare trends of dynamic topography with estimates of eustasy and regional paleogeography. Qualitative corroboration of modelled mantle structure with equivalent tomographic profiles allows us to ground-truth the models. We show that predicted dynamic topography correlates with the paleogeographic record of New Guinea from the Jurassic to the present. We find that subduction at the East Gondwana margin locally enhanced the high eustatic sea levels from the Early Cretaceous (∼145 Ma) to generate long-term regional flooding. During the Miocene, however, dynamic subsidence associated with subduction of the Maramuni Arc played a fundamental role in causing long-term inundation of New Guinea during a period of global sea level fall. © 2017 Elsevier B.VThis research was undertaken with the assistance of resources from the National Computational Infrastructure (NCI), which is supported by the Australian Government . SZ and RDM were supported by Australian Research Council grant IH130200012 and DP130101946 . NF was supported by Australian Research Council grant DE16010102

Oceanic microplate formation records the onset of India-Eurasia collision

Mapping of seafloor tectonic fabric in the Indian Ocean, using high-resolution satellite-derived vertical gravity gradient data, reveals an extinct Pacific-style oceanic microplate ('Mammerickx Microplate') west of the Ninetyeast Ridge. It is one of the first Pacific-style microplates to be mapped outside the Pacific basin, suggesting that geophysical conditions during formation probably resembled those that have dominated at eastern Pacific ridges. The microplate formed at the Indian-Antarctic ridge and is bordered by an extinct ridge in the north and pseudofault in the south, whose conjugate is located north of the Kerguelen Plateau. Independent microplate rotation is indicated by asymmetric pseudofaults and rotated abyssal hill fabric, also seen in multibeam data. Magnetic anomaly picks and age estimates calculated from published spreading rates suggest formation during chron 21o (~47.3 Ma). Plate reorganizations can trigger ridge propagation and microplate development, and we propose that Mammerickx Microplate formation is linked with the India-Eurasia collision (initial 'soft' collision). The collision altered the stress regime at the Indian-Antarctic ridge, leading to a change in segmentation and ridge propagation from an establishing transform. Fast Indian-Antarctic spreading that preceded microplate formation, and Kerguelen Plume activity, may have facilitated ridge propagation via the production of thin and weak lithosphere; however both factors had been present for tens of millions of years and are therefore unlikely to have triggered the event. Prior to the collision, the combination of fast spreading and plume activity was responsible for the production of a wide region of undulate seafloor to the north of the extinct ridge and 'W' shaped lineations that record back and forth ridge propagation. Microplate formation provides a precise means of dating the onset of the India-Eurasia collision, and is completely independent of and complementary to timing constraints derived from continental geology or convergence histories. © 2015 Elsevier B.V.K.J.M. and R.D.M. were supported by ARC Discovery Project DP130101946 . The CryoSat-2 data were provided by the European Space Agency, and NASA/CNES provided data from the Jason-1 altimeter. This research was supported by the National Science Foundation ( OCE-1128801 ), the Office of Naval Research ( N00014-12-1-0111 ), the National Geospatial-Intelligence Agency ( HM0177-13-1-0008 ) and ConocoPhillips . Version 23 of the global grids of gravity anomaly and VGG can be downloaded from the following ftp site ftp://topex.ucsd.edu/pub/global_grav_1min . All figures were produced using the Generic Mapping Tools ( GMT ) software ( Wessel et al., 2013 ). The open-source plate reconstruction software GPlates ( Boyden et al., 2011 ) was used to compute the distance from the Kerguelen Plume to the initiation point of ridge propagation using different absolute reference frames, and to produce the VGG raster reconstruction in Fig. 3 . Magnetic anomaly picks were accessed from the compilation of Seton et al. (2014) from The Global Seafloor Fabric and Magnetic Lineation Data Base Project website ( http://www.soest.hawaii.edu/PT/GSFML/ ). We thank the editor An Yin and two anonymous reviewers for their thoughtful and constructive comments that improved the manuscript. Appendix

Geodynamic reconstruction of an accreted Cretaceous back-arc basin in the Northern Andes

A complex history of subduction, back-arc basin formation, terrane accretion and transpressional shearing characterizes the evolution of the Caribbean and northern South American margin since Jurassic times. Quantitative plate tectonic reconstructions of the area do not include Jurassic-Cretaceous back-arc terranes of which there are both geological and geophysical observations. We developed a revised plate tectonic reconstruction based on geological observations and seismic tomography models to constrain the Jurassic-Cretaceous subduction history of eastern Panthalassa, along the western margin of the Caribbean region. This reconstruction considers the opening of a Northern Andean back-arc basin at 145 Ma, the Quebradagrande back-arc, closing at 120 Ma and followed by terrane accretion and northward translation along the South American margin starting at 100 Ma. This kinematic reconstruction is tested against two previously published tectonic reconstructions via coupling with global numerical mantle convection models using CitcomS. A comparison of modelled versus tomographically imaged mantle structure reveals that subduction outboard of the South American margin, lacking in previous tectonic models, is required to reproduce mid-mantle positive seismic anomalies imaged in P- and S-wave seismic tomography beneath South America, 500–2000 km in depth. Furthermore, we show that this subduction zone is likely produced by a back-arc basin that developed along the northern Andes during the Cretaceous via trench roll-back from 145 Ma and was closed at 100 Ma. The contemporaneous opening of the Quebradagrande back-arc basin with the Rocas Verdes back-arc basin in the southern Andes is consistent with a model that invokes return flow of mantle material behind a retreating slab and may explain why extension along the Peruvian and Chilean sections of the Andean margin did not experience full crustal break-up and back-arc opening during the late Jurassic-early Cretaceous Period.This research was undertaken with the assistance of the Sydney Informatics Hub in accessing resources from the National Computational Infrastructure (NCI), which is supported by the Australian Government. The authors were supported by Australian Research Council grants IH130200012 (RDM, NF), FT130101564 (MS) and DE160101020 (NF). We would like to thank Wouter Schellart and an anonymous reviewer for their helpful comments, which improved the manuscript

Oceanic Residual Topography Agrees With Mantle Flow Predictions at Long Wavelengths

Dynamic topography, the surface deflection induced by sublithosheric mantle flow, is an important prediction made by geodynamic models, but there is an apparent disparity between geodynamic model predictions and estimates of residual topography (total topography minus lithospheric and crustal contributions). We generate synthetic global topography fields with different power spectral slopes and spatial patterns to investigate how well the long-wavelength (spherical degrees 1 to 3) components can be recovered from a discrete set of samples where residual topography has been recently estimated. An analysis of synthetic topography, along with observed geoid and gravity anomalies, demonstrates the reliability of signal recovery. Appropriate damping factors, which depend on the maximum degree in the spherical harmonic expansion that is used to fit the samples, must be applied to recover the long-wavelength topography correctly; large damping factors smooth the model excessively and suppress residual topography amplitude and power spectra unrealistically. Recovered long-wavelength residual topographies based on recent oceanic point-wise estimates with different spherical expansion degrees agree with each other and with the predicted dynamic topography from mantle flow models. The peak amplitude of the long-wavelength residual topography from oceanic observations is about 1 km, suggesting an important influence of large-scale deep mantle flow. ©2017. American Geophysical Union. All Rights Reserved.T.Y. benefitted from the discussion with Judith Sippel on residual topography. The authors thank Malcolm Sambridge and two anonymous reviewers for reading the original manuscript and providing insightful suggestions. M.G. has been supported by the National Science Foundation through EAR-1358646, EAR-1600956, and EAR-1645775 and by Statoil ASA. L.M. and R. D.M. were supported by Australian Research Council grants DP130101946 and IH130200012. Dynamic topography and the recovered long-wavelength residual topography data are listed in the supporting information