Simple mass balance model for lakes based on the Vollenweider equation:
dMw/dt = Min - sMw + pMs - Mout
The model was first used in the 1960s to determine the phosphorus concentration in lakes and reservoirs for eutrophication assessment.
This version considers mercury, and adds diagenesis, using an extra state variable (mercury in the sediment), and incorporates desorption processes that release mercury trapped in the sediment back to the water column.
The temporal dynamics of the model simulate the typical development of pollution in time.
1. Low loading, low Hg concentration in lake
2. High loading, increasing Hg concentration in lake
3. Desorption rate is low, Hg in sediment increases
4. Measures implemented for source control, loading reduces
5. Hg in lake gradually decreases, but below a certain point, desorption increases, and lake Hg concentration does not improve
6. Recovery only occurs when the secondary load in the sediment is strongly reduced.