4.1 Marine Ecosystems

The Living Ocean

Marine ecosystems encompass the full range of life in the ocean, from microscopic bacteria and archaea to the largest animals that have ever lived. The ocean produces approximately 50 Gt C per year through photosynthesis—roughly half of global primary production—despite containing only ~1-2% of the planet's photosynthetic biomass. This extraordinary turnover reflects the rapid growth rates of phytoplankton compared to terrestrial plants.

Marine biomes span coastal wetlands, coral reefs, open ocean gyres, upwelling systems, polar seas, and the deep sea. Each is shaped by the interplay of light, nutrients, temperature, and circulation, creating distinct ecological provinces with characteristic community structures and productivity regimes.

Marine Biomes

Coastal & Shelf Systems

High productivity (150-300 g C/m²/yr) from nutrient inputs via rivers, upwelling, and tidal mixing. Estuaries, mangroves, seagrass beds, and kelp forests provide nursery habitats. Only 10% of ocean area but >30% of marine production.

Open Ocean (Oligotrophic Gyres)

Low productivity (30-60 g C/m²/yr) due to strong stratification limiting nutrient supply. Dominated by picophytoplankton (Prochlorococcus, Synechococcus). "Marine deserts" covering ~40% of Earth's surface. Deep chlorophyll maximum at 80-150 m.

Polar Ecosystems

Extreme seasonality: 24-hr light/dark cycles. Intense spring blooms when sea ice retreats and stratification develops. High productivity (100-300 g C/m²/yr in bloom season). Krill-based food webs supporting whales, seals, and penguins.

Tropical Coral Reefs

Extremely high biodiversity (>25% of marine species). Productivity: 1000-5000 g C/m²/yr (ecosystem basis). Symbiosis between corals and zooxanthellae. Calcium carbonate framework builds three-dimensional habitat. Threatened by warming and acidification.

Primary Production

Marine primary production is the fixation of inorganic carbon into organic matter through photosynthesis. The photosynthesis equation, including the Redfield stoichiometry for nutrients consumed, is:

$$106\text{CO}_2 + 16\text{NO}_3^- + \text{HPO}_4^{2-} + 122\text{H}_2\text{O} + 18\text{H}^+ \xrightarrow{h\nu} (\text{CH}_2\text{O})_{106}(\text{NH}_3)_{16}\text{H}_3\text{PO}_4 + 138\text{O}_2$$

The rate of primary production depends on light and nutrient availability. A commonly used formulation combines light limitation (Webb et al.) with nutrient limitation (Monod):

$$P = P_{\max} \left(1 - e^{-\alpha I / P_{\max}}\right) \cdot \frac{N}{K_N + N}$$

$P_{\max}$ = maximum photosynthetic rate, $\alpha$ = initial slope,$I$ = light intensity, $K_N$ = half-saturation for nutrients

~50 Gt C/yr

Total marine NPP

~50%

Of global photosynthesis

~1-2 Gt C

Total phytoplankton biomass

Compensation Depth & Sverdrup's Critical Depth

The compensation depth $z_c$ is where photosynthesis equals respiration for a single cell. Below this depth, a phytoplankton cell cannot sustain itself:

$$P(z_c) = R \quad \Rightarrow \quad I(z_c) = I_c$$

Typically at 1% of surface light level, i.e., $z_c \approx 4.6/K_d$ where $K_d$ is the diffuse attenuation coefficient

Sverdrup (1953) showed that a net spring bloom can only develop when the mixed layer depth is shallower than the critical depth $z_{cr}$, defined as the depth where depth-integrated production equals depth-integrated respiration:

$$\int_0^{z_{cr}} P(z)\,dz = \int_0^{z_{cr}} R\,dz$$

Spring bloom initiates when mixed layer depth < critical depth (shoaling due to warming/calming)

Oligotrophic Systems

Subtropical gyres. Strong stratification year-round. Low nutrients, low biomass but steady production. Dominated by small cells (pico/nano). f-ratio ~ 0.1.

Eutrophic Systems

Upwelling regions, coastal waters. Abundant nutrients. Large diatom blooms. High biomass, high export. f-ratio ~ 0.3-0.5. Strong seasonality at high latitudes.

Ecological Provinces (Longhurst)

Longhurst (1998, 2007) divided the ocean into 56 biogeochemical provinces based on physical forcing, nutrient regimes, and biological response. Major province types include:

Trade Wind Biome

Oligotrophic gyres. Permanent stratification. Deep nutricline. Low chlorophyll.

Westerlies Biome

Seasonal mixed layer dynamics. Spring blooms. Moderate productivity.

Polar Biome

Ice-influenced. Extreme seasonality. Short intense blooms. High export.

Coastal Biome

Upwelling-driven. Riverine nutrient input. Highest areal productivity.

Derivation: Net Ecosystem Production & Critical Depth

Step 1: Define Production and Respiration Components

Gross Primary Production (GPP) is the total carbon fixed by photosynthesis. Net Primary Production (NPP) subtracts autotrophic respiration R_a. Net Ecosystem Production (NEP) further subtracts heterotrophic respiration R_h:

$$\text{GPP} = \text{NPP} + R_a, \qquad \text{NEP} = \text{NPP} - R_h = \text{GPP} - R_a - R_h$$

Step 2: Light-Dependent Production with Depth

Light attenuates exponentially with depth following Beer-Lambert law. The photosynthesis rate at depth z using the Webb (1974) formulation is:

$$P(z) = P_{\max}\!\left(1 - e^{-\alpha I_0 e^{-K_d z}/P_{\max}}\right)$$

where I(z) = I₀ exp(-K_d z) is the irradiance at depth z, K_d is the diffuse attenuation coefficient, alpha is the photosynthetic efficiency, and P_max is the maximum photosynthetic rate.

Step 3: Compensation Depth from P(z) = R

At the compensation depth z_c, photosynthesis equals respiration for a single cell. Setting P(z_c) = R and solving for z_c:

$$P_{\max}\!\left(1 - e^{-\alpha I_0 e^{-K_d z_c}/P_{\max}}\right) = R$$

$$z_c = -\frac{1}{K_d}\ln\!\left[-\frac{P_{\max}}{\alpha I_0}\ln\!\left(1 - \frac{R}{P_{\max}}\right)\right]$$

Step 4: Sverdrup Critical Depth Integral

Sverdrup (1953) defined the critical depth z_cr as the depth where depth-integrated production equals depth-integrated community respiration R_c over the mixed layer. A bloom can develop only when the mixed layer depth (MLD) is shallower than z_cr:

$$\int_0^{z_{cr}} P(z)\,dz = \int_0^{z_{cr}} R_c\,dz = R_c \cdot z_{cr}$$

Step 5: Evaluate the Integral for Critical Depth

For the simplified linear P-I model P(z) = alpha I₀ exp(-K_d z) (valid when light is not saturating), the integral evaluates analytically:

$$\int_0^{z_{cr}} \alpha I_0 e^{-K_d z}\,dz = \frac{\alpha I_0}{K_d}\left(1 - e^{-K_d z_{cr}}\right) = R_c \cdot z_{cr}$$

This transcendental equation is solved numerically for z_cr. The bloom condition is then: MLD < z_cr, which is satisfied when the mixed layer shoals in spring due to warming and reduced wind mixing.

Derivation: Energy Flow Through Marine Ecosystems

Step 1: Trophic Transfer Efficiency

Energy flows from one trophic level to the next with an efficiency epsilon. At each level, energy is lost to respiration, excretion, and egestion. The transfer efficiency from trophic level n to n+1 is:

$$\epsilon = \frac{P_{n+1}}{P_n} = \frac{\text{Production at level } n{+}1}{\text{Production at level } n} \approx 0.10$$

Step 2: Production at Trophic Level n

Starting from NPP at the base, the production at trophic level n is a geometric decay. This is the "10 percent rule" of Lindeman (1942):

$$P_n = \text{NPP} \cdot \epsilon^{n-1}$$

Step 3: Energy Budget for a Single Trophic Level

The energy budget for any organism partitions ingested energy (I) among assimilation (A), egestion (F), production (P), and respiration (R):

$$I = F + A, \qquad A = P + R$$

$$\Rightarrow \quad I = P + R + F, \qquad \text{Gross growth efficiency} = \frac{P}{I} = \frac{A - R}{I}$$

Step 4: Export Production and the Biological Pump

The fraction of NPP that sinks below the euphotic zone as export production (EP) is related to the f-ratio. The efficiency of the biological pump in sequestering carbon to depth z is:

$$\text{EP} = f \cdot \text{NPP}, \qquad \text{Transfer efficiency}(z) = \frac{F(z)}{F(z_0)} = \left(\frac{z}{z_0}\right)^{-b}$$

Only ~1% of surface NPP reaches the deep sea floor (4000 m), with most organic matter remineralized in the upper 1000 m (mesopelagic zone).

Python: Sverdrup Critical Depth & Spring Bloom NPZ

Python: Sverdrup Critical Depth & Spring Bloom NPZ

Python

!/usr/bin/env python3

script.py83 lines

Click Run to execute the Python code

Code will be executed with Python 3 on the server

Fortran: Primary Production Model

Fortran: Primary Production Model

Fortran

-------------------------------------------------------

program.f9077 lines

Click Run to execute the Fortran code

Code will be compiled with gfortran and executed on the server

Key Takeaways

  • Marine primary production (~50 Gt C/yr) accounts for roughly half of global photosynthesis.
  • Sverdrup's critical depth hypothesis explains the onset of spring blooms at high latitudes.
  • Production is co-limited by light and nutrients, with the dominant limiter varying by region and season.
  • Longhurst ecological provinces classify the ocean into regions with similar biological responses.
  • Coastal upwelling regions have the highest areal productivity despite covering a small fraction of the ocean.