Dan J. Bower

Lower mantle structure from paleogeographically constrained dynamic Earth models

Seismic tomography reveals two large, low-shear velocity provinces (LLSVPs) beneath Africa and the Pacific Ocean. These structures may have existed for several 100 Myr and are likely compositionally distinct based on observed seismic characteristics interpreted in light of geodynamic models and mineral physics constraints. We investigate the dynamics of the LLSVPs through the use of evolutionary models of thermochemical convection from 250 Ma to present day. We use a spherical convection model in which the anomalous structures have a high bulk modulus, consistent with seismic interpretation. A new progressive assimilation method incorporates constraints from paleogeography using a refined plate history model (with 1 Myr time spacing) to guide the thermal structure of the lithosphere and steer the thermal evolution of slabs in the uppermost mantle. The thermochemical structures deform and migrate along the core-mantle boundary (CMB) through coupling to plate motions and in response to slab stresses. The models produce a ridge-like anomaly beneath Africa and a rounded pile beneath the Pacific Ocean, which at present day agrees with tomography, waveform modeling, and other geodynamic studies. Plumes emanate from the margins of the domes and ridges of thickened boundary layer between the domes. Dense and viscous slabs can undermine the stability of high bulk modulus structures at the CMB. High bulk modulus structures are not necessarily required to satisfy dynamic constraints on the LLSVPs.

Ridge subduction sparked reorganization of the Pacific plate-mantle system 60-50 million years ago

A reorganization centered on the Pacific plate occurred ~53–47 million years ago. A “top-down” plate tectonic mechanism, complete subduction of the Izanagi plate, as opposed to a “bottom-up” mantle flow mechanism, has been proposed as the main driver. Verification based on marine geophysical observations is impossible as most ocean crust recording this event has been subducted. Using a forward modeling approach, which assimilates surface plate velocities and shallow thermal structure of slabs into mantle flow models, we show that complete Izanagi plate subduction and margin-wide slab detachment induced a major change in sub-Pacific mantle flow, from dominantly southward before 60 Ma to north-northeastward after 50 Ma. Our results agree with onshore geology, mantle tomography, and the inferred motion of the Hawaiian hot spot and are consistent with a plate tectonic process driving the rapid plate-mantle reorganization in the Pacific hemisphere between 60 and 50 Ma. This reorganization is reflected in tectonic changes in the Pacific and surrounding ocean basins.

A geodynamic and mineral physics model of a solid-state ultralow-velocity zone

article i nfo Article history: Recent results (Wicks et al., 2010) suggest that a mixture of iron-enriched (Mg,Fe)O and ambient mantle is consistent with wavespeed reductions and density increases inferred for ultralow-velocity zones (ULVZs). We explore this hypothesis by simulating convection to deduce the stability and morphology of such chemically- distinct structures. The buoyancy number, or chemical density anomaly, largely dictates ULVZ shape, and the prescribed initial thickness (proxy for volume) of the chemically-distinct layer controls its size. We synthesize our dynamic results with a Voigt-Reuss-Hill mixing model to provide insight into the inherent seismic tradeoff between ULVZ thickness and wavespeed reduction. Seismic data are compatible with a solid-state origin for ULVZs, and a suite of these structures may scatter seismic energy to produce broadband PKP precursors.

Enhanced convection and fast plumes in the lower mantle induced by the spin transition in ferropericlase

Using a numerical model we explore the consequences of the intrinsic density change (Δρ/ρ ≈ 2–4%) caused by the Fe^(2+) spin transition in ferropericlase on the style and vigor of mantle convection. The effective Clapeyron slope of the transition from high to low spin is strongly positive in pressure-temperature space and broadens with high temperature. This introduces a net spin-state driving density difference for both upwellings and downwellings. In 2-D cylindrical geometry spin-buoyancy dominantly enhances the positive thermal buoyancy of plumes. Although the additional buoyancy does not fundamentally alter large-scale dynamics, the Nusselt number increases by 5–10%, and vertical velocities by 10–40% in the lower mantle. Advective heat transport is more effective and temperatures in the core-mantle boundary region are reduced by up to 12%. Our findings are relevant to the stability of lowermost mantle structures.

Direct measures of lateral velocity variation in the deep Earth

Current tomographic models of the Earth display perturbations to a radial stratified reference model. However, structures in the deep mantle that are chemically dense with low Rayleigh numbers can develop enormous relief, perhaps with boundaries closer to vertical than to radial. Such features are hard to detect with present tomographic modeling techniques because the timing anomalies are based on long-period filtered waveforms with complexity removed. Here we develop a new tool for processing array data on the basis of a decomposition referred to as a multipath detector, which can be used to distinguish between horizontal structure (in-plane multipathing) and vertical (out-of-plane multipathing) directly from processing array waveforms. A lateral gradient coefficient based on this detector provides a direct constraint on the sharpness of the boundaries and material properties. We demonstrate the usefulness of this approach by processing samples of both P and S data from the Kaapvaal array in southern Africa, which are compared with synthetic predictions from a metastable dynamic model containing sharp edges. Both data and simulations produce timing gradients larger than 2 s/deg in azimuthal changes for S waves, where only minor effects are obtained for P waves. These results further validate the case for distinct chemistry inside the African Low Shear Velocity Province.

Provenance of plumes in global convection models

In global convection models constrained by plume motions and subduction history over the last 230 Myr, plumes emerge preferentially from the edges of thermochemical structures that resemble present-day large low shear velocity provinces (LLSVPs) beneath Africa and the Pacific Ocean. It has been argued that large igneous provinces (LIPs) erupting since 200 Ma may originate from plumes that emerged from the edges of the LLSVPs and numerical models have been devised to validate this hypothesis. Although qualitative assessments that are broadly in agreement with this hypothesis have been derived from numerical models, a quantitative assessment has been lacking. We present a novel plume detection scheme and derive Monte Carlo-based statistical correlations of model plume eruption sites and reconstructed LIP eruption sites. We show that models with a chemically anomalous lower mantle are highly correlated to reconstructed LIP eruption sites, whereas the confidence level obtained for a model with purely thermal plumes falls just short of 95%. A network of embayments separated by steep ridges form in the deep lower mantle in models with a chemically anomalous lower mantle. Plumes become anchored to the peaks of the chemical ridges and the network of ridges acts as a floating anchor, adjusting to slab push forces through time. The network of ridges imposes a characteristic separation between conduits that can extend into the interior of the thermochemical structures. This may explain the observed clustering of reconstructed LIP eruption sites that mostly but not exclusively occur around the present-day LLSVPs.

Origin and evolution of the deep thermochemical structure beneath Eurasia

A unique structure in the Earths lowermost mantle, the Perm Anomaly, was recently identified beneath Eurasia. It seismologically resembles the large low-shear velocity provinces (LLSVPs) under Africa and the Pacific, but is much smaller. This challenges the current understanding of the evolution of the plate–mantle system in which plumes rise from the edges of the two LLSVPs, spatially fixed in time. New models of mantle flow over the last 230 million years reproduce the present-day structure of the lower mantle, and show a Perm-like anomaly. The anomaly formed in isolation within a closed subduction network ∼22,000 km in circumference prior to 150 million years ago before migrating ∼1,500 km westward at an average rate of 1 cm year−1, indicating a greater mobility of deep mantle structures than previously recognized. We hypothesize that the mobile Perm Anomaly could be linked to the Emeishan volcanics, in contrast to the previously proposed Siberian Traps.