No abstract

Grain size plays a key role in controlling the mechanical properties of the Earth's mantle, affecting both long‐time‐scale flow patterns and anelasticity on the time scales of seismic wave propagation. However, dynamic models of Earth's convecting mantle usually implement flow laws with constant grain size, stress‐independent viscosity, and a limited treatment of changes in mineral assemblage. We study grain size evolution, its interplay with stress and strain rate in the convecting mantle, and its influence on seismic velocities and attenuation. Our geodynamic models include the simultaneous and competing effects of dynamic recrystallization resulting from dislocation creep, grain growth in multiphase assemblages, and recrystallization at phase transitions. They show that grain size evolution drastically affects the dynamics of mantle convection and the rheology of the mantle, leading to lateral viscosity variations of 6 orders of magnitude due to grain size alone, and controlling the shape of upwellings and downwellings. Using laboratory‐derived scaling relationships, we convert model output to seismologically observable parameters (velocity and attenuation) facilitating comparison to Earth structure. Reproducing the fundamental features of the Earth's attenuation profile requires reduced activation volume and relaxed shear moduli in the lower mantle compared to the upper mantle, in agreement with geodynamic constraints. Faster lower mantle grain growth yields best fit to seismic observations, consistent with our reexamination of high‐pressure grain growth parameters. We also show that ignoring grain size in interpretations of seismic anomalies may underestimate the Earth's true temperature variations.

A large data set of surface wave phase anomalies, body wave travel times, normal‐mode splitting functions, and long‐period waveforms is used to investigate the scaling between shear velocity, density, and compressional velocity in the Earth's mantle. We introduce a methodology that allows construction of joint models with various levels of scaling complexity (ϱ = d lnρ/d lnvS, ν = d lnvS/d lnvP), in order to detect seismological signatures of chemical heterogeneity. We demonstrate that the data sets considered cannot be fit concurrently with a uniform ν or a positive and uniform ϱ throughout the mantle. The variance reductions to P wave travel times and vP‐sensitive modes are up to 40% higher with our preferred model of anisotropic shear and compressional velocity than the recent anisotropic shear velocity model S362ANI+M, which was constructed assuming a uniform ν throughout the mantle. Several features reported in earlier tomographic studies persist after the inclusion of new and larger data sets; anticorrelation between bulk sound and shear velocities in the lowermost mantle as well as an increase in ν with depth in the lower mantle are largely independent of the regularization scheme. When correlations between density and shear velocity variations are imposed in the lowermost mantle, variance reductions of several spheroidal and toroidal modes deteriorate by as much as 40%. Recent measurements of the splitting of 0S2, in particular, are largely incompatible with perfectly correlated shear velocity and density heterogeneity throughout the mantle. A way to significantly improve the fits to various data sets is by allowing independent density perturbations in the lowermost mantle. Our preferred joint model consists of denser‐than‐average anomalies (∼1% peak to peak) at the base of the mantle roughly coincident with the low‐velocity superplumes. The relative variation of shear velocity, density, and compressional velocity in our study disfavors a purely thermal contribution to heterogeneity in the lowermost mantle, with implications for the long‐term stability and evolution of superplumes.