Sources and References

Primary textbook — White, W.M. (2013). Geochemistry. Wiley-Blackwell, Oxford. An authoritative graduate-level treatment covering cosmochemistry, isotope geochemistry, and aqueous systems in rigorous mathematical detail.

Supplementary texts — Rollinson, H.R. (1993). Using Geochemical Data: Evaluation, Presentation, Interpretation. Longman, London. Faure, G. and Mensing, T.M. (2005). Isotopes: Principles and Applications (3rd ed.). Wiley. Langmuir, C.H. and Broecker, W. (2012). How to Build a Habitable Planet. Princeton University Press.

Online resources — EarthChem Portal (earthchem.org) for geochemical databases; GEOROC (geochemical data repository, Max Planck Institute); USGS National Geochemical Survey (pubs.usgs.gov); NASA Earth Observatory for planetary and remote sensing context; NOAA for water chemistry and ocean data.

---

Chapter 1: Origin of the Elements and the Periodic Table

The Cosmic Context of Geochemistry

Geochemistry is fundamentally the application of chemical principles to the study of the Earth and other planetary bodies. It asks questions of breathtaking scope: why does the planet have the composition it does, and not some other composition? Why are some elements abundant and others vanishingly rare? Why do particular elements concentrate into specific minerals, rocks, fluids, and biological organisms? The answers to all of these questions begin not with the Earth itself but with the birth of the universe and the nuclear physics that determined which atoms were forged in stellar interiors and scattered across the cosmos before our solar system ever condensed.

The discipline traces its modern origins to the early twentieth century, when Victor Goldschmidt pioneered the systematic study of elemental distributions in natural materials. Goldschmidt recognized that the periodic table is not simply a chemical catalogue but a map of cosmic history — each element’s abundance in the solar system reflects the specific nuclear reactions that created it, and each element’s behaviour in terrestrial systems reflects the electronic architecture of its atoms. His classification scheme, discussed in detail below, remains foundational to how geochemists think about which elements go where in the Earth.

Understanding geochemistry therefore requires starting at the largest possible scale. The universe began approximately 13.8 billion years ago in the Big Bang, an event that produced primarily hydrogen (roughly 75% by mass), helium (roughly 25%), and trace amounts of lithium and beryllium through primordial nucleosynthesis. Every element heavier than these was subsequently manufactured inside stars, through a series of nuclear fusion and neutron-capture reactions collectively called stellar nucleosynthesis. The composition of the Earth, and indeed of all rocky planets, is ultimately a record of this nuclear history filtered through the chemistry of the solar nebula from which the solar system condensed roughly 4.57 billion years ago. The extraordinary fact that we can read this history in the chemistry of rocks collected at the surface is what makes geochemistry one of the most intellectually rewarding of the Earth sciences.

The study of element abundances in the solar system reveals patterns that carry information about which nuclear processes were active and how efficient they were. The Oddo-Harkins rule — that elements with even atomic numbers are systematically more abundant than those with odd atomic numbers — reflects the stability of alpha-particle configurations in nuclear physics. The abundance peaks at iron and nickel mark the maximum of nuclear binding energy per nucleon. The twin humps in the abundance pattern at barium and lead mark the sites of closed neutron shells (magic numbers) that inhibit further neutron capture during the s-process. Reading these patterns quantitatively is the business of nuclear astrophysics, but recognizing them qualitatively is essential background for every geochemist.

Stellar Nucleosynthesis: Building the Elements

The production of elements in stars proceeds through several distinct processes, each operating under different physical conditions and producing characteristic sets of nuclides. Understanding these processes is essential for interpreting the patterns of elemental abundance observed in meteorites, in the bulk Earth, and in the solar photosphere — all of which provide geochemists with crucial constraints on the planet’s starting composition.

Hydrogen burning (the proton-proton chain and CNO cycle) operates in main-sequence stars like the Sun and fuses four protons into a helium-4 nucleus, releasing energy through the mass defect. This is the dominant energy source for most of a star’s lifetime. When hydrogen in the stellar core is exhausted, the core contracts and heats, igniting helium burning (the triple-alpha process), in which three helium-4 nuclei fuse to produce carbon-12. The addition of another alpha particle yields oxygen-16. These reactions explain why carbon and oxygen are among the most abundant elements in the universe beyond hydrogen and helium. Subsequent stages of stellar evolution in massive stars ignite carbon burning, neon burning, oxygen burning, and silicon burning at progressively higher temperatures, synthesizing elements up to iron (atomic number 26) and nickel. Iron-56 represents the peak of nuclear binding energy; no energy can be liberated by fusing nuclei heavier than iron, which is why iron is the ultimate endpoint of charged-particle fusion in stellar interiors.

Elements heavier than iron are produced almost entirely through neutron-capture processes, because neutrons carry no charge and can therefore be absorbed by heavy nuclei without needing to overcome a Coulomb barrier. The slow neutron-capture process (s-process) operates in the interiors of asymptotic giant branch (AGB) stars over timescales of thousands to millions of years. Neutrons are captured one at a time, and between captures, beta decay converts neutrons to protons, building elements heavier than iron in a path that stays close to the valley of beta stability in the chart of nuclides. The s-process produces roughly half of the elements heavier than iron, including barium, strontium, and lead. The rapid neutron-capture process (r-process) requires an extraordinary neutron flux and occurs during core-collapse supernovae or neutron-star mergers. Neutrons are captured so quickly that the resulting nuclei are extremely neutron-rich; they subsequently decay toward stability, producing many of the heavy nuclides not made by the s-process, including the majority of the rare earth elements on the neutron-rich side of stability, as well as uranium and thorium. In 2017, the detection of gravitational waves from the neutron-star merger GW170817 followed by observation of a kilonova — a short-duration optical and infrared transient powered by rapid neutron-capture nucleosynthesis — provided direct observational confirmation that neutron-star mergers are a major site of r-process element production, a hypothesis that had been debated for decades.

A third process, the proton-capture or p-process, is responsible for a small set of proton-rich nuclides that cannot be produced by neutron capture. These p-nuclides are typically rare and their production sites (possibly Type Ia supernovae or the inner regions of massive star explosions) remain an active area of research. Their geochemical importance is limited by their scarcity, but the isotopic anomalies they can produce in meteorites are powerful tracers of the heterogeneity of the presolar material incorporated into the solar nebula.

The Periodic Table as a Geochemical Map

The modern periodic table arranges elements in order of increasing atomic number (the number of protons in the nucleus), with elements grouped in vertical columns (groups) sharing similar electronic configurations and therefore similar chemical behaviour. For geochemists, the periodic table is indispensable because it predicts which elements will substitute for one another in crystal lattices, which elements will partition into silicate melts versus metallic phases, and which elements will form soluble aqueous species under oxidizing versus reducing conditions.

The chemical properties of an element are governed primarily by the configuration of its outermost (valence) electrons. The principal quantum number \( n \) describes the electron shell, the angular momentum quantum number \( l \) describes the subshell (s, p, d, or f), the magnetic quantum number \( m_l \) describes the orbital orientation, and the spin quantum number \( m_s \) (±1/2) describes electron spin. The Pauli exclusion principle requires that no two electrons in an atom share the same four quantum numbers, which limits each orbital to two electrons. These quantum mechanical rules determine the filling order of orbitals and thus the structure of the periodic table: the s-block (Groups 1–2), the p-block (Groups 13–18), the d-block (transition metals, Groups 3–12), and the f-block (lanthanides and actinides).

Several properties that vary systematically across the periodic table are of direct geochemical importance. Atomic radius decreases across a period (left to right) as the nuclear charge increases and pulls electrons inward, and increases down a group as electrons occupy successively higher shells. Ionization energy — the energy required to remove an electron from a gaseous atom — increases across a period and decreases down a group. Electronegativity, which measures the tendency of an atom to attract electrons in a bond, increases across a period and decreases down a group; Linus Pauling’s electronegativity scale is the most commonly used in geochemistry. These trends control how elements bond with oxygen (the dominant anion in most geological materials), how stable particular oxidation states are, and whether an element tends to form ionic or covalent bonds in silicate structures.

The lanthanide contraction is a phenomenon of major geochemical importance that arises in the f-block elements (atomic numbers 57–71, the rare earth elements). As the 4f electron subshell is progressively filled across the lanthanide series, the nuclear charge increases but the 4f electrons are poor at shielding one another from the nucleus. The result is that the ionic radii of the REE decrease steadily from La³⁺ (1.032 Å in 6-coordination) to Lu³⁺ (0.861 Å in 6-coordination) — a contraction of about 16% across 15 elements. Because all REE carry the same 3+ charge, this progressive size change is the primary variable controlling their distribution among minerals, and the systematic change in compatibility across the REE series makes them extraordinarily powerful geochemical tracers, as discussed in Chapter 8.

Goldschmidt’s Geochemical Classification

Victor Goldschmidt developed his classification in the 1920s by studying the distribution of elements in meteorites, which he used as natural experiments in planetary differentiation. When he examined chondrite meteorites (primitive meteorites that retain much of the original solar nebula composition), he found that metal and sulphide phases had segregated from silicate phases, and that particular elements were consistently concentrated in one phase or another. This observation led him to recognize that the behaviour of elements during planetary differentiation — specifically, during the separation of a metallic core from a silicate mantle — could be predicted from basic chemical principles, particularly from the relative stabilities of the oxide and sulphide forms of each element.

Lithophile elements (from the Greek lithos, stone) have a strong affinity for oxygen and therefore concentrate in silicate phases. They are typically elements whose ions have a high charge-to-radius ratio or that form very stable oxides and silicates. Major lithophile elements include silicon, aluminium, calcium, magnesium, sodium, potassium, and titanium — the dominant constituents of the continental crust and mantle. Trace lithophile elements include the rare earth elements (REE), barium, strontium, rubidium, zirconium, hafnium, niobium, tantalum, uranium, and thorium. Because the crust and mantle are silicate systems, the lithophile elements are the primary focus of most igneous and metamorphic geochemistry. The lithophile elements are further subdivided by geochemists into the large-ion lithophile elements (LILE: K, Rb, Cs, Ba, Sr) and the high-field-strength elements (HFSE: Zr, Hf, Nb, Ta, Ti, U, Th), which have very different behaviour during fluid-rock interaction and magma genesis, as discussed in Chapter 8.

Siderophile elements (from the Greek sideros, iron) prefer to associate with metallic iron. They include iron itself, nickel, cobalt, and the platinum-group elements (PGE: osmium, iridium, ruthenium, rhodium, palladium, platinum), as well as gold, rhenium, and molybdenum. These elements are strongly depleted in the Earth’s mantle and crust relative to primitive chondrite meteorites because they partitioned into the metallic core during Earth’s early differentiation. The extraordinary depletion of PGE in the mantle (concentrations typically at the parts-per-billion level) is direct evidence that a metal-silicate separation event occurred when the proto-Earth was largely or entirely molten — the magma ocean stage. The small residual abundances of PGE in the mantle are attributed to the “late veneer” — accretion of chondritic material after core formation was complete — though this interpretation has been disputed and alternative explanations involving incomplete metal-silicate equilibration or inefficient PGE extraction remain under discussion.

Chalcophile elements (from the Greek chalkos, copper) preferentially associate with sulphur and concentrate in sulphide phases. The most important chalcophile elements are copper, zinc, lead, silver, arsenic, antimony, bismuth, selenium, tellurium, and cadmium. These are the elements that form economically significant sulphide ore deposits — the ores of copper (chalcopyrite, bornite), lead (galena), and zinc (sphalerite) are all sulphides. The chalcophile character of an element reflects the relative stability of its sulphide compared to its oxide; elements that form more stable sulphides than oxides under reducing conditions are chalcophile. The boundary between siderophile and chalcophile behaviour is somewhat blurred in some elements — platinum-group elements can also concentrate in sulphide phases under some conditions, which is the basis for magmatic PGE ore deposits (discussed in EARTH 471).

Atmophile elements (from the Greek atmos, vapour) are those concentrated in the atmosphere and hydrosphere. These are primarily the noble gases (helium, neon, argon, krypton, xenon), as well as nitrogen, hydrogen, and (under some definitions) oxygen and carbon when they exist in gaseous form. In practice, geochemists rarely treat oxygen or carbon as atmophile since they form abundant oxide and carbonate minerals, but the noble gases are essentially entirely atmophile because they form no chemical bonds and cannot be incorporated into crystal structures. The noble gas geochemistry of the mantle — particularly the ratios of different helium, neon, and argon isotopes — is one of the most powerful tools for tracing the deep degassing history of the Earth and the mixing of primitive (undegassed) and depleted (degassed) mantle reservoirs.

Crystal Chemistry and Goldschmidt’s Rules of Ionic Substitution

The incorporation of trace elements into crystal lattices is governed by the principles of crystal chemistry, which Goldschmidt codified into a set of rules that remain central to igneous and metamorphic petrology. These rules explain which trace elements substitute for which major elements in minerals, thereby controlling the distribution of trace elements between minerals and melts during igneous processes.

The first rule states that if two ions have similar ionic radii (within roughly 15% of each other) and the same charge, they will substitute freely for each other in a crystal. The classic example is the substitution of Sr²⁺ (ionic radius 1.18 Å in 6-coordination) for Ca²⁺ (ionic radius 1.00 Å in 6-coordination) in calcite and feldspars. Although the radii differ by about 15%, both ions carry a 2+ charge, and Sr readily enters calcium sites in these minerals. Similarly, Rb⁺ (1.52 Å) substitutes for K⁺ (1.38 Å) in potassium-bearing minerals such as K-feldspar, phlogopite, and muscovite. This substitution is of enormous geochronological importance because Rb decays to Sr (see Chapter 4), and the K-rich minerals that preferentially incorporate Rb therefore accumulate radiogenic Sr over geological time.

The second rule concerns the competition between ions of similar radius but different charge: if two ions have similar ionic radii but different charges, the ion of higher charge is preferentially incorporated into the crystal when the site has a high bond strength requirement (i.e., in highly charged crystal sites). Conversely, the ion of lower charge tends to be excluded from such sites. This explains why Ba²⁺ (1.35 Å) prefers to substitute into K⁺ (1.38 Å) sites in feldspars only under specific conditions, while the higher-charge Ce³⁺ may be preferentially incorporated relative to La³⁺ in some minerals despite the charge mismatch with divalent sites, because the crystal site geometry may favour the smaller, more highly charged ion.

The third rule recognises that if two ions of different charge compete for the same site, the ion that forms the stronger bond with the anion (typically oxygen in silicate systems) will be preferentially incorporated. The strength of the ionic bond is proportional to the charge and inversely proportional to the square of the bond distance. This is why highly charged, small ions (such as Zr⁴⁺, Hf⁴⁺, or Ti⁴⁺) have very different crystal chemical behaviours from similarly sized but lower-charge ions, and why the high-field-strength elements (HFSE) are generally very insoluble in aqueous fluids and resistant to mobility during metamorphic fluid-rock interaction.

The ionic radius plays a central role in all of these rules, and it is critical to remember that ionic radii are coordination-dependent — an ion surrounded by four oxygen neighbours (tetrahedral coordination) has a different effective radius than the same ion in octahedral (six oxygen) or eight-fold coordination. The Shannon (1976) compilation of ionic radii is the standard reference used by geochemists, and it must always be specified whether a radius refers to tetrahedral, octahedral, or some other coordination. The commonly cited radii for geochemically important ions are those in octahedral coordination (VI), unless the specific coordination environment of the crystal site is known.

---

Chapter 2: Earth’s Formation and Planetary Geochemistry

From Solar Nebula to Planet

The formation of the Earth is inseparable from the formation of the solar system as a whole, which began when a region of the interstellar medium — enriched with heavy elements from previous generations of stars — collapsed under its own gravity approximately 4.568 billion years ago. The collapsing cloud formed a central proto-Sun surrounded by a rotating disk of gas and dust, the solar nebula. Within this disk, temperature and pressure decreased with increasing distance from the proto-Sun, creating a temperature gradient that profoundly influenced which materials condensed from the gas phase at different heliocentric distances.

The condensation sequence — the order in which different minerals and metals condense from a cooling gas of solar composition — is one of the most fundamental concepts in cosmochemistry. At the highest temperatures (above about 1,800 K at nebular pressures), refractory oxides and silicates condense first: corundum (Al₂O₃), calcium aluminates, and melilite. These refractory materials are preserved as calcium-aluminium-rich inclusions (CAIs) in chondritic meteorites and represent the oldest solid objects in the solar system, dated at 4.5682 ± 0.0001 Ga by U-Pb isotope geochronology. At lower temperatures, magnesium silicates (forsterite, enstatite) and then metallic iron-nickel alloy condense. Below about 700 K, iron sulphide (troilite, FeS) forms from the reaction of iron metal with nebular hydrogen sulphide, and at even lower temperatures, phyllosilicates and carbonates form on the surfaces of existing grains. The inner solar system, being hot, was dominated by refractory and metal-rich condensates, which is why the terrestrial planets (Mercury, Venus, Earth, Mars) are dense, silicate-metal bodies rather than ice-rich like the outer planets and their satellites.

The small dust grains in the solar nebula aggregated into successively larger bodies through a process of accretion. Grains first formed centimetre-to-metre-scale aggregates (the so-called metre-scale barrier is a major unsolved problem in planet formation theory because objects of this size are most susceptible to orbital decay by gas drag). Somehow, and the mechanism remains debated, planetesimals of kilometre scale formed, and from that scale gravitational attraction dominates accretion dynamics. Planetesimals sweep up material on crossing orbits, growing into planetary embryos (Moon- to Mars-mass bodies), which then collide with one another in a final stage of giant impacts over tens to hundreds of millions of years. The Earth’s Moon is thought to have formed from the debris ejected when a Mars-sized body (sometimes called Theia) collided with the proto-Earth approximately 4.5 billion years ago — the Giant Impact hypothesis, which is supported by the Moon’s nearly identical oxygen isotope composition to Earth and its depletion in volatile elements.

The volatile element budget of the Earth — its inventory of H, C, N, S, and the noble gases — was determined largely by the processes of accretion and differentiation. Volatile elements are those that condense from a cooling gas at relatively low temperatures; they are found in much lower abundances on the Earth and other inner planets than in CI chondrites. The degree of volatile depletion increases from Mars to Earth to the Moon and to Mercury, broadly consistent with the idea that inner planets formed from more thermally processed material and/or lost their volatiles during giant impacts. The late addition of volatile-rich material (the “late veneer” mentioned earlier) may have contributed some fraction of Earth’s water, carbon, and nitrogen, but the dominant volatile budget was likely delivered earlier in the accretion process. This is an area of very active research, with isotopic constraints on the sources of Earth’s hydrogen (D/H ratios of carbonaceous chondrites vs. comets) and nitrogen (\( \delta^{15}N \) systematics) providing key clues.

Meteorites as Windows into Planetary Formation

Meteorites are fragments of solar system bodies — asteroids, and occasionally Mars and the Moon — that survive passage through Earth’s atmosphere and land on the surface. They represent the best samples of undifferentiated and partially differentiated planetary materials available to geochemists, and they have been studied intensively since the nineteenth century as records of conditions in the early solar system.

The most primitive meteorites are the chondrites, so named because they contain millimetre-scale spherical structures called chondrules — droplets of silicate melt that solidified rapidly, probably during transient heating events in the solar nebula such as lightning discharges or bow shocks ahead of large bodies. The most chemically primitive chondrites are the CI carbonaceous chondrites (named after the Ivuna meteorite), which have elemental compositions that closely match the photospheric composition of the Sun for all non-volatile elements. This agreement is remarkable and confirms that both the Sun and the CI chondrites sampled the same well-mixed starting material. Geochemists therefore use CI chondrite composition as the reference standard for “primitive” or “chondritic” elemental abundances, and deviations from this reference in rocks from differentiated bodies (including the Earth’s mantle and crust) reveal the extent and nature of planetary differentiation.

Among the differentiated meteorites, the iron meteorites are particularly informative. These are fragments of the metallic cores of asteroidal bodies that were melted, differentiated, and subsequently disrupted by impact. Their elemental compositions — principally iron and nickel with trace amounts of siderophile elements — provide direct samples of what a small planetary core looks like. The patterns of Ni, Co, Ir, Au, Ge, and other siderophile elements within iron meteorites are diagnostic of crystallization in a slowly cooling metallic magma, and the parent bodies of the different iron meteorite groups can be distinguished by their distinct chemical compositions and crystallization ages.

The pallasites — meteorites consisting of olivine crystals embedded in a metallic matrix — are interpreted as samples from the core-mantle boundary of small differentiated asteroids. They beautifully illustrate the physical separation between mantle silicates (olivine) and core metal (iron-nickel) that occurs during planetary differentiation, and they have been used to constrain the timing of core formation in small bodies (generally within the first 5–20 million years of solar system history, based on short-lived isotope systems like ¹⁸²Hf-¹⁸²W). The Hf-W chronometer is especially powerful for dating metal-silicate separation events: ¹⁸²Hf (half-life 8.9 Ma) decays to ¹⁸²W, and because Hf is lithophile while W is siderophile, core formation transfers W into the metal phase, leaving the silicate mantle with an elevated Hf/W ratio that subsequently drives the mantle ¹⁸²W/¹⁸⁴W ratio above the chondritic value. The excess ¹⁸²W in the Earth’s mantle indicates that core formation was largely complete within the first 30–50 Ma of solar system history.

Core Formation and the Bulk Silicate Earth

The differentiation of the Earth — the separation of the metallic core from the silicate mantle — was a process of profound geochemical significance that set the stage for all subsequent geological activity. The driving force was gravitational: liquid metal, being denser than silicate, sank toward the planet’s centre while buoyant silicate rose to form the mantle and crust. For this process to occur on a planetary scale, the Earth must have been partially or completely molten at some point — the magma ocean stage — and this is supported by models of accretional heating (kinetic energy of giant impacts), radiogenic heating from short-lived isotopes like ²⁶Al (half-life 0.7 Ma) and ⁶⁰Fe (half-life 2.6 Ma), and the thermal energy released by core formation itself (which is substantial, potentially several hundred degrees of mantle heating).

The geochemical consequences of core formation are dramatic. All siderophile elements were removed from the silicate portion of the Earth during this event, partitioning into the metallic core at extremely high efficiencies determined by the relevant metal-silicate partition coefficients. The degree of depletion of siderophile elements in the primitive upper mantle (PUM, also called the primitive mantle or bulk silicate Earth, BSE) relative to CI chondrites constrains the conditions of metal-silicate equilibration. Early studies assumed that metal-silicate separation occurred at low pressures and relatively low temperatures, but this would predict far more extreme depletion of moderately siderophile elements (Ni, Co) than is actually observed — the “excess siderophile problem” of the 1980s. The resolution came from recognizing that the partition coefficients of Ni and Co between metal and silicate melt become less extreme at the very high pressures and temperatures expected in a deep magma ocean (possibly 30–60 GPa and 3,000–4,500 K), so that their residual concentrations in the mantle are consistent with metal-silicate equilibration at depth.

---

Chapter 3: Thermodynamics of Geochemical Systems

Why Thermodynamics Matters in Geochemistry

Thermodynamics is the branch of physics and chemistry concerned with the relationship between heat, work, and energy, and with the direction and extent of chemical reactions. In geochemistry, thermodynamic principles allow us to predict what minerals are stable at a given temperature and pressure, what the composition of a fluid in equilibrium with a rock should be, in which direction a reaction will proceed, and what the equilibrium partitioning of elements between phases should be. Without thermodynamics, geochemistry would be merely descriptive; with it, geochemists can make quantitative predictions about processes they cannot directly observe, such as the conditions deep in the crust or in subducting slabs.

The First Law of Thermodynamics states that energy is conserved: the change in internal energy \( \Delta U \) of a system equals the heat \( q \) absorbed by the system minus the work \( w \) done by the system:

\[ \Delta U = q - w \]

For a system at constant pressure (the most common situation in geological processes), the relevant energy function is the enthalpy \( H = U + PV \), so that \( \Delta H = q_P \), the heat exchanged at constant pressure. Exothermic reactions (burning, hydration, crystallization of many minerals) release heat (\( \Delta H < 0 \)), while endothermic reactions (melting, some dehydration reactions) absorb heat (\( \Delta H > 0 \)).

The Second Law of Thermodynamics introduces the concept of entropy \( S \), a measure of the disorder or dispersal of energy in a system. The Second Law states that the entropy of the universe always increases in any spontaneous process. For a system at constant temperature, the entropy change of the surroundings equals the heat transferred from the system divided by temperature: \( \Delta S_{surroundings} = -\Delta H_{system} / T \). The combined criterion for spontaneity is therefore:

\[ \Delta S_{universe} = \Delta S_{system} - \frac{\Delta H_{system}}{T} \geq 0 \]

The thermodynamic function that combines both enthalpy and entropy and provides the most useful criterion for chemical equilibrium at constant temperature and pressure is the Gibbs free energy:

\[ G = H - TS \]

For a reaction to be spontaneous at constant \( T \) and \( P \), the Gibbs free energy must decrease:

\[ \Delta G = \Delta H - T\Delta S < 0 \]

At equilibrium, \( \Delta G = 0 \). These simple equations contain enormous predictive power for geochemical systems. A reaction can be spontaneous even if it absorbs heat (\( \Delta H > 0 \)) provided that the entropy increase is sufficiently large (\( T\Delta S > \Delta H \)); conversely, a reaction that releases heat may not be spontaneous at high temperatures if it involves a large decrease in entropy.

Chemical Potential and the Activity Concept

For systems containing multiple components (different chemical species) distributed among multiple phases (solid minerals, melt, aqueous fluid), the full thermodynamic treatment requires the concept of chemical potential. The chemical potential of component \( i \) in a phase is:

\[ \mu_i = \left(\frac{\partial G}{\partial n_i}\right)_{T,P,n_{j\neq i}} \]

This is the partial molar Gibbs free energy — the change in Gibbs free energy per mole of component \( i \) added to the phase, holding all other quantities constant. At equilibrium, the chemical potential of each component must be equal in all phases in which it occurs. This equality-of-chemical-potential condition is the thermodynamic basis for all phase equilibria and partitioning relationships in geochemistry. For a mineral in equilibrium with a melt, for example, the chemical potential of each component in the mineral must equal its chemical potential in the melt.

The chemical potential of a component in a non-ideal mixture is related to its activity \( a_i \) through:

\[ \mu_i = \mu_i^\circ + RT\ln a_i \]

where \( \mu_i^\circ \) is the standard chemical potential (in the pure end-member state at the same T and P) and \( R \) is the gas constant (8.314 J mol^-1 K^-1). The activity is related to composition through an activity coefficient \( \gamma_i \):

\[ a_i = \gamma_i x_i \]

where \( x_i \) is the mole fraction of component \( i \). For ideal solutions, \( \gamma_i = 1 \) and activity equals mole fraction. For real solutions, activity coefficients may be greater or less than 1, depending on whether the interactions between unlike molecules are more or less favourable than those between like molecules. In dilute aqueous solutions of electrolytes, the Debye-Hückel theory provides a means of calculating activity coefficients from ionic strength, which is particularly important for predicting the solubility of minerals in natural waters.

The Equilibrium Constant and Le Chatelier’s Principle

For a chemical reaction of the form \( aA + bB \rightleftharpoons cC + dD \), the equilibrium constant \( K \) is defined as:

\[ K = \frac{a_C^c \cdot a_D^d}{a_A^a \cdot a_B^b} \]

The equilibrium constant is related to the standard Gibbs free energy change of the reaction by:

\[ \Delta G^\circ = -RT\ln K \]

This equation is one of the most important in all of thermodynamics. It shows that if \( \Delta G^\circ \) is very negative (the products are much more stable than the reactants), then \( K \) is very large and the reaction goes essentially to completion. If \( \Delta G^\circ \) is positive, \( K \) is small and the reactants are favoured. The temperature dependence of \( K \) is given by the van’t Hoff equation:

\[ \frac{d\ln K}{dT} = \frac{\Delta H^\circ}{RT^2} \]

For an exothermic reaction (\( \Delta H^\circ < 0 \)), \( K \) decreases with increasing temperature — heating shifts equilibrium toward reactants. For an endothermic reaction (\( \Delta H^\circ > 0 \)), \( K \) increases with temperature — heating favours products. This is a quantitative expression of Le Chatelier’s principle: a system at equilibrium shifts to counteract any perturbation. The integrated form of the van’t Hoff equation, assuming \( \Delta H^\circ \) is constant over the temperature range of interest, is:

\[ \ln\frac{K_2}{K_1} = \frac{\Delta H^\circ}{R}\left(\frac{1}{T_1} - \frac{1}{T_2}\right) \]

The pressure dependence of \( K \) is given by:

\[ \frac{d\ln K}{dP} = -\frac{\Delta V^\circ}{RT} \]

where \( \Delta V^\circ \) is the volume change of the reaction. Reactions that produce a net decrease in volume (\( \Delta V^\circ < 0 \)) are favoured by increasing pressure. This is profoundly important in the deep Earth: the transformation of minerals from low-density to high-density polymorphs (e.g., olivine to wadsleyite at the 410-km seismic discontinuity, or calcite to aragonite in subducting slabs) is driven by the large negative volume change of the high-pressure phases.

---

Chapter 4: Geochronology — Dating Rocks and Events

Principles of Radioactive Decay

Geochronology — the science of measuring the ages of rocks, minerals, and geological events — is built on the physics of radioactive decay. Radioactive isotopes (radionuclides) are unstable nuclides that spontaneously transform into more stable daughter nuclides by emitting particles or energy. The rate of this decay is governed by a fundamental law: the number of atoms decaying per unit time is proportional to the number of radioactive atoms present. This leads to an exponential decay equation:

\[ N(t) = N_0 e^{-\lambda t} \]

where \( N(t) \) is the number of parent atoms remaining at time \( t \), \( N_0 \) is the initial number of parent atoms, and \( \lambda \) is the decay constant (units of time^-1), which is specific to each radionuclide. The half-life \( t_{1/2} \) — the time required for half of the original parent atoms to decay — is related to the decay constant by:

\[ t_{1/2} = \frac{\ln 2}{\lambda} = \frac{0.6931}{\lambda} \]

The daughter isotope \( D \) accumulates according to:

\[ D(t) = D_0 + N(e^{\lambda t} - 1) \]

where \( D_0 \) is the initial number of daughter atoms present when the rock or mineral formed. This is the fundamental geochronological equation. It contains three unknowns (\( D_0 \), \( N \), and \( t \)) and two observables (present-day \( D \) and \( N \)), so one additional constraint — either an assumption about the initial ratio or multiple mineral analyses — is needed to determine age. The power of geochronology rests on the invariance of decay constants: unlike virtually every other rate constant in chemistry, radioactive decay rates are completely unaffected by temperature, pressure, chemical form, crystal structure, or any other physical or chemical condition within the ranges encountered on Earth or in the solar system. This invariance has been tested experimentally and theoretically and is one of the most robustly established facts in physics.

Decay Modes and Geochemically Important Systems

Several modes of radioactive decay are geochemically important. Alpha decay involves the emission of an alpha particle (helium-4 nucleus), which reduces the parent atomic number by 2 and mass number by 4. Beta-minus decay converts a neutron to a proton, emitting an electron and an antineutrino, increasing the atomic number by 1. Beta-plus decay (positron emission) converts a proton to a neutron, decreasing the atomic number by 1. Electron capture (EC) accomplishes the same transformation as positron emission by capturing an inner-shell electron. Each of these processes produces a characteristic daughter nuclide and operates at a rate determined solely by the nuclear properties of the parent.

The major radioactive decay systems used in geochronology are summarized below:

Uranium-Lead (U-Pb): Two decay chains are relevant: ²³⁸U decays to ²⁰⁶Pb through a series of eight alpha and six beta decays (\( \lambda_{238} = 1.5512 \times 10^{-10} \) yr^-1, \( t_{1/2} = 4.468 \) Ga), and ²³⁵U decays to ²⁰⁷Pb through a series of seven alpha and four beta decays (\( \lambda_{235} = 9.8485 \times 10^{-10} \) yr^-1, \( t_{1/2} = 0.7038 \) Ga). The present-day ²³⁸U/²³⁵U ratio is 137.88. The simultaneous use of both systems (the concordia method) provides a powerful internal consistency check and allows detection of open-system behaviour.

Rubidium-Strontium (Rb-Sr): ⁸⁷Rb decays by beta emission to ⁸⁷Sr with \( \lambda = 1.402 \times 10^{-11} \) yr^-1 (\( t_{1/2} = 49.44 \) Ga). The long half-life makes this system ideal for dating old Precambrian rocks and for isotope tracing in mantle geochemistry.

Samarium-Neodymium (Sm-Nd): ¹⁴⁷Sm decays by alpha emission to ¹⁴³Nd with \( \lambda = 6.54 \times 10^{-12} \) yr^-1 (\( t_{1/2} = 106 \) Ga). The extremely long half-life means that the ¹⁴³Nd/¹⁴⁴Nd ratio changes very slowly over geological time, making Sm-Nd a powerful tracer of ancient mantle and crustal reservoirs. The epsilon-Nd (\( \varepsilon_{Nd} \)) notation expresses deviations of the measured ¹⁴³Nd/¹⁴⁴Nd ratio from the chondritic value in parts per 10,000: positive \( \varepsilon_{Nd} \) indicates derivation from a depleted mantle source (that has had melt extracted from it), while negative \( \varepsilon_{Nd} \) indicates derivation from enriched continental crust.

Potassium-Argon (K-Ar): ⁴⁰K decays by a branched scheme to both ⁴⁰Ca (89.5%) and ⁴⁰Ar (10.5%), with a total decay constant of \( \lambda = 5.543 \times 10^{-10} \) yr^-1. The ⁴⁰Ar/³⁹Ar variant, in which a nuclear reactor converts ³⁹K to ³⁹Ar and both isotopes are measured simultaneously, eliminates the need to measure potassium separately and allows step-heating experiments that reveal the thermal history of the sample. The closure temperature — the temperature below which the system behaves as closed to argon diffusion — varies by mineral: approximately 350°C for hornblende, 300°C for muscovite, 250°C for biotite, and 150°C for K-feldspar (multi-diffusion domain), allowing multi-mineral thermochronology to reconstruct the cooling history of a rock through different temperature windows.

The Isochron Method

The isochron method is a graphical and mathematical technique for determining the age of a rock or suite of rocks from the Rb-Sr, Sm-Nd, or other decay systems. The fundamental insight is that if a suite of rocks or minerals formed at the same time from a common reservoir (so they all started with the same isotope ratio of the daughter element), they will plot on a straight line (an isochron) in a diagram of \( {}^{87}Sr/{}^{86}Sr \) versus \( {}^{87}Rb/{}^{86}Sr \). The slope of this line is \( e^{\lambda t} - 1 \) and the y-intercept gives the initial \( {}^{87}Sr/{}^{86}Sr \) ratio at the time of formation.

The isochron equation is:

\[ \left(\frac{{}^{87}Sr}{{}^{86}Sr}\right)_{measured} = \left(\frac{{}^{87}Sr}{{}^{86}Sr}\right)_0 + \left(\frac{{}^{87}Rb}{{}^{86}Sr}\right) \left(e^{\lambda t} - 1\right) \]

This has the form of a straight line \( y = b + mx \), where \( y = {}^{87}Sr/{}^{86}Sr \), \( x = {}^{87}Rb/{}^{86}Sr \), the slope \( m = e^{\lambda t} - 1 \), and the intercept \( b \) is the initial ratio. A suite of co-genetic rocks with different Rb/Sr ratios (e.g., a mafic rock with low Rb/Sr and a granitic rock with high Rb/Sr, both derived from the same source by differentiation at the same time) will plot on such a line, and its slope yields the age \( t \) while its intercept yields the initial ratio.

The U-Pb Concordia Diagram

The concordia diagram is one of the most powerful tools in geochronology. Because uranium undergoes two simultaneous decay chains to lead (producing ²⁰⁶Pb from ²³⁸U and ²⁰⁷Pb from ²³⁵U), and because the two systems evolve at very different rates, the two daughter/parent ratios provide independent age information. A mineral that has remained as a closed system since crystallization will have equal ages from the two systems — it is concordant — and plots on the concordia curve in the \( {}^{207}Pb/{}^{235}U \) versus \( {}^{206}Pb/{}^{238}U \) diagram (Wetherill concordia) or in the Tera-Wasserburg form (\( {}^{207}Pb/{}^{206}Pb \) versus \( {}^{238}U/{}^{206}Pb \)). The age can be read directly from the position on the concordia curve.

If a zircon experienced lead loss (Pb is more mobile than U due to radiation damage of the crystal lattice, or during high-temperature metamorphic events), it will plot below the concordia curve — the two systems give discordant ages. A family of zircons that all experienced lead loss at the same time will define a chord (a discordia line) connecting the original crystallization age on the concordia curve to the age of the lead-loss event. The upper intercept gives the crystallization age, and the lower intercept gives the age of the disturbance. This capacity to deconvolve two-stage histories is one of the great advantages of the U-Pb system and is why zircon U-Pb geochronology is the workhorse of Precambrian geology worldwide.

---

Chapter 5: Earth’s Crust and Natural Resources

The Continental Crust: Composition and Significance

The continental crust is the geochemically diverse outer shell of the Earth that forms the foundation of all landmasses and continental shelves. It ranges from about 25–35 km thick in stable cratons to 60–70 km beneath the Himalayas and other young mountain belts. Its bulk composition is andesitic to dacitic — intermediate between the mafic composition of the oceanic crust (basaltic) and the felsic composition of granite. The average composition of the upper continental crust is well constrained through decades of geochemical studies combining direct sampling (major element analysis of rocks, trace element analysis by ICP-MS, isotope ratio measurements by TIMS or MC-ICP-MS), sediment geochemistry, and mass balance arguments.

The continental crust is strongly enriched in lithophile elements relative to the mantle. Silicon (SiO₂ ≈ 66 wt%), aluminium (Al₂O₃ ≈ 15 wt%), and alkali elements (Na₂O ≈ 3.3 wt%, K₂O ≈ 3.3 wt%) are all far more concentrated in the crust than in the primitive mantle. The trace elements K, Rb, Ba, Th, U, La, and Ce — the large-ion lithophile elements (LILE) and light rare earth elements (LREE) — are extraordinarily enriched, reflecting billions of years of partial melting of the mantle, which causes these highly incompatible elements to preferentially enter the melt and accumulate in the crust. The complementary depletion of these elements in the depleted mantle — the source of mid-ocean ridge basalts (MORB) — provides geochemical evidence that the two reservoirs evolved together as complementary differentiation products over Earth history.

The lower continental crust is more mafic in composition than the upper crust, and in many places it has granulite facies metamorphic assemblages (pyroxene-bearing, water-poor rocks that have lost their volatiles and incompatible elements during high-temperature metamorphism). The boundary between the crust and the mantle is the Mohorovicic discontinuity (Moho), marked by a sharp increase in seismic P-wave velocity from about 6–7 km/s in the lower crust to 8 km/s in the uppermost mantle, representing a compositional boundary between felsic-to-intermediate crustal rocks and peridotitic mantle.

---

Chapter 6: Aqueous Geochemistry — Water Quality and the Hydrological Cycle

The Hydrological Cycle and Natural Water Chemistry

The second half of EARTH 221 shifts from the solid Earth to the aquatic realm, beginning with the hydrological cycle — the continuous movement of water through the Earth system via evaporation, precipitation, surface runoff, groundwater flow, and ultimately return to the ocean. The global hydrological cycle transfers approximately 500,000 km³ of water per year through evaporation from the oceans alone, with roughly 90% of this returning directly to the oceans as precipitation and 10% being transported by winds to the continents. The net continental water balance drives river flow and groundwater recharge, which ultimately returns to the ocean through streams, rivers, and submarine groundwater discharge.

Natural waters acquire their dissolved loads through chemical weathering of rocks — the dissolution and alteration of minerals by water, carbonic acid (produced when CO₂ dissolves in rainwater), organic acids, and other agents. The primary reactions involve the hydrolysis of silicate minerals, the dissolution of carbonates, the oxidation of sulphide minerals, and cation exchange on clay mineral surfaces. These reactions collectively transfer elements from the solid to the aqueous phase and are the primary mechanism by which river water acquires its dissolved burden of Ca²⁺, Mg²⁺, Na⁺, K⁺, HCO₃^-, SO₄^2-, and SiO₂. The rate of weathering is controlled by temperature, rainfall, the reactivity of the minerals exposed, the presence of organic acids from soil decomposition, and the degree of biological activity. Tropical regions with warm temperatures, abundant rainfall, and dense vegetation weather rocks far faster than polar or arid regions, producing dramatically different water chemistries and different rates of solute export to the oceans.

The major ions in freshwater — and their concentrations relative to each other — can be displayed on a Piper diagram, which is a trilinear diagram that plots the relative proportions of major cations (Ca²⁺, Mg²⁺, Na⁺ + K⁺) and major anions (HCO₃^- + CO₃^2-, SO₄^2-, Cl^-) in ternary diagrams, with a central diamond diagram showing the overall water type. Piper diagrams are extensively used in hydrogeology to classify water types, identify mixing relationships between different water sources, and track chemical evolution along flow paths. A water dominated by Ca²⁺ and HCO₃^- (calcium bicarbonate type, characteristic of limestone aquifers) falls in a very different position from Na⁺-Cl^- water (saline brine type), and the evolution from one type to another along a groundwater flow path records the progressive dissolution of carbonate minerals, sulphates, or halite, or cation exchange reactions on clay minerals. Understanding which water type is present is essential for assessing suitability for drinking, irrigation, or industrial use, and for identifying potential contamination sources.

The Carbonate System and pH Controls in Natural Waters

The most important buffer system controlling the pH of natural waters is the carbonate system, which involves the equilibrium between dissolved CO₂, carbonic acid, bicarbonate ion, and carbonate ion. The relevant reactions and their equilibrium constants at 25°C are:

\[ \text{CO}_{2(g)} \rightleftharpoons \text{CO}_{2(aq)} \quad K_H = 3.4 \times 10^{-2} \text{ mol/L/atm} \]\[ \text{CO}_{2(aq)} + \text{H}_2\text{O} \rightleftharpoons \text{H}_2\text{CO}_3 \quad \text{(essentially instantaneous, combined with above)} \]\[ \text{H}_2\text{CO}_3 \rightleftharpoons \text{H}^+ + \text{HCO}_3^- \quad K_1 = 4.47 \times 10^{-7} \text{ (pK}_1 = 6.35) \]\[ \text{HCO}_3^- \rightleftharpoons \text{H}^+ + \text{CO}_3^{2-} \quad K_2 = 4.68 \times 10^{-11} \text{ (pK}_2 = 10.33) \]

In equilibrium with the pre-industrial atmosphere (\( p\text{CO}_2 = 10^{-3.5} \) atm), the pH of rain water is approximately 5.6 — slightly acidic even in the absence of pollution because of dissolved CO₂. When such water contacts carbonate minerals (calcite, dolomite), dissolution occurs:

\[ \text{CaCO}_3 + \text{CO}_2 + \text{H}_2\text{O} \rightleftharpoons \text{Ca}^{2+} + 2\text{HCO}_3^- \]

This reaction releases Ca²⁺ and HCO₃^- into solution, raising pH and increasing alkalinity. The reverse reaction (calcite precipitation) occurs when CO₂ degasses from water or when temperature increases, which is why spring waters emerging from limestone aquifers commonly deposit calcium carbonate travertine at the surface. Ocean acidification — the decrease in seawater pH caused by uptake of anthropogenic CO₂ — is a direct consequence of these equilibria: rising atmospheric CO₂ drives more CO₂ into solution, increasing H⁺ concentration and reducing the carbonate ion concentration, making it harder for marine organisms (corals, molluscs, echinoderms) to build calcium carbonate shells and skeletons. Since the industrial revolution, ocean surface pH has decreased from approximately 8.2 to 8.1 — a seemingly small change that represents a 25% increase in hydrogen ion concentration.

---

Chapter 7: Electrochemistry and Redox Geochemistry

Oxidation-Reduction Reactions in Natural Systems

The redox state of a natural environment — whether it is oxidizing or reducing — is one of the most important controls on element behaviour in geochemical systems. Oxidation is the loss of electrons, and reduction is the gain of electrons; these processes must always occur together (one species cannot gain electrons unless another loses them). In natural waters and sediments, the dominant redox couples involve oxygen, sulphur, nitrogen, iron, and manganese, and the predominance of particular oxidized or reduced species determines which minerals are stable, which metals are soluble or precipitated, and what microorganisms can metabolize.

The sequence in which oxidized species are consumed as an aquatic system becomes more reducing follows a predictable thermodynamic order dictated by the relative free energies of the relevant reactions. In aerobic respiration, microorganisms use O₂ as the terminal electron acceptor, oxidizing organic matter and yielding the most energy per mole of carbon. Once O₂ is depleted, nitrate (NO₃^-) is the next electron acceptor (denitrification), followed by Mn⁴⁺ in manganese oxyhydroxides (manganese reduction), then Fe³⁺ in iron oxyhydroxides (iron reduction), then SO₄^2- (sulphate reduction, producing H₂S), and finally CO₂ (methanogenesis, producing CH₄). Each step corresponds to a less negative Gibbs free energy yield for the microorganism involved, so the sequence is both thermodynamically predicted and observed in natural anoxic environments such as stratified lake bottoms, wetland soils, marine sediments, and poorly ventilated aquifers.

The quantitative framework for redox geochemistry is provided by electrochemistry. The tendency of a chemical species to accept or donate electrons is measured by its reduction potential \( E^\circ \), expressed in volts. For the reduction half-reaction \( \text{Ox} + ne^- \rightarrow \text{Red} \), the Nernst equation gives the actual electrode potential at non-standard conditions:

\[ E = E^\circ - \frac{RT}{nF} \ln \frac{[\text{Red}]}{[\text{Ox}]} = E^\circ - \frac{0.05916}{n} \log \frac{[\text{Red}]}{[\text{Ox}]} \text{ at 25°C} \]

where \( n \) is the number of electrons transferred and \( F = 96,485 \) C/mol is the Faraday constant. In geochemistry, the redox state of a system is described by the parameter Eh — the electrode potential of the system relative to the standard hydrogen electrode (SHE). Strongly oxidizing systems (e.g., acidic mine drainage in contact with oxygen) have high positive Eh values; strongly reducing systems (e.g., deep aquifers with H₂S) have large negative Eh values.

Eh-pH Diagrams (Pourbaix Diagrams)

Eh-pH diagrams (also called Pourbaix diagrams) display the thermodynamic stability fields of different species of a given element as a function of pH (horizontal axis) and Eh (vertical axis). These diagrams are constructed by solving the Nernst equation for the equilibria between each pair of adjacent species and plotting the resulting lines in Eh-pH space. The boundaries between stability fields are lines along which two species are in equilibrium; crossing a boundary represents a change in which species is thermodynamically dominant.

Consider the iron Eh-pH diagram in detail. Under oxidizing, near-neutral to basic conditions (high Eh, pH 5–10), ferric oxyhydroxide minerals (goethite FeOOH, ferrihydrite Fe(OH)₃, hematite Fe₂O₃) are stable and iron concentrations in solution are extremely low (often below 1 μg/L). Under reducing conditions (low Eh), Fe²⁺ becomes the dominant dissolved form, and iron concentrations can rise to tens or hundreds of mg/L. Under acidic and oxidizing conditions (low pH, high Eh) — the conditions found in acid mine drainage — Fe³⁺ can be the dominant dissolved species. The sharp transition from soluble Fe²⁺ to insoluble Fe³⁺ oxyhydroxide as Eh increases (or as pH increases) is why the mixing of anoxic groundwater with oxygenated surface water produces immediate precipitation of orange-red iron oxide — a common sight at spring seeps and iron-rich bogs in areas underlain by iron-bearing rocks or sulphide deposits.

The sulphur Eh-pH diagram reveals the large stability field of sulphate (SO₄^2-) under oxidizing conditions and sulphide (H₂S, HS^-) under reducing conditions, separated by a boundary that traverses the entire Eh-pH space with a steep slope. The transition between sulphate and sulphide involves an 8-electron change (sulphur goes from +6 to -2), which means the boundary is thermodynamically sharp but kinetically slow: sulphate and sulphide can coexist out of equilibrium in many natural systems because the rate of abiotic interconversion between them is very slow at low temperatures. This kinetic persistence is exploited by sulphate-reducing bacteria, which catalyse the thermodynamically favourable but kinetically sluggish reaction and play an essential role in the deep carbon cycle and in the geochemical cycling of sulphur through geological history.

---

Chapter 8: Stable Isotope Geochemistry

Principles of Isotopic Fractionation

Stable isotopes are isotopes that do not undergo radioactive decay. The most commonly studied stable isotope systems in geochemistry involve light elements: hydrogen (H/D = ¹H/²H), carbon (¹²C/¹³C), nitrogen (¹⁴N/¹⁵N), oxygen (¹⁶O/¹⁸O and ¹⁷O), and sulphur (³²S/³⁴S and higher masses). More recently, isotope systems of heavier elements — iron, magnesium, calcium, zinc, copper, molybdenum — have become important as “non-traditional” stable isotope tracers of redox processes, biological cycling, and weathering.

The key principle of stable isotope geochemistry is isotopic fractionation: the preferential partitioning of heavier or lighter isotopes between coexisting phases or between reactants and products in a chemical reaction. Fractionation occurs because isotopes of the same element have slightly different masses, and this mass difference causes them to behave slightly differently in chemical reactions, during physical processes (evaporation, condensation, diffusion), and during biological processes (metabolism, biosynthesis). There are two end-member types of fractionation: equilibrium fractionation (resulting from the thermodynamic preference for isotopes of a particular mass in a specific bonding environment) and kinetic fractionation (resulting from the different rates at which light and heavy isotopes react or diffuse).

Equilibrium isotopic fractionation decreases with increasing temperature — at very high temperatures (deep in the mantle), isotope fractionation between coexisting phases approaches zero. This temperature dependence provides a thermometer: measuring the isotope fractionation between two minerals in equilibrium (such as calcite and dolomite, or quartz and magnetite) gives the temperature at which they equilibrated. Kinetic fractionation, by contrast, produces larger fractionations than equilibrium and does not have the same temperature dependence. Biological fractionation of carbon isotopes during photosynthesis (where ¹²C is preferentially fixed relative to ¹³C, producing organic matter depleted in ¹³C by 20–30‰ relative to the CO₂ source) is a kinetic fractionation that has profoundly influenced the carbon isotope record throughout Earth history.

The magnitude of isotopic fractionation is expressed through the fractionation factor \( \alpha \), defined for the exchange of an isotope between phases A and B as:

\[ \alpha_{A-B} = \frac{R_A}{R_B} \]

where \( R \) is the ratio of the heavy to light isotope (e.g., ¹⁸O/¹⁶O) in each phase. The fractionation is also expressed through the delta notation, which reports isotope ratios as parts per thousand (‰) deviations from a standard:

\[ \delta^{18}O = \left(\frac{R_{sample} - R_{standard}}{R_{standard}}\right) \times 1000 \text{ ‰} \]

The relationship between \( \alpha \) and \( \delta \) values is approximately:

\[ \delta^{18}O_A - \delta^{18}O_B \approx 1000 \ln \alpha_{A-B} \equiv \Delta_{A-B} \]

for small fractionations (less than ~50‰). The \( \Delta \) notation thus represents the isotope difference between coexisting phases and is directly used as a thermometer or process indicator.

Oxygen Isotopes and the Meteoric Water Line

The oxygen isotope composition of water is a powerful tool in hydrology, climatology, and aqueous geochemistry. The key process is equilibrium isotopic fractionation between water vapour and liquid water during evaporation and condensation. Because the heavier ¹⁸O prefers slightly to remain in the liquid phase relative to the vapour (fractionation factor \( \alpha_{liq-vap} \approx 1.0092 \) at 25°C), evaporating water vapour becomes isotopically lighter (more negative \( \delta^{18}O \)) and successive precipitation events become progressively more depleted in ¹⁸O as the vapour mass moves inland and upward (Rayleigh distillation).

The meteoric water line (MWL), first described by Harmon Craig in 1961, is the global relationship between \( \delta^{18}O \) and \( \delta D \) (deuterium) in precipitation:

\[ \delta D = 8 \delta^{18}O + 10 \]

The slope of 8 arises from the ratio of fractionation factors for hydrogen and oxygen isotopes during condensation at equilibrium. The intercept of +10‰ (the global weighted average deuterium excess) results from kinetic fractionation during evaporation from the ocean surface under non-equilibrium conditions (relative humidity less than 100%). Continental rainfall becomes progressively depleted in heavy isotopes (more negative \( \delta^{18}O \)) with increasing distance from the coast (the continental effect), increasing altitude (the altitude effect, approximately −2.8‰ per 1000 m for \( \delta^{18}O \) in mid-latitudes), and decreasing temperature (the temperature effect).

---

Chapter 9: Partition Coefficients and Trace Element Geochemistry

The Concept and Significance of Partition Coefficients

The distribution of trace elements between coexisting phases — minerals, melts, and fluids — is one of the most powerful tools in igneous geochemistry for deciphering the conditions of magma genesis and differentiation. The partition coefficient (also called the distribution coefficient or Nernst distribution coefficient) \( D_i^{min/melt} \) for element \( i \) between a mineral and a coexisting melt is defined as:

\[ D_i^{min/melt} = \frac{C_i^{min}}{C_i^{melt}} \]

where \( C_i^{min} \) and \( C_i^{melt} \) are the concentrations of element \( i \) in the mineral and in the melt respectively. If \( D > 1 \), the element is said to be compatible with the mineral — it prefers to enter the mineral rather than remain in the melt. If \( D < 1 \), the element is incompatible and prefers the melt. Highly incompatible elements have \( D \ll 1 \) and are strongly concentrated in melts; compatible elements have \( D > 1 \) and are retained in the residual solid.

Partition coefficients are controlled by the crystal chemistry of the mineral (the size and charge of the available crystal site), the composition of the melt (which affects both the activity of the trace element and the structure of the silicate network), temperature, and pressure. They are not constants — they vary with all of these conditions — but for particular mineral-melt pairs under well-constrained conditions, they provide powerful quantitative models of trace element behaviour during igneous processes.

The lattice strain model (Blundy and Wood, 1994) provides a quantitative framework for predicting how partition coefficients vary with ionic radius for a given crystal site. The model treats the substitution of a trace element into a crystal site as an elastic deformation of the surrounding lattice, and the energy cost of this deformation increases with the square of the mismatch between the trace element radius and the ideal site radius. The resulting Onuma diagram — a plot of \( \ln D_i \) against ionic radius — shows a parabolic relationship with a maximum at the “optimal” radius \( r_0 \) for the site. This model has been extraordinarily successful at predicting partition coefficients for the REE (which all carry the same 3+ charge and vary only in ionic radius) in a wide range of minerals including clinopyroxene, garnet, zircon, and amphibole.

Batch Melting, Fractional Melting, and Crystallisation

When a rock begins to melt, trace elements partition themselves between the residual solid and the melt according to their partition coefficients. The two end-member models for partial melting are batch (equilibrium) melting, in which the melt remains in contact with the residue until it is fully segregated, and fractional melting, in which each infinitesimal increment of melt is immediately removed from the system before it can re-equilibrate.

For batch melting, the concentration of a trace element in the melt is:

\[ C_L = \frac{C_0}{D_0 + F(1 - D_0)} \]

where \( C_L \) is the concentration in the liquid, \( C_0 \) is the initial concentration in the source rock, \( F \) is the melt fraction, and \( D_0 = \sum_j p_j D_i^{j/melt} \) is the bulk distribution coefficient (sum over all mineral phases weighted by their proportions). For highly incompatible elements (\( D_0 \ll F \)), \( C_L \approx C_0/F \): the melt concentration is simply the source concentration divided by the melt fraction. For compatible elements (\( D_0 \gg F \)), \( C_L \approx C_0/D_0 \): the melt concentration is controlled by the mineral, nearly independent of melt fraction.

For fractional melting, the instantaneous melt composition is:

\[ C_L = \frac{C_0}{D_0} (1 - F)^{(1/D_0 - 1)} \]

Fractional melting produces more extreme fractionation of incompatible elements than batch melting for the same melt fraction. In practice, melting in the Earth’s mantle probably lies between these two end-members, but the batch melting model is often a good approximation for trace element modelling of MORB and OIB magmas.

---

Chapter 10: Environmental Geochemistry and Resource Applications

Water Quality, Contamination, and Drinking Water Advisories

The intersection of geochemistry with public health is nowhere more apparent than in the study of drinking water quality. The natural geochemical cycle produces waters that are suitable for consumption across most of the Earth, but in specific geological settings, naturally occurring elements can be present at concentrations that pose risks to human health. The two most important naturally occurring drinking water contaminants globally are arsenic and fluoride, both of which occur at elevated concentrations in groundwaters derived from specific rock types.

Arsenic occurs naturally in many rock types, particularly sedimentary rocks containing pyrite (FeS₂), sulphide ore deposits, and volcanic rocks. The most common geochemical mechanism for arsenic mobilisation in groundwater is reductive dissolution of iron oxyhydroxides: in reducing, organic-rich sediments, Fe³⁺ in iron oxyhydroxide minerals is reduced to Fe²⁺ by microbially-mediated organic matter oxidation, and as the oxyhydroxide dissolves, the arsenate (As⁵⁺) adsorbed onto its surface is released into solution. This mechanism is responsible for the catastrophic arsenic crisis in the Bengal Basin, where 30–80 million people are exposed to arsenic concentrations exceeding the WHO guideline of 10 μg/L (some wells exceed 1,000 μg/L). In Canada, elevated arsenic in First Nations community water systems in Ontario and Manitoba has similar origins — reducing shallow aquifers in glacial sediments derived from arsenic-bearing bedrock.

Fluoride provides a contrasting example: it is mobilised from volcanic rocks and fluorite (CaF₂) under high-pH, alkaline conditions, particularly in the East African Rift Valley and parts of China and India. At concentrations below 1 mg/L, fluoride strengthens tooth enamel and reduces dental caries; above 1.5 mg/L (WHO guideline), it causes dental fluorosis (mottled enamel); above 4 mg/L, it causes skeletal fluorosis, a debilitating condition affecting millions of people in India and East Africa. The geochemical control on fluoride concentration in groundwater is the solubility product of fluorite, which constrains dissolved fluoride if calcium concentrations are high, but in low-calcium waters (common in highly weathered tropical soils), fluoride can accumulate to very high levels. Understanding the geochemical controls on fluoride concentrations is essential for identifying affected communities and designing appropriate treatment systems.

Ore-Forming Processes: A Geochemical Perspective

The formation of ore deposits — concentrations of economically valuable elements or minerals that are sufficiently rich to be extracted profitably — is fundamentally a geochemical problem. Nature must concentrate an element by a factor of tens to thousands of times its average crustal abundance to produce an ore deposit. Understanding how this concentration happens requires understanding the thermodynamic and kinetic controls on element solubility and transport in hydrothermal fluids, on mineral saturation and precipitation, and on the coupling of geochemical systems to tectonic and magmatic processes.

Gold provides a striking example of extreme geochemical concentration. The average crustal abundance of gold is approximately 2 ppb (parts per billion by weight). A gold ore deposit must typically contain at least 1–3 grams per tonne (1–3 ppm) to be economic, representing a concentration factor of 500 to 1500 relative to average crust. In some epithermal gold deposits, gold concentrations may reach 10–100 g/tonne — concentration factors of 5,000–50,000. Gold is transported in hydrothermal fluids primarily as gold-chloride complexes (AuCl₂^-) at high temperatures (above 350°C) and as gold-bisulphide complexes (Au(HS)₂^-) at lower temperatures (150–350°C). Precipitation occurs in response to cooling, boiling (which removes H₂S as steam), oxidation (which destroys bisulphide complexes), or mixing with more oxidized or dilute fluids. The geochemical environments where these changes occur efficiently — epithermal systems at shallow depths, orogenic gold systems along faults, Carlin-type deposits in calcareous sediments — are the exploration targets of the gold mining industry.

Environmental Geochemistry of Mining and Acid Mine Drainage

Mining operations release elements from solid minerals into the environment, often at rates and concentrations that far exceed natural weathering fluxes. The most widespread environmental impact of mining is acid mine drainage (AMD), which occurs when sulphide minerals (primarily pyrite, FeS₂) are exposed to oxygen and water. The oxidation of pyrite is thermodynamically extremely favourable (\( \Delta G^\circ \approx -1,440 \) kJ/mol) but kinetically slow under purely abiotic conditions; the process is dramatically accelerated by acidophilic, iron-oxidizing bacteria (primarily Acidithiobacillus ferrooxidans) once the pH drops below about 4. The overall oxidation reaction can be written as:

\[ \text{FeS}_2 + \frac{15}{4}\text{O}_2 + \frac{7}{2}\text{H}_2\text{O} \rightarrow \text{Fe(OH)}_3 + 2\text{H}_2\text{SO}_4 \]

The sulphuric acid produced dissolves carbonate minerals (if present), lowering pH further; once the neutralising capacity of the rock is exhausted, pH plummets to 1–3, mobilising a wide range of metals (Fe, Zn, Cu, As, Pb, Cd, Mn) from the waste rock and tailings. These metals precipitate as oxyhydroxides when the acidic water mixes with more neutral stream water, producing the characteristic orange “yellow boy” that is visible in streams draining abandoned metal mines throughout Canada, the USA, Spain, and elsewhere. The long-term legacy of AMD can persist for centuries after mining operations cease — lake sediment records show elevated lead concentrations from Roman-era silver mining at Rio Tinto in Spain and from medieval copper mining in Falun, Sweden. Addressing this legacy requires a combination of source control (preventing oxidation by water exclusion or subaqueous tailings disposal), passive treatment (alkaline substrates that neutralise acidity), and active treatment (lime neutralisation of effluent, sulphate-reducing bioreactors), all of which require a thorough understanding of the underlying geochemistry.

The geochemical principles covered in EARTH 221 — from elemental abundances set by stellar nucleosynthesis, through crystal chemistry governing mineral composition, through thermodynamics controlling mineral stability, through radioactive decay enabling geochronology, through aqueous chemistry determining water quality, through redox geochemistry controlling metal mobility, through stable isotopes tracing water sources, and through partition coefficients illuminating magmatic processes — form an integrated framework for understanding the chemistry of the Earth system at all scales. This framework is essential preparation for advanced courses in igneous and metamorphic petrology, hydrogeology, ore deposit geology, and environmental earth sciences.