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.

The adjoint equations for thermochemical compressible mantle convection: derivation and verification by twin experiments

The adjoint method is an efficient way to obtain gradient information in a mantle convection model relative to past flow structure, allowing one to retrodict mantle flow from observations of the present-day mantle state. While adjoint equations for isochemical mantle flow have been derived for both incompressible and compressible flows, here we extend the method to thermochemical mantle flow models, and present thermochemical adjoint equations in the elastic-liquid approximation. We verify the method with twin experiments, and retrodict the flow history of a thermochemical reference model (reference twin) assuming for the final state, either a consistent thermochemical interpretation, using the thermochemical adjoint equations, or an inconsistent purely thermal interpretation, using the isochemical adjoint equations. The consistent simulation correctly retrodicts the flow evolution of the reference twin. The inconsistent case, instead, restores a false flow history whereby internal buoyancy forces and convectively maintained topography are overestimated. Because the cost function is reduced in either case, our results suggest that the adjoint method can be used to link assumptions on the role of chemical mantle heterogeneity to geologic inferences of dynamic topography, thus providing additional means to test hypotheses on mantle composition and dynamics.

On the consistency of seismically imaged lower mantle slabs

The geoscience community is increasingly utilizing seismic tomography to interpret mantle heterogeneity and its links to past tectonic and geodynamic processes. To assess the robustness and distribution of positive seismic anomalies, inferred as subducted slabs, we create a set of vote maps for the lower mantle with 14 global P-wave or S-wave tomography models. Based on a depth-dependent threshold metric, an average of 20% of any given tomography model depth is identified as a potential slab. However, upon combining the 14 models, the most consistent positive wavespeed features are identified by an increasing vote count. An overall peak in the most robust anomalies is found between 1000-1400 km depth, followed by a decline to a minimum around 2000 km. While this trend could reflect reduced tomographic resolution in the middle mantle, we show that it may alternatively relate to real changes in the time-dependent subduction flux and/or a mid-lower mantle viscosity increase. An apparent secondary peak in agreement below 2500 km depth may reflect the degree-two lower mantle slow seismic structures. Vote maps illustrate the potential shortcomings of using a limited number or type of tomography models and slab threshold criteria.

Effect of cation substitution on bridgmanite elasticity: A key to interpret seismic anomalies in the lower mantle

Seismological observations show that, in some regions of the lower mantle, an increase in bulk sound velocity, interestingly, occurs in the same volume where there is a decrease in shear velocity. We show that this anti-correlated behavior occurs on cation substitution in bridgmanite by making single crystal elasticity measurements of MgSiO₃ and (Mg,Fe,Al)(Si,Al)O₃ using inelastic x-ray scattering in the ambient conditions. Cation substitution of ferrous iron and aluminum may explain large low shear velocity provinces in the lower mantle.

Stability of active mantle upwelling revealed by net characteristics of plate tectonics

Viscous convection within the mantle is linked to tectonic plate motions and deforms Earth's surface across wide areas. Such close links between surface geology and deep mantle dynamics presumably operated throughout Earth's history, but are difficult to investigate for past times because the history of mantle flow is poorly known. Here we show that the time dependence of global-scale mantle flow can be deduced from the net behaviour of surface plate motions. In particular, we tracked the geographic locations of net convergence and divergence for harmonic degrees 1 and 2 by computing the dipole and quadrupole moments of plate motions from tectonic reconstructions extended back to the early Mesozoic era. For present-day plate motions, we find dipole convergence in eastern Asia and quadrupole divergence in both central Africa and the central Pacific. These orientations are nearly identical to the dipole and quadrupole orientations of underlying mantle flow, which indicates that these 'net characteristics' of plate motions reveal deeper flow patterns. The positions of quadrupole divergence have not moved significantly during the past 250 million years, which suggests long-term stability of mantle upwelling beneath Africa and the Pacific Ocean. These upwelling locations are positioned above two compositionally and seismologically distinct regions of the lowermost mantle, which may organize global mantle flow as they remain stationary over geologic time.

Indoor seismology by probing the Earth's interior by using sound velocity measurements at high pressures and temperatures

The adiabatic bulk (K(S)) and shear (G) moduli of mantle materials at high pressure and temperature can be obtained directly by measuring compressional and shear wave velocities in the laboratory with experimental techniques based on physical acoustics. We present the application of the current state-of-the-art experimental techniques by using ultrasonic interferometry in conjunction with synchrotron x radiation to study the elasticity of olivine and pyroxenes and their high-pressure phases. By using these updated thermoelasticity data for these phases, velocity and density profiles for a pyrolite model are constructed and compared with radial seismic models. We conclude that pyrolite provides an adequate explanation of the major seismic discontinuities at 410- and 660-km depths, the gradient in the transition zone, as well as the velocities in the lower mantle, if the uncertainties in the modeling and the variations in different seismic models are considered. The characteristics of the seismic scaling factors in response to thermal anomalies suggest that anticorrelations between bulk sound and shear wave velocities, as well as the large positive density anomalies observed in the lower mantle, cannot be explained fully without invoking chemical variations.

The physical and chemical composition of the lower mantle

This article reviews some of the recent advances made within the field of mineral physics. In order to link the observed seismic and density structures of the lower mantle with a particular mineral composition, knowledge of the thermodynamic properties of the candidate materials is required. Determining which compositional model best matches the observed data is difficult because of the wide variety of possible mineral structures and compositions. State-of-the-art experimental and analytical techniques have pushed forward our knowledge of mineral physics, yet certain properties, such as the elastic properties of lower mantle minerals at high pressures and temperatures, are difficult to determine experimentally and remain elusive. Fortunately, computational techniques are now sufficiently advanced to enable the prediction of these properties in a self-consistent manner, but more results are required.A fundamental question is whether or not the upper and lower mantles are mixing. Traditional models that involve chemically separate upper and lower mantles cannot yet be ruled out despite recent conflicting seismological evidence showing that subducting slabs penetrate deep into the lower mantle and that chemically distinct layers are, therefore, unlikely.Recent seismic tomography studies giving three-dimensional models of the seismic wave velocities in the Earth also base their interpretations on the thermodynamic properties of minerals. These studies reveal heterogeneous velocity and density anomalies in the lower mantle, which are difficult to reconcile with mineral physics data.

Probabilistic Tomography Maps Chemical Heterogeneities Throughout the Lower Mantle

We obtained likelihoods in the lower mantle for long-wavelength models of bulk sound and shear wave speed, density, and boundary topography, compatible with gravity constraints, from normal mode splitting functions and surface wave data. Taking into account the large uncertainties in Earth's thermodynamic reference state and the published range of mineral physics data, we converted the tomographic likelihoods into probability density functions for temperature, perovskite, and iron variations. Temperature and composition can be separated, showing that chemical variations contribute to the overall buoyancy and are dominant in the lower 1000 kilometers of the mantle.

Geodynamic and seismic constraints on the thermochemical structure and dynamics of convection in the deep mantle

We revisit a recent study by Forte & Mitrovica in which global geophysical observables associated with mantle convection were inverted and the existence of a strong increase in viscosity near a depth of 2000 km was inferred. Employing mineral-physics data and theory we also showed that, although there are chemical anomalies in the lowermost mantle, they are unable to inhibit the dominant thermal buoyancy of the deep-mantle mega-plumes below Africa and the Pacific Ocean. New Monte Carlo simulations are employed to explore the impact of uncertainties in current mineral-physics constraints on inferences of deep-mantle thermochemical structure. To explore the impact of the high-viscosity peak at a depth of 2000 km on the evolution of lower-mantle structure, we carried out time-dependent convection simulations. The latter show that the stability and longevity of the dominant long-wavelength heterogeneity in the lowermost mantle are controlled by this viscosity peak.

Computational mineral physics and the physical properties of perovskite

The inherent uncertainties in modern first-principles calculations are reviewed using geophysically relevant examples. The elastic constants of perovskite at lower-mantle temperatures and pressures are calculated using ab initio molecular dynamics. These are used in conjunction with seismic tomographic models to estimate that the lateral temperature contrasts in the Earth's lower mantle are 800 K at a depth of 1000 km, and 1500 K at a depth of 2000 km. The effect of Al(3+) on the compressibility of MgSiO(3) perovskite is calculated using three different pseudopotentials. The results confirm earlier work and show that the compressibility of perovskites with Al(3+) substituted for both Si(4+) and Mg(2+) is very similar to the compressibility of Al(3+)-free perovskite. Even when 100% of the Si(4+) and Mg(2+) ions are replaced with Al(3+), the bulk modulus is only 7% less than that for Al(3+)-free perovskite. In contrast, perovskites where Al(3+) substitutes for Si(4+) only and that are charge balanced by oxygen vacancies do show higher compressibilities. When corrected to similar concentrations of Al(3+), the calculated compressibilities of the oxygen-vacancy-rich perovskites are in agreement with experimental results.

Zoned mantle convection

We review the present state of our understanding of mantle convection with respect to geochemical and geophysical evidence and we suggest a model for mantle convection and its evolution over the Earth's history that can reconcile this evidence. Whole-mantle convection, even with material segregated within the D" region just above the core-mantle boundary, is incompatible with the budget of argon and helium and with the inventory of heat sources required by the thermal evolution of the Earth. We show that the deep-mantle composition in lithophilic incompatible elements is inconsistent with the storage of old plates of ordinary oceanic lithosphere, i.e. with the concept of a plate graveyard. Isotopic inventories indicate that the deep-mantle composition is not correctly accounted for by continental debris, primitive material or subducted slabs containing normal oceanic crust. Seismological observations have begun to hint at compositional heterogeneity in the bottom 1000 km or so of the mantle, but there is no compelling evidence in support of an interface between deep and shallow mantle at mid-depth. We suggest that in a system of thermochemical convection, lithospheric plates subduct to a depth that depends - in a complicated fashion - on their composition and thermal structure. The thermal structure of the sinking plates is primarily determined by the direction and rate of convergence, the age of the lithosphere at the trench, the sinking rate and the variation of these parameters over time (i.e. plate-tectonic history) and is not the same for all subduction systems. The sinking rate in the mantle is determined by a combination of thermal (negative) and compositional buoyancy and as regards the latter we consider in particular the effect of the loading of plates with basaltic plateaux produced by plume heads. Barren oceanic plates are relatively buoyant and may be recycled preferentially in the shallow mantle. Oceanic plateau-laden plates have a more pronounced negative buoyancy and can more easily founder to the very base of the mantle. Plateau segregation remains statistical and no sharp compositional interface is expected from the multiple fate of the plates. We show that the variable depth subduction of heavily laden plates can prevent full vertical mixing and preserve a vertical concentration gradient in the mantle. In addition, it can account for the preservation of scattered remnants of primitive material in the deep mantle and therefore for the Ar and (3)He observations in ocean-island basalts.

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.

The elastic constants of MgSiO3 perovskite at pressures and temperatures of the Earth's mantle

The temperature anomalies in the Earth's mantle associated with thermal convection can be inferred from seismic tomography, provided that the elastic properties of mantle minerals are known as a function of temperature at mantle pressures. At present, however, such information is difficult to obtain directly through laboratory experiments. We have therefore taken advantage of recent advances in computer technology, and have performed finite-temperature ab initio molecular dynamics simulations of the elastic properties of MgSiO3 perovskite, the major mineral of the lower mantle, at relevant thermodynamic conditions. When combined with the results from tomographic images of the mantle, our results indicate that the lower mantle is either significantly anelastic or compositionally heterogeneous on large scales. We found the temperature contrast between the coldest and hottest regions of the mantle, at a given depth, to be about 800 K at 1,000 km, 1,500 K at 2,000 km, and possibly over 2,000 K at the core-mantle boundary.

Deep-mantle high-viscosity flow and thermochemical structure inferred from seismic and geodynamic data

Surface geophysical data that are related to the process of thermal convection in the Earth's mantle provide constraints on the rheological properties and density structure of the mantle. We show that these convection-related data imply the existence of a region of very high effective viscosity near 2,000 km depth. This inference is obtained using a viscous-flow model based on recent high-resolution seismic models of three-dimensional structure in the mantle. The high-viscosity layer near 2,000 km depth results in a re-organization of flow from short to long horizontal length scales, which agrees with seismic tomographic observations of very long wavelength structures in the deep mantle. The high-viscosity region also strongly suppresses flow-induced deformation and convective mixing in the deep mantle. Here we predict compositional and thermal heterogeneity in this region, using viscous-flow calculations based on the new viscosity profile, together with independent mineral physics data. These maps are consistent with the anti-correlation of anomalies in seismic shear and bulk sound velocity in the deep mantle. The maps also show that mega-plumes in the lower mantle below the central Pacific and Africa are, despite the presence of compositional heterogeneity, buoyant and actively upwelling structures.

Mantle convection and plate tectonics: toward an integrated physical and chemical theory

Plate tectonics and convection of the solid, rocky mantle are responsible for transporting heat out of Earth. However, the physics of plate tectonics is poorly understood; other planets do not exhibit it. Recent seismic evidence for convection and mixing throughout the mantle seems at odds with the chemical composition of erupted magmas requiring the presence of several chemically distinct reservoirs within the mantle. There has been rapid progress on these two problems, with the emergence of the first self-consistent models of plate tectonics and mantle convection, along with new geochemical models that may be consistent with seismic and dynamical constraints on mantle structure.

Compositional stratification in the deep mantle

A boundary between compositionally distinct regions at a depth of about 1600 kilometers may explain the seismological observations pertaining to Earth's lower mantle, produce the isotopic signatures of mid-ocean ridge basalts and oceanic island basalts, and reconcile the discrepancy between the observed heat flux and the heat production of the mid-ocean ridge basalt source region. Numerical models of thermochemical convection imply that a layer of material that is intrinsically about 4 percent more dense than the overlying mantle is dynamically stable. Because the deep layer is hot, its net density is only slightly greater than adiabatic and its surface develops substantial topography.

Compositional heterogeneity in the bottom 1000 kilometers of Earth's mantle: toward a hybrid convection model

Tomographic imaging indicates that slabs of subducted lithosphere can sink deep into Earth's lower mantle. The view that convective flow is stratified at 660-kilometer depth and preserves a relatively pristine lower mantle is therefore not tenable. However, a range of geophysical evidence indicates that compositionally distinct, hence convectively isolated, mantle domains may exist in the bottom 1000 kilometers of the mantle. Survival of these domains, which are perhaps related to local iron enrichment and silicate-to-oxide transformations, implies that mantle convection is more complex than envisaged by conventional end-member flow models.