Module 0

Ocean Zones & Physical Chemistry

From sunlit surface to the hadal deep: the vertical structure of the ocean and its driving physics

0.1 Depth Zonation of the Ocean

The ocean is divided into five major depth zones, each with distinct physical, chemical, and biological characteristics. These zones are defined primarily by light penetration, temperature, pressure, and nutrient availability.

Epipelagic (Sunlit Zone)

0–200 m. Receives sufficient sunlight for photosynthesis. Temperature 15–30 Β°C at surface, declining through the seasonal thermocline. Contains >90% of all marine life biomass. Wind-driven mixing maintains a relatively uniform mixed layer (50–200 m). O\(_2\) near saturation due to air-sea exchange and photosynthesis.

Mesopelagic (Twilight Zone)

200–1000 m. Light insufficient for photosynthesis but detectable by organisms. Temperature drops from ~10 Β°C to ~4 Β°C through the permanent thermocline. Contains the oxygen minimum zone (OMZ) where bacterial respiration exceeds O\(_2\) supply. Home to the world's largest daily migration (see Module 2).

Bathypelagic (Midnight Zone)

1000–4000 m. Completely dark except for bioluminescence. Temperature 1–4 Β°C, nearly uniform. Pressure 100–400 atm. Sparse life depends entirely on sinking organic matter (β€œmarine snow”) and vertical migrants. O\(_2\) levels recover due to thermohaline circulation bringing oxygenated polar water.

Abyssopelagic (Abyssal Zone)

4000–6000 m. Covers ~75% of the ocean floor. Temperature 1–2 Β°C, near freezing. Pressure 400–600 atm. Nutrient concentrations high (decomposition exceeds biological uptake). Extremely low organic carbon flux: <1% of surface production reaches this depth.

Hadopelagic (Hadal Zone)

6000–11,000 m. Found only in oceanic trenches (Mariana, Tonga, Kermadec, Philippine, etc.). Pressure 600–1,100 atm. Temperature 1–4 Β°C (some trenches show slight warming from adiabatic compression). Despite extreme conditions, hosts amphipods, polychaetes, foraminifera, and xenophyophores. The Mariana Trench's Challenger Deep reaches 10,994 m.

Zone Boundaries as Ecological Transitions

Zone boundaries are not sharp lines but gradients. The euphotic zone boundary (base of epipelagic) is defined by the 1% light level β€” the depth where photosynthetically active radiation (PAR) falls to 1% of its surface value. This varies with water clarity: ~20 m in turbid coastal waters to ~200 m in the clearest subtropical gyres. The compensation depth, where gross photosynthesis exactly balances algal respiration, typically lies at 1–5% of surface PAR.

Understanding these zones is foundational for ocean biogeochemistry and connects directly toclimate-biodiversity interactions andclimate science, since the ocean's vertical structure controls the biological carbon pump that sequesters atmospheric CO\(_2\).

0.2 Light Attenuation: Beer-Lambert Law

Light intensity decreases exponentially with depth due to absorption and scattering by water molecules, dissolved organic matter (CDOM), phytoplankton pigments, and suspended particles. The governing law is the Beer-Lambert law of exponential attenuation:

\( I(z) = I_0 \cdot e^{-K_d \cdot z} \)

Downwelling irradiance at depth z (m)

where \(I_0\) is the surface irradiance (typically 400–2000 \(\mu\)mol photons/m\(^2\)/s for PAR), and \(K_d\) is the diffuse attenuation coefficient (m\(^{-1}\)).

Derivation: Wavelength Dependence of \(K_d\)

The attenuation coefficient is the sum of absorption and scattering contributions:

\( K_d(\lambda) = a_w(\lambda) + a_{ph}(\lambda) + a_{CDOM}(\lambda) + b_b(\lambda) \)

Pure water absorption + phytoplankton + CDOM + backscattering

Pure water absorption \(a_w(\lambda)\) dominates the wavelength structure. Red light (\(\lambda \approx 700\) nm) is strongly absorbed: \(a_w(700) \approx 0.65\;\text{m}^{-1}\), meaning 50% absorption in just 1.1 m. Blue light (\(\lambda \approx 475\) nm) penetrates deepest: \(a_w(475) \approx 0.015\;\text{m}^{-1}\), reaching 1% at ~300 m in clear ocean water. This is why the deep ocean appears blue.

Typical \(K_d\) values for PAR (broadband):

  • Clear open ocean (oligotrophic): \(K_d \approx 0.02\text{--}0.04\;\text{m}^{-1}\)
  • Productive open ocean (mesotrophic): \(K_d \approx 0.06\text{--}0.10\;\text{m}^{-1}\)
  • Coastal waters (eutrophic): \(K_d \approx 0.10\text{--}0.50\;\text{m}^{-1}\)
  • Turbid estuaries: \(K_d \approx 1.0\text{--}5.0\;\text{m}^{-1}\)

Euphotic Zone and Compensation Depth

The euphotic zone depth \(z_{eu}\) is where irradiance falls to 1% of the surface value. Setting \(I(z_{eu}) = 0.01 \cdot I_0\):

\( 0.01 = e^{-K_d \cdot z_{eu}} \)

\( z_{eu} = \frac{\ln(100)}{K_d} = \frac{4.605}{K_d} \)

For clear ocean water (\(K_d = 0.03\)): \(z_{eu} \approx 154\) m. For coastal water (\(K_d = 0.20\)): \(z_{eu} \approx 23\) m.

The compensation depth \(z_c\) is where gross photosynthesis\(P_g\) equals phytoplankton respiration \(R\). Below \(z_c\), cells consume more carbon than they fix. For a P-I curve \(P_g = P_{\max} \tanh(\alpha I / P_{\max})\), the compensation irradiance \(I_c\) satisfies \(P_{\max} \tanh(\alpha I_c / P_{\max}) = R\), typically at 1–5 \(\mu\)mol photons/m\(^2\)/s.

0.3 Hydrostatic Pressure

Pressure increases linearly with depth in the ocean, governed by hydrostatic equilibrium. The water column exerts a weight proportional to its height and density:

\( P(z) = P_0 + \rho g z \)

\(P_0 = 1\;\text{atm} = 101{,}325\;\text{Pa}\), \(\rho \approx 1025\;\text{kg/m}^3\), \(g = 9.81\;\text{m/s}^2\)

Derivation: Pressure at Key Depths

One atmosphere of additional pressure is added for every 10.33 m of seawater:

\( \Delta z_{1\,\text{atm}} = \frac{P_0}{\rho g} = \frac{101{,}325}{1025 \times 9.81} \approx 10.07\;\text{m} \)

(Rounded to ~10 m for the β€œ1 atm per 10 m” rule of thumb)

At the bottom of the Mariana Trench (Challenger Deep, \(z = 10{,}994\) m):

\( P = 1 + \frac{1025 \times 9.81 \times 10{,}994}{101{,}325} \approx 1 + 1{,}091 = 1{,}092\;\text{atm} \)

\( \approx 1{,}100\;\text{atm} \approx 110\;\text{MPa} \)

Biochemical Effects of Pressure

High pressure profoundly affects macromolecular structure. By Le Chatelier's principle, pressure favours states with smaller molar volume. The volume change of a reaction determines its pressure sensitivity:

\( \left(\frac{\partial \ln K}{\partial P}\right)_T = -\frac{\Delta V}{RT} \)

Pressure dependence of equilibrium constant

Key effects on biomolecules:

  • Protein unfolding: The hydrophobic core of proteins has a larger molar volume than the unfolded state (water fills cavities). At >200 atm, many mesophilic enzymes lose activity. \(\Delta V_{\text{unfold}} \approx -50\;\text{to}\;-100\;\text{mL/mol}\)
  • Membrane compression: Lipid bilayers become more ordered and rigid. Deep-sea organisms incorporate more unsaturated fatty acids to maintain fluidity (homeoviscous adaptation).
  • TMAO accumulation: Trimethylamine N-oxide increases linearly with depth in fish (approximately 0.04 mol/atm). It stabilises proteins by strengthening the hydration shell around hydrophobic groups.

0.4 Temperature Structure: Mixed Layer, Thermocline, Deep Ocean

The ocean's vertical temperature structure consists of three layers:

  • Mixed layer (0–100 m): Warm, well-mixed by wind, waves, and convection. Temperature ~20–28 Β°C in tropics, ~5–15 Β°C at mid-latitudes. Nearly isothermal due to turbulent mixing.
  • Thermocline (100–1000 m): Rapid temperature decrease, acting as a barrier to vertical mixing. The permanent thermocline exists year-round; a seasonal thermocline forms and erodes above it.
  • Deep ocean (below ~1000 m): Cold (1–4 Β°C) and nearly uniform. Filled by thermohaline circulation β€” cold, dense water sinking at high latitudes (North Atlantic Deep Water, Antarctic Bottom Water).

Heat Diffusion in the Thermocline

The thermocline can be modelled as a balance between downward diffusion of heat from the warm surface and upwelling of cold deep water. At steady state, the advection-diffusion equation gives:

\( w \frac{\partial T}{\partial z} = \kappa \frac{\partial^2 T}{\partial z^2} \)

Munk (1966) abyssal recipe: w = upwelling velocity, \(\kappa\) = turbulent diffusivity

The solution is an exponential profile:

\( T(z) = T_{\text{deep}} + (T_{\text{surface}} - T_{\text{deep}}) \cdot \exp\!\left(-\frac{w}{\kappa}\,z\right) \)

Scale height \(h = \kappa/w \approx\) 500–1000 m defines thermocline thickness

With \(\kappa \approx 10^{-4}\;\text{m}^2/\text{s}\) (diapycnal diffusivity) and\(w \approx 10^{-7}\;\text{m/s}\) (\(\sim 3\) m/yr upwelling), the thermocline scale height is \(h = 10^{-4}/10^{-7} = 1000\) m, consistent with observations. This connects directly to ocean circulation discussed in Climate Science Module 2.

0.5 Dissolved Oxygen Profile & the Oxygen Minimum Zone

Dissolved oxygen shows a distinctive vertical profile with three regimes:

  • Surface (0–100 m): Near saturation (~6–8 mL/L) due to air-sea gas exchange and photosynthetic production. Saturation follows Henry's law, decreasing with temperature (warmer water holds less O\(_2\)).
  • OMZ (200–1000 m): O\(_2\) drops to <0.5 mL/L in the most intense OMZs (eastern tropical Pacific, Arabian Sea, Bay of Bengal). Biological oxygen demand from bacterial decomposition of sinking organic matter exceeds the rate of O\(_2\) supply by lateral advection and mixing.
  • Deep ocean (>1000 m): O\(_2\) gradually increases to 4–6 mL/L, replenished by thermohaline circulation carrying cold, O\(_2\)-saturated water from the polar surface to the deep ocean interior.

O\(_2\) Consumption Rate Derivation

At steady state, the vertical O\(_2\) distribution balances advection, diffusion, and biological consumption:

\( w \frac{\partial [\text{O}_2]}{\partial z} + \kappa \frac{\partial^2 [\text{O}_2]}{\partial z^2} = R_{\text{bio}}(z) \)

Biological consumption rate \(R_{\text{bio}}\) (\(\mu\)mol/kg/yr)

The bacterial decomposition follows sinking organic matter flux, which decays as a power law (Martin curve, see Module 2):

\( R_{\text{bio}}(z) = R_0 \left(\frac{z}{z_0}\right)^{-b} \quad \text{with } b \approx 0.86 \)

Where \(R_0 \approx 5\text{--}20\;\mu\text{mol O}_2/\text{kg/yr}\) at the base of the euphotic zone (\(z_0 = 100\) m). The stoichiometry of organic matter remineralisation follows the Redfield ratio:\( (\text{CH}_2\text{O})_{106}(\text{NH}_3)_{16}(\text{H}_3\text{PO}_4) + 138\text{O}_2 \to 106\text{CO}_2 + 16\text{HNO}_3 + \text{H}_3\text{PO}_4 + 122\text{H}_2\text{O} \). Thus 138 moles of O\(_2\) are consumed per 106 moles of organic carbon remineralised.

0.6 Salinity, pH, and Nutrient Profiles

Pycnocline and Density Stratification

The pycnocline is the zone of rapid density increase with depth, controlled by both temperature (thermocline) and salinity (halocline). Seawater density is computed from the equation of state:

\( \rho = \rho(T, S, P) \approx \rho_0 \left[1 - \alpha(T - T_0) + \beta(S - S_0)\right] \)

Linearised: \(\alpha \approx 2 \times 10^{-4}\;\text{K}^{-1}\) (thermal expansion), \(\beta \approx 7.5 \times 10^{-4}\;\text{psu}^{-1}\) (haline contraction)

Ocean pH and the Carbonate System

Surface ocean pH is ~8.1, decreasing with depth to ~7.6–7.8 due to CO\(_2\) released by remineralisation. The carbonate chemistry equilibria:

\( \text{CO}_2 + \text{H}_2\text{O} \rightleftharpoons \text{H}_2\text{CO}_3 \rightleftharpoons \text{H}^+ + \text{HCO}_3^- \rightleftharpoons 2\text{H}^+ + \text{CO}_3^{2-} \)

\(pK_1 \approx 6.1\), \(pK_2 \approx 9.3\) (at 25 Β°C, S = 35)

As CO\(_2\) increases from remineralisation, pH drops and carbonate ion concentration decreases. Below the carbonate compensation depth (CCD, ~4500 m), the water is sufficiently undersaturated that CaCO\(_3\) dissolves faster than it accumulates β€” no carbonate sediment is preserved. This has implications forcarbon cycle feedbacks.

Redfield Ratio: C:N:P = 106:16:1

Alfred Redfield (1934) discovered that the elemental composition of marine plankton and the dissolved nutrients in deep water follow a remarkably consistent ratio:

\( \text{C} : \text{N} : \text{P} = 106 : 16 : 1 \)

Redfield-Ketchum-Richards ratio (extended: C:N:P:Si = 106:16:1:15 for diatoms)

This ratio emerges because plankton assimilate nutrients in fixed proportions during growth, and deep-water nutrient concentrations reflect cumulative remineralisation. The deviation from Redfield reveals nutrient limitation:

  • If \(\text{N:P} < 16\) in surface water: nitrogen-limited (most of the ocean)
  • If \(\text{N:P} > 16\): phosphorus-limited (Mediterranean, some coastal)
  • Iron limitation in HNLC regions (Southern Ocean, subarctic Pacific) overrides N and P

The nutricline β€” the depth of rapid nutrient increase β€” typically coincides with the thermocline (~100–500 m), since nutrients accumulate where remineralisation occurs and vertical mixing is suppressed by stratification.

0.7 Seawater Equation of State & Stability

The full UNESCO equation of state (EOS-80) relates seawater density to temperature, salinity, and pressure through a polynomial with 41 coefficients. For most oceanographic applications, the simplified linear approximation suffices for small perturbations from a reference state:

\( \sigma_t = (\rho - 1000)\;\text{kg/m}^3 \)

Sigma-t notation: typical values 24–28 kg/m\(^3\)

Potential density \(\sigma_\theta\) removes the effect of adiabatic compression by referencing all parcels to surface pressure. This is essential for assessing vertical stability: a water column is stably stratified when \(\partial \sigma_\theta / \partial z > 0\)(density increases with depth).

Brunt-VΓ€isΓ€lΓ€ Frequency

The buoyancy frequency \(N\) quantifies the strength of stratification and the frequency of internal gravity waves:

\( N^2 = -\frac{g}{\rho_0} \frac{\partial \sigma_\theta}{\partial z} \)

Typical values: \(N \approx 10^{-2}\;\text{s}^{-1}\) in thermocline, \(N \approx 10^{-3}\;\text{s}^{-1}\) in deep ocean

When \(N^2 > 0\), the water column is stable and resists vertical mixing. A displaced parcel oscillates about its equilibrium depth with period \(T = 2\pi/N\). In the pycnocline,\(T \approx 10\) minutes; in the deep ocean, \(T \approx 1\) hour.

Strong stratification (\(N^2\) large) in the tropics inhibits vertical nutrient supply to the surface, explaining why tropical gyres are oligotrophic despite high solar radiation. Conversely, weak stratification at high latitudes allows deep winter convection and nutrient replenishment, fueling productive spring blooms (see Module 1).

Thermohaline Circulation and Ventilation

The deep ocean is ventilated by the thermohaline (meridional overturning) circulation. Cold, saline surface water at high latitudes (principally the North Atlantic and around Antarctica) becomes dense enough to sink to the abyss:

  • North Atlantic Deep Water (NADW): forms in the Greenland-Iceland-Norwegian seas, sinks to ~2000–4000 m, flows southward. Temperature ~2–4 Β°C, salinity ~34.9 psu.
  • Antarctic Bottom Water (AABW): forms on Antarctic continental shelf (brine rejection from sea ice formation), sinks to the very bottom (>4000 m). Temperature <0 Β°C, salinity ~34.7 psu. The densest water mass in the ocean.
  • Antarctic Intermediate Water (AAIW): subducts at the Antarctic Convergence (~50Β°S), spreads northward at ~700–1200 m. Low salinity (~34.3 psu), high O\(_2\).

The global overturning circulation takes approximately 1000–1500 years for a full cycle. This sets the timescale for deep-ocean ventilation and is crucial for understanding oceanic carbon sequestration and climate dynamics. Water parcels that sink in the North Atlantic today will not return to the surface for a millennium, carrying with them any dissolved CO\(_2\) and nutrients they contain.

0.8 Sound Propagation & the SOFAR Channel

Sound speed in seawater depends on temperature, salinity, and pressure through an empirical relationship:

\( c \approx 1449 + 4.6T - 0.055T^2 + 0.0003T^3 + (1.39 - 0.012T)(S - 35) + 0.017z \)

Mackenzie (1981): c in m/s, T in Β°C, S in psu, z in m

Temperature dominates near the surface (warmer = faster), while pressure dominates at depth (deeper = faster). This creates a sound speed minimum at ~700–1200 m depth (the SOFAR channel), where sound is trapped by refraction and can propagate for thousands of kilometres with minimal attenuation.

Marine mammals exploit this: blue whale calls at the SOFAR channel axis can be detected across entire ocean basins. The minimum sound speed \(c_{\min} \approx 1480\;\text{m/s}\) occurs where\(\partial c / \partial z = 0\):

\( \frac{\partial c}{\partial z} = \underbrace{\frac{\partial c}{\partial T}\frac{\partial T}{\partial z}}_{< 0\;\text{(cooling)}} + \underbrace{\frac{\partial c}{\partial P}\frac{\partial P}{\partial z}}_{> 0\;\text{(compression)}} = 0 \)

By Snell's law, sound rays bend toward regions of lower sound speed, trapping acoustic energy in the channel. This phenomenon is exploited in ocean acoustic tomography to measure large-scale temperature changes, directly relevant to climate monitoring.

The SOFAR channel also has biological significance: marine mammals use low-frequency sound for long-range communication. Blue whales produce calls at 10–20 Hz that can travel over 1000 km through the SOFAR channel. Anthropogenic noise from shipping, sonar, and seismic surveys is increasingly disrupting these communication channels, with documented effects on whale behaviour, stress hormones, and stranding events. The sound intensity from large container ships reaches 190 dB re 1 \(\mu\)Pa at 1 m, significantly raising the ambient noise floor in the 20–200 Hz frequency band throughout the world's oceans.

Ocean Depth Profile: Multi-Variable Cross-Section

The following diagram shows the vertical distribution of temperature, dissolved oxygen, light intensity, pH, and pressure across the five ocean depth zones:

DEPTH ZONESEpipelagic0–200 mSunlit ZoneMesopelagic200–1000 mTwilight ZoneBathypelagic1000–4000 mMidnight ZoneAbyssopelagic4000–6000 mHadopelagic6000–11000 mTrenches onlyTemperature (Β°C)2541.5MixedlayerThermoclineOβ‚‚ (mL/L)70.55OMZLight (% surface)100%1%0%pH8.17.7Pressure (atm)11004001100

0.9 Dissolved Gas Solubility & Henry's Law

The solubility of gases in seawater is governed by Henry's law, which states that at equilibrium, the dissolved gas concentration is proportional to its partial pressure in the overlying atmosphere:

\( [\text{Gas}]_{\text{eq}} = K_H(T, S) \cdot p_{\text{gas}} \)

\(K_H\): Henry's law constant (mol/L/atm), decreasing with temperature

For O\(_2\), solubility decreases approximately 2% per Β°C of warming: at 0 Β°C, saturated seawater holds ~8 mL/L; at 30 Β°C, only ~4.5 mL/L. This temperature dependence is a key driver of ocean deoxygenation under climate change, as discussed inClimate Science Module 6.

For CO\(_2\), the solubility is enhanced by chemical reactions with water (hydration to carbonic acid), making the ocean the largest active carbon reservoir (~38,000 Gt C dissolved inorganic carbon). The temperature dependence of CO\(_2\) solubility follows:

\( \ln K_H = A_1 + \frac{A_2}{T} + A_3 \ln T + A_4 T + S\left(B_1 + \frac{B_2}{T} + B_3 T^2\right) \)

Weiss (1974) parameterisation: T in Kelvin, S in permil

The Revelle factor \(R_f \approx 10\) describes the buffer capacity: a 1% increase in atmospheric CO\(_2\) produces only a ~0.1% increase in dissolved inorganic carbon, because most added CO\(_2\) is neutralised by reaction with carbonate ions. As ocean CO\(_2\) uptake continues, \(R_f\)increases (buffer capacity decreases), reducing the ocean's ability to absorb future emissions.

Simulation: Ocean Depth Profiles

Quantitative depth profiles for temperature, dissolved oxygen, light intensity (Beer-Lambert), pressure, and pH, with compensation depth calculation and thermocline modelling:

Python
script.py167 lines

Click Run to execute the Python code

Code will be executed with Python 3 on the server

Simulation: Sound Speed Profile & SOFAR Channel

The sound speed profile showing the SOFAR channel minimum, with temperature and pressure contributions decomposed:

Python
script.py46 lines

Click Run to execute the Python code

Code will be executed with Python 3 on the server

Simulation: Thermocline Advection-Diffusion Model

Testing the Munk (1966) model: the balance between upwelling \(w\) and turbulent diffusion\(\kappa\) sets the thermocline structure. We also explore how changes in diffusivity and upwelling rate affect the temperature profile:

\( T(z) = T_{\text{deep}} + \Delta T \cdot \exp\!\left(-\frac{w}{\kappa}\,z\right) \)

Python
script.py92 lines

Click Run to execute the Python code

Code will be executed with Python 3 on the server

References

  1. Munk, W. H. (1966). Abyssal recipes. Deep-Sea Research, 13(4), 707–730.
  2. Redfield, A. C. (1934). On the proportions of organic derivatives in sea water and their relation to the composition of plankton. James Johnstone Memorial Volume, 176–192.
  3. Kirk, J. T. O. (2011). Light and Photosynthesis in Aquatic Ecosystems (3rd ed.). Cambridge University Press.
  4. Sarmiento, J. L. & Gruber, N. (2006). Ocean Biogeochemical Dynamics. Princeton University Press.
  5. Yancey, P. H., Gerringer, M. E., Drazen, J. C., Rowden, A. A., & Jamieson, A. (2014). Marine fish may be biochemically constrained from inhabiting the deepest ocean depths. PNAS, 111(12), 4461–4465.
  6. Jamieson, A. J. (2015). The Hadal Zone: Life in the Deepest Oceans. Cambridge University Press.
  7. Emerson, S. & Hedges, J. (2008). Chemical Oceanography and the Marine Carbon Cycle. Cambridge University Press.