Mesozoic plate-motion history below the northeast Pacific Ocean from seismic images of the subducted Farallon slab

Assessing plate reconstruction models using plate driving force consistency tests

Plate reconstruction models are constructed to fit constraints such as magnetic anomalies, fracture zones, paleomagnetic poles, geological observations and seismic tomography. However, these models do not consider the physical equations of plate driving forces when reconstructing plate motion. This can potentially result in geodynamically-implausible plate motions, which has implications for a range of work based on plate reconstruction models. We present a new algorithm that calculates time-dependent slab pull, ridge push (GPE force) and mantle drag resistance for any topologically closed reconstruction, and evaluates the residuals-or missing components-required for torques to balance given our assumed plate driving force relationships. In all analyzed models, residual torques for the present-day are three orders of magnitude smaller than the typical driving torques for oceanic plates, but can be of the same order of magnitude back in time-particularly from 90 to 50 Ma. Using the Pacific plate as an example, we show how our algorithm can be used to identify areas and times with high residual torques, where either plate reconstructions have a high degree of geodynamic implausibility or our understanding of the underlying geodynamic forces is incomplete. We suggest strategies for plate model improvements and also identify times when other forces such as active mantle flow were likely important contributors. Our algorithm is intended as a tool to help assess and improve plate reconstruction models based on a transparent and expandable set of a priori dynamic constraints.

Glacial isostatic adjustment: physical models and observational constraints

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.

Subduction history of the Caribbean from upper-mantle seismic imaging and plate reconstruction

The margins of the Caribbean and associated hazards and resources have been shaped by a poorly understood history of subduction. Using new data, we improve teleseismic P-wave imaging of the eastern Caribbean upper mantle and compare identified subducted-plate fragments with trench locations predicted from plate reconstruction. This shows that material at 700-1200 km depth below South America derives from 90-115 Myr old westward subduction, initiated prior to Caribbean Large-Igneous-Province volcanism. At shallower depths, an accumulation of subducted material is attributed to Great Arc of the Caribbean subduction as it evolved over the past 70 Ma. We interpret gaps in these subducted-plate anomalies as: a plate window and tear along the subducted Proto-Caribbean ridge; tearing along subducted fracture zones, and subduction of a volatile-rich boundary between Proto-Caribbean and Atlantic domains. Phases of back-arc spreading and arc jumps correlate with changes in age, and hence buoyancy, of the subducting plate.

Imaging paleoslabs in the D″ layer beneath Central America and the Caribbean using seismic waveform inversion

D″ (Dee double prime), the lowermost layer of the Earth's mantle, is the thermal boundary layer (TBL) of mantle convection immediately above the Earth's liquid outer core. As the origin of upwelling of hot material and the destination of paleoslabs (downwelling cold slab remnants), D″ plays a major role in the Earth's evolution. D″ beneath Central America and the Caribbean is of particular geodynamical interest, because the paleo- and present Pacific plates have been subducting beneath the western margin of Pangaea since ~250 million years ago, which implies that paleoslabs could have reached the lowermost mantle. We conduct waveform inversion using a data set of ~7700 transverse component records to infer the detailed three-dimensional S-velocity structure in the lowermost 400 km of the mantle in the study region so that we can investigate how cold paleoslabs interact with the hot TBL above the core-mantle boundary (CMB). We can obtain high-resolution images because the lowermost mantle here is densely sampled by seismic waves due to the full deployment of the USArray broadband seismic stations during 2004-2015. We find two distinct strong high-velocity anomalies, which we interpret as paleoslabs, just above the CMB beneath Central America and Venezuela, respectively, surrounded by low-velocity regions. Strong low-velocity anomalies concentrated in the lowermost 100 km of the mantle suggest the existence of chemically distinct denser material connected to low-velocity anomalies in the lower mantle inferred by previous studies, suggesting that plate tectonics on the Earth's surface might control the modality of convection in the lower mantle.

Intra-oceanic subduction shaped the assembly of Cordilleran North America

The western quarter of North America consists of accreted terranes--crustal blocks added over the past 200 million years--but the reason for this is unclear. The widely accepted explanation posits that the oceanic Farallon plate acted as a conveyor belt, sweeping terranes into the continental margin while subducting under it. Here we show that this hypothesis, which fails to explain many terrane complexities, is also inconsistent with new tomographic images of lower-mantle slabs, and with their locations relative to plate reconstructions. We offer a reinterpretation of North American palaeogeography and test it quantitatively: collision events are clearly recorded by slab geometry, and can be time calibrated and reconciled with plate reconstructions and surface geology. The seas west of Cretaceous North America must have resembled today's western Pacific, strung with island arcs. All proto-Pacific plates initially subducted into almost stationary, intra-oceanic trenches, and accumulated below as massive vertical slab walls. Above the slabs, long-lived volcanic archipelagos and subduction complexes grew. Crustal accretion occurred when North America overrode the archipelagos, causing major episodes of Cordilleran mountain building.

Origin of Columbia River flood basalt controlled by propagating rupture of the Farallon slab

The origin of the Steens-Columbia River (SCR) flood basalts, which is presumed to be the onset of Yellowstone volcanism, has remained controversial, with the proposed conceptual models involving either a mantle plume or back-arc processes. Recent tomographic inversions based on the USArray data reveal unprecedented detail of upper-mantle structures of the western USA and tightly constrain geodynamic models simulating Farallon subduction, which has been proposed to influence the Yellowstone volcanism. Here we show that the best-fitting geodynamic model depicts an episode of slab tearing about 17 million years ago under eastern Oregon, where an associated sub-slab asthenospheric upwelling thermally erodes the Farallon slab, leading to formation of a slab gap at shallow depth. Driven by a gradient of dynamic pressure, the tear ruptured quickly north and south and within about two million years covering a distance of around 900 kilometres along all of eastern Oregon and northern Nevada. This tear would be consistent with the occurrence of major volcanic dikes during the SCR-Northern Nevada Rift flood basalt event both in space and time. The model predicts a petrogenetic sequence for the flood basalt with sources of melt starting from the base of the slab, at first remelting oceanic lithosphere and then evolving upwards, ending with remelting of oceanic crust. Such a progression helps to reconcile the existing controversies on the interpretation of SCR geochemistry and the involvement of the putative Yellowstone plume. Our study suggests a new mechanism for the formation of large igneous provinces.

The Emperor Seamounts: southward motion of the Hawaiian hotspot plume in Earth's mantle

The Hawaiian-Emperor hotspot track has a prominent bend, which has served as the basis for the theory that the Hawaiian hotspot, fixed in the deep mantle, traced a change in plate motion. However, paleomagnetic and radiometric age data from samples recovered by ocean drilling define an age-progressive paleolatitude history, indicating that the Emperor Seamount trend was principally formed by the rapid motion (over 40 millimeters per year) of the Hawaiian hotspot plume during Late Cretaceous to early-Tertiary times (81 to 47 million years ago). Evidence for motion of the Hawaiian plume affects models of mantle convection and plate tectonics, changing our understanding of terrestrial dynamics.

Mantle-circulation models with sequential data assimilation: inferring present-day mantle structure from plate-motion histories

Data assimilation is an approach to studying geodynamic models consistent simultaneously with observables and the governing equations of mantle flow. Such an approach is essential in mantle circulation models, where we seek to constrain an unknown initial condition some time in the past, and thus cannot hope to use first-principles convection calculations to infer the flow history of the mantle. One of the most important observables for mantle-flow history comes from models of Mesozoic and Cenozoic plate motion that provide constraints not only on the surface velocity of the mantle but also on the evolution of internal mantle-buoyancy forces due to subducted oceanic slabs. Here we present five mantle circulation models with an assimilated plate-motion history spanning the past 120 Myr, a time period for which reliable plate-motion reconstructions are available. All models agree well with upper- and mid-mantle heterogeneity imaged by seismic tomography. A simple standard model of whole-mantle convection, including a factor 40 viscosity increase from the upper to the lower mantle and predominantly internal heat generation, reveals downwellings related to Farallon and Tethys subduction. Adding 35% bottom heating from the core has the predictable effect of producing prominent high-temperature anomalies and a strong thermal boundary layer at the base of the mantle. Significantly delaying mantle flow through the transition zone either by modelling the dynamic effects of an endothermic phase reaction or by including a steep, factor 100, viscosity rise from the upper to the lower mantle results in substantial transition-zone heterogeneity, enhanced by the effects of trench migration implicit in the assimilated plate-motion history. An expected result is the failure to account for heterogeneity structure in the deepest mantle below 1500 km, which is influenced by Jurassic plate motions and thus cannot be modelled from sequential assimilation of plate motion histories limited in age to the Cretaceous. This result implies that sequential assimilation of past plate-motion models is ineffective in studying the temporal evolution of core-mantle-boundary heterogeneity, and that a method for extrapolating present-day information backwards in time is required. For short time periods (of the order of perhaps a few tens of Myr) such a method exists in the form of crude 'backward' convection calculations. For longer time periods (of the order of a mantle overturn), a rigorous approach to extrapolating information back in time exists in the form of iterative nonlinear optimization methods that carry assimilated information into the past through the use of an adjoint mantle convection model.