Isotope Geochemistry
Isotope Geochemistry
Isotope Geochemistry
◊ISOTOPES
DEFINITION
Two or more nuclides having the same atomic number, thus constituting
the same element, but differing in the mass number. Isotopes of a given
element have the same number of nuclear protons but differing numbers of
neutrons.
TYPES:
STABLE ISOTOPES
The atomic nuclei of these elements do not change to nuclei of other
elements.
RADIOACTIVE ISOTOPES
The atomic nuclei of these elements give out radiation spontaneously and
thereby change to nuclei of other elements.
◊ RADIOACTIVE ISOTOPES
►NUCLIDE
Stable nuclides - maintain their atomic configuration over long periods of time.
Unstable nuclides - spontaneously change of an unstable nuclide into another
nuclide.
►RADIOACTIVE DECAY
□Radioactive Isotopes
Parent nuclide (unstable) – before decay the atom containing the radioactive
nuclide
1. Alpha decay
Alpha emission results in releasing an alpha particle. An alpha particle has two
protons and two neutrons, so it has a positive charge. (Since it has two protons it
is a helium nucleus.) It is written in equations like this:
2. Beta decay
Beta emission is when a high speed electron (negative charge) leaves the nucleus.
Beta emission occurs in elements with more neutons than protons, so a neutron splits
into a proton and an electron. The proton stays in the nucleus and the electron is
emitted. Negative electrons are represented as follows:
3. Gamma Emission
Gamma Emission is when an excited nucleus gives off a ray in the gamma
part of the spectrum. A gamma ray has no mass and no charge. This often
occurs in radioactive elements because the other types of emission can result in
an excited nucleus. Gamma rays are represented with the following symbol.
The two types of artificial radiation are positron emission and electron capture.
Positron emission
Positron emission involves a particle that has the same mass as an electron but a
positive charge. The particle is released from the nucleus.
Electron capture
Electron capture is when an unstable nucleus grabs an electron from its inner shell to
help stabilize the nucleus. The electrons combine with a proton to form a neutron which
stays in the nucleus.
□Principle of Radioactive Decay
The graph also shows the half-life concept. The half-life is the amount of time
necessary to reduce the number of parent atoms by 50% from the original
number.
► The Basic Equation of Radioactive Decay
The number of decays you will measure each second from a sample depends
on the number of atoms in the sample, N.
2 kg 1 kg
Here are two blocks of exactly the same radioisotope. The chance of an atom
decaying from one is exactly the same as in the other but there are twice as
many atoms in the 2 kg block so there will be twice as many decays per second
in the 2 kg block.
where
λ = the constant of proportionality, called the Decay Constant.
The decay constant is the proportion of atoms that decay in an interval of time
The decay constant gives you an idea of how quickly or slowly a material will
decay.
A large λ value means that the sample will decay more quickly.
Rate of decay of a radioactive nuclide is proportional to the number of atoms
of that nuclide remaining at any time.
- dN/dt = λN (1)
where
λ is the proportionality constant known as the decay constant
N is the number of atoms remaining/present
and the minus sign indicates that the rate of decay decreases with time.
Rearrangement
- dN/N = λdt (2)
-In N = λt + C (4)
When t=0,
No= number of nuclides at t=0
-In No = (0) + C
C = - In No
In N = λt - In No (5)
Rearrangement
In N - In No = -λt (6)
Present Original
In (N/No) = -λt (7)
The equation above is the basic relationship that describes all radioactive decay
processes.
With it, we can calculate the number of parent atoms (N) that remain at any time t
from the original number of atoms (No) present at time t=0.
Rearrangement
Present Original
N = No e-λt
(9)
The graph plots the number of radioactive nuclei at any time, Nt, against time,
t. We can see that the number of radioactive nuclei decreases from N0 that is
the number at t = 0 in a rapid fashion initially and then more slowly in the classic
exponential manner.
All three curves here are exponential in nature, only the Decay Constant is
different.
• When the Decay Constant has a low value the curve decreases relatively
slowly
• When the Decay Constant is large the curve decreases very quickly.
Original Present
No = N e λt (10)
Note:
N = N0e-kt (exponential decay)
where
• N0 is the initial quantity
• t is time
• N(t) is the quantity after time t
• k is the decay constant and
• ex is the exponential function (e is the base of the natural logarithm)
www.earth.northwestern.edu/people/seth/202/DECAY/decay.pennies.slow.html
Decay of parent produces daughter or radiogenic nuclides.
original present
D = No – N (11)
Since in general there will be some atoms of the daughter nuclide around to
begin with, i.e. when t = 0, a more general expression is:
D = Do + N (e λt
– 1) (13)
Where Do is the number of daughters originally present.
Rearrangement
This is the time during which an amount of the daughter represented by D has
accumulated, leaving undecayed an amount of the parent represented by P.
Values of D and P are found by analyzing the rock or mineral in which the
radioactive isotope occurs.
If we can also find values for λ and Do the equation will give us the age of the
rock or mineral in years.
A half-life is the time it takes for half of the parent radioactive element to
decay to a daughter product.
So if you have 10 grams of a radioactive element
After one half-life there will be 5 grams of the radioactive element left.
After another half-life, there will be 2.5 g of the original element left.
After another half-life, 1.25 g will be left.
Radioactive decay occurs at a constant exponential or geometric rate.
The rate of decay is proportional to the number of parent atoms present.
Approach
o Compare amount of daughter isotope to amount of parent
originally there
Example:
Problem: The 235U: 207Pb ratio in a mineral is 1:7.
What is the age of the mineral?
Given: Half-life of 235U is 0.7 billion years (b. y.)
The half-life of an isotope equals the number of years it takes for an initial
number of parent atoms to be reduced to half that number by radioactive decay.
The half-life figure enables us to relatively quickly understand the useful age
range of a particular isotopic system.
For instance, the half life of the C-14 system is 5,730 years - you would never
use C-14 to determine the age of material older than 40 000 years which is the
practical upper limit; all of the radioactivity would be gone.
XXXXXXXXXX
□Half-life Equations:
1) Suppose the amount of time for the number of parent atoms to decrease to
half the original number i.e. t when N/No =1/2 is required to be determined.
to finally get
t1/2 = In2/λ (6)
An example problem is if you originally had 157 grams of carbon-14 and the half-
life of carbon-14 is 5730 years, how much would there be after 2000 years?
Once this rate is known, the length of time over which decay has been occurring
can be estimated by measuring the amount of radioactive parent element
and the amount of stable daughter elements.
The course examines K-Ar, U-Th-Pb, Rb-Sr decay systems and Carbon-14.
►RUBIDIUM-STRONTIUM
λ= 1.42x10-11 y-1.
□ Background
Rubidium: Strontium:
• Univalent • Divalent
• Not very common in the Earth's crust • Occurs as four stable
isotopes
(88Sr, 87Sr, 86Sr and 84Sr).
The table below lists the naturally occurring isotopes of both Rb and Sr along
with their isotopic abundances (in atom %) and their nuclide weights in atomic
mass units (a. m. u.).
iv) Of these terms, (87Sr/86Sr)t , which is the total 87Sr/86Sr, is measured in the
laboratory; 87Rb/86Sr is calculated from the measured Rb and Sr
concentrations in the sample; and (87Sr/86Sr)o and t are unknowns.
However, if neither the initial ratio nor the age of the sample is known,
then neither can be computed using the equation above.
If the body of rock under study contains rocks with different Rb/Sr
ratios and the rocks are known, based on geological observations, to
have formed at the same time from the same source, then an equation
like the one above can be written for each sample.
• If only two samples are available, the two equations may be solved
simultaneously to give both the initial ratio and the age of the samples.
• If more than two samples are available, then all of the equations are
solved simultaneously using least squares methods to give best fit
values for the initial ratio and the age of the samples. The latter approach is
preferred and is called the Isochron Method.
□ Methods
Method #1: Direct comparison
(Since no 87Rb in this sample all 87Sr must have been present to start with-- it
is not radiogenic).
In this method minerals with varying amounts of Rb are analyzed that are the
same age.
• At time of crystallization 87Sr/86Sr ratio is the same for all minerals of the same
rock.
• The amount of 87Sr that you measure is equal to the original amount PLUS
what has been generated by radioactive decay of rubidium.
• All four of these samples have the same initial 87Sr/86Sr ratio (shown by the
black dashed line) but different 87Rb content so different 87Rb/86Sr ratios.
• With time, some of the 87Rb in the samples decays to 87Sr. The red arrows
show how the locations of the samples move as a function of time (note that
one Sr is produced by each Rb that decays).
• The 87Rb decreases while the 87Sr increases. As Rb decays to form Sr and
the samples evolve, they remain colinear.
• You can think of the horizontal line originally defined by the initial ratio of the
samples rotating with its fixed point located at the initial 87Sr/86Sr ratio and
an 87Rb/86Sr value of zero.
Now consider the form of the Sr isotope evolution equation from above:
(λ = decay constant)
• For the variables in the diagram above, this equation is the equation of
a straight line (y = mx + b), where y = (87Sr/86Sr)t, x = (87Rb/86Sr), b =
(87Sr/86Sr)o and the slope of the line (m) is e λt -1.
• A straight line is then fit to the data using linear regression (most spread
sheets and hand calculators have linear regression functions).
• The intercept of the best fit line gives the initial ratio [(87Sr/86Sr)o] for the
sample suite.
The use of this method is based on the validity of the following assumptions:
□ Principle
• Potassium (K) is one of the most abundant elements in the Earth's crust
(2.4% by mass).
• 11 become Ar-40.
• The other 89% of the 40K atoms decay to 40Ca.
• For the K/Ar dating system, this decay scheme to calcium isotopes is
ignored.
40
• The buildup of radiogenic Ar (40Ar*) in a closed system can be
expressed by the equation:
• To check for this it is necessary to analyze also for 36Ar which has
a known ratio to 40Ar in the atmosphere (40Ar/36Ar = 296).
• As the K-40 in the rock decays into Ar-40, the gas is trapped in the
rock.
►URANIUM-THORIUM-LEAD
Since 238U is much more abundant than 235U, let’s look at the decay
sequence for 238U which is given below.
(235U and 232Th have their unique decay sequences as well).
The ages should agree. If they do, they are called concordant; if not,
they are discordant.
• Discordant ages indicate that there has been some type of
disturbance which has reset the "time clock" at some time in the
past, such that the steady decrease in U has been interrupted
and reset (e.g., metamorphism).
A good data set will usually produce a linear plot which may be
extrapolated back to intersect the CONCORDIA curve.
Use
This dating method is used for zircon, uranite, pitchblende - any
uranium or thorium bearing mineral.
►CARBON-14
b. (i) Living things are in equilibrium with the atmosphere, and the
radioactive carbon dioxide is absorbed and used by plants.
(ii) The radioactive carbon dioxide gets into the food chain and the
carbon cycle.
Progressively through time, the carbon-14 atoms decay and once again become
nitrogen-14.
The half-life is so short (5730 years) that this method can only be practically used
on materials less than 40,000 years old.
Beyond 40-50,000 years, there usually is not enough carbon-14 left to measure
with conventional laboratory methods.
Half-Lives
Years Past C-14 Atoms C-12 Atoms
0 0 1N 1N
1 5,730 1/2 N 1N
2 11,460 1/4 N 1N
3 17,190 1/8 N 1N
4 22,920 1/16 N 1N
5 28,650 1/32 N 1N
6 34,380 1/64 N 1N
7 40,110 1/128 N 1N
Some other things to realize, think about, etc. about "absolute" age
determination:
• Even the best instruments carry an instrumental error of at least 0.2 -
2%, meaning that on a 100 Ma date: range is between 98-102 Ma.
• Records crystallization time (detrital zircon mineral in a sedimentary
rock records date of rock from which the zircon was eroded)
• Contamination problems abound (e.g., volcanic ash flow may pick up
older/detrital feldspar or biotite minerals on its excursion over a rock.
• Metamorphism may "reset" clock
• Loss of either parent or daughter rocks due to heat or weathering. Need
FRESH samples
• Have to have enough of an isotope to measure (low crustal abundance
of some elements, or low abundance in rocks being studied).
• Half-life has to be short enough to have produced enough daughter
isotopes for analysis (limits use of "long-term" isotopes such Rb-Sr, U-
Pb, Th-Pb to rocks over a 100Ma old).
• Half-life can't be so short such that the original parent is beyond
detection.
◊OTHER USES OF RADIOACTIVE ISOTOPES
► GEOCHEMICAL PROCESSES
□Introduction
• The UR model is indicated. For this model, the 87Sr/86Sr ratio of the Earth evolves from
0.69897 (BABI) 4.6 Ga ago to 0.7045 today.
□Generation of Crust
Now consider what happens when the Uniform Reservoir is partially melted to generate
melt, which becomes crust.
• During a partial melting event, the mantle that is partially melted, the
partial melt and the mantle that is left behind (residual mantle) all have
the same 87Sr/86Sr ratio (shown by the dot in the figure).
• However, partial melting does fractionate Rb from Sr in such a way that
the partial melt (= crust) has a higher Rb/Sr ratio than the original
mantle that was melted and the residual mantle has a lower Rb/Sr ratio
than the original mantle.
• Thus, the 87Sr/86Sr ratio of the partial melt (= crust) grows more quickly
than that of the original (UR) mantle (melt + crust line in the figure) and
the 87Sr/86Sr ratio of the residual (depleted) mantle grows more slowly
than the UR mantle (residual mantle line in the figure).
• The result of this type of fractionation event is that crustal rocks, in
general, have higher 87Sr/86Sr ratios than the UR model and depleted
mantle rocks (from which partial melts were derived) have lower
87
Sr/86Sr ratios than the UR model.
□Melting of Crust
This type of modelling can be carried further by considering the isotopic
composition of a partial melt of the crust (i.e., a granite magma).
• During melting of the crust to form a granite magma, the magma will
again have a higher Rb/Sr ratio than the crust that is being melted and
the residual (unmelted) crust will have a lower Rb/Sr ratio than the
original crust.
• Thus, over time the (rock) granite's 87Sr/86Sr ratio will evolve to higher
values than the original crust that was melted and the 87Sr/86Sr of the
residual crust in the granite source region will evolve more slowly than
that of the original crust.
• Consequently, granites are expected to have quite high 87Sr/86Sr ratios (as
high as 0.715) compared to the UR model.
In both of these cases, the initial ratios of the basalt and the granite (the
87
Sr/86Sr composition of the melts at the time they were formed) will reflect
the isotopic compositions of the source regions of those melts. These initial
ratios are very useful to the isotope geochemistry.