Mathematical modeling of phytoplankton populations evolution in the Azov Sea

Abstract

The paper describes the construction of a three-dimensional mathematical model of biogeochemical processes, considering the salinity and temperature influence on the phytoplankton populations development. The paper proposes a new difference scheme for solving convection-diffusion-reaction problems at large values of the Peclet grid number (2<Pe≤20), which is a linear combination of the central and the upwind leapfrog difference schemes. The three-layer difference scheme is more accurate than the traditional upwind leapfrog difference scheme for problems where convection prevails over diffusion. The construction of discrete equations for solving the problem of biogeochemical cycles on the basis of the scheme considering the filling of cells is described. The stationary regimes of phytoplankton dynamics problem were researched considering the transformation of phosphorus, nitrogen and silicon forms. Results of software complex, which allows to simulate biogeochemical processes in the Azov Sea, were described. The software package allows to forecast the dynamics of the Azov Sea ecosystem development in the conditions of modern salinization.