The minimum mantle viscosity of an accreting Earth

SUMMARY Previous studies have shown that a low viscosity upper mantle can impact the wavelength of mantle flow and the balance of plate driving to resisting forces. Those studies assumed that mantle viscosity is independent of mantle flow. We explore the potential that mantle flow is not only influenced by viscosity but can also feedback and alter mantle viscosity structure owing to a non-Newtonian upper-mantle rheology. Our results indicate that the average viscosity of the upper mantle, and viscosity variations within it, are affected by the depth to which a non-Newtonian rheology holds. Changes in the wavelength of mantle flow, that occur when upper-mantle viscosity drops below a critical value, alter flow velocities which, in turn, alter mantle viscosity. Those changes also affect flow profiles in the mantle and the degree to which mantle flow drives the motion of a plate analogue above it. Enhanced upper-mantle flow, due to an increasing degree of non-Newtonian behaviour, decreases the ratio of upper- to lower-mantle viscosity. Whole layer mantle convection is maintained but upper- and lower-mantle flow take on different dynamic forms: fast and concentrated upper-mantle flow; slow and diffuse lower-mantle flow. Collectively, mantle viscosity, mantle flow wavelengths, upper- to lower-mantle velocities and the degree to which the mantle can drive plate motions become connected to one another through coupled feedback loops. Under this view of mantle dynamics, depth-variable mantle viscosity is an emergent flow feature that both affects and is affected by the configuration of mantle and plate flow.

Download Full-text

Dislocation interactions in olivine control postseismic creep of the upper mantle

Nature Communications ◽

10.1038/s41467-021-23633-8 ◽

2021 ◽

Vol 12 (1) ◽

Author(s):

David Wallis ◽

Lars N. Hansen ◽

Angus J. Wilkinson ◽

Ricardo A. Lebensohn

Keyword(s):

Upper Mantle ◽

Stress Transfer ◽

Transient Creep ◽

Slip Systems ◽

Stress Distributions ◽

Correlation Lengths ◽

Mantle Viscosity ◽

Dislocation Interactions ◽

Different Temperatures ◽

New Generation

AbstractChanges in stress applied to mantle rocks, such as those imposed by earthquakes, commonly induce a period of transient creep, which is often modelled based on stress transfer among slip systems due to grain interactions. However, recent experiments have demonstrated that the accumulation of stresses among dislocations is the dominant cause of strain hardening in olivine at temperatures ≤600 °C, raising the question of whether the same process contributes to transient creep at higher temperatures. Here, we demonstrate that olivine samples deformed at 25 °C or 1150–1250 °C both preserve stress heterogeneities of ~1 GPa that are imparted by dislocations and have correlation lengths of ~1 μm. The similar stress distributions formed at these different temperatures indicate that accumulation of stresses among dislocations also provides a contribution to transient creep at high temperatures. The results motivate a new generation of models that capture these intragranular processes and may refine predictions of evolving mantle viscosity over the earthquake cycle.

Download Full-text

Mass balance of the Antarctic ice sheet 1992–2016: reconciling results from GRACE gravimetry with ICESat, ERS1/2 and Envisat altimetry

Journal of Glaciology ◽

10.1017/jog.2021.8 ◽

2021 ◽

pp. 1-27

Author(s):

H. Jay Zwally ◽

John W. Robbins ◽

Scott B. Luthcke ◽

Bryant D. Loomis ◽

Frédérique Rémy

Keyword(s):

Mass Balance ◽

Antarctic Peninsula ◽

East Antarctica ◽

Mantle Material ◽

West Antarctica ◽

Mantle Viscosity ◽

Grounding Line ◽

The Antarctic ◽

Total Gain

Abstract GRACE and ICESat Antarctic mass-balance differences are resolved utilizing their dependencies on corrections for changes in mass and volume of the same underlying mantle material forced by ice-loading changes. Modeled gravimetry corrections are 5.22 times altimetry corrections over East Antarctica (EA) and 4.51 times over West Antarctica (WA), with inferred mantle densities 4.75 and 4.11 g cm−3. Derived sensitivities (Sg, Sa) to bedrock motion enable calculation of motion (δB0) needed to equalize GRACE and ICESat mass changes during 2003–08. For EA, δB0 is −2.2 mm a−1 subsidence with mass matching at 150 Gt a−1, inland WA is −3.5 mm a−1 at 66 Gt a−1, and coastal WA is only −0.35 mm a−1 at −95 Gt a−1. WA subsidence is attributed to low mantle viscosity with faster responses to post-LGM deglaciation and to ice growth during Holocene grounding-line readvance. EA subsidence is attributed to Holocene dynamic thickening. With Antarctic Peninsula loss of −26 Gt a−1, the Antarctic total gain is 95 ± 25 Gt a−1 during 2003–08, compared to 144 ± 61 Gt a−1 from ERS1/2 during 1992–2001. Beginning in 2009, large increases in coastal WA dynamic losses overcame long-term EA and inland WA gains bringing Antarctica close to balance at −12 ± 64 Gt a−1 by 2012–16.

Download Full-text

The heat flow constraint on mantle tomography-based convection models: Towards a geodynamically self-consistent inference of mantle viscosity

Journal of Geophysical Research Atmospheres ◽

10.1029/95jb01078 ◽

1995 ◽

Vol 100 (B7) ◽

pp. 12731-12751 ◽

Cited By ~ 42

Author(s):

Giovanni Pari ◽

W. Richard Peltier

Keyword(s):

Heat Flow ◽

Mantle Viscosity ◽

Self Consistent

Download Full-text

A new inference of mantle viscosity based upon joint inversion of convection and glacial isostatic adjustment data

Earth and Planetary Science Letters ◽

10.1016/j.epsl.2004.06.005 ◽

2004 ◽

Vol 225 (1-2) ◽

pp. 177-189 ◽

Cited By ~ 434

Author(s):

J.X. Mitrovica ◽

A.M. Forte

Keyword(s):

Joint Inversion ◽

Glacial Isostatic Adjustment ◽

Isostatic Adjustment ◽

Mantle Viscosity

Download Full-text

Lithosphere and upper-mantle structure of the southern Baltic Sea estimated from modelling relative sea-level data with glacial isostatic adjustment

Solid Earth ◽

10.5194/se-5-447-2014 ◽

2014 ◽

Vol 5 (1) ◽

pp. 447-459 ◽

Cited By ~ 10

Author(s):

H. Steffen ◽

G. Kaufmann ◽

R. Lampe

Keyword(s):

Baltic Sea ◽

Sea Level ◽

Upper Mantle ◽

Glacial Isostatic Adjustment ◽

Relative Sea Level ◽

Southern Baltic Sea ◽

Lithospheric Thickness ◽

Isostatic Adjustment ◽

Mantle Viscosity ◽

Ice Model

Abstract. During the last glacial maximum, a large ice sheet covered Scandinavia, which depressed the earth's surface by several 100 m. In northern central Europe, mass redistribution in the upper mantle led to the development of a peripheral bulge. It has been subsiding since the begin of deglaciation due to the viscoelastic behaviour of the mantle. We analyse relative sea-level (RSL) data of southern Sweden, Denmark, Germany, Poland and Lithuania to determine the lithospheric thickness and radial mantle viscosity structure for distinct regional RSL subsets. We load a 1-D Maxwell-viscoelastic earth model with a global ice-load history model of the last glaciation. We test two commonly used ice histories, RSES from the Australian National University and ICE-5G from the University of Toronto. Our results indicate that the lithospheric thickness varies, depending on the ice model used, between 60 and 160 km. The lowest values are found in the Oslo Graben area and the western German Baltic Sea coast. In between, thickness increases by at least 30 km tracing the Ringkøbing-Fyn High. In Poland and Lithuania, lithospheric thickness reaches up to 160 km. However, the latter values are not well constrained as the confidence regions are large. Upper-mantle viscosity is found to bracket [2–7] × 1020 Pa s when using ICE-5G. Employing RSES much higher values of 2 × 1021 Pa s are obtained for the southern Baltic Sea. Further investigations should evaluate whether this ice-model version and/or the RSL data need revision. We confirm that the lower-mantle viscosity in Fennoscandia can only be poorly resolved. The lithospheric structure inferred from RSES partly supports structural features of regional and global lithosphere models based on thermal or seismological data. While there is agreement in eastern Europe and southwest Sweden, the structure in an area from south of Norway to northern Germany shows large discrepancies for two of the tested lithosphere models. The lithospheric thickness as determined with ICE-5G does not agree with the lithosphere models. Hence, more investigations have to be undertaken to sufficiently determine structures such as the Ringkøbing-Fyn High as seen with seismics with the help of glacial isostatic adjustment modelling.

Download Full-text