LectureNote Basics IsotopeHydro BNOC 21sept14
LectureNote Basics IsotopeHydro BNOC 21sept14
------------------------------------------------------------------------------------
---------------------------------------------------------------------------------------------------------------------------------------
Nasir Ahmed
CE and Head, Isotope Hydrology Division
Institute of Nuclear Science & Technology (INST)
Atomic Energy Research Establishment (AERE)
Ganakbari, Savar, Dhaka
Page 2 of 14
Environmental tracers
A tracer is a substance, which is present only in trace concentrations and marks a trace in a
natural system (marker). In hydrology, it ideally marks the water itself or dissolved
substances therein, and moves passively with the water or the solutes.
Isotopes are often nearly ideal tracers, e.g.:
- isotopes of H and O mark the water molecule
- isotopes of C, N, S etc. mark dissolved substances
- Conservative solutes (e.g. noble gases) are good tracers
Environmental tracers are isotopes or other tracers that are widely distributed in the
environment, especially the hydrosphere. They can be both of natural or anthropogenic
origin. They can be further classified into
The term "environmental tracers" is mainly used to distinguish them from "artificial tracers",
both of which are often simply called tracers but are used in quite different settings and by
largely different researchers. Artificial tracers are released deliberately into the investigated
system, in local tracer experiments (e.g., dyes, chemicals, also isotopes, SF6 and so on).
The essential difference is that environmental tracers can be applied for studies on
large spatial and temporal scales, whereas artificial are suitable for targeted studies
of small scale processes.
2. Introduction to isotopes
Nuclei of atoms consist of protons and neutrons (= nucleons)
Z = number of protons (atomic number)
N = number of neutrons
A = Z + N = number of nucleons (mass number)
Z determines to which element a nucleus belongs and is equal to the number of electrons in
a neutral atom. The number of electrons determines the chemical properties of the elements.
As a general term for the type of nucleus determined by Z and N we use the word nuclide.
Isotopes are nuclides of the same element (same Z), with a different number of neutrons
(thus varying N and A). The term is derived from the Greek iso topos = same place (in the
periodic system).
There are two forms of notation for isotopes:
Full notation:
Short notation:
A
Z
N , e.g.
14
6
14
, e.g.
C8
The full notation is rarely used, because the short notation contains all information, if one
knows the number of protons (Z) of the elements.
The Chart of the Nuclides
The chart of the nuclide chart provides an overview of the existing isotopes and is an
important source of basic information for isotope sciences. The isotopes are mapped in plot
of neutron number N versus proton number Z. Isotopes of a single element stand in rows
Nuclear Methods in Hydrology
Page 3 of 14
(Fig. 1). Stable isotopes of the light elements align roughly along the 1:1 diagonal, for
heavier elements, there is always an excess of neutrons over protons.
Fig. 1: The lowest part of the chart of the nuclides. Adapted from Mook (2001).
Page 4 of 14
Isotopic abundance (ratio of the abundance of a given isotope to that of all isotopes of the
element) of the stable isotopes of the light elements. There is always one dominating light
isotope, and one or two rare, heavier isotope(s) with higher neutron number.
There are several different modes of radioactive decay. Particularly important for the
isotopes that are applied in hydrology are the and - decay. In an -decay an -particle,
which is nothing else than the very stable 4He nucleus, is emitted. The proton and neutron
numbers of the decaying nuclide are thus both reduced by 2, and the mass number is
reduced by 4. In a --decay a --particle, which is nothing else than an electron, is emitted,
and a neutron is converted to a proton. Consequently, the proton number is enlarged by 1
while the neutron number is reduced by 1, and the mass number remains the same. The
effects of radioactive decays on the position of an isotope in the chart of nuclides is shown in
Fig. 2.
Nuclear Methods in Hydrology
Page 5 of 14
Radioactive decay is a purely statistical process, with constant decay probability. The
number of decays in a given time interval is proportional to the total number of nuclides:
dN
N
dt
(4.1)
with [T-1] is the decay constant, i.e. the decay probability per unit time. The decay
constants of the different radioisotopes are characteristic material constants. They are in
principle well-known, although the precision of the values for long-lived isotopes is not
always very high.
Integration of eq. (4.1) with the initial condition N(0) = N0 yields
N (t ) N 0 e t
(4.2)
This is the well-known equation of an exponential decay (Fig. 3). The time dependence of
this decay is entirely characterised by the decay constant . Instead of , it may be more
convenient to use its inverse = 1/, which has the dimension of a time. corresponds to the
time when the number of nuclides has declined by a factor of 1/e and it equals the
expectation value for the lifetime of the nuclides. It is thus called the mean life of the isotopes.
Although mathematically less elegant, a much more customary measure of the characteristic
time of a radioactive decay process is the half-life T1/2, which is the time when the number of
nuclides has declined to half of its initial value. It can easily be seen that the half-life is linked
to t by a factor of ln2:
T1 / 2
ln 2
ln 2
N0
)
2
(4.3)
Eq. (4.3) is a simple but very useful relationship. It may even be worthwhile to remember that
ln2 is about 0.7, thus the half-life is about 70% of the mean life.
Radioactive isotopes can yield time information (age), if N(t) and N0 are known (the s are
known constants). The age is obtained simply by solving eq. (4.2) for the time:
Page 6 of 14
N (t )
N
ln 0
ln
N0
N (t )
1
(4.4)
In applications of this equation for dating, time zero is some unknown date in the past,
whereas time t is usually the present time or more precisely the time of analysis. N(t) is thus
directly measured, the problem is to know N0! We will address this problem and possible
solutions in some detail later on in the discussion of individual methods. The basic idea is
usually that N0 is assumed to be a known constant value. Another very important approach
is to measure the product of the decay (the so-called daughter isotope) in addition to the
decaying radioisotope itself (the mother isotope).
Radioisotopes are usually rather rare (because they constantly disappear), but the fact that
they decay offers an opportunity to detect them even at very low abundance. Thus, in
practice, often not N (number of atoms, abundance) is measured, but the activity, which is
the number of decays per unit time:
Definition of activity:
dN
N
dt
(4.5)
Page 7 of 14
O + p 10Be + 2n + 5p
16
O (p,2n5p) 10Be
or short
Another example is thermal neutron capture, where low-energy neutrons react less violently
with nuclides. The most commonly used environmental radioisotopes for groundwater
studies are tritium (3H) and radiocarbon (14C). Tritium is produced naturally in the upper
atmosphere through interaction of cosmic ray produced neutrons with nitrogen atoms:
14
14
C + p (proton)
10,000
Vienna
KAITOKE,
NEW ZEALAND
TRITIUM [TU]
1,000
VIENNA,
AUSTRIA
100
10
1955
1960
1965
1970
1975
1980
1985
1990
YEAR
Fig. 4:
Page 8 of 14
5. Isotope Fractionation
Isotope fractionation is defined as the phenomenon that the isotopic composition of an
element in a certain compound changes by the transition of the compound from one physical
state or chemical composition to another. One distinguishes three mass-dependent isotope
fractionation processes, namely, thermodynamic (in physical or chemical equilibrium
systems), kinetic (in one-way (bio)chemical reactions) and transport fractionation during
diffusive processes.
Within the hydrologic cycle the variability in the isotope composition results primarily from
mass-dependent isotope fractionation accompanying the phase transitions and transport
processes in the cycle. A special case of mass-independent fractionation is the phenomenon
that 17O16O in the atmospheric oxygen is affected by collisions of atmospheric gases with
cosmic particles and photochemical reactions in the higher atmosphere.
Isotope fractionation (isotopic enrichment of lighter isotopes):
Stable isotopes of water varies in time and space due to equilibrium fractionation or
kinetic fractionation.
Equilibrium isotope fractionation occurs during phase changes, transportation, oxidation
Fractionation refers to the partial separation of two isotopes of the same element,
producing reservoirs with different ratios of the isotopes. These differences are called
isotopic anomalies.
Two classes of basic mechanisms and two strange ones exist for fractionating
isotopes:
The two basic mechanisms are:
equilibrium isotope fractionation, which is due to differences in bond
energies of isotopes in compounds
kinetic isotope fractionation, which is due to differences in average
velocity or reaction rates of different isotopes
Both depend only on the mass of the isotope and are called mass
dependent fractionation; both will fractionate, say 18O/16O about twice
as much as 17O/16O
The strange mechanisms are mass-independent fractionation; one depends
on quantum mechanical symmetry terms that occur for molecules like O3 and
SO2 and reflects only the identity (not mass) of atoms, so fractionates 17O and
18
O equally. The other depends on isotope-specific wavelengths of photodissociation reactions and is thought perhaps to have been relevant in the
solar nebula.
Page 9 of 14
Since hydrologic processes involve only mass fractionation, the hydrogen isotopes
and oxygen isotopes co-vary quite closely, defining the meteoric water line D = 8 *
18O + 10
On a global scale, the transfers of water to the atmosphere from the oceans, and the
return via rainout and runoff, are close to a dynamic equilibrium (steady state).
Because of this, there is a linear relationship between 18O and 2H in precipitation, as
discovered by Craig (1961). The slope of the line relating 18O and 2H has a value
of 8, because 2H is 8 times more depleted in water vapor than is 18O (due to
differences in vapor pressures).
(6.1)
This Global Meteoric Water Line (GMWL) provides a useful benchmark against which
regional or local waters can be compared and their isotopic composition interpreted. The
slope of this curve represents Rayleigh fractionation due to repeated evaporation and
precipitation, the intercept (termed the deuterium excess) is largely a function of the mean
relative humidity of the atmosphere above the ocean water (Clark and Fritz,, 1997). Local
meteoric water lines can be established from isotopic analysis of local precipitation events.
The relationship has been confirmed in many later studies, including that of Rozanski et al.
(1993). The deuterium excess (d-excess) value is defined as:
d = D 8* 18O VSMOW (Dansgaard, 1964)
(6.2)
A d-excess of 10 is the average value for global precipitation and derived waters formed
from the vapour that was evaporated in the global average humidity of 85% over oceans
(Merlivat & Jouzel 1979; Gonfiantini 1986). Essentially, the d-excess is a measure that
characterizes the meteoric conditions prevailing in the initial moisture source (Merlivat &
Jouzel 1979). In addition, it is affected by secondary processes such as the evaporation of
raindrops beneath a cloud or evaporation of continental water bodies adding extra moisture
to the air mass (Rozanski 1984; Froehlich et al. 2002). The latter would increase the dexcess in precipitation, whereas sub-cloud evaporation, induced for example by the lowered
Nuclear Methods in Hydrology
Page 10 of 14
relative humidity of terrestrial areas, would decrease the d-excess. Seasonal variations in dexcess values are distinctive, especially for the precipitation of the northern hemisphere.
Through the combination of all primary and secondary processes, the d-excess in global
atmospheric precipitation generally varies from -2 to 15 (Froehlich et al. 2002). Waters
with low 18O and D values are fairly systematically encountered in cold regions and those
with high values in warm regions (Craig 1961).
AQUIFER CHARACTERISTICS
Porosity, dispersivity, transmissivity
H-18O, 3H(-3He),
Kr, CFCs, 14C-13C
85
GROUNDWATER QUANTITY
Area and rate of recharge
H-18O, 3H(-3He)
Kr, CFCs, 14C-13C
3
H(-3He), 85Kr, CFCs
14
C-13C, 36Cl, ...
2
H-18O, 3H, 222Rn
85
Flow regime
Interaction surface groundwater
Interaction between aquifers
Verification of conceptual models
All
GROUNDWATER QUALITY
Pollution sources
Transport of contaminants
Chemical and biological processes
isotopes of B, C, N, S, Cl
3
H(-3He), 85Kr, CFCs,
14
C-13C, U/Th series RN
isotopes of B, C, N, S
Page 11 of 14
Often isotope techniques are required as a complementary approach, and often they are
considered to be more cost-effective than conventional approaches. But there are a number
of problems where the use of isotope techniques is indispensable (Table 4).
Table 4: Fields in which the Application of Isotope Techniques is Indispensable
FP fingerprint (qualitative)
2)
FP1)
QE2)
QE
FP
FP
FP
FP
QE
QE
QE
FP
FP/QE
QE
FP/QE
FP/QE
FP/QE
QE quantitative evaluation
Page 12 of 14
Half-Life
Radiation &
[Energy of Interest
in Mev]
Chemical
Form
Bromine-82 (82Br)
35.4 h
[0.55-1.48]
NH4Br
Gold-198 (198Au)
2.7 d
[0.41]
AuCl3
Iodine-131 (131I)
8.1 d
[0.36]
NH4I
Chromium-51 (51Cr)
27.8 d
[0.32]
Cr-EDTA
Scandium-46 (46Sc)
84.0 d
[0.9, 1.1]
Sc2O3
Cobalt-60 (60Co)
5.2 a
[1.17, 1.33]
K3Co (CN)6
6h
[0.140]
Pertechnate
12.4 a
[0.018]
HTO
Isotopes
Technetium-99m (99mTc)
Tritium (3H)
Page 13 of 14
Official authorization for the acquisition and use of the radioactive material is required.
In the case of short-lived radioisotopes, they have to be acquired and used at a
specific time.
In some cases, heavy shielding lead containers have to be used for their transport,
such as for bromine-82 (82Br), scandium-46 (46Sc) and cobalt-60 (60Co).
In some cases, the use of radioactive tracers may result in some kind of risk for the
population and for the field team involved in their use.
Measurement of Radiotracers
Gamma emitting radioactive tracers are measured with the help of a radiation probe
incorporating a NaI (T1) scintillation detector. Field measurements are possible using a
portable rate meter/scaler.
Tritium being a beta emitter is measured using a liquid scintillation counter in the
laboratory.
References:
Clark, I. D. and Fritz, P., 1997. Environmental Isotopes in Hydrogeology, Lewis Publishers,
Boca Raton, FL, New York, USA.
Fritz, P. and Fontes, J. Ch. (Eds), 1980. Handbook of Environmental Isotope Geochemistry,
Volume 1, The Terrestrial Environment, A. Elsevier Scientific Publishing Company, New
York, USA.
International Atomic Energy Agency, 1981. Stable Isotope Hydrology: Deuterium and
Oxygen-18 in the water cycle, Technical Reports Series, No. 210, IAEA, Vienna, Austria.
International Atomic Energy Agency, 1983. Guidebook on Nuclear Techniques in Hydrology,
Technical Reports Series, No. 91, IAEA, Vienna, Austria.
Kendall, C. and Caldwell, E.A., 1998. Fundamentals of Isotope Geochemistry. Chapter 2 in:
Isotope Tracers in Catchment Hydrology (C. Kendall and J.J. McDonnell, editors), Elsevier,
New York, USA, pp. 51-86.
Mook, W. G. (Ed.), 2000. Environmental Isotope in the Hydrological Cycle, Principles and
Applications, Volume 1, Introduction: Theory, Methods and Review. International
Hydrological Programme, IHP-V, Technical Documents in Hydrology, No. 39, Vol. 1,
UNESCO and IAEA, Paris, France.
Nuclear Methods in Hydrology
Page 14 of 14