Timing of water plume eruptions on Enceladus explained by interior viscosity structure
|
|
- Octavia Long
- 5 years ago
- Views:
Transcription
1 SUPPLEMENTARY INFORMATION DOI: /NGEO2475 Timing of water plume eruptions on Enceladus explained by interior viscosity structure Marie Běhounková 1, Gabriel Tobie 2, Ondřej Čadek 1, Gaël Choblet 2, Carolyn Porco 3 & Francis Nimmo 4 1 Charles University in Prague, Faculty of Mathematics and Physics, Department of Geophysics, V Holešovičkách 2, Praha 8, Czech Republic. 2 Université de Nantes, CNRS, Laboratoire de Planétologie et Géodynamique de Nantes, UMR-6112, 2 rue de la Houssinière, Nantes Cedex. 3 CICLOPS, Space Science Institute, Boulder CO Dept. Earth and Planetary Sciences, University of California Santa Cruz, CA NATURE GEOSCIENCE 1
2 1 Interior model and viscosity structure Enceladus s mass M kg ice density ρ ice 925 kg m 3 water density ρ water 1000 kg m 3 core density ρ core varying outer radius of the ice shell r t km ocean width ice shell thickness L km far zone L lithospheric thickness L 20 km reduced lithospheric thickness around south pole L SP km extent of L SP (characteristic width) δ 20 width of high-to-low viscosity transition σ 5 eccentricity e viscosity η Pa s shear modulus µ Pa Andrade parameter α 0.3 Andrade parameter ζ 1 time steps per period N step number of periods (length of simulation) N per 10 horizontal resolution (maximum degree) j max 60 vertical resolution (number of layers) n l 96 Table S1: List of reference parameters. 2
3 2 Tensile stress along the stripes and the stress representations Figure S1 displays a comparison between the NAS and FOF representations, for a given interior structure, and for three different values of the threshold parameter, σ t. This shows that the modeled stress curves are more sensitive to the threshold value for a 120 regional sea than for a global ocean, especially for the FOF representation. This is explained by the fact that the amplitude of the stress along the faults is smaller for a regional sea than for a global ocean. For the largest value of the threshold parameter, the initiation of the activity can be offset. However, the time (true anomaly) at which the maximum is reached remains identical irrespective of the stress representation and ocean extent. 3
4 a NAS representation, global ocean b FOF representation, global ocean c NAS representation, = 120 d FOF representation, = 120 Figure S1: Sensitivity to the threshold value σ t ; η = Pa s, L = 60 km, L SP = 2.5 km, L = 20 km, δ = 20. 4
5 3 Comparison with other approaches In the case of a purely elastic shell with a global ocean, the numerical solution has been compared with the approach by Nimmo et al. 1 and with the results of the SatStress code based on the method of Wahr et al. 2 (Figure S2 and S3). Figures S2a and S2c show the time evolution of stress tensor components at two specific locations 1,3. Note that the sign of component σ ϑϕ differs from the one in Nimmo et al. 1 and Smith-Konter and Pappalardo 3 due to different bases (east-positive longitude is used here). Figures S2b and S2d display the evolution of the normal stress. The evolution of all stress components and normal stresses agrees very well for all three approaches. The Love numbers h 2 and l 2 needed in Nimmo et al. 1 were obtained by SatStress and the two approaches are thus not independent. Figure S3 shows the average stress along the faults as a function of time for ice shell thickness of 60 km and 40 km. 5
6 a b c d Figure S2: Comparison of the stress tensors for different approaches corresponding to 60 km thick ice layer and h 2 = and l 2 = at Alexandria tiger stripe 70 S, 165 W with fault orientation 43 1 (a b) and at Damascus tiger stripe 80 S, 315 W with fault orientation 40 3 (c d). a b Figure S3: Comparison between analytical solution 1 (green) and numerical solution for elastic ice shell and global ocean (blue) and threshold stress σ t = 0 Pa; a analytical solution for h 2 = , l 2 = and numerical solution for ice shell thickness of 60 km, b analytical solution for h 2 = and l 2 = and numerical solution for ice shell thickness of 40 km. 6
7 4 Sensitivity tests We performed a systematic exploration of the model parameter space to determine which parameters primarily control the time evolution of the stress along the SPT faults. These sensitivity tests showed that the main controlling parameters are the ocean angular width,, the minimum viscosity, η, and the lithospheric thickness beneath the SPT, L SP, and to a lesser extent the total ice shell thickness, L. The lithospheric thickness away from the SPT plays a role only in the case of a global ocean. Other details on the viscosity structure, such as the width of the transition between the low/high viscosity regions and the maximum viscosity in conductive regions, are found to have a minor influence. The main results of these sensitivity tests are detailed below. Combined effect of the ocean extent and the extent of the low viscosity region Both the ocean width and the viscosity structure influence the response of the ice shell and the time evolution of the stress along the faults. In order to separate the interplay between the viscosity structure and the boundary conditions (ocean/ice: force equilibrium or silicate/ice: no-slip), we performed a first set of simulations for a purely elastic shell and various ocean widths (Figure S4a-b). These simulations showed that the elastic shell response with a hemispheric ocean ( = 180 ) is similar to the global ocean case, and that the stress curves are delayed compared to the global ocean solution when the ocean width is further reduced. The delay reaches a maximum value of 2.5 hours for the smallest ocean width tested ( = 60 ), and depends slightly on the ice shell thickness. However, this delay is not sufficient to explain the observed time lag and an additional delay due to viscous response is required. For a global ocean, only the viscous response of the ice shell can explain the delay. As illustrated on Figure S4c-d, the delay is proportional to the extent and thickness of the viscosity region. Maximum delay of 4 hours is obtained for a global viscosity structure ( v = 360 ) and a thick ice shell (60 km). The delay is also more sensitive to the extent of the viscosity region for L = 60 km than for L = 40 km. For models with a regional sea, the low viscosity region can exist only above the ocean, as high temperature can be maintained only there. The combined effect between the ocean extent and the extent of the low viscosity zone makes the 7
8 Elastic interior a L = 40 km b L = 60 km Viscoelastic interior with a global ocean ( o = 360 ) and varying viscosity structure ( v = ) c L = 40 km d L = 60 km Viscoelastic interior with varying ocean and viscosity structure ( = o = v ) e L = 40 km f L = 60 km Figure S4: Separating the influence of the weak interior and the localized ocean, the NAS representation for the viscoelastic (η = Pa s, L SP = 2.5 km, L = 20 km, δ = 20 ) and the elastic interior; a b elastic ice shell, changing extent of the internal ocean o, note the overlapping curves for global ocean and o = 180 c d viscoelastic ice shell, viscosity structure defined by v, global ocean assumed ( o = 360 ), e f viscoelastic ice shell, viscosity structure and ocean defined by = o = v. 8
9 a NAS representation b FOF representation Figure S5: Sensitivity to the minimum viscosity η for both representations; = 120, L = 60 km, L SP = 2.5 km, L = 20 km, δ = 20. a NAS representation b FOF representation Figure S6: Sensitivity to the south pole lithospheric thickness L SP for both representations; η = Pa s, = 120, L = 60 km, L = 20 km, δ = 20. evolution of the delay less straightforward than in the case of a global ocean. The delay reaches a minimum value for = 180 (Figure S4e-f) due to the decreasing viscous response. For smaller ocean widths ( < 180 ) the delay again increases due to increasing influences of the restricted ocean boundary. For the combined ocean and viscosity effects, a maximum delay of 5.5 hours is obtained for = 60. For both global ocean and regional sea models, a delay greater or equal to 4 hours can be obtained if a low viscosity zone is considered. However, for global ocean models, a thicker and wider low viscosity region is required, which has very different implications for the thermal evolution of Enceladus. 9
10 a NAS representation, global ocean b FOF representation, global ocean c NAS representation, = 120 d FOF representation, = 120 Figure S7: Sensitivity to the extent of the south polar lithosphere δ; η = Pa s, L = 60 km and L = 20 km. Influence of the viscosity structure beneath the SPT The viscosity structure beneath the SPT is controlled by the minimum viscosity η, the thickness and extent of the thinned lithosphere, L SP and δ, and the total ice shell thickness L. As illustrated with Figure S5 and Figure S6 for a 120 regional sea, the delay increases with decreasing minimum viscosity and SPT lithospheric thickness. A delay 4 hours can be obtained only for a very thin SPT lithosphere (L SP < 5 km) and a low viscosity (< Pa s). We also tested how the extent of the thin lithosphere around the South Pole, defined with the width parameter δ, affects the viscoelastic solution (Figure S7). These tests showed that as long as the area with a minimal lithospheric thickness corresponds to a significant portion of the SPT (δ > 20 ), the solution remains unchanged. For a lower value (δ = 10 ), the delay is significantly reduced as the thickness on the external part of the SPT is larger. This effect is equivalent to an increase of the averaged thickness illustrated in Figure S6. This outcome is independent of the extent of the ocean (global ocean or = 120 ) and of the representation 10
11 a NAS representation, global ocean b FOF representation, global ocean c NAS representation, = 120 d FOF representation, = 120 Figure S8: Sensitivity to the ice shell thickness L for both representations; η = Pa s, L = 60 km, L SP = 2.5 km, L = 20 km, δ = 20. (NAS or FOF). Finally, we tested the influence of the total ice shell thickness (Figure S8). For a 120 regional sea, the influence is minor, confirming that the delay in this case is not very sensitive to the volume of dissipative ice. In contrast, for models with a global ocean, the delay is very sensitive to the thickness, which appears to be one of the main controlling parameters. In summary, global ocean models require a relatively thick ice shell (L 60 km) to produce significant time lag, while models with a regional sea can lead to significant time lag independently of the ice shell thickness, as long as the average viscosity and lithospheric thickness in the SPT are sufficiently small (η Pa s, L SP 5 km). Influence of the viscosity structure outside the SPT Besides the parameters defining the ocean and the viscosity structure beneath the SPT, i.e., δ, η, L and L SP, we tested the parameters defining the viscosity structure away from the SPT: the far zone lithospheric thickness 11
12 L, the maximum viscosity η max and the width of the transition between minimum/maximum viscosity σ. All these parameters appears to be of minor importance, except the far zone lithospheric thickness in the case of a global ocean. Far zone lithospheric thickness L Figure S9 depicts the dependence of the predicted curve on the far zone lithospheric thickness L. For models with a regional sea (S9c d), the influence is negligible as both the global viscosity structure and the lithospheric thickness below SPT remain identical. In the case of a global ocean (S9a b), the increase of the far zone lithospheric thickness leads to the stiffening of the entire ice shell and the predicted curve is thus advanced compared to curves predicted for low L. The case L = L (S9a b, magenta curve) is a limit of our parameterization corresponding to a case where convection occurs only below SPT the conduction prevails elsewhere. This example shows that the convection in the case of a global ocean cannot be limited only to the SPT in order to explain the observed shift in the plume activity. Maximum viscosity η max Figure S10 shows that predicted curves are not sensitive to the chosen value of maximum viscosity as long as η max Pa s. Width of the lateral transition σ between the low/high viscosity regions We also explored the influence of the width of the lateral transition σ between the low and the high viscosity regions (see structures in Figure S11). This effect is shown to be very limited unless the ocean is strongly localized ( = 60 ). Again, we interpret this last result as an effective increase of the lithospheric thickness beneath the SPT owing to a wider transition zone (σ = 20, Figure S12). 12
13 a NAS representation, global ocean b FOF representation, global ocean c NAS representation, = 120 d FOF representation, = 120 Figure S9: Sensitivity to the far zone lithospheric thickness L ; η = Pa s, L = 60 km, L SP = 2.5 km, δ =
14 a NAS representation, global ocean b FOF representation, global ocean c NAS representation, = 120 d FOF representation, = 120 Figure S10: Sensitivity to the maximum viscosity η max ; η = Pa s, L = 60 km, L SP = 2.5 km, δ =
15 = 120 = 60 σ = 5 σ = 20 σ = 5 σ = 20 Figure S11: Example of two different structures corresponding to localized ocean and different width of transition zone; L = 60 km, L SP = 2.5 km, η = Pa s and L = 20 km. a NAS representation b FOF representation Figure S12: Sensitivity to the width of the transition zone σ; η = Pa s, L = 60 km, L SP = 2.5 km, δ = 20 ; solid line σ = 5, dashed line σ =
16 5 Fitting procedure, successful models Tables S2, S3 and S4 show the overview successful models. When a zero threshold value is considered for stress (Table S3), the lowest misfit is observed for the NAS representation and the largest set of successful models is found for VIMS without 2005 and ISS datasets. Additionally, whereas for the NAS representation the possibility of a global ocean cannot be excluded, we found no succesfull model with a global ocean for the FOF representation. Threshold dependence (Table S4) shows that low threshold (< 1000 Pa) values are favored for the NAS representation. For the FOF representation higher threshold values (< Pa) produce lower misfits than for cases without any threshold especially for VIMS dataset. Note, however, that the misfit defining very good models is lower for the NAS representation rather than for the FOF representation, independent of the threshold considered. 16
17 a) VIMS data with year 2005 b) VIMS data without year 2005 c) ISS data Figure S13: Acceptable models for all datasets and FOF representation, red, orange and yellow denotes 5%, 10% and 20% best models. 17
18 dataset m ν C ν M F CM χ 2 νc χ 2 νm,max χ 2 νm,min NAS FOF VIMS VIMS wo ISS Table S2: Datasets and misfits; m number of data points; ν C and ν M degrees of freedom for constant model C and time dependent model M; F CM critical value of F-distribution at the confidence level of 95%; χ 2 νc misfit for constant model C, χ2 νm,max maximum possible value of misfit for acceptable models, χ 2 νm,min minimum misfit for a given representation. datasets VIMS/ISS VIMS wo 2005/ISS representation NAS FOF NAS FOF No. acceptable models No. successful models No. successful models (global ocean) No. successful models (localized ocean) Table S3: Number of acceptable and successful models for both representations. Both cases where year 2005 is included or removed from the VIMS data set are considered. Threshold value is zero. 18
19 No. No. No. No. succ. No. succ. σ t model accep. succ. global local 0 Pa 1000 Pa 2000 Pa 3000 Pa 4000 Pa 5000 Pa NAS FOF NAS FOF NAS FOF NAS FOF NAS FOF NAS FOF Table S4: Overview of the threshold dependence, acceptable and successful models for VIMS without year 2005 and ISS datasets and both representations. Whatever the threshold value, successful models are defined as χ 2 ν,cew < and χ2 ν,i/f < 9.12 for the NAS and χ2 ν,cew < and χ 2 ν,i/f < 9.39 the FOF representations, respectively. 19
20 6 Additional Figures 20
21 elastic case, global ocean, L = 60 km global ocean, L = 60 km, η = Pa s LSP = 2.5 km, L = 20 km, δ = 20 = 120, L = 60 km, η = Pa s, LSP = 2.5 km, L = 20 km, δ = 20 Figure S14: Example of internal structure (first column), maximum tensile stress along the tiger stripes (second column), mean anomaly for which maximum tensile stress is reached (third column). Ice shell thickness 60 km, η = Pa s, variable ocean width, L = 20 km, δ =
22 = 90, L = 60 km, η = Pa s, LSP = 2.5 km, L = 20 km, δ = 20 = 60, L = 60 km, η = Pa s, LSP = 2.5 km, L = 20 km, δ = 20 Figure S14: Continuation 22
23 VIMS without 2005 = 110, L = 60 km, η = Pa s, L SP = 2.5 km, L = 20 km, δ = 20 global ocean, L = 60 km, η = Pa s, L SP = 2.5 km, L = 20 km, δ = 20 Figure S15: Example of the most successful models with global and localized oceans for both representations. 23
24 VIMS with 2005 = 120, L = 80 km, η = Pa s, L SP = 2.5 km, L = 20 km, δ = 20 global ocean, L = 50 km, η = Pa s, L SP = 2.5 km, L = 20 km, δ = 20 Figure S15: Continuation 24
25 VIMS without 2005, FOF representation = 80, L = 50 km, η = Pa s, L SP = 2.5 km, L = 20 km, δ = 20 VIMS WITH 2005, FOF representation = 90, L = 50 km, η = Pa s, L SP = 2.5 km, L = 20 km, δ = 20 Figure S15: Continuation 25
26 VIMS without 2005, NAS representation L SP = 2.5 km, L = 20 km, δ = 20 L SP = 5 km, L = 20 km, δ = 20 L SP = 7.5 km, L = 20 km, δ = 20 Figure S16: The quality of the fit for the main parameters, different datasets and both representations. Successful models are shown as models with both halves of circle in red and/or orange, gray and white colors denote acceptable and non-acceptable models, respectively. 26
27 VIMS without 2005, NAS representation L SP = 2.5 km, L = 20 km, δ = 20 VIMS WITH 2005, NAS representation L SP = 2.5 km, L = 20 km, δ = 20 Figure S16: Continuation 27
28 VIMS without 2005, FOF representation L SP = 2.5 km, L = 20 km, δ = 20 VIMS WITH 2005, FOF representation L SP = 2.5 km, L = 20 km, δ = 20 Figure S16: Continuation 28
29 References [1] Nimmo, F., Spencer, J.R., Pappalardo, R.T. & Mullen, M.E. Shear heating as the origin of the plumes and heat flux on Enceladus. Nature 447, (2007). [2] Wahr, J. et al.. Modeling stresses on satellites due to nonsynchronous rotation and orbital eccentricity using gravitational potential theory, Icarus 200, (2008). [3] Smith-Konter, B. & Pappalardo, R.T. Tidally driven stress accumulation and shear failure of Enceladus s tiger stripes. Icarus 198, (2008). 29
Limits of Enceladus s ice shell thickness from tidally driven tiger stripe shear failure
GEOPHYSICAL RESEARCH LETTERS, VOL. 38,, doi:10.1029/2010gl044950, 2011 Limits of Enceladus s ice shell thickness from tidally driven tiger stripe shear failure John G. Olgin, 1 Bridget R. Smith Konter,
More informationThe effect of an asymmetric core on convection in Enceladus ice shell: Implications for south polar tectonics and heat flux
GEOPHYSICAL RESEARCH LETTERS, VOL. 40, 1 5, doi:10.1002/2013gl057149, 2013 The effect of an asymmetric core on convection in Enceladus ice shell: Implications for south polar tectonics and heat flux Adam
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 informationGeodynamics of icy satellites
Geodynamics of icy satellites Francis Nimmo (U. C. Santa Cruz) With help from: Bruce Bills, Eric Gaidos, Bernd Giese, Michael Manga, Isamu Matsuyama, Bill Moore, McCall Mullen, Bob Pappalardo, Louise Prockter,
More informationMobile lid convection beneath Enceladus south polar terrain
Click Here for Full Article JOURNAL OF GEOPHYSICAL RESEARCH, VOL. 113,, doi:10.1029/2008je003114, 2008 Mobile lid convection beneath Enceladus south polar terrain Amy C. Barr 1 Received 12 February 2008;
More informationarxiv: v1 [astro-ph.ep] 12 Dec 2016
Astronomy & Astrophysics manuscript no. manuscript_arxiv c ESO 2018 December 20, 2018 Frequency-dependent tidal dissipation in a viscoelastic Saturnian core and expansion of Mimas semi-major axis arxiv:1612.03664v1
More informationEART193 Planetary Capstone. Francis Nimmo
EART193 Planetary Capstone Francis Nimmo Last Time silicate melting How and why are melts generated? Increase in mantle potential temperature; or Reduction in solidus temperature (e.g. water); or Thinning
More informationPreviously Documented Basin-localized Extension on Mercury
GSA DATA REPOSITORY 2012323 Watters et al. Previously Documented Basin-localized Extension on Mercury The most widely distributed extensional landforms on Mercury seen in images from the Mariner 10 and
More informationStress equilibrium in southern California from Maxwell stress function models fit to both earthquake data and a quasi-static dynamic simulation
Stress equilibrium in southern California from Maxwell stress function models fit to both earthquake data and a quasi-static dynamic simulation Peter Bird Dept. of Earth, Planetary, and Space Sciences
More informationEART162: PLANETARY INTERIORS
EART162: PLANETARY INTERIORS Francis Nimmo Last Week Applications of fluid dynamics to geophysical problems Navier-Stokes equation describes fluid flow: Convection requires solving the coupled equations
More informationSUPPLEMENTARY INFORMATION
SUPPLEMENTARY INFORMATION DOI: 10.1038/NGEO1887 Diverse calving patterns linked to glacier geometry J. N. Bassis and S. Jacobs 1. Supplementary Figures (a) (b) (c) Supplementary Figure S1 Schematic of
More informationMimas and Enceladus: Formation and interior structure from astrometric reduction of Cassini images.
Mimas and Enceladus: Formation and interior structure from astrometric reduction of Cassini images. R. Tajeddine, V. Lainey, N. Rambaux, N.Cooper, S. Charnoz, C.D. Murray Center of Mass/Center of figure
More informationThermal Evolution and Overpressurization of Europa s Subsurface Ocean. Austin Green. November 4, Advisor: Laurent Montesi GEOL394H
Thermal Evolution and Overpressurization of Europa s Subsurface Ocean Austin Green November 4, 2012 Advisor: Laurent Montesi GEOL394H 1. Abstract Scientific consensus on the mechanics of formation of ridge
More informationThe thermal and orbital evolution of Enceladus: observational constraints and models
The thermal and orbital evolution of Enceladus: observational constraints and models Francis Nimmo UCSC Amy C. Barr PSI Marie Běhounková Charles University William B. McKinnon Washington University Enceladus
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 informationarxiv: v1 [astro-ph.ep] 29 Sep 2014
Astronomy & Astrophysics manuscript no. RMZL14_arxiv ESO 018 September 8, 018 The surface signature of the tidal dissipation of the core in a two-layer planet Research Note) F. Remus 1,,, S. Mathis,4,
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 information22. Enceladus: An Active Cryovolcanic Satellite
22. Enceladus: An Active Cryovolcanic Satellite John R. Spencer 1, Amy C. Barr 1, Larry W. Esposito 2, Paul Helfenstein 3, Andrew P. Ingersoll 4, Ralf Jaumann 5, Christopher P. McKay 6, Francis Nimmo 7,
More informationAppendix: Nomenclature
Appendix: Nomenclature 209 Appendix: Nomenclature Chapter 2: Roman letters A amplitude of surface wave m A j (complex) amplitude of tidal component [various] a right ascension rad or a magnitude of tidal
More informationSUPPLEMENTARY INFORMATION
SUPPLEMENTARY INFORMATION doi: 10.1038/ngeo739 Supplementary Information to variability and distributed deformation in the Marmara Sea fault system Tobias Hergert 1 and Oliver Heidbach 1,* 1 Geophysical
More informationIcarus. The orbital thermal evolution and global expansion of Ganymede
Icarus 200 (2009) 207 221 Contents lists available at ScienceDirect Icarus www.elsevier.com/locate/icarus The orbital thermal evolution and global expansion of Ganymede MichaelT.Bland a,,adamp.showman
More informationEmergence of Strong Nonlinear Viscoelastic Response of Semifluorinated Alkane Monolayers
Emergence of Strong Nonlinear Viscoelastic Response of Semifluorinated Alkane Monolayers Salomé Mielke 1, Taichi Habe 1,#, Mariam Veschgini 1, Xianhe Liu 2, Kenichi Yoshikawa 3, Marie Pierre Krafft 2,
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 informationThe global shape of Europa: Constraints on lateral shell thickness variations
Icarus 191 (2007) 183 192 www.elsevier.com/locate/icarus The global shape of Europa: Constraints on lateral shell thickness variations F. Nimmo a,, P.C. Thomas b, R.T. Pappalardo c, W.B. Moore d a Department
More informationI. Comparison between rheological data obtained by using the shear rate sweep and frequency
Electronic Supplementary Material (ESI) for Soft Matter. This journal is The Royal Society of Chemistry 2017 ELECTRONIC SUPPLEMENTARY INFORMATION (ESI) I. Comparison between rheological data obtained by
More informationSUPPLEMENTARY INFORMATION
SUPPLEMENTARY INFORMATION DOI: 10.1038/NGEO1639 Importance of density-compensated temperature change for deep North Atlantic Ocean heat uptake C. Mauritzen 1,2, A. Melsom 1, R. T. Sutton 3 1 Norwegian
More informationWhen 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 informationChapter 21 Enceladus: An Active Cryovolcanic Satellite
Chapter 21 Enceladus: An Active Cryovolcanic Satellite John R. Spencer, Amy C. Barr, Larry W. Esposito, Paul Helfenstein, Andrew P. Ingersoll, Ralf Jaumann, Christopher P. McKay, Francis Nimmo, and J.
More informationData Repository Hampel et al., page 1/5
GSA DATA REPOSITORY 2138 Data Repositor Hampel et al., page 1/5 SETUP OF THE FINITE-ELEMENT MODEL The finite-element models were created with the software ABAQUS and consist of a 1-km-thick lithosphere,
More informationTalks. December 2011.
Talks Oliveira, J.S., Langlais, B., Amit, H., Pais, M.A., Modeling the magnetic field of Mercury using the Time Dependent Equivalent Source Dipole method. European Geosciences Union, Vienna, Austria, April
More informationConvection in Enceladus ice shell: Conditions for initiation
Click Here for Full Article GEOPHYSICAL RESEARCH LETTERS, VOL. 34, L09202, doi:10.1029/2006gl028799, 2007 Convection in Enceladus ice shell: Conditions for initiation Amy C. Barr 1 and William B. McKinnon
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 informationCollisions and transport phenomena
Collisions and transport phenomena Collisions in partly and fully ionized plasmas Typical collision parameters Conductivity and transport coefficients Conductivity tensor Formation of the ionosphere and
More informationGeophysics Departmental Exam: 2004 Part 1
2004 Geophysics Departmental Exam: 2004 Part 1 This section is 90 minutes, closed book, and consists of questions designed to test your knowledge of facts and figures in the geosciences. The focus will
More informationStrain-dependent strength profiles Implication of planetary tectonics
Strain-dependent strength profiles Implication of planetary tectonics Laurent G.J. Montési 1 Frederic Gueydan 2, Jacques Précigout 3 1 University of Maryland 2 Université de Montpellier 2, 3 Université
More informationLOCAL STRESSES FROM THE GEOID
RESIF - Montpellier 30 January 2018 LOCAL STRESSES FROM THE GEOID Thierry CAMELBEECK Royal Observatory of Belgium, BE-1180 Uccle, Belgium thierry.camelbeeck@oma.be Stresses in Western Europe Classical
More informationAuthor's personal copy
Chapter 39 Francis Nimmo Dept. Earth and Planetary Sciences, University of California Santa Cruz, CA, USA Carolyn Porco CICLOPS, Space Science Institute, Boulder, CO, USA Chapter Outline 1. Introduction
More informationcentrifugal acceleration, whose magnitude is r cos, is zero at the poles and maximum at the equator. This distribution of the centrifugal acceleration
Lecture 10. Equations of Motion Centripetal Acceleration, Gravitation and Gravity The centripetal acceleration of a body located on the Earth's surface at a distance from the center is the force (per unit
More informationInformation on internal structure from shape, gravity field and rotation
Information on internal structure from shape, gravity field and rotation Seismological information is available only for the Earth and in limited amounts for the Moon. Various geodetic data put constraints
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 informationESS314. Basics of Geophysical Fluid Dynamics by John Booker and Gerard Roe. Conservation Laws
ESS314 Basics of Geophysical Fluid Dynamics by John Booker and Gerard Roe Conservation Laws The big differences between fluids and other forms of matter are that they are continuous and they deform internally
More informationNormal stress causes normal strain σ 22
Normal stress causes normal strain blue box = before yellow box = after x 2 = Eɛ 22 ɛ 22 = E x 3 x 1 force is acting on the x2 face force is acting in the x2 direction Why do I draw this with equal stresses
More informationPlate Tectonics and the cycling of Earth materials
Plate Tectonics and the cycling of Earth materials Plate tectonics drives the rock cycle: the movement of rocks (and the minerals that comprise them, and the chemical elements that comprise them) from
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 informationTectonic patterns on reoriented and despun planetary bodies
Icarus 195 2008) 459 473 www.elsevier.com/locate/icarus Tectonic patterns on reoriented and despun planetary bodies Isamu Matsuyama a,, Francis Nimmo b a Department of Terrestrial Magnetism, Carnegie Institution
More informationAccretion of Phobos and Deimos in an extended debris disc stirred by transient moons
Accretion of Phobos and Deimos in an extended debris disc stirred by transient moons Pascal Rosenblatt * Royal Observatory of Belgium, Av. circulaire 3, B-1180 Uccle, Belgium Sébastien Charnoz Institut
More informationarxiv: v1 [astro-ph.ep] 11 Mar 2019
Draft version March 13, 2019 Typeset using LATEX default style in AASTeX62 Tides between the TRAPPIST-1 planets Hamish C. F. C. Hay 1 and Isamu Matsuyama 1 1 Lunar and Planetary Laboratory, University
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 informationSATELLITES: ACTIVE WORLDS AND EXTREME ENVIRONMENTS. Jessica Bolda Chris Gonzalez Crystal Painter Natalie Innocenzi Tyler Vasquez.
SATELLITES: ACTIVE WORLDS AND EXTREME ENVIRONMENTS Jessica Bolda Chris Gonzalez Crystal Painter Natalie Innocenzi Tyler Vasquez. Areas of interest! How did the Satellites of the outer solar system form
More informationTHE ANELASTIC EQUILIBRIUM TIDE IN EXOPLANETARY SYSTEMS
SFA 0 S. Boissier, P. de Laverny, N. Nardetto, R. Samadi, D. Valls-Gabaud and H. Wozniak (eds THE ANELASTIC EQUILIBRIUM TIDE IN EXOPLANETARY SYSTEMS F. Remus,,, S. Mathis, J.-P. Zahn and V. Lainey Abstract.
More informationSUPPLEMENTARY INFORMATION
doi:10.1038/nature12400 Contents S1 Data 2 S1.1GravityModel... 2 S1.2TopographyModel... 2 S2 Ice Shell Model 6 S2.1 Origin of long-wavelength topography............... 6 S2.2Iceshellstructure... 6 S2.3Lateralflowintheiceshell...
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 informationIntroduction 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 informationLong-term Crustal Deformation in and around Japan, Simulated by a 3-D Plate Subduction Model
Long-term Crustal Deformation in and around Japan, Simulated by a 3-D Plate Subduction Model Chihiro Hashimoto (1) and Mitsuhiro Matsu ura (2) (1) Institute of Frontier Research for Earth Evolution, Japan
More informationPhysical Processes in Acoustic Wave Heating of the Thermosphere
Publications 4-9-2005 Physical Processes in Acoustic Wave Heating of the Thermosphere G. Schubert Institute of Geophysics and Planetary Physics, University of California Michael P. Hickey Ph.D. Embry-Riddle
More informationarxiv: v2 [astro-ph.ep] 10 Aug 2017
Noname manuscript No. will be inserted by the editor Hugo Folonier Sylvio Ferraz-Mello Tidal synchronization of an anelastic multi-layered body: Titan s synchronous rotation. August 11, 017 arxiv:1706.08603v
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 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 informationEnceladus. Michelle Wenz
Enceladus Michelle Wenz Enceladus General Information Diameter 314 miles (505 km) Fit in Arizona (Pozio 1995) 1/7 th of Earth s moon diameter Herschel discovered in 1789 During the equinox of Saturn (reduced
More information1. Introduction, fluid properties (1.1, 2.8, 4.1, and handouts)
1. Introduction, fluid properties (1.1, 2.8, 4.1, and handouts) Introduction, general information Course overview Fluids as a continuum Density Compressibility Viscosity Exercises: A1 Fluid mechanics Fluid
More informationCourse Business. Today: isostasy and Earth rheology, paper discussion
Course Business Today: isostasy and Earth rheology, paper discussion Next week: sea level and glacial isostatic adjustment Email did you get my email today? Class notes, website Your presentations: November
More informationMelting-induced stratification above the Earth s inner core due to convective translation
Melting-induced stratification above the Earth s inner core due to convective translation Thierry Alboussiere, Renaud Deguen, Mickael Melzani To cite this version: Thierry Alboussiere, Renaud Deguen, Mickael
More informationImplications of shear heating and fracture zones for ridge formation on Europa
Click Here for Full Article GEOPHYSICAL RESEARCH LETTERS, VOL. 35, L03202, doi:10.1029/2007gl031957, 2008 Implications of shear heating and fracture zones for ridge formation on Europa Lijie Han 1 and
More informationSupplementary Information. Text S1:
Supplementary Information Text S1: In order to characterize the change in visco-elastic response in the course of a shear thickening transition in a controlled shear stress flow, on a fresh sample of for
More informationarxiv: v1 [physics.geo-ph] 16 May 2012
Obliquity of the Galilean satellites: The influence of a global internal liquid layer arxiv:125.3628v1 [physics.geo-ph] 16 May 212 Abstract Rose-Marie Baland a,b, Marie Yseboodt b, Tim Van Hoolst b a Université
More informationHigh heat flow from Enceladus south polar region measured using cm 1 Cassini/CIRS data
JOURNAL OF GEOPHYSICAL RESEARCH, VOL. 116,, doi:10.1029/2010je003718, 2011 High heat flow from Enceladus south polar region measured using 10 600 cm 1 Cassini/CIRS data C. J. A. Howett, 1 J. R. Spencer,
More informationDoes the normal stress parallel to the sliding plane affect the friction of ice upon ice?
Journal of Glaciology, Vol. 52, No. 205, 2011 949 Does the normal stress parallel to the sliding plane affect the friction of ice upon ice? Andrew L. FORTT, Erland M. SCHULSON Ice Research Laboratory,
More informationSupplementary Figure 1 - Maximum tidal amplitude in French Polynesia. The global tide model
Supplementary Figure 1 - Maximum tidal amplitude in French Polynesia. The global tide model TPXO7.2 1 is calibrated with satellite observations and used to predict tidal amplitudes in open ocean locations.
More informationImportance of Solar System Objects discussed thus far. Interiors of Terrestrial Planets. The Terrestrial Planets
Importance of Solar System Objects discussed thus far Interiors of Terrestrial Planets Chapter 9 Sun: Major source of heat for the surfaces of planets Asteroids: Provide possible insight to the composition
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 informationLife in the outer Solar System. AST 309 part 2: Extraterrestrial Life
Life in the outer Solar System AST 309 part 2: Extraterrestrial Life Prospects for life on: Overview: 1. Europa (Jupiter moon) 2. Titan (Saturn s moon) 3. Enceladus (Saturn s moon) Life on Europa? Europa
More informationMaterials and Methods The deformation within the process zone of a propagating fault can be modeled using an elastic approximation.
Materials and Methods The deformation within the process zone of a propagating fault can be modeled using an elastic approximation. In the process zone, stress amplitudes are poorly determined and much
More informationEarth. Interior Crust Hydrosphere Atmosphere Magnetosphere Tides
Earth Interior Crust Hydrosphere Atmosphere Magnetosphere Tides Semi-major Axis 1 A.U. Inclination 0 Orbital period 1.000 tropical year Orbital eccentricity 0.017 Rotational period 23 h 56 min 4.1 s Tilt
More information12 Gravity and Topography
12. Gravity and Topography 126 12 Gravity and Topography In the limit of perfect hydrostatic equilibrium, there is an exact relationship between gravity and topography... and no new information emerges:
More information1. What is Wegener s theory of continental drift? 2. What were the 4 evidences supporting his theory? 3. Why wasn t Wegener s theory excepted?
Notebook 7a and 7b Objective (left-side): Students will be able to explain the theory of plate tectonics and be able to make predictions about plate interactions. Mapping Reflection 1. What is Wegener
More informationMimas, moon of Saturn and Death Star impersonator responsible for several gaps in Saturn s ring system
Last time: Gravitational signs of large outer moons in the rings Ring shepherding/gap maintenance Longer lived structures due to mean motion resonances with large satellites Example: 2:1 resonance with
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 informationLecture 8. Lecture 1. Wind-driven gyres. Ekman transport and Ekman pumping in a typical ocean basin. VEk
Lecture 8 Lecture 1 Wind-driven gyres Ekman transport and Ekman pumping in a typical ocean basin. VEk wek > 0 VEk wek < 0 VEk 1 8.1 Vorticity and circulation The vorticity of a parcel is a measure of its
More informationTightly linked zonal and meridional sea surface temperature gradients over the past five million years
SUPPLEMENTARY INFORMATION DOI: 10.1038/NGEO2577 Tightly linked zonal and meridional sea surface temperature gradients over the past five million years Alexey V. Fedorov 1*, Natalie J. Burls 1,4, Kira T.
More informationSupplementary material to On the rheology of pendular gels and morphological developments in paste- like ternary systems based on capillary attraction
Electronic Supplementary Material (ESI) for Soft Matter. This journal is The Royal Society of Chemistry 214 Supplementary material to On the rheology of pendular gels and morphological developments in
More informationSUPPLEMENTARY INFORMATION
In the format provided by the authors and unedited. SUPPLEMENTARY INFORMATION DOI: 10.1038/NGEO3053 1 2 Contribution of topographically-generated submesoscale turbulence to Southern Ocean overturning 3
More informationProbing planetary interiors by spacecraft orbital observations
Probing planetary interiors by spacecraft orbital observations Alexander Stark, Jürgen Oberst, Frank Preusker, Klaus Gwinner, Gregor Steinbrügge, Hauke Hussmann Funded by Deutsche Forschungsgemeinschaft
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 informationAPPENDIX. A.1. Sensitivity analysis of the non-dimensional equations for stretch growth
335 336 337 338 339 340 34 342 343 344 APPENDIX A.. Sensitivity analysis of the non-dimensional equations for stretch growth Our goal in this section of the appendix is to show how the solutions of the
More informationQuestion 9.1: A steel wire of length 4.7 m and cross-sectional area 3.0 10 5 m 2 stretches by the same amount as a copper wire of length 3.5 m and cross-sectional area of 4.0 10 5 m 2 under a given load.
More informationFundamentals 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 informationPractice Test 3. Name: Date: ID: A. Multiple Choice Identify the choice that best completes the statement or answers the question.
Name: Date: _ Practice Test 3 Multiple Choice Identify the choice that best completes the statement or answers the question. 1. A wheel rotates about a fixed axis with an initial angular velocity of 20
More informationTidal heating in Enceladus
Icarus 188 (2007) 535 539 Note Tidal heating in Enceladus wwwelseviercom/locate/icarus Jennifer Meyer, Jack Wisdom Massachusetts Institute of Technology, Cambridge, MA 02139, USA Received 20 February 2007;
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 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 informationHYBRID DECADE-MEAN GLOBAL SEA LEVEL WITH MESOSCALE RESOLUTION. University of Hawaii, Honolulu, Hawaii, U.S.A.
HYBRID DECADE-MEAN GLOBAL SEA LEVEL WITH MESOSCALE RESOLUTION Nikolai A. Maximenko 1 and Pearn P. Niiler 2 1 International Pacific Research Center, School of Ocean and Earth Science and Technology, University
More informationClass XI Physics. Ch. 9: Mechanical Properties of solids. NCERT Solutions
Downloaded from Class XI Physics Ch. 9: Mechanical Properties of solids NCERT Solutions Page 242 Question 9.1: A steel wire of length 4.7 m and cross-sectional area 3.0 10 5 m 2 stretches by the same amount
More informationDynamics of Saturn s South Polar Vortex
Dynamics of Saturn s South Polar Vortex Ulyana A. Dyudina 1, Andrew P. Ingersoll 1,Shawn P. Ewald 1, Ashwin R. Vasavada 2, Robert A. West 2, Anthony D. Del Genio 3, John M. Barbara 3, Carolyn C. Porco
More informationImage: solarsystem.nasa.gov. Orbital trajectory of Cassini spacecraft ( ).
1 Enceladus This image, acquired by the Cassini spacecraft, captures Saturn, its rings (edge on), and the moon Enceladus. It was discovered that this moon emits jets of ice from possible underground seas.
More informationAntarctic sea-ice expansion between 2000 and 2014 driven by tropical Pacific decadal climate variability
Antarctic sea-ice expansion between 2000 and 2014 driven by tropical Pacific decadal climate variability Gerald A. Meehl 1, Julie M. Arblaster 1,2, Cecilia M. Bitz 3, Christine T.Y. Chung 4, and Haiyan
More informationWhat 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 informationSupporting Information for The origin of water-vapor rings in tropical oceanic cold pools
GEOPHYSICAL RESEARCH LETTERS Supporting Information for The origin of water-vapor rings in tropical oceanic cold pools Wolfgang Langhans 1 and David M. Romps 1,2 Contents of this file 1. Texts S1 to S2
More informationSection Forces Within Earth. 8 th Grade Earth & Space Science - Class Notes
Section 19.1 - Forces Within Earth 8 th Grade Earth & Space Science - Class Notes Stress and Strain Stress - is the total force acting on crustal rocks per unit of area (cause) Strain deformation of materials
More informationEART162: PLANETARY INTERIORS
EART16: PLANETARY INTERIORS Francis Nimmo F.Nimmo EART16 Spring 10 Last Week Elasticity: G E (1 v) xx Flexural equation gives deflection w in response to load 4 d w D gw q( x) 4 m w dx The flexural parameter
More informationTEACHING PLATE TECTONICS FROM THE EVIDENCE. Stephen T. Allard Associate Professor of Geoscience Winona State University
TEACHING PLATE TECTONICS FROM THE EVIDENCE Stephen T. Allard Associate Professor of Geoscience Winona State University 1 TEACHING PLATE TECTONICS FROM THE EVIDENCE Stephen T. Allard Associate Professor
More information