The role of thermal effect on mantle seismic anomalies from observations of GIA

NASA Astrophysics Data System (ADS)

Wu, P.; Wang, H. S.; Steffen, H.

2012-04-01

Recent advance in seismic tomography reveals the structure inside the mantle. An outstanding issue is the role of thermal versus non-thermal (e.g. compositional, partial melting) contribution to seismic velocity anomalies. Here we use observations of Glacial Isostatic Adjustment (GIA), e.g. global relative sea level data, GRACE observations (with recent hydrology contributions removed) and GPS crustal uplift rates in combination with 3D GIA models to address this issue. Both ICE-4G and ICE-5G models are tested, but ICE-4G gives much better overall fit to these observations. Also, several 1-D background viscosity profiles, with different viscosity contrast at 670 km depth have also been tested and the one that gives consistent results is model RF3 which has a moderate viscosity increase across 670 km. Lateral mantle viscosity variation is inferred from Ekstrom & Dziewonski's S20A seismic tomography model using a scaling law that includes both the effect of anharmonicity and anelasticity. Thermal contribution to seismic tomography appears as the beta factor in the scaling law. The values of beta in the upper mantle, shallow part of the lower mantle and the deep part of the lower mantle are allowed to be different and the solution space of the beta values is searched to find the best combination that gives the best fit to the GIA observations simultaneously. The result of our best model (RF3 with lateral heterogeneity) shows that thermal effect increases from about 65% in the upper mantle to 80% in the shallow part of the lower mantle and to about 100% in the deep lower mantle above the D" layer. This is consistent with temperature excess in the lower mantle from high core heating. However, the uncertainty increases from < 1% in the upper mantle to 20% in the shallow lower mantle and is not very well constrained in the deep lower mantle.

Quantifying mixing and age variations of heterogeneities in models of mantle convection: Role of depth-dependent viscosity

NASA Astrophysics Data System (ADS)

Hunt, D. L.; Kellogg, L. H.

2001-04-01

Using a two-dimensional finite element model of mantle convection containing over a million tracer particles, we examine the effects of depth-dependent viscosity on the rates and patterns of mixing. We simulate the processes of recycling crust at subduction zones and the homogenization of recycled material (by dispersion and by melting at mid-ocean ridges). Particles are continually introduced at downwellings and destroyed when they either are so thoroughly dispersed that it would be impossible to measure their presence in the geochemical signature of mid-ocean ridges or oceanic islands, or when they are close to spreading centers, at which point melting would "reset" the geochemical clock. A large number of factors influence the flow pattern and thus the rate at which heterogeneities are dispersed by convection. We examine the effect of increasing the viscosity with depth, and determine how both the residence time of heterogeneities and the extent of lateral mixing and exchange between the upper and lower mantle vary with the viscosity profile of the mantle. We determine the particle distribution resulting from convection models with three viscosity profiles: uniform viscosity, a smooth increase of viscosity with depth, and an abrupt jump in viscosity between the upper and lower mantle. We characterize the resulting distribution of heterogeneities in space and time by examining the age distribution of particles and their locations relative to others introduced into the flow at separate downwellings. Mixing rates in the three models are calculated as a function of the number of particles removed from the flow through time. We found that an increase of viscosity at depth does not induce age stratification in which older particles stagnate in the lover mantle, and it does not produce an upper layer (the source of mid-ocean ridge basalt) that is well-mixed compared to the deeper regions. However, pronounced lateral heterogeneity is evident in the distribution of particles of different ages and starting locations that is not apparent from the particle positions alone.

Using GIA observables to constrain the thermal contribution to lateral variations in mantle viscosity

NASA Astrophysics Data System (ADS)

Wu, P.; Wang, H.; van der Wal, W.; Shum, C.; Lee, H.; Braun, A.

2009-05-01

Lateral heterogeneities in the mantle can be caused by thermal, chemical and non-isotropic pre-stress effects. Here, observations of the glacial isostatic adjustment (GIA) process are used to constrain the thermal contribution to lateral variations in upper and lower mantle viscosities. The Coupled Laplace-Finite Element method is used to predict the GIA response on a spherical, self-gravitating, compressible, viscoelastic earth with self-gravitating oceans, induced by either the ICE-5G or ICE-4G deglaciation models. GIA observations include global historic relative sea level data, GPS uplift rates in Laurentide and Fennoscandia, altimetry together with tide-gauge data in the Great Lakes area, and GRACE data in Laurentide. The lateral viscosity perturbations are inferred from the high resolution seismic tomography model of Grand (2002) by using a conversion relation that takes into account both anelastic and anharmonic effects (Karato 2008). To determine the contribution of thermal effects in the upper and lower mantle, the scaling factor b is also inserted into the conversion relation: For b = 1, lateral velocity variations are caused by thermal effects alone; while b < 1 indicates a decreasing contribution of thermal effects; eventually when b = 0, there is no lateral viscosity variations exist and the Earth is laterally homogeneous. The value of b in the upper mantle is b1 while that in the lower mantle is b2. The lateral viscosity variations computed this way are superposed on a reference model that is able to give a reasonably good fit to the GIA observations. The parameter space for (b1, b2) is then searched to find the combination that yields the best improvement in fitting the GIA data in Laurentide, Fennoscandia or globally.

Can lower mantle slab-like seismic anomalies be explained by thermal coupling between the upper and lower mantles?

NASA Astrophysics Data System (ADS)

Čížková, Hana; Čadek, Ondřej; van den Berg, Arie P.; Vlaar, Nicolaas J.

Below subduction zones, high resolution seismic tomographic models resolve fast anomalies that often extend into the deep lower mantle. These anomalies are generally interpreted as slabs penetrating through the 660-km seismic discontinuity, evidence in support of whole-mantle convection. However, thermal coupling between two flow systems separated by an impermeable interface might provide an alternative explanation of the tomographic results. We have tested this hypothesis within the context of an axisymmetric model of mantle convection in which an impermeable boundary is imposed at a depth of 660 km. When an increase in viscosity alone is imposed across the impermeable interface, our results demonstrate the dominant role of mechanical coupling between shells, producing lower mantle upwellings (downwellings) below upper mantle downwellings (upwellings). However, we find that the effect of mechanical coupling can be significantly weakened if a narrow low viscosity zone exists beneath the 660-km discontinuity. In such a case, both thermally induced ‘slabs’ in the lower mantle and thermally activated plumes that rise from the upper/lower mantle boundary are observed even though mass transfer between the shells does not exist.

Upper Mantle Anisotropy Under Fast Spreading Mid-ocean Ridges: 2-D Whole Mantle Convection Model With Subduction

NASA Astrophysics Data System (ADS)

Lee, C.; Zhou, Y.; King, S. D.

2008-12-01

Analyses of seismic anisotropy caused by spatial alignments of anisotropic minerals (e.g., olivine) have been widely used to infer mantle flow directions in the upper mantle. Deep seismic anisotropy beneath fast spreading mid-ocean ridges (e.g., East Pacific Rise) has been recently observed at depths of 200-300 km and even down to the transition zone, with polarization changes in radial anisotropy from VSH < VSV (shallow) to VSH < VSV (deep). We investigate the origin of the observed deep seismic anisotropy and polarization changes beneath the EPR in 2-D Cartesian numerical models using both kinematically (prescribed velocity) and dynamically (negative buoyancy) driven ridge spreading. Because subduction is thought to be an important controlling factor in the style of ridge spreading and mantle convection, we consider a subduction zone developing at the prescribed weak zone. A whole mantle domain expressed by a one by four box (2890 by 11560 km) is used to minimize the boundary effects on the subducting slab. For the upper mantle rheology, we consider composite viscosity of diffusion and dislocation creep for dry olivine to evaluate the effects of lateral variation of mantle viscosity and the rheological changes from dislocation to diffusion creep under the mid-ocean ridge. For the lower mantle rheology, we use diffusion creep for dry olivine by increasing grain size to match relevant lower mantle viscosity. We also consider the 660 km phase transition with density and viscosity jump as well as Clapeyron slope. Anisotropy is evaluated using finite-strain ellipses based on the assumption that a-axes of olivine crystals are parallel to the major axes of the finite-strain ellipses. Our preliminary results show 1) in general, the development of VSH < VSV anisotropy is confined only in a narrow region under the ridge axis at depths of 200- 300 km; 2) strong VSH > VSV anisotropy can be found in the 'asthenosphere' beneath the entire spreading oceanic lithosphere; and 3) the dominate creep mechanism changes from dislocation creep to diffusion creep at depths of 300-400 km; indicating a more isotropic lower upper mantle. We conclude that our geodynamical modeling in a passive ridge spreading system does not produce the deep seismic anisotropy recently observed beneath the EPR. However, we do not consider partial melting, dynamic recrystallization and anisotropic viscosity which would change seismic interpretation and mantle flow, and thus further study is required.

On the resolution of shallow mantle viscosity structure using post-earthquake relaxation data: Application to the 1999 Hector Mine, California, earthquake

USGS Publications Warehouse

Pollitz, Fred F.; Thatcher, Wayne R.

2010-01-01

Most models of lower crust/mantle viscosity inferred from postearthquake relaxation assume one or two uniform-viscosity layers. A few existing models possess apparently significant radially variable viscosity structure in the shallow mantle (e.g., the upper 200 km), but the resolution of such variations is not clear. We use a geophysical inverse procedure to address the resolving power of inferred shallow mantle viscosity structure using postearthquake relaxation data. We apply this methodology to 9 years of GPS-constrained crustal motions after the 16 October 1999 M = 7.1 Hector Mine earthquake. After application of a differencing method to isolate the postearthquake signal from the “background” crustal velocity field, we find that surface velocities diminish from ∼20 mm/yr in the first few months to ≲2 mm/yr after 2 years. Viscoelastic relaxation of the mantle, with a time-dependent effective viscosity prescribed by a Burgers body, provides a good explanation for the postseismic crustal deformation, capturing both the spatial and temporal pattern. In the context of the Burgers body model (which involves a transient viscosity and steady state viscosity), a resolution analysis based on the singular value decomposition reveals that at most, two constraints on depth-dependent steady state mantle viscosity are provided by the present data set. Uppermost mantle viscosity (depth ≲ 60 km) is moderately resolved, but deeper viscosity structure is poorly resolved. The simplest model that explains the data better than that of uniform steady state mantle viscosity involves a linear gradient in logarithmic viscosity with depth, with a small increase from the Moho to 220 km depth. However, the viscosity increase is not statistically significant. This suggests that the depth-dependent steady state viscosity is not resolvably different from uniformity in the uppermost mantle.

Atlas of the underworld: Slab remnants in the mantle, their sinking history, and a new outlook on lower mantle viscosity

NASA Astrophysics Data System (ADS)

van der Meer, Douwe G.; van Hinsbergen, Douwe J. J.; Spakman, Wim

2018-01-01

Across the entire mantle we interpret 94 positive seismic wave-speed anomalies as subducted lithosphere and associate these slabs with their geological record. We document this as the Atlas of the Underworld, also accessible online at www.atlas-of-the-underworld.org, a compilation comprising subduction systems active in the past 300 Myr. Deeper slabs are correlated to older geological records, assuming no relative horizontal motions between adjacent slabs following break-off, using knowledge of global plate circuits, but without assuming a mantle reference frame. The longest actively subducting slabs identified reach the depth of 2500 km and some slabs have impinged on Large Low Shear Velocity Provinces in the deepest mantle. Anomously fast sinking of some slabs occurs in regions affected by long-term plume rising. We conclude that slab remnants eventually sink from the upper mantle to the core-mantle boundary. The range in subduction-age versus - depth in the lower mantle is largely inherited from the upper mantle history of subduction. We find a significant depth variation in average sinking speed of slabs. At the top of the lower mantle average slab sinking speeds are between 10 and 40 mm/yr, followed by a deceleration to 10-15 mm/yr down to depths around 1600-1700 km. In this interval, in situ time-stationary sinking rates suggest deceleration from 20 to 30 mm/yr to 4-8 mm/yr, increasing to 12-15 mm/yr below 2000 km. This corroborates the existence of a slab deceleration zone but we do not observe long-term (> 60 My) slab stagnation, excluding long-term stagnation due to compositional effects. Conversion of slab sinking profiles to viscosity profiles shows the general trend that mantle viscosity increases in the slab deceleration zone below which viscosity slowly decreases in the deep mantle. This is at variance with most published viscosity profiles that are derived from different observations, but agrees qualitatively with recent viscosity profiles suggested from material experiments.

Continent-Wide Maps of Lg Coda Q Variation and Rayleigh-wave Attenuation Variation for Eurasia

DTIC Science & Technology

2007-01-30

lithosphere and crustal strain lead us to infer that fluids, originating by hydrothermal release from subducting lithosphere or other upper mantle heat...relatively low Qo values in the Arabian Peninsula are produced by fluids that have been released in the upper mantle by hydrothermal processes and have...Advection of plumes in mantle flow: Implications for hotspot motion, mantle viscosity and plume distribution, Geophys. J. Int., 132, 412–434. Talebian, M

A mantle plume model for the Equatorial Highlands of Venus

NASA Technical Reports Server (NTRS)

Kiefer, Walter S.; Hager, Bradford H.

1991-01-01

The possibility that the Equatorial Highlands are the surface expressions of hot upwelling mantle plumes is considered via a series of mantle plume models developed using a cylindrical axisymmetric finite element code and depth-dependent Newtonian rheology. The results are scaled by assuming whole mantle convection and that Venus and the earth have similar mantle heat flows. The best model fits are for Beta and Atla. The common feature of the allowed viscosity models is that they lack a pronounced low-viscosity zone in the upper mantle. The shape of Venus's long-wavelength admittance spectrum and the slope of its geoid spectrum are also consistent with the lack of a low-viscosity zone. It is argued that the lack of an asthenosphere on Venus is due to the mantle of Venus being drier than the earth's mantle. Mantle plumes may also have contributed to the formation of some smaller highland swells, such as the Bell and Eistla regions and the Hathor/Innini/Ushas region.

Dislocation-accommodated grain boundary sliding as the major deformation mechanism of olivine in the Earth’s upper mantle

PubMed Central

Ohuchi, Tomohiro; Kawazoe, Takaaki; Higo, Yuji; Funakoshi, Ken-ichi; Suzuki, Akio; Kikegawa, Takumi; Irifune, Tetsuo

2015-01-01

Understanding the deformation mechanisms of olivine is important for addressing the dynamic processes in Earth’s upper mantle. It has been thought that dislocation creep is the dominant mechanism because of extrapolated laboratory data on the plasticity of olivine at pressures below 0.5 GPa. However, we found that dislocation-accommodated grain boundary sliding (DisGBS), rather than dislocation creep, dominates the deformation of olivine under middle and deep upper mantle conditions. We used a deformation-DIA apparatus combined with synchrotron in situ x-ray observations to study the plasticity of olivine aggregates at pressures up to 6.7 GPa (that is, ~200-km depth) and at temperatures between 1273 and 1473 K, which is equivalent to the conditions in the middle region of the upper mantle. The creep strength of olivine deforming by DisGBS is apparently less sensitive to pressure because of the competing pressure-hardening effect of the activation volume and pressure-softening effect of water fugacity. The estimated viscosity of olivine controlled by DisGBS is independent of depth and ranges from 1019.6 to 1020.7 Pa·s throughout the asthenospheric upper mantle with a representative water content (50 to 1000 parts per million H/Si), which is consistent with geophysical viscosity profiles. Because DisGBS is a grain size–sensitive creep mechanism, the evolution of the grain size of olivine is an important process controlling the dynamics of the upper mantle. PMID:26601281

Dislocation-accommodated grain boundary sliding as the major deformation mechanism of olivine in the Earth's upper mantle.

PubMed

Ohuchi, Tomohiro; Kawazoe, Takaaki; Higo, Yuji; Funakoshi, Ken-Ichi; Suzuki, Akio; Kikegawa, Takumi; Irifune, Tetsuo

2015-10-01

Understanding the deformation mechanisms of olivine is important for addressing the dynamic processes in Earth's upper mantle. It has been thought that dislocation creep is the dominant mechanism because of extrapolated laboratory data on the plasticity of olivine at pressures below 0.5 GPa. However, we found that dislocation-accommodated grain boundary sliding (DisGBS), rather than dislocation creep, dominates the deformation of olivine under middle and deep upper mantle conditions. We used a deformation-DIA apparatus combined with synchrotron in situ x-ray observations to study the plasticity of olivine aggregates at pressures up to 6.7 GPa (that is, ~200-km depth) and at temperatures between 1273 and 1473 K, which is equivalent to the conditions in the middle region of the upper mantle. The creep strength of olivine deforming by DisGBS is apparently less sensitive to pressure because of the competing pressure-hardening effect of the activation volume and pressure-softening effect of water fugacity. The estimated viscosity of olivine controlled by DisGBS is independent of depth and ranges from 10(19.6) to 10(20.7) Pa·s throughout the asthenospheric upper mantle with a representative water content (50 to 1000 parts per million H/Si), which is consistent with geophysical viscosity profiles. Because DisGBS is a grain size-sensitive creep mechanism, the evolution of the grain size of olivine is an important process controlling the dynamics of the upper mantle.

Length-scales of Slab-induced Asthenospheric Deformation from Geodynamic Modeling, Mantle Deformation Fabric, and Synthetic Shear Wave Splitting

NASA Astrophysics Data System (ADS)

Jadamec, M. A.; MacDougall, J.; Fischer, K. M.

2017-12-01

The viscosity structure of the Earth's interior is critically important, because it places a first order constraint on plate motion and mantle flow rates. Geodynamic models using a composite viscosity based on experimentally derived flow laws for olivine aggregates show that lateral viscosity variations emerge in the upper mantle due to the subduction dynamics. However, the length-scale of this transition is still not well understood. Two-dimensional numerical models of subduction are presented that investigate the effect of initial slab dip, maximum yield stress (slab strength), and viscosity formulation (Newtonian versus composite) on the emergent lateral viscosity variations in the upper-mantle and magnitude of slab-driven mantle flow velocity. Significant viscosity reductions occur in regions of large flow velocity gradients due to the weakening effect of the dislocation creep deformation mechanism. The dynamic reductions in asthenospheric viscosity (less than 1018 Pa s) occur within approximately 500 km from driving force of the slab, with peak flow velocities occurring in models with a lower yield stress (weaker slab) and higher stress exponent. This leads to a sharper definition of the rheological base of the lithosphere and implies lateral variability in tractions along the base of the lithosphere. As the dislocation creep mechanism also leads to mantle deformation fabric, we then examine the spatial variation in the LPO development in the asthenosphere and calculate synthetic shear wave splitting. The models show that olivine LPO fabric in the asthenosphere generally increases in alignment strength with increased proximity to the slab, but can be transient and spatially variable on small length scales. The vertical flow fields surrounding the slab tip can produce shear-wave splitting variations with back-azimuth that deviate from the predictions of uniform trench-normal anisotropy, a result that bears on the interpretation of complexity in shear-wave splitting observed in real subduction zones.

Mantle dynamics in super-Earths: Post-perovskite rheology and self-regulation of viscosity

NASA Astrophysics Data System (ADS)

Tackley, P. J.; Ammann, M.; Brodholt, J. P.; Dobson, D. P.; Valencia, D.

2013-07-01

The discovery of extra-solar "super-Earth" planets with sizes up to twice that of Earth has prompted interest in their possible lithosphere and mantle dynamics and evolution. Simple scalings suggest that super-Earths are more likely than an equivalent Earth-sized planet to be undergoing plate tectonics. Generally, viscosity and thermal conductivity increase with pressure while thermal expansivity decreases, resulting in lower convective vigour in the deep mantle, which, if extralopated to the largest super-Earths might, according to conventional thinking, result in no convection in their deep mantles due to the very low effective Rayleigh number. Here we evaluate this. First, as the mantle of a super-Earth is made mostly of post-perovskite we here extend the density functional theory (DFT) calculations of post-perovskite activation enthalpy of to a pressure of 1 TPa, for both slowest diffusion (upper-bound rheology) and fastest diffusion (lower-bound rheology) directions. Along a 1600 K adiabat the upper-bound rheology would lead to a post-perovskite layer of a very high (˜1030 Pa s) but relatively uniform viscosity, whereas the lower-bound rheology leads to a post-perovskite viscosity increase of ˜7 orders of magnitude with depth; in both cases the deep mantle viscosity would be too high for convection. Second, we use these DFT-calculated values in statistically steady-state numerical simulations of mantle convection and lithosphere dynamics of planets with up to ten Earth masses. The models assume a compressible mantle including depth-dependence of material properties and plastic yielding induced plate-like lithospheric behaviour. Results confirm the likelihood of plate tectonics for planets with Earth-like surface conditions (temperature and water) and show a self-regulation of deep mantle temperature. The deep mantle is not adiabatic; instead feedback between internal heating, temperature and viscosity regulates the temperature such that the viscosity has the value needed to facilitate convective loss of the radiogenic heat, which results in a very hot perovskite layer for the upper-bound rheology, a super-adiabatic perovskite layer for the lower-bound rheology, and an azimuthally-averaged viscosity of no more than 1026 Pa s. Convection in large super-Earths is characterised by large upwellings (even with zero basal heating) and small, time-dependent downwellings, which for large super-Earths merge into broad downwellings. In the context of planetary evolution, if, as is likely, a super-Earth was extremely hot/molten after its formation, it is thus likely that even after billions of years its deep interior is still extremely hot and possibly substantially molten with a "super basal magma ocean" - a larger version of the proposal of Labrosse et al. (Labrosse, S., Hernlund, J.W., Coltice, N. [2007]. Nature 450, 866-869), although this depends on presently unknown melt-solid density contrast and solidus.

Constraints on the rheology of the lower crust in a strike-slip plate boundary: evidence from the San Quintín xenoliths, Baja California, Mexico

NASA Astrophysics Data System (ADS)

van der Werf, Thomas; Chatzaras, Vasileios; Marcel Kriegsman, Leo; Kronenberg, Andreas; Tikoff, Basil; Drury, Martyn R.

2017-12-01

The rheology of lower crust and its transient behavior in active strike-slip plate boundaries remain poorly understood. To address this issue, we analyzed a suite of granulite and lherzolite xenoliths from the upper Pleistocene-Holocene San Quintín volcanic field of northern Baja California, Mexico. The San Quintín volcanic field is located 20 km east of the Baja California shear zone, which accommodates the relative movement between the Pacific plate and Baja California microplate. The development of a strong foliation in both the mafic granulites and lherzolites, suggests that a lithospheric-scale shear zone exists beneath the San Quintín volcanic field. Combining microstructural observations, geothermometry, and phase equilibria modeling, we estimated that crystal-plastic deformation took place at temperatures of 750-890 °C and pressures of 400-560 MPa, corresponding to 15-22 km depth. A hot crustal geotherm of 40 ° C km-1 is required to explain the estimated deformation conditions. Infrared spectroscopy shows that plagioclase in the mafic granulites is relatively dry. Microstructures are interpreted to show that deformation in both the uppermost lower crust and upper mantle was accommodated by a combination of dislocation creep and grain-size-sensitive creep. Recrystallized grain size paleopiezometry yields low differential stresses of 12-33 and 17 MPa for plagioclase and olivine, respectively. The lower range of stresses (12-17 MPa) in the mafic granulite and lherzolite xenoliths is interpreted to be associated with transient deformation under decreasing stress conditions, following an event of stress increase. Using flow laws for dry plagioclase, we estimated a low viscosity of 1.1-1.3×1020 Pa ṡ s for the high temperature conditions (890 °C) in the lower crust. Significantly lower viscosities in the range of 1016-1019 Pa ṡ s, were estimated using flow laws for wet plagioclase. The shallow upper mantle has a low viscosity of 5.7×1019 Pa ṡ s, which indicates the lack of an upper-mantle lid beneath northern Baja California. Our data show that during post-seismic transients, the upper mantle and the lower crust in the Pacific-Baja California plate boundary are characterized by similar and low differential stress. Transient viscosity of the lower crust is similar to the viscosity of the upper mantle.

Postglacial rebound with a non-Newtonian upper mantle and a Newtonian lower mantle rheology

NASA Technical Reports Server (NTRS)

Gasperini, Paolo; Yuen, David A.; Sabadini, Roberto

1992-01-01

A composite rheology is employed consisting of both linear and nonlinear creep mechanisms which are connected by a 'transition' stress. Background stress due to geodynamical processes is included. For models with a non-Newtonian upper-mantle overlying a Newtonian lower-mantle, the temporal responses of the displacements can reproduce those of Newtonian models. The average effective viscosity profile under the ice-load at the end of deglaciation turns out to be the crucial factor governing mantle relaxation. This can explain why simple Newtonian rheology has been successful in fitting the uplift data over formerly glaciated regions.

Impact of lithospheric rheology on surface topography

NASA Astrophysics Data System (ADS)

Liao, K.; Becker, T. W.

2017-12-01

The expression of mantle flow such as due to a buoyant plume as surface topography is a classical problem, yet the role of rheological complexities could benefit from further exploration. Here, we investigate the topographic expressions of mantle flow by means of numerical and analytical approaches. In numerical modeling, both conventional, free-slip and more realistic, stress-free boundary conditions are applied. For purely viscous rheology, a high viscosity lithosphere will lead to slight overestimates of topography for certain settings, which can be understood by effectively modified boundary conditions. Under stress-free conditions, numerical and analytical results show that the magnitude of dynamic topography decreases with increasing lithosphere thickness (L) and viscosity (ηL), as L-1 and ηL-3. The wavelength of dynamic topography increases linearly with L and (ηL/ ηM) 1/3. We also explore the time-dependent interactions of a rising plume with the lithosphere. For a layered lithosphere with a decoupling weak lower crust embedded between stronger upper crust and lithospheric mantle, dynamic topography increases with a thinner and weaker lower crust. The dynamic topography saturates when the decoupling viscosity is 3-4 orders lower than the viscosity of upper crust and lithospheric mantle. We further explore the role of visco-elastic and visco-elasto-plastic rheologies.

Rheological structure of the lithosphere in plate boundary strike-slip fault zones

NASA Astrophysics Data System (ADS)

Chatzaras, Vasileios; Tikoff, Basil; Kruckenberg, Seth C.; Newman, Julie; Titus, Sarah J.; Withers, Anthony C.; Drury, Martyn R.

2016-04-01

How well constrained is the rheological structure of the lithosphere in plate boundary strike-slip fault systems? Further, how do lithospheric layers, with rheologically distinct behaviors, interact within the strike-slip fault zones? To address these questions, we present rheological observations from the mantle sections of two lithospheric-scale, strike-slip fault zones. Xenoliths from ˜40 km depth (970-1100 ° C) beneath the San Andreas fault system (SAF) provide critical constraints on the mechanical stratification of the lithosphere in this continental transform fault. Samples from the Bogota Peninsula shear zone (BPSZ, New Caledonia), which is an exhumed oceanic transform fault, provide insights on lateral variations in mantle strength and viscosity across the fault zone at a depth corresponding to deformation temperatures of ˜900 ° C. Olivine recrystallized grain size piezometry suggests that the shear stress in the SAF upper mantle is 5-9 MPa and in the BPSZ is 4-10 MPa. Thus, the mantle strength in both fault zones is comparable to the crustal strength (˜10 MPa) of seismogenic strike-slip faults in the SAF system. Across the BPSZ, shear stress increases from 4 MPa in the surrounding rocks to 10 MPa in the mylonites, which comprise the core of the shear zone. Further, the BPSZ is characterized by at least one order of magnitude difference in the viscosity between the mylonites (1018 Paṡs) and the surrounding rocks (1019 Paṡs). Mantle viscosity in both the BPSZ mylonites and the SAF (7.0ṡ1018-3.1ṡ1020 Paṡs) is relatively low. To explain our observations from these two strike-slip fault zones, we propose the "lithospheric feedback" model in which the upper crust and lithospheric mantle act together as an integrated system. Mantle flow controls displacement and the upper crust controls the stress magnitude in the system. Our stress data combined with data that are now available for the middle and lower crustal sections of other transcurrent fault systems support the prediction for constant shear strength (˜10 MPa) throughout the lithosphere; the stress magnitude is controlled by the shear strength of the upper crustal faults. Fault rupture in the upper crust induces displacement rate loading of the upper mantle, which in turn, causes strain localization in the mantle shear zone beneath the strike-slip fault. Such forced localization leads to higher stresses and strain rates in the shear zone compared to the surrounding rocks. Low mantle viscosity within the shear zone is critical for facilitating mantle flow, which induces widespread crustal deformation and displacement loading. The lithospheric feedback model suggests that strike-slip fault zones are not mechanically stratified in terms of shear stress, and that it is the time-dependent interaction of the different lithospheric layers - rather than their relative strengths - that governs the rheological behavior of the plate boundary, strike-slip fault zones.

Abrupt Upper-Plate Tilting Upon Slab-Transition-Zone Collision

NASA Astrophysics Data System (ADS)

Crameri, F.; Lithgow-Bertelloni, C. R.

2017-12-01

During its sinking, the remnant of a surface plate crosses and interacts with multiple boundaries in Earth's interior. The most-prominent dynamic interaction arises at the upper-mantle transition zone where the sinking plate is strongly affected by the higher-viscosity lower mantle. Within our numerical model, we unravel, for the first time, that this very collision of the sinking slab with the transition zone induces a sudden, dramatic downward tilt of the upper plate towards the subduction trench. The slab-transition zone collision sets parts of the higher-viscosity lower mantle in motion. Naturally, this then induces an overall larger return flow cell that, at its onset, tilts the upper plate abruptly by around 0.05 degrees and over around 10 Millions of years. Such a significant and abrupt variation in surface topography should be clearly visible in temporal geologic records of large-scale surface elevation and might explain continental-wide tilting as observed in Australia since the Eocene or North America during the Phanerozoic. Unravelling this crucial mantle-lithosphere interaction was possible thanks to state-of-the-art numerical modelling (powered by StagYY; Tackley 2008, PEPI) and post-processing (powered by StagLab; www.fabiocrameri.ch/software). The new model that is introduced here to study the dynamically self-consistent temporal evolution of subduction features accurate subduction-zone topography, robust single-sided plate sinking, stronger plates close to laboratory values, an upper-mantle phase transition and, crucially, simple continents at a free surface. A novel, fully-automated post-processing includes physical model diagnostics like slab geometry, mantle flow pattern, upper-plate tilt angle and trench location.

Thermal evolution of flattening slab and formation of wet plume: Insight into intraplate volcanism in East Asia

NASA Astrophysics Data System (ADS)

He, L.

2016-12-01

Geophysical observations imply the intraplate volcanism in East Asia is related to dehydration of slab stagnating in the transition zone. To better understand the dynamics of such process, a thermochemical mantle convection model is constructed to simulate numerically the thermal evolution of slab and the transportation of water in the process of subduction. Equation of water transfer is explicitly included, and water effects on density and viscosity are considered. Modeling results indicate that behavior of water transport relates closely to the transient thermal state and viscosities both of the slab and the surrounding mantle. Generally, initiation of wet plume is mainly influenced by the viscosity of the wet layer in the uppermost slab, whereas the horizontal distance of water transport and its ascending rate is affected strongly by the viscosity of the big mantle wedge. Whether water can be carried successfully by slab into the mantle transition zone and trigger wet plume at the surface of flattening slab depends on the viscosity contrast between wet layer and surrounding mantle. The complex fluid flow superposed by corner flow and free thermal convection controls the water transport pattern in the upper mantle. Modeling results together with previous modeling infer three stages of water circulation in the big mantle wedge: 1) water is brought into the mantle transition zone by downward subducting slab when water layer viscosity is much higher than the wedge viscosity, otherwise water is released at shallow depth near wedge tip; 2) wet plume generates from surface of warm flattening slab if containing water, which arrives at the lithospheric base and induces melting; and 3) water spreads all over the big mantle wedge, mantle convection within the big mantle wedge becomes more active, leading to upwelling of asthenosphere and erosion of the overriding continental lithosphere. Wet plume from the flattening Pacific Plate can explain the intraplate Cenozoic volcanoes in East Asia.

Transient rheology of the upper mantle beneath central Alaska inferred from the crustal velocity field following the 2002 Denali earthquake

USGS Publications Warehouse

Pollitz, F.F.

2005-01-01

The M7.9 2002 Denali earthquake, Alaska, is one of the largest strike-slip earthquakes ever recorded. The postseismic GPS velocity field around the 300-km-long rupture is characterized by very rapid horizontal velocity up to ???300 mm/yr for the first 0.1 years and slower but still elevated horizontal velocity up to ???100 mm/yr for the succeeding 1.5 years. I find that the spatial and temporal pattern of the displacement field may be explained by a transient mantle rheology. Representing the regional upper mantle as a Burghers body, I infer steady state and transient viscosities of ??1 = 2.8 ?? 1018 Pa s and ??2 = 1.0 ?? 1017 Pa s, respectively, corresponding to material relaxation times of 1.3 and 0.05 years. The lower crustal viscosity is poorly constrained by the considered horizontal velocity field, and the quoted mantle viscosities assume a steady state lower crust viscosity that is 7??1. Systematic bias in predicted versus observed velocity vectors with respect to a fixed North America during the first 3-6 months following the earthquake is reduced when all velocity vectors are referred to a fixed site. This suggests that the post-Denali GPS time series for the first 1.63 years are shaped by a combination of a common mode noise source during the first 3-6 months plus viscoelastic relaxation controlled by a transient mantle rheology.

Reconciling Long-Wavelength Dynamic Topography, Geoid Anomalies and Mass Distribution on Earth

NASA Astrophysics Data System (ADS)

Hoggard, M.; Richards, F. D.; Ghelichkhan, S.; Austermann, J.; White, N.

2017-12-01

Since the first satellite observations in the late 1950s, we have known that that the Earth's non-hydrostatic geoid is dominated by spherical harmonic degree 2 (wavelengths of 16,000 km). Peak amplitudes are approximately ± 100 m, with highs centred on the Pacific Ocean and Africa, encircled by lows in the vicinity of the Pacific Ring of Fire and at the poles. Initial seismic tomography models revealed that the shear-wave velocity, and therefore presumably the density structure, of the lower mantle is also dominated by degree 2. Anti-correlation of slow, probably low density regions beneath geoid highs indicates that the mantle is affected by large-scale flow. Thus, buoyant features are rising and exert viscous normal stresses that act to deflect the surface and core-mantle boundary (CMB). Pioneering studies in the 1980s showed that a viscosity jump between the upper and lower mantle is required to reconcile these geoid and tomographically inferred density anomalies. These studies also predict 1-2 km of dynamic topography at the surface, dominated by degree 2. In contrast to this prediction, a global observational database of oceanic residual depth measurements indicates that degree 2 dynamic topography has peak amplitudes of only 500 m. Here, we attempt to reconcile observations of dynamic topography, geoid, gravity anomalies and CMB topography using instantaneous flow kernels. We exploit a density structure constructed from blended seismic tomography models, combining deep mantle imaging with higher resolution upper mantle features. Radial viscosity structure is discretised, and we invert for the best-fitting viscosity profile using a conjugate gradient search algorithm, subject to damping. Our results suggest that, due to strong sensitivity to radial viscosity structure, the Earth's geoid seems to be compatible with only ± 500 m of degree 2 dynamic topography.

Thermal and mechanical structure of the upper mantle: A comparison between continental and oceanic models

NASA Technical Reports Server (NTRS)

Froidevaux, C.; Schubert, G.; Yuen, D. A.

1976-01-01

Temperature, velocity, and viscosity profiles for coupled thermal and mechanical models of the upper mantle beneath continental shields and old ocean basins show that under the continents, both tectonic plates and the asthenosphere, are thicker than they are beneath the oceans. The minimum value of viscosity in the continental asthenosphere is about an order of magnitude larger than in the shear zone beneath oceans. The shear stress or drag underneath continental plates is also approximately an order of magnitude larger than the drag on oceanic plates. Effects of shear heating may account for flattening of ocean floor topography and heat flux in old ocean basins.

Viscoelastic responses of a hard transition zone - Effects on postglacial uplifts and rotational signatures

NASA Technical Reports Server (NTRS)

Spada, Giorgio; Sabadini, Roberto; Yuen, David A.

1991-01-01

A five-layer viscoelastic spherical model is used to calculate the transient displacements of postglacial rebound, the induced polar motions, and the temporal variations of the geopotential up to degree 8 of the zonal coefficients. Two models - one with two viscoelastic layers separated at 670 km, and the other with three layers in which a hard garnet layer lies between the upper and lower mantle - are compared. Forward modeling shows that it may be possible to discern the presence of a hard garnet layer with a viscosity of at least ten times greater than the upper mantle, on the basis of uplift data near the center of the former Laurentide ice-sheet and from polar wander and j2 data. Temporal variations of higher gravity harmonics, such as j6 and j8, can potentially place even tighter constraints on the rheological properties of the hard transition zone. A lower mantle viscosity between 2 and 4 x 10 to the 22nd Pa is generally preferred in models with a garnet layer which may be as large as 50 times more viscous than the upper mantle.

Episodic large-scale overturn of two-layer mantles in terrestrial planets

NASA Astrophysics Data System (ADS)

Herrick, D. L.; Parmentier, E. M.

1994-01-01

It is usually assumed that the upper and lower mantles of a chemically stratified planet are arranged so that the upper mantle is chemically less dense and that these layers convect separately. Possible buoyant overturn of the two mantle layers has not previously been considered. Such overturn would initially occur when thermal expansion of a chemically denser lower mantle more than offsets the compositional density difference between the layers, reversing the relative sense of buoyancy. Once overturn has occurred, the chemically denser, but thermally less dense upper mantle cools more efficiently than the lower mantle and loses its relative thermal buoyancy. If mixing is slow, this leads to repeated overturns that result in thermal histories that differ radically from those obtained without this large-scale overturning. Thermal evolution calculations, for a two-layer mantle over a wide range of parameter space, show that large-scale overturn occurs cyclically with a well-defined period. This period depends most strongly on the viscosity of the lower mantle, to which it is approximately proportional. Geologically interesting overturn periods on the order of 107 to 109 years result for lower mantle viscosities of 1022 to 1024 Pa s for the Earth and Venus, and 1021 to 1023 Pa s for Mars. The mantles of Mercury and the Moon are too thin to permit two-layer convection, and therefore the model is not appropriate for them. Overturn cannot occur on Earth or Venus if the compositional density difference between the layers exceeds about 4%, or on Mars if it exceeds about 2%. Large-scale mantle overturn could have significant tectonic consequences such as the initiation of a new plate tectonic cycle on the Earth or a major resurfacing event on Mars or Venus. Such episodic events in the evolution of a planet are not easily explained by whole mantle thermal convection.

Episodic large-scale overturn of two-layer mantles in terrestrial planets

NASA Technical Reports Server (NTRS)

Herrick, David L.; Parmentier, E. M.

1994-01-01

It is usually assumed that the upper and lower mantles of a chemically stratified planet are arranged so that the upper mantle is chemically less dense and that these layers convect separately. Possible buoyant overturn of the two mantle layers has not previously been considered. Such overturn would initially occur when thermal expansion of a chemically denser lower mantle more than offsets the compositional density difference between the layers, reversing the relative sense of buoyancy. Once overturn has occurred, the chemically denser, but thermally less dense upper mantle cools more efficiently than the lower mantle and loses its relative thermal buoyancy. If mixing is slow, this leads to repeated overturns that result in thermal histories that differ radically from those obtained without this large-scale overturning. Thermal evolution calculations, for a two-layer mantle over a wide range of parameter space, show that large-scale overturn occurs cyclically with a well-defined period. This period depends most strongly on the viscosity of the lower mantle, to which it is approximately proportional. Geologically interesting overturn periods on the order of 10(exp 7) to 10(exp 9) years result for lower mantle viscosities of 10(exp 22) to 10(exp 24) Pa s for the Earth and Venus, and 10(exp 21) to 10(exp 23) Pa s for Mars. The mantles of Mercury and the Moon are too thin to permit two-layer convection, and therefore the model is not appropriate for them. Overturn cannot occur on Earth or Venus if the compositional density difference between the layers exceeds about 4%, or on Mars if it exceeds about 2%. Large-scale mantle overturn could have significant tectonic consequences such as the initiation of a new plate tectonic cycle on the Earth or a major resurfacing event on Mars or Venus. Such episodic events in the evolution of a planet are not easily explained by whole mantle thermal convection.

Experimental investigation of flow-induced fabrics in rocks at upper-mantle pressures: Application to understanding mantle dynamics and seismic anisotropy

DOE Office of Scientific and Technical Information (OSTI.GOV)

Kohlstedt, David L.

2016-04-26

The goal of this collaborative research effort between W.B. Durham at the Massachusetts Institute of Technology (MIT) and D.L. Kohlstedt and S. Mei at the University of Minnesota (UMN) was to exploit a newly developed technology for high-pressure, high-temperature deformation experimentation, namely, the deformation DIA (D-DIA) to determine the deformation behavior of a number of important upper mantle rock types including olivine, garnet, enstatite, and periclase. Experiments were carried out under both hydrous and anhydrous conditions and at both lithospheric and asthenospheric stress and temperature conditions. The result was a group of flow laws for Earth’s upper mantle that quantitativelymore » describe the viscosity of mantle rocks from shallow depths (the lithosphere) to great depths (the asthenosphere). These flow laws are fundamental for modeling the geodynamic behavior and heat transport from depth to Earth’s surface.« less

Experimental investigation of flow-induced fabrics in rocks at upper-mantle pressures. Application to understanding mantle dynamics and seismic anisotropy

DOE Office of Scientific and Technical Information (OSTI.GOV)

Durham, William B.

2016-05-02

The goal of this collaborative research effort between W.B. Durham at the Massachusetts Institute of Technology (MIT) and D.L. Kohlstedt and S. Mei at the University of Minnesota (UMN) was to exploit a newly developed technology for high-pressure, high-temperature deformation experimentation, namely, the deformation DIA (D-DIA), to determine the deformation behavior of a number of important upper mantle rock types including olivine, garnet, enstatite, and periclase. Experiments were carried out under both hydrous and anhydrous conditions and at both lithospheric and asthenospheric stress and temperature conditions. The result was a group of flow laws for Earth’s upper mantle that quantitativelymore » describe the viscosity of mantle rocks from shallow depths (the lithosphere) to great depths (the asthenosphere). These flow laws are fundamental for modeling the geodynamic behavior and heat transport from depth to Earth’s surface.-« less

Resolving Discrepancies Between Observed and Predicted Dynamic Topography on Earth

NASA Astrophysics Data System (ADS)

Richards, F. D.; Hoggard, M.; White, N. J.

2017-12-01

Compilations of well-resolved oceanic residual depth measurements suggest that present-day dynamic topography differs from that predicted by geodynamic simulations in two significant respects. At short wavelengths (λ ≤ 5,000 km), much larger amplitude variations are observed, whereas at long wavelengths (λ > 5,000 km), observed dynamic topography is substantially smaller. Explaining the cause of this discrepancy with a view to reconciling these different approaches is central to constraining the structure and dynamics of the deep Earth. Here, we first convert shear wave velocity to temperature using an experimentally-derived anelasticity model. This relationship is calibrated using a pressure and temperature-dependent plate model that satisfies age-depth subsidence, heat flow measurements, and seismological constraints on the depth to the lithosphere-asthenosphere boundary. In this way, we show that, at short-wavelengths, observed dynamic topography is consistent with ±150 ºC asthenospheric temperature anomalies. These inferred thermal buoyancy variations are independently verified by temperature measurements derived from geochemical analyses of mid-ocean ridge basalts. Viscosity profiles derived from the anelasticity model suggest that the asthenosphere has an average viscosity that is two orders of magnitude lower than that of the underlying upper mantle. The base of this low-viscosity layer coincides with a peak in azimuthal anisotropy observed in recent seismic experiments. This agreement implies that lateral asthenospheric flow is rapid with respect to the underlying upper mantle. We conclude that improved density and viscosity models of the uppermost mantle, which combine a more comprehensive physical description of the lithosphere-asthenosphere system with recent seismic tomographic models, can help to resolve spectral discrepancies between observed and predicted dynamic topography. Finally, we explore possible solutions to the long-wavelength discrepancy that exploit the velocity to density conversion described above combined with radial variation of mantle viscosity.

Modelling the possible interaction between edge-driven convection and the Canary Islands mantle plume

NASA Astrophysics Data System (ADS)

Negredo, A. M.; Rodríguez-González, J.; Fullea, J.; Van Hunen, J.

2017-12-01

The close location between many hotspots and the edges of cratonic lithosphere has led to the hypothesis that these hotspots could be explained by small-scale mantle convection at the edge of cratons (Edge Driven Convection, EDC). The Canary Volcanic Province hotspot represents a paradigmatic example of this situation due to its close location to the NW edge of the African Craton. Geochemical evidence, prominent low seismic velocity anomalies in the upper and lower mantle, and the rough NE-SW age-progression of volcanic centers consistently point out to a deep-seated mantle plume as the origin of the Canary Volcanic Province. It has been hypothesized that the plume material could be affected by upper mantle convection caused by the thermal contrast between thin oceanic lithosphere and thick (cold) African craton. Deflection of upwelling blobs due to convection currents would be responsible for the broader and more irregular pattern of volcanism in the Canary Province compared to the Madeira Province. In this study we design a model setup inspired on this scenario to investigate the consequences of possible interaction between ascending mantle plumes and EDC. The Finite Element code ASPECT is used to solve convection in a 2D box. The compositional field and melt fraction distribution are also computed. Free slip along all boundaries and constant temperature at top and bottom boundaries are assumed. The initial temperature distribution assumes a small long-wavelength perturbation. The viscosity structure is based on a thick cratonic lithosphere progressively varying to a thin, or initially inexistent, oceanic lithosphere. The effects of assuming different rheologies, as well as steep or gradual changes in lithospheric thickness are tested. Modelling results show that a very thin oceanic lithosphere (< 30 km) is needed to generate partial melting by EDC. In this case partial melting can occur as far as 700 km away from the edge of the craton. The size of EDC cells is relatively small (diameter about 300 km) for lithosphere/asthenosphere viscosity contrasts of 1000. In contrast, models assuming temperature-dependent viscosity and large viscosity variations evolve to large-scale (upper mantle) convection cells, with upwelling of hot material being enhanced by cold downwellings at the edge of cratonic lithosphere.

Effects of basin-forming impacts on the thermal evolution and magnetic field of Mars

NASA Astrophysics Data System (ADS)

Roberts, J. H.; Arkani-Hamed, J.

2017-11-01

The youngest of the giant impact basins on Mars are either weakly magnetized or completely demagnetized, indicating that a global magnetic field was not present at the time those basins formed. Eight basins are sufficiently large that the impact heating associated with their formation could have penetrated below the core-mantle boundary (CMB). Here we investigate the thermal evolution of the martian interior and the fate of the global magnetic field using 3D mantle convection models coupled to a parameterized 1D core thermal evolution model. We find that the survival of the impact-induced temperature anomalies in the upper mantle is strongly controlled by the mantle viscosity. Impact heating from subsequent impacts can accumulate in stiffer mantles faster than it can be advected away, resulting in a thermal blanket that insulates an entire hemisphere. The impact heating in the core will halt dynamo activity, at least temporarily. If the mantle is initially cold, and the core initially superheated, dynamo activity may resume as quickly as a few Myr after each impact. However unless the lower mantle has either a low viscosity or a high thermal conductivity, this restored dynamo will last for only a few hundred Myr after the end of the sequence of impacts. Thus, we find that the longevity of the magnetic field is more strongly controlled by the lower mantle properties and relatively insensitive to the impact-induced temperature anomalies in the upper mantle.

Impact of the rheological layering of the lithosphere on the topography generated by sublithospheric density anomalies: Insights from analog modeling

NASA Astrophysics Data System (ADS)

Sembroni, A.; Globig, J.; Rozel, A.; Faccenna, C.; Funiciello, F.; Fernandez, M.

2013-12-01

Density anomalies located beneath the lithosphere are thought to generate dynamic topography at the surface of the Earth. Tomographic models are often used to infer the later variations of the density field in the mantle. Surface topography can then be computed using analytical solutions or numerical simulations of mantle convection. It has been shown that the viscosity profile of the upper mantle has a strong influence on the magnitude and spectral signature of surface topography and uplift rate. Here we present results from analogue modeling of the interaction between a rising ball-shaped density anomaly and the lithosphere in an isoviscous, isothermal Newtonian mantle system. Preliminary data show that surface topography is strongly influenced not only by mantle viscosity but also by density and viscosity profiles of the lithosphere. Our apparatus consists of a plexiglass square box (40x40x50 cm3) filled with glucose syrup. From the bottom a silicon ball was free to rise up until impinging a silicon plate floating on top of the syrup, mimicking the lithosphere. In order to investigate the role of lithospheric thickness and layered continental crust on stress partitioning, maximum dynamic topography, uplift rate and signal wavelength, two different configurations were tested: homogeneous lithosphere and stratified lithosphere including a low-viscosity lower crust. The topographic evolution of the surface was tracked using a laser scanning the top of the apparatus. The rise of the density anomaly was recorded by a side camera. We observe that a thick and then more resistant lithosphere makes up to 2 times lower and laterally wider topographic signatures. Layered lithospheres including a decoupling lower crust decrease the equilibrium topography and its lateral extend by ~30% to 40%. Most importantly, the uplift rate is strongly affected by the choice of lithosphere model. Both lithosphere width and the presence of a decoupling lower crust may modify the uplift rate by a factor 3. Thus, depending on the lithosphere rheology, we show that uplift rate may vary by one order of magnitude, for the same density anomaly and mantle viscosity. This result shows that surface uplift rate can be used to infer the viscosity of the upper mantle in specific Earth regions only if the rheology of the lithosphere is well constrained. With respect to previous approaches, whether numerical or analog modeling of dynamic topography, our experiments represent a new attempt to investigate the propagation of normal stresses generated by mantle flow through a rheologically stratified lithosphere and its resulting topographic signal.

Crust and Mantle Deformation Revealed from High-Resolution Radially Anisotropic Velocity Models

NASA Astrophysics Data System (ADS)

Li, A.; Dave, R.; Yao, Y.

2017-12-01

Love wave tomography, which can achieve a similar model resolution as Rayleigh wave, so far has limited applications to the USArray data. Recently, we have developed high-resolution Love wave phase velocity maps in the Wyoming craton and Texas using data at the Transportable Array stations. 3-D, radially anisotropic velocity models are obtained by jointly inverting Love and Rayleigh wave phase velocities. A high-velocity anomaly extending to about 200 km depth beneath central Wyoming correlates with negative radial anisotropy (Vsv>Vsh), suggesting that mantle downwelling develops under the cratonic lithosphere. Surprisingly, the significantly low velocity beneath the Yellowstone hotspot, which has been interpreted as partial melting and asthenospheric upwelling, is associated with the largest radial anisotropy (Vsh>Vsv) in the area. This observation does not support mantle upwelling. Instead, it indicates that the upper mantle beneath the hotspot has experienced strong shear deformation probably by the plate motion and large-scale mantle flow. In Texas, positive radial anisotropy in the lower crust extends from the coast to the Ouachita belt, which is characterized by high velocity and negative radial anisotropy. In the upper mantle, large variations of velocity and anisotropy exit under the coastal plain. A common feature in these anisotropic models is that high-velocity anomalies in the upper mantle often correlate with negative anisotropy (Vsv>Vsh) while low-velocity anomalies are associated with positive anisotropy (Vsh>Vsv). The manifestation of mantle downweling as negative radial anisotropy is largely due to the relatively high viscosity of the high-velocity mantle block, which is less affected by the surrounding large-scale horizontal flow. However, mantle upwelling, which is often associated with low-velocity anomalies, presumably low-viscosity mantle blocks, is invisible in radial anisotropy models. Such upwelling may happen too quickly to make last effects or too slow to alter the dominant shear deformation in the asthenosphere.

Mantle flow tectonics - The influence of a ductile lower crust and implications for the formation of topographic uplands on Venus

NASA Technical Reports Server (NTRS)

Bindschadler, Duane L.; Parmentier, E. Marc

1990-01-01

The crust and mantle of Venus can be represented by a model of a layered structure stratified in both density and viscosity. This structure consists of a brittle-elastic upper crustal layer; a ductile weaker crustal layer; a strong upper mantle layer, about 10 percent denser than the crust; and a weaker substrate, representing the portion of the mantle in which convective flow occurs which is a primary source of large-scale topographic and tectonic features. This paper examines the interactions between these four layers and the mantle flow driven by thermal or compositional variations. Solutions are found for a flow driven by a buoyancy-force distribution within the mantle and by relief at the surface and crust-mantle boundary. It is shown that changes in crustal thickness are driven by vertical normal stresses due to mantle flow and by shear coupling of horizontal mantle flow into the crust.

Wet plume atop of the flattening slab: Insight into intraplate volcanism in East Asia

NASA Astrophysics Data System (ADS)

He, Lijuan

2017-08-01

Geophysical observations imply the intraplate volcanism in East Asia is related to dehydration of slab stagnating in the transition zone. To better understand the dynamics of such process, a thermochemical mantle convection model is constructed to simulate numerically the thermal evolution of slab and the transportation of water in the process of slab downgoing, flattening and stagnation. Equation of water transfer is included, and water effects on density and viscosity are considered. Model results indicate the warming of slab by surrounding mantle is rather slow. Water could be successfully dragged into the transition zone if the reference viscosity of the hydrous layer (with initial water of 2 wt%) is higher than 1017 Pa s and that of mantle is 1021 Pa s. Wet plumes could then originate in the flat-lying part of the slab, relatively far from the trench. Generally, the viscosity of the hydrous layer governs the initiation of wet plume, whereas the viscosity of the overlying mantle wedge controls the activity of the ascending wet plumes - they are more active in the weaker wedge. The complex fluid flow superposed by corner flow and free thermal convection influences greatly the water transport pattern in the upper mantle. Modeling results together with previous modeling infer three stages of water circulation in the big mantle wedge: 1) water is brought into the mantle transition zone by downward subducting slab under some specific thermo-rheological conditions, otherwise water is released at shallow depth near wedge tip; 2) wet plume generates from surface of the flattening slab warmed by surrounding mantle, and 3) water spreads over the big mantle wedge. Wet plume from the flattening Pacific Plate arrives at the lithospheric base and induces melting, which can explain the intraplate Cenozoic volcanoes in East Asia.

The role of upper mantle mineral phase transitions on the current structure of large-scale Earth's mantle convection.

NASA Astrophysics Data System (ADS)

Thoraval, C.

2017-12-01

Describing the large-scale structures of mantle convection and quantifying the mass transfer between upper and lower mantle request to account for the role played by mineral phase transitions in the transition zone. We build a density distribution within the Earth mantle from velocity anomalies described by global seismic tomographic models. The density distribution includes thermal anomalies and topographies of the phase transitions at depths of 410 and 660 km. We compute the flow driven by this density distribution using a 3D spherical circulation model, which account for depth-dependent viscosity. The dynamic topographies at the surface and at the CMB and the geoid are calculated as well. Within the range of viscosity profiles allowing for a satisfying restitution of the long wavelength geoid, we perform a parametric study to decipher the role of the characteristics of phase diagrams - mainly the Clapeyron's slopes - and of the kinetics of phase transitions, which may modify phase transition topographies. Indeed, when a phase transition is delayed, the boundary between two mineral phases is both dragged by the flow and interfere with it. The results are compared to recent estimations of surface dynamic topography and to the phase transition topographies as revealed by seismic studies. The consequences are then discussed in terms of structure of mantle flow. Comparisons between various tomographic models allow us to enlighten the most robust features. At last, the role played by the phase transitions on the lateral variations of mass transfer between upper and lower mantle are quantified by comparison to cases with no phase transitions and confronted to regional tomographic models, which reflect the variability of the behaviors of the descending slabs in the transition zone.

The Importance of Lower Mantle Structure to Plate Stresses and Plate Motions

NASA Astrophysics Data System (ADS)

Holt, W. E.; Wang, X.; Ghosh, A.

2016-12-01

Plate motions and plate stresses are widely assumed as the surface expression of mantle convection. The generation of plate tectonics from mantle convection has been studied for many years. Lithospheric thickening (or ridge push) and slab pull forces are commonly accepted as the major driving forces for the plate motions. However, the importance of the lower mantle to plate stresses and plate motions remains less clear. Here, we use the joint modeling of lithosphere and mantle dynamics approach of Wang et al. (2015) to compute the tractions originating from deeper mantle convection and follow the method of Ghosh et al. (2013) to calculate gravitational potential energy per unit area (GPE) based on Crust 1.0 (Laske et al., 2013). Absolute values of deviatoric stresses are determined by the body force distributions (GPE gradients and traction magnitudes applied at the base of the lithosphere). We use the same relative viscosity model that Ghosh et al. (2013) used, and we solve for one single adjustable scaling factor that multiplies the entire relative viscosity field to provide absolute values of viscosity throughout the lithosphere. This distribution of absolute values of lithosphere viscosities defines the magnitudes of surface motions. In this procedure, the dynamic model first satisfies the internal constraint of no-net-rotation of motions. The model viscosity field is then scaled by the single factor until we achieve a root mean square (RMS) minimum between computed surface motions and the kinematic no-net-rotation (NNR) model of Kreemer et al. (2006). We compute plate stresses and plate motions from recently published global tomography models (over 70 based on Wang et al., 2015). We find that RMS misfits are significantly reduced when details of lower mantle structure from the latest tomography models are added to models that contain only upper and mid-mantle density distributions. One of the key reasons is that active upwelling from the Large Low Shear Velocity Provinces (LLSVPs) in the lower mantle in Pacific (Frost and Rost, 2014) provides important components of mantle flow affecting plate stresses and motions. We demonstrate in this paper how lower mantle density heterogeneity has a marked influence on plate stresses and plate motions.

Glacial isostatic adjustment model with composite 3-D Earth rheology for Fennoscandia

NASA Astrophysics Data System (ADS)

van der Wal, Wouter; Barnhoorn, Auke; Stocchi, Paolo; Gradmann, Sofie; Wu, Patrick; Drury, Martyn; Vermeersen, Bert

2013-07-01

Models for glacial isostatic adjustment (GIA) can provide constraints on rheology of the mantle if past ice thickness variations are assumed to be known. The Pleistocene ice loading histories that are used to obtain such constraints are based on an a priori 1-D mantle viscosity profile that assumes a single deformation mechanism for mantle rocks. Such a simplified viscosity profile makes it hard to compare the inferred mantle rheology to inferences from seismology and laboratory experiments. It is unknown what constraints GIA observations can provide on more realistic mantle rheology with an ice history that is not based on an a priori mantle viscosity profile. This paper investigates a model for GIA with a new ice history for Fennoscandia that is constrained by palaeoclimate proxies and glacial sediments. Diffusion and dislocation creep flow law data are taken from a compilation of laboratory measurements on olivine. Upper-mantle temperature data sets down to 400 km depth are derived from surface heatflow measurements, a petrochemical model for Fennoscandia and seismic velocity anomalies. Creep parameters below 400 km are taken from an earlier study and are only varying with depth. The olivine grain size and water content (a wet state, or a dry state) are used as free parameters. The solid Earth response is computed with a global spherical 3-D finite-element model for an incompressible, self-gravitating Earth. We compare predictions to sea level data and GPS uplift rates in Fennoscandia. The objective is to see if the mantle rheology and the ice model is consistent with GIA observations. We also test if the inclusion of dislocation creep gives any improvements over predictions with diffusion creep only, and whether the laterally varying temperatures result in an improved fit compared to a widely used 1-D viscosity profile (VM2). We find that sea level data can be explained with our ice model and with information on mantle rheology from laboratory experiments, heatflow and seismology and a pure olivine rheology above 400 km. Moreover, laterally heterogeneous models provide a significantly better fit to relative sea level data than the VM2 viscosity, for our ice model as well as for the ICE-5G model that is based on the VM2 profile. The new ice model gives different constraints on mantle rheology than the ICE-5G model, indicating a possible bias towards mantle viscosity in the latter or shortcomings in our ice model. Present-day uplift rates for a dry rheology are close to GPS observed uplift rate for certain combinations of grain size and temperature fields. Sea level data show a preference for a wet olivine rheology, but in that case uplift rates are too low for all grain sizes and temperature fields. The difficulty to fit sea level data and uplift rate data simultaneously can not be resolved by varying creep parameters below 400 km. Uncertainties in the flow law and the neglect of other materials in the upper mantle, as well as the neglect of flow in the crust could affect our conclusions.

Rheological properties of bridgmanite based on deformation experiments

NASA Astrophysics Data System (ADS)

Tsujino, N.; Yamazaki, D.; Yoshino, T.; Sakurai, M.; Higo, Y.; Tange, Y.

2017-12-01

The lower mantle occupies 65% of Earth's mantle. Therefore, rheology of the Earth's lower mantle is most important to understand dynamic processes in the Earth's mantle. In Tsujino et al. (2016), we developed deformation experimental technique using D-DIA apparatus as Kawai-type (6-8 type). Crystallographic-preferred-orientation (CPO) of bridgmanite at top of the Earth's lower mantle conditions was determined by shear deformation experiments under upper most lower mantle conditions (25 GPa and 1873 K). The observed seismic shear wave anisotropies near several subducted slabs (Tonga-Kermadec, Kurile, Peru and Java) can be explained in terms of the CPO of bridgmanite as induced by mantle flow parallel to the direction of subduction. On the other hands, one dimensional viscosity models of the Earth's mantle were proposed by geophysical observations while there are large inconsistencies of viscosity (2 3 order magnitude) in the lower mantle between suggested models. It is important to determine viscosity of lower mantle minerals by high pressure experiments in order to understand mantle dynamics. In this study, we conducted in-situ stress-strain measurements of MgSiO3-bridgmanite aggregate at 1473-1673 K and 24 GPa using D-DIA type apparatus as Kawai-type at Spring-8 BL04B1. Measured uniaxial stress, strain rate of bridgmanite during deformation experiments were 0.3-1.3 GPa and 4×10-6 - 3×10-5 /s with <6% strain. Creep strength of bridgmanite at 1×10-5 /s is largest in the mantle minerals and 0.5-1 order magnitude larger than those of transition minerals when only the results using D-DIA apparatus are compared.

An Inversion of Gravity and Topography for Mantle and Crustal Structure on Mars

NASA Technical Reports Server (NTRS)

Kiefer, Walter S.; Bills, Bruce G.; Nerem, R. Steven

1996-01-01

Analysis of the gravity and topography of Mars presently provides our primary quantitative constraints on the internal structure of Mars. We present an inversion of the long-wavelength (harmonic degree less than or equal to 10) gravity and topography of Mars for lateral variations of mantle temperature and crustal thickness. Our formulation incorporates both viscous mantle flow (which most prior studies have neglected) and isostatically compensated density anomalies in the crust and lithosphere. Our nominal model has a 150-km-thick high-viscosity surface layer over an isoviscous mantle, with a core radius of 1840 km. It predicts lateral temperature variations of up to a few hundred degrees Kelvin relative to the mean mantle temperature, with high temperature under Tharsis and to a lesser extent under Elysium and cool temperatures elsewhere. Surprisingly, the model predicts crustal thinning beneath Tharsis. If correct, this implies that thinning of the crust by mantle shear stresses dominates over thickening of the crust by volcanism. The major impact basins (Hellas, Argyre, Isidis, Chryse, and Utopia) are regions of crustal thinning, as expected. Utopia is also predicted to be a region of hot mantle, which is hard to reconcile with the surface geology. An alternative model for Utopia treats it as a mascon basin. The Utopia gravity anomaly is consistent with the presence of a 1.2 to 1.6 km thick layer of uncompensated basalt, in good agreement with geologic arguments about the amount of volcanic fill in this area. The mantle thermal structure is the dominant contributor to the observed geoid in our inversion. The mantle also dominates the topography at the longest wavelengths, but shorter wavelengths (harmonic degrees greater than or equal to 4) are dominated by the crustal structure. Because of the uncertainty about the appropriate numerical values for some of the model's input parameters, we have examined the sensitivity of the model results to the planetary structural model (core radius and core and mantle densities), the mantle's viscosity stratification, and the mean crustal thickness. The model results are insensitive to the specific thickness or viscosity contrast of the high-viscosity surface layer and to the mean crustal thickness in the range 25 to 100 km. Models with a large core radius or with an upper mantle low-viscosity zone require implausibly large lateral variations in mantle temperature.

Crust and Upper Mantle Structure of Antarctica from Rayleigh Wave Tomography

NASA Astrophysics Data System (ADS)

Wiens, D. A.; Heeszel, D. S.; Sun, X.; Chaput, J. A.; Aster, R. C.; Nyblade, A.; Anandakrishnan, S.; Wilson, T. J.; Huerta, A. D.

2012-12-01

We combine data from three temporary arrays of seismometers (AGAP/GAMSEIS 2007-2010, ANET/POLENET 2007-2012, TAMSEIS 2001-2003) deployed across Antarctica, along with permanent stations in the region, to produce a large scale shear velocity model of the continent extending from the Gamburtsev Subglacial Mountains (GSM) in East Antarctica, across the Transantarctic Mountains (TAM) and West Antarctic Rift System (WARS) to Marie Byrd Land (MBL) in West Antarctica. Our combined dataset consists of Rayleigh wave phase and amplitude measurements from 112 stations across the study region. We first invert for 2-D Rayleigh wave phase velocities using the two-plane wave method. These results are then inverted for shear velocity structure using crustal thicknesses derived from ambient noise tomography and teleseismic receiver functions. We refine our shear velocity model by performing a Monte Carlo simulation that explores the tradeoff between crustal thickness and upper mantle seismic velocities. The resulting model is higher resolution than previous studies (~150 km resolution length) and highlights significant differences in crustal and uppermost mantle structure between East and West Antarctica in greater detail than previously possible. East Antarctica is underlain by thick crust (reaching ~55 km beneath the GSM) and fast, cratonic lithosphere. West Antarctica is defined by thinner crust and slow upper mantle velocities indicative of its more recent tectonic activity. The observed boundary in crustal thickness closely follows the TAM front. MBL is underlain by a thicker lithosphere than that observed beneath the WARS, but slow mantle velocities persist to depths greater than 200 km, indicating a 'deep seated' (i.e. deeper than the deepest resolvable features of our model) thermal source for volcanism in the region. The slowest seismic velocities at shallow depths are observed in the Terror Rift region of the Ross Sea along an arc following the TAM front, where the most recent extension has occurred, and in another region of active volcanism. The Ellsworth-Whitmore Mountains are underlain by relatively thick crust and an intermediate thickness lithosphere, consistent with its hypothesized origin as a remnant Precambrian crustal block. We also produce upper mantle viscosity models for the study region using a temperature-dependent rheology, assuming that mantle seismic anomalies are dominated by temperature variations. Initial results closely correlate with the velocity model, with viscosities beneath West Antarctica inferred to be orders of magnitude lower than beneath East Antarctica. These viscosity results have important implications for our understanding of glacial isostatic adjustment, which is of particular interest in producing models of past and future changes in the Antarctic Ice Sheets.

Large-scale compositional heterogeneity in the Earth's mantle

NASA Astrophysics Data System (ADS)

Ballmer, M.

2017-12-01

Seismic imaging of subducted Farallon and Tethys lithosphere in the lower mantle has been taken as evidence for whole-mantle convection, and efficient mantle mixing. However, cosmochemical constraints point to a lower-mantle composition that has a lower Mg/Si compared to upper-mantle pyrolite. Moreover, geochemical signatures of magmatic rocks indicate the long-term persistence of primordial reservoirs somewhere in the mantle. In this presentation, I establish geodynamic mechanisms for sustaining large-scale (primordial) heterogeneity in the Earth's mantle using numerical models. Mantle flow is controlled by rock density and viscosity. Variations in intrinsic rock density, such as due to heterogeneity in basalt or iron content, can induce layering or partial layering in the mantle. Layering can be sustained in the presence of persistent whole mantle convection due to active "unmixing" of heterogeneity in low-viscosity domains, e.g. in the transition zone or near the core-mantle boundary [1]. On the other hand, lateral variations in intrinsic rock viscosity, such as due to heterogeneity in Mg/Si, can strongly affect the mixing timescales of the mantle. In the extreme case, intrinsically strong rocks may remain unmixed through the age of the Earth, and persist as large-scale domains in the mid-mantle due to focusing of deformation along weak conveyor belts [2]. That large-scale lateral heterogeneity and/or layering can persist in the presence of whole-mantle convection can explain the stagnation of some slabs, as well as the deflection of some plumes, in the mid-mantle. These findings indeed motivate new seismic studies for rigorous testing of model predictions. [1] Ballmer, M. D., N. C. Schmerr, T. Nakagawa, and J. Ritsema (2015), Science Advances, doi:10.1126/sciadv.1500815. [2] Ballmer, M. D., C. Houser, J. W. Hernlund, R. Wentzcovitch, and K. Hirose (2017), Nature Geoscience, doi:10.1038/ngeo2898.

Evidence of lower-mantle slab penetration phases in plate motions.

PubMed

Goes, Saskia; Capitanio, Fabio A; Morra, Gabriele

2008-02-21

It is well accepted that subduction of the cold lithosphere is a crucial component of the Earth's plate tectonic style of mantle convection. But whether and how subducting plates penetrate into the lower mantle is the subject of continuing debate, which has substantial implications for the chemical and thermal evolution of the mantle. Here we identify lower-mantle slab penetration events by comparing Cenozoic plate motions at the Earth's main subduction zones with motions predicted by fully dynamic models of the upper-mantle phase of subduction, driven solely by downgoing plate density. Whereas subduction of older, intrinsically denser, lithosphere occurs at rates consistent with the model, younger lithosphere (of ages less than about 60 Myr) often subducts up to two times faster, while trench motions are very low. We conclude that the most likely explanation is that older lithosphere, subducting under significant trench retreat, tends to lie down flat above the transition to the high-viscosity lower mantle, whereas younger lithosphere, which is less able to drive trench retreat and deforms more readily, buckles and thickens. Slab thickening enhances buoyancy (volume times density) and thereby Stokes sinking velocity, thus facilitating fast lower-mantle penetration. Such an interpretation is consistent with seismic images of the distribution of subducted material in upper and lower mantle. Thus we identify a direct expression of time-dependent flow between the upper and lower mantle.

Lithospheric radial anisotropy beneath the Gulf of Mexico

NASA Astrophysics Data System (ADS)

Chu, Risheng; Ko, Justin Yen-Ting; Wei, Shengji; Zhan, Zhongwen; Helmberger, Don

2017-05-01

The Lithosphere-Asthenosphere Boundary (LAB), where a layer of low viscosity asthenosphere decouples with the upper plate motion, plays an essential role in plate tectonics. Most dynamic modeling assumes that the shear velocity can be used as a surrogate for viscosity which provides key information about mantle flow. Here, we derive a shear velocity model for the LAB structure beneath the Gulf of Mexico allowing a detailed comparison with that beneath the Pacific (PAC) and Atlantic (ATL). Our study takes advantage of the USArray data from the March 25th, 2013 Guatemala earthquake at a depth of 200 km. Such data is unique in that we can observe a direct upward traveling lid arrival which remains the first arrival ahead of the triplications beyond 18°. This extra feature in conjunction with upper-mantle triplication sampling allows good depth control of the LAB and a new upper-mantle seismic model ATM, a modification of ATL, to be developed. ATM has a prominent low velocity zone similar to the structure beneath the western Atlantic. The model contains strong radial anisotropy in the lid where VSH is about 6% faster than VSV. This anisotropic feature ends at the bottom of the lithosphere at about the depth of 175 km in contrast to the Pacific where it extends to over 300 km. Another important feature of ATM is the weaker velocity gradient from the depth of 175 to 350 km compared to Pacific models, which may be related to differences in mantle flow.

Mantle-circulation models with sequential data assimilation: inferring present-day mantle structure from plate-motion histories.

PubMed

Bunge, Hans-Peter; Richards, M A; Baumgardner, J R

2002-11-15

Data assimilation is an approach to studying geodynamic models consistent simultaneously with observables and the governing equations of mantle flow. Such an approach is essential in mantle circulation models, where we seek to constrain an unknown initial condition some time in the past, and thus cannot hope to use first-principles convection calculations to infer the flow history of the mantle. One of the most important observables for mantle-flow history comes from models of Mesozoic and Cenozoic plate motion that provide constraints not only on the surface velocity of the mantle but also on the evolution of internal mantle-buoyancy forces due to subducted oceanic slabs. Here we present five mantle circulation models with an assimilated plate-motion history spanning the past 120 Myr, a time period for which reliable plate-motion reconstructions are available. All models agree well with upper- and mid-mantle heterogeneity imaged by seismic tomography. A simple standard model of whole-mantle convection, including a factor 40 viscosity increase from the upper to the lower mantle and predominantly internal heat generation, reveals downwellings related to Farallon and Tethys subduction. Adding 35% bottom heating from the core has the predictable effect of producing prominent high-temperature anomalies and a strong thermal boundary layer at the base of the mantle. Significantly delaying mantle flow through the transition zone either by modelling the dynamic effects of an endothermic phase reaction or by including a steep, factor 100, viscosity rise from the upper to the lower mantle results in substantial transition-zone heterogeneity, enhanced by the effects of trench migration implicit in the assimilated plate-motion history. An expected result is the failure to account for heterogeneity structure in the deepest mantle below 1500 km, which is influenced by Jurassic plate motions and thus cannot be modelled from sequential assimilation of plate motion histories limited in age to the Cretaceous. This result implies that sequential assimilation of past plate-motion models is ineffective in studying the temporal evolution of core-mantle-boundary heterogeneity, and that a method for extrapolating present-day information backwards in time is required. For short time periods (of the order of perhaps a few tens of Myr) such a method exists in the form of crude 'backward' convection calculations. For longer time periods (of the order of a mantle overturn), a rigorous approach to extrapolating information back in time exists in the form of iterative nonlinear optimization methods that carry assimilated information into the past through the use of an adjoint mantle convection model.

Postseismic relaxation process and lithospheric rheology inferred from eight years of postseismic deformation after the 2008 Mw7.9 Wenchuan earthquake

NASA Astrophysics Data System (ADS)

Zhao, B.; Burgmann, R.; Rui, X.; Wang, D.; Yu, J.; He, K.

2017-12-01

Current inferences of postseismic deformation mechanisms and lithospheric rheology in the eastern Tibetan Plateau strongly depend on spatial and temporal observations of postseismic transients following the 2008 Mw=7.9 Wenchuan earthquake. We processed regional continuously operating and survey-mode GPS data from the Crustal Movement Observation Network of China and Sichuan Continuous Operation Reference System. These data cover a broad region and time intervals of up to eight years. The determined amplitude of postseismic displacements show clear contrast between the Sichuan Basin and eastern Tibet. In addition to significant amounts of deformation in the region between the Longmen Shan and Longriba fault, reliable deformation transients are also visible in the far field, such as regions to the west of the Longriba fault and along the left-lateral Xianshuihe fault. In contrast, no more than 10 mm of postseismic transients are observed in the Sichuan Basin. Guided by previous studies, we conducted multiple-mechanism models of afterslip and viscoelastic relaxation. We first explored a series of forward viscoelastic relaxation models using a heterogeneous rheological earth structure, and then inverted corresponding afterslip distributions on the shallowly dipping detachment to explain the remaining residuals. Our preliminary results indicate the viscoelastic relaxation in the lower crust and upper mantle dominantly contributed to the mid- and far-field observations, whereas afterslip below the coseismic asperities and on small patches near the surface can explain the near-field measurements. Time-dependent slip inversions illustrate that afterslip decays more rapidly on the shallow portions of the fault interface than on the shallowly dipping detachment. Relatively long-lived right-lateral afterslip is revealed in the north segment of the Beichuan fault, suggesting variations of frictional properties along strike of the fault zone. Our results also support previous inferences of higher mantle viscosities below the Sichuan Basin and lower viscosities of the lower crust and upper mantle below eastern Tibet. The transient and steady-state viscosities of Tibet's lower crust are constrained to be 1018 and 1019 Pa s. The upper mantle viscosity is poorly resolved due to small coseismic stress change.

A global geochemical model for the evolution of the mantle

NASA Technical Reports Server (NTRS)

Anderson, D. L.

1979-01-01

It is proposed that the upper mantle transition region, 220 to 670 km, is composed of eclogite which has been derived from primitive mantle by about 20 percent partial melting and that this is the source and sink of oceanic crust. The remainder of the upper mantle is garnet peridotite which is the source of continental basalts and hotspot magmas. This region is enriched in incompatible elements by hydrous and CO2 rich metasomatic fluids which have depleted the underlying layers in the L.I.L. elements and L.R.E.E. The volatiles make this a low-velocity, high attenuation, low viscosity region. The eclogite layer is internally heated and its controls the convection pattern in the upper mantle. Plate tectonics is intermittent. The continental thermal anomaly at a depth of 150-220 km triggers kimberlite and carbonatite activity, alkali and flood basalt volcanism, vertical tectonics and continental breakup. Hot spots remain active after the continents leave and build the oceanic islands. Mantle plumes rise from a depth of about 220 km. Midocean ridge basalts rise from the depleted layer below this depth. Material from this layer can also be displaced upwards by subducted oceanic lithosphere to form back-arc basins.

Evolution of the viscosity of Earth's upper mantle: Grain-boundary sliding and the role of microstructure in olivine deformation

NASA Astrophysics Data System (ADS)

Hansen, Lars N.

Many features of plate tectonics cannot be explained with standard rheological models of the upper mantle. In particular, the localization of deformation at plate boundaries requires the viscosity of the constituent rocks to evolve spatially and temporally. Such rheological complexity may arise from changing microstructural state variables (e.g., grain size and crystallographic-fabric strength), but the degree to which microstructure contributes to the evolution of viscosity is unclear given our current understanding of deformation mechanisms in mantle minerals. Dislocation-accommodated grain-boundary sliding (GBS) is a potentially critical mechanism for localizing deformation in olivine because it imparts a sensitivity of the viscosity to the state of the microstructure while simultaneously providing mechanisms for changing the microstructure. However, many details of GBS in olivine are currently unknown including 1) the magnitude of the sensitivity of strain rate to crystallographic fabric and grain size, 2) the strength of the crystallographic fabrics produced, and 3) the anisotropy in viscosity of polycrystalline aggregates. Detailed knowledge of these unknowns is necessary to assess the importance of microstructural evolution in the operation of plate tectonics. This dissertation investigates the details of GBS in olivine through four sets of laboratory-based experiments. In Chapter 2, triaxial compressive creep experiments on aggregates of San Carlos olivine are used to develop a flow law for olivine deforming by GBS. Extrapolations of strain rate to geological conditions using the derived flow law indicate that GBS is the dominant deformation mechanism throughout the uppermost mantle. Crystallographic fabrics observed in deformed samples are consistent with upper-mantle seismic anisotropy. In Chapter 3, torsion experiments on iron-rich olivine are used to determine the rheological behavior of olivine deforming by GBS at large strains. The sensitivity of the strain rate to grain size and stress is demonstrated to be consistent with low-strain experiments. Additionally, the sensitivity of strain rate to the development of a crystallographic fabric is determined. Constitutive relationships including microstructural evolution are developed that accurately predict the observed stress as a function of strain. The results of Chapter 3 confirm that significant weakening is associated with both grain-size reduction and crystallographic-fabric development. In Chapter 4, torsion experiments on iron-rich olivine are used to determine if microstructural evolution can lead to strain localization. Experiments were conducted with either constant-strain-rate or constant-stress boundary conditions. Localization is only observed in samples deformed at constant-stress, which suggests boundary conditions affect the critical size of strength perturbation necessary for localization to occur. Strain localization is correlated with fine-grained regions, and a feedback mechanism between grain-size reduction and strain rate is proposed. In Chapter 5, both torsion and tension experiments are used to assess the mechanical anisotropy of previously deformed samples. Based on the direction of the applied stress relative to the orientation of a pre-existing crystallographic fabric, the viscosity is demonstrated to vary by over an order of magnitude. This observation suggests deformation can localize in regions that were previously deformed and retained a strong crystallographic fabric. The results of this dissertation elucidate the interplay between microstructure and deformation of olivine in the GBS regime. Because the viscosity of olivine-rich rocks deforming by GBS is dependent on both grain size and crystallographic fabric, heterogeneities in these microstructural parameters can lead to spatial and temporal variations in viscosity, possibly explaining the large-scale patterns of deformation in the upper mantle. Future numerical simulations can test the importance of microstructure in geodynamic processes by incorporating the constitutive relationships outlined in this dissertation.

Quantifying mantle structure and dynamics using plume tracing in seismic tomography

NASA Astrophysics Data System (ADS)

O'Farrell, K. A.; Eakin, C. M.; Jackson, M. G.; Jones, T. D.; Lekic, V.; Lithgow-Bertelloni, C. R.

2017-12-01

Directly linking deep mantle processes with surface features and dynamics is a complex problem. Hotspot volcanism gives us surface observables of mantle signatures, but the depth and source of the mantle plumes feeding these hotspots are highly debated. To address these issues, it is necessary to consider the entire journey of a plume through the mantle. By analyzing the behavior of mantle plumes we can constrain the vigor of mantle convection, the net rotation of the mantle and the role of thermal versus chemical anomalies as well as the bulk physical properties such as the viscosity profile. To do this, we developed a new algorithm to trace plume-like features in shear-wave (Vs) seismic tomography models based on picking local minima in the velocity and searching for continuous features with depth. We applied this method to recent tomographic models and find 60+ continuous plume conduits that are > 750 km long. Approximately a third of these can be associated with known hotspots at the surface. We analyze the morphology of these continuous conduits and infer large scale mantle flow patterns and properties. We find the largest lateral deflections in the conduits occur near the base of the lower mantle and in the upper mantle (near the thermal boundary layers). The preferred orientation of the plume deflections show large variability at all depths and indicate no net mantle rotation. Plate by plate analysis shows little agreement in deflection below particular plates, indicating these deflected features might be long lived and not caused by plate shearing. Changes in the gradient of plume deflection are inferred to correspond with viscosity contrasts in the mantle and found below the transition zone as well as at 1000 km depth. From this inferred viscosity structure, we explore the dynamics of a plume through these viscosity jumps. We also retrieve the Vs profiles for the conduits and compare with the velocity profiles predicted for different mantle adiabat temperatures. We are able to constrain the average temperature anomaly of the conduits to be around 150 K. We use these thermal anomalies in conjunction with our measured plume tilts/deflections to further explore the dynamics of plume conduits in the lower mantle and transition zone.

The dynamics of plate tectonics and mantle flow: from local to global scales.

PubMed

Stadler, Georg; Gurnis, Michael; Burstedde, Carsten; Wilcox, Lucas C; Alisic, Laura; Ghattas, Omar

2010-08-27

Plate tectonics is regulated by driving and resisting forces concentrated at plate boundaries, but observationally constrained high-resolution models of global mantle flow remain a computational challenge. We capitalized on advances in adaptive mesh refinement algorithms on parallel computers to simulate global mantle flow by incorporating plate motions, with individual plate margins resolved down to a scale of 1 kilometer. Back-arc extension and slab rollback are emergent consequences of slab descent in the upper mantle. Cold thermal anomalies within the lower mantle couple into oceanic plates through narrow high-viscosity slabs, altering the velocity of oceanic plates. Viscous dissipation within the bending lithosphere at trenches amounts to approximately 5 to 20% of the total dissipation through the entire lithosphere and mantle.

Production and recycling of oceanic crust in the early Earth

NASA Astrophysics Data System (ADS)

van Thienen, P.; van den Berg, A. P.; Vlaar, N. J.

2004-08-01

Because of the strongly different conditions in the mantle of the early Earth regarding temperature and viscosity, present-day geodynamics cannot simply be extrapolated back to the early history of the Earth. We use numerical thermochemical convection models including partial melting and a simple mechanism for melt segregation and oceanic crust production to investigate an alternative suite of dynamics which may have been in operation in the early Earth. Our modelling results show three processes that may have played an important role in the production and recycling of oceanic crust: (1) Small-scale ( x×100 km) convection involving the lower crust and shallow upper mantle. Partial melting and thus crustal production takes place in the upwelling limb and delamination of the eclogitic lower crust in the downwelling limb. (2) Large-scale resurfacing events in which (nearly) the complete crust sinks into the (eventually lower) mantle, thereby forming a stable reservoir enriched in incompatible elements in the deep mantle. New crust is simultaneously formed at the surface from segregating melt. (3) Intrusion of lower mantle diapirs with a high excess temperature (about 250 K) into the upper mantle, causing massive melting and crustal growth. This allows for plumes in the Archean upper mantle with a much higher excess temperature than previously expected from theoretical considerations.

Asthenosphere rheology inferred from observations of the 2012 Indian Ocean earthquake.

PubMed

Hu, Yan; Bürgmann, Roland; Banerjee, Paramesh; Feng, Lujia; Hill, Emma M; Ito, Takeo; Tabei, Takao; Wang, Kelin

2016-10-20

The concept of a weak asthenospheric layer underlying Earth's mobile tectonic plates is fundamental to our understanding of mantle convection and plate tectonics. However, little is known about the mechanical properties of the asthenosphere (the part of the upper mantle below the lithosphere) underlying the oceanic crust, which covers about 60 per cent of Earth's surface. Great earthquakes cause large coseismic crustal deformation in areas hundreds of kilometres away from and below the rupture area. Subsequent relaxation of the earthquake-induced stresses in the viscoelastic upper mantle leads to prolonged postseismic crustal deformation that may last several decades and can be recorded with geodetic methods. The observed postseismic deformation helps us to understand the rheological properties of the upper mantle, but so far such measurements have been limited to continental-plate boundary zones. Here we consider the postseismic deformation of the very large (moment magnitude 8.6) 2012 Indian Ocean earthquake to provide by far the most direct constraint on the structure of oceanic mantle rheology. In the first three years after the Indian Ocean earthquake, 37 continuous Global Navigation Satellite Systems stations in the region underwent horizontal northeastward displacements of up to 17 centimetres in a direction similar to that of the coseismic offsets. However, a few stations close to the rupture area that had experienced subsidence of up to about 4 centimetres during the earthquake rose by nearly 7 centimetres after the earthquake. Our three-dimensional viscoelastic finite-element models of the post-earthquake deformation show that a thin (30-200 kilometres), low-viscosity (having a steady-state Maxwell viscosity of (0.5-10) × 10 18 pascal seconds) asthenospheric layer beneath the elastic oceanic lithosphere is required to produce the observed postseismic uplift.

The temporal evolution of a subducting plate in the lower mantle

NASA Astrophysics Data System (ADS)

Loiselet, C.; Grujic, D.; Braun, J.; Fullsack, P.; Thieulot, C.; Yamato, P.

2009-04-01

It is now widely accepted that some subducting slabs may cross the lower/upper mantle boundary to ground below the 660 km discontinuity. Indeed, geophysical data underline long and narrow traces of fast materials, associated with subducting slabs, from the upper mantle transition zone to mid-mantle depths that are visible beneath North and South America and southern Asia (Li et al, 2008). Furthermore, seismic tomography data (Van der Hilst et al., 1997; Karason and van der Hilst, 2000, 2001) show a large variety of slab geometries and of mantle flow patterns around subducting plate boundaries (e.g. the slab geometry in the lower mantle in the Tonga subduction zone). However, seismic tomography does not elucidate the temporal evolution of the slab behaviour and geometry during its descent through the upper and lower mantle. In this work, we therefore propose to study the deformation of a thin plate (slab) falling in a viscous fluid (mantle) by means of both analogue and numerical modelling. The combination of both analogue and numerical experiments provides important insights into the shape and attitude evolution of subducting slabs. Models bring information into the controls exerted by the rheology of the slab and the mantle and other physical parameters such as the density contrast between the slab and the surrounding mantle, on the rate at which this deformation takes place. We show that in function of a viscosity ratios between the plate and the surrounding fluid, the plate will acquire a characteristic shape. For the isoviscous case, the plate shape tends toward a bubble with long tails: a "jellyfish" form. The time necessary for the plate to acquire this shape is a function of the viscosity and density contrast between the slab and the mantle. To complete our approach, we have developed a semi-analytical model based on the solution of the Hadamar-Rybinski equations for the problem of a dense, yet isoviscous and thus deforming sphere. This model helps to better describe flow processes around the downgoing plate and, simultaneously, to characterize its deformation. In this way, we were able to calculate the velocities in the mantle, the forces exerted by the fluid on the plate, and the dissipated energy in the surrounding fluid. Experimental results will be correlated with geophysical data.

The Temporal Evolution Of A Subducting Plate In The Lower Mantle

NASA Astrophysics Data System (ADS)

Loiselet, C.; Grujic, D.; Fullsack, P.; Thieulot, C.; Yamato, P.; Braun, J.

2008-12-01

It is now widely accepted that some subducting slabs may cross the lower/upper mantle boundary to ground below the 660 km discontinuity. Indeed, geophysical data underline long and narrow traces of fast materials, associated with subducting slabs, from the upper mantle transition zone to mid-mantle depths that are visible beneath North and South America and southern Asia (Li et al, 2008). Furthermore, seismic tomography data (Van der Hilst et al., 1997; Karason and van der Hilst, 2000, 2001)) show a large variety of slab geometries and of mantle flow patterns around subducting plate boundaries (e.g. the slab geometry in the lower mantle in the Tonga subduction zone). However, seismic tomography does not elucidate the temporal evolution of the slab behaviour and geometry during its descent through the upper and lower mantle. In this work, we therefore propose to study the deformation of a thin plate (slab) falling in a viscous fluid (mantle). The combination of both analogue and numerical experiments provides important insights into the shape and attitude evolution of subducting slabs. Models bring information into the controls exerted by the rheology of the slab and the mantle and other physical parameters such as the density contrast between the slab and the surrounding mantle, on the rate at which this deformation takes place. We show that in function of a viscosity ratios between the plate and the surrounding fluid, the plate will acquire a characteristic shape. For the isoviscous case, the plate shape tends toward a bubble with long tails: a jellyfish form. The time necessary for the plate to acquire this shape is a function of the viscosity and density contrast between the slab and the mantle. To complete our approach, we have developed a semi-analytical model based on the solution of the Hadamar-Rybinski equations for the problem of a dense, yet isoviscous and thus deforming sphere. This model helps to better describe flow processes around the downgoing plate and, simultaneously, to characterize its deformation. In this way, we were able to calculate the velocities in the mantle, the forces exerted by the fluid on the plate, and the dissipated energy in the surrounding fluid. Experimental results will be correlated with geophysical data.

Driving forces: Slab subduction and mantle convection

NASA Technical Reports Server (NTRS)

Hager, Bradford H.

1988-01-01

Mantle convection is the mechanism ultimately responsible for most geological activity at Earth's surface. To zeroth order, the lithosphere is the cold outer thermal boundary layer of the convecting mantle. Subduction of cold dense lithosphere provides tha major source of negative buoyancy driving mantle convection and, hence, surface tectonics. There are, however, importnat differences between plate tectonics and the more familiar convecting systems observed in the laboratory. Most important, the temperature dependence of the effective viscosity of mantle rocks makes the thermal boundary layer mechanically strong, leading to nearly rigid plates. This strength stabilizes the cold boundary layer against small amplitude perturbations and allows it to store substantial gravitational potential energy. Paradoxically, through going faults at subduction zones make the lithosphere there locally weak, allowing rapid convergence, unlike what is observed in laboratory experiments using fluids with temperature dependent viscosities. This bimodal strength distribution of the lithosphere distinguishes plate tectonics from simple convection experiments. In addition, Earth has a buoyant, relatively weak layer (the crust) occupying the upper part of the thermal boundary layer. Phase changes lead to extra sources of heat and bouyancy. These phenomena lead to observed richness of behavior of the plate tectonic style of mantle convection.

Viscosity structure of Earth's mantle inferred from rotational variations due to GIA process and recent melting events

NASA Astrophysics Data System (ADS)

Nakada, Masao; Okuno, Jun'ichi; Lambeck, Kurt; Purcell, Anthony

2015-08-01

We examine the geodetically derived rotational variations for the rate of change of degree-two harmonics of Earth's geopotential, skew5dot J_2, and true polar wander, combining a recent melting model of glaciers and the Greenland and Antarctic ice sheets taken from the IPCC 2013 Report (AR5) with two representative GIA ice models describing the last deglaciation, ICE5G and the ANU model developed at the Australian National University. Geodetically derived observations of skew4dot J_2 are characterized by temporal changes of -(3.7 ± 0.1) × 10-11 yr-1 for the period 1976-1990 and -(0.3 ± 0.1) × 10-11 yr-1 after ˜2000. The AR5 results make it possible to evaluate the recent melting of the major ice sheets and glaciers for three periods, 1900-1990, 1991-2001 and after 2002. The observed skew4dot J_2 and the component of skew4dot J_2 due to recent melting for different periods indicate a long-term change in skew4dot J_2-attributed to the Earth's response to the last glacial cycle-of -(6.0-6.5) × 10-11 yr-1, significantly different from the values adopted to infer the viscosity structure of the mantle in most previous studies. This is a main conclusion of this study. We next compare this estimate with the values of skew4dot J_2 predicted by GIA ice models to infer the viscosity structure of the mantle, and consequently obtain two permissible solutions for the lower mantle viscosity (ηlm), ˜1022 and (5-10) × 1022 Pa s, for both adopted ice models. These two solutions are largely insensitive to the lithospheric thickness and upper mantle viscosity as indicated by previous studies and relatively insensitive to the viscosity structure of the D″ layer. The ESL contributions from the Antarctic ice sheet since the last glacial maximum (LGM) for ICE5G and ANU are about 20 and 30 m, respectively, but glaciological reconstructions of the Antarctic LGM ice sheet have suggested that its ESL contribution may have been less than ˜10 m. The GIA-induced skew4dot J_2 for GIA ice models with an Antarctic ESL component of ˜10 m suggests two permissible lower mantle viscosity solutions of ηlm ˜ 2 × 1022 and ˜5 × 1022 Pa s or one solution with (2-5) × 1022 Pa s. These results suggest that the effective lower mantle viscosity is larger than ˜1022 Pa s regardless of the uncertainties for an Antarctic ESL component. We also examine the polar wander due to recent melting and GIA processes, suggesting that the observed polar wander may be significantly attributed to convection motions in the mantle and/or another cause, particularly for permissible lower mantle viscosity solution of (5-10) × 1022 Pa s.

Correction to “Changes in the Earth's rotation by tectonic movements”

NASA Astrophysics Data System (ADS)

Vermeersen, L. L. A.; Vlaar, N. J.

1993-06-01

Present-day true polar wander and the secular non-tidal acceleration of the earth are usually attributed to post-glacial rebound. In the models which relate this rebound to changes in rotation, the mantle is assumed to relax passively to the melted ice-loads. The lithosphere is usually modeled as a highly viscous upper layer in these models, having viscosities which exceed mantle viscosities by several orders of magnitude. We propose that lithospheric processes unrelated to post-glacial rebound and taking place under non-isostatic conditions are also able to induce non-negligible influences on the earth's rotation. Examples of such processes are mountain building and erosion, foundering flexure of oceanic basins and lithospheric snapbacking resulting from detachment of subducting slabs. Lithospheric and crustal rheologies and intraplate-stresses are the dominant factors in these mechanisms, contrary to the mantle rheologies which are assumed to dominate the process of post-glacial rebound.

Experimental Deformation of Olivine Crystals at Mantle P and T: Evidences for a Pressure-Induced Slip Transition and Implications for Upper-Mantle Seismic Anisotropy and Low Viscosity Zone

NASA Astrophysics Data System (ADS)

Raterron, P.; Chen, J.; Geenen, T.; Girard, J.

2009-04-01

Recent developments in high-pressure deformation devices coupled with synchrotron radiation allow investigating the rheology of mantle minerals and aggregates at the extreme pressure (P) and temperature (T) of their natural occurrence in the Earth. This is particularly true in the case of olivine, which rheology has been recently investigated in the Deformation-DIA apparatus (D-DIA, see Wang et al., 2003, Rev. Scientific Instr., 74, 3002) at upper-mantle P and T conditions. Olivine deforms by dislocation creep in the shallow upper-mantle, as revealed by the seismic velocity anisotropy observed in this region. The attenuation of seismic anisotropy at depth greater than 200 km is interpreted as a pressure-induced change in olivine main deformation mechanism. It was first attributed to a transition from dislocation creep to diffusion creep (Karato and Wu, 1993, Science, 260, 771). This interpretation has been challenged by deformation data obtained at high pressure (P > 3 GPa) in the dislocation creep regime (Couvy et al., 2004, EJM, 16, 877; Raterron et al., 2007, Am. Miner., 92, 1436; Raterron et al., 2009, PEPI, 72, 74), which support a second interpretation: a transition in olivine dominant dislocation slip, from [100] slip at low P to [001] slip at high P (e.g., Mainprice et al., 2005, Nature, 433, 731). Such a P -induced [100]/[001] slip transition is also supported by recent theoretical studies based on first-principle calculations of olivine dislocation slips (Durinck et al., 2005, PCM, 32, 646; Durinck et al., 2007, Eur. J. Mineral., 19, 631). In order to further constrain the effect of pressure on olivine slip system activities, deformation experiments were carried out in poor water condition at P > 5 GPa and T =1400˚ C, on pure forsterite (Fo100) and San Carlos olivine crystals, using the D-DIA at the X17B2 beamline of the NSLS (Upton, NY, USA). Crystals were oriented in order to active either [100] slip alone or [001] slip alone in (010) plane, or both [100](001) and [001](100) systems together. Constant applied stress < 300 MPa and specimen strain rates were monitored in situ using time-resolved X-ray diffraction and radiography, respectively. Run products were investigated by transmission electron microscopy (TEM) in order to verify the actual activation of the tested dislocation slip systems. The obtained data were compared with rheological data previously obtained at comparable T and conditions, but at room P (Darot and Gueguen, 1981, JGR, 86, 6219; Bai et al., 1991, JGR, 96, 2441), resulting in creep power laws which quantify the effect of P on olivine rheology. The new data confirm the occurrence of a P -induced [100]/[001] slip transition, and suggest that [001](010) system dominates olivine deformation in the deep upper mantle. Extrapolation of the obtained rheological laws to natural condition along upper-mantle geotherms, shows that the [100] / [001] slip transition should occur in the Earth at ~ 200 km depth, thus can explain the attenuation of seismic anisotropy in the deep upper mantle. The obtained rheological laws were also integrated into a straightforward olivine aggregate model, then extrapolated to mantle condition using a 2-D geodynamic modeling application (Van den Berg et al., 1993, Geophys. J. International, 115, 62), which is the simplest approach to investigate upper-mantle steady-state deformation. In the application, the velocity of the lower boundary (the transition-zone boundary at 410-km depth) was set to 0, while that at the Earth's surface was set to 2 cm/year. Results from this modeling suggest that the combine activity of [100] and [001] slips in olivine aggregates may significantly decrease mantle viscosity below the oceanic lithosphere, thus, may contribute to the low viscosity zone (LVZ) required in plate tectonics to decouple rigid plates from the more ductile asthenophere underneath.

Dynamical influences on the moment of inertia tensor from lateral viscosity variations inferred from seismic tomographic models

NASA Technical Reports Server (NTRS)

Zhang, Shuxia; Yuen, David A.

1994-01-01

We have investigated the influences of lateral variations of viscosity on the moment of inertia tensor from viscous flows due to the density anomalies in the mantle inferred from seismic tomographic models. The scaling relations between the density and the seismic anomalies is taken as either a constant or a function increasing with depth in accord with the recent high-pressure experimental studies. The viscosity is taken as an exponential function of the 3D density anomaly. In models with an isoviscous background, the effects on the perturbed moment of inertia tensor from the lateral viscosity variations are smaller than those due to variations in the radial viscosity profiles. In mantle models with a background viscosity increasing with depth, the influences of the lateral viscosity variations are significant. The most striking feature in the latter case is that the two off-diagonal elements delta I(sub xz) and delta I(sub yz) in the inertia tensor exhibit greatest sensitivity to lateral variations of the viscosity. While the other elements of the inertia change by only about a few tens of percent in the range of lateral viscosity contrast considered (less than 300), delta I(sub xz) and delta I(sub yz) can vary up to 40 times even with a change in sign, depending on the radial viscosity stratification and the location of the strongest lateral variations. The increase in the velocity-density scaling relation with depth can reduce the influences of the lateral viscosity variations, but it does not change the overall sensitive nature of delta I(sub xz) and delta I(sub yz). This study demonstrates clearly that the lateral viscosity variations, especially in the upper mantle, must be considered in the determination of long-term polar wander, since the variations in the delta I(sub xz) and delta I(sub yz) terms are directly responsible for exciting rotational movements.

Rheological Structure of Northern Tibet Lithosphere Inferred from Postseismic Deformation Modeling of 2001 Mw7.8 Kokoxili Earthquake

NASA Astrophysics Data System (ADS)

Shen, Z.; Wang, M.; He, P.

2017-12-01

It has been debated for decades about the rheological structure of Tibetan plateau lithosphere. Occurrence of the 2001 Mw7.8 Kokoxili earthquake greatly changed the tectonic stress field in northern Tibet, and offered a rare opportunity to infer the rheological structure of lithosphere by modeling postseismic deformation observed at the Earth's surface. We collect GPS data from 45 sites observed after the quake, most of them are within 100 km distance from the fault, and the data are processed to produce the station position time series. Data observed after the 2001 Kokoxili and prior to the 2008 Yutian earthquake are used, with at most 6.4 years of observation time span. A joint inversion is performed to solve for the viscous relaxation in lithosphere and afterslip on fault simultaneously. The Bayankala-Qiangtang region and the Qaidam Basin, located south and north of the fault, are assumed to have the Burgers-body rheological structure same in the lower crust and upper mantle but different across the fault. Viscosities of the two regions are inverted through a grid search procedure. Afterslip is constrained by both of the GPS observations and the a priori Coulomb stress distribution on fault. Our inversion result shows the transient viscosities of 1.5*1018 Pa·s and 5*1018 Pa·s for the regions south and north of the East Kunlun fault, respectively. Secular viscosities of 1.5*1019 Pa·s and 1.5*1020 Pa·s are found for the regions south and north of the fault. The secular viscosities of the lower crust and upper mantle are found not differentiable statistically from the above values for both regions north and south of the fault. Our result reveals that viscosities of the lower crust and upper mantle underneath the Bayankala-Qiangtang region is much lower than that underneath the Qaidam Basin; and the result is consistent with a mechanic model that the Bayankala-Qiangtang block is relatively weak and undergoes distributed deformation, and the Qaidam block is relatively strong and undergoes block-like deformation. Viscosities of the lower crust in northern Tibet is 2-3 orders of magnitude higher than that a lower crust flow model would demand, suggesting that lower crust flow may not exist within this region.

Dominant Controls of Downdip Afterslip and Viscous Relaxation on the Postseismic Displacements Following the Mw7.9 Gorkha, Nepal, Earthquake

NASA Astrophysics Data System (ADS)

Zhao, Bin; Bürgmann, Roland; Wang, Dongzhen; Tan, Kai; Du, Ruilin; Zhang, Rui

2017-10-01

We analyze three-dimensional GPS coordinate time series from continuously operating stations in Nepal and South Tibet and calculate the initial 1 year postseismic displacements. We first investigate models of poroelastic rebound, afterslip, and viscoelastic relaxation individually and then attempt to resolve the trade-offs between their contributions by evaluating the misfit between observed and simulated displacements. We compare kinematic inversions for distributed afterslip with stress-driven afterslip models. The modeling results show that no single mechanism satisfactorily explains near- and far-field postseismic deformation following the Gorkha earthquake. When considering contributions from all three mechanisms, we favor a combination of viscoelastic relaxation and afterslip alone, as poroelastic rebound always worsens the misfit. The combined model does not improve the data misfit significantly, but the inverted afterslip distribution is more physically plausible. The inverted afterslip favors slip within the brittle-ductile transition zone downdip of the coseismic rupture and fills the small gap between the mainshock and largest aftershock slip zone, releasing only 7% of the coseismic moment. Our preferred model also illuminates the laterally heterogeneous rheological structure between India and the South Tibet. The transient and steady state viscosities of the upper mantle beneath Tibet are constrained to be greater than 1018 Pa s and 1019 Pa s, whereas the Indian upper mantle has a high viscosity ≥1020 Pa s. The viscosity in the lower crust of southern Tibet shows a clear trade-off with its southward extent and thickness, suggesting an upper bound value of 8 × 1019 Pa s for its steady state viscosity.

Geoid, topography, and convection-driven crustal deformation on Venus

NASA Technical Reports Server (NTRS)

Simons, Mark; Hager, Bradford H.; Solomon, Sean C.

1992-01-01

High-resolution Magellan images and altimetry of Venus reveal a wide range of styles and scales of surface deformation that cannot readily be explained within the classical terrestrial plate tectonic paradigm. The high correlation of long-wavelength topography and gravity and the large apparent depths of compensation suggest that Venus lacks an upper-mantle low-viscosity zone. A key difference between Earth and Venus may be the degree of coupling between the convecting mantle and the overlying lithosphere. Mantle flow should then have recognizable signatures in the relationships between surface topography, crustal deformation, and the observed gravity field.

Rheologic effects of crystal preferred orientation in upper mantle flow near plate boundaries

NASA Astrophysics Data System (ADS)

Blackman, Donna; Castelnau, Olivier; Dawson, Paul; Boyce, Donald

2016-04-01

Observations of anisotropy provide insight into upper mantle processes. Flow-induced mineral alignment provides a link between mantle deformation patterns and seismic anisotropy. Our study focuses on the rheologic effects of crystal preferred orientation (CPO), which develops during mantle flow, in order to assess whether corresponding anisotropic viscosity could significantly impact the pattern of flow. We employ a coupled nonlinear numerical method to link CPO and the flow model via a local viscosity tensor field that quantifies the stress/strain-rate response of a textured mineral aggregate. For a given flow field, the CPO is computed along streamlines using a self-consistent texture model and is then used to update the viscosity tensor field. The new viscosity tensor field defines the local properties for the next flow computation. This iteration produces a coupled nonlinear model for which seismic signatures can be predicted. Results thus far confirm that CPO can impact flow pattern by altering rheology in directionally-dependent ways, particularly in regions of high flow gradient. Multiple iterations run for an initial, linear stress/strain-rate case (power law exponent n=1) converge to a flow field and CPO distribution that are modestly different from the reference, scalar viscosity case. Upwelling rates directly below the spreading axis are slightly reduced and flow is focused somewhat toward the axis. Predicted seismic anisotropy differences are modest. P-wave anisotropy is a few percent greater in the flow 'corner', near the spreading axis, below the lithosphere and extending 40-100 km off axis. Predicted S-wave splitting differences would be below seafloor measurement limits. Calculations with non-linear stress/strain-rate relation, which is more realistic for olivine, indicate that effects are stronger than for the linear case. For n=2-3, the distribution and strength of CPO for the first iteration are greater than for n=1, although the fast seismic axis directions are similar. The greatest difference in CPO for the nonlinear cases develop at the flow 'corner' at depths of 10-30 km and 20-100 km off-axis. J index values up to 10% greater than the linear case are predicted near the lithosphere base in that region. Viscosity tensor components are notably altered in the nonlinear cases. Iterations between the texture and flow calculations for the non-linear cases are underway this winter; results will be reported in the presentation.

Mantle convection patterns reveal the enigma of the Red Sea rifting

NASA Astrophysics Data System (ADS)

Petrunin, Alexey; Kaban, Mikhail; El Khrepy, Sami; Al-Arifi, Nassir

2017-04-01

Initiation and further development of the Red Sea rift (RSR) is usually associated with the Afar plume at the Oligocene-Miocene separating the Arabian plate from the rest of the continent. Usually, the RSR is divided into three parts with different geological, tectonic and geophysical characteristics, but the nature of this partitioning is still debatable. To understand origin and driving forces responsible for the tectonic partitioning of the RSR, we have developed a global mantle convection model based on the refined density model and viscosity distribution derived from tectonic, rheological and seismic data. The global density model of the upper mantle is refined for the Middle East based on the high-resolution 3D model (Kaban et al., 2016). This model based on a joint inversion of the residual gravity and residual topography provides much better constraints on the 3D density structure compared to the global model based on seismic tomography. The refined density model and the viscosity distribution based on a homologous temperature approach provide an initial setup for further numerical calculations. The present-day snapshot of the mantle convection is calculated by using the code ProSpher 3D that allows for strong lateral variations of viscosity (Petrunin et al., 2013). The setup includes weak plate boundaries, while the measured GPS velocities are used to constrain the solution. The resulting mantle flow patterns show clear distinctions among the mantle flow patterns below the three parts of the RSR. According to the modeling results, tectonics of the southern part of the Red Sea is mainly determined by the Afar plume and the Ethiopian rift opening. It is characterized by a divergent mantle flow, which is connected to the East African Rift activity. The rising mantle flow is traced down to the transition zone and continues in the lower mantle for a few thousand kilometers south-west of Afar. The hot mantle anomaly below the central part of the RSR can be explained either by the asthenospheric upwelling due to the Red Sea floor spreading or by a secondary plume rising from the transition zone. According to our model, there is no obvious evidence for a direct connection of the hot anomaly below the central part of the RSR and the Afar plume in the upper mantle. In the northern part of the RSR, we found the ridge-axis aligned downstream flow contradicting the hypothesis of the intra-continental rifting in this area. Likely, the tectonics of this area implies a complex interplay of the Dead Sea transform fault development and the Sinai and Mediterranean tectonics. Kaban, M. K., S. El Khrepy, N. Al-Arifi, M. Tesauro, and W. Stolk (2016), Three dimensional density model of the upper mantle in the Middle East: Interaction of diverse tectonic processes, J. Geophys. Res. Solid Earth, 121, doi:10.1002/2015JB012755. Petrunin, A. G.; Kaban, M. K.; Rogozhina, I.; Trubitsyn, V. (2013). Revising the spectral method as applied to modeling mantle dynamics. Geochemistry Geophysics Geosystems (G3), EDOC: 21048.

Illumination of rheological mantle heterogeneity by the M7.2 2010 El Mayor-Cucapah earthquake

USGS Publications Warehouse

Pollitz, Fred F.; Bürgmann, Roland; Thatcher, Wayne R.

2012-01-01

Major intracontinental strike-slip faults tend to mark boundaries between lithospheric blocks of contrasting mechanical properties along much of their length. Both crustal and mantle heterogeneities can form such boundaries, but the role of crustal versus mantle strength contrasts for localizing strain sufficiently to generate major faults remains unclear. Using the crustal velocity field observed through the Global Positioning System (GPS) in the epicentral area of the M7.2 2010 El Mayor-Cucapah earthquake, Baja California, we find that transient deformation observed after the event is anomalously small in areas of relatively high seismic velocity in the shallow upper mantle (∼50 km depth). This pattern is best explained with a laterally heterogeneous viscoelastic structure that mimics the seismic structure. The mantle of the Southern Colorado River Desert (SCRD) and Peninsular Ranges (PR), which bound the fault system to its east and west, respectively, have anomalously high viscosity and seismic velocity. We hypothesize that compared with the rest of the San Andreas fault (SAF) system to its north, the strike-slip fault system in northern Baja California is narrow because of the presence of the PR and SCRD high-viscosity regions which bound it.

Constraints from Earth's heat budget on mantle dynamics

NASA Astrophysics Data System (ADS)

Kellogg, L. H.; Ferrachat, S.

2006-12-01

Recent years have seen an increase in the number of proposed models to explain Earth's mantle dynamics: while two end-members, pure layered convection with the upper and lower mantle convecting separately from each other, and pure, whole mantle convection, appear not to satisfy all the observations, several addition models have been proposed. These models include and attempt to characterize least one reservoir that is enriched in radiogenic elements relative to the mid-ocean ridge basalt (MORB) source, as is required to account for most current estimates of the Earth's heat budget. This reservoir would also be responsible for the geochemical signature in some ocean island basalts (OIBs) like Hawaii, but must be rarely sampled at the surface. Our current knowledge of the mass- and heat-budget for the bulk silicate Earth from geochemical, cosmochemical and geodynamical observations and constraints enables us to quantify the radiogenic heat enrichment required to balance the heat budget. Without assuming any particular model for the structure of the reservoir, we first determine the inherent trade-off between heat production rate and mass of the reservoir. Using these constraints, we then investigate the dynamical inferences of the heat budget, assuming that the additional heat is produced within a deep layer above the core-mantle boundary. We carry out dynamical models of layered convection using four different fixed reservoir volumes, corresponding to deep layers of thicknesses 150, 500 1000 and 1600 km, respectively, and including both temperature-dependent viscosity and an instrinsic viscosity jump between upper and lower mantle. We then assess the viability of these cases against 5 criteria: stability of the deep layer through time, topography of the interface, effective density profile, intrinsic chemical density and the heat flux at the CMB.

Mantle Flow in the Western United States Constrained by Seismic Anisotropy

NASA Astrophysics Data System (ADS)

Niday, W.; Humphreys, E.

2017-12-01

Shear wave splitting, caused by the lattice preferred orientation (LPO) of olivine crystals under shear deformation, provide a useful constraint on numerical models of mantle flow. Although it is sometimes assumed that shear wave splitting fast directions correspond with mantle flow directions, this is only true in simple shear flows that do not vary strongly with space or time. Observed shear wave splitting in the western United States is complex and inconsistent with simple shear driven by North American and Pacific plate motion, suggesting that the effects of time-dependent subduction history and spatial heterogeneity are important. Liu and Stegman (2011) reproduce the pattern of fast seismic anomalies below the western US from Farallon subduction history, and Chaparro and Stegman (2017) reproduce the circular anisotropy field below the Great Basin. We extend this to consider anisotropic structure outside the Great Basin and evaluate the density and viscosity of seismic anomalies such as slabs and Yellowstone. We use the mantle convection code ASPECT to simulate 3D buoyancy-driven flow in the mantle below the western US, and predict LPO using the modeled flow fields. We present results from a suite of models varying the sub-lithospheric structures of the western US and constraints on density and viscosity variations in the upper mantle.

The importance of grain size to mantle dynamics and seismological observations

NASA Astrophysics Data System (ADS)

Gassmoeller, R.; Dannberg, J.; Eilon, Z.; Faul, U.; Moulik, P.; Myhill, R.

2017-12-01

Grain size plays a key role in controlling the mechanical properties of the Earth's mantle, affecting both long-timescale flow patterns and anelasticity on the timescales of seismic wave propagation. However, dynamic models of Earth's convecting mantle usually implement flow laws with constant grain size, stress-independent viscosity, and a limited treatment of changes in mineral assemblage. We study grain size evolution, its interplay with stress and strain rate in the convecting mantle, and its influence on seismic velocities and attenuation. Our geodynamic models include the simultaneous and competing effects of dynamic recrystallization resulting from dislocation creep, grain growth in multiphase assemblages, and recrystallization at phase transitions. They show that grain size evolution drastically affects the dynamics of mantle convection and the rheology of the mantle, leading to lateral viscosity variations of six orders of magnitude due to grain size alone, and controlling the shape of upwellings and downwellings. Using laboratory-derived scaling relationships, we convert model output to seismologically-observable parameters (velocity, attenuation) facilitating comparison to Earth structure. Reproducing the fundamental features of the Earth's attenuation profile requires reduced activation volume and relaxed shear moduli in the lower mantle compared to the upper mantle, in agreement with geodynamic constraints. Faster lower mantle grain growth yields best fit to seismic observations, consistent with our re-examination of high pressure grain growth parameters. We also show that ignoring grain size in interpretations of seismic anomalies may underestimate the Earth's true temperature variations.

Location and extent of Tertiary structures in Cook Inlet Basin, Alaska, and mantle dynamics that focus deformation and subsidence

USGS Publications Warehouse

Haeussler, Peter J.; Saltus, Richard W.

2011-01-01

Subduction of the buoyant Yakutat microplate likely caused deformation to be focused preferentially in upper Cook Inlet. The upper Cook Inlet region has both the highest degree of shortening and the deepest part of the Neogene basin. This forearc region has a long-wavelength magnetic high, a large isostatic gravity low, high conductivity in the lower mantle, low p-wave velocity (Vp), and a high p-wave to shear-wave velocity ratio (Vp/Vs). These data suggest that fluids in the mantle wedge caused serpentinization of mafic rocks, which may, at least in part, contribute to the long-wavelength magnetic anomaly. This area lies adjacent to the subducting and buoyant Yakutat microplate slab. We suggest the buoyant Yakutat slab acts much like a squeegee to focus mantle-wedge fluid flow at the margins of the buoyant slab. Such lateral flow is consistent with observed shear-wave splitting directions. The additional fluid in the adjacent mantle wedge reduces the wedge viscosity and allows greater corner flow. This results in focused subsidence, deformation, and gravity anomalies in the forearc region.

An analytic model of axisymmetric mantle plume due to thermal and chemical diffusion

NASA Technical Reports Server (NTRS)

Liu, Mian; Chase, Clement G.

1990-01-01

An analytic model of axisymmetric mantle plumes driven by either thermal diffusion or combined diffusion of both heat and chemical species from a point source is presented. The governing equations are solved numerically in cylindrical coordinates for a Newtonian fluid with constant viscosity. Instead of starting from an assumed plume source, constraints on the source parameters, such as the depth of the source regions and the total heat input from the plume sources, are deduced using the geophysical characteristics of mantle plumes inferred from modelling of hotspot swells. The Hawaiian hotspot and the Bermuda hotspot are used as examples. Narrow mantle plumes are expected for likely mantle viscosities. The temperature anomaly and the size of thermal plumes underneath the lithosphere can be sensitive indicators of plume depth. The Hawaiian plume is likely to originate at a much greater depth than the Bermuda plume. One suggestive result puts the Hawaiian plume source at a depth near the core-mantle boundary and the source of the Bermuda plume in the upper mantle, close to the 700 km discontinuity. The total thermal energy input by the source region to the Hawaiian plume is about 5 x 10(10) watts. The corresponding diameter of the source region is about 100 to 150 km. Chemical diffusion from the same source does not affect the thermal structure of the plume.

Mixing in heterogeneous internally-heated convection

NASA Astrophysics Data System (ADS)

Limare, A.; Kaminski, E. C.; Jaupart, C. P.; Farnetani, C. G.; Fourel, L.; Froment, M.

2017-12-01

Past laboratory experiments of thermo chemical convection have dealt with systems involving fluids with different intrinsic densities and viscosities in a Rayleigh-Bénard setup. Although these experiments have greatly improved our understanding of the Earth's mantle dynamics, they neglect a fundamental component of planetary convection: internal heat sources. We have developed a microwave-based method in order to study convection and mixing in systems involving two layers of fluid with different densities, viscosities, and internal heat production rates. Our innovative laboratory experiments are appropriate for the early Earth, when the lowermost mantle was likely enriched in incompatible and heat producing elements and when the heat flux from the core probably accounted for a small fraction of the mantle heat budget. They are also relevant to the present-day mantle if one considers that radioactive decay and secular cooling contribute both to internal heating. Our goal is to quantify how two fluid layers mix, which is still very difficult to resolve accurately in 3-D numerical calculations. Viscosities and microwave absorptions are tuned to achieve high values of the Rayleigh-Roberts and Prandtl numbers relevant for planetary convection. We start from a stably stratified system where the lower layer has higher internal heat production and density than the upper layer. Due to mixing, the amount of enriched material gradually decreases to zero over a finite time called the lifetime. Based on more than 30 experiments, we have derived a scaling law that relates the lifetime of an enriched reservoir to the layer thickness ratio, a, to the density and viscosity contrasts between the two layers, and to their two different internal heating rates in the form of an enrichment factor beta=1+2*a*H1/H, where H1 is the heating rate of the lower fluid and H is the average heating rate. We find that the lifetime of the lower enriched reservoir varies as beta**(-7/3) in the low viscosity contrast limit, and as beta**(-4/3) in the large viscosity contrast limit. Our state-of-the-art experimental technique thus provides insights on chemical differentiation processes and on the evolution of mantle heterogeneities on both short and long time-scales.

Deep and persistent melt layer in the Archaean mantle

NASA Astrophysics Data System (ADS)

Andrault, Denis; Pesce, Giacomo; Manthilake, Geeth; Monteux, Julien; Bolfan-Casanova, Nathalie; Chantel, Julien; Novella, Davide; Guignot, Nicolas; King, Andrew; Itié, Jean-Paul; Hennet, Louis

2018-02-01

The transition from the Archaean to the Proterozoic eon ended a period of great instability at the Earth's surface. The origin of this transition could be a change in the dynamic regime of the Earth's interior. Here we use laboratory experiments to investigate the solidus of samples representative of the Archaean upper mantle. Our two complementary in situ measurements of the melting curve reveal a solidus that is 200-250 K lower than previously reported at depths higher than about 100 km. Such a lower solidus temperature makes partial melting today easier than previously thought, particularly in the presence of volatiles (H2O and CO2). A lower solidus could also account for the early high production of melts such as komatiites. For an Archaean mantle that was 200-300 K hotter than today, significant melting is expected at depths from 100-150 km to more than 400 km. Thus, a persistent layer of melt may have existed in the Archaean upper mantle. This shell of molten material may have progressively disappeared because of secular cooling of the mantle. Crystallization would have increased the upper mantle viscosity and could have enhanced mechanical coupling between the lithosphere and the asthenosphere. Such a change might explain the transition from surface dynamics dominated by a stagnant lid on the early Earth to modern-like plate tectonics with deep slab subduction.

Mantle viscosity structure constrained by joint inversions of seismic velocities and density

NASA Astrophysics Data System (ADS)

Rudolph, M. L.; Moulik, P.; Lekic, V.

2017-12-01

The viscosity structure of Earth's deep mantle affects the thermal evolution of Earth, the ascent of mantle upwellings, sinking of subducted oceanic lithosphere, and the mixing of compositional heterogeneities in the mantle. Modeling the long-wavelength dynamic geoid allows us to constrain the radial viscosity profile of the mantle. Typically, in inversions for the mantle viscosity structure, wavespeed variations are mapped into density variations using a constant- or depth-dependent scaling factor. Here, we use a newly developed joint model of anisotropic Vs, Vp, density and transition zone topographies to generate a suite of solutions for the mantle viscosity structure directly from the seismologically constrained density structure. The density structure used to drive our forward models includes contributions from both thermal and compositional variations, including important contributions from compositionally dense material in the Large Low Velocity Provinces at the base of the mantle. These compositional variations have been neglected in the forward models used in most previous inversions and have the potential to significantly affect large-scale flow and thus the inferred viscosity structure. We use a transdimensional, hierarchical, Bayesian approach to solve the inverse problem, and our solutions for viscosity structure include an increase in viscosity below the base of the transition zone, in the shallow lower mantle. Using geoid dynamic response functions and an analysis of the correlation between the observed geoid and mantle structure, we demonstrate the underlying reason for this inference. Finally, we present a new family of solutions in which the data uncertainty is accounted for using covariance matrices associated with the mantle structure models.

Differentiation of Secular and Postseismic Deformation in the Mojave Shear Zone in Southern California and Inference of Lithospheric Rheology

NASA Astrophysics Data System (ADS)

Shen, Z.; Liu, S.; Burgmann, R.

2015-12-01

The 1992 Mw 7.3 Landers and 1999 Mw7.1 Hector Mine earthquakes struck the Eastern California Shear Zone (ECSZ) in the Mojave Desert, Southern California. Coseismic and postseismic deformation from these events affect efforts to use Global Positioning System (GPS) observations collected since these events to establish a secular surface velocity field, especially in the near field of the coseismic ruptures. We devise block motion models constrained by both historical pre-Landers triangulation and trilateration observations and post-Landers GPS measurements to recover the secular deformation field and differentiate the postseismic transients in the Mojave region. Postseismic transients are found to remain in various "interseismic" GPS velocity solutions in the form of 2-3 mm/yr excess right-lateral shear across the Landers and Hector Mine coseismic ruptures [Liu et al., 2015 JGR]. Postseismic GPS time series differentiated from the secular velocity field reveal enduring late-stage transient motions in the near field of the coseismic ruptures. Using the postseismic time series data as model constraints, we develop postseismic deformation model invoking afterlip on faults and viscoelastic relaxation in the lower crust and upper mantle. A Burgers body material and a Maxwell material are assumed for the lower crust and upper mantle respectively. Our preliminary modeling result, constrained using GPS time series data from the SCEC Crustal Motion Map 4.0 (covering the time period of 1992-2004), reveals that both the long-term viscosities for the lower crust and upper mantle are on the order of e+19 Pa-s. This finding differs significantly from the "Crème Brulee" model predictions about the rheological structure of the lower crust and upper mantle, in which the lower crust has a substantially higher viscosity. We are incorporating more GPS time series data into our model, particularly the ones from continuous sites of the Plate Boundary Observatory network with post-2004 time span, and the modeling result will be presented at the meeting.

Thermomechanical earthquake cycle simulations with rate-and-state friction and nonlinear viscoelasticity

NASA Astrophysics Data System (ADS)

Allison, K. L.; Dunham, E. M.

2017-12-01

We simulate earthquake cycles on a 2D strike-slip fault, modeling both rate-and-state fault friction and an off-fault nonlinear power-law rheology. The power-law rheology involves an effective viscosity that is a function of temperature and stress, and therefore varies both spatially and temporally. All phases of the earthquake cycle are simulated, allowing the model to spontaneously generate earthquakes, and to capture frictional afterslip and postseismic and interseismic viscous flow. We investigate the interaction between fault slip and bulk viscous flow, using experimentally-based flow laws for quartz-diorite in the crust and olivine in the mantle, representative of the Mojave Desert region in Southern California. We first consider a suite of three linear geotherms which are constant in time, with dT/dz = 20, 25, and 30 K/km. Though the simulations produce very different deformation styles in the lower crust, ranging from significant interseismc fault creep to purely bulk viscous flow, they have almost identical earthquake recurrence interval, nucleation depth, and down-dip coseismic slip limit. This indicates that bulk viscous flow and interseismic fault creep load the brittle crust similarly. The simulations also predict unrealistically high stresses in the upper crust, resulting from the fact that the lower crust and upper mantle are relatively weak far from the fault, and from the relatively small role that basal tractions on the base of the crust play in the force balance of the lithosphere. We also find that for the warmest model, the effective viscosity varies by an order of magnitude in the interseismic period, whereas for the cooler models it remains roughly constant. Because the rheology is highly sensitive to changes in temperature, in addition to the simulations with constant temperature we also consider the effect of heat generation. We capture both frictional heat generation and off-fault viscous shear heating, allowing these in turn to alter the effective viscosity. The resulting temperature changes may reduce the width of the shear zone in the lower crust and upper mantle, and reduce the effective viscosity.

Modes of mantle convection and the removal of heat from the earth's interior

NASA Technical Reports Server (NTRS)

Spohn, T.; Schubert, G.

1982-01-01

Thermal histories for two-layer and whole-mantle convection models are calculated and presented, based on a parameterization of convective heat transport. The model is composed of two concentric spherical shells surrounding a spherical core. The models were constrained to yield the observed present-day surface heat flow and mantle viscosity, in order to determine parameters. These parameters were varied to determine their effects on the results. Studies show that whole-mantle convection removes three times more primordial heat from the earth interior and six times more from the core than does two-layer convection (in 4.5 billion years). Mantle volumetric heat generation rates for both models are comparable to that of a potassium-depleted chondrite, and thus surface heat-flux balance does not require potassium in the core. Whole and two-layer mantle convection differences are primarily due to lower mantle thermal insulation and the lower heat removal efficiency of the upper mantle as compared with that of the whole mantle.

Modulation of Intraplate Deformation in Arizona due to Far-Reaching Postseismic Relaxation Following Plate Boundary Earthquakes

NASA Astrophysics Data System (ADS)

Kreemer, C.; Broermann, J.; Bennett, R. A.; Blewitt, G.

2016-12-01

The Basin and Range of southern Arizona is mostly devoid of active faults and seismicity, although M7+ earthquakes are known to have occurred in the Quaternary. Before the 2010 Mw=7.2 El-Mayor Cucapah (EMC) earthquake, we found from GPS data 2 mm/yr of extension across this region. After the EMC, we find that the extension across southern Arizona has been reduced to near-zero; this is part of a widespread change in the region's post-EMC deformation field that can be explained by viscoelastic relaxation of the lower crust and upper mantle. If the relaxation after the EMC is so wide-spread, similar affects must have occurred after other previous large earthquakes such that the deformation field is probably ever-changing, slowly. This transient deformation field obscures any long-term deformation and makes it difficult to assess the driving forces responsible for the deformation. To model the transient deformation we consider that the observed GPS-time-series record postseismic relaxation from the most recent large earthquakes: 1992 Landers, 1999 Hector Mine, 2009 Gulf of California, 2010 EMC, and 2012 Gulf of California. The postseismic model assumes an elastic upper crust to a depth of 15 km, a Maxwell visco-elastic lower crust at 15-30 km, a Burger's body rheology for the lithospheric mantle at 30-60 km, and an asthenosphere below 60 km. A best fitting viscosity was found by a random search over the parameter space of the three viscous layers and minimizing the misfit between the corrected time series and a least squares model that includes intercept, velocity, offsets at times of earthquakes and equipment offsets and annual and semi-annual terms. The overall best fitting viscosity for the lower crust 1019.9 Pa-s, lithospheric mantle 1019.6 Pa-s, and asthenospheric mantle 1018.8 Pa-s. Both the mantle viscosities are treated as Burger rheologies with a Kelvin Voigt element 1/10 of the Maxwell element. These viscosities predict very small present-day velocity gradients associated with relaxation from pre-1992 events. We will use the corrected GPS velocity field to explore the driving sources behind the velocity field; i.e., plate motion and spatial variation in gravitational potential energy.

The fate of water within Earth and super-Earths and implications for plate tectonics

PubMed Central

2017-01-01

The Earth is likely to have acquired most of its water during accretion. Internal heat of planetesimals by short-lived radioisotopes would have caused some water loss, but impacts into planetesimals were insufficiently energetic to produce further drying. Water is thought to be critical for the development of plate tectonics, because it lowers viscosities in the asthenosphere, enabling subduction. The following issue persists: if water is necessary for plate tectonics, but subduction itself hydrates the upper mantle, how is the upper mantle initially hydrated? The giant impacts of late accretion created magma lakes and oceans, which degassed during solidification to produce a heavy atmosphere. However, some water would have remained in the mantle, trapped within crystallographic defects in nominally anhydrous minerals. In this paper, we present models demonstrating that processes associated with magma ocean solidification and overturn may segregate sufficient quantities of water within the upper mantle to induce partial melting and produce a damp asthenosphere, thereby facilitating plate tectonics and, in turn, the habitability of Earth-like extrasolar planets. This article is part of the themed issue ‘The origin, history and role of water in the evolution of the inner Solar System’. PMID:28416729

The fate of water within Earth and super-Earths and implications for plate tectonics.

PubMed

Tikoo, Sonia M; Elkins-Tanton, Linda T

2017-05-28

The Earth is likely to have acquired most of its water during accretion. Internal heat of planetesimals by short-lived radioisotopes would have caused some water loss, but impacts into planetesimals were insufficiently energetic to produce further drying. Water is thought to be critical for the development of plate tectonics, because it lowers viscosities in the asthenosphere, enabling subduction. The following issue persists: if water is necessary for plate tectonics, but subduction itself hydrates the upper mantle, how is the upper mantle initially hydrated? The giant impacts of late accretion created magma lakes and oceans, which degassed during solidification to produce a heavy atmosphere. However, some water would have remained in the mantle, trapped within crystallographic defects in nominally anhydrous minerals. In this paper, we present models demonstrating that processes associated with magma ocean solidification and overturn may segregate sufficient quantities of water within the upper mantle to induce partial melting and produce a damp asthenosphere, thereby facilitating plate tectonics and, in turn, the habitability of Earth-like extrasolar planets.This article is part of the themed issue 'The origin, history and role of water in the evolution of the inner Solar System'. © 2017 The Authors.

Glacial Isostatic Adjustment with ICE-6G{_}C (VM5a) and Laterally Heterogeneous Mantle Viscosity

NASA Astrophysics Data System (ADS)

Li, Tanghua; Wu, Patrick; Steffen, Holger

2017-04-01

Recently, Peltier et al. (2015) introduced the ICE-6GC (VM5a) ice-earth model pair, which has successfully explained many observations of Glacial Isostatic Adjustment (GIA) simultaneously. However, their earth model used (VM5a) to infer the ice history (ICE-6G_C) is laterally homogeneous with viscosity profile varying in the radial direction only. Since surface geology and seismic tomography clearly indicates that the Earth's material properties also vary in the lateral direction, laterally heterogeneity must be included in GIA models. This can be achieved by using the Coupled-Laplace-Finite-Element method (Wu 2004) to model GIA in a spherical, self-gravitating, compressible viscoelastic Earth with linear rheology and lateral heterogeneity. In fact, Wu et al (2013) have used such model with GIA observations (e.g. global relative sea level data, GRACE data with recent hydrology contributions removed and GPS crustal uplift rates) to study the thermal contribution to lateral heterogeneity in the mantle. Their lateral viscosity perturbations are inferred from the seismic shear wave tomography model S20A (Ekstrom & Dziewonski 1998) by applying a scaling law, which includes both the effect of anharmonicity and anelasticity. The thermal contribution to seismic tomography, which is represented by the beta factor in the scaling relationship, is searched in the upper and lower mantle, for the best combination that gives the best fit between GIA predictions and observations. However, their study is based on ICE-4G only, and the new ice-earth model pair may give other best beta value combinations in the upper and lower mantle. Here, we follow the work of Wu et al (2013) but use the new ICE-6GC ice model instead. The higher resolution seismic tomography model by Bunge & Grand (2000) substitutes S20A. Earth model VM5a is used as the reference background viscosity model. The full viscosity model is obtained by superposing the background model with the lateral viscosity perturbations inferred from the seismic tomography model (Bunge & Grand 2000) logarithmically. The preliminary results of these and other background viscosity profiles will be presented. References: Bunge, H.-P. & Grand, S. P. (2000). Mesozoic plate-motion history below the northeast Pacific Ocean from seismic images of the subducted Farallon slab. Nature, 405(6784):337-340. Peltier, W., Argus, D., and Drummond, R. (2015). Space geodesy constrains ice age terminal deglaciation: The global ICE-6GC (VM5a) model. Journal of Geophysical Research: Solid Earth, 120(1): 450-487. Wu, P. (2004). Using commercial finite element packages for the study of earth deformations, sea levels and the state of stress. Geophysical Journal International, 158(2): 401-408. Wu, P., Wang, H.S. & Steffen, H. (2012). The role of thermal effect on mantle seismic anomalies under Laurentia and Fennoscandia from observations of Glacial Isostatic Adjustment. Geophysical Journal International, 192(1):7-17.

Influence of the Structural Dichotomy of Antarctic Lithosphere on Regional Glacial-Isostatic Adjustment

NASA Astrophysics Data System (ADS)

Klemann, V.; Rau, D.; Martinec, Z.; Wolf, D.

2009-05-01

The strong structural dichotomy between East and West Antarctica is related to the West Antarctic Rift. The rheological implications are a reduction of the elastic-lithosphere thickness by a factor of more than 2 from East to West Antarctica as well as a strongly reduced mantle viscosity below West Antarctica and the Antarctic Peninsula. For modelling glacial-isostatic adjustment, we use a global viscoelastic earth model and apply the spectral finite-element method for the solution of the field equations. Ice models ICE-5G and IJ05 are used for parameterizing the last Pleistocene deglaciation. Lateral viscosity variations in the upper mantle are derived from variations in seismic velocity by applying scaling laws. Considering also lateral variations in the lithosphere structure, we study the implications of lateral variability on the glacial-isostatic adjustment of Antarctica.

Olivine water contents in the continental lithosphere and the longevity of cratons.

PubMed

Peslier, Anne H; Woodland, Alan B; Bell, David R; Lazarov, Marina

2010-09-02

Cratons, the ancient cores of continents, contain the oldest crust and mantle on the Earth (>2 Gyr old). They extend laterally for hundreds of kilometres, and are underlain to depths of 180-250 km by mantle roots that are chemically and physically distinct from the surrounding mantle. Forming the thickest lithosphere on our planet, they act as rigid keels isolated from the flowing asthenosphere; however, it has remained an open question how these large portions of the mantle can stay isolated for so long from mantle convection. Key physical properties thought to contribute to this longevity include chemical buoyancy due to high degrees of melt-depletion and the stiffness imparted by the low temperatures of a conductive thermal gradient. Geodynamic calculations, however, suggest that these characteristics are not sufficient to prevent the lithospheric mantle from being entrained during mantle convection over billions of years. Differences in water content are a potential source of additional viscosity contrast between cratonic roots and ambient mantle owing to the well-established hydrolytic weakening effect in olivine, the most abundant mineral of the upper mantle. However, the water contents of cratonic mantle roots have to date been poorly constrained. Here we show that olivine in peridotite xenoliths from the lithosphere-asthenosphere boundary region of the Kaapvaal craton mantle root are water-poor and provide sufficient viscosity contrast with underlying asthenosphere to satisfy the stability criteria required by geodynamic calculations. Our results provide a solution to a puzzling mystery of plate tectonics, namely why the oldest continents, in contrast to short-lived oceanic plates, have resisted recycling into the interior of our tectonically dynamic planet.

The effect of ilmenite viscosity on the dynamics and evolution of an overturned lunar cumulate mantle

NASA Astrophysics Data System (ADS)

Zhang, Nan; Dygert, Nick; Liang, Yan; Parmentier, E. M.

2017-07-01

Lunar cumulate mantle overturn and the subsequent upwelling of overturned mantle cumulates provide a potential framework for understanding the first-order thermochemical evolution of the Moon. Upwelling of ilmenite-bearing cumulates (IBCs) after the overturn has a dominant influence on the dynamics and long-term thermal evolution of the lunar mantle. An important parameter determining the stability and convective behavior of the IBC is its viscosity, which was recently constrained through rock deformation experiments. To examine the effect of IBC viscosity on the upwelling of overturned lunar cumulate mantle, here we conduct three-dimensional mantle convection models with an evolving core superposed by an IBC-rich layer, which resulted from mantle overturn after magma ocean solidification. Our modeling shows that a reduction of mantle viscosity by 1 order of magnitude, due to the presence of ilmenite, can dramatically change convective planform and long-term lunar mantle evolution. Our model results suggest a relatively stable partially molten IBC layer that has surrounded the lunar core to the present day.