Subduction II Fundamentals of Mantle Dynamics

Size: px
Start display at page:

Download "Subduction II Fundamentals of Mantle Dynamics"

Transcription

1 Subduction II Fundamentals of Mantle Dynamics Thorsten W Becker University of Southern California Short course at Universita di Roma TRE April 18 20, 2011

2

3 Rheology Elasticity vs. viscous deformation η = O (1021) Pa s = viscosity µ= O (1011) Pa = shear modulus = rigidity τ = η / µ = O(1010) sec = O(103) years = Maxwell time

4 Elastic deformation In general: σij = Cijkl εkl (in 3-D 81 degrees of freedom, in general 21 independent) For isotropic body this reduces to Hooke s law: σij = λεkkδij + 2µεij with λ and µ Lame s parameters, εkk Taking shear components ( σ12 = 2µε12 = ε11 + ε22 + ε33 i j ) gives definition of rigidity: Adding the normal components ( i=j ) for all i=1,2,3 gives: σkk = (3λ + 2µ) εkk = 3κεkk with κ = λ + 2µ/3 = bulk modulus

5 Linear viscous deformation (1) Total stress field = static + dynamic part: σ ij = pδ ij + τ ij Analogous to elasticity General case: τ ij = C 'ijkl ε kl Isotropic case: τ ij = λ ' ε kk δ ij + 2ηε ij

6 Linear viscous deformation (2) Split in isotropic and deviatoric part (latter causes deformation): 1 σ 'ij = σ ij σ kk δ ij = σ ij + pδ ij 3 1 ε 'ij = ε ij ε kk δ ij 3 which gives the following stress: σ 'ij = ( p p )δ ij + ςε kkδ ij + 2ηε 'ij With compressibility term assumed 0 (Stokes condition 2 ( ς = λ ' + η = bulk viscosity) 3 τ Now τ = 2ηε or η = ij ij ij 2ε ij In general η = f (T, d, p, H 2O ) ε kk = 0 )

7 Non-linear (or non-newtonian) deformation General stress-strain relation: n = 1 : Newtonian n > 1 : non-newtonian n : pseudo-brittle Effective viscosity ηeff ε = Aτ n n 1 1 / n (1 n ) / n = A τ = A ε 2 2 Application: different viscosities under oceans with different absolute plate motion, anisotropic viscosities by means of superposition (Schmeling, 1987)

8 Microphysical theory and observations Maximum strength of materials (1) Strength is maximum stress that material can resist In principle, viscous fluid has zero strength. In reality, all materials have finite strength. (Ranalli, 1995) Elastic deformation until atom jumps to next equilibrium position. So theoretical strength σ = O( µ )

9 Microphysical theory and observations Maximum strength of materials (2) However, from laboratory measurements: Shear strength = O (10 Structural flaws Cracks Vacancies Dislocations Subgrain boundaries 4 µ) due to: (Ranalli, 1995) Stress concentration makes deformation possible under smaller stresses (compare to breaking/tearing of sheet of paper)

10 Diffusion creep (Schubert, Turcotte & Olson, 2001)

11 Dislocation creep (Ranalli, 1995)

12 Steady state creep models Theoretically many different models: only a few relevant for Earth (Climb-controlled) dislocation creep or powerlaw creep: gliding of dislocations controlled by the climb rate around n impurities/obstacles: E p V D SD D SD =D 0 exp RT Diffusion creep (Newtonian, linear creep): Nabarro-Herring creep (diffusion through grain) Coble creep (diffusion along grain boundary) D SD m d m=2-3 Peierl s stress mechanism (low-t plasticity): dislocation glide Grain-boundary sliding (superplasticity) Pressure-solution Brittle deformation, Byerlee s law

13 Strength of the Lithosphere and Mantle Experimental data for olivine

14 Strength of the Lithosphere and Mantle (Turcotte and Schubert, 2002) Collection of rheological data for different materials

15 See Hirth & Kohlstedt (2005) for olivine in the upper mantle

16 Laboratory experiments: large temperature, strain-rate & grain-size dependence.

17 Dislocation creep decreases viscosity where the strain-rate is more than the transition value.

18 Non-deforming regions remain highly viscous. Yielding concentrates deformation.

19 Deformation maps (Ranalli, 1995)

20 Strength of the Lithosphere and Mantle (3) (Kohlstedt et al., 1995) Strength curves for different materials: lithosphere

21 Slab rheology Figure courtesy of M. Billen

22 Governing equations Conservation of mass: continuity Conservation of momentum: Stokes equation Equation of state: density constitutive equation: rheology Conservation of energy: temperature

23 Continuity (from Turcotte and Schubert, 2002) This gives: u = 0

24 Derivation of Stokes equation Static force balance f = 0 ij f i =0 xj Incompressible constitutive relationship and separation into deviatoric stress ij = ij p ij =2 ij p ij Constant viscosity ij p 2 f i =0 x j xi

25 Derivation of Stokes equation (2) 2 ij vi vi v j vj 2 = = x j x j x j xi x j x j xi x j 2 ij vi 2 2 = = v i x j x j x j p v i f i =0 xi Because for incompressible vj =0 xj 2 2 vi p f i =0 x j x j xi p f = v 0 2

26 Incompressible, Newtonian flow

27

28

29 Simplified equations

30 Simple 1-D fluid dynamics examples Couette channel flow: Channel flow with horizontal pressure gradient (Hagen Poiseuille):

31 Stokes sinker solution

32 Inferences based on Stokes flow u u L L

33 Scaling: Dimensional equations: u = 0 p + η 2u = ρ 0α (T T0 ) gδ i 3 ρ 0C p dt = k 2T ρc pu T dt Scaling parameters x' = x / h t ' = t /(h 2 / κ ) T ' = (T T0 ) / T η ' = η / η0 Scaled equations (with primes left out): u = 0 p + η u = RaTδ i 3 2 dt = 2T u T dt αρ 0 g Th 3 with Rayleigh number Ra = ηκ

34 Rayleigh numbers Bottom heated Internal heating αρ0 g Th Ra = ηκ H h =k T h =k T h 2 Hh T = k 5 0 g H h Ra H = k

35 Effect of phase transitions Clapeyron slope dp dt = 0 g T h Ra = c 3 Q latent 1 2 T a g T h Rb= = 3 Buoyancy parameter Rb P 2 g h Ra = = c = 0 g h c T = e.g. Schubert et al. (1975); Christensen (1985); Schubert et al. (2001), p. 466f

36 Linear stability analysis (1) for small T (or Ra): conduction for larger Ra: convection sets in so minimum Ra = Rac exists below which no convection occurs T = Tconductive + T1 = T0 + T1 For small velocity and T1 we can linearize system: linear stability analysis Energy equation in terms of T1 and Ψ Remove small terms Solve with separation of variables: Ψ = Ψ * cos(nπz ) sin(kx ) exp(αt ) T1 = T1 cos(nπz ) cos(kx ) exp(αt ) * α<0: stable α>0: instable α=0: marginal stability Rac = (k n π ) k2 See Schubert et al. (2001) p. 288ff

37 Linear stability analysis (2) minimum Rac = for free-slip, bottom heated, c = 2.8 Rac as function of horizontal wavelength k=2πb/λ, b = box height, i.e. one wavelength of counter-rotating cells with aspect ratio 2.8 x 1 for = c (from Turcotte and Schubert, 2002)

38 Simple convection model for high Ra (from (Davies, 1999)) B = gddρα T R = 4ηv R+B=0 v~14 cm/yr relates physical quantities

39 Boundary layer theory u Ra h 2 3 Q T Ra Nu Ra See Schubert et al. (2001) p. 353ff Grigne et al. (2005)

40 Example: Boussinesq convection, isoviscous, no internal heating Ra=1 Nu=1 Figure courtesy of P. van Keken

41 Ra=103 Nu>1 Figure courtesy of P. van Keken

42 Ra=104 Nu=3.5 Figure courtesy of P. van Keken

43 Ra=105 Nu=8 Figure courtesy of P. van Keken

44 Ra=106 <Nu>=16 Figure courtesy of P. van Keken

45 Ra=2x106 <Nu>=18 Figure courtesy of P. van Keken

46 Ra=5x106 <Nu>=20 Figure courtesy of P. van Keken

47 Logistical equation (May, Nature, 1976): xn+1 = r xn ( 1 xn )

48 Ra > Rac modeling results higher Ra gives thinner thermal boundary layers For larger Ra flow usually not steady state

49 2D convection experiments, isoviscous case Ra = 107 H=0 '(T)= 1 lm/ um = 1 Courtesy of Allen McNamara

50 Temperature dependent viscosity Ra = 107 H=0 '(T) = 1000 lm/ um = 1 Courtesy of Allen McNamara

51 Thermal convection in the mantle Ra = 107 H=0 '(T) = 1000 lm/ um = 50 Courtesy of Allen McNamara

52 Increase of internal heating Ra = 2.4e5

53 Heating mode from (Davies, 1999) bottom/internal heating passive/active upwellings (MOR?) time dependence bottom/top boundary layer independent (plumes vs. plates)

54 Convection with η(t) (1) large aspect ratios large viscosity variations in top 200 km asymmetry between up- & downwellings from Bercovici et al., (1996)

55 Convection with η(t) (2) Stagnant lid regime on Earth? missing rheology! from Bercovici et al. (1996)

Chapter 1. Continuum mechanics review. 1.1 Definitions and nomenclature

Chapter 1. Continuum mechanics review. 1.1 Definitions and nomenclature Chapter 1 Continuum mechanics review We will assume some familiarity with continuum mechanics as discussed in the context of an introductory geodynamics course; a good reference for such problems is Turcotte

More information

Benchmarks for subduction zone models Subduction zone workshop, University of Michigan, July 2003

Benchmarks 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 information

Rheology of the Mantle and Plates (part 1): Deformation mechanisms and flow rules of mantle minerals

Rheology 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 information

DETAILS 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 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 information

1 Scaling analysis and non-dimensional numbers

1 Scaling analysis and non-dimensional numbers 1 Scaling analysis and non-dimensional numbers While this is a textbook on numerical analysis, it is crucial to keep the nature of the physical processes which we would like to model in mind. This will

More information

Why does the Nazca plate slow down since the Neogene? Supplemental Information

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 information

Lecture 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 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 information

Scaling Laws for Convection with. Temperature-dependent Viscosity and. Grain-damage

Scaling Laws for Convection with. Temperature-dependent Viscosity and. Grain-damage Scaling Laws for Convection with arxiv:1410.7652v1 [physics.flu-dyn] 28 Oct 2014 Temperature-dependent Viscosity and Grain-damage Bradford J. Foley 1,2 and David Bercovici 1 1 Yale University, Department

More information

Rheology: What is it?

Rheology: 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 information

Rheology. 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 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 information

Rock Rheology GEOL 5700 Physics and Chemistry of the Solid Earth

Rock Rheology GEOL 5700 Physics and Chemistry of the Solid Earth Rock Rheology GEOL 5700 Physics and Chemistry of the Solid Earth References: Turcotte and Schubert, Geodynamics, Sections 2.1,-2.4, 2.7, 3.1-3.8, 6.1, 6.2, 6.8, 7.1-7.4. Jaeger and Cook, Fundamentals of

More information

Dynamics of the Mantle and Lithosphere ETH Zürich Continuum Mechanics in Geodynamics: Equation cheat sheet

Dynamics of the Mantle and Lithosphere ETH Zürich Continuum Mechanics in Geodynamics: Equation cheat sheet Dynamics of the Mantle and Lithosphere ETH Zürich Continuum Mechanics in Geodynamics: Equation cheat sheet or all equations you will probably ever need Definitions 1. Coordinate system. x,y,z or x 1,x

More information

Seismic and flexure constraints on lithospheric rheology and their dynamic implications

Seismic 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 information

Introduction to Geology Spring 2008

Introduction 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 information

EART162: PLANETARY INTERIORS

EART162: PLANETARY INTERIORS EART162: PLANETARY INTERIORS Francis Nimmo Last Week Global gravity variations arise due to MoI difference (J 2 ) We can also determine C, the moment of inertia, either by observation (precession) or by

More information

Fundamentals of Fluid Dynamics: Elementary Viscous Flow

Fundamentals of Fluid Dynamics: Elementary Viscous Flow Fundamentals of Fluid Dynamics: Elementary Viscous Flow Introductory Course on Multiphysics Modelling TOMASZ G. ZIELIŃSKI bluebox.ippt.pan.pl/ tzielins/ Institute of Fundamental Technological Research

More information

Tutorials for Thermal Convection. Shijie Zhong Department of Physics University of Colorado Boulder, CO

Tutorials for Thermal Convection. Shijie Zhong Department of Physics University of Colorado Boulder, CO Tutorials for Thermal Convection Shijie Zhong Department of Physics University of Colorado Boulder, CO 80305 Email: szhong@colorado.edu Table of content: 1. Introduction 2 2. Governing equations 2 3. Rheology

More information

Thermal-Mechanical Behavior of Oceanic Transform Faults

Thermal-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 information

Journal of Applied Mathematics and Computation (JAMC), 2018, 2(7),

Journal 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 information

SUPPLEMENTARY INFORMATION

SUPPLEMENTARY 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 information

Gravity Tectonics Volcanism Atmosphere Water Winds Chemistry. Planetary Surfaces

Gravity 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 information

Stress, Strain, and Viscosity. San Andreas Fault Palmdale

Stress, Strain, and Viscosity. San Andreas Fault Palmdale Stress, Strain, and Viscosity San Andreas Fault Palmdale Solids and Liquids Solid Behavior: Liquid Behavior: - elastic - fluid - rebound - no rebound - retain original shape - shape changes - small deformations

More information

Fluid Mechanics II Viscosity and shear stresses

Fluid Mechanics II Viscosity and shear stresses Fluid Mechanics II Viscosity and shear stresses Shear stresses in a Newtonian fluid A fluid at rest can not resist shearing forces. Under the action of such forces it deforms continuously, however small

More information

Recycling the lid: Effects of subduction and stirring on boundary layer dynamics in bottom-heated planetary mantle convection

Recycling the lid: Effects of subduction and stirring on boundary layer dynamics in bottom-heated planetary mantle convection Click Here for Full Article GEOPHYSICAL RESEARCH LETTERS, VOL. 33, L20318, doi:10.1029/2006gl027668, 2006 Recycling the lid: Effects of subduction and stirring on boundary layer dynamics in bottom-heated

More information

Lithospheric 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 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 information

Three-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. 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 information

Lecture 8: Tissue Mechanics

Lecture 8: Tissue Mechanics Computational Biology Group (CoBi), D-BSSE, ETHZ Lecture 8: Tissue Mechanics Prof Dagmar Iber, PhD DPhil MSc Computational Biology 2015/16 7. Mai 2016 2 / 57 Contents 1 Introduction to Elastic Materials

More information

Viscoelasticity. Basic Notions & Examples. Formalism for Linear Viscoelasticity. Simple Models & Mechanical Analogies. Non-linear behavior

Viscoelasticity. Basic Notions & Examples. Formalism for Linear Viscoelasticity. Simple Models & Mechanical Analogies. Non-linear behavior Viscoelasticity Basic Notions & Examples Formalism for Linear Viscoelasticity Simple Models & Mechanical Analogies Non-linear behavior Viscoelastic Behavior Generic Viscoelasticity: exhibition of both

More information

The anisotropic and rheological structure of the oceanic upper mantle from a simple model of plate shear

The anisotropic and rheological structure of the oceanic upper mantle from a simple model of plate shear Geophys. J. Int. (24) 158, 287 296 doi: 1.1111/j.1365-246X.24.225.x The anisotropic and rheological structure of the oceanic upper mantle from a simple model of plate shear Noah S. Podolefsky, Shijie Zhong

More information

Supplementary information on the West African margin

Supplementary 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 information

Time-varying subduction and rollback velocities in slab stagnation and buckling

Time-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 information

Modeling the Thermal-Mechanical Behavior of Mid-Ocean Ridge Transform Faults

Modeling 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 information

Viscous heating, adiabatic heating and energetic consistency in compressible mantle convection

Viscous heating, adiabatic heating and energetic consistency in compressible mantle convection Geophys. J. Int. (28) 173, 693 72 doi: 1.1111/j.1365-246X.28.3745.x Viscous heating, adiabatic heating and energetic consistency in compressible mantle convection Wei Leng and Shijie Zhong Department of

More information

1 Exercise: Linear, incompressible Stokes flow with FE

1 Exercise: Linear, incompressible Stokes flow with FE Figure 1: Pressure and velocity solution for a sinking, fluid slab impinging on viscosity contrast problem. 1 Exercise: Linear, incompressible Stokes flow with FE Reading Hughes (2000), sec. 4.2-4.4 Dabrowski

More information

Lithospheric-Scale Stresses and Shear Localization Induced by Density-Driven Instabilities

Lithospheric-Scale Stresses and Shear Localization Induced by Density-Driven Instabilities Examensarbete vid Institutionen för geovetenskaper ISSN 1650-6553 Nr 187 Lithospheric-Scale Stresses and Shear Localization Induced by Density-Driven Instabilities Christiane Heinicke Copyright c Christiane

More information

Modeling the Dynamics of Subducting Slabs

Modeling 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 information

Rheology 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 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 information

10. Buoyancy-driven flow

10. Buoyancy-driven flow 10. Buoyancy-driven flow For such flows to occur, need: Gravity field Variation of density (note: not the same as variable density!) Simplest case: Viscous flow, incompressible fluid, density-variation

More information

CH.9. CONSTITUTIVE EQUATIONS IN FLUIDS. Multimedia Course on Continuum Mechanics

CH.9. CONSTITUTIVE EQUATIONS IN FLUIDS. Multimedia Course on Continuum Mechanics CH.9. CONSTITUTIVE EQUATIONS IN FLUIDS Multimedia Course on Continuum Mechanics Overview Introduction Fluid Mechanics What is a Fluid? Pressure and Pascal s Law Constitutive Equations in Fluids Fluid Models

More information

Numerical models of convection in a rheologically stratified oceanic upper mantle: Early results

Numerical models of convection in a rheologically stratified oceanic upper mantle: Early results Earth Planets Space, 50, 1047 1054, 1998 Numerical models of convection in a rheologically stratified oceanic upper mantle: Early results Mamoru Kato Department of Geophysics, Kyoto University, Kyoto 606-8502,

More information

Abstract. Because mantle viscosity is temperature-dependent, cold subducting

Abstract. Because mantle viscosity is temperature-dependent, cold subducting Chapter 6 Mantle Convection with Strong Subduction Zones Abstract. Because mantle viscosity is temperature-dependent, cold subducting lithosphere should be strong, which implies that the rapid, localized

More information

RHEOLOGY & LINEAR ELASTICITY. B Importance of fluids and fractures in deformation C Linear elasticity for homogeneous isotropic materials

RHEOLOGY & LINEAR ELASTICITY. B Importance of fluids and fractures in deformation C Linear elasticity for homogeneous isotropic materials GG303 Lecture 2 0 9/4/01 1 RHEOLOGY & LINEAR ELASTICITY I II Main Topics A Rheology: Macroscopic deformation behavior B Importance of fluids and fractures in deformation C Linear elasticity for homogeneous

More information

Constraints on Mantle Structure from Surface Observables

Constraints 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 information

Effects of elasticity on the Rayleigh Taylor instability: implications for large-scale geodynamics

Effects of elasticity on the Rayleigh Taylor instability: implications for large-scale geodynamics Geophys. J. Int. (2006) doi: 10.1111/j.1365-246X.2006.03201.x Effects of elasticity on the Rayleigh Taylor instability: implications for large-scale geodynamics Boris J. P. Kaus and Thorsten W. Becker

More information

Lecture 5. Rheology. Earth Structure (2 nd Edition), 2004 W.W. Norton & Co, New York Slide show by Ben van der Pluijm

Lecture 5. Rheology. Earth Structure (2 nd Edition), 2004 W.W. Norton & Co, New York Slide show by Ben van der Pluijm Lecture 5 Rheology Earth Structure (2 nd Edition), 2004 W.W. Norton & Co, New York Slide show by Ben van der Pluijm WW Norton; unless noted otherwise Rheology is... the study of deformation and flow of

More information

Agricultural Science 1B Principles & Processes in Agriculture. Mike Wheatland

Agricultural Science 1B Principles & Processes in Agriculture. Mike Wheatland Agricultural Science 1B Principles & Processes in Agriculture Mike Wheatland (m.wheatland@physics.usyd.edu.au) Outline - Lectures weeks 9-12 Chapter 6: Balance in nature - description of energy balance

More information

Introduction to Marine Hydrodynamics

Introduction to Marine Hydrodynamics 1896 1920 1987 2006 Introduction to Marine Hydrodynamics (NA235) Department of Naval Architecture and Ocean Engineering School of Naval Architecture, Ocean & Civil Engineering First Assignment The first

More information

Geodynamics Lecture 10 The forces driving plate tectonics

Geodynamics 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 information

- Marine Hydrodynamics. Lecture 4. Knowns Equations # Unknowns # (conservation of mass) (conservation of momentum)

- Marine Hydrodynamics. Lecture 4. Knowns Equations # Unknowns # (conservation of mass) (conservation of momentum) 2.20 - Marine Hydrodynamics, Spring 2005 Lecture 4 2.20 - Marine Hydrodynamics Lecture 4 Introduction Governing Equations so far: Knowns Equations # Unknowns # density ρ( x, t) Continuity 1 velocities

More information

Contents. I Introduction 1. Preface. xiii

Contents. I Introduction 1. Preface. xiii Contents Preface xiii I Introduction 1 1 Continuous matter 3 1.1 Molecules................................ 4 1.2 The continuum approximation.................... 6 1.3 Newtonian mechanics.........................

More information

Eects of Plate Bending and Fault. Strength at Subduction Zones on. Published in Journal of Geophysical Research by C. P. Conrad and B. H.

Eects of Plate Bending and Fault. Strength at Subduction Zones on. Published in Journal of Geophysical Research by C. P. Conrad and B. H. Chapter 4 Eects of Plate Bending and Fault Strength at Subduction Zones on Plate Dynamics Published in Journal of Geophysical Research by C. P. Conrad and B. H. Hager, 104, 17551-17571, 1999. Copyright

More information

Beall et al., 2018, Formation of cratonic lithosphere during the initiation of plate tectonics: Geology, https://doi.org/ /g

Beall 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 information

Journal of Geophysical Research: Planets

Journal of Geophysical Research: Planets RESEARCH ARTICLE Key Points: Time variability of stresses and dynamic disequilibrium affect the initiation of plate tectonics Upper thermal boundary layer dynamics are not self-determined; the mantle below

More information

When you are standing on a flat surface, what is the normal stress you exert on the ground? What is the shear stress?

When you are standing on a flat surface, what is the normal stress you exert on the ground? What is the shear stress? When you are standing on a flat surface, what is the normal stress you exert on the ground? What is the shear stress? How could you exert a non-zero shear stress on the ground? Hydrostatic Pressure (fluids)

More information

Development of Anisotropic Structure by Solid-State Convection in the Earth s Lower Mantle

Development 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 information

Why cold slabs stagnate in the transition zone

Why 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 information

Controls on plume heat flux and plume excess temperature

Controls on plume heat flux and plume excess temperature Click Here for Full Article JOURNAL OF GEOPHYSICAL RESEARCH, VOL. 113,, doi:10.1029/2007jb005155, 2008 Controls on plume heat flux and plume excess temperature Wei Leng 1 and Shijie Zhong 1 Received 3

More information

Material 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 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 information

The combined effects of continents and the 660 km depth endothermic phase boundary on the thermal regime in the mantle

The combined effects of continents and the 660 km depth endothermic phase boundary on the thermal regime in the mantle 1 2 3 The combined effects of continents and the 66 km depth endothermic phase boundary on the thermal regime in the mantle 4 5 6 G. Sinha, a S.L. Butler a, a Department of Geological Sciences, University

More information

The effect of rheological parameters on plate behaviour in a self-consistent model of mantle convection

The effect of rheological parameters on plate behaviour in a self-consistent model of mantle convection Physics of the Earth and Planetary Interiors 142 (2004) 225 255 The effect of rheological parameters on plate behaviour in a self-consistent model of mantle convection Claudia Stein, Jörg Schmalzl, Ulrich

More information

V (r,t) = i ˆ u( x, y,z,t) + ˆ j v( x, y,z,t) + k ˆ w( x, y, z,t)

V (r,t) = i ˆ u( x, y,z,t) + ˆ j v( x, y,z,t) + k ˆ w( x, y, z,t) IV. DIFFERENTIAL RELATIONS FOR A FLUID PARTICLE This chapter presents the development and application of the basic differential equations of fluid motion. Simplifications in the general equations and common

More information

Dynamics of Glaciers

Dynamics of Glaciers Dynamics of Glaciers McCarthy Summer School 01 Andy Aschwanden Arctic Region Supercomputing Center University of Alaska Fairbanks, USA June 01 Note: This script is largely based on the Physics of Glaciers

More information

RHEOLOGY Principles, Measurements, and Applications. Christopher W. Macosko

RHEOLOGY Principles, Measurements, and Applications. Christopher W. Macosko RHEOLOGY Principles, Measurements, and Applications I -56081-5'79~5 1994 VCH Publishers. Inc. New York Part I. CONSTITUTIVE RELATIONS 1 1 l Elastic Solid 5 1.1 Introduction 5 1.2 The Stress Tensor 8 1.2.1

More information

The finite difference code (fully staggered grid) includes a marker-in-cell Lagrangian marker

The finite difference code (fully staggered grid) includes a marker-in-cell Lagrangian marker GSA Data Repository 2018289 Ruh et al., 2018, Shale-related minibasins atop a massive olistostrome in an active accretionary wedge setting: Two-dimensional numerical modeling applied to the Iranian Makran:

More information

Numerical Simulation of the Thermal Convection and Subduction Process in the Mantle

Numerical 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 information

Rheology of Soft Materials. Rheology

Rheology of Soft Materials. Rheology Τ Thomas G. Mason Department of Chemistry and Biochemistry Department of Physics and Astronomy California NanoSystems Institute Τ γ 26 by Thomas G. Mason All rights reserved. γ (t) τ (t) γ τ Δt 2π t γ

More information

Viscoelasticity in mantle convection

Viscoelasticity in mantle convection Mgr. et Mgr. Vojtěch Patočka Supervised by: RNDr. Ondřej Čadek, CSc. [patocka@karel.troja.mff.cuni.cz] 5th June 2015 Content Method of implementation Model Method Testing the method Thermal convection

More information

Geodynamics Lecture 5 Basics of elasticity

Geodynamics Lecture 5 Basics of elasticity Geodynamics Lecture 5 Basics of elasticity Lecturer: David Whipp david.whipp@helsinki.fi! 16.9.2014 Geodynamics www.helsinki.fi/yliopisto 1 Goals of this lecture Introduce linear elasticity! Look at the

More information

Geophysical Journal International

Geophysical Journal International Geophysical Journal International Geophys. J. Int. (2017) 208, 491 507 Advance Access publication 2016 October 14 GJI Geodynamics and tectonics doi: 10.1093/gji/ggw392 The effect of a power-law mantle

More information

Plate tectonics. Chapter Slab pull. Dynamics of the Mantle and Lithosphere

Plate 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 information

Dynamic processes in the mantle wedge

Dynamic 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 information

Seismotectonics 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 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 information

Exercise: concepts from chapter 10

Exercise: concepts from chapter 10 Reading:, Ch 10 1) The flow of magma with a viscosity as great as 10 10 Pa s, let alone that of rock with a viscosity of 10 20 Pa s, is difficult to comprehend because our common eperience is with s like

More information

Rheology and the Lithosphere

Rheology 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 information

Chapter 9: Differential Analysis of Fluid Flow

Chapter 9: Differential Analysis of Fluid Flow of Fluid Flow Objectives 1. Understand how the differential equations of mass and momentum conservation are derived. 2. Calculate the stream function and pressure field, and plot streamlines for a known

More information

Three-dimensional instantaneous mantle flow induced by subduction. Dip. Scienze Geologiche, Universita degli Studi Roma TRE, Rome, Italy

Three-dimensional instantaneous mantle flow induced by subduction. Dip. Scienze Geologiche, Universita degli Studi Roma TRE, Rome, Italy Three-dimensional instantaneous mantle flow induced by subduction C. Piromallo (1), T.W. Becker (2), F. Funiciello (3) and C. Faccenna (3) (1) Istituto Nazionale di Geofisica e Vulcanologia, Rome, Italy

More information

PEAT SEISMOLOGY Lecture 2: Continuum mechanics

PEAT SEISMOLOGY Lecture 2: Continuum mechanics PEAT8002 - SEISMOLOGY Lecture 2: Continuum mechanics Nick Rawlinson Research School of Earth Sciences Australian National University Strain Strain is the formal description of the change in shape of a

More information

Chapter 2: Fluid Dynamics Review

Chapter 2: Fluid Dynamics Review 7 Chapter 2: Fluid Dynamics Review This chapter serves as a short review of basic fluid mechanics. We derive the relevant transport equations (or conservation equations), state Newton s viscosity law leading

More information

What we should know about mechanics of materials

What we should know about mechanics of materials What we should know about mechanics of materials 0 John Maloney Van Vliet Group / Laboratory for Material Chemomechanics Department of Materials Science and Engineering Massachusetts Institute of Technology

More information

Thermal and compositional structure of the Mantle and Lithosphere

Thermal and compositional structure of the Mantle and Lithosphere Chapter 1 Thermal and compositional structure of the Mantle and Lithosphere 1.1 Primordial heat of the Earth The most widely accepted planetary formation theory says that the solar system accreted from

More information

Modeling the interior dynamics of terrestrial planets

Modeling 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 information

Lecture 2: Constitutive Relations

Lecture 2: Constitutive Relations Lecture 2: Constitutive Relations E. J. Hinch 1 Introduction This lecture discusses equations of motion for non-newtonian fluids. Any fluid must satisfy conservation of momentum ρ Du = p + σ + ρg (1) Dt

More information

Scaling Laws and the Problem of the Prediction of Tectonic Modes

Scaling Laws and the Problem of the Prediction of Tectonic Modes Scaling Laws and the Problem of the Prediction of Tectonic Modes D. Breuer 1 and V. Stamenkovic 2 1 DLR, Institute of Planetary Research, Berlin 2 MIT, Earth, Atmospheric and Planetary Sciences PLATE TECTONICS

More information

20. Rheology & Linear Elasticity

20. Rheology & Linear Elasticity I Main Topics A Rheology: Macroscopic deformation behavior B Linear elasticity for homogeneous isotropic materials 10/29/18 GG303 1 Viscous (fluid) Behavior http://manoa.hawaii.edu/graduate/content/slide-lava

More information

Chapter 9: Differential Analysis

Chapter 9: Differential Analysis 9-1 Introduction 9-2 Conservation of Mass 9-3 The Stream Function 9-4 Conservation of Linear Momentum 9-5 Navier Stokes Equation 9-6 Differential Analysis Problems Recall 9-1 Introduction (1) Chap 5: Control

More information

GEOL540: The Mantle System

GEOL540: 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 information

Elements of Rock Mechanics

Elements of Rock Mechanics Elements of Rock Mechanics Stress and strain Creep Constitutive equation Hooke's law Empirical relations Effects of porosity and fluids Anelasticity and viscoelasticity Reading: Shearer, 3 Stress Consider

More information

Basic concepts in viscous flow

Basic concepts in viscous flow Élisabeth Guazzelli and Jeffrey F. Morris with illustrations by Sylvie Pic Adapted from Chapter 1 of Cambridge Texts in Applied Mathematics 1 The fluid dynamic equations Navier-Stokes equations Dimensionless

More information

Physics of Continuous media

Physics of Continuous media Physics of Continuous media Sourendu Gupta TIFR, Mumbai, India Classical Mechanics 2012 October 26, 2012 Deformations of continuous media If a body is deformed, we say that the point which originally had

More information

On the effect of temperature and strain-rate dependent viscosity on global mantle flow, net rotation, and plate-driving forces

On the effect of temperature and strain-rate dependent viscosity on global mantle flow, net rotation, and plate-driving forces Geophys. J. Int. (2006) 167, 943 957 doi: 10.1111/j.1365-246X.2006.03172.x On the effect of temperature and strain-rate dependent viscosity on global mantle flow, net rotation, and plate-driving forces

More information

Les Houches School of Foam: Rheology of Complex Fluids

Les Houches School of Foam: Rheology of Complex Fluids Les Houches School of Foam: Rheology of Complex Fluids Andrew Belmonte The W. G. Pritchard Laboratories Department of Mathematics, Penn State University 1 Fluid Dynamics (tossing a coin) Les Houches Winter

More information

2 Law of conservation of energy

2 Law of conservation of energy 1 Newtonian viscous Fluid 1 Newtonian fluid For a Newtonian we already have shown that σ ij = pδ ij + λd k,k δ ij + 2µD ij where λ and µ are called viscosity coefficient. For a fluid under rigid body motion

More information

6.2 Governing Equations for Natural Convection

6.2 Governing Equations for Natural Convection 6. Governing Equations for Natural Convection 6..1 Generalized Governing Equations The governing equations for natural convection are special cases of the generalized governing equations that were discussed

More information

Continuum Mechanics. Continuum Mechanics and Constitutive Equations

Continuum Mechanics. Continuum Mechanics and Constitutive Equations Continuum Mechanics Continuum Mechanics and Constitutive Equations Continuum mechanics pertains to the description of mechanical behavior of materials under the assumption that the material is a uniform

More information

Fundamentals of Linear Elasticity

Fundamentals of Linear Elasticity Fundamentals of Linear Elasticity Introductory Course on Multiphysics Modelling TOMASZ G. ZIELIŃSKI bluebox.ippt.pan.pl/ tzielins/ Institute of Fundamental Technological Research of the Polish Academy

More information

Constraints on thermochemical convection of the mantle from plume heat flux, plume excess temperature, and upper mantle temperature

Constraints on thermochemical convection of the mantle from plume heat flux, plume excess temperature, and upper mantle temperature JOURNAL OF GEOPHYSICAL RESEARCH, VOL. 111,, doi:10.1029/2005jb003972, 2006 Constraints on thermochemical convection of the mantle from plume heat flux, plume excess temperature, and upper mantle temperature

More information

Earth and Planetary Science Letters

Earth and Planetary Science Letters Earth and Planetary Science Letters 297 (2010) 154 164 Contents lists available at ScienceDirect Earth and Planetary Science Letters journal homepage: www.elsevier.com/locate/epsl Constraints on viscous

More information

Geophysical Journal International

Geophysical Journal International Geophysical Journal International Geophys. J. Int. (2010) 180, 23 33 doi: 10.1111/j.1365-246X.2009.04404.x Long wavelength convection, Poiseuille Couette flow in the low-viscosity asthenosphere and the

More information

Interfacial dynamics

Interfacial dynamics Interfacial dynamics Interfacial dynamics = dynamic processes at fluid interfaces upon their deformation Interfacial rheological properties: elasticity, viscosity, yield stress, Relation between macroscopic

More information

Chapter 7: Natural Convection

Chapter 7: Natural Convection 7-1 Introduction 7- The Grashof Number 7-3 Natural Convection over Surfaces 7-4 Natural Convection Inside Enclosures 7-5 Similarity Solution 7-6 Integral Method 7-7 Combined Natural and Forced Convection

More information