Module 6

Coral Reef Biochemistry

Zooxanthellae symbiosis, calcification, bleaching, fluorescence, and reef nitrogen cycling

Coral reefs occupy less than 0.1% of the ocean floor yet support approximately 25% of all marine species. This extraordinary productivity arises from an intimate biochemical partnership between animal coral polyps and photosynthetic dinoflagellate algae (Symbiodiniaceae), operating in some of the most nutrient-poor waters on Earth. In this module we derive the quantitative biochemistry of coral symbiosis, aragonite calcification, thermal bleaching, and the tight nutrient recycling that resolves Darwin's paradox. See also Climate & Biodiversity and Climate Science for broader context on ocean-climate interactions.

6.1 Zooxanthellae Symbiosis

Hermatypic (reef-building) corals harbour intracellular photosynthetic dinoflagellates of the family Symbiodiniaceae (formerly genus Symbiodinium, colloquially “zooxanthellae”). These single-celled algae reside within membrane-bound vesicles called symbiosomes inside the coral's gastrodermal cells, at densities of 1–2 × 10⁶ cells per cm² of coral tissue.

Photosynthate Translocation

Zooxanthellae fix carbon via the Calvin cycle, producing glycerol, glucose, and amino acids. Remarkably, up to 90–95% of photosynthate is translocated to the coral host. In return, the coral provides CO₂, inorganic nitrogen (NH₄⁺), and phosphorus from its metabolic waste, along with a stable, light-rich environment.

Net photosynthesis of the holobiont:

\( P_{\text{net}} = P_{\text{gross}} - R_{\text{coral}} - R_{\text{zoox}} \)

where \( P_{\text{gross}} \) is total photosynthetic carbon fixation, \( R_{\text{coral}} \) is coral host respiration, and \( R_{\text{zoox}} \) is zooxanthellae respiration. For a healthy reef coral at saturating irradiance, \( P_{\text{gross}} / (R_{\text{coral}} + R_{\text{zoox}}) \approx 3\text{--}5 \), meaning the symbiosis is strongly net autotrophic.

Photosynthesis-Irradiance (P-I) Curve

Zooxanthellae photosynthesis follows a saturating P-I relationship described by the hyperbolic tangent model (Jassby & Platt 1976):

\( P(I) = P_{\max} \cdot \tanh\!\left(\frac{\alpha \cdot I}{P_{\max}}\right) - R_d \)

where \( P_{\max} \) is the maximum photosynthetic rate (~5–15 μmol O₂ cm⁻² h⁻¹), \( \alpha \) is the initial light-use efficiency (~0.02–0.05 μmol O₂ μmol photons⁻¹), \( I \) is photosynthetically active radiation (PAR), and \( R_d \) is dark respiration. The compensation irradiance \( I_c \) where \( P(I_c) = 0 \)is typically 10–50 μmol photons m⁻² s⁻¹.

The saturation irradiance \( I_k = P_{\max}/\alpha \) typically falls around 200–400 μmol photons m⁻² s⁻¹. Beyond \( I_k \), excess excitation energy must be dissipated via non-photochemical quenching (NPQ) or xanthophyll cycling to prevent photodamage — a process central to the bleaching mechanism discussed in Section 6.3.

6.2 Aragonite Calcification

Scleractinian corals precipitate calcium carbonate in the aragonite polymorph, building massive skeletal structures that form the three-dimensional framework of reef ecosystems. Calcification occurs in the sub-ectodermal space between the calicoblastic ectoderm and the existing skeleton, within a thin (10–100 μm) extracellular calcifying fluid (ECF).

Aragonite Saturation State

Saturation state:

\( \Omega_{\text{arag}} = \frac{[\text{Ca}^{2+}][\text{CO}_3^{2-}]}{K_{sp}^{\text{arag}}} \)

where \( K_{sp}^{\text{arag}} \) is the stoichiometric solubility product for aragonite. When \( \Omega_{\text{arag}} > 1 \), seawater is supersaturated and thermodynamically favours CaCO₃ precipitation. When \( \Omega_{\text{arag}} < 1 \), dissolution dominates.

Current tropical surface ocean \( \Omega_{\text{arag}} \approx 3.5 \). Under the high-emissions scenario SSP5-8.5, projections show a decline to \( \Omega_{\text{arag}} \approx 2.0 \) by 2100 due to ocean acidification (CO₂ absorption lowering [CO₃²⁻]).

Burton-Walter Calcification Equation

Calcification rate:

\( G = k(\Omega_{\text{arag}} - 1)^n \)

where \( G \) is the calcification rate (mmol CaCO₃ m⁻² d⁻¹), \( k \) is a rate constant dependent on temperature and biology, and \( n \) is the reaction order (typically 1.0–1.7 for corals). As \( \Omega \to 1 \), calcification slows dramatically, and below ~2.0 many coral species can no longer maintain positive net calcification.

The calcification process also requires active biological control: Ca²⁺-ATPase pumps concentrate calcium in the ECF to 2–3 times seawater levels, while proton pumps raise pH from ~8.0 to ~8.6, elevating the local \( \Omega \) to 10–25 at the site of crystal nucleation. This “biological uplift” explains why coral calcification rates exceed abiotic precipitation rates by an order of magnitude.

6.3 Coral Bleaching Mechanism

Coral bleaching is the breakdown of the coral-zooxanthellae symbiosis under thermal stress. The mechanistic cascade proceeds as follows:

Thermal damage to PSII: Elevated SST (>1°C above maximum monthly mean for >4 weeks) damages the D1 protein in photosystem II (PSII) and impairs the Calvin cycle enzyme Rubisco activase.
Excess excitation energy: With PSII repair rate overwhelmed, excess light energy cannot be utilized photochemically.
ROS production: Excess energy is transferred to O₂, generating reactive oxygen species (superoxide O₂⁻, hydrogen peroxide H₂O₂, hydroxyl radical OH•, singlet oxygen ¹O₂).
Oxidative stress cascade: ROS overwhelm antioxidant defences (SOD, catalase, ascorbate peroxidase), causing lipid peroxidation and protein damage.
Symbiont expulsion: Coral cells expel zooxanthellae via exocytosis, apoptosis, or host-cell detachment. Loss of pigmented symbionts reveals the white aragonite skeleton — “bleaching.”

Degree Heating Weeks (DHW)

NOAA Coral Reef Watch quantifies cumulative thermal stress using Degree Heating Weeks:

\( \text{DHW} = \frac{1}{7}\sum_{i=1}^{84} \max\!\bigl(\text{SST}_i - \text{MMM}, \; 0\bigr) \)

where SSTₐ is daily sea surface temperature over the preceding 84 days (12 weeks), and MMM is the maximum monthly mean climatology (the warmest month's average SST). The factor 1/7 converts degree-days to degree-weeks.

Bleaching thresholds:

DHW > 4 °C-weeks: Significant bleaching likely (coral stress alert level 1)
DHW > 8 °C-weeks: Severe bleaching and mortality likely (alert level 2)
DHW > 12 °C-weeks: Catastrophic mortality (>50% of corals die)

Global mass bleaching events have accelerated dramatically: 1998, 2010, 2014–2017, 2020, and 2023–2024 (the most severe on record). The Great Barrier Reef has experienced seven mass bleaching events since 2016. Under 1.5°C warming, 70–90% of tropical corals are projected to decline; at 2°C, loss exceeds 99% (IPCC SR1.5).

6.4 Coral Fluorescence

Many coral species express GFP-like fluorescent proteins (FPs) that produce vivid green, cyan, red, and orange fluorescence. These proteins are homologous to GFP from the jellyfish Aequorea victoria but have evolved independently in corals over >500 Myr.

Stokes Shift & Quantum Yield

Stokes shift:

\( \Delta\bar{\nu} = \bar{\nu}_{\text{abs}} - \bar{\nu}_{\text{em}} = \frac{1}{\lambda_{\text{abs}}} - \frac{1}{\lambda_{\text{em}}} \)

Typical coral FPs have Stokes shifts of 20–80 nm (e.g., cyan FP: excitation ~450 nm, emission ~485 nm, shift ~35 nm).

Fluorescence quantum yield:

\( \Phi_F = \frac{k_r}{k_r + k_{nr}} = \frac{\text{photons emitted}}{\text{photons absorbed}} \)

where \( k_r \) is the radiative decay rate and \( k_{nr} \) is the sum of all non-radiative decay pathways. Coral FPs typically have \( \Phi_F \approx 0.3\text{--}0.8 \), comparable to commercial fluorescent dyes.

Proposed functions of coral fluorescence remain debated:

Photoprotection: Shallow-water corals express high FP concentrations that absorb harmful UV/blue light and re-emit it at longer, less energetic wavelengths, acting as a biological sunscreen.
Light enhancement: Deep-water corals (>60 m) express FPs that convert blue/green wavelengths (which penetrate deepest) to orange/red wavelengths, potentially enhancing photosynthesis for zooxanthellae in light-limited conditions.
Antioxidant activity: Some FPs can quench reactive oxygen species (ROS), potentially reducing oxidative stress during thermal events.

6.5 Reef Framework & Bioerosion

The reef framework is constructed by three main coral morphologies: massive corals (Porites, Orbicella) providing bulk structure, branching corals (Acropora, Pocillopora) creating three-dimensional habitat complexity, and encrusting corals and coralline algae cementing the framework together.

Net Reef Accretion

\( \frac{dh}{dt} = G_{\text{calc}} - E_{\text{bio}} - D_{\text{chem}} - E_{\text{phys}} \)

where \( G_{\text{calc}} \) is gross calcification (typically 4–10 kg CaCO₃ m⁻² yr⁻¹ for healthy reefs), \( E_{\text{bio}} \) is bioerosion by boring sponges (Cliona), parrotfish scraping, and sea urchin (Diadema) grazing, \( D_{\text{chem}} \) is chemical dissolution (enhanced by acidification), and \( E_{\text{phys}} \) is physical breakage from storms.

On a healthy Indo-Pacific reef, net accretion rates average 3–5 mm yr⁻¹ (vertical reef growth). However, many Caribbean and Pacific reefs are now at net-zero or negative accretion due to coral mortality (reduced \( G_{\text{calc}} \)), increased bioerosion (parrotfish shifting to more aggressive scraping of dead substrate), and ocean acidification (increasing \( D_{\text{chem}} \)).

6.6 Reef Nitrogen Cycling & Darwin's Paradox

Darwin's paradox: How do coral reefs sustain extraordinary biomass and productivity in nutrient-poor tropical waters where dissolved inorganic nitrogen (DIN) is typically <0.5 μmol L⁻¹ and phosphate <0.1 μmol L⁻¹?

The answer lies in tight nutrient recycling within the reef community:

Internal recycling: Coral hosts recycle ammonium (NH₄⁺) from protein catabolism directly to zooxanthellae, which assimilate it via glutamine synthetase / glutamate synthase (GS/GOGAT) pathway. Recycling efficiency exceeds 95%.
N₂ fixation: Diazotrophic cyanobacteria (Trichodesmium) and endolithic bacteria fix atmospheric N₂ at rates of 0.1–5.0 mmol N m⁻² d⁻¹, providing new nitrogen input.
Sponge loop: Dissolved organic matter (DOM) released by corals and algae is consumed by sponges, which shed cellular detritus that feeds detritivores — converting DOM (inaccessible to most metazoans) into particulate organic matter (POM).
Denitrification balance: Anaerobic sediments within the reef framework support denitrification (NO₃⁻ → N₂), preventing eutrophication while maintaining the nitrogen-limited state that favours corals over macroalgae.

Nutrient retention efficiency:

\( \eta_N = 1 - \frac{F_{\text{export}}}{F_{\text{input}} + F_{\text{fixation}}} > 0.95 \)

This near-complete retention of nitrogen means the reef acts as a largely closed system for nutrients, even though it sits in an open-ocean “desert.” The total nitrogen residence time within the reef is estimated at 10–100 years.

Coral Polyp Cross-Section & Bleaching Sequence

Computational Models

Four integrated simulations: (1) calcification rate as a function of aragonite saturation state, (2) degree heating weeks bleaching threshold model, (3) photosynthesis-irradiance curves for zooxanthellae, and (4) reef accretion vs dissolution under ocean acidification scenarios.

Python

script.py180 lines

import numpy as np
import matplotlib
matplotlib.use('Agg')
import matplotlib.pyplot as plt

fig, axes = plt.subplots(2, 2, figsize=(14, 11))
fig.patch.set_facecolor('#0a0a0a')
fig.suptitle('Coral Reef Biochemistry Models', fontsize=18, color='#7dd3fc', fontweight='bold', y=0.98)

# ========== Panel 1: Calcification vs Omega ==========
ax1 = axes[0, 0]

Omega = np.linspace(0.5, 5.0, 200)

# Burton-Walter: G = k * (Omega - 1)^n
# Different coral types
params = [
    ('Massive (Porites)', 8.0, 1.0, '#60a5fa'),
    ('Branching (Acropora)', 15.0, 1.4, '#22c55e'),
    ('Encrusting (coralline)', 5.0, 1.2, '#f472b6'),
]

for label, k, n, color in params:
    G = np.where(Omega > 1, k * (Omega - 1)**n, -2 * (1 - Omega))
    ax1.plot(Omega, G, color=color, linewidth=2.5, label=f'{label} (k={k}, n={n})')

# Mark current and projected Omega
ax1.axvline(3.5, color='#22c55e', linestyle='--', alpha=0.7, linewidth=1.5)
ax1.text(3.55, 25, 'Current
(Ω=3.5)', fontsize=9, color='#22c55e')
ax1.axvline(2.0, color='#ef4444', linestyle='--', alpha=0.7, linewidth=1.5)
ax1.text(2.05, 25, '2100 SSP5
(Ω=2.0)', fontsize=9, color='#ef4444')
ax1.axhline(0, color='#6b7280', linewidth=1, linestyle=':')
ax1.axvline(1.0, color='#fbbf24', linewidth=1, linestyle=':', alpha=0.5)
ax1.text(1.05, -3, 'Ω=1', fontsize=9, color='#fbbf24')

ax1.set_xlabel('Ω aragonite', fontsize=12, color='white')
ax1.set_ylabel('Calcification rate G (mmol CaCO₃ m⁻² d⁻¹)', fontsize=10, color='white')
ax1.set_title('Burton-Walter Calcification Model', fontsize=14, color='#7dd3fc', fontweight='bold')
ax1.set_facecolor('#0a0a0a')
ax1.tick_params(colors='white')
ax1.grid(True, alpha=0.15, color='white')
ax1.legend(fontsize=9, facecolor='#1a1a1a', edgecolor='#3b82f6', labelcolor='white')
ax1.set_xlim(0.5, 5.0)
ax1.set_ylim(-5, 35)

# ========== Panel 2: DHW Bleaching Model ==========
ax2 = axes[0, 1]

np.random.seed(42)
days = np.arange(365)
# Baseline seasonal SST cycle (tropical reef)
MMM = 29.0  # Maximum monthly mean
SST_base = 27.5 + 2.0 * np.sin(2 * np.pi * (days - 60) / 365)

# Add marine heatwave
heatwave_start = 50
heatwave_duration = 60
heatwave = np.zeros(365)
heatwave[heatwave_start:heatwave_start + heatwave_duration] = 1.5 * np.sin(
    np.pi * np.arange(heatwave_duration) / heatwave_duration) + 0.8
SST = SST_base + heatwave + np.random.normal(0, 0.15, 365)

ax2.plot(days, SST, color='#60a5fa', linewidth=1.5, label='SST', alpha=0.8)
ax2.axhline(MMM, color='#fbbf24', linewidth=2, linestyle='--', label=f'MMM = {MMM}°C')
ax2.axhline(MMM + 1, color='#ef4444', linewidth=1.5, linestyle=':', label='Bleaching threshold')

# Fill above MMM
ax2.fill_between(days, SST, MMM, where=SST > MMM, alpha=0.25, color='#ef4444')

# Calculate DHW on secondary axis
ax2_twin = ax2.twinx()
DHW = np.zeros(365)
for d in range(84, 365):
    hotspots = np.maximum(SST[d-84:d] - MMM, 0)
    DHW[d] = np.sum(hotspots) / 7.0
ax2_twin.plot(days, DHW, color='#f472b6', linewidth=2.5, label='DHW')
ax2_twin.axhline(4, color='#fbbf24', linewidth=1, linestyle=':', alpha=0.7)
ax2_twin.axhline(8, color='#ef4444', linewidth=1, linestyle=':', alpha=0.7)
ax2_twin.text(340, 4.3, 'Alert 1', fontsize=9, color='#fbbf24')
ax2_twin.text(340, 8.3, 'Alert 2', fontsize=9, color='#ef4444')
ax2_twin.set_ylabel('DHW (°C-weeks)', fontsize=11, color='#f472b6')
ax2_twin.tick_params(colors='#f472b6')

ax2.set_xlabel('Day of year', fontsize=12, color='white')
ax2.set_ylabel('SST (°C)', fontsize=12, color='#60a5fa')
ax2.set_title('Degree Heating Weeks (DHW) Model', fontsize=14, color='#7dd3fc', fontweight='bold')
ax2.set_facecolor('#0a0a0a')
ax2.tick_params(colors='white')
ax2.grid(True, alpha=0.15, color='white')

lines1, labels1 = ax2.get_legend_handles_labels()
lines2, labels2 = ax2_twin.get_legend_handles_labels()
ax2.legend(lines1 + lines2, labels1 + labels2, fontsize=9, facecolor='#1a1a1a',
           edgecolor='#3b82f6', labelcolor='white', loc='upper left')

# ========== Panel 3: Photosynthesis-Irradiance Curves ==========
ax3 = axes[1, 0]

I = np.linspace(0, 1200, 500)  # PAR

# P-I parameters for different zooxanthellae clades
clades = [
    ('Clade C (shallow)', 12.0, 0.04, 2.0, '#22c55e'),
    ('Clade D (heat-tolerant)', 8.0, 0.03, 1.5, '#f59e0b'),
    ('Clade A (high-light)', 15.0, 0.035, 3.0, '#60a5fa'),
]

for label, Pmax, alpha, Rd, color in clades:
    P = Pmax * np.tanh(alpha * I / Pmax) - Rd
    ax3.plot(I, P, color=color, linewidth=2.5, label=label)

# Mark Ik
    Ik = Pmax / alpha
    P_Ik = Pmax * np.tanh(alpha * Ik / Pmax) - Rd
    ax3.plot(Ik, P_Ik, 'o', color=color, markersize=6)
    ax3.annotate(f'I_k={Ik:.0f}', (Ik, P_Ik), textcoords="offset points",
                xytext=(10, -10), fontsize=8, color=color)

ax3.axhline(0, color='#6b7280', linewidth=1, linestyle=':')
ax3.fill_between(I, -4, 0, alpha=0.05, color='#ef4444')
ax3.text(600, -1.5, 'Net respiration zone', fontsize=9, color='#ef4444', alpha=0.7)

ax3.set_xlabel('PAR (μmol photons m⁻² s⁻¹)', fontsize=12, color='white')
ax3.set_ylabel('P_net (μmol O₂ cm⁻² h⁻¹)', fontsize=11, color='white')
ax3.set_title('Zooxanthellae P-I Curves by Clade', fontsize=14, color='#7dd3fc', fontweight='bold')
ax3.set_facecolor('#0a0a0a')
ax3.tick_params(colors='white')
ax3.grid(True, alpha=0.15, color='white')
ax3.legend(fontsize=9, facecolor='#1a1a1a', edgecolor='#3b82f6', labelcolor='white')
ax3.set_ylim(-4, 16)

# ========== Panel 4: Reef Accretion Under Acidification ==========
ax4 = axes[1, 1]

years = np.arange(2020, 2105)

# Omega projections under different scenarios
Omega_ssp126 = 3.5 - 0.005 * (years - 2020)
Omega_ssp245 = 3.5 - 0.012 * (years - 2020)
Omega_ssp585 = 3.5 - 0.018 * (years - 2020)

scenarios = [
    ('SSP1-2.6', Omega_ssp126, '#22c55e'),
    ('SSP2-4.5', Omega_ssp245, '#f59e0b'),
    ('SSP5-8.5', Omega_ssp585, '#ef4444'),
]

# Net accretion = calcification - bioerosion - dissolution
k_calc = 10.0
n_calc = 1.2
bioerosion = 3.5  # kg CaCO3 m-2 yr-1 (constant)
dissolution_base = 1.0

for label, Omega_t, color in scenarios:
    calc = k_calc * np.maximum(Omega_t - 1, 0)**n_calc
    dissolution = dissolution_base * np.exp(0.3 * (3.5 - Omega_t))
    net_accretion = calc - bioerosion - dissolution

ax4.plot(years, net_accretion, color=color, linewidth=2.5, label=label)
    ax4.fill_between(years, net_accretion, 0, where=net_accretion < 0, alpha=0.1, color=color)

ax4.axhline(0, color='#fbbf24', linewidth=2, linestyle='--')
ax4.text(2025, 0.5, 'Net accretion ↑', fontsize=10, color='#22c55e')
ax4.text(2025, -0.5, 'Net erosion ↓', fontsize=10, color='#ef4444')

ax4.set_xlabel('Year', fontsize=12, color='white')
ax4.set_ylabel('Net reef accretion (kg CaCO₃ m⁻² yr⁻¹)', fontsize=10, color='white')
ax4.set_title('Reef Accretion Under Acidification Scenarios', fontsize=14, color='#7dd3fc', fontweight='bold')
ax4.set_facecolor('#0a0a0a')
ax4.tick_params(colors='white')
ax4.grid(True, alpha=0.15, color='white')
ax4.legend(fontsize=9, facecolor='#1a1a1a', edgecolor='#3b82f6', labelcolor='white')

plt.tight_layout(rect=[0, 0, 1, 0.96])
plt.savefig('output.png', dpi=150, bbox_inches='tight', facecolor='#0a0a0a')
plt.close()
print("Coral reef biochemistry models generated successfully.")

Click Run to execute the Python code

Code will be executed with Python 3 on the server

References

Muscatine, L. & Porter, J.W. (1977). Reef corals: mutualistic symbioses adapted to nutrient-poor environments. BioScience, 27(7), 454–460.
LaJeunesse, T.C. et al. (2018). Systematic revision of Symbiodiniaceae highlights the antiquity and diversity of coral endosymbionts. Current Biology, 28(16), 2570–2580.
Burton, E.A. & Walter, L.M. (1987). Relative precipitation rates of aragonite and Mg calcite from seawater. Geology, 15(2), 111–114.
Hoegh-Guldberg, O. (1999). Climate change, coral bleaching and the future of the world's coral reefs. Marine and Freshwater Research, 50(8), 839–866.
Liu, G. et al. (2014). Reef-scale thermal stress monitoring of coral ecosystems: new 5-km global products from NOAA Coral Reef Watch. Remote Sensing, 6(11), 11579–11606.
Alieva, N.O. et al. (2008). Diversity and evolution of coral fluorescent proteins. PLoS ONE, 3(7), e2680.
Perry, C.T. et al. (2018). Loss of coral reef growth capacity to track future increases in sea level. Nature, 558, 396–400.
de Goeij, J.M. et al. (2013). Surviving in a marine desert: the sponge loop retains resources within coral reefs. Science, 342(6154), 108–110.
Jassby, A.D. & Platt, T. (1976). Mathematical formulation of the relationship between photosynthesis and light for phytoplankton. Limnology and Oceanography, 21(4), 540–547.
IPCC (2018). Special Report on Global Warming of 1.5°C. Cambridge University Press.

M5. Hydrothermal & Chemosynthesis M7. Polar Ocean Life