SUPPLEMENTARY INFORMATION
|
|
- Bridget Morrison
- 5 years ago
- Views:
Transcription
1 SUPPLEMENTARY INFORMATION Linking mantle upwelling with the lithosphere decent and the Japan Sea evolution: a hypothesis Alik Ismail-Zadeh 1, 2, 3 *, Satoru Honda 4, Igor Tsepelev Institut für Angewandte Geowissenschaften, Karlsruher Institut für Technologie, Karlsruhe, GERMANY Institut de Physique du Globe de Paris, Paris, FRANCE Institute of Earthquake Prediction Theory and Mathematical Geophysics, Russian Academy of Sciences, Moscow, RUSSIA Earthquake Research Institute, University of Tokyo, Tokyo, JAPAN Institute of Mathematics and Mechanics, Russian Academy of Sciences, Yekaterinburg, RUSSIA * Correspondence to Alik.Ismail-Zadeh@kit.edu or aiz@ipgp.fr Here we describe the detail of the mathematical model including the governing equations, the rheological law, the background temperature, and the boundary conditions. Numerical methods and solvers are then presented. The performance of the sensitivity analysis is illustrated by several case studies related to (i) variations in mantle viscosity, (ii) model depth variation, and (iii) the presence of phase changes in the mantle. Also here we compare qualitatively three seismic tomography models.
2 MODEL DESCRIPTION Mathematical statement In the three-dimensional (3-D) rectangular domain [0, x1 l1] [0, x2 l2] [0, x3 h] and for time interval t [0, ] we solve the regularized Stokes, the incompressibility, and the backward heat balance equations using the quasi-reversibility method 42 and the extended Boussinesq approximation 10,1 : T P ( [ u ( u) ]) ( E ) RaT a La ( ) a La ( ) e, (1) u 0, (2) * 1 2 * 2 ( E ) T u T A B Di Ra u3 T A T Di ( eij ) (3) t i, j 1 with appropriate boundary and initial conditions (see below). Here * A 1 a a Di La T 0, w1 w2 La 2 2 B 1 a a Ra w w , 1 1 tanh i i, i zi x3 i( T Ti), i = 1, 2, 2 w i x = (x 1, x 2, x 3 ), u = (u 1, u 2, u 3 ), t, T, P, and are the dimensionless Cartesian coordinates, velocity, time, temperature, pressure, and viscosity, respectively; is the present time; e ( u ) { u / x u / x} is the strain rate tensor; e (0,0,1) is the unit vector; is the ij i j j i gradient operator; and E is the unit operator. With regard to the phase changes around 410 km and 0 km, respectively, 1 and 2 are the dimensionless excess pressures; and 1 and 2 are the phase functions describing the relative fraction of the heavier phase, respectively, and varying between 0 and 1 as a function of depth and temperature. The Rayleigh (Ra), Laplace (La), and 2
3 modified dissipation (Di * ) dimensionless numbers are defined as * * 3 * 1 Ra g T h ( ), La gh ( ), and * 3 * 1 Di ( ch T ) * * * 2 * 1, respectively. The operator ( ) 2 1 E is applied to the right-hand side of the Stokes equations (1) to smooth temperature jumps at the phase boundaries and to enhance the stability of our computations. The higher dissipation term, whose magnitude is controlled by the small parameter, is added in order to regularize the heat balance equation (3). Length, temperature, and time are normalized by h, respectively. The physical parameters used in this study are listed in Table S1. T, and h 2 1, We consider a temperature- and depth-dependent Newtonian rheology ( T( x ), x3) : E gxv RT, exp a * 3 a 0 where 0 is determined so that it will give Pa s at the depth of 290 km and temperature of 198 K. Other parameters of the rheological law are listed in Table S1. We set the upper limit of the viscosity at ~10 22 Pa s, which results in the viscosity increase from the upper to the lower mantle by about two orders of magnitude. Background temperature The initial conditions for temperature are determined from the seismic tomography model as described in the main text. The background temperature used for deriving the temperature from seismic velocity anomalies is presented in Fig. S1. Note that the solidus for dry rocks (black curve) is used in this derivation. 3
4 Fig. S1. Background temperature used in the model (the light blue curve for the temperature determined from the cooling half-space model 18, and the dark blue curve for the temperature from the geothermal adiabat 19 ). The liquidus (brown curve 2 ) and solidus for dry (black curve 18 ) and wet (red and orange curves 3 ) rocks are also presented. Boundary conditions At the upper surface of the model boundary we prescribe the velocity W (, tx1, x2, x 3 0) and temperature T = T u. The velocity W is constructed as described in the main text. Due to complexities and significant uncertainties in trench migration (oceanward versus trenchward) patterns observed on the wide Pacific and Philippine Sea subduction zone 4, we do not include the trench migration in the present model. Also we do not introduce a weak (artificial) zone at the plate interfaces as it is used in forward modelling to reduce stresses at the interfaces and to promote subduction. Such a weak zone is partly accommodated in backward modelling, because in the backward (time-reverse) modeling the hotter material from below moves upward toward the divergent plate interfaces reducing the viscosity there. At the lower surface of the model boundary we set the velocity u 0 (no-slip) and fixed temperature T = T l. To allow for the flow to pass through the lateral boundaries, we prescribe at the lateral sides of the boundary the following conditions: u T 0, 0. 4
5 P Also we prescribe 0 at the model boundary. We tested several other conditions at the lateral sides and at the bottom of the model domain (see Table S2) to study how the boundary conditions influence the restoration results. We found that the major features related to the present conclusions are reserved among the above different cases. NUMERICAL METHODS AND SOLVERS The governing equations with the prescribed boundary and initial conditions are solved numerically by the finite-volume method 10 using open source computational fluid dynamics software package OpenFoam ( We use finite volumes (rectangular hexahedrons), and hence a horizontal resolution of the model is 20 km 20 km. The model domain is divided into five horizontal layers: Layer 1 (from the surface to the depth of 400 km), layer 2 ( km), layer 3 ( km), layer 4 (50-70 km), and layer 5 (70 to 800 km). Within each layer 0, 40, 35, 40, and 15 grid points are used. Therefore, a vertical resolution of the model varies from 0.5 to 8.7 km. The accuracy of the numerical solutions has been verified by several tests including grid changes, volume preservation, and principle of the maximum 10,5. Several approaches have been developed to reconstruct the past thermal state and flow in the crust and mantle: backward advection method -9, sequential filtering method 70, variational/ adjoint method 71-73, and quasi-reversibility (QRV) method 43. Among these methods, the QRV method is less susceptible to a noise (small temperature perturbations) in restoration models 74. The QRV method for data assimilation introduces a new term in the heat balance equation (the first term in Eq. 3) to regularize the equation when solving it backward in time. The additional term describes heat flux relaxation 43,75. Velocity u and pressure P are found from the equations (1) and (2) using the SIMPLE method 10,7. The regularized heat balance equation (3) is approximated by the Euler method using the implicit approximation of the advective term and the explicit approximation of the conductive term: 5
6 n 1 n 2 T T n 1 n n ( E D) CT DT f ( u, T ) 0, dt where the discrete operators T C= C and T D=D approximate the advective and conductive terms, respectively. To solve the numerical scheme we use the splitting method 77 introducing the convection/anti-diffusion and regularization parts as n 1/2 n n ( E dtc) T ( E dtd) T dt f ( u, T ), (4) 2 n 1 n 1/2 ( E D ) T T. (5) The system of the discrete equations (4) is solved by the BiConjugate Gradient method 78 using the incomplete LU-factorization as a pre-conditioner 79. The system (5) is solved by the conjugate gradient method 10. NUMERICAL TESTS Table S3 lists the case studies performed to analyse the effects of (i) the variation in mantle viscosity (case studies ), (ii) the variation in the depth of the model domain (case studies 2.1 and 2.2.), and (iii) the phase changes (case studies ). In the following results shown in Figs. S2 S8, the figures present a 3-D (south to north) view of the temperature anomaly and viscosity in the region ( 4000 km km depth ). The region presents a central part of the model domain truncated equally from the north and from the south. A transparent mode of visualization with transparency coefficient 0.5 is used. Variations in mantle viscosity To analyse a sensitivity of the model to variations in mantle viscosity, we performed three case studies varying the dimensionless viscosity from 45 to 591 by changing the activation volume V a : m 3 mol -1 (case study 1.1, Fig. S2), m 3 mol -1 (case study 1.2, Fig. S3), and m 3 mol -1 (case study 1.3, Fig. S4). We found that small-scale upwellings in the model show persistency despite the realistic variations in viscosity of the mantle and the slab. Particularly, the upwelling beneath the northern Japan Sea linking the hotter material in the subslab mantle with the mantle wedge region is clearly recognized in Figs. S3 and S4 (see areas marked by green ovals in the figures).
7 Fig. S2. Evolution of the temperature anomaly and viscosity in case study 1.1 Fig. S3. Evolution of the temperature anomaly and viscosity in case study 1.2 7
8 Fig. S4. Evolution of the temperature anomaly and viscosity in case study 1.3 Change of the depth of the model domain To analyse a sensitivity of model results to the conditions imposed on the lower boundary of the model domain (no-slip conditions), we performed two case studies by increasing the depth of the model domain to 1000 km (case study 2.1) and 2000 km (case study 2.2). Figure S5 compares temperature anomalies and viscosities in case study 1.2 with those in the case studies 2.1 and 2.2. We see that the results are rather robust to the depth of the domain. And again there exists a link between the hotter rocks in sub-slab mantle with those in the mantle wedge (see areas marked by green ovals in Fig. S5). Influence of phase changes To analyse the effect of phase changes on the data assimilation, three case studies have been performed. Figures S-S8 present the results (temperature anomalies and viscosities) of the case studies. The upwellings connecting the sub-slab mantle with the mantle wedge are observed in the models despite changing conditions at the interfaces of 410 km and 0 km depths (see areas marked by green ovals in the figures). 8
9 Fig. S5. Evolution of the temperature anomaly and viscosity for different depths of the model domain 9
10 Fig. S. Case study 3.1: As case study 1.2, but with no phase change introduced Fig. S7. Case study 3.2: As case study 1.2, but only with the 410-km phase change interface 10
11 Fig. S8. Case study 3.3: As case study 1.2, but only with the 0-km phase change interface COMPARISON OF SEISMIC TOMOGRAPHY IMAGES Figure S9 illustrates three images of P-wave seismic tomography models. The resolution of the models is as follows: WEPP2 5 : 128 (latitude) 25 (longitude) 32 (depth), GAP_P2 15 : 288 (latitude) 57 (longitude) 29 (depth), and MIT 58 : 25 (latitude) 512 (longitude) 4 (depth). The image of the lithospheric slab and the sub-slab hot anomaly in the MIT model are less pronounced although its resolution power is comparable with the resolution power of the GAP_P2. If the MIT model is considered, the fading slab image in the northern part of the studied region (see cross section CS3 in Fig. S9c) may be much favourable for the potential link between the sub-slab mantle and the mantle wedge, than the tomography models used in the study. 11
12 Fig. S9. Three seismic tomography images: WEPP2 5, MIT 58, and GAP_P
13 Movie Linkage between the hot sub-slab mantle and mantle wedge The movie (Linkage_BandC.mp4) presents the evolution of two iso-surfaces of hot 4% anomaly (red) and cold 4% anomaly (blue). It illustrates the linkage between the hot anomaly in the subslab (region c in Fig. 3) and the hot anomaly in the mantle wedge (region B in Fig. 3) in 38.9 Ma. The movie s show starts at 38.9 Ma and finishes at present. The movie can be downloaded from References 1. Christensen, U. R. & Yuen, D. A. Layered convection induced by phase transitions. J. Geophys. Res. 90, (1985). 2. McKenzie, D. & Bickle, M. J. The volume and composition of melt generated by extension of the lithosphere, J. Petrology 29, (1988). 3. Kawamoto, T. Hydrous phase stability and partial melt chemistry in H 2 O-saturated KLB-1 peridotite up to the uppermost lower mantle conditions. Phys. Earth Planet. Int., , (2004). 4. Schellart, W. P., Stegman, D. R. & Freeman, J. Global trench migration velocities and slab migration induced upper mantle volume fluxes: Constraints to find an Earth reference frame based on minimizing viscous dissipation. Earth-Science Reviews 88, (2008). 5. Samarskii, A. A. Theory of Finite Difference Schemes (in Russian). Nauka, Moscow (1977).. Steinberger, B. & O Connell, R. J. Advection of plumes in mantle flow: implications for hotspot motion, mantle viscosity and plume distribution. Geophys. J. Int. 132, (1998). 7. Ismail-Zadeh, A. T., Talbot, C. J. & Volozh, Y. A. Dynamic restoration of profiles across diapiric salt structures: numerical approach and its applications. Tectonophysics 337, 21-3 (2001). 8. Conrad, C. P. & Gurnis, M. Seismic tomography, surface uplift, and the breakup of Gondwanaland: Integrating mantle convection backwards in time. Geochem. Geophys. Geosys. 4, 1031 (2003). 9. Ismail-Zadeh, A. T., Tsepelev, I. A., Talbot, C. J. & Volozh. Y. A. Three-dimensional forward and backward modelling of diapirism: Numerical approach and its applicability to the evolution of salt structures in the Pricaspian basin. Tectonophysics 387, (2004). 13
14 70. Bunge, H.-P., Richards, M. A. & Baumgardner, J. R. Mantle circulation models with sequential data-assimilation: Inferring present-day mantle structure from plate motion histories. Phil. Trans. Royal Soc. A 30, (2002). 71. Bunge, H.-P., Hagelberg, C. R. & Travis, B. J. Mantle circulation models with variational data assimilation: Inferring past mantle flow and structure from plate motion histories and seismic tomography. Geophys. J. Int. 152, (2003). 72. Ismail-Zadeh, A. T., Korotkii, A. I. & Tsepelev, I. A. Numerical approach to solving problems of slow viscous flow backwards in time. In: Computational Fluid and Solid Mechanics, Bathe, K. J. (ed.), Elsevier Science, Amsterdam, pp (2003). 73. Ismail-Zadeh, A., Schubert, G., Tsepelev, I. & Korotkii, A. Inverse problem of thermal convection: Numerical approach and application to mantle plume restoration. Phys. Earth Planet. Inter. 145, (2004). 74. Ismail-Zadeh, A., Schubert, G., Tsepelev, I. & Korotkii, A. Three-dimensional forward and backward numerical modeling of mantle plume evolution: Effects of thermal diffusion. J. Geophys. Res. 111, B0401 (200). 75. Yu, N., Imatani, S. & Inoue, T. Characteristics of temperature field due to pulsed heat input calculated by non-fourier heat conduction hypothesis. JSME Int. J., Series A 47, (2004). 7. Patankar, S. V. Numerical Heat Transfer and Fluid Flow. McGraw-Hill, New York (1980). 77. Samarskii, A. A. & Vabischevich, P.N. Computational Heat Transfer. Vol. 1. Mathematical Modeling. John Wiley & Sons, New York (1995). 78. Van der Vorst, H. A. Bi-CGSTAB: A fast and smoothly converging variant of Bi-CG for the solution of nonsymmetric linear systems. SIAM J. Sci. Stat. Comput. 13, (1992). 79. Saad, Y. Iterative Methods for Sparse Linear Systems, PWS, Boston (199). 14
15 TABLES Table S1. Model parameters Parameter Symbol Value Dimensionless density jump at the 410-km phase boundary a Dimensionless density jump at the 0-km phase boundary a Thermal conductivity c 1250 W m -1 K -1 Activation energy E a J mol -1 Acceleration due to gravity g m s Depth h 800 km Length (in x-direction) l km Length (in y-direction) l km Universal gas constant R J mol -1 K -1 Difference between the temperatures at the lower (T l ) T * 1594 K and upper (T u ) model boundaries Dimensionless temperature at the upper model boundary T u 290 / T * Dimensionless temperature at the lower model boundary T l 1884 / T * Dimensionless temperature at the 410-km phase boundary T / T * Dimensionless temperature at the 0-km phase boundary T / T * Activation volume V a 4 10 m 3 mol -1 Dimensionless width of the 410- km phase transition w 1 10 km / h Dimensionless width of the 0- km phase transition w 2 10 km / h Dimensionless depth of the 410- km phase boundary z km / h Dimensionless depth of the 0- km phase boundary z km / h Thermal expansivity K -1 QRV regularization parameter Dimensionless Clapeyron (pressure-temperature) slope at the 410-km phase boundary Dimensionless Clapeyron slope at the 0-km phase boundary 1 2 Pa K -1 * * T ( gh) 2 10 Pa K -1 * * 1 T ( gh) 15
16 Reference viscosity * Pa s Thermal diffusivity m s Reference density * 3400 kg m -3 Phase regularization parameter Table S2. Conditions (on velocity and pressure) at the lateral sides of the model boundary Lateral sides u P 0 0 u P 0 0 Lower surface u P 0 0 u 0 P 0 u P 0 P 0 u 0 0 P u u* 0 0 u 0 P 0 Table S3. A list of case studies No. case study Activation volume, V a, m 3 mol -1 Model depth, km Phase change introduced km and 0 km km and 0 km km and 0 km km and 0 km km and 0 km No phase changes km km 1
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 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 informationDynamic Subsidence and Uplift of the Colorado Plateau. Supplementary Material
GSA DATA REPOSITORY 2010177 Liu and Gurnis Dynamic Subsidence and Uplift of the Colorado Plateau Supplementary Material Lijun Liu and Michael Gurnis Seismological Laboratory California Institute of Technology
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 informationSupplementary Material for. Mantle induced subsidence and compression in SE Asia
Supplementary Material for Mantle induced subsidence and compression in SE Asia Ting Yang 1 *, Michael Gurnis 1, Sabin Zahirovic 2 1 Seismological Laboratory, California Institute of Technology, Pasadena,
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 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 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 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 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 informationWhole Mantle Convection
Whole Mantle Convection Overview 1. Evidence for whole mantle convection 2. Model of whole mantle convection reconciling geophysical and geochemical data Transition Zone Water Filter Model 3. Evidence
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 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 informationStudies of Arc Volcanism and Mantle Behavior in Subduction Zones
1 Studies of Arc Volcanism and Mantle Behavior in Subduction Zones Using COMSOL Lee, Changyeol Faculty of fearth and denvironmental lsciences Chonnam National University, Gwangju, Republic of Korea Sponsor:
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 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 informationES265 Order of Magnitude Phys & Chem Convection
ES265 Order of Magnitude Phys & Chem Convection Convection deals with moving fluids in which there are spatial variations in temperature or chemical concentration. In forced convection, these variations
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 informationChapter 3 Constraining Dynamic Properties of Mantle 2
Chapter 3 Constraining Dynamic Properties of Mantle 2 28 3.1 Need for Assimilation of Time-Dependent Data in Real Problems Unlike atmospheric or oceanographic problems where many parameters within the
More informationNumerical 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 informationEarth as a planet: Interior and Surface layers
Earth as a planet: Interior and Surface layers Bibliographic material: Langmuir & Broecker (2012) How to build a habitable planet Internal structure of the Earth: Observational techniques Seismology Analysis
More informationEffects of geometry on the convection with core-cooling
Earth Planets Space, 50, 387 395, 1998 Effects of geometry on the convection with core-cooling Yasuyuki Iwase and Satoru Honda Department of Earth and Planetary Systems Science, Faculty of Science, Hiroshima
More informationOn the Purpose of Toroidal Motion in a Convecting Mantle
GEOPHYSICAL RESEARCH LETTERS, VOL. 22, PAGES 3107-3110, 1995 On the Purpose of Toroidal Motion in a Convecting Mantle David Bercovici Department of Geology Geophysics, SOEST, University of Hawaii, Honolulu
More informationHigh performance mantle convection modeling
High performance mantle convection modeling Jens Weismüller, Björn Gmeiner, Siavash Ghelichkhan, Markus Huber, Lorenz John, Barbara Wohlmuth, Ulrich Rüde and Hans-Peter Bunge Geophysics Department of Earth-
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 informationExercise 1: Vertical structure of the lower troposphere
EARTH SCIENCES SCIENTIFIC BACKGROUND ASSESSMENT Exercise 1: Vertical structure of the lower troposphere In this exercise we will describe the vertical thermal structure of the Earth atmosphere in its lower
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 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 informationControls 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 informationGeotherms. Reading: Fowler Ch 7. Equilibrium geotherms One layer model
Geotherms Reading: Fowler Ch 7 Equilibrium geotherms One layer model (a) Standard model: k = 2.5 W m -1 C -1 A = 1.25 x 10-6 W m -3 Q moho = 21 x 10-3 W m -2 shallow T-gradient: 30 C km -1 deep T-gradient:
More informationInverse problem of thermal convection: numerical approach and application to mantle plume restoration
Physics of the Earth and Planetary Interiors 145 (2004) 99 114 Inverse problem of thermal convection: numerical approach and application to mantle plume restoration Alik Ismail-Zadeh a,b,c,, Gerald Schubert
More informationThree-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 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 informationAuxiliary Material. Subduction of oceanic asthenosphere: evidence from sub-slab seismic anisotropy. Teh-Ru Alex Song, Hitoshi Kawakatsu
Auxiliary Material Subduction of oceanic asthenosphere: evidence from sub-slab seismic anisotropy Teh-Ru Alex Song, Hitoshi Kawakatsu correspondence to: tehrusong@gmail.com 1 1. The effect of anisotropy
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 informationStructure of the Earth and the Origin of Magmas
Page 1 of 12 EENS 2120 Petrology Tulane University Prof. Stephen A. Nelson Structure of the Earth and the Origin of Magmas This document last updated on 23-Jan-2015 Magmas do not form everywhere beneath
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 informationSimultaneous inversion of mantle properties and initial conditions using an adjoint of mantle convection
Click Here for Full Article JOURNAL OF GEOPHYSICAL RESEARCH, VOL. 113,, doi:10.1029/2008jb005594, 2008 Simultaneous inversion of mantle properties and initial conditions using an adjoint of mantle convection
More informationPLATE TECTONIC PROCESSES
Lab 9 Name Sec PLATE TECTONIC PROCESSES 1. Fill in the blank spaces on the chart with the correct answers. Refer to figures 2.3, 2.4 p.33 (2.2 and 2.3 on p. 23) as needed. 2. With your knowledge of different
More informationConvection under a lid of finite conductivity: Heat flux scaling and application to continents
JOURNAL OF GEOPHYSICAL RESEARCH, VOL. 112,, doi:10.1029/2005jb004192, 2007 Convection under a lid of finite conductivity: Heat flux scaling and application to continents C. Grigné, 1 S. Labrosse, 2 and
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 informationSubduction zones are complex plate boundaries in which variable geometry and structure can be
1 Chapter 1 Introduction Subduction zones are complex plate boundaries in which variable geometry and structure can be seismically observed. The along-strike transition from flat to normal subduction is
More informationThe solution of the discretized incompressible Navier-Stokes equations with iterative methods
The solution of the discretized incompressible Navier-Stokes equations with iterative methods Report 93-54 C. Vuik Technische Universiteit Delft Delft University of Technology Faculteit der Technische
More informationRecycling 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 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 informationThe 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 informationOn the origin and significance of subadiabatic temperature gradients in the mantle
Click Here for Full Article JOURNAL OF GEOPHYSICAL RESEARCH, VOL. 112,, doi:10.1029/2006jb004850, 2007 On the origin and significance of subadiabatic temperature gradients in the mantle G. Sinha 1 and
More informationPhysics of the Earth and Planetary Interiors
Physics of the Earth and Planetary Interiors 173 (2009) 354 364 Contents lists available at ScienceDirect Physics of the Earth and Planetary Interiors journal homepage: www.elsevier.com/locate/pepi The
More informationEdge Driven Convection and Iceland
Edge Driven Convection and Iceland Scott D. King Department of Earth and Atmospheric Sciences Purdue University, West Lafayette, Indiana One of the alternative hypotheses for hotspot volcanism is Edge-Driven
More informationChapter 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 informationEffect of an outer-rise earthquake on seismic cycle of large interplate earthquakes estimated from an instability model based on friction mechanics
Effect of an outer-rise earthquake on seismic cycle of large interplate earthquakes estimated from an instability model based on friction mechanics Naoyuki Kato (1) and Tomowo Hirasawa (2) (1) Geological
More informationG 3. AN ELECTRONIC JOURNAL OF THE EARTH SCIENCES Published by AGU and the Geochemical Society
Geosystems G 3 AN ELECTRONIC JOURNAL OF THE EARTH SCIENCES Published by AGU and the Geochemical Society Article Volume 4, Number 3 29 March 2003 1031, doi:10.1029/2001gc000299 ISSN: 1525-2027 Seismic tomography,
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 information1 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 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 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 information1. The process by which the ocean floor sinks beneath a deep-ocean trench and back into the mantle is known as
Sample questions 1. The process by which the ocean floor sinks beneath a deep-ocean trench and back into the mantle is known as b. continental drift c. subduction d. conduction 2. The transfer of thermal
More informationThermal Analysis Contents - 1
Thermal Analysis Contents - 1 TABLE OF CONTENTS 1 THERMAL ANALYSIS 1.1 Introduction... 1-1 1.2 Mathematical Model Description... 1-3 1.2.1 Conventions and Definitions... 1-3 1.2.2 Conduction... 1-4 1.2.2.1
More informationChapter 1. Governing Equations of GFD. 1.1 Mass continuity
Chapter 1 Governing Equations of GFD The fluid dynamical governing equations consist of an equation for mass continuity, one for the momentum budget, and one or more additional equations to account for
More informationA PRECONDITIONER FOR THE HELMHOLTZ EQUATION WITH PERFECTLY MATCHED LAYER
European Conference on Computational Fluid Dynamics ECCOMAS CFD 2006 P. Wesseling, E. Oñate and J. Périaux (Eds) c TU Delft, The Netherlands, 2006 A PRECONDITIONER FOR THE HELMHOLTZ EQUATION WITH PERFECTLY
More informationContinental Margin Geology of Korea : Review and constraints on the opening of the East Sea (Japan Sea)
Continental Margin Geology of Korea : Review and constraints on the opening of the East Sea (Japan Sea) Han-Joon Kim Marine Satellite & Observation Tech. Korea Ocean Research and Development Institute
More informationCOMPOSITION and PHYSICAL PROPERTIES GENERAL SUBJECTS. GEODESY and GRAVITY
COMPOSITION and PHYSICAL PROPERTIES Composition and structure of the continental crust Composition and structure of the core Composition and structure of the mantle Composition and structure of the oceanic
More informationPlume-induced topography and geoid anomalies and their implications for the Tharsis rise on Mars
JOURNAL OF GEOPHYSICAL RESEARCH, VOL. 109,, doi:10.1029/2003je002226, 2004 Plume-induced topography and geoid anomalies and their implications for the Tharsis rise on Mars James H. Roberts Department of
More informationTutorials 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 informationInfluence of continents and lithospheric plates on the shape of D layer and the spatial distribution of mantle plumes
RUSSIAN JOURNAL OF EARTH SCIENCES vol. 15, 3, 2015, doi: 10.2205/2015ES000552 Influence of continents and lithospheric plates on the shape of D layer and the spatial distribution of mantle plumes V. P.
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 informationLow plume excess temperature and high core heat flux inferred from non-adiabatic geotherms in internally heated mantle circulation models
Physics of the Earth and Planetary Interiors 153 (2005) 3 10 Low plume excess temperature and high core heat flux inferred from non-adiabatic geotherms in internally heated mantle circulation models Hans-Peter
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 information(b) What is the amplitude at the altitude of a satellite of 400 km?
Practice final quiz 2015 Geodynamics 2015 1 a) Complete the following table. parameter symbol units temperature T C or K thermal conductivity heat capacity density coefficient of thermal expansion volumetric)
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 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 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 informationAn Introduction of Aleutian Subduction Zone. Chuanmao Yang, Hong Yang, Meng Zhang, Wenzhong Wang 2016/04/29
An Introduction of Aleutian Subduction Zone Chuanmao Yang, Hong Yang, Meng Zhang, Wenzhong Wang 2016/04/29 Outline General Introduction Formation history Structure from seismic study geochemical features
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 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 informationDevelopment of a Predictive Simulation System for Crustal Activities in and around Japan - II
Development of a Predictive Simulation System for Crustal Activities in and around Japan - II Project Representative Mitsuhiro Matsu'ura Graduate School of Science, The University of Tokyo Authors Mitsuhiro
More informationConstitution of Magmas. Magmas. Gas Law. Composition. Atomic Structure of Magma. Structural Model. PV = nrt H 2 O + O -2 = 2(OH) -
Constitution of Magmas Magmas Best, Ch. 8 Hot molten rock T = 700-1200 degrees C Composed of ions or complexes Phase Homogeneous Separable part of the system With an interface Composition Most components
More informationContinent-sized anomalous zones with low seismic velocity at the base of Earth s mantle
SUPPLEMENTARY INFORMATION DOI: 10.1038/NGEO2733 Continent-sized anomalous zones with low seismic velocity at the base of Earth s mantle Edward J. Garnero 1, Allen K. McNamara 1, and Sang-Heon D. Shim 1
More informationRemote Sensing of the Earth s Interior
Remote Sensing of the Earth s Interior Earth s interior is largely inaccessible Origin and Layering of the Earth: Geochemical Perspectives Composition of Earth cannot be understood in isolation Sun and
More informationWhat Forces Drive Plate Tectonics?
What Forces Drive Plate Tectonics? The tectonic plates are moving, but with varying rates and directions. What hypotheses have been proposed to explain the plate motion? Convection Cells in the Mantle
More informationA unifying model for fluid flow and elastic solid deformation: a novel approach for fluid-structure interaction and wave propagation
A unifying model for fluid flow and elastic solid deformation: a novel approach for fluid-structure interaction and wave propagation S. Bordère a and J.-P. Caltagirone b a. CNRS, Univ. Bordeaux, ICMCB,
More informationShape of thermal plumes in a compressible mantle with depth-dependent viscosity
GEOPHYSICAL RESEARCH LETTERS, VOL. 39,, doi:10.1029/2012gl050959, 2012 Shape of thermal plumes in a compressible mantle with depth-dependent viscosity Wei Leng 1 and Michael Gurnis 1 Received 18 January
More informationGravitational constraints
Gravitational constraints Reading: Fowler p172 187 Gravity anomalies Free-air anomaly: g F = g g( λ ) + δg obs F Corrected for expected variations due to the spheroid elevation above the spheroid Bouguer
More informationForms of Energy. Energy: commonly defined as the capacity to do work (i.e. by system on its surroundings); comes in many forms
Forms of Energy Energy: commonly defined as the capacity to do work (i.e. by system on its surroundings); comes in many forms Work: defined as the product of a force (F) times times a displacement acting
More informationDifferential relations for fluid flow
Differential relations for fluid flow In this approach, we apply basic conservation laws to an infinitesimally small control volume. The differential approach provides point by point details of a flow
More information10. Paleomagnetism and Polar Wandering Curves.
Map of ocean floor Evidence in Support of the Theory of Plate Tectonics 10. Paleomagnetism and Polar Wandering Curves. The Earth's magnetic field behaves as if there were a bar magnet in the center of
More informationDocumentation of the Solutions to the SFPE Heat Transfer Verification Cases
Documentation of the Solutions to the SFPE Heat Transfer Verification Cases Prepared by a Task Group of the SFPE Standards Making Committee on Predicting the Thermal Performance of Fire Resistive Assemblies
More informationWhole Earth Structure and Plate Tectonics
Whole Earth Structure and Plate Tectonics Processes in Structural Geology & Tectonics Ben van der Pluijm WW Norton+Authors, unless noted otherwise 4/5/2017 14:45 We Discuss Whole Earth Structure and Plate
More informationDissolution and precipitation during flow in porous media
1/25 Class project for GEOS 692: Transport processes and physical properties of rocks Dissolution and precipitation during flow in porous media Gry Andrup-Henriksen Fall 2006 1 2/25 Outline Introduction
More informationSubduction zones 3 arc magmatism
5. 3 Subduction zones 3 arc magmatism Where can we observe magmatic/volcanic activities along subduction zones? Characteristics of arc magmatism (vs. mid-ocean ridge/intraplate magmatism) Structure of
More informationModels of convection-driven tectonic plates: a comparison of methods and results
Geophys. J. Int. (1992) 19, 481-487 Models of convection-driven tectonic plates: a comparison of methods and results Scott D. King, * Carl W. Gable2 and Stuart A. Weinstein3t IGPP, Scripps Institution
More informationConvection Induced by Cooling at One Side Wall in Two-Dimensional Non-Rotating Fluid Applicability to the Deep Pacific Circulation
Journal of Oceanography Vol. 52, pp. 617 to 632. 1996 Convection Induced by Cooling at One Side Wall in Two-Dimensional Non-Rotating Fluid Applicability to the Deep Pacific Circulation ICHIRO ISHIKAWA
More informationNumerical methods for the Navier- Stokes equations
Numerical methods for the Navier- Stokes equations Hans Petter Langtangen 1,2 1 Center for Biomedical Computing, Simula Research Laboratory 2 Department of Informatics, University of Oslo Dec 6, 2012 Note:
More informationSUPPLEMENTARY INFORMATION
SUPPLEMENTARY INFORMATION Miocene drainage reversal of the Amazon River driven by plate-mantle interaction. Shephard, G.E., Müller, R.D., Liu, L., Gurnis, M. Supplementary Figures SOM Figure 1. Non-dimensional
More informationSUPPLEMENTARY INFORMATION
SUPPLEMENTARY INFORMATION Seismic evidence for a global low velocity layer within the Earth s upper mantle SUPPLEMENTARY MATERIAL Benoît Tauzin 1, Eric Debayle 2 & Gérard Wittlinger 3 1 Department of Earth
More informationSummary and Conclusions
Chapter 9 Summary and Conclusions 9.1 Summary The contents of this thesis revolve around the question of what type of geodynamics was active in the Early Earth and other terrestrial planets. The geology
More informationRational derivation of the Boussinesq approximation
Rational derivation of the Boussinesq approximation Kiyoshi Maruyama Department of Earth and Ocean Sciences, National Defense Academy, Yokosuka, Kanagawa 239-8686, Japan February 22, 2019 Abstract This
More informationPlate motions, mantle convection and chaos
Ecosystems and Sustainable Development IX 43 Plate motions, mantle convection and chaos S. Quéré 1 & A. A. Berezin 2 1 University of Toronto, Toronto, Canada 2 McMaster University, Hamilton, Canada Abstract
More informationLecture Outline Friday January 12 Friday January 19, 2018
Lecture Outline Friday January 12 Friday January 19, 2018 Questions?? Key Points for today The structure of the Earth Introduction to Plate Tectonic Theory and Convection Chapter 2 Development of the Plate
More information