Allen Keith McNamara

Thermochemical structures beneath Africa and the Pacific Ocean

Large low-velocity seismic anomalies have been detected in the Earths lower mantle beneath Africa and the Pacific Ocean that are not easily explained by temperature variations alone. The African anomaly has been interpreted to be a northwest–southeast-trending structure with a sharp-edged linear, ridge-like morphology. The Pacific anomaly, on the other hand, appears to be more rounded in shape. Mantle models with heterogeneous composition have related these structures to dense thermochemical piles or superplumes. It has not been shown, however, that such models can lead to thermochemical structures that satisfy the geometrical constraints, as inferred from seismological observations. Here we present numerical models of thermochemical convection in a three-dimensional spherical geometry using plate velocities inferred for the past 119 million years. We show that Earths subduction history can lead to thermochemical structures similar in shape to the observed large, lower-mantle velocity anomalies. We find that subduction history tends to focus dense material into a ridge-like pile beneath Africa and a relatively more-rounded pile under the Pacific Ocean, consistent with seismic observations.

Structure and Dynamics of Earth's Lower Mantle

Processes within the lowest several hundred kilometers of Earths rocky mantle play a critical role in the evolution of the planet. Understanding Earths lower mantle requires putting recent seismic and mineral physics discoveries into a self-consistent, geodynamically feasible context. Two nearly antipodal large low-shear-velocity provinces in the deep mantle likely represent chemically distinct and denser material. High-resolution seismological studies have revealed laterally varying seismic velocity discontinuities in the deepest few hundred kilometers, consistent with a phase transition from perovskite to post-perovskite. In the deepest tens of kilometers of the mantle, isolated pockets of ultralow seismic velocities may denote Earths deepest magma chamber.

Deformation of (Mg,Fe)SiO3 Post-Perovskite and D'' Anisotropy

Polycrystalline (Mg0.9,Fe0.1)SiO3 post-perovskite was plastically deformed in the diamond anvil cell between 145 and 157 gigapascals. The lattice-preferred orientations obtained in the sample suggest that slip on planes near (100) and (110) dominate plastic deformation under these conditions. Assuming similar behavior at lower mantle conditions, we simulated plastic strains and the contribution of post-perovskite to anisotropy in the D″ region at the Earth core-mantle boundary using numerical convection and viscoplastic polycrystal plasticity models. We find a significant depth dependence of the anisotropy that only develops near and beyond the turning point of a downwelling slab. Our calculated anisotropies are strongly dependent on the choice of elastic moduli and remain hard to reconcile with seismic observations.

Development of anisotropic structure in the Earth's lower mantle by solid-state convection

Seismological observations reveal highly anisotropic patches at the bottom of the Earths lower mantle, whereas the bulk of the mantle has been observed to be largely isotropic. These patches have been interpreted to correspond to areas where subduction has taken place in the past or to areas where mantle plumes are upwelling, but the underlying cause for the anisotropy is unknown—both shape-preferred orientation of elastically heterogenous materials and lattice-preferred orientation of a homogeneous material have been proposed. Both of these mechanisms imply that large-strain deformation occurs within the anisotropic regions, but the geodynamic implications of the mechanisms differ. Shape-preferred orientation would imply the presence of large elastic (and hence chemical) heterogeneity whereas lattice-preferred orientation requires deformation at high stresses. Here we show, on the basis of numerical modelling incorporating mineral physics of elasticity and development of lattice-preferred orientation, that slab deformation in the deep lower mantle can account for the presence of strong anisotropy in the circum-Pacific region. In this model—where development of the mineral fabric (the alignment of mineral grains) is caused solely by solid-state deformation of chemically homogeneous mantle material—anisotropy is caused by large-strain deformation at high stresses, due to the collision of subducted slabs with the core–mantle boundary.

Tomographic filtering of geodynamic models: Implications for model interpretation and large-scale mantle structure

[1] The resolution operator R is a critical accompaniment to tomographic models of the mantle. R facilitates the comparison between conceptual three-dimensional velocity models and tomographic models because it can filter these theoretical models to the spatial resolution of the tomographic model. We compute R for the tomographic model S20RTS (Ritsema et al., 1999, 2004) and two companion models that are based on the same data but derived with different norm damping values. The three models explain (within measurement uncertainty) S-SKS and S-SKKS travel times equally well. To demonstrate how artifacts distort tomographic images and complicate model interpretation, we apply R to (1) a thermochemical and (2) an isochemical model of convection in the mantle that feature different patterns of shear velocity heterogeneity in the deep mantle if we assume that shear velocity heterogeneity is caused by temperature variations only. R suppresses short-wavelength structures, removes strong velocity gradients, and introduces artificial stretching and tilting of velocity anomalies. Temperature anomalies in the thermochemical model resemble the spatial extent of low seismic velocity anomalies and the shear velocity spectrum in the D’’ region better than the isochemical model. However, the thermochemical model overpredicts the amplitude of shear velocity variation and places the African and Pacific anomalies imperfectly. We suspect that inaccurate velocity scaling laws and uncertain initial conditions control these mismatches. Extensive hypothesis testing is required to identify successful models.

Localization of dislocation creep in the lower mantle: implications for the origin of seismic anisotropy

Recent seismological observations reveal the presence of seismic anisotropy in localized regions at the base of the mantle within an otherwise isotropic lower mantle. These regions can be placed in a tectonic context, corresponding to locations of paleosubduction and plume upwelling. This project works toward determining whether the observed seismic anisotropy may be explained by the development of a mineral fabric by lattice-preferred orientation (LPO). Numerical modeling is used to explore whether the conditions at the base of upwelling and downwelling regions are consistent with those required for fabric development. Specifically, we examine whether dislocation creep dominates these regions within a background mantle that flows primarily by diffusion creep. The key to our study is the use of a composite rheology that includes both mechanisms of diffusion and dislocation creep and is based on mineral physics experiments. Results show that it is possible to produce a localization of dislocation creep near slabs within a background mantle dominated by diffusion creep. In contrast, upwelling regions are characterized by a domination of diffusion creep. These results indicate that LPO may be the cause of lowermost mantle seismic anisotropy near paleoslabs, but other mechanisms such as shape-preferred orientation may be required to produce the anisotropy observed near upwellings. fl 2001 Elsevier Science B.V. All rights reserved.

West African proximity of the Avalon terrane in the latest Precambrian

Considerable debate surrounds the Late Neoproterozoic paleogeographic position of the Avalon terrane, specifically whether it was adjacent to West Africa or Amazonia. New paleomagnetic results from upper Neoproterozoic rocks in the Avalon terrane challenge the latter position. Samples collected from the ca. 580‐570 Ma Marystown Group in the southern part of the Burin peninsula of Newfoundland, Canada, yield high-temperature magnetic components, including dual-polarity directions, which are considered to be primary, on the basis of positive fold and agglomerate tests. The resultant tilt-corrected inclination is 538, representing a paleolatitude of deposition of 34 81 88/‐78 for the Marystown Group. Given the likelihood that Amazonia and Laurentia were still juxtaposed around 580‐570 Ma, the Marystown Group results reveal that the paleolatitude of Avalon is significantly lower than would be expected if it was part of Amazonia. In fact, Avalon was separated from northern Amazonia by at least 1100 km at ca. 580 Ma. If West Africa was juxtaposed to Amazonia, opposite Laurentia, by this time, these results place Avalon at the same paleolatitude as the northern margin of the West African craton.

Episodic entrainment of deep primordial mantle material into ocean island basalts

Chemical differences between mid-ocean ridge basalts (MORBs) and ocean island basalts (OIBs) provide critical evidence that the Earths mantle is compositionally heterogeneous. MORBs generally exhibit a relatively low and narrow range of 3He/4He ratios on a global scale, whereas OIBs display larger variability in both time and space. The primordial origin of 3He in OIBs has motivated hypotheses that high 3He/4He ratios are the product of mantle plumes sampling chemically distinct material, but do not account for lower MORB-like 3He/4He ratios in OIBs, nor their observed spatial and temporal variability. Here we perform thermochemical convection calculations which show the variable 3He/4He signature of OIBs can be reproduced by deep isolated mantle reservoirs of primordial material that are viscously entrained by thermal plumes. Entrainment is highly time-dependent, producing a wide range of 3He/4He ratios similar to that observed in OIBs worldwide and indicate MORB-like 3He/4He ratios in OIBs cannot be used to preclude deep mantle-sourced hotspots.

The Core–Mantle Boundary Region

The Earths core–mantle boundary (CMB) is the largest density discontinuity in our planet, separating the rocky oxide mantle from the liquid metal core. The structure and dynamics of the CMB region are central to major questions such as the depth extent of lithospheric subduction, the thermal and chemical evolution of the Earth, the mechanisms and energetics involved in the generation of a geomagnetic field, and the nucleation of deep-seated mantle plumes that rise up to trace out volcanic hot spot tracks at the surface. Owing to the proliferation of digital seismic networks and developments in high-pressure experimental and theoretical mineral physics, we are now learning more about the structure and physical properties of the CMB region and the myriad ways in which deep processes are linked to the surface environment. This chapter surveys recent developments in understanding the nature and dynamics of this critical region of the Earths interior.

Fine-Scale Ultra-Low Velocity Zone Layering at the Core-Mantle Boundary and Superplumes

Ultra-low velocity layering at the Earth’s core-mantle boundary (CMB) has now been detected using a variety of seismic probes. P- and S-wave velocity reductions of up to 10’s of percent have been mapped in a thin (5–50 km) layer, which commonly underlies reduced seismic shear wave speeds in the overlying few 100 km of the mantle. Ultra-low velocity zones (ULVZ) contain properties consistent with partial melt of rock at the very base of the mantle. Strong evidence now exists for a significant density increase in the layer (∼5–10% greater than reference models), which must be included in dynamical scenarios relating ULVZ partial melt to deep mantle plume genesis. 3-D geodynamical calculations involving an initially uniform dense layer in the lowermost few 100 km of the mantle result in thermo-chemical piles that are geographically well-correlated with seismic tomography low velocities, when past plate motions are imposed as a surface boundary condition. The hottest lower mantle regions underlay edges of the dense thermo-chemical piles. A scenario is put forth where these piles geographically correlate with ultra-low velocity zones, and subsequent mantle plume genesis.