Effect of the viscoelastic lithosphere on polar wander speed caused by the Late Pleistocene glacial cycles

Previous studies of the wander of the rotation pole associated with the Late Pleistocene glacial cycles indicate that the predicted polar wander speed is sensitive to the density jump at the 670 km discontinuity, the thickness of the elastic lithosphere, and the lower mantle viscosity. In particular, the M1 mode related to the density jump at 670 km depth has been shown to contribute a dominant portion of predicted polar wander speed for sufficiently small lower mantle viscosities. In this study, we examine the sensitivity of polar wander to variations in the viscosity of the viscoelastic lithosphere using simplified compressible Maxwell viscoelastic earth models. Model calculations for earth models with a viscoelastic lithosphere of finite viscosity indicate that the contribution of the M1 mode is similar to those associated with the density discontinuity at the core–mantle boundary (C0 mode) and the lithosphere (L0 mode). We speculate that this is due to the interaction between the M1 mode and the transient mode associated with the viscoelastic lithosphere, which reduces the magnitude of polar wander rates. Therefore, the M1 mode does not contribute a dominant portion of the predicted polar wander speed for earth models with a viscoelastic lithosphere of finite viscosity. In this case, predictions of polar wander speed as a function of lower mantle viscosity exhibit the qualitative form of an 'inverted parabola', as predicted for the J ˙₂ curve. We caution, however, that these results are obtained for simplified earth models, and the results for seismological earth models such as PREM may be complicated by the interaction between the M1 mode and the large set of transient modes.     

Free-surface formulation of mantle convection—II. Implication for subduction-zone observables

Viscous and viscoelastic models for a subduction zone with a faulted lithosphere and internal buoyancy can self-consistently and simultaneously predict long-wavelength geoid highs over slabs, short-wavelength gravity lows over trenches, trench-forebulge morphology, and explain the high apparent strength of oceanic lithosphere in trench environments. The models use two different free-surface formulations of buoyancy-driven flows (see, for example, Part I): Lagrangian viscoelastic and pseudo-free-surface viscous formulations. The lower mantle must be stronger than the upper in order to obtain geoid highs at long wavelengths. Trenches are a simple consequence of the negative buoyancy of slabs and a large thrust fault, decoupling the overriding from underthrusting plates. The lower oceanic lithosphere must have a viscosity of less than to²⁴ Pa s in order to be consistent with the flexural wavelength of forebulges. Forebulges are dynamically maintained by viscous flow in the lower lithosphere and mantle, and give rise to apparently stiffer oceanic lithosphere at trenches. With purely viscous models using a pseudo-free-surface formulation, we find that viscous relaxation of oceanic lithosphere, in the presence of rapid trench rollback, leads to wider and shallower back-arc basins when compared to cases without viscous relaxation. Moreover, in agreement with earlier studies, the stresses necessary to generate forebulges are small (∼ 100 bars) compared to the unrealistically high stresses needed in classic thin elastic plate models.     

Apparent polar wander of the mean-lithosphere reference frame

Apparent polar wander in the mean-lithosphere (= no-net-rotation = no-net-torque uniform drag) reference frame is compared with apparent polar wander in the hotspot reference frame over the past 100 Myr. Palaeo-magnetic poles and plate rotations previously used to determine an apparent polar wander path for the hotspot reference frame are here used to determine an apparent polar wander path in the mean-lithosphere reference frame. We find that the two paths are similar, especially for Late Cretaceous time, when a 10°–20° shift of the pole occurred. To first-order the hotspots and lithosphere (as a whole) moved in unison relative to the palaeomagnetic axis during Late Cretaceous time. A non-dipole field explanation for the apparent shift can probably be excluded. However, either motion of the time-averaged geomagnetic axis relative to the spin axis or polar wandering could have caused this shift, the latter being the more likely explanation.