GOCE-GDC: Towards a better understanding of the Earth s interior and geophysical exploration research Work Package 5 - Geophysics B
|
|
- Margaret Ryan
- 6 years ago
- Views:
Transcription
1 GOCE-GDC: Towards a better understanding of the Earth s interior and geophysical exploration research Work Package 5 - Geophysics B Zdeněk Martinec Dublin Institute for Advanced Studies 5 Merrion Square, Dublin 2, Ireland (i) GOCE solid Earth workshop (October 16-17, 212) (ii) Joint science meeting GOCE+ GeoExplore and GOCE+ GDC (October 18) (iii) GOCE+ GDC MTR (October 19) University of Twente, Enschede
2 Contents 1 Static-geoid interpretation 2 Dynamic-geoid interpretation 3 Seismic tomography of the upper mantle 4 Congo-basin negative free-air gravity anomaly 5 Mass-density Green s functions for gravity
3 What can we learn from the geoid on the solid Earth? Newton s law of gravity: ϱ(q) V (P) = G dv Q B d PQ The inverse gravimetric problem (IGP): Find the volume-mass distribution ϱ of the Earth, provided that the potential V is known outside the Earth. An inverse problem is called well-posed if each of the following three criteria is satisfied (Hadamard): a solution exists the solution is unique the solution is stable, i.e. it continuously depends on the given data (here: V )
4 Inverse gravimetric problem is ill-posed The IGP is an ill-posed inverse problem since violates all three of Hadamard s criteria. Existence: V has to be harmonic outside the Earth which may be violated by measurement erros Uniqueness: Only the harmonic part of the density function can uniquely be reconstructed, whereas the (in the sense of the L2-space) orthogonal, so-called anharmonic, part has the external potential and, therefore, does not leave any trace in gravitational measurements. This is the most serious problem of IGP. Stability: The density does not continuously depend on the gravitational potential.
5 Non-uniqueness Let ϱ 1 be a solution of the IGP and let B have a continuous normal n. Then a class of all solutions of the IGP in L 2 (B) can be represented as where ϱ = ϱ h h B = h n = B In other words, 2 h generates a null external gravitational field. The Green s third identity proves this: for P in E 3 \B it holds 2 [ h(q) h(q) 1 V (P) = G dv Q = G h(q) ( )] 1 ds Q = B d PQ B n d PQ n d PQ i.e., V = in E 3 \B and hence on B.
6 Minimum-norm solution The minimum-norm solution: B [δϱ(r, Ω)] 2 dv = min δϱ under the integral constraints for potential coefficients V jm = δϱ(r, Ω)r j Y jm (Ω)dV B The variational form using the Lagrange-multiplier method j 1 max j ] δ [δϱ(r, Ω)] 2 dv + α jm [V jm δϱ(r, Ω)r j Y jm (Ω)dV 2 B B = j= m= j where δ is the variation with respect to δϱ. harmonic density anomaly δϱ(r, Ω) = j max j= m= j j α jm r j Y jm (Ω)
7 The decomposition of L 2 space L 2 = H H { } H := r j Y jm (Ω), j =, 1,..., m = j,..., j { H := 2 h, h B = h } n = B H the set of gravitationally equivalent bodies a solution H the harmonic solution the minimum L 2 -norm solution
8 Instability Let σ = on B and let ϱ := in E 3 \B. The Laplace-Poisson equation for potential V : 2 V = 4πGϱ in E 3 Using the 3-D Fourier transform of the form f (k) = (2π) 3/2 E 3 f (r) e ik.r d 3 r The Laplace-Poisson equation in the Fourier domain: k 2 Ṽ (k) = 4πG ϱ(k) k = k Two possible views: Forward GP Ṽ (k) = A(k) ϱ(k) Inverse GP ϱ(k) = Ṽ (k)/a(k) The transfer function mapping the density ϱ onto the potential V : A(k) = 4πG k 2 For k, 1/A(k), and a solution of the IGP is unstable.
9 Do gradiometric data bring new information on mass density? Answer: NO if the gradiometric data are considered as the Stokes potential coefficients Spherical harmonic representation of V V (r, Ω) = GM a ( a ) j+1 j V jm Y jm (Ω) r j= m= j Stokes potential coefficients V jm = 4π 1 Ma j 2j + 1 B Gradiometric model (e.g. GOCO2) grad grad V (r, Ω) = GM j= a j r j ϱ(r, Ω )Y jm(ω )dv j m= j [ ] Vjm GOCE grad grad r j 1 Y jm (Ω)
10 Newton s integral: Green s function of geoid for a static mantle δv (r, Ω) = G a r =b Using the addition theorem for r > r, 1 L(r, ψ, r ) = 4π r j j= m= j δϱ(r, Ω ) Ω L(r, ψ, r ) r 2 dr dω 1 2j + 1 ( r and writing the spherical harmonic expansions: r ) j Y jm(ω )Y jm (Ω), δv (r, Ω) = jm δv jm (r)y jm (Ω), δϱ(r, Ω) = jm δϱ jm (r)y jm (Ω) Then for r = a: δv jm (a) = a r =b G sta. j (a, r )δϱ jm (r )dr where Green s function of a static geoid is Gj sta. (a, r ) = 4πGa 2j + 1 ( r a ) j+2
11 Green s functions of static geoid (j = 2, 4, 1) The smaller spatial scales in gravitational potential, the higher sensitivity to shallower small-scale density variations, or in other words, the lower sensitivity to deeper small-scale density variations.
12 Contents 1 Static-geoid interpretation 2 Dynamic-geoid interpretation 3 Seismic tomography of the upper mantle 4 Congo-basin negative free-air gravity anomaly 5 Mass-density Green s functions for gravity
13 The Earth s interior is closed to hydrostatic equilibrium Flattening of the Earth: hydrostatic: f = 1/299. observed : f = 1/ The Earth s gravity interior is slighlty perturbed from the hydrostatic equilibrium: ϱ = ϱ }{{} hydrostatic + δϱ }{{} non-hydrostatic Additionally: Seismic tomographic model of the compressional (P) and transverse (S) seismic elastic waves Conversion of seismic anomalies to density anomalies: p f... proportionality factor δϱ ϱ = p δv f v
14 The geoid response to a static mantle Induced geoid undulations hundreds of meters Earth s mantle is not static but dynamic
15 Dynamic mantle The Earth s mantle is slighlty perturbed from the hydrostatic equilibrium: where ϱ = ϱ + δϱ g = g + δg δϱ ϱ and δg g The gravitational force acting on a material partical: ϱg ϱ g }{{} + δϱ g }{{} + ϱ δg }{{} hydrostatic force buoyancy force self-gravitating force Driving forces from equilibrium = buoyancy + self-gravitation
16 Fluid dynamic model of the geoid V δϱ V t t < t > geoid surface + δϱ > δϱ < cmb V δϱ : potential induced by mantle density anomalies δϱ V t : potential induced by deformations of boundary topographies t Courtesy of N.Tosi
17 Stokes problem for present-day mantle convection Mass conservation for an incompressible mantle: divv = Momentum conservation with Pr : divτ + f = Constitutive relation for a Newtonian viscous medium: τ = pi + 2ηε ε = ( v + T v)/2 Boundary conditions at B = B top B cmb v n = τ n (n τ n)n =
18 Weak formulation Consider a suitable functional space and create the functional F (v, p, l 1, l 2 ) = η(ε : ε) dv + f v dv + p div v dv B B B + l 1 (n v) ds + l 2 (n v) ds B top B cmb Compute the variation of F and use of Green s identity δf = (div τ + f ) δv dv + [τ n (n τ n)n] δv ds B B + (n v) δl 1 ds + (n v) δl 2 ds B top B cmb At a stationary point of F, δf =, the differential and the variational problem are equivalent.
19 Spherical-harmonic, finite-element parameterization Expansion in scalar and vector spherical harmonics with respect to angular coordinates: N j p = p jm (r)y jm (Ω) v = N j j= m= j j= m= j [U jm (r)s ( 1) jm (Ω) + V jm(r)s (1) jm (Ω) + W jm(r)s () jm (Ω)] Expansion into piecewise linear-finite elements with respect to radial coordinate: K +1 U jm (r) = Ujmψ k k (r) where: k=1 ψ k (r) = r k+1 r r k+1 r k and ψ k+1 (r) = r r k+1 r k+1 r k
20 Green s function of geoid for a dynamic mantle Assumption: 1-D viscosity model η = η(r). The gravitational effect of dynamic mantle + dynamic topographies: where a δϱ(r, Ω ) δv (r, Ω) = G r =b Ω L(r, ψ, r ) r 2 dr dω + G σ i (Ω) = ϱ i t i (Ω) n i=1 i = 1,..., n σ i (Ω ) Ω L(r, ψ, r i ) dω and t i (Ω) is the i-th dynamic topography. Easy manipulations results in δv jm (a) = a r =b G sta. j (a, r )δϱ jm (r )dr + n i=1 4πGa 2j + 1 where t jm are the coefficients of topography t(ω). Formally, δv jm (a) = a r =b G dyn. j (a, r )δϱ jm (r )dr ( ri ) j+2 ϱi (t i ) jm a Green s function of a dynamic geoid G dyn. j (a, r ) can only be modelled numerically.
21 Input data for forward dynamic-mantle modelling S-wave seismic tomographic model of Becker & Boschi (22). Conversion of seismic to density anomalies: proportionality factor p f =.2 3-layer viscosity model: δϱ ϱ = p δv f v η LT = Pa s r ( ) VM1 : η UM = Pa s r ( ) VM2 : η UM = Pa s r ( ) η LM = 4 η UM r ( )
22 Green s functions of dynamic geoid (j = 2, 8, 163) Depth Amplitude
23 Green s functions of surface dynamic topography (j = 2, 8, 16, 3) TBD
24 Geoid GOCO2 (j=2 12) SMEAN model LB model Geoid [m] Geoid [m]
25 Surface dynamic topography SMEAN model LB model Surface dynamic topography [km] Surface dynamic topography [km]
26 CMB dynamic topography SMEAN model LB model CMB dynamic topography [km] CMB dynamic topography [km]
27 SMEAN dynamic-mantle prediction SMEAN (j=2 3) 4 Free air gravity anomalies 4 Surface dynamic topography 4 CMB dynamic topography mgal km km
28 Observed vs. predicted g-fields 4 Geoid Anomalies 4 GOCO2 (j=2 3) Free Air Gravity Anomalies 4 d 2 V/dr meter mgal E Geoid Anomalies SMEAN dynamic mantle prediction Free Air Gravity Anomalies 4 4 d 2 V/dr
29 SMEAN dynamic-mantle prediction: VM1 vs. VM2 4 Geoid Anomalies 4 Viscosity model VM1 Free Air Gravity Anomalies 4 d 2 V/dr meter mgal E Geoid Anomalies 4 Viscosity model VM2 Free Air Gravity Anomalies 4 d 2 V/dr
30 Contents 1 Static-geoid interpretation 2 Dynamic-geoid interpretation 3 Seismic tomography of the upper mantle 4 Congo-basin negative free-air gravity anomaly 5 Mass-density Green s functions for gravity
31 Tomography of crust & shallow UM (Lebedev & Hilst, 28) 4 depth: 7 km 4 2 km 4 36 km km 4 8 km 4 11 km Shear velocity anomaly (%)
32 Tomography of deep UM (Lebedev & Hilst, 28) 4 depth: 15 km 4 2 km 4 26 km km 4 41 km km Shear velocity anomaly (%)
33 Lithospheric thickness Shear velocity anomaly (%) Fishwick & Bastow, 211 Lebedev & Hilst, 28
34 Contents 1 Static-geoid interpretation 2 Dynamic-geoid interpretation 3 Seismic tomography of the upper mantle 4 Congo-basin negative free-air gravity anomaly 5 Mass-density Green s functions for gravity
35 Free-air gravity anomalies GOCO2 (Free Air Gravity Anomalies) 4 Degree: mgal
36 Butterworth filter Spectral amplitude Spatial weight Cut-off frequency = 7.5 deg Angular degree Spherical distance (degree)
37 Masked GOCO2 free-air gravity anomalies j max = Free air gravity anomaly (mgal)
38 Possible contributors to the (static) gravity field Surface topography (Is ϱ = 267 kg/m 3 correct for sedimentary basin?) Lower-density sediments 3-D density variations within the crust (various geological units) Compensation at the Moho discontinuity 1-D (3-D) density variations in the lithosphere Compensation at the lithosphere boundary Pull-down or push-up by mantle flow at the lithosphere boundary 3-D density variations due to mantle convection
39 ETOPO1 topography ETOPO1 topography (m)
40 Total sediment thickness of Congo basin Sediment thickness (m) Kadima et al., 211
41 Degree amplitude spectra FA Gravity Sediments Degree amplitude spectrum (mgal) Degree j Degree amplitude spectrum (meter) Degree j
42 A density contrast The external gravitation of a body can be approximated by a material surface with either surface mass density ϱ t(ω) Airy s type σ(ω) = h ϱ(ω) Pratt s type The coefficients of the gravity induced by Airy s type mass surface: g ϱ jm = (j + 1) ϱt jm T jm = Let R jm be the residual gravity coefficients: R jm = g obs jm + g ϱ jm 3 a ϱ mean t jm 2j + 1 The density contrast ϱ will be searched by minimization of R jm : R jm Rjm = min ϱ It results in (Martinec, 1994): jm jm ϱ = (j + 1)T jm gjm obs jm (j + 1)2 T jm Tjm
43 Degree correlation and density contrast 1. 3 Degree correlation.5. Confidence level=.95 Density contrast (kg/m 3 ) kg/m Degree j Degree j Conclusions: Free-air gravity and gravity induced by sediments are correlated (in sense of statistical significance) up to degree j = 13 The density of sediments is in average of 21 kg/m 3 smaller than the density of surrounding crustal material
44 Gravity signals, j max = 13 6 FA Gravity 6 6 Sediments Free air gravity anomalies (mgal) Free air gravity anomalies (mgal)
45 Residual free-air gravity, j max = FA S gravity anomalies (mgal) Kadima et al. (211) hypothesis: Uplift of the Moho surface crustal thickening fill in by heavier mantle material
46 The effect of topography, j max = 13 Bouguer gravity corrections Geological units Bouguer gravity correction (mgal) Kadima et al., 211
47 The crustal thickness over Africa Crustal thickness (km) Crustal thickness (km) Pasyanos & Nyblade, 27
48 The Moho-induced gravity, j max = ϱ Moho = 2 kg/m Free air gravity anomalies (mgal)
49 A possible scenario The negative free-air gravity could be explained by lower-density sediments (it agrees with Kadima et al., 211), except the narrow, NW-SE oriented, positive gravity anomaly. The positive Bouguer gravity correction could be balanced by the negative Moho-induced gravity. However, the narrow, NW-SE oriented, positive gravity anomaly is amplified, not compensated (it contradicts to Kadima et al., 211 hypothesis). The Congo high-seismic velocity lithospheric craton may have a large density (due to thermal reason) and induced a higher gravity could be compensated by the negative gravity induced by mantle flow (not shown here, but also within the frame of Buiter at al., 211) Open question: How to interpret the narrow, NW-SE oriented, positive gravity anomaly?
50 Contents 1 Static-geoid interpretation 2 Dynamic-geoid interpretation 3 Seismic tomography of the upper mantle 4 Congo-basin negative free-air gravity anomaly 5 Mass-density Green s functions for gravity
51 Gravitational potential and its gradients Green s functions: G( r, r ) := V ( r) = κ B B ϱ( r ) G( r, r ) dv grad V ( r) = κ ϱ( r ) G( r, r ) dv grad grad V ( r) = κ ϱ( r ) G( r, r ) dv 1 L(r, ψ, r ) G( r, r ) := grad 1 L = r r L 3 G( r, r ) := grad grad 1 L = 1 [ I 3( r r ) ( r r ) ] L 3 L 2 B
52 An example in 2-D Cartesian geometry Observer x, x L ϱ(x, z ) Green s functions at (x =, z = ): z, z G xx = G xz = G zz = G x = G z = G = 1 L ( grad 1 ) L x ) ( grad 1 L ( grad grad 1 ) L ) ( grad grad 1 L ( grad grad 1 L ) z xx xz zz = = x L 3 = z L 3 1 x 2 + z 2 = 1 L 3 ( 1 = 3x z L 5 3x 2 L 2 ) = 1 L 3 ( 1 3z 2 L 2 )
53 Scalar, vector and tensor Green s functions for mass density 2 G zz z =const. Amplitude 1 G z G x G xz G x 1 G xx Conclusion: The higher the order of the derivative of the gravitational potential as a boundary datum, the smaller the contribution from far-zone mass-density anomalies.
54 Mass-density Green s functions for gradiometric data After lengthy derivation: in spherical geometry grad grad 1 L = 1 [K r 3 rr (t, x)e rr + K rω (t, x)(cos α e rϑ sin α e rϕ) ( ) +K ΩΩ (t, x) cos 2α (e ϑϑ e ϕϕ) 2 sin 2α e ϑϕ 1 ] 2 Krr (t, x) (e ϑϑ + e ϕϕ) Three isotropic kernels: ψ... spherical distance, x = cos ψ α... azimuth t = r /r K rr (t, x) = (j + 1)(j + 2)t j P j (x) j= K rω (t, x) = 2 1 x 2 (j + 2)t j dp j (x) dx j= K ΩΩ (t, x) = 1 2 (1 x 2 ) t j d 2 P j (x) dx 2 j=
55 Closed spatial forms K rr (t, x) = 1 3(1 tx)2 + g3 g 5 K rω (t, x) = 3t(1 tx) 1 x 2 g 5 K ΩΩ (t, x) = 1 2 (1 x 2 ) 3t 2 g 5 where g = 1 + t 2 2tx
56 Isotropic gradiometric Green s functions 6 1 Amplitude K rr K rω K ΩΩ h = 25 km Amplitude K rr K ΩΩ K rω h = 1 km ψ (degree) ψ (degree) Amplitude K rr K ΩΩ h = 5 km 2 K rω ψ (degree)
57 Integral equations for a surface-mass density For a fixed source position (r = const): Ω ϱ(ω )K α(ω, Ω )dω = d α(ω) α = 1,..., GOCE The Nyström method discretizes the integral equation by a quadrature rule: n w j ϱ(ω j )K α(ω k, Ω j ) = d α(ω k ) j=1 The Galerkin method: ϱ(ω) = n ϱ j φ j (Ω) j=1 k = 1,..., n φ j (Ω) a set of n linearly indep. functions (e.g. tapers on a spherical cap) n j=1 ϱ j dω dω Ω φ k (Ω)K α(ω, Ω )φ j (Ω ) = d α(ω)φ k (Ω)dΩ A Ω Ω A A Overdetermined s. of linear algebraic equations for a surface-mass density k = 1,..., n A α ϱ = d α α = 1,..., GOCE
7. Post Glacial Rebound. Ge 163 4/16/14-
7. Post Glacial Rebound Ge 163 4/16/14- Outline Overview Order of magnitude estimate of mantle viscosity Essentials of fluid mechanics Viscosity Stokes Flow Biharmonic equation Half-space model Channel
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 informationCHAPTER X. Second Half Review 2017
CHAPTER X Second Half Review 217 Here is a quick overview of what we covered in the second half of the class. Remember that the final covers the whole course but there will naturally be a bias towards
More informationMean Vertical Gradient of Gravity
Mean Vertical Gradient of Gravity P. Vaníek, J. Janák Department of Geodesy and Geomatics Engineering, University of New Brunswick, P.O.Box 4400, Fredericton, New Brunswick, Canada, E3B 2M8 J. Huang Geodetic
More informationGetting started: CFD notation
PDE of p-th order Getting started: CFD notation f ( u,x, t, u x 1,..., u x n, u, 2 u x 1 x 2,..., p u p ) = 0 scalar unknowns u = u(x, t), x R n, t R, n = 1,2,3 vector unknowns v = v(x, t), v R m, m =
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 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 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 informationTotal gravitational field is sum of contributions from all masses.
Gravity force (acceleration) vs potential (energy) acceleration (g) => GM/r 2 Potential => - GM/r G is Newton s gravitational constant 6.67x10-11 (S.I. units) you should determine what the S.I. units are
More informationA NEW GLOBAL CRUSTAL MODEL BASED ON GOCE DATA GRIDS
A NEW GLOBAL CRUSTAL MODEL BASED ON GOCE DATA GRIDS M. Reguzzoni, D. Sampietro DIIAR Politecnico di Milano CRUSTAL MODELS SEISMIC MODEL SEISMIC AND GRAVITY Moho model from Meier et al. 2007; unit [km].
More informationEquations of linear stellar oscillations
Chapter 4 Equations of linear stellar oscillations In the present chapter the equations governing small oscillations around a spherical equilibrium state are derived. The general equations were presented
More informationLast week we obtained a general solution: 1 cos αdv
GRAVITY II Surface Gravity Anomalies Due to Buried Bodies Simple analytical solution may be derived for bodies with uniform density contrast simple shape, such as: Sphere Horizontal/vertical cylinders
More informationQuick Recapitulation of Fluid Mechanics
Quick Recapitulation of Fluid Mechanics Amey Joshi 07-Feb-018 1 Equations of ideal fluids onsider a volume element of a fluid of density ρ. If there are no sources or sinks in, the mass in it will change
More informationBoundary-Value Problems for Gravimetric Determination of a Precise Geoid
Zdenek Martinec Boundary-Value Problems for Gravimetric Determination of a Precise Geoid With 51 Figures and 3 Tables Springer List of symbols Acknowledgements ix xii Introduction 1 1 The Stokes two-boundary-value
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 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 informationRheology of the Mantle and Plates (part 1): Deformation mechanisms and flow rules of mantle minerals
(part 1): Deformation mechanisms and flow rules of mantle minerals What is rheology? Rheology is the physical property that characterizes deformation behavior of a material (solid, fluid, etc) solid mechanics
More informationIsostacy: Compensation of Topography and Isostatic Geoid Anomalies
Isostacy: Compensation of Topography and Isostatic Geoid Anomalies November 3, 2016 1 Isostacy One of the major goals of interpreting planetary gravity fields is to understand how topography positive (
More informationv It appears that under their own units the gravity gradient is a million times worse than the gravity. v In reality this is not the case.
v It appears that under their own units the gravity gradient is a million times worse than the gravity. v In reality this is not the case. v In particular, the opposite is true in terms of dealing with
More informationExercise: concepts from chapter 10
Reading:, Ch 10 1) The flow of magma with a viscosity as great as 10 10 Pa s, let alone that of rock with a viscosity of 10 20 Pa s, is difficult to comprehend because our common eperience is with s like
More informationNavier-Stokes Equation: Principle of Conservation of Momentum
Navier-tokes Equation: Principle of Conservation of Momentum R. hankar ubramanian Department of Chemical and Biomolecular Engineering Clarkson University Newton formulated the principle of conservation
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 informationGRAVITY AND ISOSTASY
GRAVITY AND ISOSTASY Gravity The Geoid is the oblate spheroid -- the sea level surface over the entire Earth's surface Physically, the Geoid is an equipotential surface for gravity (i.e. the surface for
More 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 informationConstraints on Shallow Low-Viscosity Earth Layers from Future GOCE Data
Constraints on Shallow Low-Viscosity Earth Layers from Future GOCE Data Hugo Schotman 1,2, Bert Vermeersen 2, Pieter Visser 2 1 2 3 rd International GOCE User Workshop, ESA Esrin, 7 November 2006 glacial-isostatic
More informationV (r,t) = i ˆ u( x, y,z,t) + ˆ j v( x, y,z,t) + k ˆ w( x, y, z,t)
IV. DIFFERENTIAL RELATIONS FOR A FLUID PARTICLE This chapter presents the development and application of the basic differential equations of fluid motion. Simplifications in the general equations and common
More informationBSL Transport Phenomena 2e Revised: Chapter 2 - Problem 2B.11 Page 1 of 5
BS Transport Phenomena 2e Revised: Chapter 2 - Problem 2B11 Page 1 of 5 Problem 2B11 The cone-and-plate viscometer (see Fig 2B11 A cone-and-plate viscometer consists of a flat plate and an inverted cone,
More informationGlobal surface-wave tomography
Global surface-wave tomography Lapo Boschi (lapo@erdw.ethz.ch) October 7, 2009 Love and Rayleigh waves, radial anisotropy Whenever an elastic medium is bounded by a free surface, coherent waves arise that
More informationChapter 2. General concepts. 2.1 The Navier-Stokes equations
Chapter 2 General concepts 2.1 The Navier-Stokes equations The Navier-Stokes equations model the fluid mechanics. This set of differential equations describes the motion of a fluid. In the present work
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 informationStress, Strain, and Viscosity. San Andreas Fault Palmdale
Stress, Strain, and Viscosity San Andreas Fault Palmdale Solids and Liquids Solid Behavior: Liquid Behavior: - elastic - fluid - rebound - no rebound - retain original shape - shape changes - small deformations
More information6.2 Governing Equations for Natural Convection
6. Governing Equations for Natural Convection 6..1 Generalized Governing Equations The governing equations for natural convection are special cases of the generalized governing equations that were discussed
More informationPEAT SEISMOLOGY Lecture 2: Continuum mechanics
PEAT8002 - SEISMOLOGY Lecture 2: Continuum mechanics Nick Rawlinson Research School of Earth Sciences Australian National University Strain Strain is the formal description of the change in shape of a
More informationFORWARD MODELING OF THE GEOID ANOMALY USING SPHERICAL HARMONICS: APPLICATIONS IN THE SIERRA NEVADA. Alissa C Scire
FORWARD MODELING OF THE GEOID ANOMALY USING SPHERICAL HARMONICS: APPLICATIONS IN THE SIERRA NEVADA by Alissa C Scire A Prepublication Manuscript Submitted to the Faculty of the DEPARTMENT OF GEOSCIENCES
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 informationThe Janus Face of Turbulent Pressure
CRC 963 Astrophysical Turbulence and Flow Instabilities with thanks to Christoph Federrath, Monash University Patrick Hennebelle, CEA/Saclay Alexei Kritsuk, UCSD yt-project.org Seminar über Astrophysik,
More informationTutorial School on Fluid Dynamics: Aspects of Turbulence Session I: Refresher Material Instructor: James Wallace
Tutorial School on Fluid Dynamics: Aspects of Turbulence Session I: Refresher Material Instructor: James Wallace Adapted from Publisher: John S. Wiley & Sons 2002 Center for Scientific Computation and
More informationChapter 5. The Differential Forms of the Fundamental Laws
Chapter 5 The Differential Forms of the Fundamental Laws 1 5.1 Introduction Two primary methods in deriving the differential forms of fundamental laws: Gauss s Theorem: Allows area integrals of the equations
More informationThe Earth is Red; Dominance of large wavelength features in mantle boundary layers. Adam M. Dziewonski
The Earth is Red; Dominance of large wavelength features in mantle boundary layers Adam M. Dziewonski KITP, EARTH08, July 2, 2008 Kircher, 1665 Imagining ca. Earth s 2005 Interior from E. Garnero s files
More informationCH.9. CONSTITUTIVE EQUATIONS IN FLUIDS. Multimedia Course on Continuum Mechanics
CH.9. CONSTITUTIVE EQUATIONS IN FLUIDS Multimedia Course on Continuum Mechanics Overview Introduction Fluid Mechanics What is a Fluid? Pressure and Pascal s Law Constitutive Equations in Fluids Fluid Models
More informationChapter 2: Fluid Dynamics Review
7 Chapter 2: Fluid Dynamics Review This chapter serves as a short review of basic fluid mechanics. We derive the relevant transport equations (or conservation equations), state Newton s viscosity law leading
More informationTides in Higher-Dimensional Newtonian Gravity
Tides in Higher-Dimensional Newtonian Gravity Philippe Landry Department of Physics University of Guelph 23 rd Midwest Relativity Meeting October 25, 2013 Tides: A Familiar Example Gravitational interactions
More informationContinuum mechanism: Stress and strain
Continuum mechanics deals with the relation between forces (stress, σ) and deformation (strain, ε), or deformation rate (strain rate, ε). Solid materials, rigid, usually deform elastically, that is the
More informationEKC314: TRANSPORT PHENOMENA Core Course for B.Eng.(Chemical Engineering) Semester II (2008/2009)
EKC314: TRANSPORT PHENOMENA Core Course for B.Eng.(Chemical Engineering) Semester II (2008/2009) Dr. Mohamad Hekarl Uzir-chhekarl@eng.usm.my School of Chemical Engineering Engineering Campus, Universiti
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 informationLecture 8: Tissue Mechanics
Computational Biology Group (CoBi), D-BSSE, ETHZ Lecture 8: Tissue Mechanics Prof Dagmar Iber, PhD DPhil MSc Computational Biology 2015/16 7. Mai 2016 2 / 57 Contents 1 Introduction to Elastic Materials
More informationPhysics 451/551 Theoretical Mechanics. G. A. Krafft Old Dominion University Jefferson Lab Lecture 18
Physics 451/551 Theoretical Mechanics G. A. Krafft Old Dominion University Jefferson Lab Lecture 18 Theoretical Mechanics Fall 018 Properties of Sound Sound Waves Requires medium for propagation Mainly
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 informationIntroduction to Fluid Mechanics
Introduction to Fluid Mechanics Tien-Tsan Shieh April 16, 2009 What is a Fluid? The key distinction between a fluid and a solid lies in the mode of resistance to change of shape. The fluid, unlike the
More information1. University of Ottawa 2. Dublin Institute for Advanced Studies 3. University of Texas at Austin
Supplementary Figures: Azimuthal seismic anisotropy in the Earth s upper mantle and the thickness of tectonic plates A.J. Schaeffer 1, S. Lebedev 2 and T.W. Becker 3 Geophysical Journal International July
More informationGlobal geophysics and wave propagation
Global geophysics and wave propagation Reading: Fowler p76 83 Remote sensing Geophysical methods Seismology Gravity and bathymetry Magnetics Heat flow Seismology: Directly samples the physical properties
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 informationGS388 Handout: Radial density distribution via the Adams-Williamson equation 1
GS388 Handout: Radial density distribution via the Adams-Williamson equation 1 TABLE OF CONTENTS ADIABATIC COMPRESSION: THE ADAMS WILLIAMSON EQUATION...1 EFFECT OF NON-ADIABATIC TEMPERATURE GRADIENT...3
More informationSolution Set Two. 1 Problem #1: Projectile Motion Cartesian Coordinates Polar Coordinates... 3
: Solution Set Two Northwestern University, Classical Mechanics Classical Mechanics, Third Ed.- Goldstein October 7, 2015 Contents 1 Problem #1: Projectile Motion. 2 1.1 Cartesian Coordinates....................................
More informationDynamics of the Mantle and Lithosphere ETH Zürich Continuum Mechanics in Geodynamics: Equation cheat sheet
Dynamics of the Mantle and Lithosphere ETH Zürich Continuum Mechanics in Geodynamics: Equation cheat sheet or all equations you will probably ever need Definitions 1. Coordinate system. x,y,z or x 1,x
More 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 information1. A few words about EarthScope and USArray. 3. Tomography using noise and Aki s method
1. A few words about EarthScope and USArray 2. Surface-wave studies of the crust and mantle 3. Tomography using noise and Aki s method 4. Remarkable images of US crust (and basins)! Unlocking the Secrets
More information(TRAVELLING) 1D WAVES. 1. Transversal & Longitudinal Waves
(TRAVELLING) 1D WAVES 1. Transversal & Longitudinal Waves Objectives After studying this chapter you should be able to: Derive 1D wave equation for transversal and longitudinal Relate propagation speed
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 informationChapter 1 Fluid Characteristics
Chapter 1 Fluid Characteristics 1.1 Introduction 1.1.1 Phases Solid increasing increasing spacing and intermolecular liquid latitude of cohesive Fluid gas (vapor) molecular force plasma motion 1.1.2 Fluidity
More informationESA training. Gravity, magnetics and gradients for mapping and modelling. Jörg Ebbing. Department of Geosciences Kiel University
ESA training Gravity, magnetics and gradients for mapping and modelling Jörg Ebbing Department of Geosciences Kiel University Contributions from: Eldar Baykiev (Trondheim), Des Fitzgerald (Melbourne),
More informationOPAC102. The Acoustic Wave Equation
OPAC102 The Acoustic Wave Equation Acoustic waves in fluid Acoustic waves constitute one kind of pressure fluctuation that can exist in a compressible fluid. The restoring forces responsible for propagating
More information1 POTENTIAL FLOW THEORY Formulation of the seakeeping problem
1 POTENTIAL FLOW THEORY Formulation of the seakeeping problem Objective of the Chapter: Formulation of the potential flow around the hull of a ship advancing and oscillationg in waves Results of the Chapter:
More informationIntroduction to the use of gravity measurements in Ge111A
Introduction to the use of gravity measurements in Ge111A Background & basic intuition Data reduction Use of the instrument See Reynolds for detailed quantitative discussion What and Why Gravity measures
More informationContinuum Mechanics. Continuum Mechanics and Constitutive Equations
Continuum Mechanics Continuum Mechanics and Constitutive Equations Continuum mechanics pertains to the description of mechanical behavior of materials under the assumption that the material is a uniform
More informationPhysics 451/551 Theoretical Mechanics. G. A. Krafft Old Dominion University Jefferson Lab Lecture 18
Physics 451/551 Theoretical Mechanics G. A. Krafft Old Dominion University Jefferson Lab Lecture 18 Theoretical Mechanics Fall 18 Properties of Sound Sound Waves Requires medium for propagation Mainly
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 informationContents. I Introduction 1. Preface. xiii
Contents Preface xiii I Introduction 1 1 Continuous matter 3 1.1 Molecules................................ 4 1.2 The continuum approximation.................... 6 1.3 Newtonian mechanics.........................
More informationAdvanced Newtonian gravity
Foundations of Newtonian gravity Solutions Motion of extended bodies, University of Guelph h treatment of Newtonian gravity, the book develops approximation methods to obtain weak-field solutions es the
More informationArtificial Intelligence & Neuro Cognitive Systems Fakultät für Informatik. Robot Dynamics. Dr.-Ing. John Nassour J.
Artificial Intelligence & Neuro Cognitive Systems Fakultät für Informatik Robot Dynamics Dr.-Ing. John Nassour 25.1.218 J.Nassour 1 Introduction Dynamics concerns the motion of bodies Includes Kinematics
More informationThermal and compositional structure of the Mantle and Lithosphere
Chapter 1 Thermal and compositional structure of the Mantle and Lithosphere 1.1 Primordial heat of the Earth The most widely accepted planetary formation theory says that the solar system accreted from
More informationGravity 3. Gravity 3. Gravitational Potential and the Geoid. Chuck Connor, Laura Connor. Potential Fields Geophysics: Week 2.
Gravitational Potential and the Geoid Chuck Connor, Laura Connor Potential Fields Geophysics: Week 2 Objectives for Week 1 Gravity as a vector Gravitational Potential The Geoid Gravity as a vector We can
More informationHigh-Harmonic Geoid Signatures due to Glacial Isostatic Adjustment, Subduction and Seismic Deformation
High-Harmonic Geoid Signatures due to Glacial Isostatic Adjustment, Subduction and Seismic Deformation L.L.A. Vermeersen (1), H. Schotman (1), M.-W. Jansen (1), R. Riva (1) and R. Sabadini (2) (1) DEOS,
More informationGOCE. Gravity and steady-state Ocean Circulation Explorer
GOCE Gravity and steady-state Ocean Circulation Explorer Reiner Rummel Astronomical and Physical Geodesy Technische Universität München rummel@bv.tum.de ESA Earth Observation Summerschool ESRIN/Frascati
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 informationChapter 9: Differential Analysis of Fluid Flow
of Fluid Flow Objectives 1. Understand how the differential equations of mass and momentum conservation are derived. 2. Calculate the stream function and pressure field, and plot streamlines for a known
More informationwhen viewed from the top, the objects should move as if interacting gravitationally
2 Elastic Space 2 Elastic Space The dynamics and apparent interactions of massive balls rolling on a stretched horizontal membrane are often used to illustrate gravitation. Investigate the system further.
More information1 Exercise: Linear, incompressible Stokes flow with FE
Figure 1: Pressure and velocity solution for a sinking, fluid slab impinging on viscosity contrast problem. 1 Exercise: Linear, incompressible Stokes flow with FE Reading Hughes (2000), sec. 4.2-4.4 Dabrowski
More 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 informationLecture 41: Highlights
Lecture 41: Highlights The goal of this lecture is to remind you of some of the key points that we ve covered this semester Note that this is not the complete set of topics that may appear on the final
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 informationSolutions to exam : 1FA352 Quantum Mechanics 10 hp 1
Solutions to exam 6--6: FA35 Quantum Mechanics hp Problem (4 p): (a) Define the concept of unitary operator and show that the operator e ipa/ is unitary (p is the momentum operator in one dimension) (b)
More informationCHAPTER 8 ENTROPY GENERATION AND TRANSPORT
CHAPTER 8 ENTROPY GENERATION AND TRANSPORT 8.1 CONVECTIVE FORM OF THE GIBBS EQUATION In this chapter we will address two questions. 1) How is Gibbs equation related to the energy conservation equation?
More informationEffect of Anisotropic Permeability on Band Growth
Effect of Anisotropic Permeability on Band Growth Jesse Taylor-West October 7, 24 Introduction The mechanics of partially molten regions of the mantle are not well understood- in part due to the inaccessibility
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 informationIterative 3D gravity inversion with integration of seismologic data
BOLLETTINO DI GEOFISICA TEORICA ED APPLICATA VOL. 40, N. 3-4, pp. 469-475; SEP.-DEC. 1999 Iterative 3D gravity inversion with integration of seismologic data C. BRAITENBERG and M. ZADRO Department of Earth
More informationOCN/ATM/ESS 587. The wind-driven ocean circulation. Friction and stress. The Ekman layer, top and bottom. Ekman pumping, Ekman suction
OCN/ATM/ESS 587 The wind-driven ocean circulation. Friction and stress The Ekman layer, top and bottom Ekman pumping, Ekman suction Westward intensification The wind-driven ocean. The major ocean gyres
More informationAstr 2320 Tues. May 2, 2017 Today s Topics Chapter 23: Cosmology: The Big Bang and Beyond Introduction Newtonian Cosmology Solutions to Einstein s
Astr 0 Tues. May, 07 Today s Topics Chapter : Cosmology: The Big Bang and Beyond Introduction Newtonian Cosmology Solutions to Einstein s Field Equations The Primeval Fireball Standard Big Bang Model Chapter
More informationSummary PHY101 ( 2 ) T / Hanadi Al Harbi
الكمية Physical Quantity القانون Low التعريف Definition الوحدة SI Unit Linear Momentum P = mθ be equal to the mass of an object times its velocity. Kg. m/s vector quantity Stress F \ A the external force
More information9. Density Structure. The Australian Continent: A Geophysical Synthesis Density Structure
84 The Australian Continent: A Geophysical Synthesis Density Structure 9. Density Structure Although the primary focus for the development of the AuSREM model was the creation of 3-D seismic wavespeed
More informationSummary of the Equations of Fluid Dynamics
Reference: Summary of the Equations of Fluid Dynamics Fluid Mechanics, L.D. Landau & E.M. Lifshitz 1 Introduction Emission processes give us diagnostics with which to estimate important parameters, such
More informationAE/ME 339. Computational Fluid Dynamics (CFD) K. M. Isaac. Momentum equation. Computational Fluid Dynamics (AE/ME 339) MAEEM Dept.
AE/ME 339 Computational Fluid Dynamics (CFD) 9//005 Topic7_NS_ F0 1 Momentum equation 9//005 Topic7_NS_ F0 1 Consider the moving fluid element model shown in Figure.b Basis is Newton s nd Law which says
More informationGravity observations and 3D structure of the Earth
1 Gravity observations and 3D structure of the Earth Y. Ricard 1, F. Chambat Laboratoire des Sciences de la Terre, UMR CNRS 5570, Lyon, France. C. Lithgow-Bertelloni Department of Earth Sciences, Michigan-Ann
More informationVery basic tsunami physics...
Very basic tsunami physics... Energy Bottom uplift & Waterberg formation E R 4.8 + 1.5M E T = 1 2 ρglλ(δh)2 L ~ 1 6 m; λ ~ 1 4 m; δh ~ 5m E R 1 18 J 1 2 E T Center of mass falls... Wavelength λ H ~ 4;
More informationB6 Isostacy. B6.1 Airy and Pratt hypotheses. Geophysics 210 September 2008
B6 Isostacy B6.1 Airy and Pratt hypotheses Himalayan peaks on the Tibet-Bhutan border In the 19 th century surveyors used plumblines and theodolites to map India. A plumb line was used when measuring the
More informationChapter 9: Differential Analysis
9-1 Introduction 9-2 Conservation of Mass 9-3 The Stream Function 9-4 Conservation of Linear Momentum 9-5 Navier Stokes Equation 9-6 Differential Analysis Problems Recall 9-1 Introduction (1) Chap 5: Control
More informationFinite-frequency tomography
Finite-frequency tomography Lapo Boschi (lapo@erdw.ethz.ch) October 14, 2009 Application of Born theory (scattering theory) to the ray-theory solution (the forward problem) In earlier lectures we have
More informationTracing rays through the Earth
Tracing rays through the Earth Ray parameter p: source receiv er i 1 V 1 sin i 1 = sin i 2 = = sin i n = const. = p V 1 V 2 V n p is constant for a given ray i 2 i 3 i 4 V 2 V 3 V 4 i critical If V increases
More 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 informationRANS Equations in Curvilinear Coordinates
Appendix C RANS Equations in Curvilinear Coordinates To begin with, the Reynolds-averaged Navier-Stokes RANS equations are presented in the familiar vector and Cartesian tensor forms. Each term in the
More information