5.4 Hydrothermal Vents

Hydrothermal vents are openings in the seafloor where geothermally heated water discharges into the ocean. Discovered in 1977 on the Galapagos Rift, they sustain extraordinary chemosynthetic ecosystems and profoundly influence ocean chemistry, transferring heat and dissolved metals from Earth's interior to the deep ocean.

Black Smokers and White Smokers

Black Smokers

High-temperature vents (up to 350–407 °C) that emit particle-rich fluid. The "smoke" consists of metal sulfide precipitates (FeS, CuFeS₂, ZnS) forming instantly when hot, reduced, metal-rich vent fluid meets cold, oxygenated seawater. Chimney structures grow up to 60 m tall. The fluid is acidic (pH 2–3) and rich in H₂S, Fe²⁺, Mn²⁺, Cu²⁺, Zn²⁺, and dissolved gases (H₂, CH₄, CO₂).

White Smokers

Lower-temperature vents (150–300 °C) emitting milky-white fluids containing barium sulfate (barite), anhydrite (CaSO₄), and silica precipitates. They occur on the flanks of vent fields where fluid has mixed with ambient seawater in the subsurface. The Lost City hydrothermal field (off-axis, on 1.5 Myr old crust) produces alkaline (pH 9–11), low-temperature (40–90 °C) fluids rich in H₂ and CH₄ via serpentinisation of mantle peridotite.

Vent Fluid Chemistry

Seawater percolates through fractured oceanic crust, is heated by underlying magma, reacts with basalt, and exits at the seafloor chemically transformed. Key processes:

Sulfate Removal

$\text{CaSO}_4 \rightarrow \text{Ca}^{2+} + \text{SO}_4^{2-}$ — anhydrite precipitates above 150 °C, stripping sulfate from the fluid.

Serpentinisation

$2\text{Fe}_2\text{SiO}_4 + 3\text{H}_2\text{O} \rightarrow 2\text{Fe}_2\text{O}_3 + 2\text{SiO}_2 + 3\text{H}_2$ — releases hydrogen and leaches metals from basalt.

Metal Sulfide Precipitation

$\text{Fe}^{2+} + \text{H}_2\text{S} \rightarrow \text{FeS} + 2\text{H}^+$ — forms the black "smoke" and massive sulfide ore deposits.

Heat Flux from Mid-Ocean Ridges

The total hydrothermal heat flux is estimated at $\sim 10^{13}$ W, representing ~25–30% of the total heat loss from oceanic crust. The convective heat carried by vent fluids is:

$$Q = \dot{m}\,c_p\,(T_{vent} - T_{ambient})$$

where $\dot{m}$ is mass flow rate, $c_p \approx 4200\;\text{J/(kg\cdot K)}$ is specific heat capacity, and the temperature difference can exceed 350 °C.

Hydrothermal Plumes & MTT Model

Hot vent fluid is buoyant and rises as a turbulent plume until reaching neutral buoyancy (typically 100–400 m above the seafloor). The Morton–Taylor–Turner (MTT) entrainment model describes the dynamics:

MTT Plume Equations

$$\frac{d}{dz}(b^2 w) = 2\alpha\, b\, w$$

$$\frac{d}{dz}(b^2 w^2) = b^2 g'$$

$$\frac{d}{dz}(b^2 w g') = -b^2 w N^2$$

where $b$ is plume radius, $w$ is vertical velocity,$g' = g(\rho_a - \rho)/\rho_a$ is reduced gravity, $\alpha \approx 0.1$ is the entrainment coefficient, and $N$ is the buoyancy frequency. Rise height scales as:

$$z_{max} = 3.76\,B_0^{1/4}\,N^{-3/4}$$

Vent Biomes & Chemosynthesis

$$\text{H}_2\text{S} + \text{O}_2 + \text{CO}_2 + \text{H}_2\text{O} \rightarrow \text{CH}_2\text{O} + \text{H}_2\text{SO}_4$$

Sulfide oxidation provides energy to fix CO₂ into organic matter.

Giant Tube Worms (Riftia)

Up to 2 m long. No mouth, gut, or anus. Host endosymbiotic sulfide-oxidising bacteria in a specialised organ called the trophosome.

Vent Shrimp (Rimicaris)

Swarm around vent orifices. Possess dorsal "eyes" that detect the dim infrared glow from superheated fluid. Epibiotic bacteria on their gills.

Vent Mussels & Clams

Harbour chemosynthetic bacteria in gill tissues. Can dominate at lower-temperature diffuse venting sites around the vent periphery.

Pompeii Worms (Alvinella)

Inhabit chimney walls at up to 80 °C—among the most thermotolerant animals. Symbiotic bacteria coat their dorsal surface.

Impact on Ocean Chemistry

Hydrothermal circulation is a major source of Fe, Mn, H₂, CH₄, and ³He to the ocean, and a sink for Mg and SO₄. The entire ocean volume cycles through mid-ocean ridges every ~8 Myr.

Fe flux

~10⁹ mol/yr

Mg removal

~1.5 × 10¹² mol/yr

³He flux

~10³ mol/yr

Python: Plume Rise, Mixing Model & Metal Speciation

Python

!/usr/bin/env python3

script.py61 lines

#!/usr/bin/env python3
"""hydrothermal_vents.py - Plume rise, mixing model, metal removal."""
import numpy as np
import matplotlib.pyplot as plt

# ── 1. Plume Rise Height ──
def plume_rise_height(B0, N):
    """Max rise height (Morton-Taylor-Turner). B0: buoyancy flux (m4/s3), N: buoyancy freq (1/s)."""
    return 3.76 * B0**0.25 * N**(-0.75)

B0_range = np.logspace(-3, 0, 100)
N_values = [5e-4, 1e-3, 2e-3]

fig, axes = plt.subplots(1, 3, figsize=(18, 5))
ax = axes[0]
for N in N_values:
    z_max = plume_rise_height(B0_range, N)
    ax.loglog(B0_range, z_max, lw=2, label=f'N = {N:.0e} 1/s')
ax.set_xlabel('Buoyancy Flux B0 (m4/s3)')
ax.set_ylabel('Rise Height (m)')
ax.set_title('Hydrothermal Plume Rise Height')
ax.legend(); ax.grid(True, alpha=0.3, which='both')

# ── 2. Conservative Mixing Model ──
dilution = np.logspace(0, 5, 200)
T_vent, T_amb = 350.0, 2.0
Fe_vent, H2S_vent, Mn_vent = 1600.0, 8000.0, 600.0

T_mix = T_amb + (T_vent - T_amb) / dilution
Fe_mix = Fe_vent / dilution
H2S_mix = H2S_vent / dilution

ax = axes[1]
ax.semilogx(dilution, T_mix, 'r-', lw=2, label='Temperature (C)')
ax.set_xlabel('Dilution Factor')
ax.set_ylabel('Temperature (C)', color='r')
ax2 = ax.twinx()
ax2.semilogx(dilution, Fe_mix, 'b--', lw=2, label='Fe (umol/kg)')
ax2.semilogx(dilution, H2S_mix, 'g--', lw=2, label='H2S (umol/kg)')
ax2.set_ylabel('Concentration (umol/kg)', color='b')
ax.set_title('Conservative Mixing of Vent Fluid')
lines1, labels1 = ax.get_legend_handles_labels()
lines2, labels2 = ax2.get_legend_handles_labels()
ax.legend(lines1 + lines2, labels1 + labels2, fontsize=7)

# ── 3. Metal Removal Kinetics ──
t_hours = np.linspace(0, 48, 200)
Fe_dissolved = Fe_vent / 100 * np.exp(-0.5 * t_hours)
Mn_dissolved = Mn_vent / 100 * np.exp(-0.02 * t_hours)

ax = axes[2]
ax.plot(t_hours, Fe_dissolved, 'b-', lw=2, label='Dissolved Fe')
ax.plot(t_hours, Mn_dissolved, 'orange', lw=2, label='Dissolved Mn')
ax.set_xlabel('Time After Discharge (hours)')
ax.set_ylabel('Concentration (umol/kg)')
ax.set_title('Metal Removal in Buoyant Plume')
ax.legend(); ax.grid(True, alpha=0.3)

plt.tight_layout()
plt.savefig('hydrothermal_vents.png', dpi=150)
plt.show()

Click Run to execute the Python code

Code will be executed with Python 3 on the server

Fortran: Morton–Taylor–Turner Buoyant Plume Model

This Fortran program integrates the MTT entrainment equations to compute the evolution of plume radius, velocity, and buoyancy from the vent orifice to neutral buoyancy height.

Fortran: Morton–Taylor–Turner Buoyant Plume Model

Fortran

============================================================

program.f9068 lines

program mtt_plume
  ! ============================================================
  ! Morton-Taylor-Turner buoyant plume model
  ! d(b^2 w)/dz   = 2 alpha b w        (mass/entrainment)
  ! d(b^2 w^2)/dz = b^2 g'             (momentum)
  ! d(b^2 w g')/dz = -b^2 w N^2        (buoyancy)
  ! ============================================================
  implicit none
  integer, parameter :: dp = selected_real_kind(15, 307)
  integer, parameter :: nz = 10000
  real(dp) :: alpha_e, N2
  real(dp) :: b0, w0, gp0
  real(dp) :: dz, z_max
  real(dp) :: b, w, gp
  real(dp) :: Q, M, F, Q_new, M_new, F_new
  real(dp) :: z_rise
  integer  :: i

alpha_e = 0.1_dp
  N2      = 1.0d-6              ! N ~ 1e-3 s^-1
  z_max   = 500.0_dp
  dz      = z_max / real(nz, dp)

! Initial conditions at vent orifice
  b0  = 0.1_dp                  ! radius (m)
  w0  = 2.0_dp                  ! vertical velocity (m/s)
  gp0 = 0.03_dp                 ! reduced gravity (m/s^2)

! Top-hat fluxes
  Q = b0**2 * w0
  M = b0**2 * w0**2
  F = b0**2 * w0 * gp0

open(unit=10, file='plume_profile.dat', status='replace')
  write(10,'(A)') '# z(m)  b(m)  w(m/s)  gp(m/s2)  T_anom(C)'

z_rise = 0.0_dp

do i = 1, nz
    if (Q <= 0.0_dp .or. M <= 0.0_dp) exit
    b  = sqrt(Q**2 / M)
    w  = M / Q
    gp = F / Q

write(10,'(5(F12.4,1X))') real(i,dp)*dz, b, w, gp, gp / 1.962d-3

if (gp <= 0.0_dp) then
      z_rise = real(i, dp) * dz
      exit
    end if

! Euler forward integration
    Q_new = Q + dz * 2.0_dp * alpha_e * sqrt(M)
    M_new = M + dz * Q * gp / w
    F_new = F + dz * (-Q * N2)

Q = Q_new; M = M_new; F = F_new
  end do

close(10)
  if (z_rise > 0.0_dp) then
    write(*,'(A,F8.1,A)') 'Neutral buoyancy height: ', z_rise, ' m'
  else
    write(*,'(A)') 'Plume did not reach neutral buoyancy in domain.'
  end if
  write(*,'(A)') 'Output: plume_profile.dat'

end program mtt_plume

Click Run to execute the Fortran code

Code will be compiled with gfortran and executed on the server

Off-Axis Hydrothermal Activity & Metal Oxide Precipitation

Hydrothermal activity is not limited to the ridge axis. Off-axis systems such as the Lost City hydrothermal field (discovered in 2000, 30°N Mid-Atlantic Ridge) occur on older crust (1–2 Myr) where seawater reacts with ultramafic mantle rock (peridotite) through serpentinisation. These systems produce highly alkaline (pH 9–11), hydrogen- and methane-rich fluids at moderate temperatures (40–90 °C), with towering carbonate chimneys up to 60 m tall. Lost City-type vents are thought to be far more common than previously recognised and have been proposed as potential settings for the origin of life.

Metal Oxide Precipitation

As dissolved Fe²⁺ and Mn²⁺ in hydrothermal plumes mix with oxygenated deep water, they oxidise and precipitate as Fe-oxyhydroxide (FeOOH) and Mn-oxide (MnO₂) particles. Iron oxidises rapidly (hours to days), while manganese oxidation is much slower (weeks to months), allowing dissolved Mn to be transported thousands of kilometres from the vent site. These metalliferous sediments accumulate at rates of 1–10 mm/kyr near ridges and are enriched in trace metals (Cu, Zn, Pb, Co, REE).

Axial Magma Chambers

The heat source driving hydrothermal circulation is the axial magma chamber (AMC) beneath the ridge crest. At fast-spreading ridges (EPR), the AMC is a thin (tens of metres) melt lens at 1–2 km below the seafloor, overlain by a broader region of partially molten crystal mush. At slow-spreading ridges (MAR), discrete melt bodies are ephemeral and spatially limited. The AMC depth and geometry control the temperature, chemistry, and longevity of vent fields.

Seafloor Massive Sulfide (SMS) Deposits

The accumulation of metal sulfide precipitates around hydrothermal vents creates seafloor massive sulfide (SMS) deposits that are modern analogues of ancient volcanogenic massive sulfide (VMS) ore deposits found on land. These deposits are of increasing commercial interest for deep-sea mining, though significant environmental concerns exist.

Mineral Composition

SMS deposits contain pyrite (FeS₂), chalcopyrite (CuFeS₂), sphalerite (ZnS), galena (PbS), and precious metals (Au, Ag). Copper and zinc grades can reach 5–15%, far exceeding land-based ores. The TAG mound on the MAR contains an estimated 3.9 million tonnes of sulfide ore.

Growth & Lifetime

Individual chimneys grow at rates of 1–30 cm/day but are fragile and topple frequently. Vent fields remain active for $10^3$–$10^5$ years. Radiometric dating (⁶⁶Ra/Ba) of chimney sulfides constrains individual chimney ages to years–decades. Inactive sulfide mounds persist on the seafloor for millions of years before being subducted or obducted onto continents.

Zone Refining in Chimneys

As hot fluid flows through a chimney, minerals precipitate and re-dissolve along the thermal gradient, creating a characteristic zonation: an inner lining of chalcopyrite (high-T, >300 °C), surrounded by pyrite and sphalerite (intermediate-T), with an outer layer of anhydrite and amorphous silica (low-T). This zone refining concentrates copper in the chimney interior.

Key Equations Summary

Buoyancy Flux

$$B_0 = g\frac{\Delta\rho}{\rho_0}\,Q_0 = g\,\alpha_T\,\Delta T\,Q_0$$

where $Q_0$ is the volume flux and $\alpha_T$ is the thermal expansion coefficient.

Residence Time of Vent Chemicals

$$\tau_{res} = \frac{M_{ocean}}{F_{input}}$$

Hydrothermal iron has a short residence time (~100 years) due to rapid scavenging, while manganese persists for ~1000 years. Helium-3 ($^3$He) is conservative and traces plume dispersal across entire ocean basins.

Darcy Flow in Hydrothermal Cells

$$\mathbf{q} = -\frac{k}{\mu}(\nabla p - \rho\mathbf{g})$$

Crustal permeability $k \sim 10^{-14}$ to $10^{-10}\;\text{m}^2$ controls the vigour of convective circulation.

← Marine Sediments Coastal Geology →