Maria Seton

Lower mantle structure from paleogeographically constrained dynamic Earth models

Seismic tomography reveals two large, low-shear velocity provinces (LLSVPs) beneath Africa and the Pacific Ocean. These structures may have existed for several 100 Myr and are likely compositionally distinct based on observed seismic characteristics interpreted in light of geodynamic models and mineral physics constraints. We investigate the dynamics of the LLSVPs through the use of evolutionary models of thermochemical convection from 250 Ma to present day. We use a spherical convection model in which the anomalous structures have a high bulk modulus, consistent with seismic interpretation. A new progressive assimilation method incorporates constraints from paleogeography using a refined plate history model (with 1 Myr time spacing) to guide the thermal structure of the lithosphere and steer the thermal evolution of slabs in the uppermost mantle. The thermochemical structures deform and migrate along the core-mantle boundary (CMB) through coupling to plate motions and in response to slab stresses. The models produce a ridge-like anomaly beneath Africa and a rounded pile beneath the Pacific Ocean, which at present day agrees with tomography, waveform modeling, and other geodynamic studies. Plumes emanate from the margins of the domes and ridges of thickened boundary layer between the domes. Dense and viscous slabs can undermine the stability of high bulk modulus structures at the CMB. High bulk modulus structures are not necessarily required to satisfy dynamic constraints on the LLSVPs.

Mid-Cretaceous seafloor spreading pulse: Fact or fiction?

Two main hypotheses compete to explain the mid-Cretaceous global sea-level highstand: a massive pulse of oceanic crustal production that occurred during the Cretaceous Normal Superchron (CNS) and the “supercontinent breakup effect,” which resulted in the creation of the mid-Atlantic and Indian ocean ridges at the expense of subducting old ocean floor in the Tethys and the Pacific. We have used global oceanic paleo-age grids, including now subducted ocean floor and two alternative time scales, to test these hypotheses. Our models show that a high average seafloor spreading rate of 92 mm/a in the Early Cretaceous that decreased to 60 mm/a during the Tertiary, with peaks of 86 mm/a and 70 mm/a at 105 Ma and 75 Ma ago, respectively, correspond to the two observed sea-level highstands in the Cretaceous. Calculations using GTS2004 produce lower seafloor spreading rates during the same period and diminish the mid-Cretaceous spreading pulse. Global ridge lengths increased in the earliest Cretaceous but stayed relatively constant through time. However, we find that the average age of the ocean basins through time is only weakly dependent on the choice of time scale. The expansive mid- and Late Cretaceous epicontinental seas, coupled with warm climates and oxygen-poor water masses, were ultimately driven by the younger average age of the Cretaceous seafloor and faster seafloor spreading rather than a vast increase in mid-ocean ridge length due to the breakup of Pangea or solely by higher seafloor spreading rates, as suggested previously.

Middle Miocene tectonic boundary conditions for use in climate models

Utilizing general circulation models (GCMs) for paleoclimate study requires the construction of appropriate model boundary conditions. We present a middle Miocene paleotopographic and paleobathymetric reconstruction geographically constrained at 15 Ma for use in GCMs. Paleotopography and paleogeography are reconstructed using a published global plate rotation model and published geological data. Paleobathymetry is reconstructed through application of an age-depth relationship to a middle Miocene global digital isochron map, followed by the overlay of reconstructed sediment thickness and large igneous provinces. Adjustments are subsequently made to ensure our reconstruction may be utilized in GCMs.

Ridge subduction sparked reorganization of the Pacific plate-mantle system 60-50 million years ago

A reorganization centered on the Pacific plate occurred ~53–47 million years ago. A “top-down” plate tectonic mechanism, complete subduction of the Izanagi plate, as opposed to a “bottom-up” mantle flow mechanism, has been proposed as the main driver. Verification based on marine geophysical observations is impossible as most ocean crust recording this event has been subducted. Using a forward modeling approach, which assimilates surface plate velocities and shallow thermal structure of slabs into mantle flow models, we show that complete Izanagi plate subduction and margin-wide slab detachment induced a major change in sub-Pacific mantle flow, from dominantly southward before 60 Ma to north-northeastward after 50 Ma. Our results agree with onshore geology, mantle tomography, and the inferred motion of the Hawaiian hot spot and are consistent with a plate tectonic process driving the rapid plate-mantle reorganization in the Pacific hemisphere between 60 and 50 Ma. This reorganization is reflected in tectonic changes in the Pacific and surrounding ocean basins.

Insights on the kinematics of the India‐Eurasia collision from global geodynamic models

The Eocene India-Eurasia collision is a first order tectonic event whose nature and chronology remains controversial. We test two end-member collision scenarios using coupled global plate motion-subduction models. The first, conventional model, invokes a continental collision soon after ∼60 Ma between a maximum extent Greater India and an Andean-style Eurasian margin. The alternative scenario involves a collision between a minimum extent Greater India and a NeoTethyan back-arc at ∼60 Ma that is subsequently subducted along southern Lhasa at an Andean-style margin, culminating with continent-continent contact at ∼40 Ma. Our numerical models suggest the conventional scenario does not adequately reproduce mantle structure related to Tethyan convergence. The alternative scenario better reproduces the discrete slab volumes and their lateral and vertical distribution in the mantle, and is also supported by the distribution of ophiolites indicative of Tethyan intraoceanic subduction, magmatic gaps along southern Lhasa and a two-stage slowdown of India. Our models show a strong component of southward mantle return flow for the Tethyan region, suggesting that the common assumption of near-vertical slab sinking is an oversimplification with significant consequences for interpretations of seismic tomography in the context of subduction reference frames.

Early to Middle Miocene monsoon climate in Australia

The present-day Australian monsoon delivers substantial moisture to the northern regions of a predominantly arid continent. However, the pre-Quaternary history of the Australian monsoon is poorly constrained due to sparse and often poorly dated paleoclimate proxy evidence. Sedimentological and paleontological data suggest that warm, humid, and seasonal environments prevailed in central and north Australia during the Miocene, though it is unclear whether these were products of the Australian monsoon. We perform a series of sensitivity experiments using an atmospheric general circulation model, combined with an offl ine equilibrium vegetation model, to quantitatively constrain the areal extent of the Miocene monsoon. Our results suggest a weaker than modern monsoon climate during the Miocene. This result is insensitive to atmospheric CO 2 , although somewhat sensitive to vegetation interactions and the presumed distribution of inland water bodies. None of our Miocene experiments exhibit precipitation rates greater than modern over north Australia, in disagreement with paleoclimate record interpretations. Vegetation modeling indicates that inferred precipitation values from fossil fl ora and fauna could only support Miocene vegetation patterns if atmospheric CO 2 was twice the modern concentration. This suggests that elevated CO 2 was critical for sustaining Miocene vegetation.

Community infrastructure and repository for marine magnetic identifications

Magnetic anomaly identifications underpin plate tectonic reconstructions and form the primary data set from which the age of the oceanic lithosphere and seafloor spreading regimes in the ocean basins can be determined. Although these identifications are an invaluable resource, their usefulness to the wider scientific community has been limited due to the lack of a central community infrastructure to organize, host, and update these interpretations. We have developed an open-source, community-driven online infrastructure as a repository for quality-checked magnetic anomaly identifications from all ocean basins. We provide a global sample data set that comprises 96,733 individually picked magnetic anomaly identifications organized by ocean basin and publication reference, and provide accompanying Hellingerformat files, where available. Our infrastructure is designed to facilitate research in plate tectonic reconstructions or research that relies on an assessment of plate reconstructions, for both experts and nonexperts alike. To further enhance the existing repository and strengthen its value, we encourage others in the community to contribute to this effort.

Zealandia: Earth’s Hidden Continent

A 4.9 Mkm2 region of the southwest Pacific Ocean is made up of continental crust. The region has elevated bathymetry relative to surrounding oceanic crust, diverse and silica-rich rocks, and relatively thick and low-velocity crustal structure. Its isolation from Australia and large area support its definition as a continent—Zealandia. Zealandia was formerly part of Gondwana. Today it is 94% submerged, mainly as a result of widespread Late Cretaceous crustal thinning preceding supercontinent breakup and consequent isostatic balance. The identification of Zealandia as a geological continent, rather than a collection of continental islands, fragments, and slices, more correctly represents the geology of this part of Earth. Zealandia provides a fresh context Nick Mortimer, GNS Science, Private Bag 1930, Dunedin 9054, New Zealand; Hamish J. Campbell, GNS Science, P.O. Box 30368, Lower Hutt 5040, New Zealand; Andy J. Tulloch, GNS Science, Private Bag 1930, Dunedin 9054, New Zealand; Peter R. King, Vaughan M. Stagpoole, Ray A. Wood, Mark S. Rattenbury, GNS Science, P.O. Box 30368, Lower Hutt 5040, New Zealand; Rupert Sutherland, SGEES, Victoria University of Wellington, P.O. Box 600, Wellington 6140, New Zealand; Chris J. Adams, GNS Science, Private Bag 1930, Dunedin 9054, New Zealand; Julien Collot, Service Géologique de Nouvelle Calédonie, B.P. 465, Nouméa 98845, New Caledonia; and Maria Seton, School of Geosciences, University of Sydney, NSW 2006, Australia in which to investigate processes of continental rifting, thinning, and breakup.

The GPlates Portal: Cloud-Based Interactive 3D Visualization of Global Geophysical and Geological Data in a Web Browser.

The pace of scientific discovery is being transformed by the availability of ‘big data’ and open access, open source software tools. These innovations open up new avenues for how scientists communicate and share data and ideas with each other and with the general public. Here, we describe our efforts to bring to life our studies of the Earth system, both at present day and through deep geological time. The GPlates Portal (portal.gplates.org) is a gateway to a series of virtual globes based on the Cesium Javascript library. The portal allows fast interactive visualization of global geophysical and geological data sets, draped over digital terrain models. The globes use WebGL for hardware-accelerated graphics and are cross-platform and cross-browser compatible with complete camera control. The globes include a visualization of a high-resolution global digital elevation model and the vertical gradient of the global gravity field, highlighting small-scale seafloor fabric such as abyssal hills, fracture zones and seamounts in unprecedented detail. The portal also features globes portraying seafloor geology and a global data set of marine magnetic anomaly identifications. The portal is specifically designed to visualize models of the Earth through geological time. These space-time globes include tectonic reconstructions of the Earth’s gravity and magnetic fields, and several models of long-wavelength surface dynamic topography through time, including the interactive plotting of vertical motion histories at selected locations. The globes put the on-the-fly visualization of massive data sets at the fingertips of end-users to stimulate teaching and learning and novel avenues of inquiry.

Circum‐Arctic mantle structure and long‐wavelength topography since the Jurassic

The circum-Arctic is one of the most tectonically complex regions of the world, shaped by a history of ocean basin opening and closure since the Early Jurassic. The region is characterized by contemporaneous large-scale Cenozoic exhumation extending from Alaska to the Atlantic, but its driving force is unknown. We show that the mantle flow associated with subducted slabs of the South Anuyi, Mongol-Okhotsk, and Panthalassa oceans have imparted long-wavelength deflection on overriding plates. We identify the Jurassic-Cretaceous South Anuyi slab under present-day Greenland in seismic tomography and numerical mantle flow models. Under North America, we propose the “Farallon” slab results from Andean-style ocean-continent convergence around ~30°N and from a combination of ocean-continent and intraoceanic subduction north of 50°N. We compute circum-Arctic dynamic topography through time from subduction-driven convection models and find that slabs have imparted on average <1–16 m/Myr of dynamic subsidence across the region from at least 170 Ma to ~50 Ma. With the exception of Siberia, the main phase of circum-Arctic dynamic subsidence has been followed either by slowed subsidence or by uplift of <1–6 m/Myr on average to present day. Comparing these results to geological inferences suggest that subduction-driven dynamic topography can account for rapid Middle to Late Jurassic subsidence in the Slave Craton and North Slope (respectively, <15 and 21 m/Myr, between 170 and 130 Ma) and for dynamic subsidence (<7 m/Myr, ~170–50 Ma) followed by dynamic uplift (<6 m/Myr since 50 Ma) of the Barents Sea region. Combining detailed kinematic reconstructions with geodynamic modeling and key geological observations constitutes a powerful tool to investigate the origin of vertical motion in remote regions.