Viscosity in transition zone and lower mantle: Implications for slab penetration
|
|
- Melvin Harris
- 6 years ago
- Views:
Transcription
1 Click Here for Full Article GEOPHYSICAL RESEARCH LETTERS, VOL. 37,, doi: /2010gl043140, 2010 Viscosity in transition zone and lower mantle: Implications for slab penetration J. Quinteros, 1,2 S. V. Sobolev, 1,3 and A. A. Popov 1 Received 4 March 2010; revised 3 April 2010; accepted 13 April 2010; published 13 May [1] The evolution of a subducting slab is strongly influenced by the viscosity of the mantle which it overlies. However, there is still no consensus about the viscosity in the transition zone and shallower lower mantle. We use a numerical self consistent subduction model and run a set of experiments in order to find critical patterns of viscosity that would allow the evolution of the different subduction styles that can actually be found in nature. Our results show that a jump in viscosity of 5 times from transition zone to lower mantle gives the most reasonable results. Optimal values of viscosity in the transition zone are in the range of 3 * Pa s. Higher values would produce piling up of the slab and later sinking or even slab flattening. Lower viscosities result in velocities (>30 cm/yr) that are too high, while the new slab subducts through the upper mantle and transition zone, a phenomenon that is rarely seen in nature. Reduction of the Clapeyron slope, related to the spinel perovskite transition, variations in oceanic crustal thickness and in the age of the slab do not influence much the style of subduction. When overriding velocity is applied to the upper plate, the previously penetrating slab tends to lay at the 660 km boundary but does not substantially change the subduction velocity. Citation: Quinteros, J., S. V. Sobolev, and A. A. Popov (2010), Viscosity in transition zone and lower mantle: Implications for slab penetration, Geophys. Res. Lett., 37,, doi: /2010gl Introduction [2] It is accepted that over a long timescale, mantle behaves like a viscous fluid, while lithosphere can be seen as a brittle layer. However, a viscosity profile throughout the whole mantle is far from being accurate. Different models were proposed to find the best fit viscosity distribution based on the imposition of different sets of constraints (e.g., geoid, Haskell, tomography model) [King, 1995; Steinberger and Calderwood, 2006]. There are some common factors in all the possible profiles: 1) there seems to be a sharp viscosity increment between transition zone (TZ) and lower mantle (LM) and 2) a slow increase in viscosity with depth in the first hundreds of km of LM. [3] Many authors [e.g., Mitrovica and Forte, 2004] agree that the asthenosphere should be the layer with lowest viscosity. However, others claim that the TZ maybe weaker 1 Section 2.5, Deutsches GeoForschungsZentrum, Potsdam, Germany. 2 Department of Computer Sciences, FCEN, UBA, Buenos Aires, Argentina. 3 Institute of Physics of the Earth, Moscow, Russia. Copyright 2010 by the American Geophysical Union /10/2010GL than upper mantle (UM) and LM due to the presence of water stored in wadsleyite or ringwoodite. [4] We considered a range of possible viscosity values for the TZ presented by Steinberger and Calderwood [2006] and tested their feasibility by numerical simulation of a subducting slab. As viscosity values should be compatible with the different types of subduction detected by seismic tomography, in this work we present a set of experiments designed to study and constrain the possible critical patterns of viscosity that could determine the fate of the slab. [5] In the last decades stagnant slabs have been identified in different subduction zones (for a review see Fukao et al. [2009]). The positive Clapeyron slope associated with the olivine spinel phase transition at 410 km favours that the cold slab can go through it, accelerating the subduction process. The opposite situation seems to occur on the spinelperovskite phase transition at 660 km, where the Clapeyron slope is negative and inhibits or slows the subduction [Tackley et al., 1993]. Whether the slab is able to sink beyond this boundary depends, not only on the value of the Clapeyron slope, but also on other factors like absolute velocities of the plates and viscosity of the underlying mantle. [6] The increment on viscosity related to the 660 km boundary is likely present [e.g., Hager, 1984], but its magnitude remains unclear. In numerical experiments, this varies from 0 to 50 times [Čížková et al., 2002; Pysklywec and Ishii, 2005; Kaus et al., 2008; Enns et al., 2005]. 2. Previous Studies [7] Christensen [1996] studied the relationship of trench migration and slab penetration into lower mantle. He was the first who described different types of subduction styles as a function of the kinematically imposed trench retreat. One of the main drawbacks of the kinematic conditions applied is the lack of feedback between the subducted slab and the velocity of subduction that was held constant throughout all the simulations. [8] Čížková et al. [2002] performed a parametric study on rheological weakening in slabs that interact with the 660 km discontinuity. They varied the yield stress (in order to approximate Peierls creep mechanism) and the grain size. They conclude that stress limit is more important than grain size reduction in the evolution of the slab, because it controls the deformation on the outer rigid part of the slab, while grain size only decreases viscosity in the inner part [Karato et al., 2001]. They also claim that further studies should be made on the influence of stepwise viscosity increase at 660 km, among other factors. [9] Enns et al. [2005] investigated the influence of viscosity stratification on the trench migration. They showed that phase transition dynamics and density or viscosity 1of5
2 Figure 1. (a) Numerical model setup and boundary conditions. Subducting slab is 40 My old and is pushed on the left side for 2 My at 15 cm/yr to initiate subduction. (b) Viscosity plot and temperature contour in UM at the beginning of simulation. layering can produce a variety of ponded slab shapes near the 660 km boundary. However, the effects of the Clapeyron slopes for both phase transitions and of the increment of density related to the 410 km boundary were not taken into account. [10] Mishin et al. [2008] studied the penetration of the 660 km boundary by slabs in double subduction systems. They found that penetration is strongly influenced by the convergence rate and the relative movement of the overriding plate, among other factors. 3. Numerical Method and Set Up [11] We used a coupled elasto visco plastic thermomechanical numerical model based on code SLIM 3D [Popov and Sobolev, 2008] to run all experiments. As our set of experiments is inherently 2D, the setup contains only one element on one of the coordinates. The model has true free surface and elastic deformation is included. Realistic rheology was applied to the subducting and overriding plates as well as the upper mantle. In this part of the model, viscosity is stress and temperature dependent. Three different types of creep (diffusion, dislocation and Peierls) are included, following the approach of Kameyama et al. [1999]. In particular, Peierls creep plays a key role to reduce stress in the deep cold slab. A complete list of the parameters used can be found in the auxiliary material. 1 [12] Note that in this model we do not include adiabatic compression. Therefore, temperatures in the model are potential temperatures, i.e., temperatures extrapolated to zero pressure along adiabat. [13] Volumetric deformations are followed by corrections in the calculation of density by means of the Murnaghan approach. Density is increased by 6% at the olivine spinel phase transition [Christensen, 1995] and 8% at the spinelperovskite [Dziewonski and Anderson, 1981]. [14] All phase transition boundaries (gabbro eclogite, coesite stishovite, olivine spinel and spinel perovskite) in this 1 Auxiliary materials are available in the HTML. doi: / 2009GL implementation are dynamic. This means that their position is the result of application of the P T diagram at every timestep for every particle in the domain. Clapeyron slope for olivine spinel transition is considered to be +2.0 MPa/K and 2.5 MPa/K for the spinel perovskite. These values coincide with seismological data [Fukao et al., 2009] and laboratory experiments by Litasov et al. [2005] and Irifune et al. [1998]. [15] Subduction channel was considered as a well lubricated interface with low viscosity and friction coefficient, which favours the development of realistic one sided subduction [Sobolev and Babeyko, 2005; Gerya et al., 2008]. This interface is three elements wide, defined at every time step and restricted to the contact between oceanic and continental material at a depth less than 120 km, treated similarly by Sobolev and Babeyko [2005]. [16] The structure of the overriding plate was also handled as was done by Sobolev and Babeyko [2005]. It consists of a continental crust (40 km) and lithospheric mantle up to a depth of 100 to 140 km. [17] Subducting slab is composed of 7 km of oceanic gabbro and 73 km of oceanic lithospheric mantle. It is considered to have a thermal age of 40 My, but according to Christensen [1996], slab age does not play a key role when determining the style of subduction. Trench motion is allowed and in none of our experiments was imposed neither kinematically nor by another means. [18] No overriding is prescribed in most of our experiments in order to isolate the influence produced by the tested values. However, some experiments where repeated with overriding velocity to estimate the influence on results. [19] The viscosity profile of TZ and LM was represented in a simplified way by means of two parameters: viscosity in TZ and viscosity step at the spinel perovskite phase transition. [20] We focus not only on the more realistic range of 5 to 30 times higher viscosity in LM, according to the results of Steinberger and Calderwood [2006], but also include simulations without a jump to evaluate the differences. The viscosity for TZ considered ranges from to 3 * Pa s as suggested by Steinberger and Calderwood [2006]. 2of5
3 Figure 2. Different styles of subduction detected in the experiments. Viscosity is shown in colours and temperature with contours. A vertical viscosity profile taken at x = 300 km is shown inside every plot. (a) Penetrating slab; (b) penetrating slab with some buckling; (c) piling up and slow sinking after flat slab with late penetration under subduction zone; (d) flat slab. [21] The domain is 1700 km wide and 1300 km deep. Slab is pushed with a constant velocity of 15 cm/yr at the beginning of simulation at its left edge for 2 My to develop the subduction zone and obtain a coherent thermal state. After that moment evolution is driven by gravity. A schematic view of the setup can be seen in Figure 1. [22] Left side of the domain is opened as only hydrostatic pressure boundary conditions were defined. This gives the slab the freedom to move on any direction and with any velocity and assures that subduction will evolve in a selfconsistent way. The top side of the domain is a free surface with zero stress. On the bottom and right side free slip conditions are applied. 4. Results [23] More than 40 experiments have been run and classified in order to find constraints on the viscosity pattern. Examples of the different styles of subduction that were classified based on the experiments are shown in Figure 2, from direct penetration (Figure 2a) to slabs that lay over the 660 km boundary (Figure 2d). [24] Taking Figure 2a as a reference case, one can see in Figure 2b that the increment of viscosity in TZ reduces subduction velocity by about 2 3 times (see Figure 3) and how penetration of a small part of the subducted slab is reaccommodated by slight buckling. [25] A very similar situation (not shown) is obtained if the viscosity step between TZ and LM is increased. Even if this is not enough to stop the subduction, it is actually retarded in LM. If the jump is further increased (Figure 2c) it is more difficult for the head of the slab to sink and this causes folding and piling up of the slab, which slowly reactivates the sinking process. [26] In the last case (Figure 2d) TZ and LM are more viscous and thus offer more resistance to slab penetration. Subduction velocity diminishes to 1 3 cm/yr during almost the whole simulation ( 50 My). At My slab starts folding and velocity increases slightly. Slab flattens even before being close to LM, what delays the contact until 3of5
4 Figure 3. Subduction velocity corresponding to experiments without overriding separated by viscosity in TM; (top) 3 * 10 20, (middle) 10 21, and (bottom) 3 * Pa s. Bold arrows indicate the moment when slab contacts or crosses the 660 km phase transition boundary. Codes in legends indicate the viscosity step. Subduction style is specified next to the lines. See caption in Table 1 for description. 23 My. As subduction continues, many mln. years later a pile is formed and slowly sinks. [27] In all the experiments, subduction accelerates when slab enters TZ. The main cause is the positive Clapeyron slope related to the olivine spinel phase transition, which further increases the density contrast caused by the low temperature. Velocity is reduced when the head of the slab gets closer to the 660 km boundary and after that it becomes stable for the next mln. years (see Figure 3). Depending on the case, one possibility is that subduction can go further but, if subduction velocity is faster than sinking velocity, folding of the slab can occur. In this case, subduction reactivates after folding and velocity increases causing piling up of the slab and slowly sinking of the pile. During this process of piling up and sinking, neither avalanche type nor drastic changes in the lower mantle flux were observed. [28] A summary of the cases without overriding, classified by viscosity in TZ and the jump between TZ and LM, can be seen in Table 1. These were categorized using a similar classification as suggested by Christensen [1996]. [29] It is important to mention that in comparison with Christensen [1996] much less buckling was found. The main cause seems to be the absence of kinematical boundary conditions for the slab. If the subduction velocity were kinematically imposed, slab would not be able to reduce its velocity when hitting the 660 km boundary and would generate spurious buckling. [30] Interestingly, in almost all the experiments, viscosity in TZ was higher than the one in the bottom of UM (Figure 2). [31] Additional experiments to test the influence of other variables show that substantial reduction of Clapeyron slope at 660 km has almost no influence on the results. We also ran experiments with an oceanic crust of 10 and 12 km and Table 1. Summary of Models Without Overriding a Viscosity in TM (in Pa.s) m lm /m tz * * P P P (B) 5 P P (B) RP PIL SP 10 P RP (B) FS PIL SP 20 P RP (B) RP PIL SP F PIL SP 30 P RP PIL SP FS PIL SP F PIL SP a P, straight penetration of slab; RP, retarded penetration; B, buckling; PIL, slab folding and piling up; SP, pile sinking into LM; FS, flat slab, then penetration near subduction zone; F, flat slab. 4of5
5 the subduction styles were not changed, while velocity patterns suffered only minor changes. Inclusion of overriding velocity influenced on the subduction style (more stagnant slabs) but not on the velocity pattern. 5. Conclusions [32] From our results we can conclude that the most reasonable range of viscosity values for the TZ would be (3 * ) Pa s. A value of 3 * is too high to allow a subducting slab to penetrate into the LM with a reasonable velocity if a viscosity step is present. The same situation arises if the viscosity step is greater than 10, in which case a piling up of material next to the LM or flat slabs can be seen, but no direct penetration. It should also be noted that, for higher viscosity jumps, viscosity in LM would be too large under certain combinations. A substantial reduction of the Clapeyron slope at the 660 km boundary did not help to enhance subduction and showed only a small influence over the evolution of the slab. [33] On the other hand, when viscosity in TZ is significantly lower than 3 * Pa s slab accelerates to more than 30 cm/yr while (or before) penetrating the TZ, a phenomenon that is rarely seen in nature. If viscosity step is not present, exaggerated or very high velocity values are also obtained for the most reasonable viscosity. The best results, considering the obtained subduction velocities, were obtained with a step around 5. Even if overriding velocity can change the style of subduction toward more flattened slab, it cannot reduce the exaggerated velocities obtained with the low values of both variables. [34] With a viscosity of 3 * Pa s in TZ, a step between 5 and 10 gives the best velocity patterns, considering that this is a case of fast subduction (15 cm/yr at the beginning). If viscosity is Pa s, the step between TZ and LM should be smaller than 5 in order to get reasonable subduction velocities. These viscosity values are in agreement with those suggested by Billen [2008]. [35] Acknowledgments. The authors want to thank Bernhard Steinberger, Benjamin Heit, Xiaoi Yuan and Reiner Kind for the fruitful discussions, as well as Taras Gerya and an anonymous reviewer for their constructive reviews, which helped us to enhance the quality of the paper. References Billen, M. I. (2008), Modeling the dynamics of subducting slabs, Annu. Rev. Earth Planet. Sci., 36(1), , doi: /annurev. earth Christensen, U. R. (1995), Effects of phase transitions on mantle convection, Annu. Rev. Earth Planet. Sci., 23(1), 65 87, doi: /annurev. ea Christensen, U. R. (1996), The influence of trench migration on slab penetration into the lower mantle, Earth Planet. Sci. Lett., 140, Čížková, H., J. van Hunen, A. P. van den Berg, and N. J. Vlaar (2002), The influence of rheological weakening and yield stress on the interaction of slabs with the 670 km discontinuity, Earth Planet. Sci. Lett., 199(3 4), , doi: /s x(02) Dziewonski, A. M., and D. L. Anderson (1981), Preliminary reference Earth model, Phys. Earth Planet. Inter., 25(4), , doi: / (81) Enns, A., T. W. Becker, and H. Schmeling (2005), The dynamics of subduction and trench migration for viscosity stratication, Geophys. J. Int., 160, , doi: /j x x. Fukao, Y., M. Obayashi, and T. Nakakuki (2009), Stagnant slab: A review, Annu. Rev. Earth Planet. Sci., 37, 19 46, doi: /annurev. earth Gerya, T. V., J. A. D. Connolly, and D. A. Yuen (2008), Why is terrestrial subduction one sided?, Geology, 36(1), 43 46, doi: /g24060a.1. Hager, B. H. (1984), Subducted slabs and the geoid: Constraints on mantle rheology and flow, J. Geophys. Res., 89, , doi: / JB089iB07p Irifune, T., N. Nishiyama, K. Kuroda, T. Inoue, M. Isshiki, W. Utsumi, K. ichi Funakoshi, S. Urakawa, T. Uchida, T. Katsura, and O. Ohtaka (1998), The Postspinel Phase Boundary in Mg2SiO4 Determined by in Situ X ray Diffraction, Science, 279(5357), , doi: / science Kameyama, M., D. A. Yuen, and S. Karato (1999), Thermal mechanical effects on low temperature plasticity (the Peierls mechanism) on the deformation of a viscoelastic shear zone, Earth Planet. Sci. Lett., 168, Karato, S., M. R. Riedel, and D. A. Yuen (2001), Rheological structure and deformation of subducted slabs in the mantle transition zone: Implications for mantle circulation and deep earthquakes, Phys. Earth Planet. Inter., 127(1 4), , doi: /s (01) Kaus, B. J., C. Steedman, and T. W. Becker (2008), From passive continental margin to mountain belt: Insights from analytical and numerical models and application to taiwan, Physics Earth Planet. Inter., 171(1 4), , doi: /j.pepi King, S. D. (1995), Radial models of mantle viscosity: Results from a genetic algorithm, Geophys. J. Int., 122(3), Litasov, K. D., E. Ohtani, A. Sano, A. Suzuki, and K. Funakoshi (2005), Wet subduction versus cold subduction, Geophys. Res. Lett., 32, L13312, doi: /2005gl Mishin, Y. A., T. V. Gerya, J. P. Burg, and J. A. Connolly (2008), Dynamics of double subduction: Numerical modeling, Phys. Earth Planet. Inter., 171(1 4), , doi: /j.pepi Mitrovica, J. X., and A. M. Forte (2004), A new inference of mantle viscosity based upon joint inversion of convection and glacial isostatic adjustment data, Earth Planet. Sci. Lett., 225, , doi: /j.epsl Popov, A. A., and S. V. Sobolev (2008), SLIM3D: A tool for threedimensional thermomechanical modeling of the lithospheric deformation with elasto visco plastic rheology, Phys. Earth Planet. Inter., 171, 55 75, doi: /j.pepi Pysklywec, R. N., and M. Ishii (2005), Time dependent subduction dynamics driven by the instability of stagnant slabs in the transition zone, Phys. Earth Planet. Inter., 149(1 2), , doi: /j. pepi Sobolev, S. V., and A. Y. Babeyko (2005), What drives orogeny in the Andes?, Geology, 33, Steinberger, B., and A. Calderwood (2006), Models of large scale viscous flow in the Earth s mantle with constraints from mineral physics and surface observations, Geophys. J. Int., 167, , doi: / j x x. Tackley, P. J., D. J. Stevenson, G. A. Glatzmaier, and G. Schubert (1993), Effects of an endothermic phase transition at 670 km depth in a spherical model of convection in the Earth s mantle, Nature, 361, , doi: /361699a0. A. A. Popov, J. Quinteros, and S. V. Sobolev, Section 2.5, Deutsches GeoForschungsZentrum, Telegrafenberg, D Potsdam, Germany. (anton@gfz potsdam.de; javier@gfz potsdam.de; stephan@gfz potsdam.de) 5of5
Why does the Nazca plate slow down since the Neogene? Supplemental Information
GSA DATA REPOSITORY 2013004 Why does the Nazca plate slow down since the Neogene? Supplemental Information Javier Quinteros a,b,, Stephan V. Sobolev a,c a Deutsches GeoForschungsZentrum GFZ, Telegrafenberg,
More informationNumerical Simulation of the Thermal Convection and Subduction Process in the Mantle
Chapter 1 Earth Science Numerical Simulation of the Thermal Convection and Subduction Process in the Mantle Project Representative Yoshio Fukao Institute for Research on Earth Evolution, Japan Agency for
More informationWhy cold slabs stagnate in the transition zone
GSA Data Repository 2015085 Model 1 Model 2 Model 3 Model 4 Model 5 Model 6 Why cold slabs stagnate in the transition zone Scott D. King 1,2, Daniel J. Frost 2, and David C. Rubie 2 1 Department of Geosciences,
More informationTime-varying subduction and rollback velocities in slab stagnation and buckling
Time-varying subduction and rollback velocities in slab stagnation and buckling Hana Čížková Charles University in Prague Craig Bina Northwestern University Evanston SLAB STAGNATION Fukao et al., 29 SLAB
More informationThe importance of the South-American plate motion and the Nazca Ridge subduction on flat subduction below South Peru
Chapter 7 The importance of the South-American plate motion and the Nazca Ridge subduction on flat subduction below South Peru Abstract Flat subduction near Peru occurs only where the thickened crust of
More informationGeodynamics Lecture 10 The forces driving plate tectonics
Geodynamics Lecture 10 The forces driving plate tectonics Lecturer: David Whipp! david.whipp@helsinki.fi!! 2.10.2014 Geodynamics www.helsinki.fi/yliopisto 1 Goals of this lecture Describe how thermal convection
More informationRheology of the Mantle and Plates (part 1): Deformation mechanisms and flow rules of mantle minerals
(part 1): Deformation mechanisms and flow rules of mantle minerals What is rheology? Rheology is the physical property that characterizes deformation behavior of a material (solid, fluid, etc) solid mechanics
More informationConstraints on Mantle Structure from Surface Observables
MYRES I: Heat, Helium & Whole Mantle Convection Constraints on Mantle Structure from Surface Observables Magali Billen University of California, Davis Department of Geology The Goal Use observations of
More informationSeismotectonics of intraplate oceanic regions. Thermal model Strength envelopes Plate forces Seismicity distributions
Seismotectonics of intraplate oceanic regions Thermal model Strength envelopes Plate forces Seismicity distributions Cooling of oceanic lithosphere also increases rock strength and seismic velocity. Thus
More informationThe influence of short wavelength variations in viscosity on subduction dynamics
1 Introduction Deformation within the earth, driven by mantle convection due primarily to cooling and subduction of oceanic lithosphere, is expressed at every length scale in various geophysical observations.
More informationModeling the Dynamics of Subducting Slabs
ANNUAL REVIEWS Further Click here for quick links to Annual Reviews content online, including: Other articles in this volume Top cited articles Top downloaded articles Our comprehensive search Annu. Rev.
More informationGravity Tectonics Volcanism Atmosphere Water Winds Chemistry. Planetary Surfaces
Gravity Tectonics Volcanism Atmosphere Water Winds Chemistry Planetary Surfaces Gravity & Rotation Polar flattening caused by rotation is the largest deviation from a sphere for a planet sized object (as
More informationThe in uence of rheological weakening and yield stress on the interaction of slabs with the 670 km discontinuity
Earth and Planetary Science Letters 199 (2002) 447^457 www.elsevier.com/locate/epsl The in uence of rheological weakening and yield stress on the interaction of slabs with the 670 km discontinuity Hana
More informationPlate tectonics. Chapter Slab pull. Dynamics of the Mantle and Lithosphere
Chapter 6 In a number of aspects the Earth s mantle can be regarded as a fluid that shows simple fluid-dynamical behavior. or instance, the mantle of the Earth is probably showing a form of convection.
More informationModeling the interior dynamics of terrestrial planets
Modeling the interior dynamics of terrestrial planets Paul J. Tackley, ETH Zürich Fabio Crameri, Tobias Keller, Marina Armann, Hein van Heck, Tobias Rolf Talk Plan Introduction Tectonic modes: plates,
More informationNumerical Simulation of the Mantle Convection and Subduction Process
Numerical Simulation of the Mantle Convection and Subduction Process Project Representative Yoshio Fukao Institute for Research on Earth Evolution, Japan Agency for Marine-Earth Science and Technology
More informationThe influence of phase boundary deflection on velocity anomalies of stagnant slabs in the transition zone
GEOPHYSICAL RESEARCH LETTERS, VOL. 30, NO. 18, 1965, doi:10.1029/2003gl017754, 2003 The influence of phase boundary deflection on velocity anomalies of stagnant slabs in the transition zone K. Chambers
More informationSubduction II Fundamentals of Mantle Dynamics
Subduction II Fundamentals of Mantle Dynamics Thorsten W Becker University of Southern California Short course at Universita di Roma TRE April 18 20, 2011 Rheology Elasticity vs. viscous deformation η
More informationThermal-Mechanical Behavior of Oceanic Transform Faults
Presented at the COMSOL Conference 2008 Boston Thermal-Mechanical Behavior of Oceanic Transform Faults COMSOL Conference - Boston, Massachusetts October 2008 Emily C. Roland - MIT/WHOI Joint Program Mark
More informationThe Earth s Structure from Travel Times
from Travel Times Spherically symmetric structure: PREM - Crustal Structure - Upper Mantle structure Phase transitions Anisotropy - Lower Mantle Structure D D - Structure of of the Outer and Inner Core
More informationLecture 2: Deformation in the crust and the mantle. Read KK&V chapter 2.10
Lecture 2: Deformation in the crust and the mantle Read KK&V chapter 2.10 Tectonic plates What are the structure and composi1on of tectonic plates? Crust, mantle, and lithosphere Crust relatively light
More informationC3.4.1 Vertical (radial) variations in mantle structure
C3.4 Mantle structure Mantle behaves as a solid on short time scales (seismic waves travel through it and this requires elastic behaviour). Over geological time scales the mantle behaves as a very viscous
More informationChapter 12 Lecture. Earth: An Introduction to Physical Geology. Eleventh Edition. Earth s Interior. Tarbuck and Lutgens Pearson Education, Inc.
Chapter 12 Lecture Earth: An Introduction to Physical Geology Eleventh Edition Earth s Interior Tarbuck and Lutgens Earth s Internal Structure Earth s interior can be divided into three major layers defined
More informationGlobal Tectonics. Kearey, Philip. Table of Contents ISBN-13: Historical perspective. 2. The interior of the Earth.
Global Tectonics Kearey, Philip ISBN-13: 9781405107778 Table of Contents Preface. Acknowledgments. 1. Historical perspective. 1.1 Continental drift. 1.2 Sea floor spreading and the birth of plate tectonics.
More informationPetrology: Igneous and Metamorphic (with a few sedimentary rocks)
1 Petrology: Igneous and Metamorphic (with a few sedimentary rocks) Spring 2017 Geology 2 nd ed. - Chernicoff Why Petrology? Most of the Earth is made of rock. Information about or for: Earth history rocks
More informationTopography of the 660-km discontinuity beneath northeast China: Implications for a retrograde motion of the subducting Pacific slab
GEOPHYSICAL RESEARCH LETTERS, VOL. 35, L01302, doi:10.1029/2007gl031658, 2008 Topography of the 660-km discontinuity beneath northeast China: Implications for a retrograde motion of the subducting Pacific
More informationPetrology. Petrology: the study of rocks, especially aspects such as physical, chemical, spatial and chronoligic. Associated fields include:
Petrology Petrology: the study of rocks, especially aspects such as physical, chemical, spatial and chronoligic. Associated fields include: Petrography: study of description and classification of rocks
More informationSurface changes caused by erosion and sedimentation were treated by solving: (2)
GSA DATA REPOSITORY 214279 GUY SIMPSON Model with dynamic faulting and surface processes The model used for the simulations reported in Figures 1-3 of the main text is based on two dimensional (plane strain)
More informationInfluence of surrounding plates on 3D subduction dynamics
Click Here for Full Article GEOPHYSICAL RESEARCH LETTERS, VOL. 36, L07303, doi:10.1029/2008gl036942, 2009 Influence of surrounding plates on 3D subduction dynamics P. Yamato, 1 L. Husson, 1 J. Braun, 1
More informationJournal of Applied Mathematics and Computation (JAMC), 2018, 2(7),
Journal of Applied Mathematics and Computation (JAMC), 2018, 2(7), 266-270 http://www.hillpublisher.org/journal/jamc ISSN Online:2576-0645 ISSN Print:2576-0653 The Solution Of 2D Hydrodynamic Equations
More informationDETAILS ABOUT THE TECHNIQUE. We use a global mantle convection model (Bunge et al., 1997) in conjunction with a
DETAILS ABOUT THE TECHNIQUE We use a global mantle convection model (Bunge et al., 1997) in conjunction with a global model of the lithosphere (Kong and Bird, 1995) to compute plate motions consistent
More informationPlate Tectonics. Structure of the Earth
Plate Tectonics Structure of the Earth The Earth can be considered as being made up of a series of concentric spheres, each made up of materials that differ in terms of composition and mechanical properties.
More informationPetrology. Petrology: the study of rocks, especially aspects such as physical, chemical, spatial and chronoligic. Classification:
Petrology Petrology: the study of rocks, especially aspects such as physical, chemical, spatial and chronoligic. Associated fields include: Petrography: study of description and classification of rocks
More informationPost-perovskite 1. Galley Proofs
16 September 2005 21:39 YB061180.tex McGraw Hill YB of Science & Technology Keystroked: 29/08/2005 Initial MS Page Sequence Stamp: 02350 Article Title: Post-perovskite Article ID: YB061180 1st Classification
More informationChapter 7 Plate Tectonics
Chapter 7 Plate Tectonics Earthquakes Earthquake = vibration of the Earth produced by the rapid release of energy. Seismic Waves Focus = the place within the Earth where the rock breaks, producing an earthquake.
More informationSlab pull, slab weakening, and their relation to deep intra-slab seismicity
GEOPHYSICAL RESEARCH LETTERS, VOL. 32, L14305, doi:10.1029/2005gl022922, 2005 Slab pull, slab weakening, and their relation to deep intra-slab seismicity Susan L. Bilek Earth and Environmental Science
More informationDurham Research Online
Durham Research Online Deposited in DRO: 22 January 2015 Version of attached le: Published Version Peer-review status of attached le: Peer-reviewed Citation for published item: Agrusta, R. and van Hunen,
More informationInfluence of the Ringwoodite-Perovskite transition on mantle convection in spherical geometry as a function of Clapeyron slope and Rayleigh number
Solid Earth, 2, 315 326, 2011 doi:10.5194/se-2-315-2011 Author(s) 2011. CC Attribution 3.0 License. Solid Earth Influence of the Ringwoodite-Perovskite transition on mantle convection in spherical geometry
More informationModeling the Thermal-Mechanical Behavior of Mid-Ocean Ridge Transform Faults
Excerpt from the Proceedings of the COMSOL Conference 2008 Boston Modeling the Thermal-Mechanical Behavior of Mid-Ocean Ridge Transform Faults Emily C Roland *1, Mark Behn,2 and Greg Hirth 3 1 MIT/WHOI
More informationThree-dimensional numerical simulations of thermo-chemical multiphase convection in Earth s mantle Takashi Nakagawa a, Paul J.
Three-dimensional numerical simulations of thermo-chemical multiphase convection in Earth s mantle Takashi Nakagawa a, Paul J. Tackley b a Department of Earth and Planetary Sciences, University of Tokyo,
More informationBeall et al., 2018, Formation of cratonic lithosphere during the initiation of plate tectonics: Geology, https://doi.org/ /g
GSA Data Repository 2018160 Beall et al., 2018, Formation of cratonic lithosphere during the initiation of plate tectonics: Geology, https://doi.org/10.1130/g39943.1. Numerical Modeling Methodology Stokes
More informationSupplementary information on the West African margin
Huismans and Beaumont 1 Data repository Supplementary information on the West African margin Interpreted seismic cross-sections of the north Angolan to south Gabon west African passive margins 1-3, including
More information4-D Geodynamic Modeling With Data Assimilation: Subduction and Continental Evolution
4-D Geodynamic Modeling With Data Assimilation: Subduction and Continental Evolution PI: Lijun Liu Department of Geology, University of Illinois at Urbana-Champaign Corresponding author: Lijun Liu, ljliu@illinois.edu
More informationThe Structure of the Earth and Plate Tectonics
The Structure of the Earth and Plate Tectonics Agree or Disagree? 1. The Earth if made up of 4 different layers. 2. The crust (where we live) can be made of either less dense continental crust or the more
More informationGRAVITY AND ISOSTASY
GRAVITY AND ISOSTASY Gravity The Geoid is the oblate spheroid -- the sea level surface over the entire Earth's surface Physically, the Geoid is an equipotential surface for gravity (i.e. the surface for
More informationSUPPLEMENTARY INFORMATION
SUPPLEMENTARY INFORMATION Supplementary information: Our 1080 km x 360 km model setup includes a 7 km thick oceanic crust adjacent to a 60 km thick, 250 km wide orogenic crust. Both the oceanic and the
More informationModification of the lithospheric stress field by lateral variations in plate-mantle coupling
GEOPHYSICAL RESEARCH LETTERS, VOL. 36, L22307, doi:10.1029/2009gl040484, 2009 Modification of the lithospheric stress field by lateral variations in plate-mantle coupling J. B. Naliboff, 1 C. P. Conrad,
More informationLithospheric Rheology and Stress, Dynamics of Plate Tectonics, and Long-wavelength Mantle Convection
Lithospheric Rheology and Stress, Dynamics of Plate Tectonics, and Long-wavelength Mantle Convection Shijie Zhong and Xi Liu Dept. of Physics, University of Colorado Boulder, Colorado, USA A. B. Watts,
More informationTrans-Pacific temperature field in the mantle transition region from seismic and electromagnetic tomography
Trans-Pacific temperature field in the mantle transition region from seismic and electromagnetic tomography Yoshio Fukao 1,3, Takao Koyama 2, Masayuki Obayashi 1 and Hisashi Utada 3 1 Research Program
More informationRecent tectonic plate decelerations driven by mantle convection
Click Here for Full Article GEOPHYSICAL RESEARCH LETTERS, VOL. 36, L23301, doi:10.1029/2009gl040224, 2009 Recent tectonic plate decelerations driven by mantle convection A. M. Forte, 1 R. Moucha, 1 D.
More informationRheology. What is rheology? From the root work rheo- Current: flow. Greek: rhein, to flow (river) Like rheostat flow of current
Rheology What is rheology? From the root work rheo- Current: flow Greek: rhein, to flow (river) Like rheostat flow of current Rheology What physical properties control deformation? - Rock type - Temperature
More informationConstraints on lithosphere and mantle rheology from in-situ observations
Constraints on lithosphere and mantle rheology from in-situ observations Shijie Zhong Department of Physics University of Colorado at Boulder Collaborators: Archie Paulson and John Wahr on post-glacial
More informationA mechanical model of the San Andreas fault and SAFOD Pilot Hole stress measurements
GEOPHYSICAL RESEARCH LETTERS, VOL. 31, L15S13, doi:10.1029/2004gl019521, 2004 A mechanical model of the San Andreas fault and SAFOD Pilot Hole stress measurements Jean Chéry Laboratoire Dynamique de la
More informationHigh-Harmonic Geoid Signatures due to Glacial Isostatic Adjustment, Subduction and Seismic Deformation
High-Harmonic Geoid Signatures due to Glacial Isostatic Adjustment, Subduction and Seismic Deformation L.L.A. Vermeersen (1), H. Schotman (1), M.-W. Jansen (1), R. Riva (1) and R. Sabadini (2) (1) DEOS,
More informationRheology III. Ideal materials Laboratory tests Power-law creep The strength of the lithosphere The role of micromechanical defects in power-law creep
Rheology III Ideal materials Laboratory tests Power-law creep The strength of the lithosphere The role of micromechanical defects in power-law creep Ideal materials fall into one of the following categories:
More informationThermo-chemical structure, dynamics and evolution of the deep mantle: spherical convection calculations
Thermo-chemical structure, dynamics and evolution of the deep mantle: spherical convection calculations Paul J. Tackley ETH Zürich, Switzerland With help from Takashi Nakagawa, Frédéric Deschamps, James
More informationGeodynamics of MARGINS
Geodynamics of MARGINS Scott King (Virginia Tech), Peter van Keken (Univ. Michigan), Mike Gurnis (Caltech), Marc Spiegelman (Lamont), Shijie Zhong (Univ. Colorado), Magali Billen (U.C. Davis) contact:
More informationDevelopment of Anisotropic Structure by Solid-State Convection in the Earth s Lower Mantle
1 Development of Anisotropic Structure by Solid-State Convection in the Earth s Lower Mantle Allen K. McNamara *, Peter E. van Keken, * & Shun-Ichiro Karato ** * Department of Geological Sciences, University
More informationDynamic processes in the mantle wedge
Mantle Processes and Geodynamics Dynamic processes in the mantle wedge Ikuko Wada Woods Hole Oceanographic Institution Something Amazing: Despite the apparent complexity of the mantle wedge dynamics, the
More informationMarine Geophysics. Plate tectonics. Dept. of Marine Sciences, Ocean College, Zhejiang University. Nov. 8, 2016
Marine Geophysics Plate tectonics 何小波 Dept. of Marine Sciences, Ocean College, Zhejiang University Nov. 8, 2016 Ocean College (ZJU) Plate tectonics xbhe@zju.edu.cn 1 / 1 Mantle flow and Plate tectonics
More informationCrust : wet quartzite Arc root : dry olivine mantle = 2840 kg/m km = 3300 kg/m km (Arc root thickness) 280 km (Arc width)
Crust : wet quartzite Arc root : dry olivine mantle = 2840 kg/m 3 41.5 km = 3300 kg/m 3 o 118.5 km Temperature (C) free surface z = 0 550 1350 160 km (Arc root thickness) 280 km (Arc width) sub-lithospheric
More informationRheological structure and deformation of subducted slabs in the mantle transition zone: implications for mantle circulation and deep earthquakes
Physics of the Earth and Planetary Interiors 127 (2001) 83 108 Rheological structure and deformation of subducted slabs in the mantle transition zone: implications for mantle circulation and deep earthquakes
More informationEarth and Planetary Science Letters
Earth and Planetary Science Letters 275 (2008) 43 53 Contents lists available at ScienceDirect Earth and Planetary Science Letters journal homepage: www.elsevier.com/locate/epsl Long-wavelength character
More informationSUPPLEMENTARY INFORMATION
DOI: 1.138/NGEO1521 Possible links between long-term geomagnetic variations and whole-mantle convection processes Biggin, A.J., Steinberger, B., Aubert, J., Suttie, N., Holme, R., Torsvik, T.H., van der
More informationLecture 24: Convergent boundaries November 22, 2006
Lecture 24: Convergent boundaries November 22, 2006 Convergent boundaries are characterized by consumption of oceaninc crust in subduction zones, with attendant arc volcanism, metamorphism and uplift.
More informationWhy Does Oceanic Crust Sink Beneath Continental Crust At Convergent Boundaries
Why Does Oceanic Crust Sink Beneath Continental Crust At Convergent Boundaries What is the process by which oceanic crust sinks beneath a deep-ocean Why does oceanic crust sink beneath continental crust
More informationPhysics and Chemistry of the Earth and Terrestrial Planets
MIT OpenCourseWare http://ocw.mit.edu 12.002 Physics and Chemistry of the Earth and Terrestrial Planets Fall 2008 For information about citing these materials or our Terms of Use, visit: http://ocw.mit.edu/terms.
More informationMaterial is perfectly elastic until it undergoes brittle fracture when applied stress reaches σ f
Material is perfectly elastic until it undergoes brittle fracture when applied stress reaches σ f Material undergoes plastic deformation when stress exceeds yield stress σ 0 Permanent strain results from
More informationBenchmarks for subduction zone models Subduction zone workshop, University of Michigan, July 2003
Introduction Benchmarks for subduction zone models Subduction zone workshop, University of Michigan, July 2003 In early October 2002 a group of researchers met at the University of Michigan at Ann Arbor,
More informationFORCES ON EARTH UNIT 3.2. An investigation into how Newton s Laws of Motion are applied to the tectonic activity on Earth.
FORCES ON EARTH UNIT 3.2 An investigation into how Newton s Laws of Motion are applied to the tectonic activity on Earth. USE THESE NOTES: OUR HOME PLANET EARTH: What do you know about our planet? SO.HOW
More informationThe Theory of Plate Tectonics
Plate Tectonics Objectives Describe how plates move. Explain the features of plate tectonics. Describe the types of plate boundaries and the features that can form and events that can occur at each. The
More informationThis article appeared in a journal published by Elsevier. The attached copy is furnished to the author for internal non-commercial research and
This article appeared in a journal published by Elsevier. The attached copy is furnished to the author for internal non-commercial research and education use, including for instruction at the authors institution
More informationWhat is the Rate-Limiting Step In furthering our Understanding of Subduction Dynamics?
What is the Rate-Limiting Step In furthering our Understanding of Subduction Dynamics? Magali Billen - Katrina Arredondo, - Juan Rodriguez, - Erin Burkett, --Margerete Jadamec Processes in Series are Rate-Limited
More informationEffects of compositional and rheological stratifications on smallscale
Click Here for Full Article GEOPHYSICAL RESEARCH LETTERS, VOL. 35, L20308, doi:10.1029/2008gl035419, 2008 Effects of compositional and rheological stratifications on small-scale convection under the oceans:
More informationSeismic Anisotropy and Mantle Flow in the Izu-Bonin-Mariana Subduction System
Seismic Anisotropy and Mantle Flow in the Izu-Bonin-Mariana Subduction System Matthew J. Fouch (Department of Geological Sciences, Arizona State University, Tempe, AZ 85287, email: fouch@asu.edu) INTRODUCTION
More informationTracing rays through the Earth
Tracing rays through the Earth Ray parameter p: source receiv er i 1 V 1 sin i 1 = sin i 2 = = sin i n = const. = p V 1 V 2 V n p is constant for a given ray i 2 i 3 i 4 V 2 V 3 V 4 i critical If V increases
More informationFORCES ON EARTH. An investigation into how Newton s Laws of Motion are applied to the tectonic activity on Earth.
FORCES ON EARTH An investigation into how Newton s Laws of Motion are applied to the tectonic activity on Earth. GEOLOGY Geologists scientists who study the forces that make and shape the Earth Geologists
More informationRecall Hypsometric Curve?
Structure of the Earth (Why?) Recall Hypsometric Curve? Continental crust is very different from oceanic crust. To understand this, we need to know more about the structure & composition of Earth. But
More informationLab 1: Plate Tectonics April 2, 2009
Name: Lab 1: Plate Tectonics April 2, 2009 Objective: Students will be introduced to the theory of plate tectonics and different styles of plate margins and interactions. Introduction The planet can be
More informationGEOL540: The Mantle System
GEOL540: The Mantle System Instructor Prof. Thorsten Becker (ZHS269; (213)740-8365; twb@usc.edu) Times Two 1.5 hours, 3 units Objective This graduate class discusses the structure, dynamics, and evolution
More informationIntroduction to Geology Spring 2008
MIT OpenCourseWare http://ocw.mit.edu 12.001 Introduction to Geology Spring 2008 For information about citing these materials or our Terms of Use, visit: http://ocw.mit.edu/terms. RHEOLOGICAL MODELS Rheology
More informationContinental Drift and Plate Tectonics
Continental Drift and Plate Tectonics Continental Drift Wegener s continental drift hypothesis stated that the continents had once been joined to form a single supercontinent. Wegener proposed that the
More informationIsostasy and Tectonics Lab Understanding the Nature of Mobile Floating Lithospheric Plates
Isostasy and Tectonics Lab Understanding the Nature of Mobile Floating Lithospheric Plates Crust Mantle Dynamics Introductory Geology Lab Ray Rector - Instructor Isostasy and Tectonics Laboratory Topics
More informationRecall Hypsometric Curve?
Structure of the Earth (Why?) 1 Recall Hypsometric Curve? Continental lithosphere is very different from oceanic lithosphere. To understand this, we need to know more about the structure & composition
More informationStructure of the Earth (Why?)
Structure of the Earth (Why?) 1 Recall Hypsometric Curve? Continental lithosphere is very different from oceanic lithosphere. To understand this, we need to know more about the structure & composition
More informationRheology: What is it?
Schedule Rheology basics Viscous, elastic and plastic Creep processes Flow laws Yielding mechanisms Deformation maps Yield strength envelopes Constraints on the rheology from the laboratory, geology, geophysics
More informationPlate Tectonics. entirely rock both and rock
Plate Tectonics I. Tectonics A. Tectonic Forces are forces generated from within Earth causing rock to become. B. 1. The study of the origin and arrangement of Earth surface including mountain belts, continents,
More informationEarth s Interior and Geophysical Properties. Chapter 13
Earth s Interior and Geophysical Properties Chapter 13 Introduction Can we just go there? Deep interior of the Earth must be studied indirectly Direct access only to crustal rocks and upper mantle fragments
More informationSeismic and flexure constraints on lithospheric rheology and their dynamic implications
Seismic and flexure constraints on lithospheric rheology and their dynamic implications Shijie Zhong Dept. of Physics, University of Colorado Boulder, Colorado, USA Acknowledgement: A. B. Watts Dept. of
More informationD) outer core B) 1300 C A) rigid mantle A) 2000 C B) density, temperature, and pressure increase D) stiffer mantle C) outer core
1. In which area of Earth's interior is the pressure most likely to be 2.5 million atmospheres? A) asthenosphere B) stiffer mantle C) inner core D) outer core Base your answers to questions 2 and 3 on
More informationSUPPLEMENTARY INFORMATION
The major uncertainties in our model predictions arise from the input parameters, which include mantle density models (i.e. seismic tomography and choices about scaling velocities to temperature), crustal
More informationRheology and the Lithosphere
Rheology and the Lithosphere Processes in Structural Geology & Tectonics Ben van der Pluijm WW Norton+Authors, unless noted otherwise 3/8/2017 16:51 We Discuss Rheology and the Lithosphere What is rheology?
More informationHysteresis in mantle convection: Plate tectonics systems
GEOPHYSICAL RESEARCH LETTERS, VOL. 39,, doi:10.1029/2012gl051232, 2012 Hysteresis in mantle convection: Plate tectonics systems M. B. Weller 1 and A. Lenardic 1 Received 5 February 2012; revised 3 April
More informationStress field in the subducting lithosphere and comparison with deep earthquakes in Tonga
JOURNAL OF GEOPHYSICAL RESEARCH, VOL. 108, NO. B6, 2288, doi:10.1029/2002jb002161, 2003 Stress field in the subducting lithosphere and comparison with deep earthquakes in Tonga A. Guest 1 and G. Schubert
More informationRecall Hypsometric Curve?
Structure of the Earth (Why?) Recall Hypsometric Curve? Continental lithosphere is very different from oceanic lithosphere. To understand this, we need to know more about the structure & composition of
More informationSIO15 Midterm 1, Monday Oct. 30, 2017 TEST VARIATION: 2
SIO15 Midterm 1, Monday Oct. 30, 2017 TEST VARIATION: 2 1) Over the last 20 years, which of these natural disasters killed the most people in the U.S.? a) heat waves b) hurricanes c) lightning strikes
More information4 Deforming the Earth s Crust
CHAPTER 7 4 Deforming the Earth s Crust SECTION Plate Tectonics BEFORE YOU READ After you read this section, you should be able to answer these questions: What happens when rock is placed under stress?
More informationMost mafic magmas come from the upper mantle and lower crust. This handout will address five questions:
Geology 101 Origin of Magma From our discussions of the structure of the interior of the Earth, it is clear that the upper parts of the Earth (crust and mantle) are mostly solid because s-waves penetrate
More informationInternal Layers of the Earth
Lecture #4 notes Geology 3950, Spring 2006; CR Stern Seismic waves, earthquake magnitudes and location, and internal earth structure (pages 28-95 in the 4 th edition and 28-32 and 50-106 in the 5 th edition)
More information