EARTH 331: Volcanology and Igneous Petrology

Carson Kinney

Estimated study time: 36 minutes

Table of contents

Sources and References

Primary textbook — Winter, J.D. (2010). Principles of Igneous and Metamorphic Petrology (2nd ed.). Pearson Prentice Hall. Supplementary texts — Nesse, W.D. (2012). Introduction to Optical Mineralogy (4th ed.). Oxford University Press. — Klein, C. & Philpotts, A. (2017). Earth Materials: Introduction to Mineralogy and Petrology (2nd ed.). Cambridge University Press. — Frost, B.R. & Frost, C.D. (2019). Essentials of Igneous and Metamorphic Petrology (2nd ed.). Cambridge University Press. Online resources — GEOROC database (georoc.eu) — EarthChem Portal (earthchem.org) — USGS Volcano Hazards Program (usgs.gov/programs/VHP) — PetDB: Petrological Database of the Ocean Floor

Chapter 1: Igneous Rock Classification and the Igneous System

Introduction to Igneous Petrology

Igneous rocks — those solidified from molten or partially molten material — constitute the overwhelming bulk of Earth’s mass and the foundation upon which all other rock types are built. Oceanic crust is almost entirely basaltic; the continental crust, though more compositionally diverse, owes its existence to billions of years of magmatic addition from below. Volcanology and igneous petrology together seek to understand the generation, transport, storage, and solidification of magmas: where melts originate, what governs their compositions, how they ascend through the lithosphere, and what textures and mineralogy they produce upon crystallization.

The classification of igneous rocks rests on two axes: texture (governed by the rate of cooling and crystallization, which reflects the depth and setting of emplacement) and composition (governed by the source material and the processes — partial melting, fractional crystallization, assimilation — that shape magma chemistry). Together, these axes allow geologists to read the history of a rock from its hand-sample characteristics and to link surface geology to deep Earth processes.

Texture: Phaneritic versus Aphanitic Rocks

The rate at which a magma cools exerts the dominant control on grain size. When magma crystallizes slowly — over thousands to millions of years in the deep crust or upper mantle — crystals have time to grow to millimetre to centimetre scale, producing a phaneritic (coarse-grained) texture in which individual mineral grains are visible to the naked eye. Such rocks form intrusive (or plutonic) bodies: batholiths, stocks, sills, dykes, and lopoliths. In contrast, when lava erupts at the surface or when magma is injected into shallow, rapidly cooled settings, crystallization is rapid and grain sizes are fine to microscopic, producing aphanitic (fine-grained) textures in volcanic and hypabyssal rocks.

Phaneritic texture describes igneous rocks in which all mineral grains are large enough to be identified by eye without magnification (grain size typically > 1 mm). This texture indicates slow cooling in a plutonic environment. Aphanitic texture describes rocks in which crystals are too small to be distinguished without a microscope (grain size < 1 mm), indicating rapid cooling in a volcanic or shallow intrusive setting. Porphyritic texture involves two distinct grain-size populations: larger crystals (phenocrysts) set in a finer-grained or glassy groundmass (matrix), indicating a two-stage cooling history — an initial slow-cooling phase during magma ascent or storage that grew the phenocrysts, followed by more rapid cooling at or near the surface.

Additional igneous textures document more specific cooling histories. Glassy texture (obsidian, pitchstone) results from quenching so rapid that crystallization is suppressed altogether, producing an amorphous solid. Vesicular texture records the degassing of volatile-saturated magma: gas bubbles nucleate as pressure decreases during ascent or eruption, and if solidification overtakes bubble escape, the voids (vesicles) are preserved. When vesicles are later filled by secondary minerals such as calcite, quartz, or zeolites, the structures are called amygdales. Pyroclastic textures — found in tephra, welded tuffs, and ignimbrites — reflect fragmentation of magma during explosive eruptions; they include vitric (glassy shard), crystal, and lithic components embedded in an ash matrix.

Under the petrographic microscope, igneous textures become even more informative. Ophitic texture — characteristic of dolerites and some gabbros — involves large pyroxene crystals enclosing smaller, randomly oriented plagioclase laths. Interlocking texture (granitic or granular) results from simultaneous crystallization of major phases. Poikilitic texture involves large crystals (oikocrysts) enclosing smaller crystals of other phases that nucleated earlier. Graphic texture — occurring in pegmatites and granites — shows regular, parallel intergrowths of quartz and alkali feldspar resembling cuneiform script in cross-section, interpreted as resulting from simultaneous crystallization near the cotectic minimum.

Compositional Classification

Igneous rock composition is most practically expressed in terms of silica (SiO₂) content and the relative abundances of major minerals. The IUGS (International Union of Geological Sciences) classification, which is the international standard, defines rock names primarily on the basis of modal mineralogy (the actual volume percentages of minerals present) for phaneritic rocks, using QAPF diagrams (ternary plots of Quartz – Alkali feldspar – Plagioclase – Feldspathoid abundance). For aphanitic rocks and glasses where mineral modes cannot be determined, chemical or normative classification is used.

The TAS (Total Alkali versus Silica) diagram plots the sum of Na₂O + K₂O (wt%) against SiO₂ (wt%) for volcanic rocks and is the primary chemical classification scheme for fine-grained and glassy rocks. Field boundaries define names including basalt (SiO₂ 45–52 wt%, low alkalis), andesite (SiO₂ 52–63 wt%), dacite (SiO₂ 63–68 wt%), and rhyolite (SiO₂ > 68 wt%), as well as alkaline series rocks (trachybasalt, phonolite, trachyte) that plot above the alkaline-subalkaline divide.

Broadly, igneous rocks are classified as ultramafic (SiO₂ < 45 wt%, dominated by olivine and pyroxene — peridotite, dunite, komatiite), mafic (SiO₂ 45–52 wt%, plagioclase + pyroxene ± olivine — basalt, gabbro), intermediate (SiO₂ 52–63 wt%, plagioclase + pyroxene/hornblende — andesite, diorite), felsic (SiO₂ > 63 wt%, quartz + feldspar ± muscovite — rhyolite, dacite, granite, granodiorite). This scheme reflects a fundamental truth about magma genesis: mafic melts typically originate from the mantle through high-degree partial melting, while felsic melts are produced either by low-degree melting of mafic or metasedimentary crustal sources or by extensive differentiation of mafic magmas.

Chapter 2: Magma Generation, Migration, and Emplacement

The Physics of Melt Generation

Magma is generated wherever a portion of the solid Earth is raised above its solidus — the temperature at which the first melt appears. Because temperature increases with depth (the geothermal gradient), while the solidus temperature of mantle peridotite also increases with pressure (though more steeply than typical geotherms under normal conditions), the mantle is mostly solid under steady-state conditions. Melting can be induced by three processes that shift the system across the solidus: decompression melting, flux melting, and heat addition.

Decompression melting occurs when mantle material rises faster than heat can escape, following a nearly adiabatic pressure-temperature path. As pressure decreases, the solidus moves to lower temperatures faster than the temperature of the rising mantle decreases, so the adiabat crosses the solidus and melting begins. This process operates at mid-ocean ridges (where passive mantle upwelling fills the gap left by diverging plates) and at mantle plumes (where hot, buoyant mantle rises from depth as an active diapir). The degree of partial melting — the fraction of the source rock that becomes melt — determines how much melt is produced and, crucially, what chemical composition it has.

Flux melting (or fluid-induced melting) occurs when the addition of water (or other volatile components) to hot mantle peridotite lowers the solidus temperature sufficiently to induce melting at pressures and temperatures where the rock would otherwise remain solid. Water is a powerful solidus depressant: at 3 GPa (approximately 90 km depth), water addition can lower the peridotite solidus by 100–200°C. In subduction zones, hydrated minerals in the subducting oceanic crust and sediments break down as the slab descends and heats, releasing aqueous fluids that migrate upward into the overlying mantle wedge, inducing melting above the slab.

Magma Ascent and Intrusive Geometries

Once generated, magma rises toward the surface through a combination of buoyancy (melt densities are typically 2.4–2.8 g/cm³, less than most solid mantle rocks at ~3.3 g/cm³), hydraulic head from the pressure of overlying solid rock, and the propagation of fluid-filled fractures (dykes). The dynamics of dyke propagation are analogous to hydrofracture: magma propagates as a fluid-filled crack when the magmatic excess pressure exceeds the tensile strength of the country rock. Dykes thus tend to propagate perpendicular to the minimum principal stress, which in extensional regimes is horizontal — producing vertical dykes — and in compressional regimes may be vertical — producing horizontal intrusions (sills).

The full diversity of igneous intrusive bodies reflects the interplay of magma supply rate, country rock rheology, structural controls, and the depth of emplacement. Batholiths — the largest intrusive bodies, with surface exposures exceeding 100 km² — are composite bodies constructed by multiple intrusive episodes and represent the eroded roots of ancient magmatic arcs. Laccoliths are lens-shaped intrusions where magma domed the overlying strata upward rather than pushing laterally as a sill. Ring dykes and cone sheets occur in collapsed caldera environments where magma exploits arcuate fault systems. The geometry of an intrusion thus encodes information about the stress field and the mechanical properties of the surrounding crust at the time of emplacement.

Surface Emplacement: Lava Flows and Volcanic Architecture

The morphology of volcanic edifices and the architecture of volcanic deposits depend on the physical properties of the erupted magma — above all on viscosity and volatile content, which together govern whether eruptions are effusive (lava-producing) or explosive (tephra-producing). Viscosity, in turn, depends primarily on silica content (higher SiO₂ = more polymerized melt network = higher viscosity), temperature (lower T = higher viscosity), and dissolved water content (water depolymerizes the silicate network, dramatically reducing viscosity).

Basaltic lavas, with viscosities of \( 10^1 - 10^3 \) Pa·s (comparable to honey to tar), flow readily and can travel tens to hundreds of kilometres from their vents, building broad, gently sloping shield volcanoes or extensive flood basalt provinces. Two morphological types of basaltic lava flow are common: pahoehoe (smooth, ropy, or billowy surface, formed from low-viscosity lava flowing in insulated lava tubes) and aa (rough, clinkery, jagged surface, formed from more viscous, faster-moving lava that breaks its chilling crust). Silicic lavas, with viscosities of \( 10^8 - 10^{12} \) Pa·s, flow sluggishly and tend to form short, thick lava domes that grow endogenously (internally) or exogenously (by surface extrusion), with frequent gravitational collapse that generates pyroclastic density currents.

Chapter 3: Crystallization Mechanics, Microstructures, and Phase Diagrams

Nucleation and Crystal Growth

The solidification of a magma into an igneous rock is governed by the processes of crystal nucleation and growth, which are in turn controlled by the degree of undercooling (the difference between the actual temperature of the melt and the equilibrium liquidus temperature at which the phase would first crystallize). Nucleation — the formation of a stable embryonic crystal — has an energy barrier because the creation of a new surface between crystal and melt costs energy (surface free energy), which must be overcome by the free energy gain of forming the ordered solid phase. At small undercooling, the nucleation rate is low, and the few nuclei that form have time to grow large — producing coarse-grained rocks. At large undercooling, nucleation is rapid, many small crystals form simultaneously, and grain sizes are fine. At extreme undercooling (quenching), nucleation fails entirely and glass forms.

The liquidus is the temperature above which a magma is entirely molten (for a given composition and pressure). The solidus is the temperature below which the system is entirely crystalline. Between the liquidus and solidus, crystals and melt coexist; the proportion of solid increases as temperature falls from liquidus to solidus. In multi-component systems, the liquidus and solidus are surfaces in composition-temperature-pressure space, not simple points.

Crystal growth occurs by attachment of ions or molecular units from the melt onto existing crystal surfaces. Growth rates are fastest at moderate undercoolings (a few tens of degrees below the liquidus), where diffusion in the melt is still rapid enough to supply material to the growing crystal face and the thermodynamic driving force for crystallization is substantial. The morphology of crystals — from equant (blocky) to acicular (needle-like) to dendritic — records the ratio of diffusion rate to growth rate during crystallization: dendritic forms reflect rapid growth in diffusion-limited conditions, while equant forms grow under near-equilibrium conditions with slow cooling.

Binary and Ternary Phase Diagrams

Phase equilibrium diagrams — the graphical representation of the stable phases in a chemical system as functions of temperature, pressure, and composition — are among the most powerful tools in igneous petrology for understanding crystallization paths, the evolution of melt composition, and melting behaviour.

The binary eutectic system with no solid solution (illustrated by the system diopside-anorthite at 1 atm) shows that a melt of any composition between the two end members begins crystallizing the phase on whose side of the eutectic composition it plots. As crystallization proceeds and that phase is removed from the melt, the melt composition evolves toward the eutectic point — the temperature and composition at which the two phases crystallize together until the melt is exhausted. The eutectic temperature is the lowest melting point in the binary system. The opposite process — melting — begins at the eutectic when temperature rises from below; the first melt always has the eutectic composition, regardless of the bulk composition of the starting solid mixture, so long as both phases are present.

The Lever Rule. In a two-phase region of a phase diagram (melt + crystal coexisting), the proportion of each phase at a given bulk composition and temperature is determined by the lever rule: \[ \frac{\text{mass of phase A}}{\text{mass of phase B}} = \frac{l_B}{l_A} \]

where \( l_A \) and \( l_B \) are the distances from the bulk composition to the compositions of phases B and A, respectively, along the tie line connecting the two phase compositions. This rule follows directly from mass balance.

The system plagioclase (anorthite-albite) is a classic example of complete solid solution: the minerals crystallize across a continuous compositional range from Ab₀An₁₀₀ (pure anorthite) to Ab₁₀₀An₀ (pure albite), with liquidus and solidus curves forming a lens-shaped two-phase field. Plagioclase crystallizing in equilibrium with a melt is always more anorthite-rich than the melt, so as crystallization proceeds and anorthite is preferentially removed, the melt evolves toward the albite end, and successive generations of plagioclase crystals become more albite-rich. The resulting chemical zonation of plagioclase crystals — calcic cores grading to sodic rims — is one of the most common textural features of igneous rocks and records the progressive evolution of the magma during crystallization.

Ternary phase diagrams for systems such as Di-An-Ab (diopside-anorthite-albite) allow more realistic representation of natural magma systems and introduce the concept of cotectic lines (boundaries between the stability fields of two minerals, along which two phases crystallize simultaneously as the melt evolves to lower temperatures) and ternary eutectic/peritectic points. The 1-atmosphere basalt tetrahedron — the system Fo-Di-An-Qz (forsterite-diopside-anorthite-quartz, the essential components of basaltic magmas) — illustrates how basaltic melt compositions at the cotectic evolve toward the thermal divide that separates the basalt field from the granite field: a major reason why intermediate magma compositions are less common than basaltic and granitic end members.

Volatiles and Bulk Composition Effects

Dissolved volatile components — primarily H₂O, CO₂, SO₂, and halogens — profoundly modify magma properties and phase equilibria, even at low concentrations (tenths of a percent by weight). Water is the most geochemically important volatile in most crustal and subduction-zone settings: it depresses the solidus and liquidus of silicate magmas (expanding the P-T stability field of the melt), reduces melt viscosity by depolymerizing the silicate network, and stabilizes the crystallization of hydrous minerals (hornblende, biotite, phlogopite) that would not be present in dry systems.

The water solubility in silicate melts increases approximately as the square root of water partial pressure (or as pressure to the power 0.5–0.7), meaning that at lower pressures — as magma ascends toward the surface — water becomes progressively less soluble and exsolves as a vapour phase. This volatile exsolution is the primary driver of explosive volcanic eruptions: the nucleation and rapid expansion of water-rich bubbles (and CO₂-rich bubbles at greater depths) fragments the viscous magma into glass shards and pumice, launching pyroclastic material into the atmosphere.

Carbon dioxide has a lower solubility than water in silicate melts and exsolves at greater depths (higher pressures), making it the dominant volatile species in CO₂-rich magmas at mantle depths (kimberlites, carbonatites). The interplay of CO₂ and H₂O solubilities determines the depth at which degassing begins, the rate of bubble nucleation and growth during ascent, and ultimately the eruptive style.

Chapter 4: Igneous Field Relationships

Reading the Geological Record of Igneous Activity

The field relationships of igneous bodies — their contacts with surrounding rocks, cross-cutting relationships, intrusive versus extrusive setting, relative age relationships, and geometric form — provide evidence of the order of geological events that cannot be reconstructed from geochemistry or petrology alone. The fundamental principles of stratigraphy and relative chronology apply: cross-cutting relationships establish temporal order (a dyke that cuts a batholith is younger than the batholith), inclusions establish that the rock containing them is younger than the inclusions, and the chilled margin of an intrusion against its host rock demonstrates that the intrusion post-dates the host.

Contact metamorphism in the rocks immediately surrounding an intrusion provides information about the thermal regime of emplacement: hornfels adjacent to a small, shallow stock records temperatures of 400–600°C at the contact; the development of wollastonite (from calcite + quartz reaction) in a calc-silicate hornfels indicates temperatures exceeding ~450°C at low PCO₂. The width of the contact aureole (the zone of thermally metamorphosed rock) scales with the size of the intrusion (larger bodies have more heat to deliver), the thermal conductivity of the host rock, and the emplacement depth (deeper intrusions are in a hotter environment and produce narrower aureoles relative to their size).

Volcanic field relationships add another dimension: the stratigraphy of lava flows and pyroclastic deposits records successive eruptive episodes, with younger flows overlying older ones (the principle of superposition). Unconformities within volcanic sequences record intervals of erosion or non-deposition; the presence of paleosols (ancient soils) between flow units indicates sufficiently long inter-eruptive intervals for soil formation. Ignimbrites — extensive sheets of welded or unwelded pyroclastic deposits from large explosive eruptions — can be traced over hundreds of kilometres and correlated chemically or isotopically (geochemical fingerprinting) when physical tracing fails.

Chapter 5: Earth Formation and Geochemical Context

Planetary Differentiation and the Building Blocks of Magmas

Earth formed approximately 4.56 Ga by accretion of planetesimals in the inner solar system. Within the first tens of millions of years, sufficient energy from radioactive decay (primarily ²⁶Al and ⁶⁰Fe) and accretionary impact heating caused widespread melting and differentiation: metallic iron, nickel, and siderophile elements sank to form the core, while silicate minerals segregated into the mantle. This differentiation established the compositional stratification — core, mantle, crust — that determines the source regions and compositional envelopes of all terrestrial magmas.

The primitive mantle (or bulk silicate Earth) composition — the composition of the mantle before extraction of the continental crust — serves as the reference point from which the compositions of all igneous rocks are measured. It is estimated by combining peridotite xenolith data, ophiolite studies, and mass-balance calculations accounting for the continental crust. The primitive mantle is estimated to contain approximately 45.0 wt% SiO₂, 37.8 wt% MgO, 8.0 wt% FeO, 3.5 wt% Al₂O₃, 3.1 wt% CaO, 0.35 wt% Na₂O, 0.03 wt% K₂O, plus minor TiO₂, Cr₂O₃, MnO, and NiO.

Modern mantle geochemistry recognizes that the mantle is compositionally heterogeneous, preserving distinct geochemical reservoirs with different isotopic and trace element signatures: MORB source mantle (DMM, Depleted MORB Mantle) — the depleted upper mantle from which mid-ocean ridge basalts are derived; HIMU (high-μ = high ²³⁸U/²⁰⁴Pb) — ancient recycled oceanic crust that has developed high Pb isotopic ratios; EM1 and EM2 (Enriched Mantle types 1 and 2) — mantle domains variably enriched by ancient subducted sediments and/or delaminated lower continental crust; and FOZO/C — a primordial undegassed deep mantle component.

Chapter 6: Igneous Geochemistry and Isotopes

Trace Element Behaviour During Melting and Crystallization

The behaviour of trace elements — those present at concentrations below ~1 wt% — provides some of the most sensitive fingerprints of igneous processes because, unlike major elements that are largely controlled by stoichiometry, trace elements can be dramatically fractionated between melt and coexisting minerals depending on their ionic size and charge.

The partition coefficient (or distribution coefficient) \( D_i^{\text{mineral/melt}} \) for element \( i \) between a mineral and coexisting melt is defined as: \[ D_i = \frac{C_i^{\text{mineral}}}{C_i^{\text{melt}}} \]

An element with \( D_i \gg 1 \) is compatible — it is preferentially incorporated into the crystal relative to the melt — and is depleted in the melt as crystallization proceeds. An element with \( D_i \ll 1 \) is incompatible — it is excluded from the crystal lattice and concentrates in the melt. Highly incompatible elements (D ≪ 0.01) include the large ion lithophile elements (LILE: Rb, Cs, Ba, Sr, K, Pb²⁺) and the high field strength elements (HFSE: Zr, Nb, Ta, Hf, Ti, Th, U).

For batch (equilibrium) partial melting — in which melt remains in equilibrium with the residual solid until it is extracted — the concentration of a trace element in the melt \( C_L \) relative to the source concentration \( C_0 \) is:

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

where \( F \) is the melt fraction and \( D \) is the bulk partition coefficient (weighted average of mineral partition coefficients according to modal abundances in the source). This equation shows that incompatible elements (small D) are dramatically enriched in small-degree melts: for \( D = 0.01 \) and \( F = 0.01 \) (1% melting), \( C_L/C_0 \approx 50 \), meaning the melt contains 50 times the source concentration. Compatible elements (large D) are depleted in the melt relative to the source.

Spider Diagrams and REE Patterns

Rare earth elements (REE) — the lanthanide series from La (Z=57) to Lu (Z=71) — are particularly valuable as petrogenetic tracers because they are all trivalent under most geological conditions (with the exception of Ce and Eu which can be Ce⁴⁺ and Eu²⁺ under specific redox conditions), their ionic radii decrease systematically from La³⁺ (1.03 Å) to Lu³⁺ (0.86 Å), and their partition coefficients vary systematically and predictably with ionic radius for any given mineral. The result is that REE patterns — element concentrations normalized to primitive mantle or chondrite and plotted as a function of atomic number — are smooth, diagnostic curves whose shape encodes the mineralogy of the residual source and the degree of partial melting.

A negative Eu anomaly in a REE pattern (Eu depleted relative to adjacent Sm and Gd) indicates that plagioclase was present either as a residual phase during melting (plagioclase is the main host for Eu²⁺, which substitutes for Ca²⁺ in the plagioclase structure) or was removed by fractional crystallization; a positive Eu anomaly indicates plagioclase accumulation. The absence of a Eu anomaly in oceanic basalts generated at depths greater than ~40–50 km (where plagioclase is not stable, having converted to garnet or spinel) is thus diagnostic of a deep source.

Primitive mantle-normalized spider diagrams extend trace element comparison across a wider range of elements (typically ordered by increasing compatibility from left to right: Cs, Rb, Ba, Th, U, K, Nb, Ta, La, Ce, Pb, Pr, Nd, Sr, Sm, Zr, Hf, Eu, Ti, Gd, Dy, Y, Er, Yb, Lu). These diagrams reveal characteristic geochemical signatures of different tectonic settings: arc basalts show pronounced negative Nb-Ta anomalies (depletion in high field strength elements relative to LILE, resulting from the retention of these elements in residual rutile in the subducting slab), while MORB and OIB have smooth or flat patterns without significant HFSE depletions. These anomalies serve as geochemical fingerprints that allow ancient subduction zones to be identified in Precambrian terranes.

Radiogenic Isotopes as Tracers of Source Composition

Radiogenic isotope systems — Rb-Sr, Sm-Nd, Lu-Hf, U-Pb, Re-Os — serve a dual purpose in igneous petrology: geochronology (dating when rocks crystallized) and isotopic tracing (identifying and characterizing magma source regions). The isotopic compositions of Sr, Nd, Hf, and Pb in volcanic rocks are inherited from their mantle or crustal source and are not significantly modified during magmatic differentiation (because isotopic ratios of a given element are not fractionated by melting or crystallization). They therefore provide direct windows into the compositions and histories of mantle and crustal source regions.

The Sr-Nd Isotope Diagram. Plotting ⁸⁷Sr/⁸⁶Sr against ¹⁴³Nd/¹⁴⁴Nd for a global dataset of volcanic rocks reveals a striking negative correlation (the "mantle array"): depleted reservoirs (MORB, DMM) have low ⁸⁷Sr/⁸⁶Sr (0.7022–0.7030) and high ¹⁴³Nd/¹⁴⁴Nd (reflecting ancient depletion of Rb relative to Sr, and enrichment of Nd relative to Sm), while enriched reservoirs (EM1, EM2) and crustal rocks have high ⁸⁷Sr/⁸⁶Sr and low ¹⁴³Nd/¹⁴⁴Nd. The position of a volcanic rock sample on this diagram constrains the nature and proportions of its source components. Mixing between depleted mantle and subducted sediment (high ⁸⁷Sr/⁸⁶Sr, low ¹⁴³Nd/¹⁴⁴Nd) produces hyperbolic mixing trends consistent with observed arc basalt compositions.

Chapter 7: Magmatism at Tectonic Boundaries

Mid-Ocean Ridges and MORB

Mid-ocean ridges form a continuous, 65,000 km network of divergent plate boundaries where decompression melting of upwelling mantle produces the oceanic crust that floors the world’s ocean basins. Mid-ocean ridge basalts (MORB) are the most volumetrically abundant type of volcanic rock on Earth — approximately 20 km³ are produced annually — and represent the product of relatively high-degree (10–25%) partial melting of depleted peridotite (DMM) at depths of 60–120 km.

MORB compositions fall within a restricted range: SiO₂ 48–52 wt%, MgO 7–12 wt%, FeO 8–12 wt%, CaO 10–13 wt%, Al₂O₃ 14–17 wt%, Na₂O 2–3 wt%, TiO₂ 1–2 wt%, with uniformly low K₂O, P₂O₅, and other incompatible element concentrations relative to OIB (ocean island basalts). Normal MORB (N-MORB) is strongly depleted in incompatible elements; enriched MORB (E-MORB) shows higher concentrations of incompatible elements and overlaps in some characteristics with OIB, indicating contamination by or proximity to mantle plume material. The global array of MORB compositions records variations in mantle potential temperature, extent of melting, and the heterogeneity of the depleted mantle source.

Ocean Island Basalts and Mantle Plumes

Ocean island basalts (OIB) erupt at oceanic islands — chains like Hawaii, the Canaries, Reunion, Iceland, and the Galapagos — that are interpreted as surface expressions of mantle plumes: quasi-cylindrical upwellings of anomalously hot mantle material rising from deep within the mantle (possibly from the core-mantle boundary). Plume-generated magmas are typically more enriched in incompatible elements, volatiles, and isotopically distinct components (HIMU, EM1, EM2) than MORB, reflecting a source that contains recycled oceanic lithosphere or other exotic components mixed into the lower mantle over billions of years.

The Hawaiian volcanic chain exemplifies the classic plume model: a stationary hotspot above a deep-seated plume melts through the overlying Pacific plate as it moves northwestward at ~8 cm/year, producing a systematic age-progressive chain of volcanoes, with active volcanism at the present-day position of the hotspot (Kilauea, Mauna Loa) and increasingly older, eroded seamounts and atolls extending toward the northwest (the Emperor Seamount chain turning northwest-trending, interrupted by the Hawaii-Emperor bend at ~50 Ma).

Island Arc and Continental Arc Magmatism

Island arcs and continental arcs form above subduction zones where oceanic lithosphere descends into the mantle. The compositional diversity of arc magmas — ranging from basalt through andesite to rhyolite/dacite — and their distinctive geochemical signatures (high LILE/HFSE ratios, elevated Ba/Nb, La/Nb, negative Nb-Ta anomalies in primitive-mantle-normalized diagrams) reflect the interplay of three source/process contributions: (1) partial melting of the metasomatized mantle wedge (enriched by slab-derived fluids), (2) fractionation of primary arc basalts during ascent and storage in the deep crust, and (3) assimilation of continental crust (in continental arcs).

Calc-alkaline magmatism — the characteristic magmatic series of subduction zones — is defined by increasing SiO₂ with decreasing FeO (relative to MgO), in contrast to the tholeiitic trend of increasing Fe up to a maximum before silica increases. The calc-alkaline trend is attributed to early crystallization of magnetite (a Fe³⁺-bearing oxide) in water-rich arc magmas, which removes iron from the melt before it can be significantly enriched. High water contents in arc magmas (derived from the slab) also stabilize hornblende, whose crystallization further reduces the FeO/MgO ratio of the evolving melt.

Continental Rifts, Stable Cratons, and Archean Magmatism

Where continents are being pulled apart by extensional tectonics — as in the East African Rift System — decompression melting of the underlying lithosphere and asthenosphere generates continental rift magmatism. Rift volcanics are compositionally diverse, ranging from mantle-derived basalts and picrites (in areas of thin lithosphere and high mantle potential temperature) to highly evolved, alkali-rich phonolites, trachytes, and peralkaline rhyolites (produced by extensive differentiation of the primary melts within crustal magma chambers). The unique volcanoes of the East African Rift — including the explosive carbonatite volcano Ol Doinyo Lengai in Tanzania, which erupts natrocarbonatite lavas — illustrate the extreme compositional range possible in rift settings.

Archean magmatism (> 2.5 Ga) differs fundamentally from modern magmatism in several respects. The Archean mantle was hotter than today, producing higher-degree partial melts and erupting komatiites — ultramafic volcanic rocks with MgO > 18 wt% and characteristic spinifex texture (elongate, blade-like olivine or pyroxene crystals) that record rapid quenching of very hot, low-viscosity lavas (eruption temperatures of 1550–1650°C, versus ~1200°C for modern basalts). The origin of the voluminous tonalite-trondhjemite-granodiorite (TTG) suites that dominate Archean cratons is debated: most models involve slab melting or melting of thick mafic lower crust, producing distinctive high-Al, low-Mg, steep REE patterns (high La/Yb, low Yb, indicating garnet as a residual phase) that are the geochemical hallmark of Archean TTGs worldwide.

Chapter 8: Magmatic Ore Deposits

Metals Concentrated by Magmatic Processes

The concentration of economically valuable metals in igneous rock systems — creating ore deposits — occurs through a variety of magmatic and magmatic-hydrothermal processes that selectively enrich certain elements by factors of hundreds to thousands above their average crustal abundances. Understanding these processes is central to mineral exploration and to the long-term security of supply of metals that underpin modern technology.

Orthomagmatic ore deposits form by the separation of an immiscible metallic sulphide or oxide liquid from a silicate magma. In mafic and ultramafic intrusions, when sufficient sulphur is available (either from the magma or from assimilation of sulphur-bearing country rocks), an Fe-Ni-Cu sulphide melt segregates from the silicate melt because the two liquids have limited mutual solubility. Siderophile and chalcophile metals (Ni, Cu, Co, and the platinum-group elements — Pt, Pd, Rh, Ir, Os, Ru) are strongly partitioned into the sulphide liquid, concentrating dramatically relative to the silicate magma. The Sudbury Igneous Complex in Ontario — the product of a large meteorite impact that melted the target rocks and created a unique igneous body — hosts one of the world’s largest Ni-Cu-PGE deposits, formed by this sulphide segregation mechanism.

The R-factor in magmatic sulphide ore deposit formation quantifies the ratio of the volume of silicate magma that equilibrated with the sulphide melt to the volume of sulphide melt itself: \[ C_{\text{sulphide}} = C_0 \cdot D \cdot \frac{R + 1}{R + D} \]

where \( C_{\text{sulphide}} \) is the metal concentration in the sulphide liquid, \( C_0 \) is the metal concentration in the silicate magma, \( D \) is the sulphide-silicate partition coefficient, and \( R \) is the mass ratio of silicate to sulphide magma. For highly siderophile elements with very large \( D \) values (e.g., Pt: \( D \approx 10^4 - 10^5 \)), even small R-factors produce extreme enrichments. The R-factor concept explains why small volumes of sulphide melt that equilibrated with large magma batches (high R) produce the richest PGE ores.

Chromite deposits form by the crystal settling or in-situ crystallization of chromite (\( \text{FeCr}_2\text{O}_4 \)) in layered mafic intrusions when magma mixing events trigger chromite saturation. The Bushveld Complex in South Africa — Earth’s largest layered mafic intrusion (at ~65,000 km² surface exposure) — hosts the world’s largest reserves of Cr and PGE in the UG2 chromitite and Merensky Reef layers, which are interpreted as products of repeated magma mixing events in a stratified magma chamber.

Porphyry copper deposits — the world’s primary source of copper and molybdenum, and an important source of gold — form at shallow depths (1–4 km) above porphyritic calc-alkaline intrusions in arc settings. As hydrous arc magmas crystallize in crustal magma chambers, water and other volatiles become increasingly concentrated in the residual melt until a supercritical fluid phase exsolves. This fluid is enriched in Cu, Mo, Au, and S; it migrates upward along fractures, depositing ore minerals as it cools, depressurizes, and mixes with meteoric water. The resulting deposit geometry — a central, nearly barren intrusion surrounded by concentric shells of stockwork veining with characteristic sulphide mineral zonation (chalcopyrite, bornite, molybdenite core grading outward to pyrite ± galena, sphalerite) — is the template for porphyry deposit classification and exploration worldwide.

Chapter 9: Case Studies in Volcanology and Igneous Petrology

Case Study 1: The 1980 Mount St. Helens Eruption

The 18 May 1980 eruption of Mount St. Helens (Washington, USA) remains one of the most thoroughly studied volcanic events in history and illustrates the cascade of processes from volcanic unrest to catastrophic eruption. Two months of escalating seismicity, ground deformation, and phreatic (steam-driven) explosions preceded the climactic eruption, which was triggered by a magnitude 5.1 earthquake that destabilized the bulging north flank. A massive sector collapse — the largest historical landslide in recorded history — removed the volcano’s summit and northern face, instantaneously decompressing the shallow magmatic system and triggering a lateral blast that devastated 600 km² of forest. The resulting Plinian eruption column reached 24 km into the stratosphere, depositing ash across 11 U.S. states.

The St. Helens eruption demonstrated several key principles: that dome-building eruptions can escalate to cataclysmic events with little warning; that flank instability (in this case amplified by intrusion of a cryptodome of dacitic magma into the north flank) can precede and trigger eruptions; that real-time geophysical monitoring (seismicity, deformation, gas emissions) is essential for volcanic hazard assessment; and that phreatic (non-magmatic) and magmatic eruption phases can occur in sequence during a single eruptive episode.

Case Study 2: Hawaiian Volcanism and Mantle Plumes

The Hawaiian Islands provide an unparalleled natural laboratory for studying the temporal evolution of plume-derived magmas and the interaction between a moving plate and a stationary deep-Earth heat source. Kilauea volcano on the Big Island is one of the world’s most active volcanoes, erupting almost continuously from 1983 to 2018 from its East Rift Zone (the Pu’u ‘Ō’ō–Kupaianaha eruption), producing over 4 km³ of basalt and extending the island coastline by ~230 hectares before transitioning to the Halema’uma’u summit eruption phase.

The geochemical evolution of Hawaiian volcanoes follows a systematic temporal pattern: the pre-shield stage produces small volumes of alkali basalt and basanite; the shield stage (the main volume-building phase) produces tholeiitic basalt at high eruption rates; the post-shield stage generates capping lava flows of alkali basalt, hawaiite, mugearite, and benmoreite as eruption rates decline and the volcano drifts off-axis of the plume; and the rejuvenated stage, occurring long after the main volcanic activity, produces small volumes of highly alkalic lavas (basanite, nephelinite) representing very small-degree melts of the lithospheric mantle.

Case Study 3: Large Igneous Provinces and Mass Extinctions

Large igneous provinces (LIPs) — vast accumulations of mantle-derived mafic magma erupted over geologically brief intervals (typically < 1–5 million years) — represent some of the most consequential volcanic events in Earth history. The Siberian Traps LIP, erupted at the Permian-Triassic boundary (~252 Ma), is temporally associated with the most devastating mass extinction in Earth history (~96% of marine species lost). The Deccan Traps (India, ~66 Ma) erupted close in time to the Chicxulub impact and the end-Cretaceous extinction, contributing to a debate about the relative roles of volcanism and impact in driving the extinction.

The causal mechanisms linking LIP volcanism to mass extinction involve the massive emission of CO₂ (causing greenhouse warming and ocean acidification), SO₂ (causing short-term cooling and acid rain), and halogens (HCl, HBr — causing ozone depletion) over millions of years. The Siberian Traps also intruded into evaporite sequences rich in organic carbon, sulphate, and halite, potentially releasing thermogenic greenhouse gases and halogens at rates far exceeding those estimated from magma volumes alone. LIP studies thus link deep Earth processes — mantle plume dynamics, lithospheric structure, volatile budgets — to global environmental change, demonstrating that volcanology is inseparable from the history of life on Earth.

Back to top