Mantle Rheology Research Articles

Quantitative evaluation of mantle flow traction on overlying tectonic plate: linear versus power-law mantle rheology

SUMMARY Sub-plate mantle flow traction (MFT) has been considered as a major driving force for plate motion; however, the force acting on the overlying plate is difficult to constrain. One of the reasons lies in the variable rheological flow laws of mantle rocks, for example, linear versus power-law rheology, applied in previous studies. Here, systematic numerical models are conducted to evaluate MFT under variable rheological, geometrical and kinematic conditions. The results indicate that MFT with power-law rheology is much lower than that with linear rheology under the same mantle/plate velocity contrast. In addition, existence of a lithospheric root in the overlying plate could enhance MFT, where integrated normal force acting on the walls of lithospheric root is much lower than the shear force in a large-scale domain. In the acting domain of several thousand kilometres, MFT with power-law rheology is comparable to the ridge push of about 3×1012 N m−1, whereas that with linear rheology is comparable to the slab pull of about 3×1013 N m−1. The roles of MFT in driving plate motion are further analysed for the Tethyan evolution. It indicates that MFT with power-law rheology could partially support the Wilson cycles experienced in the Tethyan system, whereas that with linear rheology could easily dominate any kinds of plate tectonic evolutions. The quantitative evaluation of MFT in this study clarifies the roles of rheological flow laws on MFT and could help us to better understand the contrasting results in previous numerical studies.

Constraining Earth's mantle rheology with Love and Shida numbers at the [formula omitted] tidal frequency

We use measurements of Earth's tidal response at the M2 frequency to constrain the rheology of its mantle. The viscoelasticity and anelasticity of the planet are modeled with an Andrade rheology that depends on two parameters: α and ζ. In this paper, we propose an improved algorithm to compute Earth's tidal deformation. Its Love and Shida numbers k2, h2, k3 and l2 as well as the tidal lag ϵ were calculated for two viscosity profiles and for a wide range of values of α and ζ. By comparing our results with geodetic measurements we obtain the range of values of α and ζ that successfully describes Earth's viscoelastic behavior. Values of ζ as high as 105 can not be excluded. For ζ=1, α should be in the range from 0.19 to 0.33, while for a ζ=105, α is most likely between 0.11 and 0.17. We believe that a similar rheology should be used in geophysical models of other rocky planets and satellites. The obtained results are mostly representative of the lower mantle.We have also shown that for several combinations of the two parameters α and ζ we could obtain nearly identical values of Earth's ℜk2, ℜh2, ℜl2 and ℜk3 with considerably different values of the associated tidal lag. This shows that the approach of always setting ζ=1 might be too simplistic and an Andrade rheology with two free parameters is needed to constrain both the real and imaginary parts of Love and Shida numbers.

Rheological Structure and Lithospheric Stress Interaction in the Alaska Subduction Zone Gleaned From the 2018 Mw 7.9 Oceanic Crustal Earthquake

AbstractPostseismic deformation at convergent margins is controlled mainly by continuous slip on the fault (afterslip) and relaxation of the earthquake‐induced stress in the viscoelastic upper mantle (viscoelastic relaxation). Study of these deformation processes provides insight into the rheological properties of upper mantle and slip behavior of the fault. We have constructed a three‐dimensional finite element model to investigate the postseismic deformation of the 2018 Mw 7.9 Kodiak earthquake. We derived the first 2‐year postseismic Global Positioning System observations to constrain afterslip and upper mantle rheology in the south‐central Alaska. The upper mantle is separated into the mantle wedge and oceanic upper mantle topped by an 80‐km thick asthenosphere layer by the subducting slab. Results show that afterslip generally occurred in areas adjacent to the rupture zone and has a small magnitude of a few tens of millimeters. The viscosities of the asthenosphere and mantle wedge are determined to be in a range of 1–4 × 1018 and 0.5–5 × 1019 Pa s with an optimal value of 2 × 1018 and 2 × 1019 Pa s, respectively. Model results reveal a localized weak mantle wedge of ∼1018 Pa s beneath Lower Cook Inlet that may be due to the fluids dehydrated from the slab. Coulomb stress changes show that the earthquake enhanced coseismic and postseismic stress loading of up to 0.9 and 0.1 bar, respectively, on the shallow subduction interface near Kodiak Island, but there is no obvious triggered seismicity, probably due to the low stress status already released by the 1964 Mw 9.2 Alaska earthquake.

Hydration State and Rheologic Stratification of the Lithospheric Mantle Beneath the North Anatolian Fault, Turkey

AbstractWe present constraints on the hydration state and rheology of the lithospheric mantle beneath the North Anatolian fault zone (NAFZ). Peridotite xenoliths from the Biyikali and Çorlu volcanic centers record deformational microstructures consistent with shearing in a lithosphere‐scale transcurrent fault system. Analysis by Fourier transform infrared spectroscopy indicates that nominally anhydrous phases retain some OH−, but bulk rock concentrations are generally restricted to <50 ppm H2O by weight. From the rock microstructure, we determined differential stress magnitude and active deformation mechanism(s); combined with estimates of hydration state, we constrained the rheology. Recrystallized grain size piezometry shows that the mantle beneath the NAFZ sustained differential stresses of 10–20 MPa, largely independent of depth. The dominant deformation mechanism(s) change with depth; xenoliths extracted from shallower depths record evidence for grain size‐sensitive creep possibly in the presence of melt. At intermediate depths, both dislocation creep and grain size‐sensitive mechanisms were active, and we did not observe evidence for deformation in the presence of melt. The deepest samples were dominated by dislocation creep. The strong temperature sensitivity of creep mechanisms, combined with the low variability in differential stress, contributes to a stratified viscosity profile ranging from 1018 Pa s for the deepest samples to >1022 Pa s at shallower depths (assuming a melt‐free rheology). Although difficult to quantify from the rock record, melt likely reduced the viscosity of the shallow lithospheric mantle. Vertical stratification in viscosity beneath the NAFZ, the result of melt‐present deformation and/or transitions in the dominant deformation mechanism, has important consequences for the seismic cycle of strike‐slip fault systems.

Sensitivity of gravity anomalies to mantle rheology at mid-ocean ridge – transform fault systems

Marine gravity data can provide information on the distribution of mass anomalies in the oceanic crust and upper mantle. Computing corresponding gravity anomalies, especially so-called ‘residual’ gravity anomalies that directly reflect variations in the crustal structure, relies on gravity corrections of both seafloor relief and lithospheric thermal structure. The lithospheric thermal gravity correction involves either a plate cooling approximation or a mantle flow model with the latter typically done using simplified assumptions on mantle rheology. However, a detailed study of how differing rheological models affect the computed gravity anomalies is still missing. Here, we systematically examine the differences in residual mantle Bouguer anomalies (RMBA) caused by differing assumptions on mantle rheology for 16 mid-ocean ridge – transform fault systems. Our calculations show that isoviscous models tend to underpredict RMBA values within the transform deformation zone and overpredict them in the far field at older plate ages, when compared to plate cooling and nonlinear viscoplastic models. This discrepancy stems from isoviscous models failing to capture plate-like deformation, as well as their inability to resolve brittle failure and the associated strain localization that leads to warm upwelling beneath the transform fault. By exploring a wide parameter range, we find that the importance of mantle rheology scales with plate tectonic parameters at the mid-ocean ridge – transform fault system such as transform age offset, spreading rate, and transform fault length. These findings suggest that gravity thermal corrections at the intrinsically three-dimensional ridge – transform systems should employ mantle flow models that resolve plate-like deformation and brittle failure.

Post-seismic gravity change modelling based on non-linear power-law upper mantle rheology

SUMMARY Post-seismic gravity change modelling is commonly based on earth model with Newtonian linear rheology. Here, we present a novel way of modelling post-seismic gravity change by using a non-linear power-law rheology earth model. The method is constructed based on the framework of spectral finite element method (SFEM). SFEM has been proven practical for the purpose of modelling gravity change occurring during megathrust earthquakes. Our method implements the strain rate expression of non-linear power-law rheology into the mathematical framework of SFEM. Using our method, simulations of geoid change caused by synthetic point source earthquakes were made. The results revealed the potential of using non-linear power-law rheology for the explanation of rapid gravity changes in the beginning of the post-seismic epoch. Further on, we computed the post-seismic gravity change of the 2011 M9.0 Tohoku-Oki earthquake based on Maxwell non-linear power-law rheology in the upper mantle. It demonstrated the potential of power-law upper mantle flow as a possible candidate for the explanation of post-seismic gravity change after the earthquake occurs. However, a Maxwell non-linear power-law rheology alone is insufficient for the explanation of the post-seismic gravity change occurred.

Selective inversion of rift basins in lithospheric-scale analogue experiments

Abstract. Basin inversion is commonly attributed to the reverse reactivation of basin-bounding normal faults. This association implies that basin uplift and inversion-related structures are mainly controlled by the frictional behaviour of pre-existing faults and associated damage zones. In this study, we use lithospheric-scale analogue experiments of orthogonal extension followed by shortening to explore how the flow behaviour of ductile layers underneath rift basins promote or suppress basin inversion. Our experiments show that the rheology of the ductile lower crust and lithospheric mantle, modulated by the imposed bulk strain rate, determine (1) basin distribution in a wide rift setting and (2) strain accommodation by fault reactivation and basin uplift during subsequent shortening. When the ductile layers deform uniformly during extension (i.e. stretching) and shortening (i.e. thickening), all of the basins are inverted. When deformation in the ductile layers is localised during extension (i.e. necking) and shortening (i.e. folding), only some basins – which are evenly spaced apart – are inverted. We interpret the latter as selective basin inversion, which may be related to the superposition of crustal-scale and lithospheric-scale boudinage during the previous basin-forming extensional phase and/or folding of the ductile layers during shortening.

Effects of Composite Rheology on Plate‐Like Behavior in Global‐Scale Mantle Convection

AbstractEarth's upper mantle rheology controls lithosphere‐asthenosphere coupling and thus surface tectonics. Rock deformation experiments and seismic anisotropy measurements indicate that composite rheology (co‐existing diffusion and dislocation creep) occurs in the Earth's uppermost mantle, potentially affecting convection and surface tectonics. Here, we investigate how the spatio‐temporal distribution of dislocation creep in an otherwise diffusion‐creep‐controlled mantle impacts the planform of convection and the planetary tectonic regime as a function of the lithospheric yield strength in numerical models of mantle convection self‐generating plate‐like tectonics. The low upper‐mantle viscosities caused by zones of substantial dislocation creep produce contrasting effects on surface dynamics. For strong lithosphere (yield strength > 35 MPa), the large lithosphere‐asthenosphere viscosity contrasts promote stagnant‐lid convection. In contrast, the increase of upper mantle convective vigor enhances plate mobility for lithospheric strength <35 MPa. For the here‐used model assumptions, composite rheology does not facilitate the onset of plate‐like behavior at large lithospheric strength.

Grain size reduction by plug flow in the wet oceanic upper mantle explains the asthenosphere's low seismic Q zone

The prominent seismic low-velocity zone (LVZ) in the oceanic low-viscosity asthenosphere is approximately coincident with the low seismic Q zone (LQZ; a zone of high seismic attenuation). Small olivine grain sizes could link these seismic and rheological properties because they reduce viscosity, seismic velocity and seismic Q. Because rock deformation reduces grain sizes, the asthenosphere's seismic properties can place constraints on asthenospheric flow. To determine dominant flow patterns, we develop a self-consistent analytical 1-D channel flow model that accounts for upper mantle rheology and its dependence on flow-modified grain-sizes, water content, and melt fraction, for flow driven by both surface plate motions (Couette flow) and/or horizontal pressure gradients (Poiseuille flow). We find that Couette flow dominates if the upper mantle is dry, and plug flow (a Poiseuille flow for power law rheology) dominates if it is weakened by wet conditions. A plug flow configuration spanning the upper 670 km of the mantle best explains the LQZ in the asthenosphere, and can be attributed to significant grain-size reduction due to extensive shearing across the asthenosphere. This flow configuration also explains high seismic Q values below the asthenosphere associated with minimal shear deformation and large grain sizes above the mantle transition zone. We suggest that the asthenospheric LQZ and LVZ seismic anomalies can be largely explained by grain-size variations associated with plug flow in the wet upper mantle.

Probabilistic Assessment of Antarctic Thermomechanical Structure: Impacts on Ice Sheet Stability

AbstractUncertainty in present‐day glacial isostatic adjustment (GIA) rates represents at least 44% of the total gravity‐based ice mass balance signal over Antarctica. Meanwhile, physical couplings between solid Earth, sea level and ice dynamics enhance the dependency of the spatiotemporally varying GIA signal on three‐dimensional variations in mantle rheology. Improved knowledge of thermomechanical mantle structure is therefore required to refine estimates of current and projected ice mass balance. Here, we present a Bayesian inverse method for self‐consistently mapping shear‐wave velocities from high‐resolution adjoint tomography into thermomechanical structure using calibrated parameterisations of anelasticity at seismic frequency. We constrain the model using regional geophysical data sets containing information on upper mantle temperature, attenuation and viscosity structure. Our treatment allows formal quantification of parameter covariances, and naturally permits propagation of material parameter uncertainties into thermomechanical structure estimates. We find that uncertainty in steady‐state viscosity structure at 150 km depth can be reduced by 4–5 orders of magnitude compared with a forward‐modeling approach neglecting covariance between viscoelastic parameters. By accounting for the dependence of apparent viscosity on loading timescale, we find good agreement between our estimates of mantle viscosity beneath West Antarctica, and those derived from satellite GPS. Direct access to temperature structure allows us to estimate lateral variations in lithosphere‐asthenosphere boundary (LAB) depth, geothermal heat flow (GHF), and associated uncertainties. We find evidence for shallow LAB depths (63 ± 13 km), and high GHF (76 ± 7 mW m−2) beneath West Antarctica that, combined with low asthenospheric viscosities, indicate a highly dynamic response to ice mass loss.

Build-up of multiple volatiles in Earth's continental keels: Implications for craton stability

Elevated contents of structurally-bound hydrogen in olivine, introduced to cratonic mantle during metasomatism, are widely believed to decrease the effective viscosity so that it is similar to the underlying asthenosphere (1018–1019 Pa⋅s) and susceptible to basal erosion. Nevertheless, large variations exist in H contents of mantle olivines from the Kaapvaal, Siberia, Slave, Tanzania and Wyoming cratons (0 to 320 ppmw) and the role of water as the main driver of rheological weakening is contentious. While recent experimental studies have shown that olivine has the capacity to host large quantities of fluorine, as well as hydrogen, the magnitude of this has not yet been established for cratonic mantle.Our new dataset for Kaapvaal craton peridotites shows that the style of metasomatism influences the addition of H2O and F to cratonic mantle. Silicic fluids derived from subducted slabs, and responsible for the pervasive orthopyroxene enrichment observed in the Kaapvaal (and many of the world's other cratons), deliver significant quantities of H2O but lesser amounts of F, whereas proto-kimberlite melts transport high quantities of both H2O and F. Kaapvaal mantle olivines are major hosts of both H2O (up to 105 ppmw) and F (up to 180 ppm), and we propose that a major increases in the bulk H2O (∼30%) and F (∼65%) of the Kaapvaal craton occurred over a short (20 Ma) time interval between the two main pulses of Cretaceous kimberlite emplacement in southern Africa.By combining the thermal structure of individual cratons with corresponding H2O data for olivine, we show that the effective viscosity contrast at the base of the mechanical boundary layer and asthenosphere below the Kaapvaal, Slave and Siberia cratons varies from 2 to 14. Some of the lowest effective viscosities (3.5×1017 Pa⋅s) occur in H2O-rich olivine in Siberia sheared peridotites from the base of the craton and are consistent with highly-localised metasomatism. More viscous peridotites (up to 1.5×1021 Pa⋅s) were entrained from the thermal boundary layer beneath the Kaapvaal and Tanzania cratons but residence times in this region are short due to convective overturn (<100 Ma) and this dehydrated mantle would offer only temporary protection to cratonic keels. Our viscosity estimates are based on the H content of olivine and, while the additional effect of structurally bound F in olivine is currently uncertain, we anticipate that it will have a similar effect on mantle rheology to H, so that the base of the mechanical boundary layer beneath cratons at the time and location of kimberlite generation is close to the tipping point for instability.The rapid build-up of volatiles associated with pervasive kimberlite activity may have been the catalyst for lithospheric thinning on the southern margin of the Kaapvaal craton, but more localised pulses of kimberlite activity occurred over longer time intervals in the Siberia and Slave cratons and so had a less profound effect on their stability. A prolonged and widespread subduction-flux of H2O to the Wyoming and N. China cratons, and subsequent sub-solidus partitioning in olivine, may have been the driver for the rheological weakening that ultimately led to delamination.

Volatiles and Intraplate Magmatism: a Variable Role for Carbonated and Altered Oceanic Lithosphere in Ocean Island Basalt Formation

Abstract Recycling of material at subduction zones has fundamental implications for melt composition and mantle rheology. Ocean island basalts (OIBs) sample parts of the mantle from variable depths that have been diversely affected by subduction zone processes and materials, including the subducted slab, metasomatising melts and fluids. Resultant geochemical differences are preserved at a variety of scales from melt inclusions to whole rocks, from individual islands to chains of islands. Here we examine a global dataset of ocean island basalt compositions with a view to understanding the connection between silica-saturation, olivine compositions, and halogens in glass and olivine-hosted melt inclusions to reveal information regarding the mantle sources of intraplate magmatism. We find that minor elements incorporated into olivine, although informative, cannot unambiguously discriminate between different source contributions, but indicate that none of the OIB analysed here are derived solely from dry peridotite melting. Nor can differences in lithospheric thickness explain trace element variability in olivine between different ocean islands. We present new halogen (F, Cl, Br/Cl, I/Cl) data along with incompatible trace element data for the global array and encourage measurement of fluorine along with heavier halogens to obtain better insight into halogen cycling. We suggest that Ti-rich silica-undersaturated melts require a contribution from carbonated lithosphere, either peridotite or eclogite and are an important component sampled by ocean island basalts, together with altered oceanic crust. These results provide new insights into our understanding of mantle-scale geochemical cycles, and also lead to the potential for the mantle transition zone as an underestimated source for observed volatile and trace-element enrichment in ocean island basalts.

Electrical Moho and upper mantle metasomatism in the Archaean Bundelkhand craton of the Indian peninsular shield

The Bundelkhand Craton (BKC), one of the Archaean nuclei in the north central part of the Indian shield, has a lesser surface exposure compared to other cratons. Recent magnetotelluric (MT) study in the Central Ganga Basin (CGB) showed NNE-SSW extension of the cratonic block as a transverse ridge and delineated a shear zone within the buried section of the BKC. In the present study, we extend the MT profile to further south with 19 broadband MT sites in the exposed BKC to complement the lithospheric electrical resistivity architecture. Two-dimensional (2-D) geoelectric structure down to 100 km is obtained from inverse modelling of the distortion corrected MT impedance tensors together with the 31 sites in the CGB. The resulting geoelectric model reveals a flat and sharp resistivity contrast at 40–42 km depth, the electrical Moho (e-Moho), separating highly resistive Bundelkhand massif and moderately conductive upper mantle. The 2-D model also brings out two shear zones represented by near-vertical moderately conductive zones within the massif, indicating zones of amalgamation of different fragments of the pre-Bundelkhand craton. In contrast to other major cratons of the Indian shield, the present MT study has delineated a moderately conductive upper mantle beneath the BKC. In view of the estimated low mantle heat flow along with relatively flat topography, we infer metasomatism of the mantle lithosphere by percolation of fluids. We hypothesize that infiltration of fraction of water (H2O) in nominally anhydrous mantle minerals as hydrogen (H) bonded to structural oxygen (O) originated from the subducted oceanic slab has altered the mantle rheology and metasomatized the lithosphere.

The history of global strain and geodynamics on Mars

The tectonic record of Mars is dominated by compressional wrinkle ridges on volcanic surfaces, and these structures have been widely used as a record of tectonic and geodynamic evolution. This study analyzes the density of compressional tectonic structures and inferred strain as functions of time, using the lengths and heights of compressional ridges together with geologic estimates of surface age. Our analyses confirm an apparent peak in compressional strain in the late Noachian and early Hesperian, and comparatively lower values before and after. The lower tectonic strain in the early and middle Noachian relative to the early Hesperian reflects the incompleteness of the ancient compressional tectonic record of strain, as older surfaces should accumulate more strain than younger surfaces. This strain deficit necessitates other means of accommodating contractional strain in the ancient crust, including distributed strain and small-scale faulting. The abrupt decrease in the accumulated tectonic strain and strain rate after the early Hesperian reflects a rapid episode of global contraction followed by much lower rates, in conflict with models of steady-state mantle evolution. The decreasing strain rate may be explained by a changing mantle rheology due to volcanic outgassing. Alternatively, the strain history may be explained by the dominance of mantle plumes in the late Noachian and early Hesperian associated with the formation of Tharsis and the early Hesperian volcanic provinces, which could have led to a pulse of rapid contraction caused by disequilibrium cooling and volcanic outpourings. The tectonic record of strain on Mars – including the incomplete ancient record and evidence for strong departures from steady-state models – may have implications for the interpretation of the tectonic record and thermal evolution of other bodies as well.

Periclase deforms more slowly than bridgmanite under mantle conditions

Transport of heat from the interior of the Earth drives convection in the mantle, which involves the deformation of solid rocks over billions of years. The lower mantle of the Earth is mostly composed of iron-bearing bridgmanite MgSiO3 and approximately 25% volume periclase MgO (also with some iron). It is commonly accepted that ferropericlase is weaker than bridgmanite1. Considerable progress has been made in recent years to study assemblages representative of the lower mantle under the relevant pressure and temperature conditions2,3. However, the natural strain rates are 8 to 10 orders of magnitude lower than in the laboratory, and are still inaccessible to us. Once the deformation mechanisms of rocks and their constituent minerals have been identified, it is possible to overcome this limitation thanks to multiscale numerical modelling, and to determine rheological properties for inaccessible strain rates. In this work we use 2.5-dimensional dislocation dynamics to model the low-stress creep of MgO periclase at lower mantle pressures and temperatures. We show that periclase deforms very slowly under these conditions, in particular, much more slowly than bridgmanite deforming by pure climb creep. This is due to slow diffusion of oxygen in periclase under pressure. In the assemblage, this secondary phase hardly participates in the deformation, so that the rheology of the lower mantle is very well described by that of bridgmanite. Our results show that drastic changes in deformation mechanisms can occur as a function of the strain rate.

On the Choice and Implications of Rheologies That Maintain Kinematic and Dynamic Consistency Over the Entire Earthquake Cycle

AbstractViscoelastic processes in the upper mantle redistribute seismically generated stresses and modulate crustal deformation throughout the earthquake cycle. Geodetic observations of these motions at the surface of the crust‐mantle system offer the possibility of constraining the rheology of the upper mantle. Parsimonious representations of viscoelastically modulated deformation through the aseismic phase of the earthquake cycle should simultaneously explain geodetic observations of (a) rapid postseismic deformation, (b) late in the earthquake cycle near‐fault strain localization. To understand how rheological formulations affect kinematics, we compare predictions from time‐dependent forward models of deformation over the entire earthquake cycle for an idealized vertical strike‐slip fault in a homogeneous elastic crust underlain by a homogeneous viscoelastic upper‐mantle. We explore three different rheologies as inferred from laboratory experiments: (a) linear Maxwell, (b) linear Burgers, (c) power‐law. The linear Burgers and power‐law rheologies are consistent with fast and slow deformation phenomenology over the entire earthquake cycle, while the single‐layer linear Maxwell model is not. The kinematic similarity of linear Burgers and power‐law models suggests that geodetic observations alone may be insufficient to distinguish between them, but indicate that one may serve as an effective proxy for the other. However, the power‐law rheology model displays a postseismic response that is non‐linearly dependent on earthquake magnitude, which may offer a partial explanation for observations of limited postseismic deformation near some magnitude 6.5–7.0 earthquakes. We discuss the role of mechanical coupling between frictional slip and viscous creep in controlling the time‐dependence of regional stress transfer following large earthquakes and how this may affect the seismic hazard and risk to communities living close to fault networks.

Kinematic Representations of Linear and Power‐Law Viscoelastic Deformation Through the Earthquake Cycle

AbstractViscoelastic deformation below the Earth's elastic crust is modulated by stresses generated by both plate tectonic and earthquake cycle processes. Rapid near‐fault deformation following large earthquakes has been interpreted as the signature of viscoelastic stress diffusion in the upper mantle and earthquake cycle models have been developed that integrated this effect throughout the duration of the earthquake cycle. Here we develop a surrogate method to approximate the upper crustal kinematics associated with time‐dependent viscoelastic deformation that does not require knowledge of nor provide constraints on subcrustal rheology. To do this we show how the effects of deformation in the viscoelastic upper mantle can be emulated by introducing a set of effective dislocations at crust‐mantle interface. The approach is shown to approximate both linear and power‐law viscoelastic rheologies and offers a way to accurately represent viscoelastic kinematics even in the case where upper mantle rheology is poorly constrained.

Simultaneous inference of plate boundary stresses and mantle rheology using adjoints: large-scale 2-D models

SUMMARY Plate motions are a primary surface constraint on plate and mantle dynamics and rheology, plate boundary stresses and the occurrence of great earthquakes. Within an optimization method, we use plate motion data to better constrain uncertain mantle parameters. For the optimization problem characterizing the maximum a posteriori rheological parameters we derive gradients using adjoints and expressions to approximate the posterior distributions for stresses within plate boundaries. We apply these methods to a 2-D cross section from the western to eastern Pacific, with temperature distributions and fault zone geometries developed primarily from seismic and plate motion data. We find that the best-fitting stress exponent, n, is about 2.8 and the yield stress about 100 MPa or less. The normal stress on the interplate fault zones is about 100 MPa and the shear stresses about 10 MPa or less.

A Seismic Tomography, Gravity, and Flexure Study of the Crust and Upper Mantle Structure of the Emperor Seamounts at Jimmu Guyot

AbstractThe intraplate Hawaiian‐Emperor Seamount Chain has long been considered a hotspot track generated by the motion of the Pacific plate over a deep mantle plume, and an ideal feature therefore for studies of volcanic structure, magma supply, plume‐crust interaction, flexural loading, and upper mantle rheology. Despite their importance as a major component of the chain, the Emperor Seamounts have been relatively little studied. In this paper, we present the results of an active‐source wide‐angle reflection and refraction experiment conducted along an ocean‐bottom‐seismograph (OBS) line oriented perpendicular to the seamount chain, crossing Jimmu guyot. The tomographic P wave velocity model, using ∼20,000 travel times from 26 OBSs, suggests that there is a high‐velocity (&gt;6.0 km/s) intrusive core within the edifice, and the extrusive‐to‐intrusive ratio is estimated to be ∼2.5, indicating that Jimmu was built mainly by extrusive processes. The total volume for magmatic material above the top of the oceanic crust is ∼5.3 × 104 km3, and the related volume flux is ∼0.96 m3/s during the formation of Jimmu. Under volcanic loading, the ∼5.3‐km‐thick oceanic crust is depressed by ∼3.8 km over a broad region. Using the standard relationships between Vp and density, the velocity model is verified by gravity modeling, and plate flexure modeling indicates an effective elastic thickness (Te) of ∼14 km. Finally, we find no evidence for large‐scale magmatic underplating beneath the pre‐existing crust.

Zoisite in cratonic eclogite xenoliths - Implications for water in the upper mantle

A significant portion of the Earth's water is found at trace levels in nominally anhydrous rock forming minerals in the crust and mantle. Its distribution has far-reaching implications in melt-generating processes, mantle rheology and the stability of the lithosphere. We report petrographic, major and trace element, Fourier transform infrared spectroscopy and secondary ion mass spectrometry data on four eclogite xenoliths from the Obnazhennaya kimberlite, Siberia that contain hydrous minerals at upper mantle conditions. Al-rich (>9.9 wt% Al2O3) omphacitic clinopyroxene, a major mineral of eclogites, contains crystallographically controlled exsolutions of CaO-rich (up to 23.6 wt%) garnet (coronae, lenses), garnet-zoisite (lamellae) or amphiboles (needle/blade-like exsolutions). The omphacite hosting the exsolved lamellae shows enrichment or depletion in major (Al, Si, Mg) and trace (e.g., HREE) elements compatible with garnet and zoisite, proportional to the abundance of exsolutions.We argue that zoisite likely formed concomitantly with garnet exsolution from a water- and Al-rich precursory omphacite, without any fluid addition. The compositional gradients are consistent with diffusion-controlled partitioning in a closed system, and their preservation indicates the exsolution likely took place shortly prior to eclogite entrapment by the kimberlite. Large, non-exsolved, omphacite grains contain minimum ~870–1500 ppm wt. H2O (determined as OH). We infer that intertwined lamellae of hydrous garnet (600–960 ppm wt. H2O) and minor (<1%) zoisite (2 wt% H2O) formed at mantle conditions due to progressive exsolution, without external fluids. Recalculated whole rock water contents range from ~320 to 970 ppm wt. H2O. These estimates exceed by far those for the surrounding peridotitic mantle and suggest that, although a minor component, eclogites may locally be an important water reservoir in the cratonic mantle.