The Dynamic Structure of the Deep Earth: An Interdisciplinary Approach

Peer Reviewedhttp://deepblue.lib.umich.edu/bitstream/2027.42/94755/1/eost14654.pd

Thermodynamics of mantle minerals – I. Physical properties

We present a theory for the computation of phase equilibria and physical properties of multicomponent assemblages relevant to the mantle of the Earth. The theory differs from previous treatments in being thermodynamically self-consistent: the theory is based on the concept of fundamental thermodynamic relations appropriately generalized to anisotropic strain and in encompassing elasticity in addition to the usual isotropic thermodynamic properties. In this first paper, we present the development of the theory, discuss its scope, and focus on its application to physical properties of mantle phases at elevated pressure and temperature including the equation of state, thermochemical properties and the elastic wave velocities. We find that the Eulerian finite strain formulation captures the variation of the elastic moduli with compression. The variation of the vibrational frequencies with compression is also cast as a Taylor series expansion in the Eulerian finite strain, the appropriate volume derivative of which leads to an expression for the GrÜneisen parameter that agrees well with results from first principles theory. For isotropic materials, the theory contains nine material-specific parameters: the values at ambient conditions of the Helmholtz free energy, volume, bulk and shear moduli, their pressure derivatives, an effective Debye temperature, its first and second logarithmic volume derivatives (Γ 0 , q 0 ) , and the shear strain derivative of Γ. We present and discuss in some detail the results of a global inversion of a wide variety of experimental data and first principles theoretical results, supplemented by systematic relations, for the values of these parameters for 31 mantle species. Among our findings is that the value of q is likely to be significantly greater than unity for most mantle species. We apply the theory to the computation of the shear wave velocity, and temperature and compositional (Fe content) derivatives at relevant mantle pressure temperature conditions. Among the patterns that emerge is that garnet is anomalous in being remarkably insensitive to iron content or temperature as compared with other mantle phases.Peer Reviewedhttp://deepblue.lib.umich.edu/bitstream/2027.42/73485/1/j.1365-246X.2005.02642.x.pd

Why is Africa rifting?

Continental rifting has a fundamental role in the tectonic behaviour of the Earth, shaping the surface we live on. Although there is not yet a consensus about the dominant mechanism for rifting, there is a general agreement that the stresses required to rift the continental lithosphere are not readily available. Here we use a global finite element model of the lithosphere to calculate the stresses acting on Africa. We consider the stresses induced by mantle flow, crustal structure and topography in two types of models: one in which flow is exclusively driven by the subducting slabs and one in which it is derived from a shear wave tomographic model. The latter predicts much larger stresses and a more realistic dynamic topography. It is therefore clear that the mantle structure beneath Africa plays a key part in providing the radial and horizontal tractions, dynamic topography and gravitational potential energy necessary for rifting. Nevertheless, the total available stress (c. 100 MPa) is much less than that needed to break thick, cold continental lithosphere. Instead, we appeal to a model of magma-assisted rifting along pre-existing weaknesses, where the strain is localized in a narrow axial region and the strength of the plate is reduced significantly. Mounting geological and geophysical observations support such a model

Influence of continental roots and asthenosphere on plate‐mantle coupling

Peer Reviewedhttp://deepblue.lib.umich.edu/bitstream/2027.42/94652/1/grl21239.pd

Inferring the thermochemical structure of the upper mantle from seismic data

We test a mineral physics model of the upper mantle against seismic observations. The model is based on current knowledge of material properties at high temperatures and pressures. In particular, elastic properties are computed with a recent self-consistent thermodynamic model, based on a six oxides (NCFMAS) system. We focus on average structure between 250 and 800 km. We invert normal modes eigenfrequencies and traveltimes to obtain best-fitting average thermal structures for various compositional profiles. The thermochemical structures are then used to predict long-period waveforms, SS precursors waveforms and radial profiles of attenuation. These examples show the potential of our procedure to refine the interpretation combining different data sets. We found that a mixture of MORB and Harzburgite, with the MORB component increasing with depth, is able to reproduce well all the seismic data for realistic thermal structures. If the proportions of MORB with depth do not change, unrealistic negative thermal gradients below 250 km would be necessary to explain the data. Equilibrium assemblages, such as pyrolite, cannot fit the seismic data. The elastic velocities predicted by the reference mineral physics model tested are too low at the top of the lower mantle, even for the fastest (and most depleted) composition, that is, harzburgite. An increase in VP of 1 per cent and in VS of 2 per cent improves the data fit significantly and is required to find models that fit both traveltimes and normal modes, indicating the need for further experimental measurements of these properties at the simultaneously elevated pressure—temperature conditions of the lower mantle. Extending our procedure to other seismic and density data and interpreting the 3-D structure holds promise to further improve our knowledge of the thermochemical structure of the upper mantle. In addition, the same database of material properties can be used in dynamic models to test whether the thermochemical structure inferred from geophysical observations is consistent with the Earth's evolutio

Inferring the thermochemical structure of the upper mantle from seismic data

We test a mineral physics model of the upper mantle against seismic observations. The model is based on current knowledge of material properties at high temperatures and pressures. In particular, elastic properties are computed with a recent self-consistent thermodynamic model, based on a six oxides (NCFMAS) system. We focus on average structure between 250 and 800 km. We invert normal modes eigenfrequencies and traveltimes to obtain best-fitting average thermal structures for various compositional profiles. The thermochemical structures are then used to predict long-period waveforms, SS precursors waveforms and radial profiles of attenuation. These examples show the potential of our procedure to refine the interpretation combining different data sets.We found that a mixture of MORB and Harzburgite, with the MORB component increasing with depth, is able to reproduce well all the seismic data for realistic thermal structures. If the proportions of MORB with depth do not change, unrealistic negative thermal gradients below 250 km would be necessary to explain the data. Equilibrium assemblages, such as pyrolite, cannot fit the seismic data.The elastic velocities predicted by the reference mineral physics model tested are too low at the top of the lower mantle, even for the fastest (and most depleted) composition, that is, harzburgite. An increase in V P of 1 per cent and in V S of 2 per cent improves the data fit significantly and is required to find models that fit both traveltimes and normal modes, indicating the need for further experimental measurements of these properties at the simultaneously elevated pressure–temperature conditions of the lower mantle.Extending our procedure to other seismic and density data and interpreting the 3-D structure holds promise to further improve our knowledge of the thermochemical structure of the upper mantle. In addition, the same database of material properties can be used in dynamic models to test whether the thermochemical structure inferred from geophysical observations is consistent with the Earth's evolution.Peer Reviewedhttp://deepblue.lib.umich.edu/bitstream/2027.42/78699/1/j.1365-246X.2009.04338.x.pd

Iceland, the Farallon slab, and dynamic topography of the North Atlantic

ABSTRACT Upwelling or downwelling flow in Earth's mantle is thought to elevate or depress Earth's surface on a continental scale. Direct observation of this ''dynamic topography'' on the seafloor, however, has remained elusive because it is obscured by isostatically supported topography caused by near-surface density variations. We calculate the nonisostatic topography of the North Atlantic by correcting seafloor depths for lithospheric cooling and sediment loading, and find that seafloor west of the Mid-Atlantic Ridge is an average of 0.5 km deeper than it is to the east. We are able to reproduce this basic observation in a model of mantle flow driven by tomographically inferred mantle densities. This model shows that the Farallon slab, currently in the lower mantle beneath the east coast of North America, induces downwelling flow that deepens the western North Atlantic relative to the east. Our model also predicts dynamic support of observed topographic highs near Iceland and the Azores, but suggests that the Icelandic high is due to local upper-mantle upwelling, while the Azores high is part of a plate-scale lower-mantle upwelling to the south. An anomalously deep area off the coast of Nova Scotia may be associated with the downwelling component of edge-driven convection at the continental boundary. Thus, several of the seafloor's topographic features can only be understood in terms of dynamic support from flow in Earth's mantle

Extrinsic elastic anisotropy in a compositionally heterogeneous Earth's mantle

Several theoretical studies indicate that a substantial fraction of the measured seismic anisotropy could be interpreted as extrinsic anisotropy associated with compositional layering in rocks, reducing the significance of strain‐induced intrinsic anisotropy. Here, we quantify the potential contribution of grain‐scale and rock‐scale compositional anisotropy to the observations by (i) combining effective medium theories with realistic estimates of mineral isotropic elastic properties, and (ii) measuring velocities of synthetic seismic waves propagating through modelled strain‐induced microstructures. It is shown that for typical mantle and oceanic crust sub‐solidus compositions, rock‐scale compositional layering does not generate any substantial extrinsic anisotropy (<1%) because of the limited contrast in isotropic elastic moduli among different rocks. Quasi‐laminated structures observed in subducting slabs using P‐ and S‐ wave scattering are often invoked as a source of extrinsic anisotropy, but our calculations show that they only generate minor seismic anisotropy (<0.1‐0.2% of Vp and Vs radial anisotropy). More generally, rock‐scale compositional layering, when present, cannot be detected with seismic anisotropy studies, but mainly with wave scattering. In contrast, when grain‐scale layering is present, significant extrinsic anisotropy could exist in vertically limited levels of the mantle such as in a MORB‐rich lower transition zone or in the uppermost lower mantle where foliated basalts and pyrolites display up to 2‐3% Vp and 3‐6% Vs radial anisotropy. Thus, seismic anisotropy observed around the 660 km discontinuity could be possibly related to grain‐scale SPO. Extrinsic anisotropy can form also in a compositionally homogeneous mantle, where velocity variations associated with major phase transitions can generate up to 1% of positive radial anisotropy

Temperature and velocity measurements of a rising thermal plume

Author Posting. © American Geophysical Union, 2015. This article is posted here by permission of American Geophysical Union for personal use, not for redistribution. The definitive version was published in Geochemistry, Geophysics, Geosystems 16 (2015): 579–599, doi:10.1002/2014GC005576.The three-dimensional velocity and temperature fields surrounding an isolated thermal plume in a fluid with temperature-dependent viscosity are measured using Particle-Image Velocimetry and thermochromatic liquid crystals, respectively. The experimental conditions are relevant to a plume rising through the mantle. It is shown that while the velocity and the isotherm surrounding the plume can be used to visualize the plume, they do not reveal the finer details of its structure. However, by computing the Finite-Time Lyapunov Exponent fields from the velocity measurements, the material lines of the flow can be found, which clearly identify the shape of the plume head and characterize the behavior of the flow along the plume stem. It is shown that the vast majority of the material in the plume head has undergone significant stretching and originates from a wide region very low in the fluid domain, which is proposed as a contributing factor to the small-scale isotopic variability observed in ocean-island basalt regions. Lastly, the Finite-Time Lyapunov Exponent fields are used to calculate the steady state rise velocity of the thermal plume, which is found to scale linearly with the Rayleigh number, in contrast to some previous work. The possible cause and the significance of these conflicting results are discussed, and it is suggested that the scaling relationship may be affected by the temperature-dependence of the fluid viscosity in the current work.This work was funded by the National Science Foundation (grant EAR-055199) and the MAPS Dean's Office at UCL.2015-09-0