James Wookey

First-principles constraints on diffusion in lower-mantle minerals and a weak D′′ layer

Post-perovskite MgSiO3 is believed to be present in the D′′ region of the Earth’s lowermost mantle. Its existence has been used to explain a number of seismic observations, such as the D′′ reflector and the high degree of seismic anisotropy within the D′′ layer. Ionic diffusion in post-perovskite controls its viscosity, which in turn controls the thermal and chemical coupling between the core and the mantle, the development of plumes and the stability of deep chemical reservoirs. Here we report the use of first-principles methods to calculate absolute diffusion rates in post-perovskite under the conditions found in the Earth’s lower mantle. We find that the diffusion of Mg2+ and Si4+ in post-perovskite is extremely anisotropic, with almost eight orders of magnitude difference between the fast and slow directions. If post-perovskite in the D′′ layer shows significant lattice-preferred orientation, the fast diffusion direction will render post-perovskite up to four orders of magnitude weaker than perovskite. The presence of weak post-perovskite strongly increases the heat flux across the core–mantle boundary and alters the geotherm. It also provides an explanation for laterally varying viscosity in the lowermost mantle, as required by long-period geoid models. Moreover, the behaviour of very weak post-perovskite can reconcile seismic observation of a D′′ reflector with recent experiments showing that the width of the perovskite-to-post-perovskite transition is too wide to cause sharp reflectors. We suggest that the observed sharp D′′ reflector is caused by a rapid change in seismic anisotropy. Once sufficient perovskite has transformed into post-perovskite, post-perovskite becomes interconnected and strain is partitioned into this weaker phase. At this point, the weaker post-perovskite will start to deform rapidly, thereby developing a strong crystallographic texture. We show that the expected seismic contrast between the deformed perovskite-plus-post-perovskite assemblage and the overlying isotropic perovskite-plus-post-perovskite assemblage is consistent with seismic observations.

Mid-mantle deformation inferred from seismic anisotropy

With time, convective processes in the Earths mantle will tend to align crystals, grains and inclusions. This mantle fabric is detectable seismologically, as it produces an anisotropy in material properties—in particular, a directional dependence in seismic-wave velocity. This alignment is enhanced at the boundaries of the mantle where there are rapid changes in the direction and magnitude of mantle flow, and therefore most observations of anisotropy are confined to the uppermost mantle or lithosphere and the lowermost-mantle analogue of the lithosphere, the D″ region. Here we present evidence from shear-wave splitting measurements for mid-mantle anisotropy in the vicinity of the 660-km discontinuity, the boundary between the upper and lower mantle. Deep-focus earthquakes in the Tonga–Kermadec and New Hebrides subduction zones recorded at Australian seismograph stations record some of the largest values of shear-wave splitting hitherto reported. The results suggest that, at least locally, there may exist a mid-mantle boundary layer, which could indicate the impediment of flow between the upper and lower mantle in this region.

Deformation of the Lowermost Mantle from Seismic Anisotropy

The lowermost part of the Earth’s mantle—known as D″—shows significant seismic anisotropy, the variation of seismic wave speed with direction. This is probably due to deformation-induced alignment of MgSiO3-post-perovskite (ppv), which is believed to be the main mineral phase present in the region. If this is the case, then previous measurements of D″ anisotropy, which are generally made in one direction only, are insufficient to distinguish candidate mechanisms of slip in ppv because the mineral is orthorhombic. Here we measure anisotropy in D″ beneath North and Central America, where material from subducting oceanic slabs impinges on the core–mantle boundary, using shallow as well as deep earthquakes to increase the azimuthal coverage in D″. We make more than 700 individual measurements of shear wave splitting in D″ in three regions from two different azimuths in each case. We show that the previously assumed case of vertical transverse isotropy (where wave speed shows no azimuthal variation) is not possible, and that more complicated mechanisms must be involved. We test the fit of different MgSiO3-ppv deformation mechanisms to our results and find that shear on (001) is most consistent with observations and the expected shear above the core–mantle boundary beneath subduction zones. With new models of mantle flow, or improved experimental determination of the dominant ppv slip systems, this method will allow us to map deformation at the core–mantle boundary and link processes in D″, such as plume initiation, to the rest of the mantle.

Precambrian plate tectonics: Seismic evidence from northern Hudson Bay, Canada

The Canadian Shield is one of the largest exposures of Precambrian rocks on Earth. It is a mosaic of several Archean terranes that were brought together during a series of Paleoproterozoic orogens culminating in the so-called Trans-Hudson orogen, which is thought to have been similar to the Himalayan orogen in scale and nature. The tectonic evolution and lithospheric subdivisions of this region are poorly understood, but new seismic networks in northern Hudson Bay provide fresh opportunity to place constraints on the Precambrian processes that formed and shaped it. We show, via a study of seismic anisotropy, that the lithosphere of the northern Hudson Bay region retains a strong signature of Archean–Paleoproterozoic tectonics. We map the boundary between the upper (Churchill) and lower (Superior) plates that collided ca. 1.8 Ga and identify back azimuth–dependent splitting parameters (φ, δ t ) on Baffin Island that indicate complex anisotropy (e.g., dipping fabric) beneath the region. Our results support the view that significant lithospheric deformation occurred during the Paleoproterozoic and that modern-day plate tectonic processes were thus in operation by at least ca. 1.8 Ga.

Inner-core shear-wave anisotropy and texture from an observation of PKJKP waves

Since the discovery of the Earth’s core a century ago, and the subsequent discovery of a solid inner core (postulated to have formed by the freezing of iron) seismologists have striven to understand this most remote part of the deep Earth. The most direct evidence for a solid inner core would be the observation of shear-mode body waves that traverse it, but these phases are extremely difficult to observe. Two reported observations in short-period data have proved controversial. Arguably more successful have been studies of longer-period data, but such averaging limits the usefulness of the observations to reported sightings. We present two observations of an inner-core shear-wave phase at higher frequencies in stacked data from the Japanese High-Sensitivity Array, Hi-Net. From an analysis of timing, amplitude and waveform of the ‘PKJKP’ phase we derive constraints on inner-core compressional-wave velocity and shear attenuation at about 0.3 Hz which differ from standard isotropic core models. We can explain waveform features and can partially reconcile the otherwise large differences between core wavespeed and attenuation models that our observations apparently suggest if we invoke shear-wave anisotropy in the inner core. A simple model of an inner core composed of hexagonal close-packed iron with its c axis aligned perpendicular to the rotation axis yields anisotropy that is compatible with both the shear-wave anisotropy that we observe and the well-established 3 per cent compressional-wave anisotropy.

Elastic anisotropy of D″ predicted from global models of mantle flow

In order to test the hypothesis that seismic anisotropy in the lowermost mantle is caused by the development of a post-perovskite lattice preferred orientation, and that anisotropy can thus be used as a probe of the dynamics of the mantles lower boundary layer, an integrated model of texture generation in D″ is developed. This is used to predict the elastic anisotropy of the lowermost mantle as probed by global anisotropic tomographic inversions. The model combines the current 3D mantle flow field with simulations of the deformation of post-perovskite polycrystalline aggregates. Different descriptions of single crystal plasticity can lead to model results which are anti-correlated to each other. In models where post-perovskite deformation is accommodated by dislocations moving on (010) or (100), patterns of anisotropy are approximately correlated with the results of tomographic inversions. On the other hand, in models where dislocations move on (001) patterns of anisotropy are nearly anti-correlated with tomographic inversions. If all the seismic anisotropy in D″ extracted from global anisotropic inversions is due to the presence of a lattice preferred orientation in post-perovskite in the lowermost mantle, and if the results of the tomographic inversions are not strongly biased by the sampling geometries, these results suggest that, in contrast to ideas based on the 1D anisotropic signal, deformation of post-perovskite in the lowermost mantle may be accommodated by dislocations moving on (010) or (100). Alternatively, a significant portion of the anisotropic signal may be caused by mechanisms other than the alignment of post-perovskite crystals.

Differentiating flow, melt, or fossil seismic anisotropy beneath Ethiopia

Ethiopia is a region where continental rifting gives way to oceanic spreading. Yet the role that pre-existing lithospheric structure, melt, mantle flow, or active upwellings may play in this process is debated. Measurements of seismic anisotropy are often used to attempt to understand the contribution that these mechanisms may play. In this study, we use new data in Afar, Ethiopia along with legacy data across Ethiopia, Djibouti, and Yemen to obtain estimates of mantle anisotropy using SKS-wave splitting. We show that two layers of anisotropy exist, and we directly invert for these. We show that fossil anisotropy with fast directions oriented northeast-southwest may be preserved in the lithosphere away from the rift. Beneath the Main Ethiopian Rift and parts of Afar, anisotropy due to shear segregated melt along sharp changes in lithospheric thickness dominates the shear-wave splitting signal in the mantle. Beneath Afar, away from regions with significant lithospheric topography, melt pockets associated with the crustal and uppermost mantle magma storage dominate the signal in localized regions. In general, little anisotropy is seen in the uppermost mantle beneath Afar suggesting melt retains no preferential alignment. These results show the important role melt plays in weakening the lithosphere and imply that as rifting evolves passive upwelling sustains extension. A dominant northeast-southwest anisotropic fast direction is observed in a deeper layer across all of Ethiopia. This suggests that a conduit like plume is lacking beneath Afar today, rather a broad flow from the southwest dominates flow in the upper mantle.

Seismic anisotropy as an indicator of reservoir quality in siliciclastic rocks

Abstract Improving the accuracy of subsurface imaging is commonly the main incentive for including the effects of anisotropy in seismic processing. However, the anisotropy itself holds valuable information about rock properties and, as such, can be viewed as a seismic attribute. Here we summarize results from an integrated project that explored the potential to use observations of seismic anisotropy to interpret lithological and fluid properties (the SAIL project). Our approach links detailed petrofabric analyses of reservoir rocks, laboratory based measurements of ultrasonic velocities in core samples, and reservoir-scale seismic observations. We present results for the Clair field, a Carboniferous–Devonian reservoir offshore Scotland, west of the Shetland Islands. The reservoir rocks are sandstones that are variable in composition and exhibit anisotropy on three length-scales: the crystal, grain and fracture scale. We have developed a methodology for assessing crystal-preferred-orientation (CPO) using a combination of electron back-scattered diffraction (EBSD), X-ray texture goniometry (XRTG) and image analysis. Modal proportions of individual minerals are measured using quantitative X-ray diffraction (QXRD). These measurements are used to calculate the intrinsic anisotropy due to CPO via Voigt-Reuss-Hill averaging of individual crystal elasticities and their orientations. The intrinsic anisotropy of the rock is controlled by the phyllosilicate content and to a lesser degree the orientation of quartz and feldspar; the latter can serve as a palaeoflow indicator. Our results show remarkable consistency in CPO throughout the reservoir and allow us to construct a mathematical model of reservoir anisotropy. A comparison of CPO-predicted velocities and those derived from laboratory measurements of ultrasonic signals allows the estimation of additional elastic compliance terms due to grain-boundary interactions. The results show that the CPO estimates are good proxies for the intrinsic anisotropy of the clean sandstones. The more micaceous rocks exhibit enhanced anisotropy due to interactions between the phyllosilicate grains. We then compare the lab-scale predictions with reservoir-scale measurements of seismic anisotropy, based on amplitude variation with offset and azimuth (AVOA) analysis and non-hyperbolic moveout. Our mathematical model provides a foundation for interpreting the reservoir-scale seismic data and improving the geological modelling of complex reservoirs. The observed increases in AVOA signal with depth can only be explained with an increase in fracturing beneath the major unit boundaries, rather than a change in intrinsic CPO properties. In general, the style and magnitude of anisotropy in the Clair field appears to be indicative of reservoir quality.

Upper mantle anisotropy beneath the Seychelles microcontinent

[1] The Seychelles plateau is a prime example of a microcontinent, yet mechanisms for its creation and evolution are poorly understood. Recently acquired teleseismic data from a deployment of 26 stations on 18 islands in the Seychelles are analyzed to study upper mantle seismic anisotropy using SKS splitting results. Strong microseismic noise is attenuated using a polarization filter. Results show significant variation in time delays (δt = 0.4–2.4 s) and smooth variations in orientation (ϕ = 15°–69°, where ϕ is the polarization of the fast shear wave). The splitting results cannot be explained by simple asthenospheric flow associated with absolute plate motions. Recent work has suggested that anisotropy measurements for oceanic islands surrounding Africa can be explained by mantle flow due to plate motion in combination with density-driven flow associated with the African superswell. Such a mechanism explains our results only if there are lateral variations in the viscosity of the mantle. It has been suggested that the Seychelles are underlain by a mantle plume. Predictions of flow-induced anisotropy from plume-lithosphere interaction can explain our results with a plume possibly impinging beneath the plateau. Finally, we consider lithospheric anisotropy associated with rifting processes that formed the Seychelles. The large variation in the magnitude of shear wave splitting over short distances suggests a shallow source of anisotropy. Fast directions align parallel to an area of transform faulting in the Amirantes. Farther from this area the orientation of anisotropy aligns in similar directions as plate motions. This supports suggestions of transpressive deformation during the opening of the Mascarene basin.

Rapid Continental Breakup and Microcontinent Formation in the Western Indian Ocean

Two of the main factors that determine the nature of a rifted continental margin are rheology and magmatism during extension. Numerical models of lithospheric extension suggest that both factors vary with extension rate; yet until now extension rates of studied margins, as indicated by the rate of initial seafloor spreading, are mostly less than -30 mm/yr on each margin. This article presents the first geophysical results from the Seychelles-Laxmi Ridge conjugate pair of rifted margins which separated at -65 mm/yr. The Seychelles, with its spectacular exposures of Precambrian granite, was the earliest scientifically recognized microcontinent and arguably remains the classic example of one [Wegener, 1924; Matthews and Davies, 1966]. However, it is still unknown whether microcontinents result from plumes, changes in plate-boundary forces, lithospheric heterogeneity, or a combination of these factors.