Lattice preferred orientation of olivine aggregates deformed in simple shear

Using the traditional microscope for mineral grain orientation determination: A prototype image analysis pipeline for optic-axis mapping (POAM)

This paper reports on the development of an open-source image analysis software 'pipeline' dedicated to petrographic microscopy. Using conventional rock thin sections and images from a standard polarising microscope, the pipeline can classify minerals and subgrains into objects and obtain information about optic-axis orientation. Five metamorphic rocks were chosen to test and illustrate the method. Thin sections were imaged using reflected and cross- and plane-polarised transmitted light. Images were taken at different angles of the polariser and analyser (360° with 10° steps), both with and without the full-lambda plate. The resulting image stacks were analysed with a modular pipeline for optic-axis mapping (POAM). POAM consists of external and internal software packages that register, segment, classify, and interpret the visible light spectra using object-based image analysis (OBIAS). The mapped fields-of-view and grain orientation stereonets of interest are presented in the context of whole-slide images. Two innovations are reported. First, we used hierarchical tree region merging on blended multimodal images to classify individual grains of rock-forming minerals into objects. Second, we assembled a new optical mineralogy algorithm chain that identifies the mineral slow axis orientation. The c-axis orientation results were verified with scanning electron microscopy electron backscattered diffraction (SEM-EBSD) data. For quartz (uniaxial) in a granite mylonite the test yielded excellent correspondence of c-axis azimuth and good agreement for inclination. For orthorhombic orthopyroxene in a deformed garnet harzburgite, POAM produced acceptable results for slow axis azimuth. In addition, the method identified slight anisotropy in garnet that would not be appreciated by traditional microscopy. We propose that our method is ideally suited for two commonly performed tasks in mineralogy. First, for mineral grain classification of entire thin sections scans on blended images to provide automated modal abundance estimates and grain size distribution. Second, for prospective fields of view of interest, POAM can rapidly generate slow axis crystal orientation maps from multiangle image stacks on conventionally prepared thin sections for targeting detailed SEM-EBSD studies.

Crustal structure and seismic anisotropy of rift basins in Somaliland

Rift margins provide insights into the processes governing the rupture of the continental lithosphere and the subsequence formation of sedimentary basins. The Proterozoic basement underlying Somaliland has been affected by multiple rifting; however, the crustal structure of these rifted basins remains unknown. This study utilized teleseismic receiver function analysis, Bayesian inversion, common conversion point imaging and 2D forward gravity modelling to examine the crust and upper mantle of Somaliland. The results indicate 36.8-38.2 km of crust in southern Somaliland, while the central and northern regions feature thinned crust (~ 21 km) with 5-6 km thick sediments. The joint analysis of radial and transverse components of receiver functions and shear wave splitting revealed fast axis directions trending to 50-56° in the upper mantle, indicating that azimuthal anisotropy is oriented in the regional Africa-Arabia plate motion. Such orientation may have resulted from lattice preferred orientation of olivine from the asthenospheric flow. Additionally, the fast polarization of the crust in central Somaliland is oriented at - 15°, indicating fossil deformation in the thinned crust related to the NW-SE trending Late Jurassic rift event. Further, the fast polarization for stations near the Gulf of Aden is oriented at 75-80°, suggesting crustal deformation associated with the Oligocene rift event. The crustal anisotropy at southern Somaliland revealed fast polarization oriented at - 85°, indicating a preserved far-field response of the WNW-ESE trending Late Cretaceous rift event. Overall, the study provides for the first-time insight into the rift-related extensional strain fabric in the crust and upper mantle anisotropy induced by asthenospheric flow in Somaliland.

Mantle flow underneath the South China Sea revealed by seismic anisotropy

It has long been established that plastic flow in the asthenosphere interacts constantly with the overlying lithosphere and plays a pivotal role in controlling the occurrence of geohazards such as earthquakes and volcanic eruptions. Unfortunately, accurately characterizing the direction and lateral extents of the mantle flow field is notoriously difficult, especially in oceanic areas where deployment of ocean bottom seismometers (OBSs) is expensive and thus rare. In this study, by applying shear wave splitting analyses to a dataset recorded by an OBS array that we deployed between mid-2019 and mid-2020 in the South China Sea (SCS), we show that the dominant mantle flow field has a NNW-SSE orientation, which can be attributed to mantle flow extruded from the Tibetan Plateau by the ongoing Indian-Eurasian collision. In addition, the results suggest that E-W oriented flow fields observed in South China and the Indochina Peninsula do not extend to the central SCS.

Influence of the asthenosphere on earth dynamics and evolution

The existence of a thin, weak asthenospheric layer beneath Earth's lithospheric plates is consistent with existing geological and geophysical constraints, including Pleistocene glacio-isostatic adjustment, modeling of gravity anomalies, studies of seismic anisotropy, and post-seismic rebound. Mantle convection models suggest that a pronounced weak zone beneath the upper thermal boundary layer (lithosphere) may be essential to the plate tectonic style of convection found on Earth. The asthenosphere is likely related to partial melting and the presence of water in the sub-lithospheric mantle, further implying that the long-term evolution of the Earth may be controlled by thermal regulation and volatile recycling that maintain a geotherm that approaches the wet mantle solidus at asthenospheric depths.

New shear wave splitting measurements using broadband seismic data recorded by the BX seismological network in Botswana

New measurements shear-wave splitting (SWS) were determined using a teleseismic dataset from 21 broadband seismic stations belonging to the GSN-BX network in Botswana archived at IRIS-DMC. The stations are distributed across major geologic terrains such as the Kaapvaal craton, Limpopo Mobile Belt and Zimbabwe Craton. Earthquakes with magnitude larger than 5.2 mb with epicentral distances in the range of 90° and 180° were analyzed for SKS (S-wave traversing the core as P-wave) and SKKS (S-wave traversing the core as P-wave with one reflection from inner side of the Core Mantle Boundary) phases, while PKS (P-wave bottoming in the core and converting to S-wave at the Core Mantle Boundary) phases were analyzed in the range of 130°-165°. SWS parameters were calculated at each station by minimizing the energy on the transverse component of the PKS, SKKS, and SKS phases (hereafter collectively called XKS). The measurement included the polarization direction ϕ, which gives a reliable indication of the alignment of olivine in the upper mantle, and, the delay time dt between the arrival times of the fast and slow travelling split shear waves. The value of dt is related to the intrinsic anisotropy and thickness of the anisotropic layer. The SWS parameters provide insights into the present and past deformation processes of upper mantle.

On the Relationship Between Oceanic Plate Speed, Tectonic Stress, and Seismic Anisotropy

Seismic radial anisotropy (the squared ratio between the speeds of horizontally and vertically polarized shear waves,ξ = V S H 2 V S V 2 ) is a powerful tool to probe the direction of mantle flow and accumulated strain. While previous studies have confirmed the dependence of azimuthal anisotropy on plate speed, the first order control on radial anisotropy is unclear. In this study, we develop 2D ridge flow models combined with mantle fabric calculations to report that faster plates generate higher tectonics stresses and strain rates which lower the dislocation creep viscosity and lead to deeper anisotropy than beneath slower plates. We apply the SGLOBE-rani tomographic filter, resulting in a flat depth-age trend and stronger anisotropy beneath faster plates, which correlates well with 3D global anisotropic mantle models. Our predictions and observations suggest that as plate speed increases from 2 to 8 cm/yr, radial anisotropy increases by ∼0.01-0.025 in the upper 100-200 km of the mantle between 10 and 60 Ma.

Mantle-flow diversion beneath the Iranian plateau induced by Zagros' lithospheric keel

Previous investigation of seismic anisotropy indicates the presence of a simple mantle flow regime beneath the Turkish-Anatolian Plateau and Arabian Plate. Numerical modeling suggests that this simple flow is a component of a large-scale global mantle flow associated with the African superplume, which plays a key role in the geodynamic framework of the Arabia-Eurasia continental collision zone. However, the extent and impact of the flow pattern farther east beneath the Iranian Plateau and Zagros remains unclear. While the relatively smoothly varying lithospheric thickness beneath the Anatolian Plateau and Arabian Plate allows progress of the simple mantle flow, the variable lithospheric thickness across the Iranian Plateau is expected to impose additional boundary conditions on the mantle flow field. In this study, for the first time, we use an unprecedented data set of seismic waveforms from a network of 245 seismic stations to examine the mantle flow pattern and lithospheric deformation over the entire region of the Iranian Plateau and Zagros by investigation of seismic anisotropy. We also examine the correlation between the pattern of seismic anisotropy, plate motion using GPS velocities and surface strain fields. Our study reveals a complex pattern of seismic anisotropy that implies a similarly complex mantle flow field. The pattern of seismic anisotropy suggests that the regional simple mantle flow beneath the Arabian Platform and eastern Turkey deflects as a circular flow around the thick Zagros lithosphere. This circular flow merges into a toroidal component beneath the NW Zagros that is likely an indicator of a lateral discontinuity in the lithosphere. Our examination also suggests that the main lithospheric deformation in the Zagros occurs as an axial shortening across the belt, whereas in the eastern Alborz and Kopeh-Dagh a belt-parallel horizontal lithospheric deformation plays a major role.

Seismic evidence for subduction-induced mantle flows underneath Middle America

Laboratory experiments and geodynamic simulations demonstrate that poloidal- and toroidal-mode mantle flows develop around subduction zones. Here, we use a new 3-D azimuthal anisotropy model constructed by full waveform inversion, to infer deep subduction-induced mantle flows underneath Middle America. At depths shallower than 150 km, poloidal-mode flow is perpendicular to the trajectory of the Middle American Trench. From 300 to 450 km depth, return flows surround the edges of the Rivera and Atlantic slabs, while escape flows are inferred through slab windows beneath Panama and central Mexico. Furthermore, at 700 km depth, the study region is dominated by the Farallon anomaly, with fast axes perpendicular to its strike, suggesting the development of lattice-preferred orientations by substantial stress. These observations provide depth-dependent seismic anisotropy for future mantle flow simulations, and call for further investigations about the deformation mechanisms and elasticity of minerals in the transition zone and uppermost lower mantle.

The motions of subducted slabs are expected to drive mantle flow around slab edges, however, evidence of deep mantle flow has so far remained elusive. Here, the authors present a Full Waveform Inversion 3-D anisotropy model which allows them to infer deep subduction-induced mantle flows underneath the Mid-Americas and the Caribbean.

Microstructural constraints on magmatic mushes under Kīlauea Volcano, Hawai'i

Distorted olivines of enigmatic origin are ubiquitous in erupted products from a wide range of volcanic systems (e.g., Hawaiʻi, Iceland, Andes). Investigation of these features at Kīlauea Volcano, Hawaiʻi, using an integrative crystallographic and chemical approach places quantitative constraints on mush pile thicknesses. Electron backscatter diffraction (EBSD) reveals that the microstructural features of distorted olivines, whose chemical composition is distinct from undistorted olivines, are remarkably similar to olivines within deformed mantle peridotites, but inconsistent with an origin from dendritic growth. This, alongside the spatial distribution of distorted grains and the absence of adcumulate textures, suggests that olivines were deformed within melt-rich mush piles accumulating within the summit reservoir. Quantitative analysis of subgrain geometry reveals that olivines experienced differential stresses of ∼3–12 MPa, consistent with their storage in mush piles with thicknesses of a few hundred metres. Overall, our microstructural analysis of erupted crystals provides novel insights into mush-rich magmatic systems.

Olivine crystals with prominent intracrystalline distortions have previously been used to quantify deformational processes within the mantle. Here, the authors show that similar techniques can be applied to deformed volcanic olivine crystals, providing quantitative constraints on the geometry of melt-rich mush piles within magmatic plumbing systems.

The Role of Oceanic Transform Faults in Seafloor Spreading: A Global Perspective From Seismic Anisotropy

Mantle anisotropy beneath mid-ocean ridges and oceanic transforms is key to our understanding of seafloor spreading and underlying dynamics of divergent plate boundaries. Observations are sparse, however, given the remoteness of the oceans and the difficulties of seismic instrumentation. To overcome this, we utilize the global distribution of seismicity along transform faults to measure shear wave splitting of over 550 direct S phases recorded at 56 carefully selected seismic stations worldwide. Applying this source-side splitting technique allows for characterization of the upper mantle seismic anisotropy, and therefore the pattern of mantle flow, directly beneath seismically active transform faults. The majority of the results (60%) return nulls (no splitting), while the non-null measurements display clear azimuthal dependency. This is best simply explained by anisotropy with a near vertical symmetry axis, consistent with mantle upwelling beneath oceanic transforms as suggested by numerical models. It appears therefore that the long-term stability of seafloor spreading may be associated with widespread mantle upwelling beneath the transforms creating warm and weak faults that localize strain to the plate boundary.

Evidence for frozen melts in the mid-lithosphere detected from active-source seismic data

The interactions of the lithospheric plates that form the Earth's outer shell provide much of the evidentiary basis for modern plate tectonic theory. Seismic discontinuities in the lithosphere arising from mantle convection and plate motion provide constraints on the physical and chemical properties of the mantle that contribute to the processes of formation and evolution of tectonic plates. Seismological studies during the past two decades have detected seismic discontinuities within the oceanic lithosphere in addition to that at the lithosphere-asthenosphere boundary (LAB). However, the depth, distribution, and physical properties of these discontinuities are not well constrained, which makes it difficult to use seismological data to examine their origin. Here we present new active-source seismic data acquired along a 1,130 km profile across an old Pacific plate (148-128 Ma) that show oceanic mid-lithosphere discontinuities (oceanic MLDs) distributed 37-59 km below the seafloor. The presence of the oceanic MLDs suggests that frozen melts that accumulated at past LABs have been preserved as low-velocity layers within the current mature lithosphere. These observations show that long-offset, high-frequency, active-source seismic data can be used to image mid-lithospheric structure, which is fundamental to understanding the formation and evolution of tectonic plates.

Olivine anisotropy suggests Gutenberg discontinuity is not the base of the lithosphere

Tectonic plates are a key feature of Earth's structure, and their behavior and dynamics are fundamental drivers in a wide range of large-scale processes. The operation of plate tectonics, in general, depends intimately on the manner in which lithospheric plates couple to the convecting interior. Current debate centers on whether the transition from rigid lithosphere to flowing asthenosphere relates to increases in temperature or to changes in composition such as the presence of a small amount of melt or an increase in water content below a specified depth. Thus, the manner in which the rigid lithosphere couples to the flowing asthenosphere is currently unclear. Here we present results from laboratory-based torsion experiments on olivine aggregates with and without melt, yielding an improved database describing the crystallographic alignment of olivine grains. We combine this database with a flow model for oceanic upper mantle to predict the structure of the seismic anisotropy beneath ocean basins. Agreement between our model and seismological observations supports the view that the base of the lithosphere is thermally controlled. This model additionally supports the idea that discontinuities in velocity and anisotropy, often assumed to be the base of the lithosphere, are, instead, intralithospheric features reflecting a compositional boundary established at midocean ridges, not a rheological boundary.

In situ analysis of texture development from sinusoidal stress at high pressure and temperature

Here, we present a new experimental protocol to investigate the relationship between texture, plastic strain, and the mechanisms of plastic deformation at high pressure and temperature. The method utilizes synchrotron X-ray radiation as the probing tool, coupled with a large-volume high pressure deformation device (D-DIA). The intensity of X-ray diffraction peaks within the spectrum of the sample is used for sampling texture development in situ. The unique feature of this study is given by the sinusoidal variation of the intensity when a sinusoidal strain is applied to the sample. For a sample of magnesium oxide at elevated pressure and temperature, we demonstrate observations that are consistent with elasto-plastic models for texture development and for diffraction-peak measurements of apparent stress. The sinusoidal strain magnitude was 3%.

Single-crystal Brillouin spectroscopy with CO2 laser heating and variable q

We describe a Brillouin spectroscopy system integrated with CO2 laser-heating and Raman spectroscopic capabilities. Temperature is determined by measurements of the grey-body thermal radiation emitted by the hot sample, with the system response calibrated relative to a standard tungsten ribbon lamp. High-pressure laser-heating Brillouin scattering measurements of acoustic velocities on liquid water and ice compressed in a diamond-anvil cell were performed at temperatures up to 2500 ± 150 K at high pressure. Single-crystal laser-heating Brillouin measurements were made on the (111) plane of San Carlos olivine at ∼13 GPa, 1300 ± 200 K. The pressure as measured by ruby fluorescence is shown to be within ±0.5 GPa of the pressure on the olivine sample during laser heating when KCl and KBr are used as pressure-transmitting media. In addition, the system is designed for continuously variable scattering angles from forward scattering (near 0° scattering angle) up to near back scattering (∼141°). This novel setup allows us to probe a wide range of wave vectors q for investigation of phonon dispersion on, for example, crystals with large unit cells (on the scale of hundreds of nm).

Crystal preferred orientation of an amphibole experimentally deformed by simple shear

Seismic anisotropy has been widely observed in crust and mantle materials and plays a key role in the understanding of structure and flow patterns. Although seismic anisotropy can be explained by the crystal preferred orientation (CPO) of highly anisotropic minerals in the crust, that is, amphibole, experimental studies on the CPO of amphibole are limited. Here we present the results of novel experiments on simple shear deformation of amphibolite at high pressure and temperatures (1 GPa, 480-700 °C). Depending on the temperature and stress, the deformed amphibole produced three types of CPOs and resulted in a strong seismic anisotropy. Our data provide a new understanding of the observed seismic anisotropy. The seismic data obtained from the amphibole CPOs revealed that anomalous seismic anisotropy observed in the deep crust, subducting slab and mantle wedge can be attributed to the CPO of amphibole.

Crust-mantle mechanical coupling in Eastern Mediterranean and eastern Turkey

Present-day crust-mantle coupling in the Eastern Mediterranean and eastern Turkey is studied using the Global Positioning System (GPS) and seismic anisotropy data. The general trend of the shear wave fast-splitting directions in NE Turkey and Lesser Caucaus align well with the geodetic velocities in an absolute plate motion frame of reference pointing to an effective coupling in this part of the region of weak surface deformation. Farther south, underneath the Bitlis Suture, however, there are significant Pn delays with E-W anisotropy axes indicating significant lateral escape. Meanwhile, the GPS reveals very little surface deformation. This mismatch possibly suggests a decoupling along the suture. In the Aegean, the shear wave anisotropy and the Pn anisotropy directions agree with the extensional component of the right-lateral shear strains except under the Crete Basin and other parts of the southern Aegean Sea. This extensional direction matches perfectly also with the southward pulling force vectors across the Hellenic trench; however, the maximum right-lateral shear directions obtained from the GPS data in the Aegean do not match either of these anisotropies. Seismic anisotropy from Rayleigh waves sampled at 15 s, corresponding to the lower crust, match the maximum right-lateral maximum shear directions from the GPS indicating decoupling between the crust and the mantle. This decoupling most likely results from the lateral variations of the gravitational potential energies and the slab-pull forces.

Global azimuthal seismic anisotropy and the unique plate-motion deformation of Australia

Differences in the thickness of the high-velocity lid underlying continents as imaged by seismic tomography, have fuelled a long debate on the origin of the 'roots' of continents. Some of these differences may be reconciled by observations of radial anisotropy between 250 and 300 km depth, with horizontally polarized shear waves travelling faster than vertically polarized ones. This azimuthally averaged anisotropy could arise from present-day deformation at the base of the plate, as has been found for shallower depths beneath ocean basins. Such deformation would also produce significant azimuthal variation, owing to the preferred alignment of highly anisotropic minerals. Here we report global observations of surface-wave azimuthal anisotropy, which indicate that only the continental portion of the Australian plate displays significant azimuthal anisotropy and strong correlation with present-day plate motion in the depth range 175-300 km. Beneath other continents, azimuthal anisotropy is only weakly correlated with plate motion and its depth location is similar to that found beneath oceans. We infer that the fast-moving Australian plate contains the only continental region with a sufficiently large deformation at its base to be transformed into azimuthal anisotropy. Simple shear leading to anisotropy with a plunging axis of symmetry may explain the smaller azimuthal anisotropy beneath other continents.

Dynamic recrystallization and strain softening of olivine aggregates in the laboratory and the lithosphere

The effects of dynamic recrystallization on the deformation mechanisms and rheology of olivine aggregates in the laboratory and the lithosphere are reviewed in this paper. The low-strain rheology of olivine is well documented; however, deformation in the lithosphere often involves large strains. Large strain experiments show that recrystallization can result in both hardening and softening during deformation. Moderate strain softening in experimental shear and torsion can be explained by the operation of dislocation-accommodated grain boundary sliding in bands of fine recrystallized grains.

Data on the temperature dependence of recrystallized grain size are needed to extrapolate the effects of dynamic recrystallization to the lithosphere. Theories of dynamic recrystallization suggest that grain size is strongly stress dependent and moderately temperature dependent. A re-analysis of experimental grain size data indicates that the recrystallized grain size is temperature independent for olivine aggregates with low water content (<300 ppm H/Si).

Rheological regime maps have been constructed for the lithospheric mantle. The maps suggest that grain size sensitive power law creep, involving both grain boundary sliding and dislocation creep, will produce strong strain softening, greater than found so far in experimental studies, in dry and wet lithosphere shear zones.

High-strain zones: laboratory perspectives on strain softening during ductile deformation

Deformation in the Earth’s outer shell is mostly localized into narrow high-strain zones. Because they can have displacements up to several hundreds or thousands of kilometres, they can affect the entire lithosphere. The properties of high-strain zones control the kinematics and dynamics of our planet, and are therefore of key importance for an understanding of plate tectonics, stress accumulation and release (e.g. earthquakes), mountain building, etc.

One of the requirements of shear zone formation in ductile rocks is localized strain softening (Hobbs et al. 1990). In this paper we review the strain softening mechanisms that were identified and proposed 25 years ago and analyse their relevance in light of recent experimental results conducted to large strains. For this purpose, some of the newer developments in experimental deformation techniques that permit high strain in torsion are summarized and recent results are reviewed. Using these results we discuss mechanisms, processes and conditions that lead to localization.

Natural examples of olivine lattice preferred orientation patterns with a flow-normal a-axis maximum

Tectonic plate motion is thought to cause solid-state plastic flow within the underlying upper mantle and accordingly lead to the development of a lattice preferred orientation of the constituent olivine crystals. The mechanical anisotropy that results from such preferred orientation typically produces a direction of maximum seismic wave velocity parallel to the plate motion direction. This has been explained by the existence of an olivine preferred orientation with an 'a-axis' maximum parallel to the induced mantle flow direction. In subduction zones, however, the olivine a axes have been inferred to be arranged roughly perpendicular to plate motion, which has usually been ascribed to localized complex mantle flow patterns. Recent experimental work suggests an alternative explanation: under conditions of high water activity, a 'B-type' olivine preferred orientation may form, with the a-axis maximum perpendicular to the flow direction. Natural examples of such B-type preferred orientation are, however, almost entirely unknown. Here we document widespread B-type olivine preferred orientation patterns from a subduction-type metamorphic belt in southwest Japan and show that these patterns developed in the presence of water. Our discovery implies that mantle flow above subduction zones may be much simpler than has generally been thought.

Boundary-layer mantle flow under the Dead Sea transform fault inferred from seismic anisotropy

Lithospheric-scale transform faults play an important role in the dynamics of global plate motion. Near-surface deformation fields for such faults are relatively well documented by satellite geodesy, strain measurements and earthquake source studies, and deeper crustal structure has been imaged by seismic profiling. Relatively little is known, however, about deformation taking place in the subcrustal lithosphere--that is, the width and depth of the region associated with the deformation, the transition between deformed and undeformed lithosphere and the interaction between lithospheric and asthenospheric mantle flow at the plate boundary. Here we present evidence for a narrow, approximately 20-km-wide, subcrustal anisotropic zone of fault-parallel mineral alignment beneath the Dead Sea transform, obtained from an inversion of shear-wave splitting observations along a dense receiver profile. The geometry of this zone and the contrast between distinct anisotropic domains suggest subhorizontal mantle flow within a vertical boundary layer that extends through the entire lithosphere and accommodates the transform motion between the African and Arabian plates within this relatively narrow zone.

Melt segregation and strain partitioning: implications for seismic anisotropy and mantle flow

One of the principal means of understanding upper mantle dynamics involves inferring mantle flow directions from seismic anisotropy under the assumption that the seismic fast direction (olivine a axis) parallels the regional flow direction. We demonstrate that (i) the presence of melt weakens the alignment of a axes and (ii) when melt segregates and forms networks of weak shear zones, strain partitions between weak and strong zones, resulting in an alignment of a axes 90 degrees from the shear direction in three-dimensional deformation. This orientation of a axes provides a new means of interpreting mantle flow from seismic anisotropy in partially molten deforming regions of Earth.

Seismic anisotropy: tracing plate dynamics in the mantle

Elastic anisotropy is present where the speed of a seismic wave depends on its direction. In Earth's mantle, elastic anisotropy is induced by minerals that are preferentially oriented in a directional flow or deformation. Earthquakes generate two seismic wave types: compressional (P) and shear (S) waves, whose coupling in anisotropic rocks leads to scattering, birefringence, and waves with hybrid polarizations. This varied behavior is helping geophysicists explore rock textures within Earth's mantle and crust, map present-day upper-mantle convection, and study the formation of lithospheric plates and the accretion of continents in Earth history.

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 Earth's 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 heterogeneous 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.

The mantle flow field beneath western North America

Although motions at the surface of tectonic plates are well determined, the accompanying horizontal mantle flow is not. We have combined observations of surface deformation and upper mantle seismic anisotropy to estimate this flow field for western North America. We find that the mantle velocity is 5.5 +/- 1.5 centimeters per year due east in a hot spot reference frame, nearly opposite to the direction of North American plate motion (west-southwest). The flow is only weakly coupled to the motion of the surface plate, producing a small drag force. This flow field is probably due to heterogeneity in mantle density associated with the former Farallon oceanic plate beneath North America.

Seismic evidence for hotspot-induced buoyant flow beneath the Reykjanes Ridge

Volcanic hotspots and mid-ocean ridge spreading centers are the surface expressions of upwelling in Earth's mantle convection system, and their interaction provides unique information on upwelling dynamics. I investigated the influence of the Iceland hotspot on the adjacent mid-Atlantic spreading center using phase-delay times of seismic surface waves, which show anomalous polarization anisotropy-a delay-time discrepancy between waves with different polarizations. This anisotropy implies that the hotspot induces buoyancy-driven upwelling in the mantle beneath the ridge.

Water-induced fabric transitions in olivine

The interpretation of seismic anisotropy in Earth's upper mantle has traditionally been based on the fabrics (lattice-preferred orientation) of relatively water-poor olivine. Here we show that when a large amount of water is added to olivine, the relation between flow geometry and seismic anisotropy undergoes marked changes. Some of the puzzling observations of seismic anisotropy in the upper mantle, including the anomalous anisotropy in the central Pacific and the complicated anisotropy in subduction zones, can be attributed to the enrichment of water in these regions.

A complex pattern of mantle flow in the Lau backarc

Shear-wave splitting analysis of local events recorded on land and on the ocean floor in the Tonga arc and Lau backarc indicate a complex pattern of azimuthal anisotropy that cannot be explained by mantle flow coupled to the downgoing plate. These observations suggest that the direction of mantle flow rotates from convergence-parallel in the Fiji plateau to north-south beneath the Lau basin and arc-parallel beneath the Tonga arc. These results correlate with helium isotopes that map mantle flow of the Samoan plume into the Lau basin through an opening tear in the Pacific plate.