Cenozoic plate driving forces

Past studies of plate driving forces have concluded that the forces due to subducted slabs in the upper mantle and those due to the thickening of the oceanic lithosphere are the principal driving forces. We reexamine the balance of driving forces for the present-day and extend our analysis through the Cenozoic, using an analytical torque balance method which accounts for interactions between plates via viscous coupling to the induced mantle flow, We use an evolving mantle density heterogeneity field based on the last 200 Myr. of subduction to drive plate motions, an approach which has proven successful in predicting the present-day mantle heterogeneity field. We find that for plausible upper mantle viscosities the forces due to subducted slabs in the Cenozoic and Mesozoic account for in excess of 90% of plate driving forces and those due to lithospheric thickening for less than 10%

The Dynamic Structure of the Deep Earth: An Interdisciplinary Approach

Peer Reviewedhttp://deepblue.lib.umich.edu/bitstream/2027.42/94755/1/eost14654.pd

Triggered seismicity associated with the 1990 Nicoya, Costa Rica, M-w=7.0 earthquake

The 25 March 1990 (M-w = 7.0) subduction megathrust earthquake that occurred offshore the Nicoya Peninsula, Costa Rica, produced a large number of aftershocks on the subduction plate interface as expected and preceded an unusual sequence of earthquakes 75 km inland that had two periods of significant increase, one at 60-90 days and one near 270 days, following the main shock. This inland sequence of events would not typically fall within the classification of aftershocks given their spatial and temporal distance, and we show here that this sequence was likely triggered by the 25 March main shock. We compute stress changes on representative faults within this inland region using both a simple half-space model as well as with a 2-D finite element model that incorporates variable rheologic properties. The half-space model predicts a minor increase in Coulomb stress changes and a large amount of unclamping in this region, likely enough to cause triggering on the inland right-lateral strike-slip faults. Models that include a viscoelastic response also indicate stress increases that may link to triggering, particularly related to the time delay. Earthquakes on the subduction zone thrust along Costa Rica should be considered in hazard assessments for the inland populated region as several sets of strike-slip faults have been mapped in the fore-arc region

Dynamics and excess temperature of a plume throughout its life cycle

Measurements of the velocity field associated with plumes rising through a viscous fluid are performed using stereoscopic Particle-Image Velocimetry in the Rayleigh number range 4.4 × 105 − 6.4 × 105. The experimental model is analogous to a mantle plume rising from the core-mantle boundary to the base of the lithosphere. The behaviour of the plume is studied throughout its life cycle, which is broken up into four stages; the Formation Stage, when the plume forms; the Rising Stage, when the plume rises through the fluid; the Spreading Stage, when the plume reaches the surface and spreads; and finally the Declining Stage, when the heat source has been removed and the plume weakens. The latter three stages are examined in terms of the Finite-Time Lyapunov Exponent fields and the advection of passive tracers throughout the flow. The temperature at the heater and near the fluid surface are measured using thermocouples to infer how the presence of a mantle plume would produce excess temperature near the lithosphere throughout the various stages of its life cycle. In all experiments a time lag is observed between the removal of the heat source and the decline in the excess temperature near the surface, which is proportional to the rise time. A simple analytical model is presented, which suggests that under mantle conditions (i.e. negligible thermal diffusion), the relationship between the time lag and the rise time is robust and independent of the Rayleigh number; however, the constant of proportionality is closer to unity in the absence of diffusion. Once the heat source is removed, the excess temperature near the surface declines exponentially at a rate that is inversely proportional to the rise time. The implications of this result are discussed in terms of the decline in volcanism in the Louisville hotspot chain over the past 20 Ma. The rise velocity of material in the plume is examined; the rise velocity is found to vary significantly with the plume height in a manner that is inconsistent with many of the common semi-analytical models of thermal plumes in the literature. It is also argued that this height-dependency will cause estimates of the rise velocity based on the decay series of Uranium isotopes to significantly underestimate the true value

Thermodynamics of mantle minerals – I. Physical properties

We present a theory for the computation of phase equilibria and physical properties of multicomponent assemblages relevant to the mantle of the Earth. The theory differs from previous treatments in being thermodynamically self-consistent: the theory is based on the concept of fundamental thermodynamic relations appropriately generalized to anisotropic strain and in encompassing elasticity in addition to the usual isotropic thermodynamic properties. In this first paper, we present the development of the theory, discuss its scope, and focus on its application to physical properties of mantle phases at elevated pressure and temperature including the equation of state, thermochemical properties and the elastic wave velocities. We find that the Eulerian finite strain formulation captures the variation of the elastic moduli with compression. The variation of the vibrational frequencies with compression is also cast as a Taylor series expansion in the Eulerian finite strain, the appropriate volume derivative of which leads to an expression for the GrÜneisen parameter that agrees well with results from first principles theory. For isotropic materials, the theory contains nine material-specific parameters: the values at ambient conditions of the Helmholtz free energy, volume, bulk and shear moduli, their pressure derivatives, an effective Debye temperature, its first and second logarithmic volume derivatives (Γ 0 , q 0 ) , and the shear strain derivative of Γ. We present and discuss in some detail the results of a global inversion of a wide variety of experimental data and first principles theoretical results, supplemented by systematic relations, for the values of these parameters for 31 mantle species. Among our findings is that the value of q is likely to be significantly greater than unity for most mantle species. We apply the theory to the computation of the shear wave velocity, and temperature and compositional (Fe content) derivatives at relevant mantle pressure temperature conditions. Among the patterns that emerge is that garnet is anomalous in being remarkably insensitive to iron content or temperature as compared with other mantle phases.Peer Reviewedhttp://deepblue.lib.umich.edu/bitstream/2027.42/73485/1/j.1365-246X.2005.02642.x.pd

Why is Africa rifting?

Continental rifting has a fundamental role in the tectonic behaviour of the Earth, shaping the surface we live on. Although there is not yet a consensus about the dominant mechanism for rifting, there is a general agreement that the stresses required to rift the continental lithosphere are not readily available. Here we use a global finite element model of the lithosphere to calculate the stresses acting on Africa. We consider the stresses induced by mantle flow, crustal structure and topography in two types of models: one in which flow is exclusively driven by the subducting slabs and one in which it is derived from a shear wave tomographic model. The latter predicts much larger stresses and a more realistic dynamic topography. It is therefore clear that the mantle structure beneath Africa plays a key part in providing the radial and horizontal tractions, dynamic topography and gravitational potential energy necessary for rifting. Nevertheless, the total available stress (c. 100 MPa) is much less than that needed to break thick, cold continental lithosphere. Instead, we appeal to a model of magma-assisted rifting along pre-existing weaknesses, where the strain is localized in a narrow axial region and the strength of the plate is reduced significantly. Mounting geological and geophysical observations support such a model

Why is Africa rifting?

Continental rifting has a fundamental role in the tectonic behaviour of the Earth, shaping the surface we live on. Although there is not yet a consensus about the dominant mechanism for rifting, there is a general agreement that the stresses required to rift the continental lithosphere are not readily available. Here we use a global finite element model of the lithosphere to calculate the stresses acting on Africa. We consider the stresses induced by mantle flow, crustal structure and topography in two types of models: one in which flow is exclusively driven by the subducting slabs and one in which it is derived from a shear wave tomographic model. The latter predicts much larger stresses and a more realistic dynamic topography. It is therefore clear that the mantle structure beneath Africa plays a key part in providing the radial and horizontal tractions, dynamic topography and gravitational potential energy necessary for rifting. Nevertheless, the total available stress (c. 100 MPa) is much less than that needed to break thick, cold continental lithosphere. Instead, we appeal to a model of magma-assisted rifting along pre-existing weaknesses, where the strain is localized in a narrow axial region and the strength of the plate is reduced significantly. Mounting geological and geophysical observations support such a model

Orphaning Regimes: The Missing Link Between Flattened and Penetrating Slab Morphologies

Slab orphaning is a newly discovered phenomenological behavior, where the slab tip breaks off at the top of the lower mantle (~660 km depth) and is abandoned by its parent slab. Upon orphaning, subduction continues uninterrupted through the lateral motion of the parent slab above 660 km depth. In this work, we present a regime diagram for the range of conditions under which slabs can orphan at the top of the lower mantle. Our models show that a viscosity jump at 1,000 km depth not coincident with the endothermic phase change responsible for the 660 km seismic discontinuity, is necessary for orphaning as is the presence of a low viscosity channel between 660 and 1,000 km depth. We show that orphan slabs, similar to other deep slab morphologies, can be the end result for a wide range of physical parameters governing slab dynamics: slab orphaning persists across wide variations in slab dip, slab yield stress/strength, Clapeyron slope values, and overriding plate nature. The diversity in orphan slab sizes and orphaning periods is tied to the orphaning regime space, which describes a hitherto unexplored region between deflected and penetrating deep-subduction modes. Orphaning provides a simple dynamic link between the well-known deflection and penetration, and provides one possible way for slabs to switch from direct penetration to deflection, littering the mantle with abandoned fragments. Orphan slabs are therefore the intermediary between these two extensively studied slab morphologies

Influence of continental roots and asthenosphere on plate‐mantle coupling

Peer Reviewedhttp://deepblue.lib.umich.edu/bitstream/2027.42/94652/1/grl21239.pd

Surface exposure constraints on the mantle water budget

Mantle water content estimates range from 0.5 to 15 oceans of water. Its evolution is even more unclear. Rapid degassing during mantle solidification likely released much of the water to the surface, initially flooding Earth. However, evidence for subaerial land from at least 3.5 Ga means that much of this water must have been rapidly cycled back into the mantle. Here, we used a parameterized convection model and hypsometric curve to assess how much water could have been taken into the mantle and still satisfy evidence for subaerial land. Even if only the highest peaks were exposed, the initial ocean must have been less than 1.5 current oceans to explain subaerial exposure throughout most of Earth history. Today, this implies any water in the mantle?>0.5 oceans must be primordial and has been isolated from the convecting mantle for most of Earth’s history