A dynamic mass-balance model for phosphorus in lakes with a focus on criteria for applicability and boundary conditions

This paper presents an improved version of a general, process-based mass-balance model (LakeMab/LEEDS) for phosphorus in entire lakes (the ecosystem scale). The focus in this work is set on the boundary conditions, i.e., the domain of the model, and critical tests to reveal those boundary conditions using data from a wide limnological range. The basic structure of the model, and many key equations have been presented and motivated before, but this work presents several new developments. The LakeMab-model is based on ordinary differential equations regulating inflow, outflow and internal fluxes and the temporal resolution is one month to reflect seasonal variations. The model consists of four compartments: surface water, deep water, sediment on accumulation areas and sediment on areas of erosion and transportation. The separation between the surface-water layer and the deep-water layer is not done from water temperature data, but from sedimentological criteria (from the theoretical wave base, which regulates where wind/wave-induced resuspension of fine sediments occurs). There are algorithms for processes regulating internal fluxes and internal loading, e.g., sedimentation, resuspension, diffusion, mixing and burial. Critical model tests were made using data from 41 lakes of very different character and the results show that the model could predict mean monthly TP-concentrations in water very well (generally within the uncertainty bands given by the empirical data). The model is even easier to apply than the well-known OECD and Vollenweider models due to more easily accessed driving variables