What are the main sources of nutrients that cause eutrophication?

The primary nutrients driving eutrophication are nitrogen and phosphorus. Agricultural fertilizers, especially nitrates and phosphates, leach into streams and groundwater during rain events. Animal waste from livestock operations contributes large nitrogen loads. Untreated or partially treated sewage releases both nutrients. Atmospheric nitrogen deposition from vehicle and industrial emissions also contributes to coastal eutrophication. Urban lawns, golf courses, and stormwater runoff are additional point and non-point sources.

Why is phosphorus considered the main limiting nutrient in freshwater?

In most freshwater lakes, phosphorus is the nutrient in shortest supply relative to what algae need (Redfield ratio), making it the primary growth-limiting factor. Adding small amounts of phosphorus can dramatically increase algal production. Nitrogen is more often limiting in estuarine and coastal marine systems. This distinction matters for management: targeting phosphorus removal from wastewater is most effective for protecting lakes, while nitrogen reduction is critical for coastal dead zones.

What are hypoxic dead zones and how do they form?

Dead zones are bottom-water regions where dissolved oxygen drops below 2 mg/L—too low for most fish and invertebrates to survive. They form when large algal blooms die and sink to the bottom, where bacterial decomposition consumes oxygen faster than mixing can replenish it. Seasonal stratification (warm surface water floating on cold bottom water) prevents reoxygenation from the surface, allowing hypoxia to persist for months. When stratification breaks down in autumn, oxygen is restored but ecological damage has already occurred.

Home▸Nature & Climate▸Eutrophication

🦠 Eutrophication

Interactive lake eutrophication simulator: phosphorus loading, algae blooms, dissolved oxygen crash and fish kill. Explore regime shifts and hysteresis.

Nature & Climate3DModerate60 FPS💧 Water

About this simulation

This model runs a five-variable system of ordinary differential equations — phosphorus (P), algae biomass (A), dissolved oxygen (O), fish biomass (F) and dead organic matter (D) — integrated step by step with fourth-order Runge-Kutta. Algal uptake of phosphorus follows Monod kinetics, whilst oxygen is replenished by reaeration and drawn down by respiration and decomposition. Push phosphorus loading past a threshold and the lake tips from a clear, oxygenated state into a turbid, oxygen-starved one, and winding the loading back down does not retrace the same path, producing a hysteresis loop you can watch trace out on the phase plot.

🔬 What it shows

A lake tank sits alongside live bar readouts for phosphorus, algae, oxygen and fish, plus a loading-versus-algae phase plot. Four regimes — clear, transition, blooming and anoxic — are colour-coded and update in real time as the underlying equations are integrated.

🎮 How to use

Choose a starting preset (Oligotrophic, Mesotrophic, Eutrophic, Collapsing), then drag Phosphorus loading, Fish stocking (grazing) and Reaeration to see how each shifts the balance. Heavy rain and Drought inject short nutrient shocks, Sim speed rescales time, and Pause/Reset control the run.

💡 Did you know?

Real dead zones, such as the one in the Gulf of Mexico, can exceed ten thousand square kilometres in a single summer — driven by exactly the phosphorus loading this model condenses into a single slider.

Frequently asked questions

Why doesn't lowering phosphorus loading immediately clear the lake again?

Dead algae settle as organic matter, and once oxygen falls low enough that matter recycles phosphorus back out of the sediment faster than it settles, so the lake's internal phosphorus supply keeps feeding blooms even after external loading drops. This delayed feedback is what produces the hysteresis loop traced on the phase plot.

What determines when fish start dying in the model?

Fish biomass loses mass through a constant background mortality plus an oxygen-linked term that switches on once dissolved oxygen drops below 3 mg/L; the further oxygen falls beneath that threshold, the faster fish die off, which is why the regime label flips to anoxic once oxygen or fish biomass crosses its critical value.

What does the Reaeration slider actually change?

It scales how quickly oxygen exchanges with the atmosphere at the lake surface, pulling dissolved oxygen back towards its saturation value of 10 mg/L. Turning reaeration up can partly offset a deficit caused by algal respiration and decomposition, but cannot by itself stop a bloom from forming.

Why do Heavy rain and Drought behave differently from the loading slider?

Heavy rain adds a temporary pulse to phosphorus input that decays away over time, while Drought temporarily raises the concentration multiplier on existing loading. Both effects fade exponentially, modelling short-lived weather events rather than a permanent change in nutrient management.

What is the hysteresis area statistic tracking?

It is a running shoelace-formula estimate of the area enclosed by the loop traced in the loading-versus-algae phase plot. A larger area means the clear and turbid states are more strongly separated, so more effort is needed to push the lake back to clear once it has tipped into bloom.

⚙ Under the hood

Slowly raise phosphorus loading and a clear lake tips into algal blooms, anoxia and fish kill — and may not recover when loading is reduced. A regime-shift ecosystem with hysteresis.

Canvas 2DEcologyAlgal BloomAnoxiaRegime Shift

3D · Three.js / WebGL renderer · 60 FPS target · runs fully client-side, no install