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.

High-resolution seismic constraints on flow dynamics in the oceanic asthenosphere

Convective flow in the mantle and the motions of tectonic plates produce deformation of Earth's interior, and the rock fabric produced by this deformation can be discerned using the anisotropy of the seismic wave speed. This deformation is commonly inferred close to lithospheric boundaries beneath the ocean in the uppermost mantle, including near seafloor-spreading centres as new plates are formed via corner flow, and within a weak asthenosphere that lubricates large-scale plate-driven flow and accommodates smaller scale convection. Seismic models of oceanic upper mantle differ as to the relative importance of these deformation processes: seafloor spreading fabric is very strong just beneath the crust-mantle boundary (the Mohorovičić discontinuity, or Moho) at relatively local scales, but at the global and ocean-basin scales, oceanic lithosphere typically appears weakly anisotropic when compared to the asthenosphere. Here we use Rayleigh waves, recorded across an ocean-bottom seismograph array in the central Pacific Ocean (the NoMelt Experiment), to provide unique localized constraints on seismic anisotropy within the oceanic lithosphere-asthenosphere system in the middle of a plate. We find that azimuthal anisotropy is strongest within the high-seismic-velocity lid, with the fast direction coincident with seafloor spreading. A minimum in the magnitude of azimuthal anisotropy occurs within the middle of the seismic low-velocity zone, and then increases with depth below the weakest portion of the asthenosphere. At no depth does the fast direction correlate with the apparent plate motion. Our results suggest that the highest strain deformation in the shallow oceanic mantle occurs during corner flow at the ridge axis, and via pressure-driven or buoyancy-driven flow within the asthenosphere. Shear associated with motion of the plate over the underlying asthenosphere, if present, is weak compared to these other processes.

Inferring nonlinear mantle rheology from the shape of the Hawaiian swell

The convective circulation generated within the Earth's mantle by buoyancy forces of thermal and compositional origin is intimately controlled by the rheology of the rocks that compose it. These can deform either by the diffusion of point defects (diffusion creep, with a linear relationship between strain rate and stress) or by the movement of intracrystalline dislocations (nonlinear dislocation creep). However, there is still no reliable map showing where in the mantle each of these mechanisms is dominant, and so it is important to identify regions where the operative mechanism can be inferred directly from surface geophysical observations. Here we identify a new observable quantity--the rate of downstream decay of the anomalous seafloor topography (swell) produced by a mantle plume--which depends only on the value of the exponent in the strain rate versus stress relationship that defines the difference between diffusion and dislocation creep. Comparison of the Hawaiian swell topography with the predictions of a simple fluid mechanical model shows that the swell shape is poorly explained by diffusion creep, and requires a dislocation creep rheology. The rheology predicted by the model is reasonably consistent with laboratory deformation data for both olivine and clinopyroxene, suggesting that the source of Hawaiian lavas could contain either or both of these components.

Rheology of the upper mantle: a synthesis

Rheological properties of the upper mantle of the Earth play an important role in the dynamics of the lithosphere and asthenosphere. However, such fundamental issues as the dominant mechanisms of flow have not been well resolved. A synthesis of laboratory studies and geophysical and geological observations shows that transitions between diffusion and dislocation creep likely occur in the Earth's upper mantle. The hot and shallow upper mantle flows by dislocation creep, whereas cold and shallow or deep upper mantle may flow by diffusion creep. When the stress increases, grain size is reduced and the upper mantle near the transition between these two regimes is weakened. Consequently, deformation is localized and the upper mantle is decoupled mechanically near these depths.

Mineral phases of the Earth´s mantle

Our knowledge of the structure of the Earth´s interior has been obtained by analysing seismic waves that travel in the Earth, and the reference Earth global models used by geophysicists are essentially seismological. Depth profiles of the seismic waves velocities reveal that the deep Earth is divided in several shells, separated by velocity and density discontinuities. The main discontinuity located at a depth of 2900 km corresponds to the transition between the mantle and the core. The Earth´s mantle can be further divided into the upper mantle and the lower mantle, with a transition zone characterised by two prominent increases in velocities observed at 410- and 660-km depths. This article will be focused on the mineral phases of the Earth´s mantle. The interpretation of seismological models in terms of chemical composition and temperature relies on the knowledge of the nature, structure and elastic properties of the candidate materials. We will describe to what extent recent advances in experimental mineral physics and X-ray diffraction have yielded essential knowledge on the structure and high-pressure high-temperature behaviour of pertinent materials, and major improvements in our understanding of the chemical and mineralogical composition of the Earth´s mantle.

The depth distribution of azimuthal anisotropy in the continental upper mantle

The most likely cause of seismic anisotropy in the Earth's upper mantle is the lattice preferred orientation of anisotropic minerals such as olivine. Its presence reflects dynamic processes related to formation of the lithosphere as well as to present-day tectonic motions. A powerful tool for detecting and characterizing upper-mantle anisotropy is the analysis of shear-wave splitting measurements. Because of the poor vertical resolution afforded by this type of data, however, it has remained controversial whether the splitting has a lithospheric origin that is 'frozen-in' at the time of formation of the craton, or whether the anisotropy originates primarily in the asthenosphere, and is induced by shear owing to present-day absolute plate motions. In addition, predictions from surface-wave-derived models are largely incompatible with shear-wave splitting observations. Here we show that this disagreement can be resolved by simultaneously inverting surface waveforms and shear-wave splitting data. We present evidence for the presence of two layers of anisotropy with different fast-axis orientations in the cratonic part of the North American upper mantle. At asthenospheric depths (200-400 km) the fast axis is sub-parallel to the absolute plate motion, confirming the presence of shear related to current tectonic processes, whereas in the lithosphere (80-200 km), the orientation is significantly more northerly. In the western, tectonically active, part of North America, the fast-axis direction is consistent with the absolute plate motion throughout the depth range considered, in agreement with a much thinner lithosphere.

Shear velocity structure in the Aegean region obtained by joint inversion of Rayleigh and Love waves

We present a shear velocity model of the crust and uppermost mantle under the Aegean region by simultaneous inversion of Rayleigh and Love waves. The database consists of regional earthquakes recorded by portable broadband three-component digital stations that were installed for a period of 6 months in the broader Aegean region. For each epicentre–station ray path group velocity dispersion curves are measured using appropriate frequency time analysis (FTAN). The dispersion measurements for more than 600 Love wave paths have been used. We have also incorporated previous results forc. 700 Rayleigh wave paths for the study area. The single-path dispersion curves of both waves were inverted to regional group velocity maps for different values of period (6–32 s) via a tomographic method. The local dispersion curves of discrete grid points for both surface waves were inverted nonlinearly to construct 1D models of shear-wave velocity v. depth. In most cases the joint inversion of Rayleigh and Love waves resulted in a single model (from the multiple models compatible with the data) that could interpret both Rayleigh and Love wave data. Around 60 local dispersion curves for both Rayleigh and Love waves were finally jointly inverted. As expected, because of the complex tectonic environment of the Aegean region the results show strong lateral variations of the S-wave velocities for the crust and uppermost mantle. Our results confirm the presence of a thin crust typically less than 28–30 km in the whole Aegean Sea, which in some parts of the southern and central Aegean Sea becomes significantly thinner (20–22 km). In contrast, a large crustal thickness of about 40–45 km exists in western Greece, and the remaining part of continental Greece is characterized by a mean crustal thickness of about 35 km. A significant sub-Moho upper mantle low-velocity zone (LVLmantle) with velocities as low as 3.7 km s−1, is clearly identified in the southern and central Aegean Sea, correlated with the high heat flow in the mantle wedge above the subducted slab and the related active volcanism in the region. The results obtained results are compared with independent body-wave tomographic information on the velocity structure of the study area and exhibit a generally good agreement, although significant small-scale differences are also identified.

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.

Inferences on Flow at the Base of Earth's Mantle Based on Seismic Anisotropy

We applied global waveform tomography to model radial anisotropy in the whole mantle. We found that in the last few hundred kilometers near the core-mantle boundary, horizontally polarized S-wave velocities (VSH) are, on average, faster (by approximately 1%) than vertically polarized S-wave velocities (VSV), suggesting a large-scale predominance of horizontal shear. This confirms that the D" region at the base of the mantle is also a mechanical boundary layer for mantle convection. A notable exception to this average signature can be found at the base of the two broad low-velocity regions under the Pacific Ocean and under Africa, often referred to as "superplumes," where the anisotropic pattern indicates the onset of vertical flow.

The evolution of mantle mixing

We present a geochemical overview of the canonic model that suggests a two-layer mantle for most of the Earth's history. A change in the Rayleigh number may have modified the convection and now allows the subduction of oceanic plates into the lower mantle, which was not the case in the past. The measurement of stirring time in the source of mid-ocean-ridge basalt, together with Xe- and Pb-isotopic ratios in the mid-ocean-ridge-basalt source, suggests that the upper mantle is separated into two domains, one above the 400 km discontinuity (asthenosphere) with rapid mixing and short residence time, and another between the 400 and 670 km discontinuities with sluggish mixing and a residence time of ca. 1.5 x 10(9) yr.

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.

Superplumes from the core-mantle boundary to the lithosphere: implications for heat flux

Three-dimensional modeling of upper-mantle anelastic structure reveals that thermal upwellings associated with the two superplumes, imaged by seismic elastic tomography at the base of the mantle, persist through the upper-mantle transition zone and are deflected horizontally beneath the lithosphere. This explains the unique transverse shear wave isotropy in the central Pacific. We infer that the two superplumes may play a major and stable role in supplying heat and horizontal flow to the low-viscosity asthenospheric channel, lubricating plate motions and feeding hot spots. We suggest that more heat may be carried through the core-mantle boundary than is accounted for by hot spot fluxes alone.

Electrical anisotropy below slow- and fast-moving plates: paleoflow in the upper mantle?

Upper mantle electrical conductivities can be explained by hydrogen diffusivity in hydrous olivine. Diffusivity enhances the conductivity of olivine anisotropically, making the a axis the most conductive of the three axes. Therefore, the hypothesis that plate motion induces lattice-preferred orientation of olivine can be tested with the use of long-period electromagnetic array measurements. Here, we compared electrical anisotropies below the slow-moving Fennoscandian and fast-moving Australian plates. The degree of olivine alignment is greater in the mantle below the Fennoscandian plate than below the Australian plate. This finding may indicate that convection rather than plate motion is the dominant deformation mechanism.

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.

Phase velocities of rayleigh waves in the MELT experiment on the east pacific rise

The phase velocities of Rayleigh waves increase more rapidly with distance from the East Pacific Rise (EPR) axis than is predicted by models of conductive cooling of the lithosphere. Low velocities near the axis are probably caused by partial melt at depths of 20 to 70 kilometers in a zone several hundred kilometers wide. The lowest velocities are offset to the west of the EPR. Wave propagation is anisotropic; the fast direction is approximately perpendicular to the ridge, parallel to the spreading direction. Anisotropy increases from a minimum near the axis to 3 percent or more on the flanks.