the Creative Commons Attribution 4.0 License.
the Creative Commons Attribution 4.0 License.
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
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.
- Article
(4934 KB) - Full-text XML
- BibTeX
- EndNote
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 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 AQUAPHY 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 hydro-biodynamic 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 programming 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.
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 (HD1 and HP1), (2) slowly biodegradable in 45 d (HD2 and HP2), and (3) refractory (HD3 and HP3).
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 (bmax20,hb) and a bacterial growth yield (Yhb) are used to calculate the growth rate of heterotrophic bacteria () (Eq. 3). The fraction of uptake not used for growth (1−Yhb) is respired.
Here, is the effective substrate uptake rate of the ith species of heterotrophic bacteria [h−1], is the maximal substrate uptake rate of the ith species of heterotrophic bacteria at 20 ∘C [h−1], [SMS] is the small monomeric substrate concentration [mgC L−1], is the half-saturation constant for small monomeric substrate of the ith species of heterotrophic bacteria [mgC L−1], is the water temperature weight of the ith species of heterotrophic bacteria at T ∘C [–], is the optimal temperature of the ith species of heterotrophic bacteria for its growth [∘C], is the range of temperature for the ith species of heterotrophic bacteria [∘C], is the bacterial growth yield of the ith species of heterotrophic bacteria [–], and is the effective growth rate of the ith species of heterotrophic bacteria [h−1]
A sinking velocity (vshb) 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 (ϵhd and ϵhp, Fig. 1).
Here, is the effective growth rate of the ith species of heterotrophic bacteria [h−1], is the mortality rate of the ith species of heterotrophic bacteria at 20 ∘C [h−1], is the sinking rate of the ith species of heterotrophic bacteria [h−1], [HBi] is the biomass concentration of the ith species of heterotrophic bacteria [mgC L−1], is the water temperature weight at T ∘C defined by Eq. (2) [–], is the sinking velocity of the ith species of heterotrophic bacteria [m h−1], and “depth” represents the water depth [m].
2.1.2 Hydrolysis of high-weight polymer
The particulate biodegradable high-weight polymer (HP1 and HP2) is firstly hydrolyzed to the dissolved biodegradable high-weight polymer (HD1 and HD2). The dissolved biodegradable high-weight polymer is then hydrolyzed exoenzymatically 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).
Here, [HPi] is the concentration of particulate high-weight polymer [mgC L−1], is the hydrolysis rate of HPi () [h−1], f(T)j is the water temperature weight of the jth living species at T ∘C defined as in Eq. 2 [–], kd20,j is the mortality rate of the jth living species (such as phytoplankton, zooplankton, bacteria) at 20 ∘C [h−1], [LS]j is the concentration of the jth living species [mgC L−1], is the proportion to convert the dead biomass to HPi () [–], and is the sinking rate for HPi () [h−1].
Here, [HDi] is the concentration of dissolved high-weight polymer () [mgC L−1], is the maximum hydrolysis rate of HDi at 20 ∘C related to HBk and [h−1], is the water temperature weight of the kth species of heterotrophic bacteria at T ∘C (Eq. 2) [–], is the half-saturation constant for HDi related to HBk and [mgC L−1], [HBk] is the concentration of the kth species of heterotrophic bacteria [mgC L−1], kd20,j is the mortality rate of the jth living species (such as phytoplankton, zooplankton, bacteria) at 20 ∘C [h−1], f(T)j is the water temperature weight of the jth living species at T ∘C defined as in Eq. (2) [–], [LS]j is the concentration of the jth living species [mgC L−1], and is the proportion to convert the dead biomass to HDi and [–].
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 () and nitrate (). The nitrifying bacteria get energy by oxidizing NH (ammonium) and NO (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.
Here, μaob and μnob are the effective growth rates of AOB and NOB [h−1], μmax20,aob and μmax20,nob are the maximal growth rates of AOB and NOB at 20 ∘C, respectively [h−1], f(T)aob and f(T)nob are the water temperature weight at T ∘C defined as in Eq. (2) [–], and are the half-saturation constants for NH (AOB) and for NO (NOB) [mgN L−1], and and are the half-saturation constants for oxygen (AOB and NOB) [mgO2 L−1].
The mortality and sinking of nitrifying bacteria are simulated the same way as for other living species.
Here, μaob and μnob are the effective growth rates of AOB and NOB defined by Eqs. (7) and (8) [h−1], kd20,aob and kd20,nob are the mortality rates of AOB and NOB at 20 ∘C [h−1], ksink,aob and ksink,nob are the sinking rates of AOB and NOB [h−1], f(T)aob and f(T)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.
At any time, the biomass of the jth phytoplankton species (mgC L−1), [PHY]j, is equal to the sum of the three internal constituents (Eq. 11), [F]j, [R]j, and [S]j :
The most common way of measuring phytoplankton biomass is using the chlorophyll a concentration (µg chl a L−1). A carbon to chlorophyll a ratio of 35 mgC µg chl a−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−2 s−1).
Here, is the photosynthesis rate of the jth phytoplankton species at water depth z m [h−1], is the maximal photosynthesis rate of the jth phytoplankton species at 20 ∘C [h−1], is the water temperature weight of the jth phytoplankton species at T ∘C defined as in Eq. (2) [–], is the photosynthetic efficiency of the jth phytoplankton species [h−1 (µE m−2 s−1)−1], is the photoinhibition capacity of the jth phytoplankton species [h−1 (µE m−2 s−1)−1], and I(z) is photosynthetically active radiation (PAR) or active irradiance in the water column at depth z m [µE m−2 s−1] or [W m−2].
The averaged photosynthesis rate of the jth phytoplankton species over water column is obtained by integrating P(z):
where depth is the water height [m] and is the averaged photosynthesis rate over the water column [h−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 (η). The extinction coefficient is composed of three parts: pure water (ηbase), suspended solid (ηss), and algal self-shading (ηchl a).
Here, I(z) is the active irradiance at water depth z [µE m−2 s−1 or W m−2], I0 is photosynthetically active radiation (PAR) or active irradiance at the water surface measured by the photosynthetic photon1 flux density (PPFD) [µE m−2 s−1 or W m−2], η is the light extinction coefficient [m−1], ηbase is the light extinction coefficient related to pure water [m−1], ηchl a is the linear algal self-shading light extinction coefficient [m−1 (µg chl a L−1)−1], ηss is the light extinction coefficient related to suspended solid [m−1 (mg L−1)−1], [chl a] is the total chlorophyll a concentration [µg chl a L−1], and [SS] is the suspended solid concentration [mg L−1].
2.3.2 Growth
The growth of phytoplankton involves the transformation of small precursors (S) into structural and functional macromolecules (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 () 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).
Here, is the effective growth rate of functional macromolecules for the jth phytoplankton species [h−1], is the maximal growth rate of functional macromolecules for the jth phytoplankton species at 20 ∘C [h−1], is the water temperature weight of the jth phytoplankton species at T ∘C defined as in Eq. (2) [–], [Sj] and [Fj] are the concentrations of small precursors and functional macromolecules for the jth phytoplankton species [mgC L−1], 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 (, mgN L−1), dissolved inorganic phosphorus (, mgP L−1), and dissolved silica (DSi, mgSi L−1), and are the half-saturation constants for dissolved inorganic nitrogen and phosphorus of the jth phytoplankton species [mgN L−1 and mgP L−1], and is the half-saturation constant for dissolved silica in the case of diatoms [mgSi L−1].
2.3.3 Respiration
The respiration rate of phytoplankton (rphy) is divided into two components (Eq. 17): one (Rm,phy) ensuring the survival of the cell (maintenance process) and the other (Rμ,phy) corresponding to the energetic cost of growth.
Here, is the respiration rate of the jth phytoplankton species [h−1], is the maintenance respiration rate of the jth phytoplankton species at 20 ∘C [h−1], is the water temperature weight of the jth phytoplankton species at T ∘C defined as in Eq. (2) [–], is the respiration for the energetic cost of the jth phytoplankton species [–], and is the effective growth rate of the jth phytoplankton species (Eq. 15) [h−1].
2.3.4 Excretion
Included later by Garnier et al. (1998), the phytoplankton excretion (ephy) includes two terms: a constant excretion rate (Ecst,phy) and another that depends on the photosynthesis rate (Ephot,phy). The product of excretion is the small monomeric substrate (SMS), assimilated directly by heterotrophic bacteria for their growth and respiration (Fig. 1).
Here, is the excretion rate of the jth phytoplankton species [h−1], is the basic excretion rate of the jth phytoplankton species [h−1], is the excretion of the jth phytoplankton species related to photosynthesis [–], and is the photosynthesis rate of the jth phytoplankton species (Eq. 13) [h−1].
The variation of small monomeric substrate (SMS) can then be established (Eq. 19).
Here, hydr is the hydrolysis of the dissolved high-weight polymers HD1 and HD2 (Eq. 6) [mgC L−1 h−1], is the effective rate of substrate uptake by the ith heterotrophic bacteria species (Eq. 3) [h−1], [HBi] is the biomass concentration of the ith species of heterotrophic bacteria [mgC L−1], is the effective excretion rate of the jth phytoplankton species (Eq. 18) [h−1], and [Fj] is the concentration of functional macromolecules for the jth phytoplankton species [mgC L−1].
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 ratio based on a Michaelis–Menten-like function (Eq. 20).
Here, is the synthesis rate of reserve products of jth phytoplankton species [h−1], is the maximal synthesis rate of reserve products of jth phytoplankton species at 20 ∘C [h−1], is the water temperature weight of jth phytoplankton species at T ∘C defined as in Eq. (2) [–], [Sj] and [Fj] are the concentrations of small precursors and functional macromolecules for the jth phytoplankton species [mgC L−1], and is the half-saturation constant for small precursors of the jth phytoplankton species [–].
Reserve products (R) are likely to be catabolized to produce small precursors (S). A first-order kinetic (cR,phy, h−1) 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 zooplankton (Sect. 2.4.1). The phytoplankton lysis is represented by first-order kinetics using a mortality rate (kd20,phy, h−1). For ease of presentation, all three processes are assumed in an overall extinction rate dphy (h−1).
Here, is the extinction rate of the jth phytoplankton species [h−1], is the mortality rate of the jth phytoplankton species at 20 ∘C [h−1], is the water temperature weight of jth phytoplankton species at T ∘C defined as in Eq. (2) [–], is the sinking rate of the jth phytoplankton species [h−1], is the grazing rate of the ith zooplankton species (Eq. 26, Sect. 2.4.1) [h−1], [ZOOi] is the zooplankton concentration of the ith zooplankton species [mgC L−1], and is the total phytoplankton concentration (with NS the number of phytoplankton species grazed by zooplankton) [mgC L−1].
2.3.7 Phytoplankton budgets
According to the processes related to phytoplankton (photosynthesis, growth, mortality, etc.), different budgets can be established for the jth phytoplankton species as follows.
Here, is the photosynthesis rate of the jth phytoplankton species (Eq. 13) [h−1], is the respiration rate of the jth phytoplankton species (Eq. 17) [h−1], is the growth rate of the jth phytoplankton species (Eq. 15) [h−1], is the synthesis rate of reserve products of the jth phytoplankton species (Eq. 20) [h−1], is the catabolysis rate of reserve products of the jth phytoplankton species [h−1], is the excretion rate of the jth phytoplankton species (Eq. 18) [h−1], is the extinction rate of the jth phytoplankton species (Eq. 21) [h−1], [Sj], [Fj], and [Rj] are the concentrations of [Sj], [Fj], and [Rj] for the jth phytoplankton species [mgC L−1], and [PHYj] is the biomass concentration of the jth phytoplankton species [mgC L−1].
2.4 Zooplankton dynamics
The zooplankton dynamics include grazing on phytoplankton 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 (bmax20,zoo) limited by the phytoplankton biomass based on a Monod function (Eq. 26). The grazing of zooplankton takes place only when the total phytoplankton biomass exceeds a certain threshold ([PHY0]). No specific preference for grazing on particular phytoplankton species is considered among zooplankton species. Instead, the phytoplankton biomass grazed by the ith species of zooplankton is divided proportionally among each species of phytoplankton (Eq. 21). The growth rate of zooplankton is considered proportional to the grazing rate using a growth yield factor (Eq. 27).
Here, is the effective grazing rate of the ith zooplankton species [h−1], is the maximal grazing rate of the ith zooplankton species at 20 ∘C [h−1], is the water temperature weight of the ith zooplankton species at T ∘C defined as in Eq. (2) [–], is the total phytoplankton biomass with NS the number of phytoplankton species grazed by zooplankton [mgC L−1], is the phytoplankton biomass threshold above which grazing takes place for the ith zooplankton species [mgC L−1], is the half-saturation constant for phytoplankton biomass of the ith zooplankton species [mgC L−1], is the growth rate of the ith zooplankton species [h−1], and 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 . The mortality of zooplankton is simulated by first-order kinetics (kd20,zoo).
Here, is the respiration rate of the ith zooplankton species [h−1], is the growth yield of the ith zooplankton species [–], is the effective grazing rate of the ith zooplankton species (Eq. 26) [h−1], is the mortality rate of the ith zooplankton species at 20 ∘C [h−1], is the water temperature weight of the ith zooplankton species at T ∘C defined as in Eq. (2) [–], is the sinking rate of the ith zooplankton species [h−1], and [ZOOi] is the biomass concentration of the ith zooplankton species [mgC L−1].
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 ( ; – Redfield et al., 1963; Conley et al., 1989) is used to determine the composition of carbon, nitrogen, and phosphorus in organic matter. Constant , , and mass ratios are considered to calculate the uptake of nutrients associated with phytoplankton growth.
Here, is the effective growth rate of the ith phytoplankton species (Eq. 15) [h−1], is the excretion rate of the ith phytoplankton species (Eq. 18) [h−1], [Fi] represents the functional macromolecule concentration of the ith phytoplankton species [mgC L−1], and are the concentrations of ammonium and nitrate [mgN L−1], uptN is the uptake of nitrogen for phytoplankton growth [mgN L−1], is the uptake of NH for phytoplankton growth [mgN L−1], is the uptake of NO for phytoplankton growth [mgN L−1], uptP is the uptake of phosphorus for phytoplankton growth [mgP L−1], is the carbon to nitrogen mass ratio [mgC mgN−1], is the carbon to phosphorus mass ratio [mgC mgP−1], is the carbon to silica mass ratio [mgC mgSi−1], uptSi is the uptake of silica for diatom growth [mgSi L−1], μF,dia is the effective growth rate of diatoms [h−1], and [Fdia] is the functional macromolecule (F) concentration of diatoms [mgC L−1]
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.
Here, respHB is the respiration of heterotrophic bacteria species [mgC L−1 h−1], is the growth yield of the ith heterotrophic bacteria species [–], is the effective rate of substrate uptake by the ith heterotrophic bacteria species (Eq. 1) [h−1], [HBi] is the concentration of the ith heterotrophic bacteria species [mgC L−1], respZOO is the respiration of zooplankton species [mgC L−1 h−1], is the growth yield of the jth zooplankton species [–], is the effective grazing rate of the jth zooplankton species (Eq. 26) [h−1], [ZOOj] is the concentration of the jth zooplankton species [mgC L−1], relN is the release of nitrogen [mgN L−1 h−1], is the carbon to nitrogen mass ratio [mgC mgN−1], relP is the release of phosphorus [mgP L−1 h−1], and is the carbon to phosphorus mass ratio [mgC mgP−1].
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 ( and ) 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).
Here, and are the effective growth rates of the ith AOB species (Eq. 7) and the jth NOB species (Eq. 8) [h−1], and are the growth yields of the ith AOB species and the jth NOB species [mgC mgN−1], nitraob is nitrification calculated as [mgN L−1 h−1], nitrnob is nitrification calculated as [mgN L−1 h−1], and [AOBi] and [NOBj] are the biomass concentrations of the ith AOB species and the jth NOB species [mgC L−1].
The budgets of , , and can then be established.
Here, denit is denitrification [mgN L−1 h−1], nitrnob is nitrification by NOB (Eq. 40) [mgN L−1 h−1], nitraob is nitrification by AOB (Eq. 39) [mgN L−1 h−1], is the uptake of NO by phytoplankton growth (Eq. 32) [mgN L−1 h−1], relN is the release of nitrogen by respiration of heterotrophic bacteria and zooplankton (Eq. 37) [mgN L−1 h−1], and is the uptake of NH by phytoplankton growth (Eq. 31) [mgN L−1 h−1].
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
with the inorganic P content of MSS [mgP mgMSS−1], [PIP] and [MSS] the concentrations of PIP and MSS [mgP L−1 and mgMSS L−1], Pac the maximum adsorption capacity of MSS [mgP mgMSS−1], the concentration of orthophosphate [mgP L−1], and Kps the half-saturation adsorption constant [mgP L−1].
Considering this equilibrium to be instantaneously reached implies that a relationship exists between the variables PIP, MSS, PO, and TIP (total inorganic phosphorus):
This equilibrium relationship can be written as
with the concentration of orthophosphate [mgP L−1], [TIP] and [MSS] the concentrations of TIP and MSS [mgP L−1 and mgMSS L−1], Pac the maximal adsorption capacity of MSS [mgP mgMSS−1], and Kps the half-saturation adsorption constant [mgP L−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):
with KbSi20 the dissolution rate of biogenic silica at 20 ∘C [h−1], [BSi] the concentration of biogenic silica [mgSi L−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).
Here, krea is the re-aeration coefficient [m h−1], “depth” indicates water height [m], [O2]sat is the saturated concentration of dissolved oxygen in water [mgO2 L−1], [O2] is the concentration of dissolved oxygen in water [mgO2 L−1], is the molar mass ratio between dissolved oxygen and carbon [mgO2 mgC−1], is the photosynthesis rate of the ith phytoplankton species (Eq. 13) [h−1], is the respiration rate of the ith phytoplankton species (Eq. 17) [h−1], [Fi] is the functional biomass concentration of the ith phytoplankton species [mgC L−1], respHB is the respiration of all heterotrophic bacteria species (Eq. 35) [mgC L−1 h−1], is the respiration rate of the jth zooplankton species (Eq. 28) [h−1], [ZOOj] is the biomass concentration of the jth zooplankton species [mgC L−1], is the molar mass ratio between dissolved oxygen and nitrogen [mgO2 mgN−1], nitraob is nitrification to produce nitrite by oxidizing (Eq. 39, with the stoichiometric coefficient) [mgN L−1 h−1], and nitrnob is nitrification to produce nitrate by oxidizing (Eq. 40, with the stoichiometric coefficient) [mgN L−1 h−1].
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 CO2 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 – CO2((aq)).
2.7.1 CO2 flux at the air–water interface
The DIC is defined as the sum of three dissolved carbonate species:
The calculation of pH is derived from Culberson (1980) using TA and DIC. Then the aqueous carbon dioxide (CO2 (aq)) is derived from the carbonate chemical equilibrium using DIC and pH (Marescaux et al., 2020; Yan et al., 2022a).
Here, K1 and K2 are the equilibrium constants of carbonate equilibrium reactions (Stumm and Morgan, 1996) [mol L−1], [H+] is the concentration of hydrogen ions with pH [mol L−1], and [DIC] is the concentration of dissolved inorganic carbon [mgC L−1].
The flux of CO2 at the water–air interface (, gC m−2 h−1) is calculated based on Fick's first law (Fick, 1855) with a gas transfer velocity of CO2 ().
Here, kco2 is the gas transfer velocity of CO2 [m h−1], [CO2(sat)] is the solubility of CO2 in water calculated based on Henry's law (Weiss, 1974) [mgC L−1], and [CO2(aq)] is the aqueous carbon dioxide concentration [mgC L−1].
The gas transfer velocity of CO2 (kco2) depends on water temperature and k600 (gas transfer velocity of CO2 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 CO2 (kco2) at water temperature T (∘C) can be calculated as
where k600 (m h−1) is the gas transfer velocity of CO2 for a Schmidt number of 600, and is the Schmidt number (dimensionless) calculated with the water temperature in degrees Celsius (∘C). The can be determined as,=
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.
Here, TANet_input (µmol L−1 h−1) and DICNet_input (mgC L−1 h−1) are the net input fluxes. The respiration of all phytoplankton, bacteria, and zooplankton species (respPHY, respHB, respZOO; mgC L−1 h−1) transforms organic carbon to CO2 by full oxidization. The denitrification (denit, mgN L−1 h−1) is also considered in the calculation of TA and DIC. (gC m−2 h−1) is the CO2 flux at the air–water interface, and “depth” is the water depth (m). , and 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.
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 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 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).
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 (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 m3 s−1, 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 software packages are (pyRIVE 1.0 and C-RIVE 0.32).
Here, Ci represents the concentrations simulated by C-RIVE 0.32 (in ANSI C) and Pyi 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 mgO2 L−1, which is relatively low, is obtained for dissolved oxygen concentration. The mean absolute difference (MAD) for dissolved oxygen concentration is 0.00678 mgO2 L−1 (Table 3) and the corresponding percent bias (PBIAS) is 0 %. The MaAD of 0.0307 mgO2 L−1 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−1 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 ). 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 %.
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−1, 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−1. 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 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 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., 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.
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.
The 120 parameter values necessary for running unified RIVE v1.0 are provided hereafter.
* Parameters depend on water temperature and are multiplied by , where T is water temperature in ∘C.
* Parameters depend on water temperature and are multiplied by , where T is water temperature in ∘C.
* Parameters depend on water temperature and are multiplied by , where T is water temperature in ∘C.
* Parameters depend on water temperature and are multiplied by , where T is water temperature in ∘C.
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).
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.
The contact author has declared that none of the authors has any competing interests.
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.
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.
This paper was edited by Vassilios Vervatis and reviewed by two anonymous referees.
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 modeling algal growth, Environ. Sci. Pollut. Res., 25, 23515–23528, https://doi.org/10.1007/s11356-016-7555-7, 2016. a
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. a
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. a
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. a
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. a, b
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. a, b
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. a
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. a
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. a, b, c, d, e, f, g, h
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. a, b
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. a, b, c
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. a, b
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é H14CO et son application pour estimer des biomasses de bactéries nitrifiantes, Revue des Sciences de l'Eau, 11, 283–302, 1998. a
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. a
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. a
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. a, b
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. a
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. a, b, c
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. a, b, c, d, e
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. a, b, c
Fick, A.: Ueber diffusion, Annalen der Physik und Chemie, J.A. Barth, https://doi.org/10.1002/andp.18551700105, 1855. a
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. a, b, c, d, e
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. a, b
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. a
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. a
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. a, b, c, d
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. a
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. a, b, c
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. a
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. a, b, c, d, e
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. a
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. a, b
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. a
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. a, b, c, d, e
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. a, b, c
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. a, b
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. a, b
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. a, b
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. a, b
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 Approaches to High-Frequency Data Analysis, Handbook of Environmental Chemistry, Springer, Berlin, Heidelberg, 1–28, https://doi.org/10.1007/698_2019_379, 2020. a, b
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. a, b, c
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. a
Gypens, N., Lancelot, C., and Borges, A. V.: Carbon dynamics and CO2 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. a
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. a
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. a, b
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. a
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. a, b
Lancelot, C., Veth, C., and Mathot, S.: Modelling ice-edge phytoplankton bloom in the Scotia-Weddell sea sector of the southern ocean during spring 1998, J. Mar. Syst., 2, 333–346, 1991. a, b, c, d
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. a
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.jseaes.2009.08.010, 2010. a, b
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. a, b
Lindenschmidt, K.-E., Carr, M. K., Sadeghian, A., and Morales-Marin, L.: CE-QUAL-W2 model of dam outflow elevation impact on temperature, dissolved oxygen and nutrients in a reservoir, Sci. Data, 6, 312, https://doi.org/10.1038/s41597-019-0316-y, 2019. a
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, ScienceAsia, 47, 86–95, https://doi.org/10.2306/scienceasia1513-1874.2021.014, 2021. a, b
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. a
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. a, b
Marescaux, A., Thieu, V., Gypens, N., Silvestre, M., and Garnier, 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. a, b, c, d, e, f
Michaelis, L. and Menten, M. L.: Die kinetik der invertinwirkung - The kinetics of invertin action, Biochem. Z., 49, 333–369, 1913. a
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. a, b
Monod, J.: The growth of bacterial cultures, Annu. Rev. Microbiol., 3, 371–394, https://doi.org/10.1146/annurev.mi.03.100149.002103, 1949. a
Nguyen, H. T. M., Billen, G., Garnier, J., Rochelle-Newall, E., Ribolzi, O., Servais, P., and Le, T. P. Q.: Modelling of faecal indicator bacteria (FIB) in the Red River basin (Vietnam), Environ. Monit. Assess., 188, 517, https://doi.org/10.1007/s10661-016-5528-4, 2016. a, b
Odum, H. T.: Primary production in Flowing waters, Limnol. Oceanogr., 1, 795–801, 1956. a
Platt, T., Gallegos, C., and Harrison, W.: Photoinhibition of photosynthesis in natural assemblages of marine phytoplankton, J. Mar. Res., 38, 687–701, 1980. a
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. a
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. a
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 sea-water, Interscience Publishers, 26–77, ISBN 9780674017283, 1963. a
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. a
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. a
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. a
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. a, b
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. a
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. a
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. a
Servais, P., Billen, G., and Vives-Rego, J.: Rate of Bacterial Mortality in Aquatic Environments, Appl. Environ. Microbiol., 49, 1448–1454, 1985. a
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. a
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. a
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. a
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. a
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. a, b
Somville, M. and Billen, G.: A method for determining exoproteolytic activity in natural water, Limnol. Oceanogr., 28, 190–193, 1983. a
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. a, b
Stumm, W. and Morgan, J.: Aquatic Chemistry: Chemical Equilibria and Rates in Natural Waters, A Wiley-interscience publication, Wiley, ISBN 9780471511854, 1996. a
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. a, b, c
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. a, b, c, d, e, f
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. a, b, c, d
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. a, b, c
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. a, b
Thouvenot, M., Billen, G., and Garnier, J.: Modelling nutrient exchange at the sediment-water interface of river systems, J. Hydrol., 341, 55–78, 2007. a, b
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. a
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. a, b, c, d
Vilmin, L., Flipo, N., de Fouquet, C., and Poulin, M.: Pluri-annual 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. a, b, c, d
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. a, b, c, d, e
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. a, b
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. a
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. a
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. a, b, c
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. a, b, c
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. a, b
Wang, S., Vilmin, L., Hasanyar, M., and Flipo, N.: C-RIVE, Zenodo [code], https://doi.org/10.5281/zenodo.7849609, 2023b. a, b
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. a
Wanninkhof, R.: Relationship between wind speed and gas exchange over the ocean, J. Geophys. Res.-Oceans, 97, 7373–7382, https://doi.org/10.1029/92JC00188, 1992. a
Weiss, R.: Carbon dioxide in water and seawater: the solubility of a non-ideal gas, Mar. Chem., 2, 203–215, https://doi.org/10.1016/0304-4203(74)90015-2, 1974. a
Wilke, C. R. and Chang, P.: Correlation of diffusion coefficients in dilute solutions, AIChE J., 1, 264–270, https://doi.org/10.1002/aic.690010222, 1955. a
Yan, X., Garnier, J., Billen, G., Wang, S., and Thieu, V.: Unravelling 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. a, b, c, d, e
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. a, b, c
Photons within the range of visible light between 400 and 700 nm.
- Abstract
- Introduction
- Model description
- Results
- Discussion
- Conclusions
- Appendix A: Various implementations of the RIVE model and its applications in different freshwater systems
- Appendix B: Parameter values for unified RIVE v1.0
- Code and data availability
- Author contributions
- Competing interests
- Disclaimer
- Acknowledgements
- Review statement
- References
- Abstract
- Introduction
- Model description
- Results
- Discussion
- Conclusions
- Appendix A: Various implementations of the RIVE model and its applications in different freshwater systems
- Appendix B: Parameter values for unified RIVE v1.0
- Code and data availability
- Author contributions
- Competing interests
- Disclaimer
- Acknowledgements
- Review statement
- References