# The community-centered freshwater biogeochemistry model unified RIVE v1.0: a unified version for water column

Shuaitao Wang, Vincent Thieu, Gilles Billen, Josette Garnier, Marie Silvestre, Audrey Marescaux, Xingcheng Yan, Nicolas Flipo

► To cite this version:

Shuaitao Wang, Vincent Thieu, Gilles Billen, Josette Garnier, Marie Silvestre, et al.. The community-centered freshwater biogeochemistry model unified RIVE v1.0: a unified version for water column. Geoscientific Model Development, 2024, 17 (1), pp.449-476. 10.5194/gmd-17-449-2024. hal-04456280

HAL Id: hal-04456280

https://hal.science/hal-04456280v1

Submitted on 14 Feb 2024

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
<figure><img src="image_2.png" /></figure>

# The community-centered freshwater biogeochemistry model unified RIVE v1.0: a unified version for water column

Shuaitao Wang<fnref n="1" />, Vincent Thieu<fnref n="1" />, Gilles Billen<fnref n="1" />, Josette Garnier<fnref n="1" />, Marie Silvestre<fnref n="2" />, Audrey Marescaux<fnref n="1" />, Xingcheng Yan<fnref n="1" />, and Nicolas Flipo<fnref n="3" />

**Correspondence:** Shuaitao Wang (shuaitao.wang@sorbonne-universite.fr, shuaitaowang@outlook.com)

Received: 27 June 2023 – Discussion started: 15 August 2023

Revised: 7 November 2023 – Accepted: 9 November 2023 – Published: 16 January 2024

**Abstract.** Research on mechanisms of organic matter degradation, bacterial activities, phytoplankton dynamics, and other processes has led to the development of numerous sophisticated water quality models. The earliest model, dating back to 1925, was based on first-order kinetics for organic matter degradation. The community-centered freshwater biogeochemistry model RIVE was initially developed in 1994 and has subsequently been integrated into several software programs such as Seneque-Riverstrahler, pyNuts-Riverstrahler, PROSE/PROSE-PA, and Barman. After 30 years of research, the use of different programming languages including QBasic, Visual Basic, Fortran, ANSI C, and Python, as well as parallel evolution and the addition of new formalisms, raises questions about their comparability.

This paper presents a unified version of the RIVE model for the water column, including formalisms for bacterial communities (heterotrophic and nitrifying), primary producers, zooplankton, nutrients, inorganic carbon, and dissolved oxygen cycles. The unified RIVE model is open-source and implemented in Python 3 to create pyRIVE 1.0 and in ANSI C to create C-RIVE 0.32. The organic matter degradation module is validated by simulating batch experiments. The comparability of the pyRIVE 1.0 and C-RIVE 0.32 software is verified by modeling a river stretch case study. The case study considers the full biogeochemical cycles (microorganisms, nutrients, carbon, and oxygen) in the water column, as well as the effects of light and water temperature. The results show that the simulated concentrations of all state variables, including microorganisms and chemical species, are very similar for pyRIVE 1.0 and C-RIVE 0.32. This open-source

project highly encourages contributions from the freshwater biogeochemistry community to further advance the project and achieve common objectives.

## 1 Introduction

Modeling the water quality of a freshwater system (river, lake, or reservoir) is critical to understanding and managing its functioning. The functioning of a freshwater system is the result of complex interrelated biogeochemical processes. The first water quality model developed by Streeter and Phelps (1925) describes the degradation of organic matter (OM) in a river. The organic matter, measured globally by biochemical oxygen demand in 5 d (BOD5), is considered to be degraded according to first-order kinetics. Although dating back more than a century (the study was completed in 1915, but publication was delayed to 1925 due to World War I; Hellweger, 2015), this model is still widely used to represent the dynamics of organic matter in water quality modeling (Hellweger, 2015).

While the role of microorganisms in the degradation of organic matter has been acknowledged since the end of the 19th century, there is an important limitation of this type of representation. The microbiological nature of the organic matter degradation process and the bacterial population dynamics intrinsically involved are completely obscured. They are implicitly taken into account only through a biodegradability constant of organic matter and its dependence on temperature. Microbial biogeochemical work in the 1980s–1990s led

<footnote n="1">¹CNRS, EPHE, UMR Metis, Sorbonne Université, 75005 Paris, France</footnote>

<footnote n="2">²CNRS, FR3020 FIRE, Sorbonne Université, 75005 Paris, France</footnote>

<footnote n="3">³Mines Paris, PSL University, Center for geosciences and geoengineering, 77300 Fontainebleau, France</footnote>
to the elucidation of the detailed mechanisms of the organic matter degradation process and the associated heterotrophic bacterial activities (Fuhrman and Azam, 1982; Azam et al., 1983; Somville and Billen, 1983; Servais et al., 1985; Rego et al., 1985; Fontigny et al., 1987; Servais et al., 1987; Billen et al., 1988; Servais et al., 1989; Billen et al., 1990; Garnier et al., 1992a, b). This new corpus of knowledge led to the development and the formulation of the biogeochemical model RIVE (Billen et al., 1994; Garnier and Billen, 1994). It is capable of simulating the degradation of OM in freshwater systems and the associated oxygen consumption by bacterial activities, which is more realistic than the model of Streeter and Phelps (1925). In the RIVE model, the HSB model (Billen and Servais, 1989; Billen, 1991) is used to represent the degradation of organic matter and heterotrophic bacterial activities. This model simulates the exoenzymatic hydrolysis of particulate and dissolved organic matter (split into biodegradable and refractory pools), including high-weight polymers, into small monomeric substrates. These substrates are subsequently assimilated by bacteria for their growth and respiration.

Apart from the degradation of organic matter, the AQUA-PHY model (Lancelot et al., 1991) is used to simulate the dynamics of phytoplankton in the RIVE model (Billen et al., 1994). The model explicitly simulates photosynthesis of phytoplankton, growth, mortality, and respiration processes. In addition to water temperature, the photosynthesis depends on the light intensity, while growth is controlled by nutrient availability and small organic metabolites. The small organic metabolites are formed either directly by photosynthesis or by catabolysis of reserve products. This conceptualization allows for the growth of phytoplankton during dark periods. In addition, the model also introduces a limiting factor of nutrients in the growth of phytoplankton and considers the cycling of nutrients during the life cycle of phytoplankton.

Since its initial development by Billen et al. (1994), the RIVE model has co-existed within several software packages (Table A1) developed for different aquatic compartments and supported by the PIREN-Seine program (https://www.piren-seine.fr/, last access: 20 December 2023). The RIVE model was firstly applied in river systems using the Riverstrahler drainage network approach (Billen et al., 1994; Garnier et al., 1995). It was initially coded in QBasic and later on piloted by a GIS graphical interface, Seneque-Riverstrahler (Visual Basic; Ruelland et al., 2007). And it is now fully integrated within the pyNuts-Riverstrahler (https://gitlab.in2p3.fr/rive/pynuts/, last access: 20 December 2023). Python framework (Thieu et al., 2017). It can model the biogeochemical functioning of hydrographic networks at scales ranging from local to continental. The RIVE model was also applied to lentic freshwater systems like regulated reservoirs (Barman software, Garnier et al., 2000; Thieu et al., 2006; Yan et al., 2022a – Table A1) or to simulate hydrobiodynamic functioning of highly human-impacted river systems (PROSE software – Even et al., 1998, 2004, 2007;

Flipo et al., 2004; Vilmin et al., 2015b; PROSE-PA software – Wang et al., 2019, 2023a, https://gitlab.com/prose-pa/prose-pa last access: 20 December 2023; developed in ANSI C coupled with a self-developed lex and yacc parser; Table A1). The RIVE model is also coupled with the Soil & Water Assessment Tool (SWAT) to simulate the water quality of the Vienne basin, France (Manteaux et al., 2023), and incorporated into the QUAL-NET model (Minaudo et al., 2018) to simulate river eutrophication in the drainage network of the middle Loire River corridor, France. Moreover, the RIVE model is implemented into the VEMALA V3 model to simulate phosphorus and nitrogen loading in Finnish watersheds (Korppoo et al., 2017).

Based on the above implementations, different versions of the RIVE model code have successfully simulated a large variety of freshwater systems (lake or reservoirs, river systems) across the world. The parameter values were determined through laboratory experiments or calibrated with observation data (Garnier et al., 1992a; Servais and Garnier, 1993; Garnier and Billen, 1994; Billen et al., 1994; Garnier et al., 1995). These applications (Table A1) were carried out for different networks and scales as well as various degrees of anthropogenic impacts in a wide climatic gradient using either Riverstrahler (possibly with its Seneque or pyNuts environments) or PROSE/PROSE-PA, such as the Seine River (France) (Billen et al., 1994; Garnier et al., 1995; Even et al., 1998, 2004, 2007; Billen et al., 2007; Servais et al., 2007; Thieu et al., 2009, 2010; Vilmin et al., 2015b, a; Aissa-Grouz et al., 2016; Vilmin et al., 2016; Desmit et al., 2018; Vilmin et al., 2018; Romero et al., 2019; Marescaux et al., 2020; Wang et al., 2022), the Danube River (Romania and Bulgaria) (Garnier et al., 2002), the Red River (China and Vietnam) (Le et al., 2010, 2015; Nguyen et al., 2016) and its distributary the Day-Nhue River (Luu et al., 2021), the Lule and Kalix rivers (Sweden) (Sferratore et al., 2008), the Scheldt river (Belgium and Netherlands) (Billen et al., 2005; Thieu et al., 2009), the Zenne River (Belgium) (Garnier et al., 2013), the Mosel River (Germany) (Garnier et al., 1999a), the Somme River (France) (Thieu et al., 2009, 2010), the Loire River (France) (Garnier et al., 2018a), the Lot River (France) (Garnier et al., 2018b), and the Orgeval watershed (France) (Flipo et al., 2004, 2007; Garnier et al., 2014). Moreover, the RIVE model has also been applied to stagnant systems (e.g., sand-pit lake – Lake Crétail, France, Garnier and Billen, 1994; reservoirs – Marne, Aube, Seine, France, Garnier et al., 2000; Yan et al., 2022a).

After 30 years of research, the parallel evolutions of these codes, the numerical adaptations inherent in programming languages (QBasic, Visual Basic, Fortran, Python, and ANSI C), and the addition of new formalisms raise the question of their comparability. The identification of a unified version of the RIVE model is then necessary. A project aiming at unifying these RIVE implementations was undertaken. The unified version brings together all recent developments, especially the ones achieved with Python and ANSI C pro-
gramming languages. This action will strengthen the collaboration of the research teams involved in the development of the model. This paper presents a unified version of RIVE for the water column (called unified RIVE v1.0) with a presentation of the formalisms for the biogeochemical cycles. That integrates the bacterial communities (heterotrophic and nitrifying), primary producers, zooplankton, and fate of detritic organic matter either particulate or dissolved and as biodegradable and refractory, as well as the associated nutrients and dissolved oxygen cycles. The most recent developments in modeling inorganic forms of carbon are also presented. The unified RIVE v1.0 included in pyRIVE 1.0 (tested with Python 3 versions up to 3.10 release) and C-RIVE 0.32 is open-source and therefore available to the scientific community. A numerical experiment is then introduced to evaluate the comparability of the pyRIVE 1.0 and C-RIVE 0.32 through a systematic comparison of simulations produced under controlled conditions. We thus establish a reference framework to evaluate different implementations (programming languages, performance – comparability) of the unified RIVE v1.0 formulation that continues to evolve in several water quality models.

## 2 Model description

The unified RIVE v1.0 model simulates the cycling of carbon, nutrients, and oxygen within a freshwater system (river, lake, reservoir). Biogeochemical cycles are simulated with a community-centered or agent-based model. That means that the freshwater system functioning is explicitly modeled, taking into consideration the activities of microorganisms such as phytoplankton, zooplankton, heterotrophic bacteria, and nitrifying bacteria. Additionally, it accounts for physical processes like oxygen re-aeration and dilution. This modeling approach is developed in relation to water temperature, macronutrients, and organic matter, (particulate, dissolved, and biodegradable fractions). The organic matter degradation, nitrifying bacteria dynamics, primary producer dynamics, zooplankton dynamics, nutrients, and inorganic carbon cycling are described subsequently. A high number of model parameters are used to characterize the microorganisms’ properties and most of them have been determined through field or laboratory experiments under controlled conditions. This paper presents a focus on the conceptualization of the unified RIVE v1.0 model in the water column exclusively. While the RIVE model does have applications for sediment dynamics and its interaction with the water column (Even et al., 2004; Thouvenot et al., 2007; Billen et al., 2015; Vilmin et al., 2015a, 2016; Yan et al., 2022b), relevant community-centered efforts need to be made in future work, which is not the focus of this study.

### 2.1 Organic matter degradation

The mechanisms of organic matter degradation by the activity of heterotrophic bacteria are represented using the HSB model (Billen and Servais, 1989; Billen, 1991). It contains three variables: H represents high-weight polymers (large molecules), which form the majority of dissolved and particulate organic matter but must be exoenzymatically hydrolyzed to be accessible to heterotrophic bacteria; S represents small monomeric substrates (SMSs), directly accessible to microbial uptake; and B represents heterotrophic bacteria that absorb the substrates for their growth and respiration (Fig. 1). In diagrams of this paper (for instance the HSB model, Fig. 1), the state variables are shown as circles and represent either concentrations or stocks entering and leaving the (biogeochemical) processes. The biogeochemical processes are represented by squares.

The high-weight polymer (total organic carbon) is conceptually divided for each phase (dissolved – HD; particulate – HP) into three pools. Each pool is characterized by a specific biodegradability: (1) rapidly biodegradable in 5 d (HD₁ and HP₁), (2) slowly biodegradable in 45 d (HD₂ and HP₂), and (3) refractory (HD₃ and HP₃).

#### 2.1.1 Heterotrophic bacteria dynamics

The dynamic of heterotrophic bacteria is explicitly simulated, including growth, mortality, and respiration. The growth of heterotrophic bacteria depends on water temperature and the availability of small monomeric substrate (SMS). The dependence is represented by the Monod equation (Monod, 1949). A maximal substrate uptake rate at 20°C ($b_{max20,hb}$) and a bacterial growth yield ($Y_{hb}$) are used to calculate the growth rate of heterotrophic bacteria ($\mu_{hb_i}$) (Eq. 3). The fraction of uptake not used for growth ($1 - Y_{hb}$) is respired.

$$b_{hb_i} = b_{max20,hb_i} f(T)_{hb_i} \frac{[SMS]}{[SMS] + K_{sms,hb_i}} \quad (1)$$

$$f(T)_{hb_i} = \frac{e^{-\frac{(T-T_{opt,hb_i})^2}{\sigma_{hb_i}^2}}}{-\frac{(20-T_{opt,hb_i})^2}{\sigma_{hb_i}^2}} \quad (2)$$

$$\mu_{hb_i} = Y_{hb_i} b_{hb_i} \quad (3)$$

Here, $b_{hb_i}$ is the effective substrate uptake rate of the ith species of heterotrophic bacteria [h⁻¹], $b_{max20,hb_i}$ is the maximal substrate uptake rate of the ith species of heterotrophic bacteria at 20°C [h⁻¹], [SMS] is the small monomeric substrate concentration [mgCL⁻¹], $K_{sms,hb_i}$ is the half-saturation constant for small monomeric substrate of the ith species of heterotrophic bacteria [mgCL⁻¹], $f(T)_{hb_i}$ is the water temperature weight of the ith species of heterotrophic bacteria at T °C [–], $T_{opt,hb_i}$ is the optimal temperature of
<figure><img src="image_3.png" /><figcaption>Figure 1. Flowchart of the HSB model. HD: dissolved high-weight polymer; HP: particulate high-weight polymer; SMS: small monomeric substrate; HB: heterotrophic bacteria; PHY: phytoplankton; ZOO: zooplankton; nitr. bact.: nitrifying bacteria; extr. excretion of phytoplankton; sink.: sinking. respi.: respiration; $\epsilon_{hp1,2,3}$ and $\epsilon_{hd1,2,3}$: proportion to convert dead biomass to HP and HD.</figcaption></figure>

the ith species of heterotrophic bacteria for its growth [°C], $\sigma_{hb_i}$ is the range of temperature for the ith species of heterotrophic bacteria [°C], $Y_{hb_i}$ is the bacterial growth yield of the ith species of heterotrophic bacteria [–], and $\mu_{hb_i}$ is the effective growth rate of the ith species of heterotrophic bacteria [h⁻¹]

A sinking velocity ($v_{sb}$) is associated with each particulate species to represent particulate sinking by gravity. The mortality of heterotrophic bacteria is simulated by first-order kinetics (Eq. 4). The dead biomass of living species is converted into varying types of organic matter content, including both dissolved and particulate forms, based on specified proportions ($\epsilon_{hd}$ and $\epsilon_{hp}$, Fig. 1).

$$ \frac{d[HB_i]}{dt} = (\mu_{hb_i} - k_{d20,hb_i} f(T)_{hb_i} - k_{sink,hb_i})[HB_i] \\ k_{sink,hb_i} = \frac{v_{sbhb_i}}{\text{depth}} \quad (4) $$

Here, $\mu_{hb_i}$ is the effective growth rate of the ith species of heterotrophic bacteria [h⁻¹], $k_{d20,hb_i}$ is the mortality rate of the ith species of heterotrophic bacteria at 20 °C [h⁻¹], $k_{sink,hb_i}$ is the sinking rate of the ith species of heterotrophic bacteria [h⁻¹], $[HB_i]$ is the biomass concentration of the ith species of heterotrophic bacteria [mgC L⁻¹], $f(T)_{hb_i}$ is the water temperature weight at T °C defined by Eq. (2) [–], $v_{sbhb_i}$ is the sinking velocity of the ith species of heterotrophic bacteria [m h⁻¹], and “depth” represents the water depth [m].

### 2.1.2 Hydrolysis of high-weight polymer

The particulate biodegradable high-weight polymer ($HP_1$ and $HP_2$) is firstly hydrolyzed to the dissolved biodegradable high-weight polymer ($HD_1$ and $HD_2$). The dissolved biodegradable high-weight polymer is then hydrolyzed enzymatically to small monomeric substrate (Fig. 1). The hydrolysis of HP is represented by first-order kinetics (Eq. 5), while a Michaelis–Menten function (Michaelis and Menten, 1913) is used to express the exoenzymatic hydrolysis of HD (Eq. 6).

$$ \frac{d[HP_i]}{dt} = -k_{hp_i} \times [HP_i] + \left( \sum_j k_{d20,j} f(T)_j [LS]_j \right) \\ \epsilon_{hp_i} - k_{sink,hp_i} [HP_i] \quad (5) $$

Here, $[HP_i]$ is the concentration of particulate high-weight polymer $i \in \{1, 2\}$ [mgC L⁻¹], $k_{hp_i}$ is the hydrolysis rate of $HP_i$ ($i \in \{1, 2\}$) [h⁻¹], $f(T)_j$ is the water temperature weight of the jth living species at T °C defined as in Eq. 2 [–], $k_{d20,j}$ is the mortality rate of the jth living species (such as phytoplankton, zooplankton, bacteria) at 20 °C [h⁻¹], $[LS]_j$ is the concentration of the jth living species [mgC L⁻¹], $\epsilon_{hp_j}$ is the proportion to convert the dead biomass to $HP_j$ ($i \in \{1, 2\}$) [–],
<figure><img src="image_4.png" /><figcaption>**Figure 2.** Nitrifying bacteria dynamics. AOB: ammonia-oxidizing bacteria; NOB: nitrite-oxidizing bacteria; mort.: mortality; sink.: sinking.</figcaption></figure>

and $k_{\text{sink, hp}_i}$ is the sinking rate for HP$_i$ ($i \in \{1, 2\}$) [h$^{-1}$].

$$ \frac{\mathrm{d}[\mathrm{HD}_i]}{\mathrm{d}t} = -\sum_k \left( e_{\mathrm{max20,hd}_i, \mathrm{hb}_k} f(T)_{\mathrm{hb}_k} \frac{[\mathrm{HD}_i]}{[\mathrm{HD}_i] + K_{\mathrm{hd}_i, \mathrm{hb}_k}} [\mathrm{HB}_k] + k_{\mathrm{hp}_i} \times [\mathrm{HP}_i] + \left( \sum_j k_{\mathrm{d20}, j} f(T)_j [\mathrm{LS}]_j \right) \epsilon_{\mathrm{hd}_i} \right) \quad (6) $$

Here, $[\mathrm{HD}_i]$ is the concentration of dissolved high-weight polymer ($i \in \{1, 2\}$) [mgC L$^{-1}$], $e_{\mathrm{max20,hd}_i, \mathrm{hb}_k}$ is the maximum hydrolysis rate of HD$_i$ at 20°C related to HB$_k$ and $i \in \{1, 2\}$ [h$^{-1}$], $f(T)_{\mathrm{hb}_k}$ is the water temperature weight of the $k$th species of heterotrophic bacteria at T °C (Eq. 2) [–], $K_{\mathrm{hd}_i, \mathrm{hb}_k}$ is the half-saturation constant for HD$_i$ related to HB$_k$ and $i \in \{1, 2\}$ [mgC L$^{-1}$], $[\mathrm{HB}_k]$ is the concentration of the $k$th species of heterotrophic bacteria [mgC L$^{-1}$], $k_{\mathrm{d20}, j}$ is the mortality rate of the $j$th living species (such as phytoplankton, zooplankton, bacteria) at 20°C [h$^{-1}$], $f(T)_j$ is the water temperature weight of the $j$th living species at T °C defined as in Eq. (2) [–], $[\mathrm{LS}]_j$ is the concentration of the $j$th living species [mgC L$^{-1}$], and $\epsilon_{\mathrm{hd}_i}$ is the proportion to convert the dead biomass to HD$_i$ and $i \in \{1, 2\}$ [–].

## 2.2 Nitrifying bacteria dynamics

The unified RIVE v1.0 model includes the description of the nitrification microbial process, mediated by two types of nitrifying bacteria. They are respectively responsible for the production of nitrite ($NH_4^+$ + $\frac{3}{2}O_2 \rightarrow NO_2^-$ + 2H$^+$ + H$_2O$) and nitrate ($NO_2^-$ + $\frac{1}{2}O_2 \rightarrow NO_3^-$). The nitrifying bacteria get energy by oxidizing $NH_4^+$ (ammonium) and $NO_2^-$ (nitrite) for their growth. These two bacteria are named AOB

(ammonia-oxidizing bacteria) and NOB (nitrite-oxidizing bacteria) (Brion and Billen, 1998). The growth of nitrifying bacteria is limited by the availability of ammonium, nitrite, and oxygen, which is represented with Monod functions (Eq. 7). The effect of water temperature is also taken into account.

$$ \mu_{\mathrm{aob}} = \mu_{\mathrm{max20,aob}} f(T)_{\mathrm{aob}} \left( \frac{[\mathrm{NH}_4^+]}{[\mathrm{NH}_4^+] + K_{\mathrm{nh4,aob}}} \right) \left( \frac{[\mathrm{O}_2]}{[\mathrm{O}_2] + K_{\mathrm{o2,aob}}} \right) \quad (7) $$

$$ \mu_{\mathrm{nob}} = \mu_{\mathrm{max20,nob}} f(T)_{\mathrm{nob}} \left( \frac{[\mathrm{NO}_2^-]}{[\mathrm{NO}_2^-] + K_{\mathrm{no2,nob}}} \right) \left( \frac{[\mathrm{O}_2]}{[\mathrm{O}_2] + K_{\mathrm{o2,nob}}} \right) \quad (8) $$

Here, $\mu_{\mathrm{aob}}$ and $\mu_{\mathrm{nob}}$ are the effective growth rates of AOB and NOB [h$^{-1}$], $\mu_{\mathrm{max20,aob}}$ and $\mu_{\mathrm{max20,nob}}$ are the maximal growth rates of AOB and NOB at 20°C, respectively [h$^{-1}$], $f(T)_{\mathrm{aob}}$ and $f(T)_{\mathrm{nob}}$ are the water temperature weight at T °C defined as in Eq. (2) [–], $K_{\mathrm{nh4,aob}}$ and $K_{\mathrm{no2,nob}}$ are the half-saturation constants for $NH_4^+$ (AOB) and for $NO_2^-$ (NOB) [mgN L$^{-1}$], and $K_{\mathrm{o2,aob}}$ and $K_{\mathrm{o2,nob}}$ are the half-saturation constants for oxygen (AOB and NOB) [mgO$_2$ L$^{-1}$].

The mortality and sinking of nitrifying bacteria are simulated the same way as for other living species.

$$ \frac{\mathrm{d}[\mathrm{AOB}]}{\mathrm{d}t} = (\mu_{\mathrm{aob}} - k_{\mathrm{d20,aob}} f(T)_{\mathrm{aob}} - k_{\mathrm{sink,aob}})[\mathrm{AOB}] \quad (9) $$

$$ \frac{\mathrm{d}[\mathrm{NOB}]}{\mathrm{d}t} = (\mu_{\mathrm{nob}} - k_{\mathrm{d20,nob}} f(T)_{\mathrm{nob}} - k_{\mathrm{sink,nob}})[\mathrm{NOB}] \quad (10) $$

Here, $\mu_{\mathrm{aob}}$ and $\mu_{\mathrm{nob}}$ are the effective growth rates of AOB and NOB defined by Eqs. (7) and (8) [h$^{-1}$], $k_{\mathrm{d20,aob}}$ and $k_{\mathrm{d20,nob}}$ are the mortality rates of AOB and NOB at 20°C [h$^{-1}$], $k_{\mathrm{sink,aob}}$ and $k_{\mathrm{sink,nob}}$ are the sinking rates of AOB and NOB [h$^{-1}$], $f(T)_{\mathrm{aob}}$ and $f(T)_{\mathrm{nob}}$ are the water temperature weights at T °C defined as in Eq. (2) [–], and [AOB] and [NOB] are the concentrations of AOB and NOB [mgC L$^{-1}$].

## 2.3 Primary producer dynamics

The behavior of primary producers is represented using the AQUAPHY model (Lancelot et al., 1991). Biomass of a phytoplankton species is composed of three different cellular constituents (Fig. 3):

a. the structural and functional macromolecules of the cell, F, mainly proteins, chlorophyll, and structural lipids (such as membranes);

b. polysaccharides playing the role of reserve products, R;

c. monomeric (amino acids) and oligomeric precursors for macromolecular synthesis, S.
<figure><img src="image_5.png" /><figcaption>**Figure 3.** Description of the AQUAPHY model. F: functional macromolecules of the cell; R: reserve products; S: monomeric (amino acids) and oligomeric precursors for macromolecular synthesis. Phytoplankton biomass equals the sum of the three cellular constituents (F, R, S). SMS: small monomeric substrate. photos.: photosynthesis; respi.: respiration; excr.: excretion; synth.: synthesis; catab.: catabolysis; sink.: sinking</figcaption></figure>

At any time, the biomass of the *j*th phytoplankton species (mgC L⁻<fnref n="1" />), [PHY]j, is equal to the sum of the three internal constituents (Eq. 11), [F]j, [R]j, and [S]j :

$$[\text{PHY}]_j = [\text{F}]_j + [\text{R}]_j + [\text{S}]_j. \qquad (11)$$

The most common way of measuring phytoplankton biomass is using the chlorophyll *a* concentration (µg chl *a* L⁻<fnref n="1" />). A carbon to chlorophyll *a* ratio of 35 mgC µg chl *a*⁻<fnref n="1" /> is therefore considered to convert experimental data into phytoplankton biomass. The initial proportions of different constituents (F, R, S) are fixed (Lancelot et al., 1991). They are only used to determine the initial concentrations of the three cellular constituents and their concentrations in incoming water fluxes for each phytoplankton species. According to Lancelot et al. (1991), the structural and functional macromolecules of the cell ([F]j) account for about 85% of the phytoplankton biomass ([PHY]j), while the reserve products ([R]j) account for about 10% of the biomass. The remainder (5%) of the biomass constitutes the small precursors for macromolecules synthesis ([S]j). These proportions of F, S, and R are updated at each time step for each phytoplankton species.

### 2.3.1 Photosynthesis

The photosynthesis process forms small precursors (S) by fixing carbon dioxide. Its rate is determined by the photosynthesis-irradiance relationship (Platt et al., 1980) including three parameters (Eq. 12) and the active irradiance (I(z), µE m⁻² s⁻<fnref n="1" />).

$$P(z)_{\text{phy}_j} = P_{\text{max20,phy}_j} f(T)_{\text{phy}_j} \left(1 - e^{-\frac{\alpha_{\text{phy}_j} I(z)}{P_{\text{max20,phy}_j} f(T)}}\right) e^{-\frac{\beta_{\text{phy}_j} I(z)}{P_{\text{max20,phy}_j} f(T)}} \quad (12)$$

Here, $P(z)_{\text{phy}_j}$ is the photosynthesis rate of the *j*th phytoplankton species at water depth *z* m [h⁻<fnref n="1" />], $P_{\text{max20,phy}_j}$ is the maximal photosynthesis rate of the *j*th phytoplankton species at 20 °C [h⁻<fnref n="1" />], $f(T)_{\text{phy}_j}$ is the water temperature weight of the *j*th phytoplankton species at T °C defined as in Eq. (2) [-], $\alpha_{\text{phy}_j}$ is the photosynthetic efficiency of the *j*th phytoplankton species [h⁻<fnref n="1" /> (µE m⁻² s⁻<fnref n="1" />)⁻<fnref n="1" />], $\beta_{\text{phy}_j}$ is the photoinhibition capacity of the *j*th phytoplankton species [h⁻<fnref n="1" /> (µE m⁻² s⁻<fnref n="1" />)⁻<fnref n="1" />], and $I(z)$ is photosynthetically active radiation (PAR) or active irradiance in the water column at depth *z* m [µE m⁻² s⁻<fnref n="1" />] or [W m⁻²].

The averaged photosynthesis rate of the *j*th phytoplankton species over water column is obtained by integrating $P(z)$:

$$P_{\text{phy}_j} = \frac{\int_{0}^{\text{depth}} P(z)_{\text{phy}_j} dz}{\text{depth}}, \qquad (13)$$

where depth is the water height [m] and $p_{\text{phy}_j}$ is the averaged photosynthesis rate over the water column [h⁻<fnref n="1" />].

The active irradiance at water depth *z* m ($I(z)$) follows the Beer-Lambert law (Eq. 14). The decrease in active irradiance from the water surface to water bottom is represented by the light extinction coefficient ($\eta$). The extinction coefficient is composed of three parts: pure water ($\eta_{\text{base}}$), suspended solid ($\eta_{\text{ss}}$), and algal self-shading ($\eta_{\text{chl a}}$).

$$I(z) = I_0 e^{-\eta z} \\ \eta = \eta_{\text{base}} + \eta_{\text{chl a}} [\text{chl a}] + \eta_{\text{ss}} [\text{SS}] \qquad (14)$$

Here, $I(z)$ is the active irradiance at water depth *z* [µE m⁻² s⁻<fnref n="1" /> or W m⁻²], $I_0$ is photosynthetically active radiation (PAR) or active irradiance at the water surface measured by the photosynthetic photon flux density (PPFD) [µE m⁻² s⁻<fnref n="1" /> or W m⁻²], $\eta$ is the light extinction coefficient [m⁻<fnref n="1" />], $\eta_{\text{base}}$ is the light extinction coefficient related to pure water [m⁻<fnref n="1" />], $\eta_{\text{chl a}}$ is the linear algal self-shading light extinction coefficient [m⁻<fnref n="1" /> (µg chl *a* L⁻<fnref n="1" />)⁻<fnref n="1" />], $\eta_{\text{ss}}$ is the light extinction coefficient related to suspended solid [m⁻<fnref n="1" /> (mg L⁻<fnref n="1" />)⁻<fnref n="1" />], [chl *a*] is the total chlorophyll *a* concentration [µg chl *a* L⁻<fnref n="1" />], and [SS] is the suspended solid concentration [mg L⁻<fnref n="1" />].

### 2.3.2 Growth

The growth of phytoplankton involves the transformation of small precursors (S) into structural and functional macro-

<footnote n="1">¹Photons within the range of visible light between 400 and 700 nm.</footnote>
molecules (F), which also requires the uptake of dissolved inorganic nutrients (nitrogen – N, phosphorus – P, silicon – Si) from the environment (Fig. 3). The nutrients can potentially limit phytoplankton growth if their quantities are insufficient. The limitation of nutrients is represented using multiple Monod functions (Eq. 15). The maximum growth rate ($\mu_{max,F_j}$) is itself weighted by a limitation based on the availability of small precursors (S). The limitation by dissolved silica (DSi) is applied only for diatoms (DIA).

$$ \mu_{F_j} = \mu_{max20,F_j} f(T)_{phy_j} \left( \frac{\frac{[S_j]}{[F_j]} + K_{S,phy_j}}{\frac{[S_j]}{[F_j]} + K_{F_j}} \right) \text{Nut\_lim} \quad (15) $$

$$ \text{Nut\_lim} = \min \left( \frac{[\text{DIN}]}{[\text{DIN}] + K_{\text{N,phy}_j}}, \frac{[\text{DIP}]}{[\text{DIP}] + K_{\text{P,phy}_j}}, \frac{[\text{DSi}]}{[\text{DSi}] + K_{\text{Si,phy}_j}} \right) $$

$$ \text{or Nut\_lim} = \min \left( \frac{[\text{DIN}]}{[\text{DIN}] + K_{\text{N,phy}_j}}, \frac{[\text{DIP}]}{[\text{DIP}] + K_{\text{P,phy}_j}} \right) \quad (16) $$

Here, $\mu_{F_j}$ is the effective growth rate of functional macromolecules for the jth phytoplankton species [h⁻¹], $\mu_{max20,F_j}$ is the maximal growth rate of functional macromolecules for the jth phytoplankton species at 20 °C [h⁻¹], $f(T)_{phy_j}$ is the water temperature weight of the jth phytoplankton species at T °C defined as in Eq. (2) [-], $[S_j]$ and $[F_j]$ are the concentrations of small precursors and functional macromolecules for the jth phytoplankton species [mgC L⁻¹], $K_{S,phy_j}$ is the half-saturation constant for small precursors of the jth phytoplankton species [-], Nut_lim is the nutrient-limiting factor [-], [DIN], [DIP], and [DSi] are the concentrations of dissolved inorganic nitrogen ([DIN] = [NO₃⁻] + [NH₄⁺], mgN L⁻¹), dissolved inorganic phosphorus ([DIP] = [PO₄³⁻], mgP L⁻¹), and dissolved silica (DSi, mgSi L⁻¹), $K_{N,phy_j}$ and $K_{P,phy_j}$ are the half-saturation constants for dissolved inorganic nitrogen and phosphorus of the jth phytoplankton species [mgN L⁻¹ and mgP L⁻¹], and $K_{Si,phy_j}$ is the half-saturation constant for dissolved silica in the case of diatoms [mgSi L⁻¹].

### 2.3.3 Respiration

The respiration rate of phytoplankton ($r_{phy}$) is divided into two components (Eq. 17): one ($R_{m,phy}$) ensuring the survival of the cell (maintenance process) and the other ($R_{\mu,phy}$) corresponding to the energetic cost of growth.

$$ r_{phy_j} = R_{m20,phy_j} f(T)_{phy_j} + \mu_{F_j} R_{\mu,phy_j} \quad (17) $$

Here, $r_{phy_j}$ is the respiration rate of the jth phytoplankton species [h⁻¹], $R_{m20,phy_j}$ is the maintenance respiration rate of the jth phytoplankton species at 20 °C [h⁻¹], $f(T)_{phy_j}$ is the water temperature weight of the jth phytoplankton species at T °C defined as in Eq. (2) [-], $R_{\mu,phy_j}$ is the respiration rate of the jth phytoplankton species [h⁻¹], $\mu_{\mu,phy_j}$ is the effective growth rate of functional macromolecules for the jth phytoplankton species [h⁻¹], $\mu_{max20,phy_j}$ is the maximal growth rate of functional macromolecules for the jth phytoplankton species [h⁻¹], $E_{cst,phy_j}$ is the basic excretion rate of the jth phytoplankton species [h⁻¹], $E_{phot,phy_j}$ is the excretion of the jth phytoplankton species related to photosynthesis [-], and $p_{phy_j}$ is the photosynthesis rate of the jth phytoplankton species (Eq. 13) [h⁻¹].

species at T °C defined as in Eq. (2) [-], $R_{\mu,phy_j}$ is the respiration for the energetic cost of the jth phytoplankton species [-], and $\mu_{F_j}$ is the effective growth rate of the jth phytoplankton species (Eq. 15) [h⁻¹].

### 2.3.4 Excretion

Included later by Garnier et al. (1998), the phytoplankton excretion ($e_{phy}$) includes two terms: a constant excretion rate ($E_{cst,phy}$) and another that depends on the photosynthesis rate ($E_{phot,phy}$). The product of excretion is the small monomeric substrate (SMS), assimilated directly by heterotrophic bacteria for their growth and respiration (Fig. 1).

$$ e_{phy_j} = E_{cst,phy_j} + p_{phy_j} E_{phot,phy_j} \quad (18) $$

Here, $e_{phy_j}$ is the excretion rate of the jth phytoplankton species [h⁻¹], $E_{cst,phy_j}$ is the basic excretion rate of the jth phytoplankton species [h⁻¹], $E_{phot,phy_j}$ is the excretion of the jth phytoplankton species related to photosynthesis [-], and $p_{phy_j}$ is the photosynthesis rate of the jth phytoplankton species (Eq. 13) [h⁻¹].

The variation of small monomeric substrate (SMS) can then be established (Eq. 19).

$$ \frac{d[\text{SMS}]}{dt} = \text{hydr} - \sum_{i} b_{\text{hb}_i} [\text{HB}_i] + \sum_{j} e_{\text{phy}_j} [\text{F}_j] \quad (19) $$

Here, hydr is the hydrolysis of the dissolved high-weight polymers HD₁ and HD₂ (Eq. 6) [mgC L⁻¹ h⁻¹], $b_{hb_i}$ is the effective rate of substrate uptake by the ith heterotrophic bacteria species (Eq. 3) [h⁻¹], $[HB_i]$ is the biomass concentration of the ith species of heterotrophic bacteria [mgC L⁻¹], $e_{phy_j}$ is the effective excretion rate of the jth phytoplankton species (Eq. 18) [h⁻¹], and $[F_j]$ is the concentration of functional macromolecules for the jth phytoplankton species [mgC L⁻¹].

### 2.3.5 Synthesis and catabolysis of reserve products

The carbon fixed in the cell by photosynthesis forms small precursors (S) that can be transformed, either into functional macromolecules (F) or into reserve products (R). The synthesis of reserve products is limited by the $\frac{[S]}{[F]}$ ratio based on a Michaelis–Menten-like function (Eq. 20).

$$ s_{R,phy_j} = s_{R,max20,phy_j} f(T)_{phy_j} \frac{\frac{[S_j]}{[F_j]} + K_{S,phy_j}}{\frac{[S_j]}{[F_j]} + K_{F_j}} \quad (20) $$

Here, $s_{R,phy_j}$ is the synthesis rate of reserve products of jth phytoplankton species [h⁻¹], $s_{R,max20,phy_j}$ is the maximal synthesis rate of reserve products of jth phytoplankton species at 20 °C [h⁻¹], $f(T)_{phy_j}$ is the water temperature weight of jth phytoplankton species at T °C defined as in Eq. (2) [-], $R_{\mu,phy_j}$ is the respiration rate of the jth phytoplankton species [h⁻¹], $\mu_{\mu,phy_j}$ is the effective growth rate of functional macromolecules for the jth phytoplankton species [h⁻¹], $\mu_{max20,phy_j}$ is the maximal growth rate of functional macromolecules for the jth phytoplankton species [h⁻¹], $E_{cst,phy_j}$ is the basic excretion rate of the jth phytoplankton species [h⁻¹], $E_{phot,phy_j}$ is the excretion of the jth phytoplankton species related to photosynthesis [-], and $p_{phy_j}$ is the photosynthesis rate of the jth phytoplankton species (Eq. 13) [h⁻¹].
as in Eq. (2) [-], [S<sub>j</sub>] and [F<sub>j</sub>] are the concentrations of small precursors and functional macromolecules for the jth phytoplankton species [mgCL<sup>-1</sup>], and K<sub>S,phyj</sub> is the half-saturation constant for small precursors of the jth phytoplankton species [-].

Reserve products (R) are likely to be catabolized to pro-
duce small precursors (S). A first-order kinetic (c<sub>R,phy</sub>, h<sup>-1</sup>)
is used to represent catabolysis of a reserve product.

2.3.6 Extinction of phytoplankton

Three ways of phytoplankton extinction are implemented in
the unified RIVE v1.0: lysis, sinking, and grazing by zoo-
plankton (Sect. 2.4.1). The phytoplankton lysis is represented
by first-order kinetics using a mortality rate (k<sub>d20,phy</sub>, h<sup>-1</sup>).
For ease of presentation, all three processes are assumed in
an overall extinction rate d<sub>phy</sub> (h<sup>-1</sup>).

$$
d_{\text{phy}_j} = k_{\text{d20,phy}_j} f(T)_{\text{phy}_j} + k_{\text{sink,phy}_j} \\
+ \sum_i \frac{b_{\text{zoo}_i} [\text{ZOO}_i]}{\sum_{k=1}^{\text{NS}} [\text{PHY}_k]} \quad (21)
$$

Here, $d_{\text{phy}_j}$ is the extinction rate of the $j$th phytoplankton species [h$^{-1}$], $k_{\text{d20,phy}_j}$ is the mortality rate of the $j$th phytoplankton species at 20 °C [h$^{-1}$], $f(T)_{\text{phy}_j}$ is the water temperature weight of $j$th phytoplankton species at $T$ °C defined as in Eq. (2) [-], $k_{\text{sink,phy}_j}$ is the sinking rate of the $j$th phytoplankton species [h$^{-1}$], $b_{\text{zoo}_i}$ is the grazing rate of the $i$th zooplankton species (Eq. 26, Sect. 2.4.1) [h$^{-1}$], [ZOO$_i$] is the zooplankton concentration of the $i$th zooplankton species [mgCL$^{-1}$], and $\sum_{k=1}^{\text{NS}} [\text{PHY}_k]$ is the total phytoplankton concentration (with NS the number of phytoplankton species grazed by zooplankton) [mgCL$^{-1}$].

2.3.7 Phytoplankton budgets

According to the processes related to phytoplankton (photo-
synthesis, growth, mortality, etc.), different budgets can be
established for the jth phytoplankton species as follows.

$$
\begin{align}
\frac{d[\mathrm{S}_j]}{dt} ={}& (p_{\mathrm{phy}_j} - r_{\mathrm{phy}_j} - \mu_{\mathrm{F}_j} - s_{\mathrm{R,phy}_j})[\mathrm{F}_j] \nonumber \\
& + c_{\mathrm{R,phy}_j}[\mathrm{R}_j] - e_{\mathrm{phy}_j}[\mathrm{F}_j] - d_{\mathrm{phy}_j}[\mathrm{S}_j] \tag{22}
\end{align}
$$

$$
\frac{d[\mathrm{R}_j]}{dt} = s_{\mathrm{R,phy}_j}[\mathrm{F}_j] - c_{\mathrm{R,phy}_j}[\mathrm{R}_j] - d_{\mathrm{phy}_j}[\mathrm{R}_j] \quad (23)
$$

$$
\frac{d[\mathrm{F}_j]}{dt} = (\mu_{\mathrm{F}_j} - d_{\mathrm{phy}_j})[\mathrm{F}_j] \quad (24)
$$

$$
\frac{d[\text{PHY}_j]}{dt} = (p_{\text{phy}_j} - r_{\text{phy}_j} - e_{\text{phy}_j})[F_j] \\ - d_{\text{phy}_j}[\text{PHY}_j] \qquad (25)
$$

Here, *p*<sub>phy<sub>j</sub></sub> is the photosynthesis rate of the <sub>j</sub>th phytoplank-
ton species (Eq. 13) [h<sup>-1</sup>], *r*<sub>phy<sub>j</sub></sub> is the respiration rate of the
<sub>j</sub>th phytoplankton species (Eq. 17) [h<sup>-1</sup>], *μ*<sub>F<sub>j</sub></sub> is the growth

<figure><img src="image_6.png" /><figcaption>**Figure 4.** Dynamics of zooplankton. PHY: phytoplankton species; ZOO: zooplankton species; respi.: respiration; sink.: sinking.</figcaption></figure>

rate of the jth phytoplankton species (Eq. 15) [h⁻¹], s<sub>R,phyj</sub> is the synthesis rate of reserve products of the jth phytoplankton species (Eq. 20) [h⁻¹], c<sub>R,phyj</sub> is the catabolysis rate of reserve products of the jth phytoplankton species [h⁻¹], e<sub>phyj</sub> is the excretion rate of the jth phytoplankton species (Eq. 18) [h⁻¹], d<sub>phyj</sub> is the extinction rate of the jth phytoplankton species (Eq. 21) [h⁻¹], [S<sub>j</sub>], [F<sub>j</sub>], and [R<sub>j</sub>] are the concentrations of [S<sub>j</sub>], [F<sub>j</sub>], and [R<sub>j</sub>] for the jth phytoplankton species [mgCL⁻¹], and [PHY<sub>j</sub>] is the biomass concentration of the jth phytoplankton species [mgCL⁻¹].

2.4 Zooplankton dynamics

The zooplankton dynamics include grazing on phytoplank-
ton as well as growth, respiration, mortality, and sinking
(Fig. 4).

2.4.1 Grazing and growth

Grazing on phytoplankton by zooplankton and growth of
zooplankton are expressed based on a maximal grazing rate
at 20 °C (b<sub>max20,zoo</sub>) limited by the phytoplankton biomass
based on a Monod function (Eq. 26). The grazing of zoo-
plankton takes place only when the total phytoplankton
biomass exceeds a certain threshold ([PHY<sub>0</sub>]). No specific
preference for grazing on particular phytoplankton species is
considered among zooplankton species. Instead, the phyto-
plankton biomass grazed by the <i>t</i>th species of zooplankton
is divided proportionally among each species of phytoplank-
ton (Eq. 21). The growth rate of zooplankton is considered
proportional to the grazing rate using a growth yield factor
(Eq. 27).

$$b_{zoo_i} = b_{max20,zoo_i} f(T)_{zoo_i} \frac{\left(\sum_{j}^{NS}[PHY_j] - [PHY_0]_{zoo_i}\right)}{\left(\sum_{j}^{NS}[PHY_j] - [PHY_0]_{zoo_i}\right) + K_{phy,zoo_i}} \quad (26)$$

$$\mu_{zoo_i} = Y_{zoo_i} b_{zoo_i} \quad (27)$$

Here, $b_{zoo_i}$ is the effective grazing rate of the ith zooplankton species [h⁻¹], $b_{max20,zoo_i}$ is the maximal grazing rate of the ith zooplankton species at 20°C [h⁻¹], $f(T)_{zoo_i}$ is the water temperature weight of the ith zooplankton species at T °C defined as in Eq. (2) [–], $\sum_{j}^{NS}[PHY_j]$ is the total phytoplankton biomass with NS the number of phytoplankton species grazed by zooplankton [mgC L⁻¹], $[PHY_0]_{zoo_i}$ is the phytoplankton biomass threshold above which grazing takes place for the ith zooplankton species [mgC L⁻¹], $K_{phy,zoo_i}$ is the half-saturation constant for phytoplankton biomass of the ith zooplankton species [mgC L⁻¹], $\mu_{zoo_i}$ is the growth rate of the ith zooplankton species [h⁻¹], and $Y_{zoo_i}$ is the growth yield of the ith zooplankton species [–].

## 2.4.2 Respiration and mortality

Grazed phytoplankton not used for zooplankton growth is respired (Fig. 4). The rate of respiration is then obtained by $(1 - Y_{zoo}) \times b_{zoo}$. The mortality of zooplankton is simulated by first-order kinetics ($k_{d20,zoo}$).

$$r_{zoo_i} = (1 - Y_{zoo_i}) \times b_{zoo_i} \quad (28)$$

$$\frac{d[ZOO_i]}{dt} = (Y_{zoo_i} b_{zoo_i} - k_{d20,zoo_i} f(T)_{zoo_i} - k_{sink,zoo_i}) [ZOO_i] \quad (29)$$

Here, $r_{zoo_i}$ is the respiration rate of the ith zooplankton species [h⁻¹], $Y_{zoo_i}$ is the growth yield of the ith zooplankton species [–], $b_{zoo_i}$ is the effective grazing rate of the ith zooplankton species (Eq. 26) [h⁻¹], $k_{d20,zoo_i}$ is the mortality rate of the ith zooplankton species at 20°C [h⁻¹], $f(T)_{zoo_i}$ is the water temperature weight of the ith zooplankton species at T °C defined as in Eq. (2) [–], $k_{sink,zoo_i}$ is the sinking rate of the ith zooplankton species [h⁻¹], and $[ZOO_i]$ is the biomass concentration of the ith zooplankton species [mgC L⁻¹].

## 2.5 Nutrient cycling

As shown above, several processes related to nutrients are taken into account: uptake by phytoplankton, mineralization, nitrification, and denitrification (Figs. 5 and 6).

### 2.5.1 Uptake of nutrients (N, P, Si) by phytoplankton

The Redfield–Conley stoichiometry (C : N : P : Si = 106 : 16 : 1 : 42; – Redfield et al., 1963; Conley et al., 1989) is used to determine the composition of carbon, nitrogen, and phosphorus in organic matter. Constant C/N, C/P, and C/Si mass

ratios are considered to calculate the uptake of nutrients associated with phytoplankton growth.

$$\frac{d[uptN]}{dt} = \sum_{i} \frac{(\mu_{F_i} + e_{phy_i})[F_i]}{C/N} \quad (30)$$

$$uptN_{4}^{+} = \min\left([NH_{4}^{+}], uptN\left(\frac{[NH_{4}^{+}]}{[NH_{4}^{+}] + [NO_{3}^{-}]}\right)\right)^{0.025} \quad (31)$$

$$uptNO_{3}^{-} = uptN - uptNH_{4}^{+} \quad (32)$$

$$\frac{d[uptP]}{dt} = \sum_{i} \frac{(\mu_{F_i} + e_{phy_i})[F_i]}{C/P} \quad (33)$$

$$\frac{d[uptSi]}{dt} = \frac{\mu_{F,dia}[F_{dia}]}{C/Si} \quad (34)$$

Here, $\mu_{F_i}$ is the effective growth rate of the ith phytoplankton species (Eq. 15) [h⁻¹], $e_{phy_i}$ is the excretion rate of the ith phytoplankton species (Eq. 18) [h⁻¹], $[F_i]$ represents the functional macromolecule concentration of the ith phytoplankton species [mgC L⁻¹], $[NH_4^+]$ and $[NO_3^-]$ are the concentrations of ammonium and nitrate [mgN L⁻¹], $uptN$ is the uptake of nitrogen for phytoplankton growth [mgN L⁻¹], $uptNH_4^+$ is the uptake of $NH_4^+$ for phytoplankton growth [mgN L⁻¹], $uptNO_3^-$ is the uptake of $NO_3^-$ for phytoplankton growth [mgN L⁻¹], $uptP$ is the uptake of phosphorus for phytoplankton growth [mgP L⁻¹], $C/N$ is the carbon to nitrogen mass ratio [mgC mgN⁻¹], $C/P$ is the carbon to phosphorus mass ratio [mgC mgP⁻¹], $C/Si$ is the carbon to silica mass ratio [mgC mgSi⁻¹], $uptSi$ is the uptake of silica for diatom growth [mgSi L⁻¹], $\mu_{F,dia}$ is the effective growth rate of diatoms [h⁻¹], and $[F_{dia}]$ is the functional macromolecule (F) concentration of diatoms [mgC L⁻¹]

### 2.5.2 Release of nutrients by mineralization

The mineralization of organic matter by heterotrophic bacteria and zooplankton is achieved by its oxidation through respiration (Fig. 5). The process consumes organic matter and releases nitrogen and phosphorus from the fraction that is not assimilated for growth of heterotrophic bacteria and zooplankton.

$$\text{respHB} = \sum_{i} (1 - Y_{hb_i}) b_{hb_i} [HB_i] \quad (35)$$

$$\text{respZOO} = \sum_{j} (1 - Y_{zoo_j}) b_{zoo_j} [ZOO_j] \quad (36)$$

$$\text{reln} = \frac{\text{respHB}}{C/N} + \frac{\text{respZOO}}{C/N} \quad (37)$$

$$\text{relP} = \frac{\text{respHB}}{C/P} + \frac{\text{respZOO}}{C/P} \quad (38)$$

Here, respHB is the respiration of heterotrophic bacteria species [mgC L⁻¹ h⁻¹], $Y_{hb_i}$ is the growth yield of the ith heterotrophic bacteria species [–], $b_{hb_i}$ is the effective rate of substrate uptake by the ith heterotrophic bacteria species (Eq. 1) [h⁻¹], $[HB_i]$ is the concentration of
<figure><img src="image_7.png" /><figcaption>Figure 5. Cycling of nitrogen. PHY: phytoplankton species; HB: heterotrophic bacteria; ZOO: zooplankton species; AOB: ammonia-oxidizing bacteria; NOB: nitrite-oxidizing bacteria; respi.: respiration; excr.: excretion; denit: denitrification.</figcaption></figure>

the ith heterotrophic bacteria species [mgC L⁻¹], respZOO is the respiration of zooplankton species [mgC L⁻¹ h⁻¹], Yzooᵢ is the growth yield of the jth zooplankton species [–], bzooᵢ is the effective grazing rate of the jth zooplankton species (Eq. 26) [h⁻¹], [ZOOᵢ] is the concentration of the jth zooplankton species [mgC L⁻¹], relN is the release of nitrogen [mgN L⁻¹ h⁻¹], C/N is the carbon to nitrogen mass ratio [mgC mgN⁻¹], relP is the release of phosphorus [mgP L⁻¹ h⁻¹], and C/P is the carbon to phosphorus mass ratio [mgC mgP⁻¹].

### 2.5.3 Nitrification and denitrification

As mentioned in Sect. 2.2, the nitrification process (Figs. 2 and 5) is related to the growth of AOB (ammonia-oxidizing bacteria) and NOB (nitrite-oxidizing bacteria). Growth yields (Y$_{aob,i}$ and Y$_{nob,i}$) are used to describe the amount of nitrogen consumed by nitrifying bacteria (Eqs. 39 and 40). Denitrification occurs when dissolved oxygen is not present in sufficient quantity (Fig. 5).

$$ \text{nitr}_{\text{aob}} = \sum_{i} \frac{\mu_{\text{aob}_{i}}}{Y_{\text{aob}_{i}}} [\text{AOB}_{i}] \quad (39) $$

$$ \text{nitr}_{\text{nob}} = \sum_{j} \frac{\mu_{\text{nob}_{j}}}{Y_{\text{nob}_{j}}} [\text{NOB}_{j}] \quad (40) $$

Here, μ$_{aob,i}$ and μ$_{nob,i}$ are the effective growth rates of the ith AOB species (Eq. 7) and the jth NOB species (Eq. 8)

[h⁻¹], Y$_{aob,i}$ and Y$_{nob,i}$ are the growth yields of the ith AOB species and the jth NOB species [mgC mgN⁻¹], nitr$_{aob}$ is nitrification calculated as NH₄⁺ + 3/2O₂ → NO₂⁻ + H₂O + 2H⁺ [mgN L⁻¹ h⁻¹], nitr$_{nob}$ is nitrification calculated as NO₂⁻ + 1/2O₂ → NO₃⁻ [mgN L⁻¹ h⁻¹], and [AOBᵢ] and [NOBᵢ] are the biomass concentrations of the ith AOB species and the jth NOB species [mgC L⁻¹].

The budgets of NO₃⁻, NH₄⁺, and NO₂⁻ can then be established.

$$ \frac{d[\text{NO}_3^-]}{dt} = -\text{denit} + \text{nitr}_{\text{nob}} - \frac{\text{uptNO}_3^-}{dt} \quad (41) $$

$$ \frac{d[\text{NH}_4^+]}{dt} = \text{relN} - \text{nitr}_{\text{aob}} - \frac{\text{uptNH}_4^+}{dt} \quad (42) $$

$$ \frac{d[\text{NO}_2^-]}{dt} = \text{nitr}_{\text{aob}} - \text{nitr}_{\text{nob}} \quad (43) $$

Here, denit is denitrification [mgN L⁻¹ h⁻¹], nitr$_{nob}$ is nitrification by NOB (Eq. 40) [mgN L⁻¹ h⁻¹], nitr$_{aob}$ is nitrification by AOB (Eq. 39) [mgN L⁻¹ h⁻¹], uptNO₃⁻/dt is the uptake of NO₃⁻ by phytoplankton growth (Eq. 32) [mgN L⁻¹ h⁻¹], relN is the release of nitrogen by respiration of heterotrophic bacteria and zooplankton (Eq. 37) [mgN L⁻¹ h⁻¹], and uptNH₄⁺/dt is the uptake of NH₄⁺ by phytoplankton growth (Eq. 31) [mgN L⁻¹ h⁻¹].

### 2.5.4 Phosphate adsorption desorption

Orthophosphate (PO₄³⁻) is released by mineralization and taken up by phytoplankton exactly as inorganic nitrogen (Fig. 6). Once released in the water column, however, orthophosphates are subject to a process of adsorption-desorption on mineral suspended solids (MSSs) to form PIP (particulate inorganic phosphorus). In addition, the impact of sediment dynamics on P fluxes should be considered in future work with unified RIVE. Vilmin et al. (2015a) showed that P fluxes are mainly driven by hydrological conditions and sediment-related processes in the Seine River system.

The process is represented according to an instantaneous hyperbolic equilibrium relationship of the form

$$ \frac{[\text{PIP}]}{[\text{MSS}]} = P_{\text{ac}} \times \frac{[\text{PO}_4^{3-}]}{[\text{PO}_4^{3-}] + K_{\text{ps}}}, \quad (44) $$

with [PIP]/[MSS] the inorganic P content of MSS [mgP mgMSS⁻¹], [PIP] and [MSS] the concentrations of PIP and MSS [mgP L⁻¹ and mgMSS L⁻¹], P$_{ac}$ the maximum adsorption capacity of MSS [mgP mgMSS⁻¹], [PO₄³⁻] the concentration of orthophosphate [mgP L⁻¹], and K$_{ps}$ the half-saturation adsorption constant [mgP L⁻¹].

Considering this equilibrium to be instantaneously reached implies that a relationship exists between the variables PIP, MSS, PO₄³⁻, and TIP (total inorganic phosphorus):

$$ [\text{TIP}] = [\text{PO}_4^{3-}] + [\text{PIP}]. \quad (45) $$
<figure><img src="image_8.png" /><figcaption>Figure 6. Phosphorus and silica dynamics. HD: dissolved high-weight polymer; PHY: phytoplankton; PIP: particulate inorganic phosphorus; MSS: mineral suspended solids; BSi: biogenic silica; DSi: dissolved silica; adsorp.: adsorption; desorp.: desorption; sink.: sinking; disso. dissolution.</figcaption></figure>

This equilibrium relationship can be written as

$$[PO_4^{3-}] = \frac{([TIP] - P_{ac} \times [MSS] - K_{ps}) + \sqrt{(-[TIP] + P_{ac} \times [MSS] + K_{ps})^2 + 4 \times [TIP] \times K_{ps}}}{2}, \quad (46)$$

with $[PO_4^{3-}]$ the concentration of orthophosphate [mgPL$^{-1}$], [TIP] and [MSS] the concentrations of TIP and MSS [mgPL$^{-1}$ and mgMSSL$^{-1}$], $P_{ac}$ the maximal adsorption capacity of MSS [mgPmgMSS$^{-1}$], and $K_{ps}$ the half-saturation adsorption constant [mgPL$^{-1}$].

### 2.5.5 Silica dynamics

Dissolved silica (DSi) is produced by the dissolution of dead frustules of diatoms (designated as biogenic silica, BSi). Rock weathering contributes also dissolved silica, while it is considered null in unified RIVE v1.0. DSi is taken up by the growth of diatoms (Fig. 6). Biogenic silica is produced by the lysis and grazing of diatoms; it settles down and dissolves according to first-order kinetics, dependent on water temperature (Rickert et al., 2002):

$$BSi_{disso.} = K_{bSi} \times [BSi], \quad (47)$$

$$K_{bSi} = K_{bSi20} \times ftp_{Si}(T), \quad (48)$$

$$ftp_{Si}(T) = \exp\left(\frac{60000}{8.314} \times \left(\frac{1}{275} - \frac{1}{273 + T}\right)\right), \quad (49)$$

with $K_{bSi20}$ the dissolution rate of biogenic silica at 20°C [h$^{-1}$], [BSi] the concentration of biogenic silica [mgSiL$^{-1}$], and T the water temperature [°C].

## 2.6 Dissolved oxygen

Dissolved oxygen is especially influenced by photosynthesis and respiration. Re-aeration at the water-air interface is also included in the unified RIVE v1.0 model. Sediment dynamics are important for sediment oxygen demand (not shown here). Vilmin et al. (2016) showed that benthic respiration accounts for one-third of the total Seine River respiration. Relevant efforts toward sediment dynamics need to be made in future work, which is not the focus of this study. An oxygen budget can then be established (Eq. 50).

$$\frac{d[O_2]}{dt} = \text{rea} + \frac{32}{12} \left( \sum_i (p_{\text{phy}_i} - r_{\text{phy}_i}) [F_i] - \text{respHB} - \sum_j r_{\text{zoo}_j} [\text{ZOO}_j] \right) - \frac{32}{14} \left( \frac{3}{2} \text{nitr}_{\text{aob}} + \frac{1}{2} \text{nitr}_{\text{nob}} \right) \quad (50)$$

$$\text{rea} = \frac{k_{\text{rea}}}{\text{depth}} ([\text{O}_2]_{\text{sat}} - [\text{O}_2]) \quad (51)$$

Here, $k_{rea}$ is the re-aeration coefficient [m h$^{-1}$], “depth” indicates water height [m], $[O_2]_{sat}$ is the saturated concentration of dissolved oxygen in water [mgO$_2$ L$^{-1}$], $[O_2]$ is the concentration of dissolved oxygen in water [mgO$_2$ L$^{-1}$], $\frac{32}{12}$ is the molar mass ratio between dissolved oxygen and carbon [mgO$_2$ mgC$^{-1}$], $p_{phy_i}$ is the photosynthesis rate of the $i$th phytoplankton species (Eq. 13) [h$^{-1}$], $r_{phy_i}$ is the respiration rate of the $i$th phytoplankton species (Eq. 17) [h$^{-1}$], $[F_i]$ is the functional biomass concentration of the $i$th phytoplankton species [mgCL$^{-1}$], respHB is the respiration of all heterotrophic bacteria species (Eq. 35) [mgCL$^{-1}$ h$^{-1}$],
$r_{zoo_j}$ is the respiration rate of the jth zooplankton species (Eq. 28) [h⁻¹], [ZOO_j] is the biomass concentration of the jth zooplankton species [mgC L⁻¹], 32/14 is the molecular mass ratio between dissolved oxygen and nitrogen [mgO₂ mgN⁻¹], $nitr_{aob}$ is nitrification to produce nitrite by oxidizing NH₄⁺ (Eq. 39, with 3/2 the stoichiometric coefficient) [mgN L⁻¹ h⁻¹], and $nitr_{nob}$ is nitrification to produce nitrate by oxidizing NO₂⁻ (Eq. 40, with 1/2 the stoichiometric coefficient) [mgN L⁻¹ h⁻¹].

## 2.7 Inorganic carbon

An inorganic carbon module is implemented in unified RIVE v1.0. The carbonate system is described by a set of equations (named the CO₂ module) based on a previous representation provided by Gypens et al. (2004) and adapted for freshwater environments (Marescaux et al., 2020). In this module, four state variables are defined: dissolved inorganic carbon – DIC, total alkalinity – TA, acidity – pH, and aqueous carbon dioxide – CO₂((aq)).

### 2.7.1 CO₂ flux at the air-water interface

The DIC is defined as the sum of three dissolved carbonate species:

$$[\text{DIC}] = [\text{H}_2\text{CO}_3] + [\text{HCO}_3^-] + [\text{CO}_3^{2-}]. \quad (52)$$

The calculation of pH is derived from Culberson (1980) using TA and DIC. Then the aqueous carbon dioxide (CO₂ (aq)) is derived from the carbonate chemical equilibrium using DIC and pH (Marescaux et al., 2020; Yan et al., 2022a).

$$[\text{CO}_2(\text{aq})] = \frac{[\text{DIC}]_{\frac{[\text{H}^+]}{K_1}}}{\left(1 + \frac{[\text{H}^+]}{K_1} + \frac{K_2}{[\text{H}^+]}\right)} \quad (53)$$

Here, $K_1$ and $K_2$ are the equilibrium constants of carbonate equilibrium reactions (Stumm and Morgan, 1996) [mol L⁻¹], [H⁺] is the concentration of hydrogen ions with pH = -log([H⁺]) [mol L⁻¹], and [DIC] is the concentration of dissolved inorganic carbon [mgC L⁻¹].

The flux of CO₂ at the water-air interface ($F_{CO_2}$, gC m⁻² h⁻¹) is calculated based on Fick's first law (Fick, 1855) with a gas transfer velocity of CO₂ ($k_{CO_2}$).

$$F_{\text{CO}_2} = k_{\text{CO}_2} ([\text{CO}_2(\text{sat})] - [\text{CO}_2(\text{aq})]) \quad (54)$$

Here, $k_{CO_2}$ is the gas transfer velocity of CO₂ [m h⁻¹], [CO₂(sat)] is the solubility of CO₂ in water calculated based on Henry's law (Weiss, 1974) [mgC L⁻¹], and [CO₂(aq)] is the aqueous carbon dioxide concentration [mgC L⁻¹].

The gas transfer velocity of CO₂ ($k_{CO_2}$) depends on water temperature and $k_{600}$ (gas transfer velocity of CO₂ for a Schmidt number of 600, corresponding to a temperature of 20 °C in freshwater). According to Wilke and Chang (1955), Jähne et al. (1987), and Wanninkhof (1992), the gas transfer

velocity of CO₂ ($k_{CO_2}$) at water temperature T (°C) can be calculated as

$$k_{\text{CO}_2} = k_{600} \sqrt{\frac{600}{S_{\text{CO}_2}(T)}}, \quad (55)$$

where $k_{600}$ (m h⁻¹) is the gas transfer velocity of CO₂ for a Schmidt number of 600, and $S_{CO_2}(T)$ is the Schmidt number (dimensionless) calculated with the water temperature in degrees Celsius (°C). The $S_{CO_2}(T)$ can be determined as,=

$$S_{\text{CO}_2}(T) = 1911.1 - 118.11T + 3.4527T^2 - 0.04132T^3. \quad (56)$$

### 2.7.2 Budgets of TA and DIC

Processes such as respiration, photosynthesis, nitrification, denitrification, and input flows affect TA and DIC. The unified RIVE v1.0 considers these processes explicitly.

$$\frac{\mathrm{d} \mathrm{TA}}{\mathrm{dt}} = \left[ \left( \frac{14}{106} \times \frac{(\mathrm{respPHY} + \mathrm{respHB} + \mathrm{respZOO})}{12} \right) \right. \\ + \frac{(\mathrm{denit} - 2 \cdot \mathrm{nitr}_{\mathrm{aob}})}{14} + \left( \frac{17}{106} \times \frac{\mathrm{uptNO}_3^-}{\mathrm{uptN}} - \frac{15}{106} \right) \\ \left. \times \frac{\mathrm{uptNH}_4^+}{\mathrm{uptN}} \right) \times \frac{\sum (\mu_{\mathrm{F}_i} + e_{\mathrm{phy}_i})[\mathrm{F}_i]}{12} \\ \times 1000 + \mathrm{TA}_{\mathrm{Net\_Input}} \quad (57)$$

$$\frac{\mathrm{d} \mathrm{DIC}}{\mathrm{dt}} = (\mathrm{respPHY} + \mathrm{respHB} + \mathrm{respZoo}) + \mathrm{denit} \\ \times \frac{12}{14} \times \frac{5}{4} - \sum p_{\mathrm{phy}_i} [\mathrm{F}_i] + \frac{F_{\mathrm{CO}_2}}{\mathrm{depth}} \\ + \mathrm{DIC}_{\mathrm{Net\_Input}} \quad (58)$$

Here, $TA_{Net\_Input}$ (µmol L⁻¹ h⁻¹) and $DIC_{Net\_Input}$ (mgC L⁻¹ h⁻¹) are the net input fluxes. The respiration of all phytoplankton, bacteria, and zooplankton species (respPHY, respHB, respZOO; mgC L⁻¹ h⁻¹) transforms organic carbon to CO₂ by full oxidization. The denitrification (denit, mgN L⁻¹ h⁻¹) is also considered in the calculation of TA and DIC. $F_{CO_2}$ (gC m⁻² h⁻¹) is the CO₂ flux at the air-water interface, and “depth” is the water depth (m). 14/106, 17/106, 15/106, and 5/4 are the stoichiometry coefficients of biogeochemical processes (Marescaux et al., 2020).

## 2.8 Kinetic parameters in the unified RIVE model

A total of 120 parameters are used to describe the aforementioned processes considering three phytoplankton species and two heterotrophic bacteria species. Some of them depend on water temperature and are calculated with a water temperature function (Eq. 2). Their definitions and reference values are provided in Appendix B.
<figure><img src="image_9.png" /><figcaption>**Figure 7.** Simulation by the HSB model (unified RIVE v1.0) of the dynamics of heterotrophic bacteria in a filtered and reinoculated sample from drainage pond water (Seine basin, France) taken in February 2021 (Garnier et al., 2021). DOM: dissolved organic matter; SHB: small heterotrophic bacteria.</figcaption></figure>

# 3 Results

## 3.1 Digital implementation with Python 3 (pyRIVE 1.0) or ANSI C (C-RIVE 0.32)

The above unified governing equations are implemented in Python 3 to create pyRIVE 1.0 (https://doi.org/10.48579/PRO/Z9ACP1; Thieu et al., 2023) and in ANSI C to create C-RIVE 0.32 (https://doi.org/10.5281/zenodo.7849609; Wang et al., 2023b), respectively. A Jupyter Notebook is used for pedagogical exercises with pyRIVE 1.0, while C-RIVE 0.32 needs to be compiled with gcc under a Linux or a MAC OS operating system. In addition, the user interface of C-RIVE 0.32 uses its own parser based on flex and bison, which allows the software to read ASCII files.

In practice, the number of living species is predefined in pyRIVE 1.0, while we have the ability to define as many species as desired in C-RIVE 0.32 (Table 1). For instance, three communities of phytoplankton (DIA: diatoms; GRA: green algae; CYA: cyanobacteria), two populations of heterotrophic bacteria distinct in their growth rate and size (small one SHB and large one LHB; Garnier et al., 1992a), and two zooplankton communities (ZOR: rotifer and ZOC, micro-crustaceans; Billen et al., 1994; Garnier et al., 1995, 2000) are predefined in pyRIVE 1.0.

In addition, the TIP (total inorganic phosphorus) is considered to be a state variable in pyRIVE 1.0. $PO_4^{3-}$ and PIP are derived from it according to Eqs. (45) and (46). TIP is subject to release by heterotrophic bacteria and zooplankton respiration (Eq. 38), uptake by phytoplankton, and settling of PIP together with MSS. However, the $PO_4^{3-}$ is treated as a state variable and released by respiration (Eq. 38) in C-RIVE 0.32 and only PIP (particulate inorganic phosphorus) is derived from the Eq. (44).

<figure><img src="image_10.png" /><figcaption>**Figure 8.** Simulation by the HSB model (unified RIVE v1.0) of the dynamics of heterotrophic bacteria in a filtered and reinoculated sample from urban sewage water (Rosny-sur-Seine, France) taken in February 2021 (Garnier et al., 2021). DOM: dissolved organic matter; SHB: small heterotrophic bacteria; LHB: large heterotrophic bacteria.</figcaption></figure>

<figure><img src="image_11.png" /><figcaption>**Figure 9.** Geometric and hydraulic description of a river stretch.</figcaption></figure>

## 3.2 Modeling the organic matter degradation by unified RIVE v1.0 (HSB model)

The ability of the HSB model (Fig. 1) to simulate organic matter degradation has been verified by modeling two batch experiments conducted by Garnier et al. (2021). Two water samples were used in the study. One sample was obtained from a drainage pond in the Seine basin, France, in February 2021 (Fig. 7). The other sample was taken from an urban sewage collector in Rosny-sur-Seine, France, in February 2021 (Fig. 8). These samples were incubated in the dark at a temperature of 21 °C for a period of 45 d, during which aerobic bacteria consumed organic matter (Servais et al., 1995). Only DOM and bacterial biomass are measured during batch experiments and then used to show validation. The HSB model is able to effectively reproduce the concentrations of dissolved organic matter and bacterial biomass with a trial-error adjustment of its parameter values (Figs. 7 and 8). The parameter values are kept the same for both water samples.

## 3.3 A river stretch simulated with unified RIVE v1.0: pyRIVE 1.0 vs. C-RIVE 0.32

A river stretch with a Strahler order of 8 (Fig. 9) is designed to compare the results simulated by two versions of unified RIVE v1.0 implemented in pyRIVE 1.0 and C-RIVE 0.32. The case study allows us to compare the two versions of unified RIVE v1.0 under transient contrasting conditions
<table><caption>**Table 1.** Number of living species defined in pyRIVE 1.0 and C-RIVE 0.32 which implement the unified RIVE v1.0.</caption><thead><tr><th>Species</th><th>PHY</th><th>HB</th><th>AOB</th><th>NOB</th><th>ZOO</th></tr></thead><tbody><tr><td>pyRIVE 1.0</td><td>3</td><td>2</td><td>1</td><td>1</td><td>2</td></tr><tr><td>C-RIVE 0.32</td><td>User-defined</td><td>User-defined</td><td>User-defined</td><td>User-defined</td><td>User-defined</td></tr></tbody></table>

<table><caption>**Table 2.** Initial concentrations and boundary conditions.</caption><thead><tr><th>Species</th><th>Description</th><th>C<sub>init</sub></th><th>C<sub>boundary</sub></th><th>Unit</th></tr></thead><tbody><tr><td>SHB</td><td>Small heterotrophic bacteria</td><td>0.005</td><td>0.005</td><td>[mgC L<sup>-1</sup>]</td></tr><tr><td>LHB</td><td>Large heterotrophic bacteria</td><td>0.004</td><td>0.004</td><td>[mgC L<sup>-1</sup>]</td></tr><tr><td>AOB</td><td>Ammonia-oxidizing bacteria</td><td>0.001</td><td>0.001</td><td>[mgC L<sup>-1</sup>]</td></tr><tr><td>NOB</td><td>Nitrite-oxidizing bacteria</td><td>0.0002</td><td>0.0002</td><td>[mgC L<sup>-1</sup>]</td></tr><tr><td>DIA</td><td>Diatoms</td><td>0.447</td><td>0.447</td><td>[mgC L<sup>-1</sup>]</td></tr><tr><td>GRA</td><td>Green algae</td><td>0.539</td><td>0.539</td><td>[mgC L<sup>-1</sup>]</td></tr><tr><td>CYA</td><td>Cyanobacteria</td><td>0.662</td><td>0.662</td><td>[mgC L<sup>-1</sup>]</td></tr><tr><td>ZOR</td><td>Rotifer</td><td>9.33 × 10<sup>-5</sup></td><td>9.33 × 10<sup>-5</sup></td><td>[mgC L<sup>-1</sup>]</td></tr><tr><td>ZOC</td><td>Micro-crustaceans</td><td>9.33 × 10<sup>-6</sup></td><td>9.33 × 10<sup>-6</sup></td><td>[mgC L<sup>-1</sup>]</td></tr><tr><td>SMS</td><td>Small monomeric substrate</td><td>0.036</td><td>0.036</td><td>[mgC L<sup>-1</sup>]</td></tr><tr><td>DOM<sub>1</sub></td><td>Rapidly biodegradable dissolved organic matter</td><td>0.022</td><td>0.022</td><td>[mgC L<sup>-1</sup>]</td></tr><tr><td>DOM<sub>2</sub></td><td>Slowly biodegradable dissolved organic matter</td><td>0.174</td><td>0.174</td><td>[mgC L<sup>-1</sup>]</td></tr><tr><td>DOM<sub>3</sub></td><td>Dissolved refractory organic matter</td><td>1.625</td><td>1.625</td><td>[mgC L<sup>-1</sup>]</td></tr><tr><td>POM<sub>1</sub></td><td>Rapidly biodegradable particulate organic matter</td><td>0.005</td><td>0.005</td><td>[mgC L<sup>-1</sup>]</td></tr><tr><td>POM<sub>2</sub></td><td>Slowly biodegradable particulate organic matter</td><td>0.021</td><td>0.021</td><td>[mgC L<sup>-1</sup>]</td></tr><tr><td>POM<sub>3</sub></td><td>Particulate refractory organic matter</td><td>0.107</td><td>0.107</td><td>[mgC L<sup>-1</sup>]</td></tr><tr><td>NH4</td><td>Ammonium</td><td>1.5</td><td>1.5</td><td>[mgN L<sup>-1</sup>]</td></tr><tr><td>NO2</td><td>Nitrite</td><td>0.016</td><td>0.016</td><td>[mgN L<sup>-1</sup>]</td></tr><tr><td>NO3</td><td>Nitrate</td><td>0.941</td><td>0.941</td><td>[mgN L<sup>-1</sup>]</td></tr><tr><td>TIP</td><td>Total inorganic phosphorus</td><td>0.2</td><td>0.2</td><td>[mgP L<sup>-1</sup>]</td></tr><tr><td>DSi</td><td>Dissolved silica</td><td>3.090</td><td>3.090</td><td>[mgSi L<sup>-1</sup>]</td></tr><tr><td>MSS</td><td>Mineral suspended solids</td><td>2.611</td><td>2.611</td><td>[mgMSS L<sup>-1</sup>]</td></tr><tr><td>OXY</td><td>Dissolved oxygen</td><td>9.446</td><td>9.446</td><td>[mgO<sub>2</sub> L<sup>-1</sup>]</td></tr><tr><td>TA</td><td>Total alkalinity</td><td>5291</td><td>5291</td><td>[µmol L<sup>-1</sup>]</td></tr><tr><td>DIC</td><td>Dissolved inorganic carbon</td><td>62.728</td><td>62.728</td><td>[mgC L<sup>-1</sup>]</td></tr><tr><td>CO<sub>2</sub>(aq)</td><td>Aqueous carbon dioxide</td><td>0.343</td><td>0.343</td><td>[mgC L<sup>-1</sup>]</td></tr><tr><td>pH</td><td>Acidity</td><td>8.659</td><td>8.695</td><td>[-]</td></tr></tbody></table>

(i) between species communities and (ii) temporally for each species community.

### 3.3.1 River stretch morphology and hydraulic conditions

The stretch measures 10 000 m long and 300 m in width. To simplify the boundary conditions, the upstream inflow and downstream outflow are fixed at 25 m³ s⁻¹, which corresponds to a residence time of 7 d. The water height is fixed at 5 m.

### 3.3.2 Simulation settings and evaluation strategy

The concentrations of all water quality variables of inflow are defined as their initial concentrations in the stretch and remain constant during the simulation (Table 2). Since this paper focuses on the conceptualization of the unified RIVE v1.0 in the water column, no exchange between the benthic layer and water column is considered. The time step of the simulation is 6 min and a simulation period of 365 d is considered. To compare the results of the two digital implementations of unified RIVE v1.0, daily concentrations at 00:00 are plotted. Three statistical criteria (PBIAS: percent bias, %; MAD: mean absolute difference; MaAD: maximum absolute difference) are calculated to evaluate the similarity of the two sets of results. The closer the criteria are to 0, the more similar the concentrations simulated by the two soft-
ware packages are (pyRIVE 1.0 and C-RIVE 0.32).

$$ \text{PBIAS} = 100 \frac{\sum_{i=1}^{N} (C_i - \text{Py}_i)}{\sum_{i=1}^{N} \text{Py}_i} \qquad (59) $$

$$ \text{MAD} = \frac{\sum_{i=1}^{N} |C_i - \text{Py}_i|}{N} \qquad (60) $$

$$ \text{MaAD} = \max(|C_i - \text{Py}_i|) \qquad (61) $$

Here, $C_i$ represents the concentrations simulated by C-RIVE 0.32 (in ANSI C) and $\text{Py}_i$ those simulated by pyRIVE 1.0 (in Python 3). $N$ is the number of values.

### 3.3.3 Simulated concentrations of water quality variables

The concentrations simulated by pyRIVE 1.0 and C-RIVE 0.32 are very similar (and superimposed) for all water quality variables (Fig. 10). A maximum absolute difference (MaAD) of 0.0307 mgO₂ L⁻¹, which is relatively low, is obtained for dissolved oxygen concentration. The mean absolute difference (MAD) for dissolved oxygen concentration is 0.00678 mgO₂ L⁻¹ (Table 3) and the corresponding percent bias (PBIAS) is 0%. The MaAD of 0.0307 mgO₂ L⁻¹ for dissolved oxygen is the cause of the depletion of CYA S (small precursors S of cyanobacteria, Fig. 3) at the beginning of the simulation (not shown here). To correct this depletion of CYA S, the growth of functional macromolecules (CYA F) is reduced according to the availability of CYA S in C-RIVE 0.32. This is why the simulated concentrations of CYA (cyanobacteria) depict a MaAD of 0.0321 mgC L⁻¹ between pyRIVE 1.0 and C-RIVE 0.32. Due to this auto-correction in C-RIVE 0.32, the simulated concentrations of CYA by C-RIVE 0.32 are slightly smaller than those simulated by pyRIVE 1.0 (PBIAS = -1.2%). The values of PBIAS also indicate the similarity between the concentrations simulated by pyRIVE 1.0 and C-RIVE 0.32. Except for CYA, the discrepancies of other variables are extremely low compared to their concentrations (PBIAS ≤ 0.6%). More than half of simulated variables have a PBIAS of 0%.

## 4 Discussion

The results show the ability of the unified RIVE v1.0 to correctly simulate the organic matter degradation and the similarity of its two digital implementations (pyRIVE 1.0 and C-RIVE 0.32). Here, we discuss the biogeochemical cycling simulated by unified RIVE v1.0 in the water column (Sect. 4.1), the model limitations, and the future developments (Sect. 4.3), as well as its benefits for the scientific community (Sect. 4.4).

### 4.1 Biogeochemical cycling in water column simulated by unified RIVE v1.0

The unified RIVE v1.0 simulates the dynamics of microorganisms involving biogeochemical cycling, although the boundary conditions are defined as constant for modeling a river stretch (Fig. 9). Here we interpret the dynamics of diatoms (DIA) and large heterotrophic bacteria (LHB). For this purpose, the budget fluxes of DIA and LHB are calculated.

The decrease in DIA biomass from day 1 to day 15 is related to the low water temperature and low active irradiance, which limit its photosynthesis (Figs. 10, 12). The optimal temperature for the growth of DIA is 21°C, while the water temperature is lower than 3°C (Fig. 12). The low photosynthesis rate leads to negative net production (Fig. 11, green line), which is the difference between the fluxes of photosynthesis and the combined fluxes of respiration, mortality, and excretion. During this period, while the input factors play a positive role, the net change in DIA is still negative (Fig. 11, black line). Over the following days, as the water temperature and active irradiance increase (Fig. 12), the net production shows an increase. However, it still remains negative. The net change in DIA shifts to a positive direction due to a combination of net input and net production, leading to a simulated increase in DIA biomass. This trend continues until day 130 when the maximum DIA biomass is reached (Fig. 10). The decline in DIA biomass simulated from day 130 onwards is due to a combination of factors. Firstly, the input factor could be contributing to the decrease when the DIA biomass exceeds the concentration of DIA in input flow (0.447 mgC L⁻¹, Table 2, Fig. 10). Additionally, the net production rate also plays a role (Fig. 11). Although the photosynthesis rate increases with water temperature and active irradiance until day 179 (not shown here), it is not enough to compensate for the other processes occurring in the diatom population (days 150–170), resulting in an overall decrease in biomass. Despite the positive contributions of net input and net production to DIA biomass around day 175, a significant decrease in biomass occurred due to zooplankton grazing (Fig. 11, red line). Two factors impact the zooplankton dynamics: water temperature and the half-saturation constant of grazing (Eq. 26). The optimal temperature for zooplankton is 25°C and the half-saturation constant of grazing for zooplankton is set to 0.4 mgC L⁻¹. Then, equilibrium of DIA biomass is simulated until day 260 (Fig. 10), which means that the net production and net input in DIA biomass are balanced by the grazing of zooplankton. The input in DIA biomass primarily contributes to the increase in DIA biomass from day 260 (Fig. 11). As the water temperature and active irradiance decrease during this time, the net production of DIA decreases and changes to negative by day 292.

The growth rate of large heterotrophic bacteria (LHB) increases (Fig. 13) with the increase in water temperature, causing a rise in LHB biomass until day 170 (Fig. 10). The fast decrease in small monomeric substrate (SMS) around
<figure><img src="image_12.png" /><figcaption>**Figure 10.** Simulated concentrations of main species by pyRIVE 1.0 and C-RIVE 0.32. See Table 2 for their definitions.</figcaption></figure>

<table><caption>**Table 3.** Statistical criteria for comparing the simulated variables by pyRIVE 1.0 and C-RIVE 0.32 which implement the unified RIVE v1.0. PBIAS: percent bias [%]; MAD: mean absolute difference; MaAD: maximum absolute difference. The unit of MAD and MaAD depends on the species. It can be [mgC L⁻¹], [mgN L⁻¹], [mgP L⁻¹], [mgSi L⁻¹], or [µmol L⁻¹].</caption><thead><tr><th>Species</th><th>PBIAS</th><th>MAD</th><th>MaAD</th><th>Species</th><th>PBIAS</th><th>MAD</th><th>MaAD</th></tr></thead><tbody><tr><td>SHB</td><td>-0.4</td><td>2.07 × 10<sup>-5</sup></td><td>1.40 × 10<sup>-4</sup></td><td>LHB</td><td>-0.4</td><td>1.10 × 10<sup>-4</sup></td><td>4.27 × 10<sup>-4</sup></td></tr><tr><td>AOB</td><td>0</td><td>1.78 × 10<sup>-5</sup></td><td>6.95 × 10<sup>-5</sup></td><td>NOB</td><td>0</td><td>4.85 × 10<sup>-6</sup></td><td>1.98 × 10<sup>-5</sup></td></tr><tr><td>DIA</td><td>-0.1</td><td>1.07 × 10<sup>-3</sup></td><td>3.63 × 10<sup>-3</sup></td><td>GRA</td><td>-0.1</td><td>4.98 × 10<sup>-4</sup></td><td>1.87 × 10<sup>-3</sup></td></tr><tr><td>CYA</td><td>-1.2</td><td>7.26 × 10<sup>-3</sup></td><td>3.15 × 10<sup>-2</sup></td><td>ZOR</td><td>-0.2</td><td>2.26 × 10<sup>-4</sup></td><td>2.89 × 10<sup>-3</sup></td></tr><tr><td>ZOC</td><td>0</td><td>5.25 × 10<sup>-9</sup></td><td>1.69 × 10<sup>-7</sup></td><td>SMS</td><td>-0.6</td><td>6.65 × 10<sup>-4</sup></td><td>4.31 × 10<sup>-3</sup></td></tr><tr><td>DOM<sub>1</sub></td><td>0.1</td><td>6.54 × 10<sup>-5</sup></td><td>3.61 × 10<sup>-4</sup></td><td>DOM<sub>2</sub></td><td>-0.1</td><td>3.94 × 10<sup>-4</sup></td><td>1.34 × 10<sup>-3</sup></td></tr><tr><td>DOM<sub>3</sub></td><td>0</td><td>4.06 × 10<sup>-4</sup></td><td>2.06 × 10<sup>-3</sup></td><td>OXY</td><td>0</td><td>6.78 × 10<sup>-3</sup></td><td>3.07 × 10<sup>-2</sup></td></tr><tr><td>POM<sub>1</sub></td><td>-0.3</td><td>4.68 × 10<sup>-4</sup></td><td>2.27 × 10<sup>-3</sup></td><td>POM<sub>2</sub></td><td>-0.3</td><td>7.80 × 10<sup>-4</sup></td><td>3.97 × 10<sup>-3</sup></td></tr><tr><td>POM<sub>3</sub></td><td>-0.2</td><td>4.06 × 10<sup>-4</sup></td><td>2.06 × 10<sup>-3</sup></td><td>NH4</td><td>0</td><td>4.66 × 10<sup>-5</sup></td><td>5.28 × 10<sup>-4</sup></td></tr><tr><td>NO2</td><td>0</td><td>2.13 × 10<sup>-5</sup></td><td>3.42 × 10<sup>-4</sup></td><td>NO3</td><td>0</td><td>4.04 × 10<sup>-4</sup></td><td>1.90 × 10<sup>-3</sup></td></tr><tr><td>TIP</td><td>0</td><td>1.45 × 10<sup>-4</sup></td><td>3.62 × 10<sup>-4</sup></td><td>DSi</td><td>0</td><td>1.81 × 10<sup>-3</sup></td><td>8.46 × 10<sup>-3</sup></td></tr><tr><td>TA</td><td>0</td><td>4.39 × 10<sup>-2</sup></td><td>2.65 × 10<sup>-1</sup></td><td>pH</td><td>0</td><td>1.14 × 10<sup>-3</sup></td><td>7.22 × 10<sup>-3</sup></td></tr><tr><td>DIC</td><td>0</td><td>8.49 × 10<sup>-3</sup></td><td>6.54 × 10<sup>-2</sup></td><td>CO<sub>2</sub>(aq)</td><td>0</td><td>2.24 × 10<sup>-3</sup></td><td>1.65 × 10<sup>-2</sup></td></tr></tbody></table>

day 175, synchronized with the grazing of zooplankton (Fig. 11), causes a decrease in growth rate of LHB. Its mortality rate is not impacted (not shown here). Consequently, this leads to a significant reduction in LHB biomass around day 175 (Fig. 10). The biomass of LHB remains stable until day 260. After that, it increases in conjunction with the rise in SMS concentration, which is synchronized with the increase in phytoplankton biomass.

## 4.2 Complexity and strengths of the RIVE model

Complexity can be understood in terms of the large number of variables represented and interacting with each other. The RIVE model is a multi-element, multi-form model and the kinetics it represents inevitably incorporate a large number of parameters. This is especially true as the RIVE model has opted for an explicit representation of the living communities (bacteria, phytoplankton, zooplankton, etc.) involved in
<figure><img src="image_13.png" /><figcaption>**Figure 11.** Budget fluxes of DIA (ton C d⁻¹). Grazing: zooplankton grazing fluxes; net production: Eq. (25); net input: input flux-output flux; net change: daily variation of DIA in the river stretch.</figcaption></figure>

<figure><img src="image_14.png" /><figcaption>**Figure 12.** Simulated water temperature, active irradiance, and zoom of active irradiance for days 111–114.</figcaption></figure>

the carbon and nutrient cycles. The model has thus become more complex over time and the addition of new processes (and therefore new parameters) has, as far as possible, been systematically based on experimental work in the laboratory or in the field to reduce the ranges of uncertainty around the kinetic parameters.

The RIVE model is designed as a tool for generating knowledge about the functioning of freshwater ecosystems, and therefore it documents a large number of biogeochemical processes, whether they are expressed weakly or strongly

in a given freshwater ecosystem. The underlying hypothesis is that environmental factors control the intensity with which the various processes involved in the overall functioning of a hydrosystem are expressed.

Nevertheless, some work has specifically focused on analyzing the influence of RIVE parameters, particularly those controlling oxygen levels (Wang et al., 2018). This work identified key physical and physiological parameters. Based on the result of sensitivity analysis, a continuous oxygen data assimilation scheme has been developed (Prose-PA, Wang et al., 2019, 2022). The data assimilation allows determining the physiological properties of microorganisms by integrating the associated uncertainties over time. The recent work of Hasanyar et al. (2023) has also helped to better quantify the sensitivity of oxygen to bacterial kinetics parameters as well as those relating to the composition of organic matter with the aim of parsimonious simplification of the number of parameters.

In these two examples, RIVE (C-RIVE) biogeochemical modeling is implemented in much more complex modeling platforms (particle filter, data assimilation, etc.) and the various analyses (sensitivity, uncertainties, etc.) are also supported by an overall assessment of the performance of the model applied to the Seine River.

## 4.3 Model limitations and future developments of unified RIVE

Currently, the unified RIVE v1.0 presented in this paper describes only the biogeochemical processes in the water column. Comparison of benthic processes and simulations have not been investigated yet. Previous studies showed that sediment plays an important role in the metabolism of river (Vilmin et al., 2016) and lakes (Yan et al., 2022b). A unified sediment module should be further elaborated based on existing modules (Even et al., 2004; Flipo et al., 2004; Thouvenot et al., 2007; Billen et al., 2015; Vilmin et al., 2015b) and implemented into unified RIVE. This sediment module will have to take into account not only the dissolved exchanges between the water column and the sediment but also the resuspension of particulates.

In addition, the unified RIVE v1.0 simulates phytoplankton dynamics, but periphyton or macrophyte development is not implemented in current versions. Flipo et al. (2004) showed that periphyton plays a major role in carbon cycling (primary productivity) in small rivers, not only in the carbon stock fixed at the bottom of the river but also in the carbon enrichment downstream of the river. These limitations should be considered in future developments.

## 4.4 Benefit of the unified RIVE model

The unified RIVE provides a set of governing equations of freshwater biogeochemical processes across different software platforms, such as pyNuts-Riverstrahler (Billen et al.,
<figure><img src="image_15.png" /><figcaption>**Figure 13.** Large heterotrophic bacteria (LHB) dynamics and simulated concentrations of small monomeric substrate (SMS). Net change: daily variation of LHB in river stretch; net input: input flux–output flux.</figcaption></figure>

1994; Garnier et al., 1995; Thieu et al., 2017), PROSE-PA (Wang et al., 2019, 2023a), SWAT-RIVE (Manteaux et al., 2023), QUAL-NET (Minaudo et al., 2018), VEMALA V3 (Korppoo et al., 2017), and Barman (Garnier et al., 2000; Thieu et al., 2006; Yan et al., 2022b), while incorporating the latest developments. The unity of the kinetics is important for facilitating and reinforcing the collaboration nationally or internationally within different research teams. Thanks to the unity property formerly pointed out by the river continuum concept (Vannote et al., 1980), the software packages based on unified RIVE can leverage the already identified parameter values, regardless of the location in the network (Garnier et al., 2020). This feature is of great interest to the different research teams involved in freshwater quality research, for instance river metabolism (Odum, 1956; Garnier and Billen, 2007; Escoffier et al., 2018; Gurung et al., 2019; Rodríguez-Castillo et al., 2019; Garnier et al., 2020; Segatto et al., 2020; Battin et al., 2023) or nutrient cycling (Garnier et al., 1999b; Alexander et al., 2002; Garnier et al., 2002; Billen et al., 2007; Lauerwald et al., 2013; Lindenschmidt et al., 2019; Maavara et al., 2020; Marescaux et al., 2020; Yan et al., 2022a).

Open science has become increasingly popular and even indispensable in the scientific community as it allows for easier

accessibility and the reproduction of scientific results. The unified RIVE project, as an open-source project, allows for the dissemination and wider use of the RIVE biogeochemical model by creating a public repository with different programming languages.

## 5 Conclusions

This paper presents a conceptual freshwater biogeochemistry model: unified RIVE v1.0, programmed in Python 3 and ANSI C. The degradation of organic matter by heterotrophic bacteria and the dynamics of primary producers (phytoplankton) and zooplankton including carbon cycling and nutrient cycling are described exhaustively. In unified RIVE v1.0, organic matter is degraded via bacteria activity, which is simulated by an HSB model. According to the results, the HSB model is able to model the organic matter degradation and bacterial dynamics in batch experiments. A case study is designed to compare the simulations of the two digital implementations (Python 3 for pyRIVE 1.0 and ANSI C for C-RIVE 0.32). These implementations simulate similar concentrations of all state variables including microorganisms, organic carbon, nutrients, and inorganic carbon.
The river stretch case study allows us to compare the two implementations of unified RIVE V1.0 under transient contrasting conditions involving complex biogeochemical cycles. The specific dynamics of each simulated species depend on different limiting factors. The calculation of photosynthesis of phytoplankton (diatoms, Chlorophyceae, cyanobacteria) takes into account the light that naturally presents a day-night variation. The development of diatoms specifically takes into account the dissolved silica in the simulated environment. The growth of microorganisms depends on the quantity of nutrients (primary producer, nitrifying bacteria) and the small monomeric substrate (heterotrophic bacteria). In addition, the effect of water temperature is also taken into account for the physiology of the simulated microorganism communities (photosynthesis, growth, respiration, mortality).

Finally, unified RIVE being an open-source project, contributions from the freshwater biogeochemistry community are strongly encouraged to achieve a better understanding of freshwater ecosystem functioning and further investigate the future of river systems in a changing world.
## Appendix A: Various implementations of the RIVE model and its applications in different freshwater systems

<table><caption>**Table A1.** Various implementations of the RIVE model and its applications in different freshwater systems.</caption><thead><tr><td>Freshwater systems</td><td>Climates</td><td>Software platforms</td><td>References</td></tr></thead><tbody><tr><td>Danube River (Romania and Bulgaria)</td><td>Continental</td><td>Riverstrahler</td><td>Garnier et al. (2002)</td></tr><tr><td>Day-Nhue River (Vietnam)</td><td>Tropical</td><td>Seneque-Riverstrahler</td><td>Luu et al. (2021)</td></tr><tr><td>Grand Morin River (France)</td><td>Temperate</td><td>PROSE</td><td>Flipo et al. (2004, 2007)</td></tr><tr><td>Loire River (France)</td><td>Temperate</td><td>Grafs-Seneque/Riverstrahler</td><td>Garnier et al. (2018a)</td></tr><tr><td>Lot River (France)</td><td>Temperate</td><td>Grafs-Seneque/Riverstrahler</td><td>Garnier et al. (2018b)</td></tr><tr><td>Lule and Kalix rivers (Sweden)</td><td>Subarctic</td><td>Riverstrahler</td><td>Sferratore et al. (2008)</td></tr><tr><td>Mosel River (Germany)</td><td>Temperate</td><td>Riverstrahler</td><td>Garnier et al. (1999a)</td></tr><tr><td>Orgeval watershed (France)</td><td>Temperate</td><td>Seneque-Riverstrahler</td><td>Garnier et al. (2014)</td></tr><tr><td>Red River (China and Vietnam)</td><td>Tropical</td><td>Seneque-Riverstrahler</td><td>Le et al. (2010, 2015); Nguyen et al. (2016)</td></tr><tr><td>Scheldt River (Belgium and Netherlands)</td><td>Temperate</td><td>Seneque-Riverstrahler</td><td>Billen et al. (2005); Thieu et al. (2009)</td></tr><tr><td>Somme River (France)</td><td>Temperate</td><td>Seneque-Riverstrahler</td><td>Thieu et al. (2009, 2010)</td></tr><tr><td>Seine River (France)</td><td>Temperate</td><td>Seneque-Riverstrahler</td><td>Billen et al. (2007); Thieu et al. (2009, 2010); Romero et al. (2019)</td></tr><tr><td>Seine River (France)</td><td>Temperate</td><td>pyNuts-Riverstrahler</td><td>Thieu et al. (2017); Desmit et al. (2018); Raimonet et al. (2018); Marescaux et al. (2020)</td></tr><tr><td>Seine River (France)</td><td>Temperate</td><td>PROSE/PROSE-PA</td><td>Even et al. (1998, 2004, 2007); Raimonet et al. (2015); Vilmin et al. (2015b, a, 2016, 2018); Wang (2019); Wang et al. (2022)</td></tr><tr><td>Zenne River (Belgium)</td><td>Temperate</td><td>Seneque-Riverstrahler</td><td>Garnier et al. (2013)</td></tr><tr><td>Sand-pit lake, reservoirs (France)</td><td>Temperate</td><td>Barman</td><td>Garnier and Billen (1994); Garnier et al. (2000); Thieu et al. (2006); Yan et al. (2022a)</td></tr></tbody></table>
## Appendix B: Parameter values for unified RIVE v1.0

The 120 parameter values necessary for running unified RIVE v1.0 are provided hereafter.

<table><caption>**Table B1.** Parameters related to heterotrophic bacteria.</caption><thead><tr><th>Parameter</th><th>Description</th><th>Value</th><th>Unit</th></tr></thead><tbody><tr><td>e<sup>*</sup><sub>max20 dom1 lhb</sub></td><td>Maximal hydrolysis rate of DOM<sub>1</sub> at 20 °C related to LHB</td><td>0.75</td><td>[h<sup>-1</sup>]</td></tr><tr><td>e<sup>*</sup><sub>max20 dom1 shb</sub></td><td>Maximal hydrolysis rate of DOM<sub>1</sub> at 20 °C related to SHB</td><td>0.75</td><td>[h<sup>-1</sup>]</td></tr><tr><td>e<sup>*</sup><sub>max20 dom2 lhb</sub></td><td>Maximal hydrolysis rate of DOM<sub>2</sub> at 20 °C related to LHB</td><td>0.25</td><td>[h<sup>-1</sup>]</td></tr><tr><td>e<sup>*</sup><sub>max20 dom2 shb</sub></td><td>Maximal hydrolysis rate of DOM<sub>2</sub> at 20 °C related to SHB</td><td>0.25</td><td>[h<sup>-1</sup>]</td></tr><tr><td>b<sup>*</sup><sub>max20 lhb</sub></td><td>Maximal substrate (SMS) uptake rate of LHB at 20 °C</td><td>0.6</td><td>[h<sup>-1</sup>]</td></tr><tr><td>b<sup>*</sup><sub>max20 shb</sub></td><td>Maximal substrate (SMS) uptake rate of SHB at 20 °C</td><td>0.16</td><td>[h<sup>-1</sup>]</td></tr><tr><td>Y<sub>lhb</sub></td><td>Growth yield of LHB</td><td>0.25</td><td>[-]</td></tr><tr><td>Y<sub>shb</sub></td><td>Growth yield of SHB</td><td>0.25</td><td>[-]</td></tr><tr><td>k<sup>*</sup><sub>d20 lhb</sub></td><td>Mortality rate of LHB at 20 °C</td><td>0.05</td><td>[h<sup>-1</sup>]</td></tr><tr><td>k<sup>*</sup><sub>d20 shb</sub></td><td>Mortality rate of SHB at 20 °C</td><td>0.02</td><td>[h<sup>-1</sup>]</td></tr><tr><td>v<sub>shb</sub></td><td>Sinking velocity of LHB</td><td>0.0</td><td>[m h<sup>-1</sup>]</td></tr><tr><td>v<sub>lhb</sub></td><td>Sinking velocity of LHB</td><td>0.02</td><td>[m h<sup>-1</sup>]</td></tr><tr><td>T<sub>opt shb</sub></td><td>Optimal temperature of SHB</td><td>20</td><td>°C</td></tr><tr><td>T<sub>opt lhb</sub></td><td>Optimal temperature of LHB</td><td>22</td><td>°C</td></tr><tr><td>σ<sub>shb</sub></td><td>Range of temperature for SHB</td><td>17</td><td>°C</td></tr><tr><td>σ<sub>lhb</sub></td><td>Range of temperature for LHB</td><td>12</td><td>°C</td></tr><tr><td>K<sub>sms lhb</sub></td><td>Half-saturation constant of LHB for small monomeric substrate</td><td>0.1</td><td>[mgC L<sup>-1</sup>]</td></tr><tr><td>K<sub>sms shb</sub></td><td>Half-saturation constant of SHB for small monomeric substrate</td><td>0.1</td><td>[mgC L<sup>-1</sup>]</td></tr><tr><td>K<sub>dom1 shb</sub></td><td>Half-saturation constant for DOM<sub>1</sub> hydrolysis related to SHB</td><td>0.25</td><td>[mgC L<sup>-1</sup>]</td></tr><tr><td>K<sub>dom1 lhb</sub></td><td>Half-saturation constant for DOM<sub>1</sub> hydrolysis related to LHB</td><td>0.25</td><td>[mgC L<sup>-1</sup>]</td></tr><tr><td>K<sub>dom2 shb</sub></td><td>Half-saturation constant for DOM<sub>2</sub> hydrolysis related to SHB</td><td>2.5</td><td>[mgC L<sup>-1</sup>]</td></tr><tr><td>K<sub>dom2 lhb</sub></td><td>Half-saturation constant for DOM<sub>2</sub> hydrolysis related to LHB</td><td>2.5</td><td>[mgC L<sup>-1</sup>]</td></tr></tbody></table>

<table><caption>**Table B2.** Parameters related to nitrifying bacteria.</caption><thead><tr><th>Parameter</th><th>Description</th><th>Value</th><th>Unit</th></tr></thead><tbody><tr><td>μ<sup>*</sup><sub>max20 aob</sub></td><td>Maximal growth rate of AOB at 20 °C</td><td>0.07</td><td>[h<sup>-1</sup>]</td></tr><tr><td>μ<sup>*</sup><sub>max20 nob</sub></td><td>Maximal growth rate of NOB at 20 °C</td><td>0.05</td><td>[h<sup>-1</sup>]</td></tr><tr><td>K<sub>O2 aob</sub></td><td>Half-saturation constant of AOB for O<sub>2</sub></td><td>0.64</td><td>[mgO<sub>2</sub> L<sup>-1</sup>]</td></tr><tr><td>K<sub>O2 nob</sub></td><td>Half-saturation constant of NOB for O<sub>2</sub></td><td>1.088</td><td>[mgO<sub>2</sub> L<sup>-1</sup>]</td></tr><tr><td>K<sub>nh4 aob</sub></td><td>Half-saturation constant of AOB for NH<sub>4</sub><sup>+</sup></td><td>0.75</td><td>[mgN L<sup>-1</sup>]</td></tr><tr><td>K<sub>NO2 nob</sub></td><td>Half-saturation constant of NOB for NO<sub>2</sub></td><td>0.05</td><td>[mgN L<sup>-1</sup>]</td></tr><tr><td>Y<sub>aob</sub></td><td>Growth yield of AOB</td><td>0.07</td><td>[mgC mgN<sup>-1</sup>]</td></tr><tr><td>Y<sub>nob</sub></td><td>Growth yield of AOB</td><td>0.02</td><td>[mgC mgN<sup>-1</sup>]</td></tr><tr><td>k<sup>*</sup><sub>d20 aob</sub></td><td>Mortality rate of AOB at 20 °C</td><td>0.005</td><td>[h<sup>-1</sup>]</td></tr><tr><td>k<sup>*</sup><sub>d20 nob</sub></td><td>Mortality rate of NOB at 20 °C</td><td>0.005</td><td>[h<sup>-1</sup>]</td></tr><tr><td>v<sub>aob</sub></td><td>Sinking velocity of AOB</td><td>0.005</td><td>[m h<sup>-1</sup>]</td></tr><tr><td>v<sub>nob</sub></td><td>Sinking velocity of NOB</td><td>0.005</td><td>[m h<sup>-1</sup>]</td></tr><tr><td>T<sub>opt aob</sub></td><td>Optimal temperature of AOB</td><td>23</td><td>[°C]</td></tr><tr><td>σ<sub>aob</sub></td><td>Range of temperature for NOB</td><td>18</td><td>[°C]</td></tr><tr><td>T<sub>opt nob</sub></td><td>Optimal temperature of NOB</td><td>23</td><td>[°C]</td></tr><tr><td>σ<sub>nob</sub></td><td>Range of temperature for NOB</td><td>18</td><td>[°C]</td></tr></tbody></table>

<footnote>* Parameters depend on water temperature and are multiplied by $f(T) = \frac{(T-T_{opt})^2}{\frac{e}{\sigma^2}}$, where T is water temperature in °C.</footnote>

<footnote>* Parameters depend on water temperature and are multiplied by $f(T) = \frac{(T-T_{opt})^2}{\frac{e}{\sigma^2}}$, where T is water temperature in °C.</footnote>
<table><caption>Table B3. Parameters related to primary producer dynamics.</caption><thead><tr><th>Parameter</th><th>Description</th><th>Value</th><th>Unit</th></tr></thead><tbody><tr><td>P<sup>*</sup><sub>max20,dia</sub></td><td>Maximal photosynthesis rate of diatoms at 20 °C</td><td>0.2</td><td>[h<sup>-1</sup>]</td></tr><tr><td>P<sup>*</sup><sub>max20,gra</sub></td><td>Maximal photosynthesis rate of green algae at 20 °C</td><td>0.25</td><td>[h<sup>-1</sup>]</td></tr><tr><td>P<sup>*</sup><sub>max20,cya</sub></td><td>Maximal photosynthesis rate of cyanobacteria at 20 °C</td><td>0.1</td><td>[h<sup>-1</sup>]</td></tr><tr><td>α<sub>dia</sub></td><td>Photosynthetic efficiency of diatoms</td><td>0.0012</td><td>[h<sup>-1</sup> (µE m<sup>-2</sup> s<sup>-1</sup>)<sup>-1</sup>]</td></tr><tr><td>α<sub>gra</sub></td><td>Photosynthetic efficiency of green algae</td><td>0.0012</td><td>[h<sup>-1</sup> (µE m<sup>-2</sup> s<sup>-1</sup>)<sup>-1</sup>]</td></tr><tr><td>α<sub>cya</sub></td><td>Photosynthetic efficiency of cyanobacteria</td><td>0.0012</td><td>[h<sup>-1</sup> (µE m<sup>-2</sup> s<sup>-1</sup>)<sup>-1</sup>]</td></tr><tr><td>β<sub>dia</sub></td><td>Photoinhibition capacity of diatoms</td><td>0.0</td><td>[h<sup>-1</sup> (µE m<sup>-2</sup> s<sup>-1</sup>)<sup>-1</sup>]</td></tr><tr><td>β<sub>gra</sub></td><td>Photoinhibition capacity of green algae</td><td>0.0</td><td>[h<sup>-1</sup> (µE m<sup>-2</sup> s<sup>-1</sup>)<sup>-1</sup>]</td></tr><tr><td>β<sub>cya</sub></td><td>Photoinhibition capacity of cyanobacteria</td><td>0.0</td><td>[h<sup>-1</sup> (µE m<sup>-2</sup> s<sup>-1</sup>)<sup>-1</sup>]</td></tr><tr><td>η<sub>base</sub></td><td>Light-extinction-related coefficient for pure water</td><td>0.2</td><td>[m<sup>-1</sup>]</td></tr><tr><td>η<sub>chl a</sub></td><td>Light algal self-shading light extinction coefficient</td><td>0.02</td><td>[m<sup>-1</sup> (µg chl a L<sup>-1</sup>)<sup>-1</sup>]</td></tr><tr><td>η<sub>ss</sub></td><td>Light extinction coefficient related to suspended solid</td><td>0.042</td><td>[m<sup>-1</sup> (mg L<sup>-1</sup>)<sup>-1</sup>]</td></tr><tr><td>μ<sup>*</sup><sub>max20,dia</sub></td><td>Maximal growth rate of diatoms at 20 °C</td><td>0.05</td><td>[h<sup>-1</sup>]</td></tr><tr><td>μ<sup>*</sup><sub>max20,gra</sub></td><td>Maximal growth rate of green algae at 20 °C</td><td>0.05</td><td>[h<sup>-1</sup>]</td></tr><tr><td>μ<sup>*</sup><sub>max20,cya</sub></td><td>Maximal growth rate of cyanobacteria at 20 °C</td><td>0.025</td><td>[h<sup>-1</sup>]</td></tr><tr><td>K<sub>S,dia</sub></td><td>Half-saturation constant for small precursors of diatoms</td><td>0.06</td><td>[-]</td></tr><tr><td>K<sub>S,gra</sub></td><td>Half-saturation constant for small precursors of green algae</td><td>0.06</td><td>[-]</td></tr><tr><td>K<sub>S,cya</sub></td><td>Half-saturation constant for small precursors of cyanobacteria</td><td>0.06</td><td>[-]</td></tr><tr><td>K<sub>N,dia</sub></td><td>Half-saturation constant for nitrogen of diatoms</td><td>0.014</td><td>[mg N L<sup>-1</sup>]</td></tr><tr><td>K<sub>N,gra</sub></td><td>Half-saturation constant for nitrogen of green algae</td><td>0.014</td><td>[mg N L<sup>-1</sup>]</td></tr><tr><td>K<sub>N,cya</sub></td><td>Half-saturation constant for nitrogen of cyanobacteria</td><td>0.014</td><td>[mg N L<sup>-1</sup>]</td></tr><tr><td>K<sub>P,dia</sub></td><td>Half-saturation constant for phosphorus of diatoms</td><td>0.0155</td><td>[mg P L<sup>-1</sup>]</td></tr><tr><td>K<sub>P,gra</sub></td><td>Half-saturation constant for phosphorus of green algae</td><td>0.062</td><td>[mg P L<sup>-1</sup>]</td></tr><tr><td>K<sub>P,cya</sub></td><td>Half-saturation constant for phosphorus of cyanobacteria</td><td>0.062</td><td>[mg P I<sup>-1</sup>]</td></tr><tr><td>K<sub>Si,dia</sub></td><td>Half-saturation constant for silica of diatoms</td><td>0.196</td><td>[mg Si L<sup>-1</sup>]</td></tr><tr><td>R<sup>*</sup><sub>m20,dia</sub></td><td>Maintenance respiration coefficient of diatoms at 20 °C</td><td>0.002</td><td>[h<sup>-1</sup>]</td></tr><tr><td>R<sup>*</sup><sub>m20,gra</sub></td><td>Maintenance respiration coefficient of green algae at 20 °C</td><td>0.002</td><td>[h<sup>-1</sup>]</td></tr><tr><td>R<sup>*</sup><sub>m20,cya</sub></td><td>Maintenance respiration coefficient of cyanobacteria at 20 °C</td><td>0.002</td><td>[h<sup>-1</sup>]</td></tr><tr><td>R<sub>μ,dia</sub></td><td>Energetic cost of growth of diatoms</td><td>0.5</td><td>[-]</td></tr><tr><td>R<sub>μ,gra</sub></td><td>Energetic cost of growth of green algae</td><td>0.5</td><td>[-]</td></tr><tr><td>R<sub>μ,cya</sub></td><td>Energetic cost of growth of cyanobacteria</td><td>0.5</td><td>[-]</td></tr><tr><td>E<sub>cst,dia</sub></td><td>Basic excretion rate of diatoms</td><td>0.006</td><td>[h<sup>-1</sup>]</td></tr><tr><td>E<sub>cst,gra</sub></td><td>Basic excretion rate of green algae</td><td>0.006</td><td>[h<sup>-1</sup>]</td></tr><tr><td>E<sub>cst,cya</sub></td><td>Basic excretion rate of cyanobacteria</td><td>0.006</td><td>[h<sup>-1</sup>]</td></tr><tr><td>E<sub>phot,dia</sub></td><td>Excretion constant of diatoms related to photosynthesis</td><td>0.001</td><td>[-]</td></tr><tr><td>E<sub>phot,gra</sub></td><td>Excretion constant of green algae related to photosynthesis</td><td>0.001</td><td>[-]</td></tr><tr><td>E<sub>phot,cya</sub></td><td>Excretion constant of cyanobacteria related to photosynthesis</td><td>0.001</td><td>[-]</td></tr><tr><td>S<sup>*</sup><sub>R,max20,dia</sub></td><td>Maximal rate of reserve product synthesis for diatoms at 20 °C</td><td>0.15</td><td>[h<sup>-1</sup>]</td></tr><tr><td>S<sup>*</sup><sub>R,max20,gra</sub></td><td>Maximal rate of reserve product synthesis for green algae at 20 °C</td><td>0.2</td><td>[h<sup>-1</sup>]</td></tr><tr><td>S<sup>*</sup><sub>R,max20,cya</sub></td><td>Maximal rate of reserve product synthesis for cyanobacteria at 20 °C</td><td>0.075</td><td>[h<sup>-1</sup>]</td></tr><tr><td>C<sup>*</sup><sub>R,max20,dia</sub></td><td>Maximal rate of reserve product catabolism for diatoms at 20 °C</td><td>0.2</td><td>[h<sup>-1</sup>]</td></tr><tr><td>C<sup>*</sup><sub>R,max20,gra</sub></td><td>Maximal rate of reserve product catabolism for green algae at 20 °C</td><td>0.2</td><td>[h<sup>-1</sup>]</td></tr><tr><td>C<sup>*</sup><sub>R,max20,cya</sub></td><td>Maximal rate of reserve product catabolism for cyanobacteria at 20 °C</td><td>0.2</td><td>[h<sup>-1</sup>]</td></tr><tr><td>k<sup>*</sup><sub>d20,dia</sub></td><td>Rate of diatom mortality at 20 °C</td><td>0.025</td><td>[h<sup>-1</sup>]</td></tr><tr><td>k<sup>*</sup><sub>d20,gra</sub></td><td>Rate of green algae mortality at 20 °C</td><td>0.025</td><td>[h<sup>-1</sup>]</td></tr><tr><td>k<sup>*</sup><sub>d20,cya</sub></td><td>Rate of cyanobacteria mortality at 20 °C</td><td>0.015</td><td>[h<sup>-1</sup>]</td></tr><tr><td>vs<sub>dia</sub></td><td>Sinking velocity of diatoms</td><td>0.006</td><td>[m h<sup>-1</sup>]</td></tr><tr><td>vs<sub>gra</sub></td><td>Sinking velocity of green algae</td><td>0.001</td><td>[m h<sup>-1</sup>]</td></tr><tr><td>vs<sub>cya</sub></td><td>Sinking velocity of cyanobacteria</td><td>0.006</td><td>[m h<sup>-1</sup>]</td></tr><tr><td>T<sub>opt,dia</sub></td><td>Optimal temperature of diatoms</td><td>21</td><td>[°C]</td></tr><tr><td>T<sub>opt,gra</sub></td><td>Optimal temperature of green algae</td><td>37</td><td>[°C]</td></tr><tr><td>T<sub>opt,cya</sub></td><td>Optimal temperature of cyanobacteria</td><td>37</td><td>[°C]</td></tr><tr><td>σ<sub>dia</sub></td><td>Range of temperature for diatoms</td><td>13</td><td>[°C]</td></tr><tr><td>σ<sub>gra</sub></td><td>Range of temperature for green algae</td><td>15</td><td>[°C]</td></tr><tr><td>σ<sub>cya</sub></td><td>Range of temperature for cyanobacteria</td><td>12</td><td>[°C]</td></tr></tbody></table>

<footnote>* Parameters depend on water temperature and are multiplied by $f(T) = \frac{e^{-\frac{(T-T_{opt})^2}{\sigma^2}}}{e^{-\frac{(20-T_{opt})^2}{\sigma^2}}}$, where T is water temperature in °C.</footnote>
<table><caption>**Table B4.** Parameters for organic matter dynamics.</caption><thead><tr><th>Parameter</th><th>Description</th><th>Value</th><th>Unit</th></tr></thead><tbody><tr><td>ε<sub>dom1</sub></td><td>DOM<sub>1</sub> fraction in lysis products</td><td>0.2</td><td>[-]</td></tr><tr><td>ε<sub>dom2</sub></td><td>DOM<sub>2</sub> fraction in lysis products</td><td>0.2</td><td>[-]</td></tr><tr><td>ε<sub>dom3</sub></td><td>DOM<sub>3</sub> fraction in lysis products</td><td>0.1</td><td>[-]</td></tr><tr><td>ε<sub>pom1</sub></td><td>POM<sub>1</sub> fraction in lysis products</td><td>0.2</td><td>[-]</td></tr><tr><td>ε<sub>pom2</sub></td><td>POM<sub>2</sub> fraction in lysis products</td><td>0.2</td><td>[-]</td></tr><tr><td>ε<sub>pom3</sub></td><td>POM<sub>3</sub> fraction in lysis products</td><td>0.1</td><td>[-]</td></tr><tr><td>k<sup>*</sup><sub>pom1,20</sub></td><td>POM<sub>1</sub> hydrolysis rate constant at 20 °C</td><td>0.005</td><td>[h<sup>-1</sup>]</td></tr><tr><td>k<sub>pom2</sub></td><td>POM<sub>2</sub> hydrolysis rate constant</td><td>0.00025</td><td>[h<sup>-1</sup>]</td></tr></tbody></table>

* Parameters depend on water temperature and are multiplied by $f(T) = \frac{e^{\frac{(T-T_{opt})^2}{\sigma^2}}}{e^{\frac{(20-T_{opt})^2}{\sigma^2}}}$, where T is water temperature in °C.

<table><caption>**Table B5.** Zooplankton parameters.</caption><thead><tr><th>Parameter</th><th>Description</th><th>Value</th><th>Unit</th></tr></thead><tbody><tr><td>μ<sup>*</sup><sub>max20,zor</sub></td><td>Maximal growth rate of ZOR at 20 °C</td><td>0.025</td><td>[h<sup>-1</sup>]</td></tr><tr><td>μ<sup>*</sup><sub>max20,zoc</sub></td><td>Maximal growth rate of ZOC at 20 °C</td><td>0.015</td><td>[h<sup>-1</sup>]</td></tr><tr><td>b<sup>*</sup><sub>max20,zor</sub></td><td>Maximal grazing rate of ZOR at 20 °C</td><td>0.1</td><td>[h<sup>-1</sup>]</td></tr><tr><td>b<sup>*</sup><sub>max20,zoc</sub></td><td>Maximal grazing rate of ZOC at 20 °C</td><td>0.05</td><td>[h<sup>-1</sup>]</td></tr><tr><td>K<sub>phy,zor</sub></td><td>Half-saturation constant for grazing phytoplankton of ZOR</td><td>0.1</td><td>[mgC L<sup>-1</sup>]</td></tr><tr><td>K<sub>phy,zoc</sub></td><td>Half-saturation constant for grazing phytoplankton of ZOC</td><td>0.1</td><td>[mgC L<sup>-1</sup>]</td></tr><tr><td>PHY<sub>0,zor</sub></td><td>Threshold phytoplankton concentration for grazing of ZOR</td><td>0.1</td><td>[mgC L<sup>-1</sup>]</td></tr><tr><td>PHY<sub>0,zoc</sub></td><td>Threshold phytoplankton concentration for grazing of ZOC</td><td>0.1</td><td>[mgC L<sup>-1</sup>]</td></tr><tr><td>*k<sub>d20,zor</sub></td><td>Mortality rate of ZOR at 20 °C</td><td>0.007</td><td>[h<sup>-1</sup>]</td></tr><tr><td>*k<sub>d20,zoc</sub></td><td>Mortality rate of ZOC at 20 °C</td><td>0.007</td><td>[h<sup>-1</sup>]</td></tr><tr><td>T<sub>opt,zor</sub></td><td>Optimal temperature of ZOR</td><td>25</td><td>[°C]</td></tr><tr><td>T<sub>opt,zoc</sub></td><td>Optimal temperature of ZOC</td><td>25</td><td>[°C]</td></tr><tr><td>σ<sub>zor</sub></td><td>Range of temperature for ZOR</td><td>10</td><td>[°C]</td></tr><tr><td>σ<sub>zoc</sub></td><td>Range of temperature for ZOC</td><td>10</td><td>[°C]</td></tr><tr><td>vs<sub>zor</sub></td><td>Sinking velocity of ZOR</td><td>0.02</td><td>[m h<sup>-1</sup>]</td></tr><tr><td>vs<sub>zoc</sub></td><td>Sinking velocity of ZOC</td><td>0.02</td><td>[m h<sup>-1</sup>]</td></tr></tbody></table>

* Parameters depend on water temperature and are multiplied by $f(T) = \frac{e^{\frac{(T-T_{opt})^2}{\sigma^2}}}{e^{\frac{(20-T_{opt})^2}{\sigma^2}}}$, where T is water temperature in °C.

<table><caption>**Table B6.** Phosphate- and silica-related parameters.</caption><thead><tr><th>Parameter</th><th>Description</th><th>Value</th><th>Unit</th></tr></thead><tbody><tr><td colspan="4" style="text-align: center;">Phosphate adsorption desorption</td></tr><tr><td>P<sub>ac</sub></td><td>Maximal adsorption capacity of mineral suspended solids (MSSs)</td><td>0.00558</td><td>[mgP mgMSS<sup>-1</sup>]</td></tr><tr><td>K<sub>ps</sub></td><td>Half-saturation adsorption constant</td><td>0.682</td><td>[mgP L<sup>-1</sup>]</td></tr><tr><td colspan="4" style="text-align: center;">Silica dynamics</td></tr><tr><td>K<sub>bSi20</sub></td><td>Biogenic silica dissolution rate at 20 °C</td><td>0.0001</td><td>[h<sup>-1</sup>]</td></tr></tbody></table>
*Code and data availability.* C-RIVE 0.32 implements the unified RIVE v1.0 in ANSI C. It is available under Eclipse Public License 2.0 in the following Zenodo repository: https://doi.org/10.5281/zenodo.7849609 (Wang et al., 2023b). pyRIVE 1.0 implements the unified RIVE v1.0 in Python 3 and is available under Eclipse Public License 2.0 in the InDoRES repository: https://doi.org/10.48579/PRO/Z9ACP1 (Thieu et al., 2023). The dataset used in this paper is available in the following Zenodo repository: https://doi.org/10.5281/zenodo.10490669 (Wang et al., 2024).

*Author contributions.* SW: conceptualization, methodology, software, validation, formal analysis, investigation, visualization, writing (original draft). VT: conceptualization, methodology, software, formal analysis, writing (review and editing), supervision. GB: conceptualization, methodology, data curation, formal analysis, writing (review and editing). JG: conceptualization, data curation, formal analysis, writing (review and editing). MS: software, writing (review and editing). AM: conceptualization, software. XY: software. NF: conceptualization, methodology, formal analysis, software, writing (review and editing), supervision, funding acquisition.

*Competing interests.* The contact author has declared that none of the authors has any competing interests.

*Disclaimer.* Publisher's note: Copernicus Publications remains neutral with regard to jurisdictional claims made in the text, published maps, institutional affiliations, or any other geographical representation in this paper. While Copernicus Publications makes every effort to include appropriate place names, the final responsibility lies with the authors.

*Acknowledgements.* This work was carried out under the PIREN-Seine program (https://www.piren-seine.fr/, last access: 20 December 2023), part of French Long-Term Socio-Ecological Research (LTSER) (in French: "Zone Atelier Seine" - INEE, CNRS). We thank Lou Weidenfeld, Benjamin Mercier, and Anun Marninez for carrying out organic matter batch experiments.

Particularly, we would like to dedicate this paper to the memory of Pierre Servais: his knowledge, expertise, and passion for this field of study will continue to guide us in the pursuit of our scientific work.

*Review statement.* This paper was edited by Vassilios Vervatis and reviewed by two anonymous referees.

## References

Aissa-Grouz, N., Garnier, J., and Billen, G.: Long trend reduction of phosphorus wastewater loading in the Seine: determination of phosphorus speciation and sorption for model-

ing algal growth, Environ. Sci. Pollut. Res., 25, 23515-23528, https://doi.org/10.1007/s11356-016-7555-7, 2016.

Alexander, R. B., Elliott, A. H., Shankar, U., and McBride, G. B.: Estimating the sources and transport of nutrients in the Waikato River Basin, New Zealand, Water Resour. Res., 38, 4-1-4-23, https://doi.org/10.1029/2001WR000878, 2002.

Azam, F., Fenchel, T., Gray, J., Meyer, L., and Thingstad, F.: The Ecological Role of Water-Column Microbes in the Sea, Mar. Ecol. Prog. Ser., 10, 257-263, 1983.

Battin, T. J., Lauerwald, R., Bernhardt, E. S., Bertuzzo, E., Gener, L. G., Hall, R. O., Hotchkiss, E. R., Maavara, T., Pavelsky, T. M., Ran, L., Raymond, P., Rosentreter, J. A., and Regnier, P.: River Ecosystem Metabolism and Carbon Biogeochemistry in a Changing World, Nature, 613, 449-459, https://doi.org/10.1038/s41586-022-05500-8, 2023.

Billen, G.: Protein degradation in Aquatic Environments, in: Microbial Enzyme in Aquatic Environments, edited by: Chrost, R., Springer Verlag, Berlin, https://doi.org/10.1007/978-1-4612-3090-8_7, 1991.

Billen, G. and Servais, P.: Modélisation des processus de dégradation bactérienne de la matière organique en milieu aquatique, in: Micro-organismes dans les écosystèmes océaniques, edited by: Bianchi, M. et al., 219-245, Masson Paris, ISBN 9782225815522, 1989.

Billen, G., Servais, P., and Fontigny, A.: Growth and mortality in bacterial populations dynamics of aquatic environments, Arch. Hydrobiol. Beih. Ergebn. Limnol., 31, 173-183, 1988.

Billen, G., Servais, P., and Becquevort, S.: Dynamics of bacterioplankton in oligotrophic and eutrophic aquatic environments: bottom-up or top-down control?, in: Fluxes Between Trophic Levels and Through the Water-Sediments Interface, edited by: Bonin, D. and Golterman, H., 37-42, Kluwer Academic Publishers, https://doi.org/10.1007/BF00041438, 1990.

Billen, G., Garnier, J., and Hanset, P.: Modelling phytoplankton development in whole drainage networks: the RIVERSTRAHLER Model applied to the Seine river system, Hydrobiologia, 289, 119-137, 1994.

Billen, G., Garnier, J., and Rousseau, V.: Nutrient fluxes and water quality in the drainage network of the Scheldt basin over the last 50 year, Hydrobiologia, 540, 47-67, 2005.

Billen, G., Garnier, J., Mouchel, J.-M., and Silvestre, M.: The Seine system: Introduction to a multidisciplinary approach of the functioning of a regional river system, Sci. Total Environ., 375, 1-12, 2007.

Billen, G., Garnier, J., and Silvestre, M.: A simplified algorithm for calculating benthic nutrient fluxes in river systems, Ann. Limnol.-Int. J. Lim., 51, 37-47, https://doi.org/10.1051/limn/2014030, 2015.

Brion, N. and Billen, G.: Une réevaluation de la méthode de mesure de l'activité nitrifiante autotrophe par la méthode d'incorporation de bicarbonate marqué $H^{14}CO_3$ et son application pour estimer des biomasses de bactéries nitrifiantes, Revue des Sciences de l'Eau, 11, 283-302, 1998.

Conley, D. J., Kilham, S. S., and Theriot, E.: Differences in silica content between marine and freshwater diatoms, Limnol. Oceanogr., 34, 205-212, https://doi.org/10.4319/lo.1989.34.1.0205, 1989.
Culberson, C. H.: Calculation of the in situ pH of seawater1, Limnol. Oceanogr., 25, 150–152, https://doi.org/10.4319/lo.1980.25.1.0150, 1980.

Desmit, X., Thieu, V., Billen, G., Campuzano, F., Dulière, V., Garnier, J., Lassaletta, L., Ménesguen, A., Neves, R., Pinto, L., Silvestre, M., Sobrinho, J., and Lacroix, G.: Reducing marine eutrophication may require a paradigmatic change, Sci. Total Environ., 635, 1444–1466, https://doi.org/10.1016/j.scitotenv.2018.04.181, 2018.

Escoffier, N., Bensoussan, N., Vilmin, L., Flipo, N., Rocher, V., David, A., Métivier, F., and Groleau, A.: Estimating ecosystem metabolism from continuous multi-sensor measurements in the Seine River, Environ. Sci. Pollut. Res., 25, 23451–23467, https://doi.org/10.1007/s11356-016-7096-0, 2018.

Even, S., Poulin, M., Garnier, J., Billen, G., Servais, P., Chesterikoff, A., and Coste, M.: River ecosystem modelling: Application of the PROSE model to the Seine river (France), Hydrobiologia, 373, 27–45, https://doi.org/10.1023/A:1017045522336, 1998.

Even, S., Poulin, M., Mouchel, J. M., Seidl, M., and Servais, P.: Modelling oxygen deficits in the Seine river downstream of combined sewer overflows, Ecol. Model., 173, 177–196, 2004.

Even, S., Mouchel, J. M., Servais, P., Flipo, N., Poulin, M., Blanc, S., Chabanel, M., and Paffoni, C.: Modeling the impacts of Combined Sewer Overflows on the river Seine water quality, Sci. Total Environ., 375, 140–151, https://doi.org/10.1016/j.scitotenv.2006.12.007, 2007.

Fick, A.: Ueber diffusion, Annalen der Physik und Chemie, J.A. Barth, https://doi.org/10.1002/andp.18551700105, 1855.

Flipo, N., Even, S., Poulin, M., Tusseau-Vuillemin, M. H., Améziane, T., and Dauta, A.: Biogeochemical Modelling at the River Scale: Plankton and Periphyton Dynamics – Grand Morin case study, France, Ecol. Model., 176, 333–347, 2004.

Flipo, N., Rabouille, C., Poulin, M., Even, S., Tusseau-Vuillemin, M. H., and Lalande, M.: Primary production in headwater streams of the Seine basin: the Grand Morin case study, Sci. Total Environ., 375, 98–109, https://doi.org/10.1016/j.scitotenv.2006.12.015, 2007.

Fontigny, A., Billen, G., and Vives-Rego, J.: Some kinetic characteristics of exoproteolytic activity in coastal seawater, Estuar. Coast. Shelf Sci., 25, 127–133, https://doi.org/10.1016/0272-7714(87)90030-8, 1987.

Fuhrman, J. and Azam, F.: Thymidine Incorporation as a Measure of Heterotrophic Bacterioplankton Production in Marine Surface Waters: Evaluation and Field Results, Mar. Biol., 66, 109–120, 1982.

Garnier, J. and Billen, G.: Ecological interactions in a shallow sand-pit lake (Lake Créteil, Parisian Basin, France): a modelling approach, Hydrobiologia, 275/276, 97–114, 1994.

Garnier, J. and Billen, G.: Production vs. Respiration in river systems: An indicator of an “ecological status”, Sci. Total Environ., 375, 110–124, https://doi.org/10.1016/j.scitotenv.2006.12.006, 2007.

Garnier, J., Billen, G., and Servais, P.: Physiological characteristics and ecological role of small- and large-sized bacteria in a polluted river (Seine River, France), Arch. Hydrobiol. Beih., 37, 83–94, 1992a.

Garnier, J., Servais, P., and Billen, G.: Bacterioplankton in the Seine river (France): impact of the Parisian urban effluent, Can. J. Microbiol., 38, 56–64, 1992b.

Garnier, J., Billen, G., and Coste, M.: Seasonal succession of diatoms and chlorophycae in the drainage network of the river Seine: Observations and modelling, Limnol. Oceanogr., 40, 750–765, 1995.

Garnier, J., Billen, G., Hanset, P., Testard, P., and Coste, M.: Développement algual et eutrophisation dans le réseau hydrographique de la Seine, in: La Seine en son bassin-Fonctionnement écologique d'un système fluvial anthropisé, edited by: Meybeck, M., de Marsily, G., and Fustec, E., 593–626, Elsevier, ISBN 9782842990589, 1998.

Garnier, J., Billen, G., and Palfner, L.: Understanding the Oxygen Budget and Related Ecological Processes in the River Mosel: The RIVERSTRAHLER Approach, in: Man and River Systems, edited by: Garnier, J. and Mouchel, J.-M., Springer Netherlands, Dordrecht, 151–166, ISBN 978-90-481-5393-0 978-94-017-2163-9, https://doi.org/10.1007/978-94-017-2163-9_17, 1999a.

Garnier, J., Leporcq, B., Sanchez, N., and Philippon, X.: Biogeochemical mass-balances (C, N, P, Si) in three large reservoirs of the Seine Basin (France), Biogeochemistry, 47, 119–146, https://doi.org/10.1023/A:1006101318417, 1999b.

Garnier, J., Billen, G., Sanchez, N., and Leporcq, B.: Ecological functioning of the Marne Reservoir (upper Seine basin, France), Regul. River., 16, 51–71, 2000.

Garnier, J., Billen, G., Hannon, E., Fonbonne, S., Videnina, Y., and Soulie, M.: Modeling transfer and retention of nutrients in the drainage network of the Danube River, Estuar. Coast. Shelf, 54, 285–308, 2002.

Garnier, J., Brion, N., Callens, J., Passy, P., Deligne, C., Billen, G., Servais, P., and Billen, C.: Modeling historical changes in nutrient delivery and water quality of the Zenne River (1790s–2010): The role of land use, waterscape and urban wastewater management, J. Marine Syst., 128, 62–76, https://doi.org/10.1016/j.jmarsys.2012.04.001, 2013.

Garnier, J., Billen, G., Vilain, G., Benoit, M., Passy, P., Tallec, G., Tournebize, J., Anglade, J., Billy, C., Mercier, B., Ansart, P., Azougui, A., Sebilo, M., and Kao, C.: Curative vs. preventive management of nitrogen transfers in rural areas: Lessons from the case of the Orgeval watershed (Seine River basin, France), J. Environ. Manage., 144, 125–134, https://doi.org/10.1016/j.jenvman.2014.04.030, 2014.

Garnier, J., Ramarson, A., Billen, G., Théry, S., Thiéry, D., Thieu, V., Minaudo, C., and Moatar, F.: Nutrient inputs and hydrology together determine biogeochemical status of the Loire River (France): Current situation and possible future scenarios, Sci. Total Environ., 637–638, 609–624, https://doi.org/10.1016/j.scitotenv.2018.05.045, 2018a.

Garnier, J., Ramarson, A., Thieu, V., Némery, J., Théry, S., Billen, G., and Coynel, A.: How can water quality be improved when the urban waste water directive has been fulfilled? A case study of the Lot river (France), Environ. Sci. Pollut. Res., 25, 11924–11939, https://doi.org/10.1007/s11356-018-1428-1, 2018b.

Garnier, J., Marescaux, A., Guillon, S., Vilmin, L., Rocher, V., Billen, G., Thieu, V., Silvestre, M., Passy, P., Raimonet, M., Groleau, A., Théry, S., Tallec, G., and Flipo, N.: The Handbook of Environmental Chemistry, chap. Ecological Functioning of the Seine River: From Long-Term Modelling Ap
proaches to High-Frequency Data Analysis, Handbook of En-
vironmental Chemistry, Springer, Berlin, Heidelberg, 1–28,
https://doi.org/10.1007/698_2019_379, 2020.

Garnier, J., Weidenfeld, L., Billen, G., Martinez, A., Mercier,
B., Rocher, V., Tabuchi, J.-P., and Azimi, S.: La matière
organique dans le continuum terrestre-aquatique du bassin
de la Seine, Rapport annuel PIREN-Seine, PIREN-Seine,
https://doi.org/10.26047/PIREN.rapp.ann.2021.vol21, 2021.

Gurung, A., Iwata, T., Nakano, D., and Urabe, J.: River Metabolism along a Latitudinal Gradient across Japan and in a global scale, Sci. Rep., 9, 4932, https://doi.org/10.1038/s41598-019-41427-3, 2019.

Gypens, N., Lancelot, C., and Borges, A. V.: Carbon dynamics and CO₂ air-sea exchanges in the eutrophied coastal waters of the Southern Bight of the North Sea: a modelling study, Biogeosciences, 1, 147–157, https://doi.org/10.5194/bg-1-147-2004, 2004.

Hasanyar, M., Romary, T., Wang, S., and Flipo, N.: How much do bacterial growth properties and biodegradable dissolved organic matter control water quality at low flow?, Biogeosciences, 20, 1621–1633, https://doi.org/10.5194/bg-20-1621-2023, 2023.

Hellweger, F. L.: 100 Years since Streeter and Phelps: It Is Time To Update the Biology in Our Water Quality Models, Environ. Sci. Technol., 49, 6372–6373, https://doi.org/10.1021/acs.est.5b02130, 2015.

Jähne, B., Heinz, G., and Dietrich, W.: Measurement of the diffusion coefficients of sparingly soluble gases in water, J. Geophys. Res.-Oceans, 92, 10767–10776, https://doi.org/10.1029/JC092iC10p10767, 1987.

Korppoo, M., Huttunen, M., Huttunen, I., Piirainen, V., and Vehviläinen, B.: Simulation of bioavailable phosphorus and nitrogen loading in an agricultural river basin in Finland using VEMALA v.3, J. Hydrol., 549, 363–373, https://doi.org/10.1016/j.jhydrol.2017.03.050, 2017.

Lancelot, C., Veth, C., and Mathot, S.: Modelling ice-edge phyto-
plankton bloom in the Scotia-Weddell sea sector of the southern
ocean during spring 1998, J. Mar. Syst., 2, 333–346, 1991.

Lauerwald, R., Hartmann, J., Moosdorf, N., Dürr, H. H., and Kempe, S.: Retention of dissolved silica within the fluvial system of the conterminous USA, Biogeochemistry, 112, 637–659, 2013.

Le, T. P. Q., Billen, G., Garnier, J., Théry, S., Ruelland, D., Anh, N. X., and Van, M. C.: Nutrient (N, P, Si) transfers in the subtropical Red River system (China and Vietnam): Modelling and budget of nutrient sources and sinks, J. Asian Earth Sci., 37, 259–274, https://doi.org/10.1016/j.jseas.2009.08.010, 2010.

Le, T. P. Q., Billen, G., Garnier, J., and Chau, V.: Long-term biogeochemical functioning of the Red River (Vietnam): past and present situations, Reg. Environ. Change, 15, 329–339, https://doi.org/10.1007/s10113-014-0646-4, 2015.

Lindenschmidt, K.-E., Carr, M. K., Sadeghian, A., and Morales-
Marin, L.: CE-QUAL-W2 model of dam outflow elevation im-
pact on temperature, dissolved oxygen and nutrients in a reser-
voir, Sci. Data, 6, 312, https://doi.org/10.1038/s41597-019-0316-
y, 2019.

Luu, M. T., Dinh, T. D., Trinh, D. A., and Doc, N. T.: Water Quality in an Urbanized River Basin Impacted by Multi-Pollution Sources: From Comprehensive Surveys to Modelling, Sci-

enceAsia, 47, 86-95, https://doi.org/10.2306/scienceasia1513-1874.2021.014, 2021.

Maavara, T., Chen, Q., Van Meter, K., Brown, L. E., Zhang, J., Ni, J., and Zarfl, C.: River dam impacts on biogeochemical cycling, Nat. Rev. Earth Environ., 1, 103–116, https://doi.org/10.1038/s43017-019-0019-0, 2020.

Manteaux, S., Sauvage, S., Samie, R., Monteil, C., Garnier, J., Thieu, V., Cakir, R., and Sánchez-Pérez, J.-M.: Modeling in-stream biogeochemical processes at catchment scale: Coupling SWAT and RIVE models, Environ. Modell. Softw., 170, 105856, https://doi.org/10.1016/j.envsoft.2023.105856, 2023.

Marescaux, A., Thieu, V., Gypens, N., Silvestre, M., and Gar-
nier, J.: Modeling inorganic carbon dynamics in the Seine River
continuum in France, Hydrol. Earth Syst. Sci., 24, 2379–2398,
https://doi.org/10.5194/hess-24-2379-2020, 2020.

Michaelis, L. and Menten, M. L.: Die kinetik der invertinwirkung - The kinetics of invertin action, Biochem. Z., 49, 333–369, 1913.

Minaudo, C., Curie, F., Jullian, Y., Gassama, N., and Moatar, F.: QUAL-NET, a high temporal-resolution eutrophication model for large hydrographic networks, Biogeosciences, 15, 2251–2269, https://doi.org/10.5194/bg-15-2251-2018, 2018.

Monod, J.: The growth of bacterial cultures, Annu. Rev. Microbiol., 3, 371–394, https://doi.org/10.1146/annurev.mi.03.100149.002103, 1949.

Nguyen, H. T. M., Billen, G., Garnier, J., Rochelle-Newall, E., Ri-
bolzi, O., Servais, P., and Le, T. P. Q.: Modelling of faecal in-
dicator bacteria (FIB) in the Red River basin (Vietnam), Env-
iron. Monit. Assess., 188, 517, https://doi.org/10.1007/s10661-
016-5528-4, 2016.

Odum, H. T.: Primary production in Flowing waters, Limnol. Oceanogr., 1, 795–801, 1956.

Platt, T., Gallegos, C., and Harrison, W.: Photoinhibition of pho-
tosynthesis in natural assemblages of marine phytoplankton, J.
Mar. Res., 38, 687–701, 1980.

Raimonet, M., Vilmin, L., Flipo, N., Rocher, V., and Laverman, A.: Modelling the fate of nitrite in an urbanized river using experimentally obtained nitrifier growth parameters, Water Res., 73, 373–387, https://doi.org/10.1016/j.watres.2015.01.026, 2015.

Raimonet, M., Thieu, V., Silvestre, M., Oudin, L., Rabouille, C., Vautard, R., and Garnier, J.: Landward perspective of coastal eutrophication potential under future climate change: The Seine River case (France), Front. Mar. Sci., 5, 136, https://doi.org/10.3389/fmars.2018.00136, 2018.

Redfield, A., Ketchum, B., and Richards, F.: The Sea. Ideas and Observations on Progress in the Study of the Seas. The Composition of the Sea-Water Comparative and Descriptive Oceanography, vol. 2, in: The influence of organisms on the composition of seawater, Interscience Publishers, 26–77, ISBN 9780674017283, 1963.

Rego, J. V., Billen, G., Fontigny, A., and Somville, M.: Free and Attached proteolytic activity in water environments, Mar. Ecol. Prog. Ser., 21, 245–249, 1985.

Rickert, D., Schlüter, M., and Wallmann, K.: Dissolution kinetics of biogenic silica from the water column to the sediments, Geochim. Cosmochim. Ac., 66, 439–455, https://doi.org/10.1016/S0016-7037(01)00757-8, 2002.

Rodríguez-Castillo, T., Estévez, E., González-Ferreras, A., and Barquín, J.: Estimating Ecosystem Metabolism
to Entire River Networks, Ecosystems, 22, 892–911, https://doi.org/10.1007/s10021-018-0311-8, 2019.

Romero, E., Garnier, J., Billen, G., Ramarson, A., Riou, P., and Le Gendre, R.: Modeling the biogeochemical functioning of the Seine estuary and its coastal zone: Export, retention, and transformations, Limnol. Oceanogr., 64, 895–912, https://doi.org/10.1002/lno.11082, 2019.

Ruelland, D., Billen, G., Brunstein, D., and Garnier, J.: SENEQUE: a multi-scaling GIS interface to the Riverstrahler model of the biogeochemical functioning of river systems, Sci. Total Environ., 375, 257–73, https://doi.org/10.1016/j.scitotenv.2006.12.014, 2007.

Segatto, P., Battin, T., and Bertuzzo, E.: Modeling the coupled dynamics of stream metabolism and microbial biomass, Limnol. Oceanogr., 65, 1573–1593, https://doi.org/10.1002/lno.11407, 2020.

Servais, P. and Garnier, J.: Contribution of Heterotrophic Bacterial Production to the Carbon Budget of the River Seine (France), Microb. Ecol., 25, 19–33, 1993.

Servais, P., Billen, G., and Vives-Rego, J.: Rate of Bacterial Mortality in Aquatic Environments, Appl. Environ. Microbiol., 49, 1448–1454, 1985.

Servais, P., Billen, G., and Hascoët, M.-C.: Determination of the biodegradable fraction of dissolved organic matter in waters, Water Res., 21, 445–450, 1987.

Servais, P., Billen, G., Martinez, J., and Vives-Rego, J.: Estimating bacterial mortality by the disappearance of 3H-labeled intracellular DNA, FEMS Microbiol. Lett., 62, 119–125, https://doi.org/10.1016/0378-1097(89)90021-9, 1989.

Servais, P., Barillier, A., and Garnier, J.: Determination of the biodegradable fraction of dissolved and particulate organic carbon in waters, Ann. Limnol.-Int. J. Lim., 31, 75–80, https://doi.org/10.1051/limn/1995005, 1995.

Servais, P., Billen, G., Goncalves, A., and Garcia-Armisen, T.: Modelling microbiological water quality in the Seine river drainage network: past, present and future situations, Hydrol. Earth Syst. Sci., 11, 1581–1592, https://doi.org/10.5194/hess-11-1581-2007, 2007.

Sferratore, A., Billen, G., Garnier, J., Smedberg, E., Humborg, C., and Rahm, L.: Modelling nutrient fluxes from sub-arctic basins: Comparison of pristine vs. dammed rivers, J. Marine Syst., 73, 236–249, https://doi.org/10.1016/j.jmarsys.2007.10.012, 2008.

Somville, M. and Billen, G.: A method for determining exoproteolytic activity in natural water, Limnol. Oceanogr., 28, 190–193, 1983.

Streeter, H. and Phelps, E. B.: A study of the pollution and natural purification of the Ohio River, Tech. Rep. 146, U.S. Public Health Service, Treasury Department, Washington DC, 1925.

Stumm, W. and Morgan, J.: Aquatic Chemistry: Chemical Equilibrium and Rates in Natural Waters, A Wiley-interscience publication, Wiley, ISBN 9780471511854, 1996.

Thieu, V., Guillon, T., Billen, G., Garnier, J., and Thouvenot, M.: Applicatif BARMAN – Notice d'utilisation., Rapport annuel PIREN-Seine, PIREN-Seine, https://doi.org/10.26047/PIREN.rapp.ann.2006.vol07, 2006.

Thieu, V., Billen, G., and Garnier, J.: Nutrient transfer in three contrasting NW European watersheds: The Seine, Somme, and Scheldt Rivers. A comparative application of

the Seneque/Riverstrahler model, Water Res., 43, 1740–1754, https://doi.org/10.1016/j.watres.2009.01.014, 2009.

Thieu, V., Mayorga, E., Billen, G., and Garnier, J.: Subregional and downscaled global scenarios of nutrient transfer in river basins: Seine-Somme-Scheldt case study, Global Biogeochem. Cy., 24, GB0A10, https://doi.org/10.1029/2009GB003561, 2010.

Thieu, V., Silvestre, M., Garnier, J., and Billen, G.: Introducing the biogeochemical pynuts-riverstrahler model to assess prospective scenario impact along the aquatic continuum in western eu-rivers, in: Proceedings of ASLO Aquatic Sciences Meeting: "Moutains to the Sea", 27 February 2017, Honolulu, USA, 28904, 2017.

Thieu, V., Silvestre, M., Wang, S., Marescaux, A., Yan, X., Garnier, J., and Billen, G.: pyRIVE, InDoRES [code], https://doi.org/10.48579/PRO/Z9ACP1, 2023.

Thouvenot, M., Billen, G., and Garnier, J.: Modelling nutrient exchange at the sediment-water interface of river systems, J. Hydrol., 341, 55–78, 2007.

Vannote, R., Minshall, G. W., Cummins, K., Sedell, J., and Cushing, C.: The River Continuum Concept, Can. J. Fish. Aquat. Sci., 37, 130–137, https://doi.org/10.1139/f80-017, 1980.

Vilmin, L., Aissa-Grouz, N., Garnier, J., Billen, G., Mouchel, J. M., Poulin, M., and Flipo, N.: Impact of hydro-sedimentary processes on the dynamics of soluble reactive phosphorus in the Seine River, Biogeochemistry, 122, 229–251, https://doi.org/10.1007/s10533-014-0038-3, 2015a.

Vilmin, L., Flipo, N., de Fouquet, C., and Poulin, M.: Pluriannual sediment budget in a navigated river system: The Seine River (France), Sci. Total Environ., 502, 48–59, https://doi.org/10.1016/j.scitotenv.2014.08.110, 2015b.

Vilmin, L., Flipo, N., Escoffier, N., Rocher, V., and Groleau, A.: Carbon fate in a large temperate human-impacted river system: Focus on benthic dynamics, Global Biogeochem. Cy., 30, 1086–1104, https://doi.org/10.1002/2015GB005271, 2016.

Vilmin, L., Flipo, N., Escoffier, N., and Groleau, A.: Estimation of the water quality of a large urbanized river as defined by the European WFD: what is the optimal sampling frequency?, Environ. Sci. Pollut. Res., 25, 23485–23501, https://doi.org/10.1007/s11356-016-7109-z, 2018.

Wang, S.: Simulation du métabolisme de la Seine par assimilation de données en continu, Theses, Mines Paris – PSL, https://pastel.archives-ouvertes.fr/tel-02388690 (last access: 20 December 2023), 2019.

Wang, S., Flipo, N., and Romary, T.: Time-dependent global sensitivity analysis of the C-RIVE biogeochemical model in contrasted hydrological and trophic contexts, Water Res., 144, 341–355, https://doi.org/10.1016/j.watres.2018.07.033, 2018.

Wang, S., Flipo, N., and Romary, T.: Oxygen data assimilation for estimating micro-organism communities' parameters in river systems, Water Res., 165, 115021, https://doi.org/10.1016/j.watres.2019.115021, 2019.

Wang, S., Flipo, N., Romary, T., and Hasanyar, M.: Particle filter for high frequency oxygen data assimilation in river systems, Environ. Modell. Softw., 151, 105382, https://doi.org/10.1016/j.envsoft.2022.105382, 2022.

Wang, S., Flipo, N., and Romary, T.: Which filter for data assimilation in water quality models? Focus on oxygen reaeration and heterotrophic bacteria activity, J. Hydrol., 620, 129423, https://doi.org/10.1016/j.jhydrol.2023.129423, 2023a.
Wang, S., Vilmin, L., Hasanyar, M., and Flipo, N.: C-RIVE, Zenodo [code], https://doi.org/10.5281/zenodo.7849609, 2023b.

Wang, S., Thieu, V., Billen, G., Garnier, J., Silvestre, M., Marescaux, A., Yan, X., and Flipo, N.: Dataset for the paper ”The community-centered freshwater biogeochemistry model unified RIVE v1.0: a unified version for water column”, Zenodo, [data set], https://doi.org/10.5281/zenodo.10490669, 2024.

Wanninkhof, R.: Relationship between wind speed and gas ex-
change over the ocean, J. Geophys. Res.-Oceans, 97, 7373–7382,
https://doi.org/10.1029/92JC00188, 1992.

Weiss, R.: Carbon dioxide in water and seawater: the sol-
ubility of a non-ideal gas, Mar. Chem., 2, 203–215,
https://doi.org/10.1016/0304-4203(74)90015-2, 1974.

Wilke, C. R. and Chang, P.: Correlation of diffusion co-
efficients in dilute solutions, AIChE J., 1, 264–270,
https://doi.org/10.1002/aic.690010222, 1955.

Yan, X., Garnier, J., Billen, G., Wang, S., and Thieu, V.: Unrav-
elling nutrient fate and CO2 concentrations in the reservoirs of
the Seine Basin using a modelling approach, Water Res., 225,
119135, https://doi.org/10.1016/j.watres.2022.119135, 2022a.

Yan, X., Thieu, V., Wu, S., and Garnier, J.: Reservoirs change pCO2 and water quality of downstream rivers: Evidence from three reservoirs in the Seine Basin, Water Res., 213, 118158, https://doi.org/10.1016/j.watres.2022.118158, 2022b.