Geology · Fluid Dynamics · Earth Science
📅 April 2026 ⏱ ≈ 12 min read 🎯 Beginner–Intermediate

How Volcanoes Work — Magma, Pressure & Eruption Physics

A volcano is a pressure relief valve for the planet's interior. Molten rock forms deep in the mantle, rises through buoyancy, pools in magma chambers, and — when pressure exceeds the strength of surrounding rock — erupts explosively or effusively. This article traces that journey from partial melting at depth to the pyroclastic flows and lava fields at the surface.

1. Magma Formation in the Mantle

The mantle is mostly solid rock — peridotite — kept solid by pressure despite temperatures of 1 300–1 500 °C. Melting occurs when this stability is disrupted by one of three mechanisms:

Partial melt fractions are typically 1–30%. The silica content of the resulting magma determines its viscosity and eruption style. Basaltic magmas (low SiO₂, ~50%) are hot and runny; rhyolitic magmas (high SiO₂, ~75%) are viscous and gas-rich, prone to explosive eruptions.

η_basalt ≈ 10²–10³ Pa·s (honey-like at eruption temp ~1200 °C) η_rhyolite ≈ 10⁶–10¹² Pa·s (near-solid at surface conditions) Viscosity doubles roughly every 10% SiO₂ increase, and rises exponentially as temperature drops.

2. Buoyancy-Driven Ascent

Magma (density ≈ 2 500–2 700 kg/m³) is less dense than surrounding peridotite (≈ 3 300 kg/m³) and rises through buoyancy. The driving pressure gradient is:

ΔP/Δz = (ρ_rock − ρ_melt) · g ≈ (3300 − 2600) · 9.8 ≈ 6 860 Pa/m For a 10 km column: ΔP ≈ 68 MPa — enough to fracture crustal rock.

The melt collects into pockets, then migrates through permeable rock at low melt fraction, and eventually into larger conduits or dikes — vertical sheet-like fractures through which magma rapidly ascends. Dike propagation is driven by the internal magma pressure exceeding the fracture toughness of host rock.

As magma rises it decompresses, allowing dissolved gases (H₂O, CO₂, SO₂) to exsolve and form bubbles — a process called vesiculation. This is the critical step distinguishing explosive from effusive eruptions.

Volatile solubility: at 200 MPa (≈8 km depth), basaltic melt dissolves ~4 wt% H₂O. As pressure drops during ascent, water solubility decreases dramatically: at 10 MPa it holds only ~0.5 wt%. Exsolution accelerates exponentially as magma nears the surface.

3. Magma Chambers and Pressure Buildup

In many volcanic systems, ascending magma pools in a magma chamber at mid-crustal depths (5–20 km). The chamber is a partially molten zone — not a simple liquid "lake," but a crystal mush with 40–70% solid crystals and a silicate melt interstitial phase.

Pressure in the chamber increases when:

Eruption is triggered when the excess pressure ΔP exceeds the tensile strength of the overlying rock (~10–40 MPa for typical crust). Seismic tomography and ground deformation measurements (InSAR) allow volcanologists to monitor chamber inflation in real time.

ΔP_critical ≈ T_f + ρ_melt · g · h_chamber where T_f = fracture toughness of roof rock (~10–30 MPa) Ground uplift δ ≈ (ΔP · V_chamber) / (4π μ_crust d²) (Mogi point-source model for surface deformation)

4. Eruption Styles

The key variable controlling eruption style is magma viscosity combined with volatile content. High-viscosity, gas-rich magma cannot outgas gently — pressure builds until the melt fragments explosively.

Hawaiian Eruptions

Low-viscosity basaltic lava effuses quietly as lava fountains and flows. Gas escapes gradually. Flow velocities can reach 10 km/h on steep slopes. Classic examples: Kīlauea (Hawaii), Stromboli (Italy).

Strombolian Eruptions

Intermittent mild explosions every few minutes to hours. Large gas pockets rise through the conduit and burst at the surface, throwing incandescent lava bombs up to 300 m. Characteristic of basaltic to basaltic-andesitic magmas.

Vulcanian Eruptions

Short violent explosions when a solidified plug sealing the vent is blown out by overpressured gas. Produces thick ash clouds and ballistic blocks up to several km range.

Plinian Eruptions

The most violent style — named after Pliny the Younger's account of the 79 CE Vesuvius eruption. A sustained jet of gas, ash and pumice forms an eruption column that can reach 40–50 km into the stratosphere. Mass discharge rate can exceed 10⁸ kg/s.

Column height H ≈ k · Q^(1/4) (Morton-Taylor-Turner plume equation) where Q = mass flux (kg/s), k ≈ 0.24 km/(kg/s)^(1/4) Q = 10⁸ kg/s → H ≈ 0.24 × (10⁸)^(0.25) ≈ 24 km

Caldera-Forming "Super-eruptions"

When a very large magma chamber partially empties, the roof can collapse inward forming a caldera depression (Yellowstone, Toba, Campi Flegrei). These events release >1 000 km³ of material (VEI ≥ 8) and can trigger volcanic winters.

5. Lava Flow Rheology

Lava behaves as a Bingham plastic with a yield stress τ_y — it only flows when shear stress exceeds the yield threshold (due to crystals and vesicles forming a network). Above τ_y it follows a Newtonian-like relationship:

τ = τ_y + η_plastic · (du/dy) (Bingham constitutive equation) For a channel flow of width W and depth h on slope θ: U_max = (ρg sin θ / 2η) · (h − h_y)² where h_y = τ_y / (ρg sin θ) is the rigid plug thickness at surface.

As lava cools it develops a solid crust even while the interior remains mobile. This creates:

Simulation note: SPH (Smoothed Particle Hydrodynamics) and LBM (Lattice Boltzmann Method) are the primary numerical methods used to simulate lava flows, including the non-Newtonian rheology and thermal coupling. The tectonic plates simulation uses a simplified viscous model.

6. Pyroclastic Flows and Tephra

When a Plinian eruption column collapses, or when lava domes shatter, the mixture of hot gas and fragmented rock (pyroclasts) becomes denser than air and flows as a gravity current — a pyroclastic density current (PDC).

Tephra (airborne volcanic fragments) is classified by grain size: volcanic blocks (>64 mm), lapilli (2–64 mm), volcanic ash (<2 mm). Fine ash can remain suspended in the stratosphere for 1–3 years, scattering sunlight and reducing global temperatures by up to 0.5°C (Pinatubo 1991 cooled Earth by ~0.4°C for 2 years).

Stokes settling velocity for a tephra particle: v_s = (2r² · (ρ_particle − ρ_air) · g) / (9 η_air) For r = 25 μm ash: v_s ≈ 0.03 m/s → stays airborne for days For r = 1 mm lapilli: v_s ≈ 50 m/s → settles within minutes

7. Connection to Simulation

The Tectonic Plates simulation models the large-scale mantle convection that drives plate motion and creates the conditions for volcanism. Subduction zones — where one plate dives under another — are the source of 75% of Earth's volcanic eruptions.

The simulation uses a simplified viscous fluid model on a 2-D spherical surface, with temperature-dependent viscosity. As plates converge, cooler denser material sinks, driving return flow that can bring hot plumes to the surface — a visual analogue of the flux melting mechanism described above.

🌋 Open Tectonic Plates Simulation →