28.08.2013 Views

Earth and Planetary Physics ? Problems ? 2 Planetary interiors

Earth and Planetary Physics ? Problems ? 2 Planetary interiors

Earth and Planetary Physics ? Problems ? 2 Planetary interiors

SHOW MORE
SHOW LESS

You also want an ePaper? Increase the reach of your titles

YUMPU automatically turns print PDFs into web optimized ePapers that Google loves.

Jacobs University Bremen<br />

Joachim Vogt<br />

Spring 2009<br />

2 <strong>Planetary</strong> <strong>interiors</strong><br />

<strong>Earth</strong> <strong>and</strong> <strong>Planetary</strong> <strong>Physics</strong><br />

— <strong>Problems</strong> —<br />

2.1 Heat production of a granite layer (Q)<br />

What is the thickness of a granite layer (density ρ = 2.6 · 10 3 kg m −3 <strong>and</strong> heat production rate<br />

1.07 · 10 −9 J kg −1 s −1 ) that yields the observed heat flow Q = 75 · 10 −3 W m −2 from the <strong>Earth</strong>’s<br />

interior through the surface ?<br />

2.2 Central pressure for uniform density <strong>and</strong> constant gravity (Q)<br />

To obtain a first crude estimate for the pressure Pc at the center of a spherical body, assume not<br />

only the density ρ to be uniform, but also the gravity g to be constant throughout the body (<strong>and</strong><br />

equal to its surface value). Use the hydrostatic equilibrium condition to compute Pc.<br />

2.3 Central pressure for uniform density <strong>and</strong> linear gravity profile (E)<br />

Derive the formula given in class for the pressure Pc at the center of a spherical uniform body<br />

in hydrostatic equilibrium, <strong>and</strong> verify the estimates for the minimum radii of spherical bodies.<br />

(a) Show that the radial gravity profile g = g(r) is linear in the radial coordinate r, i.e.,<br />

g(r) = gs<br />

where gs = −GM/R 2 denotes the surface value of the gravity.<br />

Hint: Note that in general ∇ · g = −4πGρ where ρ is the mass density, <strong>and</strong> G is the<br />

gravitational constant.<br />

(b) The resulting expression for g(r) will be used to integrate the hydrostatic equilibrium<br />

condition, <strong>and</strong> to obtain the radial pressure profile P(r). Express the uniform mass density<br />

ρ in terms of the total mass M of the planet.<br />

1<br />

r<br />

R


(c) To estimate the minimum radii of (quasi-)spherical objects, insert numerical values of<br />

the material strengths for Pc. For rocky bodies, use ∼ 2 · 10 8 Pa <strong>and</strong> an (uncompressed)<br />

density ρ = 3.5 · 10 3 kg/m 3 . For iron bodies, the material strength is ∼ 4 · 10 8 Pa <strong>and</strong><br />

ρ = 7.9 · 10 3 kg/m 3 .<br />

2.4 <strong>Earth</strong>’s moment of inertia <strong>and</strong> core density (E)<br />

Use the values RE <strong>and</strong> ME of the spherical <strong>Earth</strong> model <strong>and</strong> the core radius Rc = 3480 km (infered<br />

from seismological recordings) to determine the core density ρc <strong>and</strong> the mantle density<br />

ρm of a simplified two-layer model of the <strong>Earth</strong>’s interior. The model comprises only the mantle<br />

<strong>and</strong> the core. Suppose that both regions are characterized by homogeneous mass density<br />

functions ρm <strong>and</strong> ρc.<br />

Obviously, you need additional information to solve the problem.<br />

(a) Use a typical surface rock density value as a proxy for the mantle density (ρm = 3 g/cm 3 )<br />

<strong>and</strong> then determine ρc.<br />

(b) Compute the moment of inertia for the spherical two-layer <strong>Earth</strong> model <strong>and</strong> identiy this<br />

theoretical value with the measured value of the <strong>Earth</strong>’s (main) moment of inertia C given<br />

in class. Now determine, both, the mantle density ρm <strong>and</strong> the core density ρc from the two<br />

equations for the total mass <strong>and</strong> the total angular momentum.<br />

Hint: The equations become simpler if normalized variables are used. For normalization of the<br />

densities you may use the mean density of the whole <strong>Earth</strong>:<br />

ρ E =<br />

mass of the <strong>Earth</strong><br />

volume of the <strong>Earth</strong><br />

<strong>and</strong> the core radius can be normalized by RE. This yields<br />

ˆρm = ρm/ρ E ,<br />

ˆρc = ρc/ρ E ,<br />

ˆRc = Rc/RE .<br />

3ME<br />

=<br />

4πR3 ,<br />

E<br />

2.5 Temperature, pressure, <strong>and</strong> density in the <strong>Earth</strong>’s mantle [H3]<br />

In this problem you are supposed to determine approximately the profiles of temperature T,<br />

pressure P, <strong>and</strong> mass density ρ in the <strong>Earth</strong>’s mantle.<br />

(a) In order to find the temperature profile T = T(r), assume that the gradient is adiabatic:<br />

<br />

dT <br />

<br />

α g<br />

= − T .<br />

dr<br />

cp<br />

adiab<br />

2


Use constant values for the parameters in this equation: gravity g = 9.81 m s −2 = go<br />

(surface value), specific heat (at constant pressure) cp = 10 3 J kg −1 K −1 , <strong>and</strong> thermal expansion<br />

coefficient α = 3 · 10 −5 K −1 . The temperature boundary condition at the surface<br />

is To = 300 K.<br />

(b) Compute P = P(r) <strong>and</strong> ρ = ρ(r) using the hydrostatic equilibrium condition <strong>and</strong> a fixed<br />

value of the adiabatic compressibility: β = 4.3 · 10−12 Pa −1 . Note that<br />

β = 1<br />

<br />

∂ρ <br />

.<br />

ρ ∂P<br />

adiab<br />

The boundary values are ρo = 3.3 · 10 3 kg m −3 <strong>and</strong> Po = 0 Pa.<br />

(c) Plot ρ, P <strong>and</strong> T as functions of radial distance r.<br />

2.6 Interior of gas giants [H3]<br />

This is problem 6.28 in the textbook <strong>Planetary</strong> Sciences of Imke de Pater <strong>and</strong> Jack J. Lissauer.<br />

Approximate Jupiter <strong>and</strong> Saturn by pure hydrogen spheres, with an equation of state P = Kρ 2<br />

<strong>and</strong> K = 2.7 · 10 12 cm 5 g −1 s −2 .<br />

(a) Determine the moment of inertia for the planets.<br />

(b) Assume the planets each have a core of density 10 g cm −3 , <strong>and</strong> the moment of inertia<br />

I/MR 2 = 0.253 for Jupiter <strong>and</strong> 0.227 for Saturn. Determine the mass of the cores.<br />

Submission due date for homework problems: 9 March 2009<br />

3

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!