Clinton P. Conrad

Global mantle flow and the development of seismic anisotropy : differences between the oceanic and continental upper mantle

[1] Viscous shear in the asthenosphere accommodates relative motion between Earth’s surface plates and underlying mantle, generating lattice-preferred orientation (LPO) in olivine aggregates and a seismically anisotropic fabric. Because this fabric develops with the evolving mantle flow field, observations of seismic anisotropy can constrain asthenospheric flow patterns if the contribution of fossil lithospheric anisotropy is small. We use global viscous mantle flow models to characterize the relationship between asthenospheric deformation and LPO and compare the predicted pattern of anisotropy to a global compilation of observed shear wave splitting measurements. For asthenosphere >500 km from plate boundaries, simple shear rotates the LPO toward the infinite strain axis (ISA, the LPO after infinite deformation) faster than the ISA changes along flow lines. Thus we expect the ISA to approximate LPO throughout most of the asthenosphere, greatly simplifying LPO predictions because strain integration along flow lines is unnecessary. Approximating LPO with the ISA and assuming A-type fabric (olivine a axis parallel to ISA), we find that mantle flow driven by both plate motions and mantle density heterogeneity successfully predicts oceanic anisotropy (average misfit 13). Continental anisotropy is less well fit (average misfit 41), but lateral variations in lithospheric thickness improve the fit in some continental areas. This suggests that asthenospheric anisotropy contributes to shear wave splitting for both continents and oceans but is overlain by a stronger layer of lithospheric anisotropy for continents. The contribution of the oceanic lithosphere is likely smaller because it is thinner, younger, and less deformed than its continental counterpart.

Seismic tomography, surface uplift, and the breakup of Gondwanaland: Integrating mantle convection backwards in time

Mantle density heterogeneities, imaged using seismic tomography, contain information about time-dependent mantle flow and mantle structures that existed in the past. We model the history of mantle flow using a tomographic image of the mantle beneath southern Africa as an initial condition while reversing the direction of flow and analytically incorporating cooling plates as a boundary condition. If the resulting (backwards integrated) model for structures is used as a starting point for a forwards convection model, todays mantle can be adequately reconstructed if we do not integrate backwards more than than about 50–75 Ma. Flow can also be reliably reversed through the Mesozoic, but only if instability of the lower boundary layer can be suppressed. Our model predicts that the large seismically-slow and presumably hot structure beneath southern Africa produced 500–700 m of dynamic topography throughout the Cenozoic. Since ~30 Ma, uplift has moved from eastern to southern Africa, where uplift rates are ~10 m/Myr, consistent with observations. During the Mesozoic, the modeled topographic high is situated near Gondwanaland rifting, raising the possibility that this buoyant structure may have been involved with this breakup.

Constraints on lithosphere net rotation and asthenospheric viscosity from global mantle flow models and seismic anisotropy

Although an average westward rotation of the Earths lithosphere is indicated by global analyses of surface features tied to the deep mantle (e.g., hot spot tracks), the rate of lithospheric drift is uncertain despite its importance to global geodynamics. We use a global viscous flow model to predict asthenospheric anisotropy computed from linear combinations of mantle flow fields driven by relative plate motions, mantle density heterogeneity, and westward lithosphere rotation. By comparing predictions of lattice preferred orientation to asthenospheric anisotropy in oceanic regions inferred from SKS splitting observations and surface wave tomography, we constrain absolute upper mantle viscosity (to 0.5–1.0 × 1021 Pa s, consistent with other constraints) simultaneously with net rotation rate and the decrease in the viscosity of the asthenosphere relative to that of the upper mantle. For an asthenosphere 10 times less viscous than the upper mantle, we find that global net rotation must be <0.26°/Myr (<60% of net rotation in the HS3 (Pacific hot spot) reference frame); larger viscosity drops amplify asthenospheric shear associated with net rotation and thus require slower net rotation to fit observed anisotropy. The magnitude of westward net rotation is consistent with lithospheric drift relative to Indo-Atlantic hot spots but is slower than drift in the Pacific hot spot frame (HS3 ≈ 0.44°/Myr). The latter may instead express net rotation relative to the deep mantle beneath the Pacific plate, which is moving rapidly eastward in our models.

Influence of continental roots and asthenosphere on plate‐mantle coupling

The shear tractions that mantle flow exerts on the base of Earths lithosphere contribute to plate-driving forces and lithospheric stresses. We investigate the sensitivity of these tractions to sub-lithospheric viscosity variations by comparing shear tractions computed from a mantle flow model featuring laterally-varying lithosphere and asthenosphere viscosity with those from a model with layered viscosity. Lateral viscosity variations generally do not change the direction of shear tractions, but deeply penetrating continental roots increase traction magnitudes by a factor of 2-5 compared to 100 km thick lithosphere. A low-viscosity asthenosphere decreases traction magnitudes by a smaller amount, and is important only if >100 km thick. Increased plate-mantle coupling beneath thick continental lithosphere may increase plate-driving forces, surface deformation, and mantle-derived lithospheric stresses in these regions. By contrast, a low-viscosity asthenosphere does not decouple the lithosphere from mantle flow, highlighting the geological importance of mantle tractions on the lithosphere.

The temporal evolution of plate driving forces: Importance of “slab suction” versus “slab pull” during the Cenozoic

Although mantle slabs ultimately drive plate motions, the mechanism by which they do so remains unclear. A detached slab descending through the mantle will excite mantle flow that exerts shear tractions on the base of the surface plates. This slab suction"" force drives subducting and overriding plates symmetrically toward subduction zones. Alternatively, cold, strong slabs may effectively transmit stresses to subducting surface plates, exerting a direct slab pull"" force on these plates, drawing them rapidly toward subduction zones. This motion induces mantle flow that pushes overriding plates away from subduction zones. We constrain the relative importance of slab suction and slab pull by comparing Cenozoic plate motions to model predictions that include viscous mantle flow and a proxy for slab strength. We find that slab pull from upper mantle slabs combined with slab suction from lower mantle slabs explains the observation that subducting plates currently move similar to4 times faster than nonsubducting plates. This implies that upper mantle slabs are strong enough to support their own weight. Slab suction and slab pull presently account for about 40 and 60% of the forces on plates, but slab suction only similar to30% if a low-viscosity asthenosphere decouples plates from mantle flow. The importance slab pull has been increasing steadily through the Cenozoic because the mass and length of upper mantle slabs has been increasing. This causes subducting plates to double their speed relative to nonsubducting plates during this time period. Our model explains this temporal evolution of plate motions for the first time.

Influence of dynamic topography on sea level and its rate of change

Mantle fl ow likely supports up to 2 km of long-wavelength topographic relief over Earth’s surface. Although the average of this dynamic support must be zero, a net defl ection of the ocean basins can change their volume and induce sea-level change. By calculating dynamic topography using a global mantle fl ow model, we fithat continents preferentially conceal depressed topography associated with mantle downwelling, leading to net seafl oor uplift and ~90 ± 20 m of positive sea-level offset. Upwelling mantle fl ow is currently amplifying positive dynamic topography and causing up to 1.0 m/Ma of sea-level rise, depending on mantle viscosity. Continental motions across dynamic topography gradients also affect sea level, but uncertainty over the plate motion reference frame permits sea-level rise or fall by ±0.3 m/Ma, depending on net lithosphere rotation. During a complete Wilson cycle, sea level should fall during supercontinent stability and rise during periods of dispersal as mantle fl ow pushes continents down dynamic topography gradients toward areas of mantle downwelling. We estimate that a maximum of ~1 m/Ma of sea-level rise may have occurred during the most recent continental dispersal. Because this rate is comparable in magnitude to other primary sea-level change mechanisms, dynamic offset of sea level by mantle fl ow should be considered a potentially signifi cant contributor to long-term sea-level change.

The thermal evolution of an Earth with strong subduction zones

It is commonly supposed that plate tectonic rates are controlled by the temperature-dependent viscosity of Earths deep interior. If this were so, a small decrease in mantle temperature would lead to a large decrease in global heat transport. This negative feedback mechanism would prevent mantle temperatures from changing rapidly with time. We propose alternatively that convection is primarily resisted by the bending of oceanic lithosphere at subduction zones. Because lithospheric strength should not depend strongly on interior mantle temperature, this relationship decreases the sensitivity of heat flow to changes in interior mantle viscosity, and thus permits more rapid temperature changes there. The bending resistance is large enough to limit heat flow rates for effective viscosities of the lithosphere greater than about 1023 Pa s, and increases with the cube of plate thickness. As a result, processes that affect plate thickness, such as small-scale convection or subduction initiation, could profoundly influence Earths thermal history.

The solid Earth’s influence on sea level

Because it lies at the intersection of Earth’s solid, liquid, and gaseous components, sea level links the dynamics of the fluid part of the planet with those of the solid part of the planet. Here, I review the past quarter century of sea-level research and show that the solid components of Earth exert a controlling influence on the amplitudes and patterns of sea-level change across time scales ranging from years to billions of years. On the shortest time scales (10 0 –10 2 yr), elastic deformation causes the ground surface to uplift instantaneously near deglaciating areas while the sea surface depresses due to diminished gravitational attraction. This produces spatial variations in rates of relative sea-level change (measured relative to the ground surface), with amplitudes of several millimeters per year. These sea-level “fingerprints” are characteristic of (and may help identify) the deglaciation source, and they can have significant societal importance because they will control rates of coastal inundation in the coming century. On time scales of 10 3 –10 5 yr, the solid Earth’s time-dependent viscous response to deglaciation also produces spatially varying patterns of relative sea-level change, with centimeters-per-year amplitude, that depend on the time-history of deglaciation. These variations, on average, cause net seafloor subsidence and therefore global sea-level drop. On time scales of 10 6 –10 8 yr, convection of Earth’s mantle also supports long-wavelength topographic relief that changes as continents migrate and mantle flow patterns evolve. This changing “dynamic topography” causes meters per millions of years of relative sea-level change, even along seemingly “stable” continental margins, which affects all stratigraphic records of Phanerozoic sea level. Nevertheless, several such records indicate sea-level drop of ∼230 m since a mid-Cretaceous highstand, when continental transgressions were occurring worldwide. This global drop results from several factors that combine to expand the “container” volume of the ocean basins. Most importantly, ridge volume decrease since the mid-Cretaceous, caused by an ∼50% slowdown in seafloor spreading rate documented by tectonic reconstructions, explains ∼250 m of sea-level fall. These tectonic changes have been accompanied by a decline in the volume of volcanic edifices on Pacific seafloor, continental convergence above the former Tethys Ocean, and the onset of glaciation, which dropped sea level by ∼40, ∼20, and ∼60 m, respectively. These drops were approximately offset by an increase in the volume of Atlantic sediments and net seafloor uplift by dynamic topography, which each elevated sea level by ∼60 m. Across supercontinental cycles, expected variations in ridge volume, dynamic topography, and continental compression together roughly explain observed sea-level variations throughout Pangean assembly and dispersal. On the longest time scales (10 9 yr), sea level may change as ocean water is exchanged with reservoirs stored by hydrous minerals within the mantle interior. Mantle cooling during the past few billion years may have accelerated drainage down subduction zones and decreased degassing at mid-ocean ridges, causing enough sea-level drop to impact the Phanerozoic sea-level budget. For all time scales, future advances in the study of sea-level change will result from improved observations of lateral variations in sea-level change, and a better understanding of the solid Earth deformations that cause them.

Faster seafloor spreading and lithosphere production during the mid-Cenozoic

Concurrent changes in seawater chemistry, sea level, and climate since the mid-Cretaceous are thought to result from an ongoing decrease in the global rate of lithosphere production at ridges. The present-day area distribution of seafloor ages, however, is most easily explained if lithosphere production rates were nearly constant during the past 180 m.y. We examined spatial gradients of present-day seafloor ages and inferred ages for the subducted Farallon plate to construct a history of spreading rates in each major ocean basin since ca. 140 Ma, revealing dramatic Cenozoic events. Globally, seafloor spreading rates increased by similar to 20% during the early Cenozoic due to an increase in plate speeds in the Pacific basin. Since then, subduction of the fast-spreading Pacific-Farallon ridge system has led to a 12% decrease in average global spreading rate and an 18% or more decrease in the total rate of lithosphere production by the most conservative estimates. These rapid changes during the Cenozoic defy models of steady-state seafloor formation, and demonstrate the time-dependent and evolving nature of plate tectonics on Earth.

Evidence of earthquake triggering by the solid earth tides

article i nfo Clearevidence forearthquake triggeringbythe earth tides hasremained elusive for morethan a century.Using the largest global earthquake catalog available (the NEIC catalog with 442412 events), we observe a clear correlation (with ∼99% confidence) between the phase of the solid Earth tide and the timing of seismic events: earthquakes occur slightly more often at the time of ground uplift by the Earth tide, i.e. when normal stresses are reduced within the lithosphere. We observe that this phase distribution anomaly is larger for smaller and shallowerearthquakes.Although earthquakes inregionswithdominantlynormalandstrike-slip faulting seem to exhibit more tidal triggering than regions dominated by thrust faulting, there is no statistically significant evidence for a focal mechanism-dependence on earthquake triggering. Finally, we show here that it is highly probablethat the observedtriggeringiscaused bythe solid Earth tide, ratherthan byloading fromthe oceanor atmospheric tides. Although an additional impact due to loading from ocean tides is possible and probable, we cannot detect it here because the earthquake database is not sufficiently complete and homogeneous (more smallmagnitudeearthquakesinoceanicareasareneeded).Ourresultsareconsistentwiththeideaofadamped