Module 3

Bathypelagic (1–4 km) — The Midnight Zone

Complete darkness, crushing pressure, and life's ingenious adaptations at the frontier of the deep ocean

3.1 Complete Darkness & the Marine Snow Economy

Below 1,000 metres, the ocean is in perpetual darkness. Not a single photon of sunlight penetrates to the bathypelagic zone. Photosynthesis is impossible. Every organism here depends on energy that originates at the surface and sinks downward as marine snow — a continuous rain of dead phytoplankton, faecal pellets, mucous aggregates, and other organic detritus.

The flux of particulate organic carbon (POC) decreases dramatically with depth. The empirical relationship discovered by Martin et al. (1987) is known as the Martin curve:

$ F(z) = F_{100} \left(\frac{z}{100}\right)^{-b} $

where $F_{100}$ is the POC flux at the base of the euphotic zone (100 m),$z$ is depth in metres, and $b \approx 0.858$ is the Martin exponent. At 1,000 m depth:

$ \frac{F(1000)}{F(100)} = \left(\frac{1000}{100}\right)^{-0.858} = 10^{-0.858} \approx 0.139 $

So only about 14% of the surface export production reaches 1,000 m. By 4,000 m, just $\sim 3\%$ remains. The rest has been consumed by mesopelagic organisms on the way down.

Carbon Residence Time in the Deep Ocean

Carbon that reaches the bathypelagic zone is effectively sequestered from the atmosphere for extremely long periods. The residence time of dissolved inorganic carbon (DIC) in the deep ocean is determined by the thermohaline circulation overturning time:

$ \tau_{\text{deep}} = \frac{M_{\text{deep}}}{Q} \approx \frac{1.0 \times 10^{18}\;\text{m}^3}{30 \times 10^{6}\;\text{m}^3/\text{s}} \approx 1000\;\text{years} $

where $M_{\text{deep}}$ is the volume of the deep ocean below 1 km and $Q$is the volumetric flow rate of the global thermohaline overturning circulation ($\sim 30$ Sv). This 1,000-year timescale is why the biological carbon pump is so important for climate regulation: carbon exported to the deep ocean is locked away for a millennium.

Alternative Energy: Chemosynthesis

Some bathypelagic organisms do not rely on surface-derived organic matter at all. Near hydrothermal vents and cold seeps, chemosynthetic bacteria fix carbon using chemical energy from reduced compounds (H₂S, CH₄, H₂). The free energy available from sulfide oxidation is:

$ \text{H}_2\text{S} + 2\text{O}_2 \rightarrow \text{SO}_4^{2-} + 2\text{H}^+ \quad \Delta G^\circ = -798\;\text{kJ/mol} $

This is sufficient to drive carbon fixation and support dense communities of specialised organisms in an otherwise food-poor environment (detailed in Module 5).

3.1b Marine Snow Sinking Dynamics

The sinking rate of marine snow particles determines how much organic matter reaches the bathypelagic zone before being consumed. Sinking velocity follows a modified Stokes law for porous aggregates:

$ w_s = \frac{2}{9} \frac{(\rho_p - \rho_w) g r^2}{\mu} \cdot (1 - \phi)^{n} $

where $\rho_p$ is the particle density, $\phi$ is the porosity (marine snow is 97–99% water by volume), $r$ is aggregate radius, and$n \approx 4.65$ is the Richardson–Zaki exponent. Typical sinking rates are 50–200 m/day for large aggregates ($>0.5\;\text{mm}$) and 10–50 m/day for smaller particles. At 100 m/day, it takes 10–40 days for a particle to reach the bathypelagic zone — during which time bacteria attached to the particle consume a significant fraction of its carbon content.

Faecal pellets from zooplankton sink fastest (100–1,000 m/day) because they are dense and compact. Salps produce particularly large, dense faecal pellets that can reach the deep ocean in just 1–2 days, bypassing mesopelagic remineralisation. This salp shunt can increase bathypelagic carbon flux by 2–5-fold during salp bloom events.

3.2 Deep-Sea Fish Adaptations

Deep-sea fishes display a suite of remarkable adaptations driven by three relentless selection pressures: darkness, food scarcity, and high pressure.

The Visual-Interaction Hypothesis & Reduced Metabolism

Childress (1995) proposed the visual-interaction hypothesis: in dark environments, organisms have fewer visually-mediated predator–prey encounters, reducing the selective pressure for high locomotory capacity. Consequently, bathypelagic fish have dramatically lower metabolic rates than their shallow-water counterparts.

The depth-dependent scaling of basal metabolic rate can be approximated as:

$ \text{BMR}_{\text{deep}} = \text{BMR}_{\text{surface}} \times \left(\frac{z}{z_{\text{ref}}}\right)^{-0.4} $

where $z_{\text{ref}} = 100\;\text{m}$ is a reference depth. At 2,000 m:

$ \frac{\text{BMR}_{2000}}{\text{BMR}_{\text{surface}}} = \left(\frac{2000}{100}\right)^{-0.4} = 20^{-0.4} \approx 0.30 $

Deep-sea fish at 2 km depth require only $\sim 30\%$ of the metabolic energy of a surface fish of equivalent mass. This is achieved through multiple morphological adaptations:

Watery muscles: reduced protein content (sometimes $<5\%$ wet weight vs $\sim 20\%$ in surface fish), replaced by water and gelatinous tissue
Reduced ossification: bones are poorly mineralised, reducing the metabolic cost of maintaining the skeleton
Large mouths and expandable stomachs: food encounters are rare, so the ability to swallow prey larger than oneself is essential
Slow growth rates: some species grow less than 1 mm per year, with lifespans exceeding 100 years

Bioluminescence as Communication

In the absence of sunlight, over 90% of bathypelagic organisms produce their own light via bioluminescence. The anglerfish (Melanocetus johnsonii) uses a bioluminescent lure (esca) containing symbiotic bacteria (Photobacterium) to attract prey. The light is produced by the luciferin–luciferase reaction:

$ \text{Luciferin} + \text{O}_2 \xrightarrow{\text{luciferase}} \text{Oxyluciferin} + h\nu $

The peak emission wavelength is $\lambda \approx 470\text{--}490\;\text{nm}$ (blue-green), matching the wavelength that travels farthest through seawater. The quantum yield of bacterial bioluminescence is remarkably high: $\Phi \approx 0.1\text{--}0.3$, meaning 10–30% of the chemical energy is converted to photons.

3.3 Pressure-Adapted Biochemistry

At 2,000 m depth, the hydrostatic pressure is approximately 200 atm ($\sim 20\;\text{MPa}$). This extreme pressure would denature proteins and collapse cell membranes of surface organisms. Deep-sea life has evolved elegant molecular solutions.

Pressure Effects on Protein Stability

The effect of pressure on a chemical equilibrium is governed by the volume change of reaction $\Delta V$. For a protein folding/unfolding equilibrium with equilibrium constant $K$:

$ \ln\!\left(\frac{K_P}{K_0}\right) = -\frac{\Delta V \cdot P}{RT} $

where $K_P$ is the equilibrium constant at pressure $P$,$K_0$ at atmospheric pressure, $R = 8.314\;\text{J/(mol\cdot K)}$, and $T$ is temperature. Protein unfolding typically has $\Delta V < 0$(the unfolded state is more compact due to hydration of exposed residues), so high pressure promotes unfolding.

For a typical deep-sea protein at 200 atm with $\Delta V = -20\;\text{mL/mol}$and $T = 275\;\text{K}$:

$ \ln\!\left(\frac{K_P}{K_0}\right) = -\frac{(-20 \times 10^{-6})(200 \times 101325)}{8.314 \times 275} = +0.177 $

So $K_P/K_0 = e^{0.177} \approx 1.19$ — a 19% increase in the unfolding equilibrium constant at 200 atm. Deep-sea proteins have evolved to counteract this through smaller cavities in their folded structures, reducing $|\Delta V|$.

TMAO: The Deep-Sea Osmolyte

Trimethylamine N-oxide (TMAO) is a small organic molecule that stabilises proteins against pressure-induced denaturation. Yancey et al. (2014) discovered a remarkable linear relationship between TMAO concentration and depth:

$ [\text{TMAO}] \approx 0.04 \times \text{depth(m)} \;\;\text{mmol/kg} $

At 2,000 m: $[\text{TMAO}] \approx 80\;\text{mmol/kg}$. At 4,000 m:$[\text{TMAO}] \approx 160\;\text{mmol/kg}$. TMAO counteracts pressure by preferentially hydrating the protein backbone, making the unfolded state thermodynamically unfavourable. The free energy of stabilisation is proportional to the osmolyte concentration:

$ \Delta G_{\text{TMAO}} = m \cdot [\text{TMAO}] $

where $m > 0$ is the transfer free energy coefficient (typically $\sim 40\text{--}100\;\text{J/(mol\cdot M)}$per backbone unit).

Membrane Fluidity Under Pressure

Cell membranes must remain in a liquid-crystalline state to function. High pressure promotes the gel (ordered) phase, reducing membrane fluidity. Deep-sea organisms compensate by increasing the proportion of unsaturated fatty acids in their membranes — each $\text{C=C}$ double bond introduces a kink that disrupts packing:

Surface fish: membrane lipids $\sim 30\%$ unsaturated
1,000 m depth: $\sim 50\%$ unsaturated
3,000 m depth: $\sim 70\%$ unsaturated, enriched in DHA (22:6$\omega$3) and EPA (20:5$\omega$3)

The phase transition temperature of a lipid bilayer shifts with pressure according to the Clausius–Clapeyron relation for phase transitions:

$ \frac{dT_m}{dP} = \frac{T_m \Delta V_m}{\Delta H_m} \approx 0.02\;\text{°C/atm} $

At 200 atm, the melting temperature shifts up by $\sim 4$°C, which would solidify membranes if not compensated by increased unsaturation.

3.4 Giant Squid: The Deep's Apex Predator

Architeuthis dux, the giant squid, inhabits the bathypelagic zone at depths of 600–1,200 m. Reaching total lengths of 13 m (females) and masses of 275 kg, it is one of the largest invertebrates on Earth.

The Largest Eyes in the Animal Kingdom

Giant squid possess eyes up to 27 cm in diameter — the size of a dinner plate and the largest of any living animal. Nilsson et al. (2012) showed these enormous eyes are specifically optimised for detecting the bioluminescent wakes of sperm whales approaching from distances of up to 120 m.

The photon detection sensitivity of a camera-type eye scales with the square of the pupil diameter$D$ and inversely with the solid angle subtended by each photoreceptor:

$ S = \frac{\pi}{4} D^2 \cdot \frac{\pi}{4} d^2 \cdot \eta \cdot \frac{1}{f^2} $

where $D$ is pupil diameter ($\approx 9\;\text{cm}$ in giant squid),$d$ is photoreceptor diameter ($\approx 5\;\mu\text{m}$),$\eta$ is quantum efficiency ($\approx 0.3$), and$f$ is focal length ($\approx 14\;\text{cm}$). The enormous pupil allows the eye to collect $\sim 100\times$ more photons than a typical fish eye ($D \approx 1\;\text{cm}$).

The detection range for a point source of bioluminescence ($I_0$ photons/s) follows from the inverse-square law and the minimum detectable photon flux:

$ r_{\max} = \sqrt{\frac{I_0}{4\pi \cdot N_{\min}}} \cdot e^{-\alpha r / 2} $

where $\alpha \approx 0.03\;\text{m}^{-1}$ is the attenuation coefficient of deep seawater at 480 nm, and $N_{\min}$ is the minimum detectable photon flux per unit area.

Ammonia Buoyancy

Unlike most squid, Architeuthis achieves near-neutral buoyancy not through a gas-filled chamber but by accumulating ammonium chloride (NH₄Cl) in its tissues and coelomic fluid. Ammonium ions are lighter than sodium ions:

$ \rho_{\text{NH}_4\text{Cl}} = 1.005\;\text{g/mL} \quad \text{vs} \quad \rho_{\text{NaCl}} = 1.025\;\text{g/mL} \quad \text{(at same molarity)} $

By replacing NaCl with NH₄Cl in tissue fluids, the squid reduces its overall density by approximately 2%. This is metabolically cheaper than actively swimming to maintain depth, though it makes the flesh taste strongly of ammonia — explaining why giant squid are not a commercial fishery target.

3.5 Whale Falls: Deep-Sea Oases

When a great whale dies and its carcass sinks to the bathypelagic or abyssal floor, it creates a whale fall — an isolated, nutrient-rich island in the otherwise barren deep-sea desert. A single 40-tonne whale delivers approximately$2 \times 10^9\;\text{kJ}$ of chemical energy to the seafloor, equivalent to 2,000 years of normal marine snow flux over the same area.

Smith & Baco (2003) described four successional stages, each dominated by different ecological processes:

Stage 1: Mobile Scavenger Stage (0–2 years)

Large scavengers arrive within hours: sleeper sharks (Somniosus pacificus), hagfish (Eptatretus), rattail fish (Macrouridae), and lysianassid amphipods. Consumption rates are extraordinary: a 35-tonne carcass can lose 40–60 kg of soft tissue per day. Hagfish can strip a whale to bone in 2–4 years.

Stage 2: Enrichment Opportunist Stage (2–4 years)

After the soft tissue is consumed, dense mats of polychaete worms (especially Ophryotrochaand Vigtorniella) colonise the bones and surrounding sediment. Densities can reach$>40{,}000$ individuals/m². The bone marrow and residual lipids provide rich organic substrates for enrichment fauna.

Stage 3: Sulfophilic Stage (4–50+ years)

Whale bones contain up to 60% lipid by dry weight. Anaerobic decomposition of these lipids by sulfate-reducing bacteria produces H₂S, which fuels chemoautotrophic bacteria (Beggiatoa mats, endosymbiont-bearing vesicomyid clams, bathymodiolin mussels). This stage is analogous to a miniature hydrothermal vent ecosystem and can persist for 50–100 years. The key reaction:

$ \text{Bone lipids} + \text{SO}_4^{2-} \xrightarrow{\text{SRB}} \text{H}_2\text{S} + \text{HCO}_3^- $

Stage 4: Reef Stage (50–100+ years)

The mineralised skeleton persists as a hard substrate — rare on the soft-sediment abyssal plain — attracting suspension feeders, sponges, and other organisms that require solid attachment points. The whale skeleton becomes a reef lasting centuries.

An estimated 690,000 great whale carcasses currently rest on the ocean floor at any given time, creating a network of deep-sea stepping stones separated by an average distance of $\sim 12\;\text{km}$ —within the dispersal range of many deep-sea larvae.

3.6 Deep-Sea Gigantism

Many bathypelagic organisms grow to sizes far exceeding their shallow-water relatives, a phenomenon known as deep-sea gigantism or abyssal gigantism. Examples include the giant isopod (Bathynomus giganteus, up to 50 cm vs 1–5 cm for coastal isopods), the giant amphipod (Alicella gigantea, 34 cm), and the giant squid itself.

Several hypotheses have been proposed for this phenomenon:

Bergmann's rule extended: lower temperatures favour larger body sizes to reduce surface-area-to-volume ratio and conserve metabolic heat. The relationship $M \propto T^{-\alpha}$ with $\alpha \approx 0.6\text{--}0.8$predicts significant size increases at bathypelagic temperatures.
Kleiber's law & starvation resistance: metabolic rate scales as $B = B_0 M^{0.75}$, while energy reserves scale as$E \propto M^{1.0}$. Therefore, starvation time scales as$\tau \propto E/B \propto M^{0.25}$ — larger animals survive longer between rare meals.
Reduced predation pressure: in the food-scarce bathypelagic, predator densities are low, relaxing size constraints imposed by predator gape-limitation.
Delayed maturation: slow growth rates but long lifespans allow organisms to reach very large sizes before reproduction. Some deep-sea species grow for decades before first reproducing.

The oxygen minimum zone (OMZ) at 500–1,000 m acts as a size bottleneck: organisms must be small to meet O₂ demands in low-oxygen water. Below the OMZ, as oxygen increases again in the cold, well-ventilated deep water, the constraints on body size are lifted.

3.7 Bathypelagic Microbial Ecology

While large organisms receive the most attention, the bathypelagic zone is dominated by microorganisms. Bacterial and archaeal cell densities are approximately $10^4\text{--}10^5$ cells per mL, compared to$10^6$ cells per mL in surface waters.

The turnover time of bathypelagic dissolved organic carbon (DOC) is remarkably slow:

$ \tau_{\text{DOC}} = \frac{[\text{DOC}]}{R_{\text{bacterial}}} \approx \frac{40\;\mu\text{mol C/L}}{0.01\;\mu\text{mol C L}^{-1}\text{d}^{-1}} \approx 4{,}000\;\text{days} \approx 11\;\text{years} $

This is because the DOC at these depths is highly refractory —it has been processed multiple times during its descent and consists mainly of humic substances and other complex molecules resistant to enzymatic degradation. The concept of$^{14}\text{C}$ dating reveals that bathypelagic DOC is, on average, 4,000–6,000 years old.

Bathypelagic archaea, particularly Thaumarchaeota (now Nitrososphaerota), play a critical role as chemoautotrophic ammonia oxidisers. They fix CO₂ using energy from ammonia oxidation:

$ \text{NH}_3 + 1.5\text{O}_2 \rightarrow \text{NO}_2^- + \text{H}^+ + \text{H}_2\text{O} \quad \Delta G^\circ = -275\;\text{kJ/mol} $

These archaea constitute up to 40% of all prokaryotic cells in the bathypelagic water column and are responsible for a significant fraction of deep-ocean carbon fixation — a process termed dark carbon fixation that was only recently recognised as quantitatively important (Herndl et al., 2005).

The Viral Shunt

Viral lysis is a major source of mortality for bathypelagic bacteria, killing an estimated 20–40% of prokaryotic standing stock per day. This viral shunt converts cellular biomass back into dissolved organic matter (DOM), short-circuiting the microbial food web:

$ \text{Bacterial cell} \xrightarrow{\text{viral lysis}} \text{DOM} + \text{viruses} \xrightarrow{\text{bacterial uptake}} \text{new cells} $

Viral abundances in the bathypelagic are $\sim 10^5\text{--}10^6$ per mL, with virus-to-bacteria ratios of 5–25. This viral recycling loop is estimated to release 0.37–0.63 Gt C per year in the deep ocean — a globally significant carbon flux that sustains heterotrophic microbial communities far from any photosynthetic source.

Bathypelagic Zone Cross-Section

An overview of the bathypelagic zone (1,000–4,000 m) showing representative organisms, marine snow, and key adaptations:

Simulation: Marine Snow & Deep-Sea Adaptations

Modelling the Martin curve POC flux attenuation, depth-dependent metabolic rate scaling, TMAO piezolyte concentration, and giant squid eye sensitivity:

Bathypelagic Zone: Carbon Flux, Metabolism, TMAO & Eye Sensitivity

Python

Four-panel model of key bathypelagic adaptations and processes

script.py150 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('#0a0a1a')
fig.suptitle('Bathypelagic Zone: Adaptations & Processes', color='white', fontsize=16, fontweight='bold', y=0.98)

# ===== Panel 1: Martin Curve - POC Flux =====
ax = axes[0, 0]
ax.set_facecolor('#0a0a2e')
z = np.linspace(100, 6000, 500)
b_values = [0.6, 0.858, 1.2]
colors_b = ['#44ddaa', '#4488ff', '#ff6644']
labels_b = ['b=0.60 (oligotrophic)', 'b=0.858 (global mean)', 'b=1.20 (productive)']

for b, c, lab in zip(b_values, colors_b, labels_b):
    flux_frac = (z / 100.0) ** (-b)
    ax.plot(flux_frac * 100, z, color=c, linewidth=2.5, label=lab)

# Bathypelagic zone shading
ax.axhspan(1000, 4000, alpha=0.15, color='#2244aa')
ax.text(55, 2500, 'Bathypelagic', color='#4488cc', fontsize=11, fontweight='bold')

ax.set_xlabel('POC Flux (% of export at 100 m)', color='white', fontsize=11)
ax.set_ylabel('Depth (m)', color='white', fontsize=11)
ax.set_title('Martin Curve: POC Flux Attenuation', color='white', fontsize=13, fontweight='bold')
ax.invert_yaxis()
ax.set_xlim(0, 100)
ax.tick_params(colors='white')
ax.grid(True, alpha=0.15, color='#4444aa')
ax.legend(facecolor='#1a1a3e', edgecolor='#4444aa', labelcolor='white', fontsize=9)

# ===== Panel 2: Metabolic Rate vs Depth =====
ax = axes[0, 1]
ax.set_facecolor('#0a0a2e')
z_met = np.linspace(100, 5000, 500)
z_ref = 100.0

# Visual-interaction hypothesis scaling
bmr_ratio = (z_met / z_ref) ** (-0.4)
ax.plot(bmr_ratio * 100, z_met, color='#ff8844', linewidth=2.5, label='BMR = BMR$_{surf}$ (z/100)$^{-0.4}$')

# Data points from Childress (1995) - representative
data_z = [200, 500, 800, 1200, 2000, 3000, 4000]
data_bmr = [(d/100)**(-0.4) * 100 for d in data_z]
noise = [2, -3, 4, -2, 1, -1, 0.5]
data_bmr_noisy = [b + n for b, n in zip(data_bmr, noise)]
ax.scatter(data_bmr_noisy, data_z, color='#ffcc44', s=50, zorder=5, edgecolors='white', linewidth=0.5, label='Observed (Childress 1995)')

ax.axhspan(1000, 4000, alpha=0.15, color='#2244aa')
ax.set_xlabel('Metabolic Rate (% of surface)', color='white', fontsize=11)
ax.set_ylabel('Depth (m)', color='white', fontsize=11)
ax.set_title('Metabolic Rate Decline with Depth', color='white', fontsize=13, fontweight='bold')
ax.invert_yaxis()
ax.tick_params(colors='white')
ax.grid(True, alpha=0.15, color='#4444aa')
ax.legend(facecolor='#1a1a3e', edgecolor='#4444aa', labelcolor='white', fontsize=9)

# ===== Panel 3: TMAO Concentration =====
ax = axes[1, 0]
ax.set_facecolor('#0a0a2e')
z_tmao = np.linspace(0, 8500, 500)
tmao = 0.04 * z_tmao  # mmol/kg

ax.plot(z_tmao, tmao, color='#44ccaa', linewidth=2.5, label='[TMAO] = 0.04 * depth')

# Literature data points (Yancey et al. 2014)
species_z = [100, 500, 1000, 1500, 2000, 2500, 3000, 4000, 5000, 6000, 7000, 8000]
species_tmao = [6, 22, 38, 55, 82, 105, 125, 158, 210, 245, 290, 325]
ax.scatter(species_z, species_tmao, color='#88eedd', s=40, zorder=5, edgecolors='white', linewidth=0.5, label='Measured in fish muscle')

# Bathypelagic zone
ax.axvspan(1000, 4000, alpha=0.15, color='#2244aa')
ax.text(2500, 300, 'Bathypelagic', color='#4488cc', fontsize=11, fontweight='bold', ha='center')

# TMAO physiological limit annotation
ax.axhline(y=300, color='#ff6644', linewidth=1, linestyle='--', alpha=0.6)
ax.text(500, 310, 'TMAO alone insufficient above ~8000 m', color='#ff8866', fontsize=9)

ax.set_xlabel('Depth (m)', color='white', fontsize=11)
ax.set_ylabel('[TMAO] (mmol/kg)', color='white', fontsize=11)
ax.set_title('TMAO Piezolyte Concentration vs Depth', color='white', fontsize=13, fontweight='bold')
ax.tick_params(colors='white')
ax.grid(True, alpha=0.15, color='#4444aa')
ax.legend(facecolor='#1a1a3e', edgecolor='#4444aa', labelcolor='white', fontsize=9)

# ===== Panel 4: Giant Squid Eye Sensitivity =====
ax = axes[1, 1]
ax.set_facecolor('#0a0a2e')

# Detection range vs pupil diameter for bioluminescent source
alpha_water = 0.03  # attenuation coefficient m^-1 at 480 nm
I0 = 1e10  # photons/s from bioluminescent disturbance
N_min_values = [1e4, 1e3, 1e2]  # min photon flux per m^2 (different sensitivities)

D_pupil = np.linspace(0.5, 12, 200)  # pupil diameter cm

for N_min, color, label in zip(N_min_values, ['#ff6644', '#ffaa44', '#44ccff'],
                                ['Low sensitivity', 'Medium sensitivity', 'High sensitivity (dark-adapted)']):
    # Simplified detection range (ignoring attenuation for analytic form)
    r_detect = np.sqrt(I0 * (D_pupil/100)**2 * 0.3 / (4 * N_min))  # metres
    # Apply attenuation correction (iterative approx)
    r_detect_atten = r_detect * np.exp(-alpha_water * r_detect / 4)  # rough correction
    r_detect_atten = np.clip(r_detect_atten, 0, 300)
    ax.plot(D_pupil, r_detect_atten, color=color, linewidth=2, label=label)

# Mark key species
ax.axvline(x=1.0, color='#666688', linewidth=1, linestyle=':', alpha=0.6)
ax.text(1.2, 10, 'Typical fish
(1 cm)', color='#8888aa', fontsize=8)
ax.axvline(x=9.0, color='#4488cc', linewidth=1.5, linestyle=':', alpha=0.8)
ax.text(6.5, 10, 'Giant squid
(9 cm pupil)', color='#4488cc', fontsize=9, fontweight='bold')
ax.axvline(x=11.0, color='#44cc88', linewidth=1, linestyle=':', alpha=0.6)
ax.text(9.5, 250, 'Colossal squid
(11 cm)', color='#44cc88', fontsize=8)

ax.set_xlabel('Pupil Diameter (cm)', color='white', fontsize=11)
ax.set_ylabel('Detection Range (m)', color='white', fontsize=11)
ax.set_title('Eye Size vs Bioluminescence Detection Range', color='white', fontsize=13, fontweight='bold')
ax.tick_params(colors='white')
ax.grid(True, alpha=0.15, color='#4444aa')
ax.legend(facecolor='#1a1a3e', edgecolor='#4444aa', labelcolor='white', fontsize=9)

plt.tight_layout(rect=[0, 0, 1, 0.96])
plt.savefig('output.png', dpi=150, bbox_inches='tight', facecolor='#0a0a1a')

# Print key results
print("=== Martin Curve POC Flux ===")
for depth in [500, 1000, 2000, 4000]:
    flux = (depth/100)**(-0.858) * 100
    print(f"  Depth {depth}m: {flux:.1f}% of surface export")

print("
=== Metabolic Rate Scaling ===")
for depth in [500, 1000, 2000, 4000]:
    bmr = (depth/100)**(-0.4) * 100
    print(f"  Depth {depth}m: {bmr:.1f}% of surface BMR")

print("
=== TMAO Concentrations ===")
for depth in [1000, 2000, 3000, 4000]:
    print(f"  Depth {depth}m: {0.04*depth:.0f} mmol/kg")

print("
=== Giant Squid Eye ===")
print(f"  Pupil area ratio (squid/fish): {(9/1)**2:.0f}x")
print(f"  Estimated detection range: ~120 m (Nilsson et al. 2012)")

Click Run to execute the Python code

Code will be executed with Python 3 on the server

References

Childress, J. J. (1995). Are there physiological and biochemical adaptations of metabolism in deep-sea animals? Trends in Ecology & Evolution, 10(1), 30–36.

Martin, J. H., Knauer, G. A., Karl, D. M., & Broenkow, W. W. (1987). VERTEX: carbon cycling in the northeast Pacific. Deep Sea Research Part A, 34(2), 267–285.

Nilsson, D.-E., Warrant, E. J., Johnsen, S., Hanlon, R., & Shashar, N. (2012). A unique advantage for giant eyes in giant squid. Current Biology, 22(8), 683–688.

Smith, C. R., & Baco, A. R. (2003). Ecology of whale falls at the deep-sea floor. Oceanography and Marine Biology: An Annual Review, 41, 311–354.

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. Proceedings of the National Academy of Sciences, 111(12), 4461–4465.

Seibel, B. A., & Drazen, J. C. (2007). The rate of metabolism in marine animals: environmental constraints, ecological demands and energetic opportunities. Philosophical Transactions of the Royal Society B, 362(1487), 2061–2078.

Somero, G. N. (1992). Adaptations to high hydrostatic pressure. Annual Review of Physiology, 54(1), 557–577.