Phosphorus stock P with sigmoidal sediment recycling: dP/dt = L − sP + r·P^q/(P^q+m^q). Two stable states: clear (low P) and turbid (high P). Restoration shows hysteresis.