Earth and Space Sciences · Book 4 · Bachelor Year 2

University Earth and Space Sciences — Year 2

University Earth and Space Sciences — Year 2 · Bachelor Year 2

2Planet Formation and Meteorites

Every year some fifteen thousand tonnes of dust from space enter the Earth’s atmosphere; most of it burns up or settles as microscopic spherules, and once in a while a stone the size of a fist falls with a fireball and a sonic boom. Most of these stones are pieces of asteroids, and many have hardly changed since they formed. Cut one open and polish it, and it shows millimetre-sized beads of once-molten rock, white irregular grains of calcium and aluminium minerals, flakes of metal, all held in a dark, fine-grained cement. The white grains are the oldest solids of the Solar System: lead isotopes date them to 4567 million years. A second clock, the decay of a radioactive isotope of aluminium that has long since vanished, dates the other components relative to them to within a fraction of a million years. This chapter follows the disc of the previous chapter as it turns into planets: which solids condense from the cooling gas, and where; how dust grows into bodies kilometres across and then into planets; what the meteorites are, and how they date these events; and why some of the first bodies melted and grew iron cores while others stayed as they had formed.

You already know

Young stars are surrounded by protoplanetary discs, which disappear within a few million years; a grain in starlight at the distance rr from a star of luminosity LL settles at T=(L/(16πσr2))1/4T = \bigl(L/(16\pi\sigma r^2)\bigr)^{1/4} (Chapter 1). Asteroids are the small rocky or metallic bodies of the Solar System, most of which orbit between Mars and Jupiter (Year 1 volume). A radioactive nucleus decays with the probability λ\lambda per unit time: N=N0e−λtN = N_0 e^{-\lambda t}; the half-life is ln⁡2/λ\ln 2/\lambda, and the mean life τ=1/λ\tau = 1/\lambda (physics). In a closed system, cogenetic samples of different parent-to-daughter ratios lie on a straight line, the isochron, whose slope gives the age and whose intercept is the initial ratio (Year 1 volume). Heating a mass mm of specific heat capacity cc by ΔT\Delta T takes the energy mc ΔTmc\,\Delta T (physics).

2.1 The solar nebula

Definition 2.1 (Solar nebula)

The solar nebula is the protoplanetary disc of gas and dust that surrounded the young Sun and from which the planets, their satellites, the asteroids and the comets formed.

How much matter did it hold? A lower limit comes from the planets themselves. Add to each planet enough hydrogen and helium to give it the composition of the Sun, spread the result over a ring centred on its orbit, and the surface density, the mass per unit area of the disc, falls off with the distance from the Sun. This reconstruction is the minimum-mass solar nebula. In a widely used version the surface density of the gas is

Σg=2200 g cm−2(r1 au)−3/2,\Sigma_g = 2200\,\mathrm{g}\,\mathrm{cm}^{-2}\left(\frac{r}{1\,\mathrm{au}}\right)^{-3/2},

and that of the solids Σp=33 g cm−2 Zrel (r/1 au)−3/2\Sigma_p = 33\,\mathrm{g}\,\mathrm{cm}^{-2}\,Z_{\mathrm{rel}}\, (r/1\,\mathrm{au})^{-3/2}, where ZrelZ_{\mathrm{rel}} is the fraction of the condensable matter actually solid: 0.330.33 where only rock and metal have condensed, about 0.780.78 where water ice has joined them, 11 beyond the place where methane freezes. The numbers are uncertain: the heavy-element content of Jupiter is known only to within a factor of four, and the planets need not have formed where they orbit now.

Proposition 2.2 (Mass of a disc)

A disc of surface density Σ(r)=Σ1(r/r1)−3/2\Sigma(r) = \Sigma_1 (r/r_1)^{-3/2} holds, inside the radius RR, the mass

M(R)=4πΣ1r12(Rr1)1/2.M(R) = 4\pi\Sigma_1 r_1^2 \left(\frac{R}{r_1}\right)^{1/2}.

Most of the mass lies far out: doubling the radius adds 41 %41\,\% to the mass.

Proof. A ring of radius rr and width dr\mathrm dr has the area 2πr dr2\pi r\,\mathrm dr, so

M(R)=∫0R2πr Σ1(rr1)−3/2dr=2πΣ1r13/2∫0Rr−1/2 dr=4πΣ1r13/2R1/2,M(R) = \int_0^R 2\pi r\,\Sigma_1 \left(\frac{r}{r_1}\right)^{-3/2}\mathrm dr = 2\pi\Sigma_1 r_1^{3/2}\int_0^R r^{-1/2}\,\mathrm dr = 4\pi\Sigma_1 r_1^{3/2} R^{1/2},

which is the formula; the integral converges at r=0r = 0. Doubling RR multiplies MM by 2=1.41\sqrt 2 = 1.41. ∎

With Σ1=2.2×104 kg m−2\Sigma_1 = 2.2 \times 10^{4}\,\mathrm{kg}\,\mathrm{m}^{-2} at r1=1 aur_1 = 1\,\mathrm{au}, the minimum-mass nebula holds 0.030.03 solar masses inside 100 au100\,\mathrm{au} (Exercise 2.5). The discs observed around stars a million years old, weighed by the millimetre glow of their dust, hold 0.001 to 0.10.001\text{ to }0.1 solar masses: the minimum-mass nebula sits comfortably among them.

Almost all of this mass is gas that will never condense. The Sun was born with a mass fraction of hydrogen X0=0.711X_0 = 0.711, of helium Y0=0.274Y_0 = 0.274 and of all heavier elements together only Z0=0.0149Z_0 = 0.0149. Cooling a gas of this composition until everything that can freeze has frozen turns only 1.5 %1.5\,\% of its mass into solids: 0.49 %0.49\,\% of rock and metal, and 1.0 %1.0\,\% of ices of water, methane and ammonia, water ice alone being 0.57 %0.57\,\%. All the rocky planets, and the cores of the giants, were built from this small fraction.

The temperature of the disc decides what is solid where. A disc heated only by its star, a passive disc, is warm at its surface, where the grains see the star directly (the proposition of Chapter 1), and much colder in its midplane, which receives only the infrared re-emitted by those surface grains. For a young star of one solar mass a model of such a disc gives the midplane temperature

T=120 K(r1 au)−3/7.T = 120\,\mathrm{K}\left(\frac{r}{1\,\mathrm{au}}\right)^{-3/7}.

Accretion of gas through the disc, which releases gravitational energy, makes the midplane hotter, especially early on and close to the star (figure below).

Two temperatures of a disc around a star of one solar luminosity. Red: a black grain in direct starlight (). Blue: the midplane of a passive disc around a young solar-mass star, shielded from direct light. Water ice is stable below 182\, K: the arrows mark where each curve crosses that line. The real snow line lay between such limits and moved as the disc evolved. Data: Chiang and Youdin (2010); Lodders (2003).
Two temperatures of a disc around a star of one solar luminosity. Red: a black grain in direct starlight (Chapter 1). Blue: the midplane of a passive disc around a young solar-mass star, shielded from direct light. Water ice is stable below 182 K182\,\mathrm{K}: the arrows mark where each curve crosses that line. The real snow line lay between such limits and moved as the disc evolved. Data: Chiang and Youdin (2010); Lodders (2003).

2.2 The condensation sequence

Imagine the gas of the nebula, close to the young Sun, hot enough for every element to be vapour, and let it cool slowly enough for gas and solids to remain in chemical equilibrium. As the temperature falls, each element leaves the gas for a solid in turn, either as a compound of its own or dissolved in a solid that has already formed.

Definition 2.3 (Condensation temperature, condensation sequence)

The condensation temperature of an element in a cooling gas of given composition and pressure is the temperature at which half of it has passed into solids at chemical equilibrium (the 50 per cent condensation temperature; the temperature at which its first solid appears is the appearance temperature). The condensation sequence is the order in which the elements and their solid compounds appear as the gas cools.

Calculations for a gas of solar composition are made at a total pressure of 10−410^{-4} bar, typical of the nebula near 1 au1\,\mathrm{au}. The first solids are oxides of aluminium and calcium (table below), which form at about 1700 K1700\,\mathrm{K}; the bulk of the rock-forming matter, magnesium silicate and iron metal, condenses near 1350 K1350\,\mathrm{K}; iron sulfide waits until 704 K704\,\mathrm{K}, and water ice until 182 K182\,\mathrm{K}.

solidformulaappears at (K)
corundumAlX2OX3\ce{Al2O3}1677
hiboniteCaAlX12OX19\ce{CaAl12O19}1659
perovskiteCaTiOX3\ce{CaTiO3}1593
spinelMgAlX2OX4\ce{MgAl2O4}1397
iron–nickel metalFe\ce{Fe}, with Ni\ce{Ni}, Co\ce{Co}, Cr\ce{Cr}1357
forsteriteMgX2SiOX4\ce{Mg2SiO4}1354
schreibersiteFeX3P\ce{Fe3P}1248
troiliteFeS\ce{FeS}704
magnetiteFeX3OX4\ce{Fe3O4}371
water iceHX2O\ce{H2O}182
ammonia hydrateNHX3 ⋅ HX2O\ce{NH3.H2O}131
methane clathrateCHX4 ⋅ 7 HX2O\ce{CH4.7H2O}78
Table 2.1. The condensation sequence of a gas of solar-system composition at 10−410^{-4} bar: the temperature at which each solid first appears. Sodium and potassium enter feldspar; nickel, cobalt and chromium dissolve in the metal. Data: Lodders (2003).

The 50 per cent condensation temperatures sort the elements into groups of volatility (figure below). Elements that condense above magnesium, silicon and iron are refractory; those that condense below them, volatile. Lodders divides each side further: the ultrarefractory elements, above 1650 K1650\,\mathrm{K}, include aluminium, zirconium, the heavy rare earths and the metals osmium, rhenium and tungsten; the highly refractory ones (1500 to 1650 K1500\text{ to }1650\,\mathrm{K}) include calcium and titanium; the refractory ones (1360 to 1500 K1360\text{ to }1500\,\mathrm{K}) include strontium, barium, vanadium and platinum; the common elements magnesium, silicon and iron condense between 12901290 and 1360 K1360\,\mathrm{K}, together with nickel, cobalt and chromium; the moderately volatile elements, among them sodium, potassium, phosphorus, copper, zinc and lead, condense between 1290 K1290\,\mathrm{K} and the appearance of troilite at 704 K704\,\mathrm{K}; the volatile ones, among them sulfur, cadmium and thallium, between 704704 and 371 K371\,\mathrm{K}; and the highly volatile elements, carbon, nitrogen, oxygen as ice, the noble gases and mercury, below that.

The 50 per cent condensation temperatures of 36 elements in a gas of solar-system composition at 10-4 bar, coloured by Lodders’ volatility groups (oxygen: as water ice; carbon: as methane clathrate and ice). Data: Lodders (2003).
The 50 per cent condensation temperatures of 36 elements in a gas of solar-system composition at 10−410^{-4} bar, coloured by Lodders’ volatility groups (oxygen: as water ice; carbon: as methane clathrate and ice). Data: Lodders (2003).

The sequence explains a pattern that runs through the whole Solar System. The oldest solids in meteorites are made of the first condensates; the Earth and the other rocky planets are poor in volatile elements; the far reaches of the disc, cold enough to keep ices, made icy bodies. The most important boundary is the one for water, because water ice is the most abundant solid of all: by mass, there is about as much water ice in a gas of solar composition as rock and metal together, 1.171.17 times more.

Definition 2.4 (Snow line)

The snow line of a protoplanetary disc is the distance from the star beyond which the midplane is cold enough for water ice to be stable, about 180 K180\,\mathrm{K} at the pressures of the disc.

Proposition 2.5 (Snow line of a passive disc)

If the midplane temperature falls as T=T1(r/r1)−3/7T = T_1 (r/r_1)^{-3/7}, the snow line of water, which condenses at TwT_w, lies at

rs=r1(T1Tw)7/3,r_s = r_1 \left(\frac{T_1}{T_w}\right)^{7/3},

and beyond it the surface density of solids jumps by the factor 1+mw/mr1 + m_w/m_r, where mwm_w and mrm_r are the masses of water ice and of rock in the gas.

Proof. Solving Tw=T1(rs/r1)−3/7T_w = T_1 (r_s/r_1)^{-3/7} for rsr_s gives rs/r1=(T1/Tw)7/3r_s/r_1 = (T_1/T_w)^{7/3}. Inside rsr_s only rock and metal are solid; beyond it the water ice adds its mass mwm_w to the rock mrm_r of the same gas, so the solid mass is multiplied by (mr+mw)/mr(m_r + m_w)/m_r. ∎

For the passive disc above, T1=120 KT_1 = 120\,\mathrm{K} and Tw=182 KT_w = 182\,\mathrm{K} give rs=0.38 aur_s = 0.38\,\mathrm{au}; methane, at 41 K41\,\mathrm{K}, freezes beyond 12 au12\,\mathrm{au} (Exercise 2.6). With 0.57 %0.57\,\% of water ice and 0.49 %0.49\,\% of rock, the solids jump by a factor 1+0.571/0.489=2.21 + 0.571/0.489 = 2.2. The model places the snow line far inside the Earth’s orbit, which shows only that a young disc was not passive: heated by accretion, its inner midplane was warmer, and the snow line lay farther out, moving inward as the disc thinned. Beyond it, the planets had more than twice as much solid matter to build from, and this, more than anything else, separates the small rocky planets from the giants.

2.3 From dust to planets

Between the micrometre grains of the interstellar medium and planets thousands of kilometres across lie some thirteen orders of magnitude of size. The growth is best described in three legs. Below about a centimetre, grains that collide gently stick, held by chemical bonds and van der Waals forces, and build fluffy aggregates. Well above a kilometre, bodies hold together by their own gravity, and the fragments of a collision fall back onto the larger body. Between the two lies the hardest leg: bodies from centimetres to metres are too large to stick and too small to attract one another, and the gas drags them inward. In numerical simulations, solids 10 to 100 cm10\text{ to }100\,\mathrm{cm} across gather into dense swarms that collapse under their own weight directly into bodies kilometres across, skipping the troublesome sizes; how the Solar System crossed this leg is still debated.

Definition 2.6 (Planetesimal, planetary embryo)

A planetesimal is a solid body of a protoplanetary disc large enough, from about a kilometre to a few hundred kilometres across, for its growth to be governed by gravity rather than by sticking. A planetary embryo is a body that has grown, by accreting planetesimals, large enough to dominate its own region of the disc and to sweep it clear of planetesimals.

Once gravity rules, growth runs away: a body slightly larger than its neighbours draws planetesimals towards itself from farther away, grows faster than they do, and pulls ahead. Its growth slows only when it has swept up the solids within its reach. Each embryo then rules its own annulus of the disc, and many embryos share the region of the rocky planets.

Proposition 2.7 (Isolation mass)

An embryo of mass MM at the distance rr from a star of mass M∗M_* attracts the solids of an annulus whose half-width is a few times its Hill radius RH=r (M/3M∗)1/3R_H = r\,(M/3M_*)^{1/3}. If it gathers all the solids of an annulus of total width bRHbR_H and surface density Σp\Sigma_p, it stops growing at the isolation mass

Miso=(2πb r2Σp)3/2(3M∗)1/2.M_{\mathrm{iso}} = \frac{(2\pi b\, r^2 \Sigma_p)^{3/2}}{(3M_*)^{1/2}} .

Proof. The annulus of radius rr and width bRHbR_H holds the mass 2πr bRH Σp2\pi r\, bR_H\,\Sigma_p. The embryo has swallowed it when M=2πbrΣp r (M/3M∗)1/3M = 2\pi b r\Sigma_p\, r\,(M/3M_*)^{1/3}, that is M2/3=2πbr2Σp/(3M∗)1/3M^{2/3} = 2\pi b r^2\Sigma_p/(3M_*)^{1/3}; raising both sides to the power 3/23/2 gives the formula. ∎

Models take b≈5b \approx 5. In the minimum-mass nebula the isolation mass is about 0.050.05 Earth masses at 1 au1\,\mathrm{au}, between the masses of the Moon (0.0120.012) and of Mars (0.110.11): the embryos of the inner Solar System were bodies of that size, and the Earth needed many of them. Beyond the snow line, at 5.2 au5.2\,\mathrm{au}, the ice-rich solids give about 0.80.8 Earth masses, and a disc with 3.53.5 times the minimum mass of solids gives some five Earth masses (Exercise 2.7): enough, in models of Jupiter’s growth, for a core to pull in its gas envelope within three million years. The giant planets had to form before the gas of the disc was gone, within a few million years; Uranus and Neptune are 80 to 90 %80\text{ to }90\,\% rock and ice by mass: they gathered little gas.

The last stage of the rocky planets was violent. Once the gas was gone, the embryos’ orbits crossed, and they merged in giant collisions; the Moon is thought to have formed from the debris of one of them, which struck the young Earth. The figure below sums up the sequence.

From dust to planets in the region of the rocky planets (schematic, not to scale). Grains stick into aggregates; swarms of pebbles collapse into planetesimals; gravity makes some planetesimals run away into embryos, which end by merging into planets. Beyond the snow line the embryos grew larger and could capture gas. Sizes from Chiang and Youdin (2010) and Lugaro and others (2018).
From dust to planets in the region of the rocky planets (schematic, not to scale). Grains stick into aggregates; swarms of pebbles collapse into planetesimals; gravity makes some planetesimals run away into embryos, which end by merging into planets. Beyond the snow line the embryos grew larger and could capture gas. Sizes from Chiang and Youdin (2010) and Lugaro and others (2018).

The planets kept little memory of this history: they melted, mixed and reworked their matter. The record is in the bodies that stopped growing, the asteroids, and in the pieces of them that reach the Earth.

2.4 Meteorites

Definition 2.8 (Meteorite, meteorite fall, meteorite find)

A meteorite is a natural solid object from space that has survived its passage through the atmosphere and reached the ground. A meteorite fall is a meteorite collected after its arrival was observed, usually as a fireball; a meteorite find is one recognised on the ground without its fall having been seen.

The Earth sweeps up far more dust than stones. Measured in the snow of the Antarctic plateau, about 5200 t5200\,\mathrm{t} of extraterrestrial dust reaches the ground each year as micrometeorites and molten spherules, out of some 15 000 t15\,000\,\mathrm{t} that enter the atmosphere, most of it shed by comets and about a fifth by asteroids. Meteorites, the stones and irons large enough to be recognised one by one, are fewer: the world’s catalogue holds about 45 00045\,000 of them, of which only 11071107 are falls.

The distinction matters. A fall is collected within days, before rain and oxygen have altered it, and every kind of meteorite that falls is collected alike; a find may have lain on the ground for thousands of years, and its chance of being found depends on how it looks and where it fell. Half of all finds have been made in Antarctica (box below), where the search is systematic. The falls are the fairer sample of what arrives (figure below): more than four in five are chondrites, three in four ordinary chondrites.

Meteorites by class, as a percentage of the falls and of the finds (logarithmic scale); HED: howardites, eucrites and diogenites. Data: NASA export of the Meteoritical Bulletin Database of The Meteoritical Society (45\,641 valid meteorites).
Meteorites by class, as a percentage of the falls and of the finds (logarithmic scale); HED: howardites, eucrites and diogenites. Data: NASA export of the Meteoritical Bulletin Database of The Meteoritical Society (45 64145\,641 valid meteorites).

In the field — Searching for meteorites on the blue ice

A line of searchers on snowmobiles crossing a field of blue ice below a rocky peak (AI-generated illustration).

Most of the Antarctic ice sheet is white snow, but on about 1 %1\,\% of its surface, where mountains block the flow of the ice, some 5 cm5\,\mathrm{cm} of ice is lost each year and bare blue ice shows. Meteorites that fell on the snow upstream, over tens of thousands of years, are buried, carried along in the ice and brought back to the surface where it wastes away: the blue ice is a trap. Since systematic collecting began in the mid-1970s, teams have recovered more than 30 00030\,000 specimens from such fields, and in the catalogue half of all finds come from south of 60∘60{}^{\circ} S. The searchers cross the ice in lines, on foot or on snowmobiles, and every stone found is collected, however ordinary it looks, so that the collection is not biased towards the striking ones. Each find is photographed where it lies, to record its context, and given a field number; later it receives its official name, a three-letter code for the ice field followed by the year of the season and a serial number. In the curation laboratory the stones are chipped with stainless-steel tools or cut dry with a band saw; carbonaceous chondrites, which may hold organic compounds, have had a cabinet of their own since 1978.

2.4.1 Chondrites

Definition 2.9 (Chondrite, carbonaceous chondrite, chondrule, CAI)

A chondrite is a stony meteorite that never melted after it was assembled, and whose composition, apart from the most volatile elements, is close to that of the Sun. Most of it is made of chondrules, rounded droplets, a tenth of a millimetre to about a centimetre across, of silicate rock that solidified from a melt, set in a fine-grained matrix with grains of iron–nickel metal and iron sulfide. Some chondrites hold calcium–aluminium-rich inclusions (CAIs), irregular white grains of the minerals of calcium, aluminium and titanium that are the first condensates of the nebula. A carbonaceous chondrite is a chondrite rich in carbon, water-bearing minerals and matrix, whose content of the volatile elements comes closest to the Sun’s.

The chondrules were molten for only minutes and cooled over hours to days: some brief, local flash heating in the nebula, whose cause is still unknown, melted clumps of dust, which froze as droplets before they were swept together with dust and metal into a planetesimal. They fill up to 90 %90\,\% of the volume of the most primitive meteorites. The CAIs make up from less than 1 %1\,\% of ordinary chondrites to over 10 %10\,\% of some carbonaceous ones, and they formed within about a hundred thousand years of one another.

A cut face of the Allende carbonaceous chondrite (class CV3), which fell in northern Mexico in 1969: about two tonnes were collected. The rounded grey grains are chondrules; the irregular white patches are CAIs; the dark cement is the matrix. The field is a few centimetres wide. James St. John, CC BY 2.0.
A cut face of the Allende carbonaceous chondrite (class CV3), which fell in northern Mexico in 1969: about two tonnes were collected. The rounded grey grains are chondrules; the irregular white patches are CAIs; the dark cement is the matrix. The field is a few centimetres wide. James St. John, CC BY 2.0.

Chondrites fall into groups that differ in the size and abundance of their chondrules, their content of metal, the iron content of their olivine and the proportion of their matrix (table below). The ordinary chondrites H, L and LL have high, low and very low contents of metal and of total iron (atomic ratios Fe/Si of 0.800.80, 0.590.59 and 0.530.53).

HLLLCICMCV
chondrules (vol. %)60–8060–8060–80≪1\ll 12045
mean chondrule diameter (mm)0.30.70.9—0.31.0
CAIs (vol. %)0.1–10.1–10.1–1≪1\ll 1510
matrix (vol. %)10–1510–1510–15>99> 997040
iron–nickel metal (vol. %)105200.10–5
fayalite in olivine (mol %)16–2023–2627–32
petrologic types3–63–63–6123
highest temperature (K)122012201220430670870
Table 2.2. Properties of the ordinary (H, L, LL) and of three carbonaceous chondrite groups. Fayalite: the iron end-member of olivine. Data: Gail and others (2014).

Definition 2.10 (Petrologic type)

The petrologic type of a chondrite is a number from 1 to 6 that records what happened to it in its parent body: type 3 is the least altered; types 4 to 6 are increasingly heated, their chondrules blurred and their minerals recrystallised and of uniform composition; types 2 and 1 are increasingly altered by liquid water, which turned their minerals into clays.

The last row of the table gives the highest temperatures reached by each group: up to about 1220 K1220\,\mathrm{K} for the most heated ordinary chondrites, never more than 430 K430\,\mathrm{K} for the CI chondrites, soaked in water but never hot. The heating, as the last section shows, came from inside the parent bodies.

2.4.2 Achondrites and irons

Definition 2.11 (Achondrite, iron meteorite, stony-iron meteorite)

An achondrite is a stony meteorite that crystallised from a melt or recrystallised at high temperature in a body that melted, and so has no chondrules. An iron meteorite is made mostly of iron–nickel metal. A stony-iron meteorite is made of metal and silicate in roughly equal amounts.

These are the pieces of bodies that melted and separated: achondrites are their crusts and mantles, irons their metal cores, and stony-irons come from where metal and silicate met. The largest family of achondrites, the howardites, eucrites and diogenites (HED), basalts and pyroxene-rich rocks, makes up 5.5 %5.5\,\% of all falls. Its source was traced, by its spectrum and its orbit, to the asteroid Vesta, about 525 km525\,\mathrm{km} across, and the Dawn spacecraft, which reached Vesta in 2011, confirmed it: the chemistry and mineralogy of Vesta’s surface match the HED meteorites, a basin 500 km500\,\mathrm{km} wide near its south pole shows where the fragments were knocked off, and its gravity field points to an iron core about 225 km225\,\mathrm{km} across. Vesta melted, nearly or entirely, early in its history. Other achondrites come from the Moon and from Mars. The meteorites in the collections appear to come from about a hundred different parent bodies.

Iron meteorites keep the record of how their cores cooled. A cut and acid-etched face shows a criss-cross of plates (figure below): at high temperature the metal is a single alloy of iron and nickel; on slow cooling, a nickel-poor alloy grows in plates along the planes of the original crystal, and the slower the cooling, the wider the plates. Measured on many irons, the plate widths give the cooling rates of the cores they came from.

The cut, polished and acid-etched face of a piece of the Gibeon iron meteorite, found in Namibia: the criss-cross of plates of iron–nickel alloy (the Widmanstätten pattern) grew as the metal core of its parent body cooled slowly. James St. John, CC BY 2.0.
The cut, polished and acid-etched face of a piece of the Gibeon iron meteorite, found in Namibia: the criss-cross of plates of iron–nickel alloy (the Widmanstätten pattern) grew as the metal core of its parent body cooled slowly. James St. John, CC BY 2.0.

2.5 Dating the Solar System

The age of the Solar System is the age of its oldest solids. Uranium–lead isochrons of CAIs, built as in the Year 1 volume, give 4567 to 4568 Ma4567\text{ to }4568\,\mathrm{Ma}; the CAIs formed within about 0.1 Myr0.1\,\mathrm{Myr} of one another, a time short enough for their age to be taken as the birth of the Sun. The events of the next few million years, the formation of chondrules and the assembly and melting of the first bodies, are best dated relative to the CAIs, with clocks that ran fast and stopped early.

Definition 2.12 (Short-lived radionuclide)

A short-lived radionuclide is a radioactive nuclide with a half-life of roughly 0.10.1 to 100 Myr100\,\mathrm{Myr}, present when the Solar System formed and now extinct, whose former presence is shown by an excess of its daughter in minerals that were rich in its element.

The best known is aluminium-26, which decays to magnesium-26 with a half-life of 0.717 Myr0.717\,\mathrm{Myr} (a mean life of 1.035 Myr1.035\,\mathrm{Myr}). When the CAIs formed, one aluminium atom in about 19 00019\,000 was aluminium-26: the ratio X26X2226Al/X27X2227Al\ce{^{26}Al}/\ce{^{27}Al} was (5.23±0.13)×10−5(5.23 \pm 0.13)\times10^{-5}, the canonical value. Other short-lived nuclides cover longer spans (table below).

nuclidedaughterreferencehalf-life (Myr)early ratio
X26X2226Al\ce{^{26}Al}X26X2226Mg\ce{^{26}Mg}X27X2227Al\ce{^{27}Al}0.7175.23×10−55.23\times10^{-5}
X60X2260Fe\ce{^{60}Fe}X60X2260Ni\ce{^{60}Ni}X56X2256Fe\ce{^{56}Fe}2.621.0×10−81.0\times10^{-8}
X53X2253Mn\ce{^{53}Mn}X53X2253Cr\ce{^{53}Cr}X55X2255Mn\ce{^{55}Mn}3.747×10−67\times10^{-6}
X107X22107Pd\ce{^{107}Pd}X107X22107Ag\ce{^{107}Ag}X108X22108Pd\ce{^{108}Pd}6.56.6×10−56.6\times10^{-5}
X182X22182Hf\ce{^{182}Hf}X182X22182W\ce{^{182}W}X180X22180Hf\ce{^{180}Hf}8.901.02×10−41.02\times10^{-4}
X129X22129I\ce{^{129}I}X129X22129Xe\ce{^{129}Xe}X127X22127I\ce{^{127}I}15.71.28×10−41.28\times10^{-4}
Table 2.3. Some short-lived radionuclides of the early Solar System, their stable daughter, the stable isotope of the same element used as reference, their half-life and their ratio to the reference when the CAIs formed. Data: Lugaro and others (2018).

Proposition 2.13 (Relative age from an extinct nuclide)

If the ratio of a short-lived nuclide to its reference isotope was R0R_0 everywhere when the CAIs formed, and is found to have been RR in an object when it formed, the object formed

Δt=τln⁡R0R\Delta t = \tau \ln\frac{R_0}{R}

after the CAIs, where τ\tau is the mean life of the nuclide. An error δR\delta R on RR gives the error τ δR/R\tau\,\delta R/R on Δt\Delta t.

Proof. The reference isotope is stable, so the ratio decays like the radioactive nuclide itself: R=R0e−Δt/τR = R_0 e^{-\Delta t/\tau}. Taking logarithms, Δt=τln⁡(R0/R)\Delta t = \tau\ln(R_0/R). Differentiating, d(Δt)=−τ dR/R\mathrm d(\Delta t) = -\tau\,\mathrm dR/R. ∎

A relative error of 10 %10\,\% on the ratio thus costs 0.1 Myr0.1\,\mathrm{Myr} with aluminium-26 but 0.54 Myr0.54\,\mathrm{Myr} with manganese-53 and 1.3 Myr1.3\,\mathrm{Myr} with hafnium-182: each clock is precise only over a few of its own mean lives, and the right clock depends on the event. The ratio RR in an object is read from its minerals.

Proposition 2.14 (Aluminium–magnesium isochron)

Let several minerals of one object form together at the time Δt\Delta t with the same isotopic composition of magnesium, (X26X2226Mg/X24X2224Mg)0(\ce{^{26}Mg}/\ce{^{24}Mg})_0, but different ratios X27X2227Al/X24X2224Mg\ce{^{27}Al}/\ce{^{24}Mg}, and stay closed since. Today they lie on the straight line

X26X2226MgX24X2224Mg=(X26X2226MgX24X2224Mg)0+R X27X2227AlX24X2224Mg,\frac{\ce{^{26}Mg}}{\ce{^{24}Mg}} = \left(\frac{\ce{^{26}Mg}}{\ce{^{24}Mg}}\right)_0 + R\,\frac{\ce{^{27}Al}}{\ce{^{24}Mg}},

whose slope is the ratio R=X26X2226Al/X27X2227AlR = \ce{^{26}Al}/\ce{^{27}Al} of the object when it formed.

Proof. Every aluminium-26 atom present at formation has since become a magnesium-26 atom, so a mineral gained X26X2226Althen=R X27X2227Al\ce{^{26}Al}_{\text{then}} = R\,\ce{^{27}Al} atoms of magnesium-26, aluminium-27 being stable. Dividing by the number of magnesium-24 atoms, which is unchanged, gives the line. ∎

Aluminium-rich, magnesium-poor minerals such as plagioclase show large excesses of magnesium-26; olivine, almost free of aluminium, shows none. This is how the time sequence of the early Solar System was built: CAIs first; chondrules over the following few million years, with the evidence from chondrules setting the lifetime of the gas disc at less than five million years; the parent bodies of the chondrites assembled after their youngest chondrules. The diogenites, rocks of Vesta’s mantle, began to crystallise within about a million years of the CAIs, and Vesta’s core formed within the first one to four million years: iron meteorites and achondrites are older than many chondrites.

2.6 Heat from aluminium-26 and the first melting

Where did the heat that melted Vesta, and metamorphosed the ordinary chondrites, come from? Models and meteorites agree that in the earliest bodies the dominant source was the decay of aluminium-26. In rock of chondritic composition it carried, at the start, 5 MJ5\,\mathrm{MJ} per kilogram ready to be released within a few million years: enough to heat the rock by thousands of kelvins, if the heat could not escape.

Proposition 2.15 (Heating of a planetesimal by aluminium-26)

Rock that receives the radioactive power H0e−t/τH_0 e^{-t/\tau} per kilogram, has the specific heat capacity cc and loses no heat, warms after its accretion at the time tat_a as

T(t)=T0+H0τc(e−ta/τ−e−t/τ),T(t) = T_0 + \frac{H_0\tau}{c}\left(e^{-t_a/\tau} - e^{-t/\tau}\right),

and tends to the peak temperature Tmax⁡=T0+(H0τ/c) e−ta/τT_{\max} = T_0 + (H_0\tau/c)\,e^{-t_a/\tau}. A critical temperature Tc>T0T_c > T_0 is reached only if the body accreted before

ta∗=τln⁡H0τc (Tc−T0).t_a^* = \tau\ln\frac{H_0\tau}{c\,(T_c - T_0)}.

Proof. With no losses, c dT/dt=H0e−t/τc\,\mathrm dT/\mathrm dt = H_0 e^{-t/\tau}. Integrating from tat_a, where T=T0T = T_0, gives c (T−T0)=H0τ (e−ta/τ−e−t/τ)c\,(T - T_0) = H_0\tau\,(e^{-t_a/\tau} - e^{-t/\tau}), which tends to H0τe−ta/τH_0\tau e^{-t_a/\tau} as t→∞t \to \infty. Setting Tmax⁡=TcT_{\max} = T_c and solving for tat_a gives ta∗t_a^*; for an earlier accretion the peak is higher. ∎

A model of the parent body of the ordinary chondrites uses H0=1.9×10−7 W kg−1H_0 = 1.9 \times 10^{-7}\,\mathrm{W}\,\mathrm{kg}^{-1}, c=837 J kg−1 K−1c = 837\,\mathrm{J}\,\mathrm{kg}^{-1}\,\mathrm{K}^{-1} and T0=250 KT_0 = 250\,\mathrm{K}. Then H0τ/cH_0\tau/c is about 7400 K7400\,\mathrm{K}. Three temperatures matter: near 750 K750\,\mathrm{K} the porous rock starts to sinter into a compact one; near 1200 K1200\,\mathrm{K} the metal and the iron sulfide begin to melt together; near 1700 K1700\,\mathrm{K} the last silicates melt. A body that accreted before 1.7 Myr1.7\,\mathrm{Myr} melted entirely; before 2.1 Myr2.1\,\mathrm{Myr} its metal and sulfide melted and could sink into a core; until about 2.8 Myr2.8\,\mathrm{Myr} it was at least sintered and metamorphosed; later bodies stayed cool (figure below).

Heating by aluminium-26 of the inside of a chondritic body with no heat loss (the proposition above; H_0 = 1.9 × 10-7\, W\, kg-1, c = 837\, J\, kg-1\, K-1, T_0 = 250\, K, = 1.035\, Myr). Left: temperature against time for bodies accreted 1.5, 2, 2.5 and 3\, Myr after the CAIs. Right: peak temperature against accretion time; bodies that accreted before 1.7\, Myr melted, before 2.1\, Myr their metal and sulfide melted. Real bodies spend part of the heat on melting, which caps the temperature, and lose heat near their surface. Model values: Gail and others (2014).
Heating by aluminium-26 of the inside of a chondritic body with no heat loss (the proposition above; H0=1.9×10−7 W kg−1H_0 = 1.9 \times 10^{-7}\,\mathrm{W}\,\mathrm{kg}^{-1}, c=837 J kg−1 K−1c = 837\,\mathrm{J}\,\mathrm{kg}^{-1}\,\mathrm{K}^{-1}, T0=250 KT_0 = 250\,\mathrm{K}, τ=1.035 Myr\tau = 1.035\,\mathrm{Myr}). Left: temperature against time for bodies accreted 1.5, 2, 2.5 and 3 Myr3\,\mathrm{Myr} after the CAIs. Right: peak temperature against accretion time; bodies that accreted before 1.7 Myr1.7\,\mathrm{Myr} melted, before 2.1 Myr2.1\,\mathrm{Myr} their metal and sulfide melted. Real bodies spend part of the heat on melting, which caps the temperature, and lose heat near their surface. Model values: Gail and others (2014).

The proposition ignores heat loss, which holds only deep inside a large body. Heat diffuses over the distance κt\sqrt{\kappa t} in the time tt, where κ\kappa is the thermal diffusivity of the rock: with the values of the model, about 5 km5\,\mathrm{km} in one mean life of aluminium-26 (Exercise 2.12). Bodies much smaller than this lose their heat as fast as they make it; detailed models show that bodies more than 10 to 20 km10\text{ to }20\,\mathrm{km} across could melt.

Definition 2.16 (Planetary differentiation)

Planetary differentiation is the separation of a body, once molten or partly molten, into layers of different composition and density: a metal core, a silicate mantle and a crust made of the last melts.

The meteorites thus record two kinds of early bodies. Those that formed within about the first two million years had enough aluminium-26 to melt and differentiate: their cores are the iron meteorites, their crusts and mantles the achondrites, Vesta the survivor. Those that formed later stayed chondrites, heated from inside to a degree set by their date of birth and by the depth of each fragment: the deepest became type 6, the shallowest stayed type 3. The Earth itself, built from embryos and planetesimals of both kinds and struck by giant impacts, differentiated into its core, mantle and crust.

2.7 Exercises

Exercise 2.1 ★

In the minimum-mass solar nebula, compute the surface density of the gas at 1 au1\,\mathrm{au} and at 5.2 au5.2\,\mathrm{au}, and the surface density of the solids at 1 au1\,\mathrm{au} (rock only, Zrel=0.33Z_{\mathrm{rel}} = 0.33) and at 5.2 au5.2\,\mathrm{au} (all ices condensed, Zrel=1Z_{\mathrm{rel}} = 1). How many times heavier than the solids is the gas at 5.2 au5.2\,\mathrm{au}?

Solution

Solution of Exercise 2.1.

Gas: 2200 g cm−22200\,\mathrm{g}\,\mathrm{cm}^{-2} at 1 au1\,\mathrm{au}; at 5.2 au5.2\,\mathrm{au}, 2200×5.2−3/2=2200/11.86=186 g cm−22200\times5.2^{-3/2} = 2200/11.86 = 186\,\mathrm{g}\,\mathrm{cm}^{-2}. Solids: 33×0.33=10.9 g cm−233\times0.33 = 10.9\,\mathrm{g}\,\mathrm{cm}^{-2} at 1 au1\,\mathrm{au}; 33/11.86=2.78 g cm−233/11.86 = 2.78\,\mathrm{g}\,\mathrm{cm}^{-2} at 5.2 au5.2\,\mathrm{au}. At 5.2 au5.2\,\mathrm{au} the gas is 186/2.78=67186/2.78 = 67 times heavier than the solids.

Exercise 2.2 ★

Using the figure of the condensation temperatures, sort aluminium, sodium, magnesium, iron, sulfur and oxygen (as water ice) from the first to condense to the last, and give the volatility group of each. Which of them would you expect to find enriched in CAIs, and which depleted in the Earth compared with the Sun?

Solution

Solution of Exercise 2.2.

Aluminium (1653 K1653\,\mathrm{K}, ultrarefractory), magnesium (1336 K1336\,\mathrm{K}, common), iron (1334 K1334\,\mathrm{K}, common), sodium (958 K958\,\mathrm{K}, moderately volatile), sulfur (664 K664\,\mathrm{K}, volatile), oxygen as water ice (180 K180\,\mathrm{K}, highly volatile). CAIs, the first condensates, are enriched in aluminium (and calcium, titanium); the Earth is depleted in the volatile ones, sodium, sulfur and water, relative to magnesium and aluminium.

Exercise 2.3 ★

The catalogue lists 11071107 falls, of which 856 are ordinary chondrites and 49 irons, and 44 53444\,534 finds, of which 38 64738\,647 are ordinary chondrites, 10201020 irons and 22 09622\,096 come from south of 60∘60{}^{\circ} S. Compute the share of ordinary chondrites and of irons among the falls and among the finds, and the share of Antarctic finds. Which sample would you use to estimate the proportions in which the classes arrive on the Earth, and why?

Solution

Solution of Exercise 2.3.

Falls: ordinary chondrites 856/1107=77.3 %856/1107 = 77.3\,\%, irons 49/1107=4.4 %49/1107 = 4.4\,\%. Finds: 38647/44534=86.8 %38647/44534 = 86.8\,\% and 1020/44534=2.3 %1020/44534 = 2.3\,\%. Antarctic finds: 22096/44534=49.6 %22096/44534 = 49.6\,\%. The falls: every observed fall is collected whatever its class, whereas a find depends on how well the meteorite survives on the ground, how it looks and where people search, and half of the finds come from one kind of place.

Exercise 2.4 ★

The minerals of an object give the initial ratio X26X2226Al/X27X2227Al=1.0×10−5\ce{^{26}Al}/\ce{^{27}Al} = 1.0 \times 10^{-5}. How long after the CAIs did it form? Take R0=5.23×10−5R_0 = 5.23 \times 10^{-5} and a mean life of 1.035 Myr1.035\,\mathrm{Myr}.

Solution

Solution of Exercise 2.4.

Δt=1.035ln⁡(5.23×10−5/1.0×10−5)=1.035ln⁡5.23=1.035×1.654=1.71 Myr\Delta t = 1.035\ln(5.23 \times 10^{-5}/1.0 \times 10^{-5}) = 1.035\ln 5.23 = 1.035\times1.654 = 1.71\,\mathrm{Myr} after the CAIs.

Exercise 2.5 ★★

Using Proposition 2.2 with Σ1=2.2×104 kg m−2\Sigma_1 = 2.2 \times 10^{4}\,\mathrm{kg}\,\mathrm{m}^{-2} at 1 au1\,\mathrm{au} (1 au1\,\mathrm{au} =1.496×1011 m= 1.496 \times 10^{11}\,\mathrm{m}, M⊙=1.989×1030 kgM_\odot = 1.989 \times 10^{30}\,\mathrm{kg}), compute the mass of the minimum-mass nebula inside 30 au30\,\mathrm{au} and inside 100 au100\,\mathrm{au}, in solar masses. Compare with the observed discs, 0.001 to 0.10.001\text{ to }0.1 solar masses.

Solution

Solution of Exercise 2.5.

4πΣ1r12=4π×2.2×104×(1.496×1011)2=6.19×1027 kg4\pi\Sigma_1 r_1^2 = 4\pi\times2.2 \times 10^{4}\times(1.496 \times 10^{11})^2 = 6.19 \times 10^{27}\,\mathrm{kg}. Inside 30 au30\,\mathrm{au}: ×30=5.48\times\sqrt{30} = 5.48, so 3.39×1028 kg3.39 \times 10^{28}\,\mathrm{kg} =0.017= 0.017 solar masses; inside 100 au100\,\mathrm{au}: ×10\times10, 6.19×1028 kg6.19 \times 10^{28}\,\mathrm{kg} =0.031= 0.031 solar masses. Both lie within the observed range, nearer its upper end: the nebula was an ordinary disc.

Exercise 2.6 ★★

In the passive disc T=120 K (r/1 au)−3/7T = 120\,\mathrm{K}\,(r/1\,\mathrm{au})^{-3/7}, find where water ice (182 K182\,\mathrm{K}) and methane ice (41 K41\,\mathrm{K}) become stable. A grain in direct sunlight reaches 182 K182\,\mathrm{K} at 2.3 au2.3\,\mathrm{au}. Why do the two estimates of the snow line differ, and which one applies to the bodies that built the planets?

Solution

Solution of Exercise 2.6.

r=(120/T)7/3r = (120/T)^{7/3} au: for water, (120/182)7/3=0.6597/3=0.38 au(120/182)^{7/3} = 0.659^{7/3} = 0.38\,\mathrm{au}; for methane, (120/41)7/3=2.937/3=12.3 au(120/41)^{7/3} = 2.93^{7/3} = 12.3\,\mathrm{au}. The grain in direct sunlight is the surface layer of the disc; the midplane, shielded by the dust above it, is colder, so ice survives closer in. The planetesimals formed in the midplane, so the midplane estimate applies, provided the disc is passive; heating by accretion moves the snow line out.

Exercise 2.7 ★★

Compute the isolation mass, with b=5b = 5, of an embryo at 1 au1\,\mathrm{au} in the minimum-mass nebula (Σp=10.9 g cm−2\Sigma_p = 10.9\,\mathrm{g}\,\mathrm{cm}^{-2}), and at 5.2 au5.2\,\mathrm{au} beyond the snow line (Σp=2.78 g cm−2\Sigma_p = 2.78\,\mathrm{g}\,\mathrm{cm}^{-2}), in Earth masses (M⊕=5.972×1024 kgM_\oplus = 5.972 \times 10^{24}\,\mathrm{kg}, M⊙=1.989×1030 kgM_\odot = 1.989 \times 10^{30}\,\mathrm{kg}). What does it become at 5.2 au5.2\,\mathrm{au} in a disc with 3.53.5 times more solids?

Solution

Solution of Exercise 2.7.

At 1 au1\,\mathrm{au}: r2=2.238×1022r^2 = 2.238 \times 10^{22}, Σp=109 kg m−2\Sigma_p = 109\,\mathrm{kg}\,\mathrm{m}^{-2}, so 2πbr2Σp=31.42×2.238×1022×109=7.66×10252\pi b r^2\Sigma_p = 31.42\times2.238 \times 10^{22}\times109 = 7.66 \times 10^{25}; (7.66×1025)3/2=6.71×1038(7.66 \times 10^{25})^{3/2} = 6.71 \times 10^{38}, divided by (3×1.989×1030)1/2=2.443×1015(3\times1.989 \times 10^{30})^{1/2} = 2.443 \times 10^{15}: 2.75×1023 kg2.75 \times 10^{23}\,\mathrm{kg} =0.046= 0.046 Earth masses. At 5.2 au5.2\,\mathrm{au}: r2=6.05×1023r^2 = 6.05 \times 10^{23}, Σp=27.8 kg m−2\Sigma_p = 27.8\,\mathrm{kg}\,\mathrm{m}^{-2}, 2πbr2Σp=5.29×10262\pi b r^2\Sigma_p = 5.29 \times 10^{26}, M=1.22×1040/2.443×1015=4.98×1024 kg=0.83M = 1.22 \times 10^{40}/2.443 \times 10^{15} = 4.98 \times 10^{24}\,\mathrm{kg} = 0.83 Earth masses. With 3.53.5 times more solids, MM grows as 3.53/2=6.553.5^{3/2} = 6.55: 5.5 Earth masses.

Exercise 2.8 ★★

Rock of CI composition held 1.01×10161.01 \times 10^{16} atoms of aluminium-26 per gram when the CAIs formed, and each decay leaves 5.00×10−13 J5.00 \times 10^{-13}\,\mathrm{J} of heat. Compute the heat released per kilogram over the whole decay, and the temperature rise it would give with c=837 J kg−1 K−1c = 837\,\mathrm{J}\,\mathrm{kg}^{-1}\,\mathrm{K}^{-1} and no losses. How many times more is it than the heat needed to warm the rock from 250250 to 1700 K1700\,\mathrm{K}?

Solution

Solution of Exercise 2.8.

Per kilogram, 1.01×10191.01 \times 10^{19} atoms, so 1.01×1019×5.00×10−13=5.05×106 J kg−11.01 \times 10^{19}\times5.00 \times 10^{-13} = 5.05 \times 10^{6}\,\mathrm{J}\,\mathrm{kg}^{-1}. The temperature rise would be 5.05×106/837=6030 K5.05 \times 10^{6}/837 = 6030\,\mathrm{K}. Warming from 250250 to 1700 K1700\,\mathrm{K} takes 837×1450=1.21×106 J kg−1837\times1450 = 1.21 \times 10^{6}\,\mathrm{J}\,\mathrm{kg}^{-1}: aluminium-26 supplied about four times that.

Exercise 2.9 ★★

Over its whole decay, aluminium-26 releases 5.07×1010 erg5.07 \times 10^{10}\,\mathrm{erg} per gram of CI rock and iron-60 7.49×107 erg7.49 \times 10^{7}\,\mathrm{erg} per gram. Compute the ratio of the two energies. The half-life of iron-60 is 2.62 Myr2.62\,\mathrm{Myr}: compute its mean life, and the ratio of the two initial heating powers, energy divided by mean life.

Solution

Solution of Exercise 2.9.

5.07×1010/7.49×107=6775.07 \times 10^{10}/7.49 \times 10^{7} = 677. The mean life of iron-60 is 2.62/ln⁡2=3.78 Myr2.62/\ln2 = 3.78\,\mathrm{Myr}. Initial powers per gram: 5.07×1010/1.035=4.90×10105.07 \times 10^{10}/1.035 = 4.90 \times 10^{10} and 7.49×107/3.78=1.98×1077.49 \times 10^{7}/3.78 = 1.98 \times 10^{7} erg per gram per million years, a ratio of about 24702470: aluminium-26 dominated by far.

Exercise 2.10 ★★★

With the results of Exercise 2.9, find the time after which iron-60 heats rock faster than aluminium-26. Could iron-60 have melted a body that accreted that late? (Compare the energy it still has to release then with the 1.2×106 J kg−11.2 \times 10^{6}\,\mathrm{J}\,\mathrm{kg}^{-1} needed to warm rock from 250250 to 1700 K1700\,\mathrm{K}.)

Solution

Solution of Exercise 2.10.

The powers are equal when 4.90×1010 e−t/1.035=1.98×107 e−t/3.784.90 \times 10^{10}\,e^{-t/1.035} = 1.98 \times 10^{7}\,e^{-t/3.78}, that is t(1/1.035−1/3.78)=ln⁡(4.90×1010/1.98×107)=ln⁡2475=7.81t(1/1.035 - 1/3.78) = \ln(4.90 \times 10^{10}/1.98 \times 10^{7}) = \ln 2475 = 7.81; with 0.966−0.265=0.7020.966 - 0.265 = 0.702 per million years, t=11.1 Myrt = 11.1\,\mathrm{Myr}. By then iron-60 has left 7.49×107 e−11.1/3.78=7.49×107×0.053=4.0×106 erg g−17.49 \times 10^{7}\,e^{-11.1/3.78} = 7.49 \times 10^{7}\times0.053 = 4.0 \times 10^{6}\,\mathrm{erg}\,\mathrm{g}^{-1}, that is 400 J kg−1400\,\mathrm{J}\,\mathrm{kg}^{-1}: thousands of times too little to melt anything.

Exercise 2.11 ★★★

With H0=1.9×10−7 W kg−1H_0 = 1.9 \times 10^{-7}\,\mathrm{W}\,\mathrm{kg}^{-1}, τ=1.035 Myr\tau = 1.035\,\mathrm{Myr}, c=837 J kg−1 K−1c = 837\,\mathrm{J}\,\mathrm{kg}^{-1}\,\mathrm{K}^{-1} and T0=250 KT_0 = 250\,\mathrm{K} (a year of 3.156×107 s3.156 \times 10^{7}\,\mathrm{s}), compute the latest accretion times for which the centre of a large body reaches 1200 K1200\,\mathrm{K} and 1700 K1700\,\mathrm{K}. Some measurements suggest that the initial aluminium-26 was only half the canonical value in parts of the disc: by how much would both times move? Show that the shift does not depend on the critical temperature.

Solution

Solution of Exercise 2.11.

τ=1.035×3.156×1013=3.267×1013 s\tau = 1.035\times3.156 \times 10^{13} = 3.267 \times 10^{13}\,\mathrm{s} and H0τ=6.21×106 J kg−1H_0\tau = 6.21 \times 10^{6}\,\mathrm{J}\,\mathrm{kg}^{-1}. For 1200 K1200\,\mathrm{K}: ta∗=1.035ln⁡(6.21×106/(837×950))=1.035ln⁡7.81=2.13 Myrt_a^* = 1.035\ln\bigl(6.21 \times 10^{6}/(837\times950)\bigr) = 1.035\ln7.81 = 2.13\,\mathrm{Myr}; for 1700 K1700\,\mathrm{K}: 1.035ln⁡(6.21×106/(837×1450))=1.035ln⁡5.12=1.69 Myr1.035\ln\bigl(6.21 \times 10^{6}/(837\times1450)\bigr) = 1.035\ln5.12 = 1.69\,\mathrm{Myr}. Halving H0H_0 subtracts τln⁡2=0.72 Myr\tau\ln2 = 0.72\,\mathrm{Myr} from both, since ta∗=τln⁡H0+τln⁡(τ/c(Tc−T0))t_a^* = \tau\ln H_0 + \tau\ln(\tau/c(T_c - T_0)) and only the first term changes: the shift is one half-life, whatever the critical temperature.

Exercise 2.12 ★★★

Rock of conductivity 2.1 W m−1 K−12.1\,\mathrm{W}\,\mathrm{m}^{-1}\,\mathrm{K}^{-1}, density 3300 kg m−33300\,\mathrm{kg}\,\mathrm{m}^{-3} and specific heat capacity 837 J kg−1 K−1837\,\mathrm{J}\,\mathrm{kg}^{-1}\,\mathrm{K}^{-1} has the thermal diffusivity κ=k/(ρc)\kappa = k/(\rho c). Compute it, and the distance κτ\sqrt{\kappa\tau} over which heat diffuses in one mean life of aluminium-26. Explain why the heating proposition applies to the centre of a body 100 km100\,\mathrm{km} across but not to a body 2 km2\,\mathrm{km} across, and compare with the smallest bodies that could melt, 10 to 20 km10\text{ to }20\,\mathrm{km} across.

Solution

Solution of Exercise 2.12.

κ=2.1/(3300×837)=7.60×10−7 m2 s−1\kappa = 2.1/(3300\times837) = 7.60 \times 10^{-7}\,\mathrm{m}^{2}\,\mathrm{s}^{-1} and κτ=7.60×10−7×3.267×1013=4980 m\sqrt{\kappa\tau} = \sqrt{7.60 \times 10^{-7}\times3.267 \times 10^{13}} = 4980\,\mathrm{m}, about 5 km5\,\mathrm{km}. The centre of a body 100 km100\,\mathrm{km} across is 50 km50\,\mathrm{km} deep, ten times this distance: its heat cannot escape while aluminium-26 decays, and the proposition applies. A body 2 km2\,\mathrm{km} across loses its heat as fast as it makes it. Melting needs a radius of a few times 5 km5\,\mathrm{km}, that is a diameter of some 10 to 20 km10\text{ to }20\,\mathrm{km} or more, as the models find.

2.8 Problem: The stone that fell on the frozen lake

Problem 2.1

Weekend problem — classifying a freshly fallen stone, dating its chondrules with aluminium-26, and finding how hot its parent body became

On a winter night a fireball crosses the sky above an invented northern lake, and a boom is heard. Next morning a black-crusted stone of 2.4 kg2.4\,\mathrm{kg} is picked up from the ice, in the middle of a fresh crater in the snow. The cut stone shows chondrules making 75 %75\,\% of its volume, 0.3 mm0.3\,\mathrm{mm} across on average, with blurred edges; metal makes 9 %9\,\% of its volume, and its olivine holds 18 mole per cent of fayalite. A laboratory analyses one chondrule (table below). Use the tables of the chapter, R0=5.23×10−5R_0 = 5.23 \times 10^{-5}, τ=1.035 Myr\tau = 1.035\,\mathrm{Myr}, a year of 3.156×107 s3.156 \times 10^{7}\,\mathrm{s}, H0=1.9×10−7 W kg−1H_0 = 1.9 \times 10^{-7}\,\mathrm{W}\,\mathrm{kg}^{-1}, c=837 J kg−1 K−1c = 837\,\mathrm{J}\,\mathrm{kg}^{-1}\,\mathrm{K}^{-1}, T0=250 KT_0 = 250\,\mathrm{K}, a conductivity of 2.1 W m−1 K−12.1\,\mathrm{W}\,\mathrm{m}^{-1}\,\mathrm{K}^{-1}, a density of 3300 kg m−33300\,\mathrm{kg}\,\mathrm{m}^{-3}, 1 au1\,\mathrm{au} =1.496×1011 m= 1.496 \times 10^{11}\,\mathrm{m}, M⊙=1.989×1030 kgM_\odot = 1.989 \times 10^{30}\,\mathrm{kg} and M⊕=5.972×1024 kgM_\oplus = 5.972 \times 10^{24}\,\mathrm{kg}.

mineralX27X2227Al/X24X2224Mg\ce{^{27}Al}/\ce{^{24}Mg}X26X2226Mg/X24X2224Mg\ce{^{26}Mg}/\ce{^{24}Mg}
olivine0.050.139 5000.139\,500
glass3.000.139 5180.139\,518
plagioclase15.00.139 5900.139\,590

Part I — The stone.

  1. Is the stone a fall or a find? Why does that make it valuable?
  2. Of the 11071107 falls of the catalogue, 856 are ordinary chondrites. What is the chance that a fall, taken at random, is one?
  3. Using the table of chondrite groups, say to which group the stone belongs.
  4. Its chondrules have blurred edges and its minerals are of uniform composition. Which petrologic types does that suggest, and what is the highest temperature its group reached, in kelvins and in degrees Celsius?
  5. Why would you expect few CAIs in it, and why is the stone nonetheless younger than any CAI?

Part II — When did the chondrule form?

  1. Explain why the ratio X26X2226Mg/X24X2224Mg\ce{^{26}Mg}/\ce{^{24}Mg} rises with X27X2227Al/X24X2224Mg\ce{^{27}Al}/\ce{^{24}Mg} from mineral to mineral.
  2. From the olivine and the plagioclase, compute the slope of the isochron. What does it measure?
  3. Check that the glass lies on the same line.
  4. How long after the CAIs did the chondrule form?
  5. Taking the CAIs to be 4567.5 Ma4567.5\,\mathrm{Ma} old, give the age of the chondrule.
  6. What would the ratio X26X2226Al/X27X2227Al\ce{^{26}Al}/\ce{^{27}Al} have been in a mineral formed 10 Myr10\,\mathrm{Myr} after the CAIs? Which nuclide of the table would date such an event better?

Part III — Heating the parent body. The stone’s parent body, 20 km20\,\mathrm{km} in radius, assembled 2.4 Myr2.4\,\mathrm{Myr} after the CAIs, soon after its chondrules.

  1. Compute the thermal diffusivity of the rock and the time R2/κR^2/\kappa for heat to escape from the centre, in millions of years. May the centre be treated as losing no heat?
  2. Compute the heating power of aluminium-26 per kilogram of rock at the time of assembly.
  3. Compute the heat still to be released per kilogram after assembly.
  4. Compute the temperature rise and the peak temperature at the centre.
  5. Compare with the three critical temperatures of the chapter. Did the centre melt? Is this consistent with question 4?
  6. What is the latest assembly time for which the metal and sulfide of the centre would have begun to melt?
  7. How long after assembly is 90 %90\,\% of the temperature rise reached? Give the date after the CAIs.

Part IV — The neighbours.

  1. An achondrite comes from a body that assembled 0.8 Myr0.8\,\mathrm{Myr} after the CAIs. Compute the peak temperature that the proposition would give, and say what really happened to the body.
  2. How many times more aluminium-26 per kilogram did that body have at its assembly than the stone’s parent body?
  3. The stone is poorer than the Sun in zinc and lead, but not in magnesium or aluminium. Use the condensation temperatures to explain why.
  4. In the minimum-mass nebula, compute the surface density of rock (Zrel=0.33Z_{\mathrm{rel}} = 0.33) at 2.5 au2.5\,\mathrm{au}, where the parent body orbited.
  5. Compute the isolation mass at 2.5 au2.5\,\mathrm{au}, with b=5b = 5, in kilograms and in Earth masses, and compare it with the mass of the parent body.
  6. State the peak temperature of the centre of the parent body, the key number of this problem, in degrees Celsius.
Solution

Solution of Problem 2.1.

1. A fall: its arrival was seen and it was picked up the next morning. It is unweathered, and it adds to the unbiased sample of what reaches the Earth.

2. 856/1107=0.77856/1107 = 0.77: about three chances in four.

3. Chondrules 0.3 mm0.3\,\mathrm{mm} across, 9 %9\,\% metal and olivine with 18 %18\,\% fayalite: an H chondrite (10 %10\,\% metal, Fa 16–20).

4. Blurred chondrules and uniform minerals mean thermal metamorphism: types 5 to 6. The ordinary chondrites reached at most about 1220 K1220\,\mathrm{K}, that is 947∘C947{}^{\circ}\mathrm{C}.

5. Ordinary chondrites hold only 0.1 to 1 %0.1\text{ to }1\,\% CAIs by volume. The stone is made of chondrules and matrix assembled after the CAIs, which are the oldest solids: any rock that contains them is younger.

6. Each mineral received at formation the same magnesium plus aluminium-26 in proportion to its aluminium; the aluminium-26 has since turned into magnesium-26, so the more aluminium a mineral had per magnesium, the more magnesium-26 it gained.

7. (0.139590−0.139500)/(15.0−0.05)=0.000 090/14.95=6.02×10−6(0.139590 - 0.139500)/(15.0 - 0.05) = 0.000\,090/14.95 = 6.02 \times 10^{-6}: the ratio X26X2226Al/X27X2227Al\ce{^{26}Al}/\ce{^{27}Al} of the chondrule when it formed.

8. 0.139500+6.02×10−6×3.00=0.1395180.139500 + 6.02 \times 10^{-6}\times3.00 = 0.139518: on the line.

9. Δt=1.035ln⁡(5.23×10−5/6.02×10−6)=1.035ln⁡8.69=1.035×2.162=2.24 Myr\Delta t = 1.035\ln(5.23 \times 10^{-5}/6.02 \times 10^{-6}) = 1.035\ln8.69 = 1.035\times2.162 = 2.24\,\mathrm{Myr}.

10. 4567.5−2.2=4565.3 Ma4567.5 - 2.2 = 4565.3\,\mathrm{Ma}.

11. 5.23×10−5 e−10/1.035=5.23×10−5×6.4×10−5=3.3×10−95.23 \times 10^{-5}\,e^{-10/1.035} = 5.23 \times 10^{-5}\times6.4 \times 10^{-5} = 3.3 \times 10^{-9}, far too little magnesium-26 to measure. Hafnium-182 (half-life 8.9 Myr8.9\,\mathrm{Myr}) or manganese-53 (3.7 Myr3.7\,\mathrm{Myr}) suit events some ten million years after the CAIs.

12. κ=2.1/(3300×837)=7.60×10−7 m2 s−1\kappa = 2.1/(3300\times837) = 7.60 \times 10^{-7}\,\mathrm{m}^{2}\,\mathrm{s}^{-1}; R2/κ=4×108/7.60×10−7=5.26×1014 s=16.7 MyrR^2/\kappa = 4 \times 10^{8}/7.60 \times 10^{-7} = 5.26 \times 10^{14}\,\mathrm{s} = 16.7\,\mathrm{Myr}, sixteen times the mean life of aluminium-26: the centre keeps its heat.

13. 1.9×10−7 e−2.4/1.035=1.9×10−7×0.0984=1.87×10−8 W kg−11.9 \times 10^{-7}\,e^{-2.4/1.035} = 1.9 \times 10^{-7}\times0.0984 = 1.87 \times 10^{-8}\,\mathrm{W}\,\mathrm{kg}^{-1}.

14. 1.87×10−8×3.267×1013=6.11×105 J kg−11.87 \times 10^{-8}\times3.267 \times 10^{13} = 6.11 \times 10^{5}\,\mathrm{J}\,\mathrm{kg}^{-1}.

15. ΔT=6.11×105/837=730 K\Delta T = 6.11 \times 10^{5}/837 = 730\,\mathrm{K}; peak 250+730=980 K250 + 730 = 980\,\mathrm{K}.

16. Above 750 K750\,\mathrm{K}, so the rock sintered and was metamorphosed; below 1200 K1200\,\mathrm{K}, so neither metal nor sulfide melted. The stone is a heated chondrite, of high petrologic type, below the 1220 K1220\,\mathrm{K} of question 4: consistent.

17. 1.035ln⁡(6.21×106/(837×950))=2.13 Myr1.035\ln\bigl(6.21 \times 10^{6}/(837\times950)\bigr) = 2.13\,\mathrm{Myr}: about 0.3 Myr0.3\,\mathrm{Myr} earlier.

18. The rise is ΔT(1−e−(t−ta)/τ)\Delta T(1 - e^{-(t - t_a)/\tau}); 90 %90\,\% needs e−(t−ta)/τ=0.1e^{-(t - t_a)/\tau} = 0.1, t−ta=1.035ln⁡10=2.38 Myrt - t_a = 1.035\ln10 = 2.38\,\mathrm{Myr}, that is 4.8 Myr4.8\,\mathrm{Myr} after the CAIs.

19. ΔT=7420 e−0.8/1.035=7420×0.462=3430 K\Delta T = 7420\,e^{-0.8/1.035} = 7420\times0.462 = 3430\,\mathrm{K}, peak 3680 K3680\,\mathrm{K}, far above the 1700 K1700\,\mathrm{K} at which the silicates are molten. The body melted (the heat of melting held its temperature down) and differentiated into a metal core and a silicate mantle and crust.

20. e(2.4−0.8)/1.035=e1.546=4.7e^{(2.4 - 0.8)/1.035} = e^{1.546} = 4.7 times more.

21. Zinc and lead condense near 727 K727\,\mathrm{K}, moderately volatile; magnesium (1336 K1336\,\mathrm{K}) and aluminium (1653 K1653\,\mathrm{K}) condense far earlier. The solids of the stone’s region gathered before the gas had cooled enough for all the zinc and lead to condense.

22. 33×0.33×2.5−3/2=10.89/3.95=2.76 g cm−233\times0.33\times2.5^{-3/2} = 10.89/3.95 = 2.76\,\mathrm{g}\,\mathrm{cm}^{-2}, 27.6 kg m−227.6\,\mathrm{kg}\,\mathrm{m}^{-2}.

23. r2=(3.74×1011)2=1.40×1023r^2 = (3.74\times10^{11})^2 = 1.40 \times 10^{23}; 2π×5×1.40×1023×27.6=1.21×10262\pi\times5\times1.40 \times 10^{23}\times27.6 = 1.21 \times 10^{26}; (1.21×1026)3/2/2.443×1015=5.45×1023 kg(1.21 \times 10^{26})^{3/2}/2.443 \times 10^{15} = 5.45 \times 10^{23}\,\mathrm{kg}, 0.0910.091 Earth masses. The parent body, 43π(2×104)3×3300=1.1×1017 kg\tfrac43\pi(2 \times 10^{4})^3\times3300 = 1.1 \times 10^{17}\,\mathrm{kg}, is five million times lighter: a planetesimal, not an embryo.

24. The centre of the parent body peaked at about 707∘C707{}^{\circ}\mathrm{C} (980−273980 - 273): hot enough to recrystallise the chondrite, too cool to melt its metal.

Terms defined in this chapter

See all 957 terms in the glossary