Know more

Our use of cookies

Cookies are a set of data stored on a user’s device when the user browses a web site. The data is in a file containing an ID number, the name of the server which deposited it and, in some cases, an expiry date. We use cookies to record information about your visit, language of preference, and other parameters on the site in order to optimise your next visit and make the site even more useful to you.

To improve your experience, we use cookies to store certain browsing information and provide secure navigation, and to collect statistics with a view to improve the site’s features. For a complete list of the cookies we use, download “Ghostery”, a free plug-in for browsers which can detect, and, in some cases, block cookies.

Ghostery is available here for free:

You can also visit the CNIL web site for instructions on how to configure your browser to manage cookie storage on your device.

In the case of third-party advertising cookies, you can also visit the following site:, offered by digital advertising professionals within the European Digital Advertising Alliance (EDAA). From the site, you can deny or accept the cookies used by advertising professionals who are members.

It is also possible to block certain third-party cookies directly via publishers:

Cookie type

Means of blocking

Analytical and performance cookies

Google Analytics

Targeted advertising cookies


The following types of cookies may be used on our websites:

Mandatory cookies

Functional cookies

Social media and advertising cookies

These cookies are needed to ensure the proper functioning of the site and cannot be disabled. They help ensure a secure connection and the basic availability of our website.

These cookies allow us to analyse site use in order to measure and optimise performance. They allow us to store your sign-in information and display the different components of our website in a more coherent way.

These cookies are used by advertising agencies such as Google and by social media sites such as LinkedIn and Facebook. Among other things, they allow pages to be shared on social media, the posting of comments, and the publication (on our site or elsewhere) of ads that reflect your centres of interest.

Our EZPublish content management system (CMS) uses CAS and PHP session cookies and the New Relic cookie for monitoring purposes (IP, response times).

These cookies are deleted at the end of the browsing session (when you log off or close your browser window)

Our EZPublish content management system (CMS) uses the XiTi cookie to measure traffic. Our service provider is AT Internet. This company stores data (IPs, date and time of access, length of the visit and pages viewed) for six months.

Our EZPublish content management system (CMS) does not use this type of cookie.

For more information about the cookies we use, contact INRA’s Data Protection Officer by email at or by post at:

24, chemin de Borde Rouge –Auzeville – CS52627
31326 Castanet Tolosan CEDEX - France

Dernière mise à jour : Mai 2018

Menu Logo Principal ANR logo of NSFC mistea niglas leesu iEES LEMON LBE


State of Knowledge

Cyanobacteria have attracted the interest of the scientific community for a long time. Since the 1990s, mathematical models have been developed in order to better understand their dynamics and the formation of blooms [31], [32]. Cyanobacteria biomass dynamics was first represented by its only growth function, whose dependence with respect to light intensity, water temperature and nutrient availability has been intensively studied (see for example [33]), leading to a catalogue of mathematical expressions. Attention has then been focused on the representation of the vertical migration of the cyanobacteria in the water column, which is part of their strategy to capture light [34], [35]. But to reproduce the dynamics of cyanobacteria blooms, models had to take into account the interaction of cyanobacteria with the environment (especially hydrodynamics) and the other organisms of the ecosystem (bacteria, zooplankton, etc.).

Ecosystem models

According to [31] , we can distinguish two main categories of models of Harmful Algal Bloom (HAB), that we will call “minimal dynamic models” [36] and “complex dynamic models”.

“Minimal dynamic models” are “heuristic, examining the likelihood of certain processes generating a HAB” [31] . These models generally take into account between 2 and 5 variables like in the well-known Nutrient-Phytoplankton-Zooplankton (NPZ) models [37] . Several hypotheses on the origin of blooms have been studied with these models. In [38] , the relaxation of the zooplankton grazing pressure compared to the phytoplankton growth rate is identified as a potential key mechanism of bloom formation. Conversely, the control of the initiation and demise of the blooms by nutrients levels (bottom-up control) is numerically investigated in [39] . Most of these models are deterministic differential equations but the stochastic nature of the events and pathways leading to cyanobacterial blooms [40], [41] has more recently advocated for the development of stochastic models [42] . 

“Complex dynamic models” are designed to be used as a virtual reality [36] for the simulation of the whole ecosystem (e.g., CAEDYM, PCLake, DELWAQ). They are therefore generally coupled with a detailed hydrodynamic model and validated on field data. They are used to better understand the cyanobacterial blooms dynamics [43] , to support the lake management strategies by forecasting the blooms occurrence at short term [44] or to anticipate the future behaviour of the ecosystem in response to different local management or global changes scenarios [45] . The main approaches used for lake ecosystem modelling are reviewed in [36], [46], [47] . For short-term (day to week) forecast, genetic models or artificial neural networks [48] can be used. But for the simulation of annual cycles, such kind of models would not be able to infer the shift in the phytoplankton succession whereas deterministic models could [46] . Individual based models have also been considered, but as the number of individuals is limited by the computational cost, it is necessary to develop some derived approaches [36] .

“Complex dynamic models” generally include numerous biogeochemical variables and processes (biological and physical). Consequently, the outcomes of simulated scenarios are hampered with great uncertainty because of the huge set of parameters (sometimes more then 100) to be calibrated. Similarly with weather forecast approach, several “ensemble approaches” have been proposed to reduce this uncertainty. In [49] , [50] an ensemble of model runs obtained with several “equivalent” set of parameters are used to improve the predictive power of the model. In [45] , the same scenario is simulated using an ensemble of models. Another effective way to improve the performances of 3D hydrodynamic-ecological models is to conduct cross-comparisons of similar models on the same study sites. This type of comparative studies has been performed recently for 1D models [51]–[53] . To our knowledge very few attempts have been made for hydrodynamic-ecological models [45] and none for 3D models. To do so, portability of the model software is required. This type of concern is just starting to be addressed by the modeller community (e.g. [54], [55] ). Defining a statistical framework to compare the simulations would therefore be very helpful [56], [57] .

Modelling of specific processes

Despite the complexity of the models, some processes are still poorly represented and the models need therefore to be improved. For example, most of the models are not able to simulate the shifts in the phytoplankton community and the variation of the stoichiometry of living organisms [58] . Recently, it has been shown that the universally used Michaelis-Menten kinetics is not able to fully represent the variable stoichiometry of biomass and of nutrient uptake [59] . Promising size- and trait-based kinetics of nutrient uptake have been proposed but not yet tested against experimental or field data [60] . Some recent process-based approaches, which take into account the major eco-physiological traits of phytoplankton [61], [62] , have also recently been proposed in order to well represent the adaptive shifts of the phytoplankton community. In particular, light utilization efficiency which is a contrasted trait between cyanobacteria and other phytoplankton groups [63] can be used to improve phytoplankton succession modelling.

Regarding the cyanobacteria specific processes, modelling their migration in the water column remains an essential research issue. The key processes driving the sinking or floating of the cells are still poorly simulated [64] . Experimental [65], [66] and numerical studies [67] have also proved that the microbial loop can have a great impact on the cyanobacteria dynamics. Yet, the process of mineralization of the organic matter by bacteria is often represented as a simple chemical reaction with a rate either constant or dependent on the bacteria concentration at best [68], [69] . Physical processes, occurring in lakes at bi- or three-dimensional scales, such as hydrodynamics, sediment resuspension and internal waves, play also a key role (e.g. [70] ). Yet, simulations with 2D or 3D ecological models and their confrontation with consistent field data set are still very rare [71] . Until now, 3D models were principally used to assess the effect of environmental changes on lake hydrodynamics [72] .

Impacts of changes in C/N/P ratio and water temperature on cyanobacteria population dynamics

The increase in the occurrence of cyanobacterial blooms in freshwater ecosystems is due for a large part to an increase in nutrient (P and N) availability, these two elements being limiting for the primary producers of these ecosystems (e.g. [73] ). Rising CO2 concentrations and water temperature also enhance the growth rate of cyanobacteria and consequently the potential occurrence of blooms in eutrophic ecosystems (e.g. [24], [74] ). Other processes, also influenced by global changes, contribute indirectly to the increase of cyanobacterial blooms. For example, climate warming enhances the release of nutrients from the sediments and the mineralization of the organic matter conducting to an increase of P and N availability [75], [76] .

However, little is known about the consequences of the changes occurring in C/N/P ratio on cyanobacterial blooms while freshwater ecosystems are subject to different trajectories concerning this ratio. Indeed, as mentioned previously, P concentrations are decreasing in Europe but increasing in China and carbonate equilibrium is changing due to the increase in CO2 concentrations. We therefore expect changes in the ecological stoichiometry of freshwater ecosystems, whose consequences are largely unknown. It has been shown [77] that rising pCO2 will change elemental stoichiometry of phytoplankton with putative consequences on zooplankton feeding on phytoplankton. In the same way, a recent study conducted on 24 freshwater bacterial strains [78] has revealed that the stoichiometric flexibility of heterotrophic bacteria depends on the species and on their cellular phosphorus quota, suggesting that changes in C/N/P ratio can lead to changes in the composition of the bacterial communities. The same authors [79] also demonstrated that the aquatic bacterial communities display very contrasting values in their P content and stoichiometry. These authors proposed that anthropogenic changes in C/N/P ratio might have consequences on aquatic ecosystem productivity and on the extent and periodicity of nutrient fluctuations.

To the best of our knowledge, no experimental study has been conducted to test the combined impact of changes occurring in P, N and CO2 concentrations and in water temperature on the mineralization of the organic matter, the nutrient release from the sediment and their consequence on the blooms of cyanobacteria.