Evolution of the slab bending radius and the bending dissipation in three-dimensional subduction models with a variable slab to upper mantle viscosity ratio

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

Moho depths and three-dimensional velocity structure of the crust and upper mantle beneath the Baikal region, from local tomography

Russian Geology and Geophysics ◽

10.1016/j.rgg.2007.02.005 ◽

2007 ◽

Vol 48 (2) ◽

pp. 204-220 ◽

Cited By ~ 12

Author(s):

A.V. Yakovlev ◽

I.Yu. Koulakov ◽

S.A. Tychkov

Keyword(s):

Upper Mantle ◽

Baikal Region ◽

Velocity Structure ◽

Three Dimensional ◽

Crust And Upper Mantle ◽

Local Tomography ◽

Moho Depths

Download Full-text

Plumbing and architecture of a crustal mineralizing system in the Archean Superior Province, Canada

10.5194/egusphere-egu21-8885 ◽

2021 ◽

Author(s):

Eric Roots ◽

Graham Hill ◽

Ben M. Frieman ◽

James A. Craven ◽

Richard S. Smith ◽

...

Keyword(s):

Upper Mantle ◽

Spatial Association ◽

Three Dimensional ◽

3D Models ◽

Superior Province ◽

Seismic Reflection Data ◽

Reflection Data ◽

Abitibi Subprovince ◽

Lower Crustal

The role of melts and magmatic/metamorphic fluids in mineralization processes is well established. However, the role of crustal architecture in defining source and sink zones in the middle to lower crust remains enigmatic. Integration of three dimensional magnetotelluric (MT) modelling and seismic reflection data across the Archean Abitibi greenstone belt of the Superior Province, Canada, reveals a &#8216;whole-of-crust&#8217; mineralizing system and highlights the controls by crustal architecture on metallogenetic processes. Electrically conductive conduits in an otherwise resistive upper crust are coincident with truncations and offsets of seismic reflections that are mostly interpreted as major brittle-ductile fault zones. The spatial association between these features and low resistivity zones imaged in the 3D models suggest that these zones acted as pathways through which fluids and melts ascended toward the surface. At mid-crustal levels, these &#8216;conduit&#8217; zones connect to ~50 km long, north-south striking conductors, and are inferred to represent graphite and/or sulphide deposited from cooling fluids. At upper mantle to lower crustal depths, east-west trending conductive zones dominate and display shallow dips. The upper mantle features are broadly coincident with the surface traces of the major deformation zones with which a large proportion of the gold endowment is associated. We suggest that these deep conductors represent interconnected graphitic zones perhaps augmented by sulphides that are relicts from metamorphic fluid and melt emplacement associated primarily with the later stages of regional deformation. &#160;Thus, from the combined MT and seismic data, we develop a crustal-scale architectural model that is consistent with existing geological and deformational models, providing constraints on the sources for and signatures of fluid and magma emplacement that resulted in widespread metallogenesis in the Abitibi Subprovince.

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

Gravity constraints on structure of the East-West Antarctic lithospheric transition zone

10.26686/wgtn.15110442.v1 ◽

2021 ◽

Author(s):

Sam Treweek

Keyword(s):

Transition Zone ◽

Upper Mantle ◽

Geophysical Methods ◽

East Antarctica ◽

Gravity Models ◽

West Antarctica ◽

Asthenospheric Mantle ◽

Surface Uplift ◽

Mantle Viscosity ◽

Taylor Valley

The differing structural evolution of cratonic East Antarctica and younger West Antarctica has resulted in contrasting lithospheric and asthenospheric mantle viscosities between the two regions. Combined with poor constraints on the upper mantle viscosity structure of the continent, estimates of surface uplift in Antarctica predicted from models of glacial isostatic adjustment (GIA) and observed by Global Satellite Navigation System (GNSS) contain large misfits. This thesis presents a gravity study ofthe lithospheric transition zone beneath the Taylor Valley, Antarctica, conducted to constrain the variation in lithological parameters such as viscosity and density of the upper mantle across this region. During this study 119 new gravity observations were collected in the ice-free regions of the Taylor Valley and amalgamated with 154 existing land-based gravity observations, analysed alongside aerogravity measurements of southern Victoria Land. Gravity data are used to construct 2D gravity models of the subsurface beneath this region. An eastward gradient in Bouguer anomalies of ~- 1.6 mGal/km is observed within the Taylor Valley. Models reveal thickening of the Moho from 23±5 km beneath the Ross Sea to 35±5 km in the Polar Plateau (dipping at 24.5±7.2°), and lithospheric mantle 100 km thicker in East Antarctica (~200±30 km) than West Antarctica (~90±30 km). Models of predicted surface uplift history are used to estimate an asthenospheric mantle viscosity of 2.1x1020 Pa.s at full surface recovery beneath the Ross Embayment, differing by ~14% from the viscosity at 50% recovery. The temperature contrast between lithospheric and asthenospheric mantle is estimated as ~400°C, equivalent to a viscosity that decreases by a factor of about 30 over the mantle boundary. Results demonstrate that the history of surface uplift in the study area may be complicated, resulting in observations of uplift, or subsidence, at GNSS stations. Future work should incorporate additional geophysical methods, such as seismicity and electrical resistivity, improving constraints on gravity models. A better understanding of the surface uplift (or subsidence) history in the Transantarctic Mountains is critical, with implications in reducing uncertainty in GIA models.

Download Full-text