The Collaborative Seismic Earth Model: Generation 1.

We present a general concept for evolutionary, collaborative, multiscale inversion of geophysical data, specifically applied to the construction of a first-generation Collaborative Seismic Earth Model. This is intended to address the limited resources of individual researchers and the often limited use of previously accumulated knowledge. Model evolution rests on a Bayesian updating scheme, simplified into a deterministic method that honors today's computational restrictions. The scheme is able to harness distributed human and computing power. It furthermore handles conflicting updates, as well as variable parameterizations of different model refinements or different inversion techniques. The first-generation Collaborative Seismic Earth Model comprises 12 refinements from full seismic waveform inversion, ranging from regional crustal- to continental-scale models. A global full-waveform inversion ensures that regional refinements translate into whole-Earth structure.

Evidence for active upper mantle flow in the Atlantic and Indo-Australian realms since the Upper Jurassic from hiatus maps and spreading rate changes.

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.

Continent-scale Hiatus Maps for the Atlantic Realm and Australia since the Upper Jurassic and links to mantle flow induced dynamic topography.

Interregional geological maps hold important information for geodynamic models. Here, we use such maps to visualize major conformable and unconformable contacts at interregional scales and at the level of geologic series from the Upper Jurassic onward across North and South America, Europe, Africa and Australia. We extract hiatus information from these paleogeological maps, which we plot in a paleogeographical reference frame to link the maps to the plate and plume modes of mantle convection. We assume that interregional patterns of hiatus surfaces are proxy records of continent-scale mantle-induced vertical motion of the lithosphere. We find significant differences in the distribution of hiatus across and between continents at the timescale of geologic series, that is ten to a few tens of millions of years (Myrs). This is smaller than the mantle transit time, which, as the timescale of convection, is about 100-200 Myrs. Our results imply that different timescales for convection and topography in convective support must be an integral component of time-dependent geodynamic Earth models, consistent with the presence of a weaker upper mantle relative to the lower mantle. Additional geological constraints together with interregional geological maps at the resolution of stages (1-2 Myrs), are needed to assist in future geodynamic interpretations of interregional geologic hiatus.

Hotspot motion caused the Hawaiian-Emperor Bend and LLSVPs are not fixed.

Controversy surrounds the fixity of both hotspots and large low shear velocity provinces (LLSVPs). Paleomagnetism, plate-circuit analyses, sediment facies, geodynamic modeling, and geochemistry suggest motion of the Hawaiian plume in Earth's mantle during formation of the Emperor seamounts. Herein, we report new paleomagnetic data from the Hawaiian chain (Midway Atoll) that indicate the Hawaiian plume arrived at its current latitude by 28 Ma. A dramatic decrease in distance between Hawaiian-Emperor and Louisville chain seamounts between 63 and 52 Ma confirms a high rate of southward Hawaiian hotspot drift (~47 mm yr-1), and excludes true polar wander as a relevant factor. These findings further indicate that the Hawaiian-Emperor chain bend morphology was caused by hotspot motion, not plate motion. Rapid plume motion was likely produced by ridge-plume interaction and deeper influence of the Pacific LLSVP. When compared to plate circuit predictions, the Midway data suggest ~13 mm yr-1 of African LLSVP motion since the Oligocene. LLSVP upwellings are not fixed, but also wander as they attract plumes and are shaped by deep mantle convection.

The bent Hawaiian-Emperor hotspot track: inheriting the mantle wind.

Bends in volcanic hotspot lineaments, best represented by the large elbow in the Hawaiian-Emperor chain, were thought to directly record changes in plate motion. Several lines of geophysical inquiry now suggest that a change in the locus of upwelling in the mantle induced by mantle dynamics causes bends in hotspot tracks. Inverse modeling suggests that although deep flow near the core-mantle boundary may have played a role in the Hawaiian-Emperor bend, capture of a plume by a ridge, followed by changes in sub-Pacific mantle flow, can better explain the observations. Thus, hotspot tracks can reveal patterns of past mantle circulation.

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.