Module 15: Ocean Currents & Marine Biodiversity

The global thermohaline circulation (THC) is the density-driven overturning that connects surface and deep oceans on timescales of 1–2 kyr. Its Atlantic limb — the Atlantic Meridional Overturning Circulation (AMOC)— transports ~1.3 PW of heat northward and sustains the mild climate of NW Europe. The RAPID-MOCHA array at 26.5°N has measured a mean overturning stream function of\(\Psi\approx 17\pm 1\;\text{Sv}\) (McCarthy et al., 2015).

The Stommel Two-Box Model

Stommel (1961) reduced the overturning to a pair of well-mixed boxes (low-latitude warm+salty, high-latitude cold+fresh) connected by a pipe whose transport responds linearly to density difference. Temperatures relax quickly to the atmosphere (\(\tau_T\ll\tau_S\)), leaving salinity as the slow variable:

The steady states \(\dot S=0\) trace out a cubic curve \(F_S(S)\)with a saddle-node bifurcation at \(S^*=\alpha\Delta T/(2\beta)\). Beyond this critical freshwater forcing, the strong-overturning branch vanishes and the system collapses onto a weak/reversed branch. The resulting hysteresis loop is the physical basis for the often-quoted “tipping point” of the AMOC (Rahmstorf, 1995; Lenton et al., 2008).

Fingerprints of AMOC Slowdown

Caesar et al. (2018) used a subpolar North Atlantic SST “cold blob” index as a proxy and concluded that the AMOC has weakened ~15% since the mid-20th century. Rahmstorf et al. (2015) arrived at similar conclusions from coral-inferred Florida Current transport. Paleo analogues such as the Younger Dryas (12.9–11.7 kyr BP) demonstrate that North Atlantic climate can undergo 5–10 °C swings within decades when freshwater hosing collapses NADW formation.

Biogeochemical & Ecological Imprint

Weakening AMOC chills the NW European shelf, redistributes nutrient supply, and alters the timing and intensity of the North Atlantic spring bloom. Drinkwater (2005) linked cod (Gadus morhua) recruitment to bottom temperature windows of 4–8 °C — precisely the range AMOC maintains on the Grand Banks. A collapsed AMOC would both cool and freshen the subpolar gyre, with cascading effects on calanoid copepod phenology (\(\textit{Calanus finmarchicus}\)) and therefore on herring, cod and right-whale prey (Head & Pepin, 2010).

Beneath the thermohaline blanket sits the vigorous wind-driven circulation: the Gulf Stream, Kuroshio, Agulhas and Brazil Current (western boundary intensification), the returning California, Canary, Benguela, Humboldt and West Australian currents (eastern boundaries), the Oyashio/Labrador subpolar returns, and the Antarctic Circumpolar Current (ACC, ~137 Sv across Drake Passage). Their shape follows from Sverdrup balance:

Ekman Pumping

The vertical velocity at the base of the surface Ekman layer — the Ekman pumping \(w_E\) — determines whether a region experiences upwelling (divergent transport, cold nutrient-rich) or downwelling (convergent, warm oligotrophic):

Along a N–S eastern boundary (e.g. Peru) the equatorward alongshore trade wind produces offshore Ekman transport of the surface layer, forced upward replacement along the coast. With \(\tau\approx 0.10\;\text{N/m}^2\) at 15°S,\(f=-3.8\times 10^{-5}\;\text{s}^{-1}\) and a coastal scale of 30 km, we get\(w_E\approx 2\;\text{m/day}\) — enough to deliver 25 µM NO\(_3\) water to the euphotic zone and sustain the most productive ecosystem on Earth. The Humboldt anchoveta(Engraulis ringens) routinely contributes ~5% of global marine fish catch from that single system (Chavez et al., 2008).

Four Eastern-Boundary Systems, One Physics

The California, Canary, Benguela and Humboldt systems all operate on the same Ekman pumping, but their biota diverge because of bathymetry, continental-shelf width and teleconnections. The Benguela is bracketed by the warm Agulhas to the south and the Angola Current to the north, giving rise to endemic sardine (Sardinops sagax) and hake (Merluccius capensis) fisheries and the iconic Spheniscus demersus (African penguin).

Stommel (1948) explained the dramatic asymmetry of ocean gyres — narrow swift western boundary currents, broad sluggish eastern returns — as the consequence of the latitudinal variation of the Coriolis parameter (\(\beta\) effect). The resulting jets (Gulf Stream 100 km wide, 1.5 m/s; Kuroshio similar) are globally conspicuous in satellite SST imagery.

Since the 1990s, altimetry has revealed that these jets shed mesoscale eddies with radii of 50–200 km and lifetimes of months. Eddies contribute an extra stirring diffusivity\(\kappa_e\sim 10^3\,{\rm m}^2/{\rm s}\) that redistributes heat, salt and phytoplankton seed communities (Chelton et al., 2011). Warm-core rings of the Gulf Stream entrain Sargassum and transport tropical larvae into temperate shelf waters, delivering invasive species to the Gulf of Maine (e.g. Hemigrapsus sanguineus). Similarly, Agulhas rings pinch off the Agulhas retroflection and cross the Atlantic carrying Indo-Pacific heat and salt — the so-called Agulhas leakage that modulates the AMOC on decadal timescales (Beal et al., 2011).

Wu et al. (2012) documented that the four principal western boundary currents have intensified and migrated poleward by ~200 km over 1900–2008, creating marine heatwave “hotspots” where warming rates are 3–4× the global average. The southeast Australian shelf under the East Australian Current has recorded the most extreme marine heatwaves (MHW), triggering kelp-forest to urchin-barren phase shifts across 500 km of coastline (Wernberg et al., 2016).

Superimposed on the mean currents is a zoo of interannual-to-decadal modes that repaint the biogeography on biological lifespans: ENSO(NINO3.4 \(\equiv\) SST anomaly 5°S–5°N, 170–120°W), the Pacific Decadal Oscillation (PDO), the Atlantic Multidecadal Oscillation(AMO), the Indian Ocean Dipole (IOD) and the Southern Annular Mode(SAM). Each is diagnosed by an empirical orthogonal function of SST or sea-level pressure.

During an El Niño the Walker circulation reverses, warm-pool water slides east, the thermocline deepens off Peru and cold nutrient-rich source water is cut off. The 1972/73 event collapsed the anchoveta from ~12 Mt catch to <2 Mt within a year, took seabird (Phalacrocorax bougainvillii, Sula variegata, Pelecanus thagus) populations from 30 M to <6 M individuals, and is widely credited with triggering the 1976–1977 Pacific regime shift (Chavez et al., 2003). The 1997/98 and 2015/16 events had similar signatures with slight shifts in the peak longitude (CP vs EP El Niño; Ashok et al., 2007).

Trenberth et al. (2014) showed that the ocean heat content budget during the 1998–2013 “warming hiatus” was dominated by increased heat uptake below 700 m in the tropical Pacific, orchestrated by the negative phase of the PDO. In short: multi-decadal ocean modes can mask or amplify radiative forcing on decadal timescales, a fact that has obvious consequences for marine ecological attribution.

Sardine–Anchovy Alternation

Chavez et al. (2003) documented a roughly 25-year alternation between sardine-dominated (warm PDO) and anchovy-dominated (cool PDO) regimes, synchronous in the California, Humboldt and Kuroshio systems. A mechanistic explanation follows from differences in thermal niche of early life stages: anchovies spawn near upwelling fronts (12–18 °C), whereas sardines spawn in warmer transition waters (16–22 °C).

Warming surface waters drive a poleward reshuffling of marine species. Cheung et al. (2013) reported a global mean shift of\(72\pm 0.36\;\text{km}\,\text{decade}^{-1}\) for marine taxa — an order of magnitude faster than terrestrial analogues. Key case studies include:

• Barents Sea: Gadus morhua and other boreal fish expanded ~160 km northward in a single decade, replacing Arctic sculpins and capelin (Fossheim et al., 2015).
• NW Atlantic: cod stocks collapsed in the 1990s and have never fully recovered under warmer shelf waters (Pershing et al., 2015 — “fastest-warming Gulf of Maine”).
• Tuna: skipjack (Katsuwonus pelamis), yellowfin (T. albacares) and bluefin (T. thynnus) are tracking warm-pool expansion; Atlantic bluefin have re-occupied Greenland waters not seen since the 1960s.
• Humpback whales (Megaptera novaeangliae) follow krill belts; acoustic monitoring shows feeding ground shifts keyed to AMOC-modulated productivity (Burtenshaw et al., 2004).
• Blue whales (Balaenoptera musculus) show phenology delays of 20–30 days in the California Current under warm PDO phases.
• Southern right whales (Eubalaena australis) shift breeding grounds following SAM-driven changes in copepod abundance (Leaper et al., 2006).
• Wandering albatross (Diomedea exulans) exploit southward-shifting westerlies to cut foraging trip length by ~20% (Weimerskirch et al., 2012).

Underlying these shifts is a “climate velocity” \(v_c=\partial T/\partial t\,\big/\,|\nabla T|\)(Loarie et al., 2009): the speed a species must travel to stay in its thermal niche. For boreal/polar species on flat continental shelves \(|\nabla T|\) is small, so\(v_c\) is large — explaining why high-latitude fish are the fastest movers in the ocean.

Currents are not only thermostats but also dispersal kernels. Cowen et al. (2006) used high-resolution Lagrangian simulations to show that reef-fish larval retention is much higher than the diffusive prediction — typically\(10\text{--}60\%\) of spawn recruit within 50 km of the natal reef. Treml et al. (2008) formalised this as a weighted directed graph: reefs are nodes, larval fluxes are edges, pelagic larval duration (PLD) is the edge cost. The resulting network reveals stepping-stone reefs that maintain genetic connectivity across the Indo-Pacific.

The dispersal kernel is approximately \(k(r)\propto \exp(-r/L_d)\) with\(L_d = u\,\mathrm{PLD}\), where \(u\) is the advecting current. For acroporid corals (\(\mathrm{PLD}\sim 14\)d) in the Coral Triangle (\(u\sim 0.2\;\mathrm{m/s}\)), \(L_d\approx 240\;\mathrm{km}\)— enough to knit Sulawesi, Halmahera and the Papuan reefs into a single metapopulation.

Mass-bleaching events recorded on the Great Barrier Reef by Hughes et al. (2018) show that\(\sim 30\%\) of coral cover was lost during the 2016 marine heatwave alone. Because the 2016 event aligned with the East Australian Current’s southward intensification, connectivity pathways that previously delivered larvae from northern GBR sources to southern recovery sinks became simultaneously stressed, collapsing the network’s resilience.

Stramma et al. (2008) documented expansion of tropical Oxygen Minimum Zones (OMZs) over 1960–2008 at the rate of \(\sim 4.5\%\) per decade in the eastern tropical Pacific. The mechanism couples warming (reduced O\(_2\) solubility,\([{\rm O}_2]\propto T^{-1}\) approximately), stratification (reduced ventilation) and enhanced export production at upwelling fronts. Billfish and tunas with high oxygen demand are squeezed into shallower, warmer layers (the “habitat compression” of Prince & Goodyear, 2006).

Brotz et al. (2012) quantified increasing jellyfish populations in 62% of 45 regional Large Marine Ecosystems. Aurelia aurita and the mauve stinger Pelagia noctiluca exploit the combination of warm surface temperatures, overfished predators and OMZ-driven fish exclusion. Mnemiopsis leidyi invasions of the Black and Caspian Seas collapsed local anchovy fisheries in the 1990s.

Gyre Accumulation & Plastics

The five subtropical gyres (N & S Pacific, N & S Atlantic, Indian) act as convergence zones. Lebreton et al. (2018) mapped the North Pacific garbage patch at\(79\pm 10\;\mathrm{kt}\) of floating plastic spread over\(1.6\times 10^6\;\mathrm{km}^2\), with 92% of the mass in fragments >5 mm. Wang et al. (2019) tracked the Great Atlantic Sargassum Belt (GASB), a \(\sim 20\)-Mt bloom fuelled by Amazon nutrient plumes coupled to the North Equatorial Recirculation.

Ocean circulation drives the biological pump that exports ~10 Gt C/yr from the euphotic zone, with 5–20% reaching the deep ocean. The Martin curve (Martin et al., 1987) describes flux attenuation with depth:

In the HNLC (high-nutrient, low-chlorophyll) regions of the Southern Ocean, equatorial Pacific and subarctic N Pacific, iron is the limiting micronutrient (Martin iron hypothesis, 1990). Antarctic krill (Euphausia superba) act as an iron recycling machine via fecal pellet production and vertical migration, and their decline under sea-ice loss (Atkinson et al., 2019) threatens the entire Southern Ocean pump.

Chlorophyll trends from MODIS-Aqua (Gregg & Rousseaux, 2019) show declines of\(\sim 1.5\%\) per decade in the subtropical gyres, consistent with enhanced stratification. These trends feed back into climate by reducing CO\(_2\)drawdown and DMS-based cloud nucleation (Charlson-Lovelock-Andreae-Warren hypothesis).

The Arctic halocline isolates cold surface water (0–50 m) from warmer Atlantic water below by a salinity step. Loss of multi-year sea ice increases surface freshwater input, strengthens the halocline, and traps heat below, thawing ice from beneath — the “Atlantification”of the Barents Sea (Polyakov et al., 2017). Ecosystem consequences include poleward advance of Atlantic copepods (Calanus finmarchicus) and fish, and decline of the Arctic-endemic Calanus glacialis — a smaller, lipid-rich species that sustains seabird chicks on limited prey loads.

In the Southern Ocean, meltwater from Antarctic ice sheets freshens the surface layer along the continental shelf, slowing AABW production (Silvano et al., 2018). The biogeochemical footprint is asymmetric with the Arctic: Southern Ocean iron limitation persists despite increased light and stratification because iron sources (dust, sea ice, melting glaciers) remain episodic. Antarctic krill response to these dual stressors is mediated through chitinase expression, lipid metabolism and reproductive phenology (Atkinson et al., 2019).

Below 2000 m the abyssal circulation is dominated by Antarctic Bottom Water (AABW) flowing northward on the western side of each basin. Benthic biodiversity peaks at bathyal depths (1000–3000 m) where lateral advection delivers particulate organic carbon from productive surface waters. Methane seep communities (Bathymodiolus mussels, siboglinid tubeworms with chemoautotrophic endosymbionts) expand as Arctic methane hydrates dissociate, providing localised oases of biomass amid oligotrophic sediments (Boetius & Wenzhöfer, 2013).

The weakening of the AABW export over the past three decades (~30% reduction of densest-class production; Purkey & Johnson, 2012) lengthens the residence time of abyssal waters and allows O\(_2\) consumption by organic matter respiration to deplete bottom-layer oxygen — a slow-motion threat to long-lived deep-sea fauna with centennial generation times (e.g. the rougheye rockfish).

Marine mammals and seabirds are ocean-current integrators on timescales of weeks to decades. Blue whales (Balaenoptera musculus) in the California Current exploit seasonal krill (Euphausia pacifica) aggregations that form where southward-flowing upwelled water meets warm subtropical water; Burtenshaw et al. (2004) combined satellite chlorophyll with acoustic calls to show that the feeding onset shifts two weeks later in warm-PDO years. Humpback whales (Megaptera novaeangliae) undertake 10\(\,000\)-km migrations from tropical calving grounds to polar feeding grounds, with both endpoints set by circulation-controlled prey abundance.

The wandering albatross (Diomedea exulans) foraging envelope extends 5\(\,000\) km from Crozet Island in a single trip, powered by dynamic soaring in the westerlies. Weimerskirch et al. (2012) showed that SAM-driven poleward strengthening of the westerlies has allowed females to cut trip length by ~20% and raise body mass, while males (which forage further south) approach the edge of their thermal niche and have stagnated. Such sex-asymmetric responses illustrate how wind-current changes can decouple reproductive physiology even within a single species.

Southern right whales (Eubalaena australis) calve off southeast Brazil, Argentina-Patagonia, South Africa and southwest Australia; calving success is tied to the Falkland (Malvinas) Current and the Agulhas Return Current productivity — both sensitive to SAM (Leaper et al., 2006). North Atlantic right whales (Eubalaena glacialis), in contrast, are tied to Calanus finmarchicuspatches shaped by the Labrador Current intrusion into the Gulf of Maine; their 2015 shift into the Gulf of St Lawrence was a direct consequence of Labrador Sea reorganisation and produced catastrophic ship-strike and entanglement mortality (Record et al., 2019).

Climate Velocity in the Sea

Combining the arguments of earlier sections, the “climate velocity” of the ocean — the speed a given isotherm moves — is given by\(v_c=\partial T/\partial t\,/\,|\nabla T|\). On the North Sea shelf, with warming of \(0.04\;{}^{\circ}\mathrm{C/yr}\) and horizontal gradients\(\sim 2\times 10^{-5}\,{}^{\circ}\mathrm{C/m}\),\(v_c\approx 65\;\mathrm{km/decade}\). For flat continental shelves like the Bering Sea, \(v_c\) exceeds 100 km/decade — faster than most fish can disperse in a generation, predicting local extirpation under even moderate emissions.

A marine heatwave (MHW) is defined by Hobday et al. (2016) as a period of at least 5 days with SST above the climatological 90th percentile. Frequency and duration of MHWs have doubled since 1982 (Oliver et al., 2018). Landmark events include:

• NE Pacific “Blob” (2013–2015): anomalies up to +3 °C sustained 2 years; triggered toxic Pseudo-nitzschia diatom bloom (domoic acid), Cassin’s auklet die-off, sea lion pup starvation, and a permanent shift of tuna and mahi distributions northward.
• GBR 2016 & 2017 bleaching: consecutive MHWs killed 50% of corals over two seasons (Hughes et al., 2018).
• Tasman Sea 2015–2016: East Australian Current intrusion; collapse of oyster farming and southward shift of East Australian kelp (Wernberg et al., 2016).
• Mediterranean 2003, 2022: mass mortality of gorgonian corals, sponges and Posidonia seagrass.

Under 2 °C warming, the IPCC AR6 projects permanent MHW conditions over\(50\%\) of the ocean by 2100 — rendering the concept of a marine heatwave meaningless relative to the pre-industrial baseline. Ecosystems that have evolved on inter-annual variability will undergo state changes: kelp forests to urchin barrens, coral reefs to algal pavements, kelp-sea-otter-urchin classic trophic cascades into novel Stellar sea-lion-kelp-urchin configurations.

Ocean currents integrate radiative forcing, wind patterns, freshwater fluxes and bathymetry into a three-dimensional delivery system for heat, salt, oxygen, nutrients and larvae. Any one of the processes we have surveyed — AMOC slowdown, expanding OMZs, intensifying eastern-boundary upwelling, Sargassum belt formation, larval connectivity collapse — would by itself reshape the marine biosphere over the next century. In combination, they constitute a compound forcing far larger than the sum of its parts.

Three cross-cutting lessons emerge:

• Bistability is generic. Stommel’s two-box argument is the simplest of an entire family of nonlinear ocean-biogeochemical oscillators (ENSO, PDO, AMO, NAO, IOD, SAM). Biology sits on top of these attractors and often bistability appears in food-webs themselves (sardine/anchovy, jellyfish/fish regimes).
• Scale-couplings are fast in the ocean. Fish and plankton disperse on timescales of days, whales on seasons, currents on years; climate velocity on the shelf commonly exceeds dispersal capacity, creating range-shift debt.
• Coastal upwelling supports disproportionate biodiversity and human protein supply. Even small changes in wind-stress curl via jet stream shifts or monsoon strengthening propagate through the biological pump to cetaceans and seabirds.

This module therefore sets up the molecular toolkit needed by Module 16: every organism riding an altered current must deploy biochemical defences — HSPs for heat, AFPs for freezing, osmolytes for salinity, photoprotective flavonoids for light — and those defences determine where the new range limits will settle.