Articles | Volume 18, issue 5
Biogeosciences, 18, 1719–1747, 2021
Biogeosciences, 18, 1719–1747, 2021

Research article 11 Mar 2021

Research article | 11 Mar 2021

Modeling silicate–nitrate–ammonium co-limitation of algal growth and the importance of bacterial remineralization based on an experimental Arctic coastal spring bloom culture study

Modeling silicate–nitrate–ammonium co-limitation of algal growth and the importance of bacterial remineralization based on an experimental Arctic coastal spring bloom culture study
Tobias R. Vonnahme1, Martial Leroy2, Silke Thoms3, Dick van Oevelen4, H. Rodger Harvey5, Svein Kristiansen1, Rolf Gradinger1, Ulrike Dietrich1, and Christoph Völker3 Tobias R. Vonnahme et al.
  • 1Department of Arctic and Marine Biology, UiT – The Arctic University of Norway, Tromsø, Norway
  • 2Université Grenoble Alpes, Grenoble, France
  • 3Department of Biogeosciences, Alfred Wegener Institute Helmholtz Centre for Polar and Marine Research, Bremerhaven, Germany
  • 4Department of Estuarine and Delta Systems, Royal Netherlands Institute for Sea Research (NIOZ), and Utrecht University, Texel, Yerseke, the Netherlands
  • 5Department of Ocean and Earth Sciences, Old Dominion University, Norfolk, USA

Correspondence: Tobias R. Vonnahme ( and Christoph Voelker (


Arctic coastal ecosystems are rapidly changing due to climate warming. This makes modeling their productivity crucially important to better understand future changes. System primary production in these systems is highest during the pronounced spring bloom, typically dominated by diatoms. Eventually the spring blooms terminate due to silicon or nitrogen limitation. Bacteria can play an important role for extending bloom duration and total CO2 fixation through ammonium regeneration. Current ecosystem models often simplify the effects of nutrient co-limitations on algal physiology and cellular ratios and simplify nutrient regeneration. These simplifications may lead to underestimations of primary production. Detailed biochemistry- and cell-based models can represent these dynamics but are difficult to tune in the environment. We performed a cultivation experiment that showed typical spring bloom dynamics, such as extended algal growth via bacterial ammonium remineralization, reduced algal growth and inhibited chlorophyll synthesis under silicate limitation, and gradually reduced nitrogen assimilation and chlorophyll synthesis under nitrogen limitation. We developed a simplified dynamic model to represent these processes. Overall, model complexity in terms of the number of parameters is comparable to the phytoplankton growth and nutrient biogeochemistry formulations in common ecosystem models used in the Arctic while improving the representation of nutrient-co-limitation-related processes. Such model enhancements that now incorporate increased nutrient inputs and higher mineralization rates in a warmer climate will improve future predictions in this vulnerable system.

1 Introduction

Marine phytoplankton are responsible for half of the CO2 fixation on Earth (Field et al., 1998; Westberry et al., 2008). In high-latitude oceans, diatoms are an important group contributing 20 %–40 % of the global CO2 fixation (Nelson et al., 1995; Uitz et al., 2010). Marine primary production can be bottom-up limited by light and/or nutrients like nitrogen (N), phosphorous (P), silicon (Si), and iron (Fe). Their availability is affected by pronounced geographical and seasonal variations (Eilertsen et al., 1989; Loebl et al., 2009; Iversen and Seuthe, 2011; Moore et al., 2013). Arctic coasts are one of the fastest changing systems due to climate change. Thus, modeling their dynamics is difficult but crucial for predictions of primary production with climate change (e.g., Slagstad et al., 2015; Fritz et al., 2017; Lannuzel et al., 2020). In Arctic coastal ecosystems, primary production is typically highest in spring. In spring, previous winter mixing supplied fresh nutrients, and a stratified surface layer with sufficient light is facilitated by increasing temperatures and potentially sea ice melt (Sverdrup, 1953; Eilertsen et al., 1989; Eilertsen and Frantzen, 2007; Iversen and Seuthe, 2011). With increasing temperatures and runoff, stratification in coastal Arctic systems is expected to increase (Tremblay and Gagnon, 2009). This will lead to decreased mixing and nutrient upwelling in autumn and winter and an earlier stratified surface layer in spring, which may lead to an earlier spring bloom (Tremblay and Gagnon, 2009). However, at the same time, brownification and increased sediment resuspension is already leading to light inhibition in spring, which may lead to a delayed spring bloom (Opdal et al., 2019). The spring bloom typically consists of chain-forming diatoms and is terminated by Si or N limitation (Eilertsen et al., 1989; Iversen and Seuthe, 2011). Zooplankton grazing is typically of low importance for terminating blooms (e.g., Saiz et al., 2013), while inorganic nutrients are considered to drive bloom termination (Krause et al., 2019, Mills et al., 2018). Heterotrophic bacteria remineralization of organic matter may supply additional N and Si (Legendre and Rassoulzadegan, 1995; Bidle and Azam, 1999; Johnson et al., 2007). N regeneration has been described as a mostly bacteria related process (Legendre and Rassoulzadegan, 1995), while Si dissolution is mainly controlled by abiotic dissolution of silica (Bidle and Azam, 1999). Zooplankton may also release some ammonium and urea after feeding on phytoplankton, but we suggest that this process is likely far less important than bacterial regeneration (e.g., Saiz et al., 2013). Previously measured ammonium excretion of Arctic mesozooplankton is typically low compared to bacterial remineralization (Conover and Gustavson, 1999), with the exception of one study in summer in a more open ocean setting (Alcaraz et al., 2010). In some Arctic systems urea excreted by zooplankton may be an important N source for regenerated algae production (Conover and Gustavson, 1999). A warmer climate will increase both bacteria-related remineralization rates (Legendre and Rassoulzadegan, 1995; Lannuzel et al., 2020) and abiotic silica dissolution (Bidle and Azam, 1999). However, the magnitude is not well understood.

Phytoplankton blooms may be dominated by a single or a few algal species, often with a similar physiology during certain phases of the bloom (e.g., Eilertsen et al., 1989; Degerlund and Eilertsen, 2010; Iversen and Seuthe, 2011). Chain-forming centric diatoms share physiological needs and responses to nutrient limitations (e.g., Eilertsen et al., 1989; von Quillfeldt, 2005) and typically dominate these blooms. In some Arctic and sub-Arctic areas the Arctic phytoplankton species chosen for this model, Chaetoceros socialis, can be dominant during spring blooms (Rey and Skjoldal, 1987; Eilertsen et al., 1989; Booth et al., 2002; Ratkova and Wassmann, 2002; von Quillfeldt, 2005; Degerlund and Eilertsen, 2010). Such spring phytoplankton blooms are accompanied by heterotrophic bacterioplankton blooms also showing typical succession patterns and distinct re-occurring taxa that dominate the community (Teeling et al., 2012, 2016). The importance of bacterial nutrient recycling for regenerated production has been recognized in several ecosystem models (e.g., van der Meersche et al., 2004; Vichi et al., 2007; Weitz et al., 2015) and algae bioreactor models focusing on nutrient conversions (e.g., Zambrano et al., 2016). However, these models are typically highly simplified or omitted in more sophisticated dynamic multi-nutrient, quota-based models (e.g., Flynn and Fasham, 1997; Wassmann et al., 2006; Ross and Geider, 2009). These latter models have been often developed and tuned based on cultivation experiments in which microbial remineralization reactions were assumed to be absent (e.g., Geider et al., 1998; Flynn, 2001) despite the fact that most algae cultures, likely including Geider et al. (1998) and Flynn (2001), are not axenic. Parameters estimated by fitting axenic models on non-axenic experiments may be misleading, mostly due to an inflated efficiency of dissolved inorganic nitrogen (DIN) uptake. Additional positive effects of bacteria include vitamin synthesis (Amin et al., 2012), trace metal chelation (Amin et al., 2012), the scavenging of oxidative stressors (Hünken et al., 2008), and exchange of growth factors (Amin et al., 2015). Especially in the stationary algal growth phase, Christie-Oleza et al. (2017) found that marine phototrophic cyanobacteria cultures are dependent on heterotrophic bacteria contaminants mainly due to their importance in degrading potentially toxic dissolved organic matter (DOM) exudates and regenerating ammonium. The current study aimed to bridge the gap between detailed representations of algae physiology and the role of microbial activity in an accurate way while keeping model complexity low.

Most ecosystem models consider only a single limiting nutrient to control primary production following Liebig's law of the minimum (Wassmann et al., 2006; Vichi et al., 2007). Yet we know that nutrient co-limitation is more complex. For example, ammonium and glutamate can inhibit nitrate uptake (Morris, 1974; Dortch, 1990; Flynn et al., 1997); C and N uptake is reduced under Fe limitation, while Si uptake continues (Werner, 1977; Firme et al., 2003); and the effects on photosynthesis differ for nitrogen and silicon limitations and for different algal groups (Werner, 1977; Flynn, 2003; Hohn et al., 2009). Complex interaction models considering intracellular biochemistry (NH4–NO3 co-limitation, Flynn et al., 1997), transporter densities and mobility (Flynn et al., 2018), and cell cycles (Si limitation, Flynn, 2001) can accurately describe these dynamics (Flynn, 2003) but are ultimately too computationally expensive to be integrated and parameterized in large-scale ecosystem models. Some models (Hohn et al., 2009; Le Quéré et al., 2016) implemented multi-nutrient (Hohn et al., 2009) and heterotrophic bacterial dynamics (Le Quéré et al., 2016) in Southern Ocean ecosystem models but have their limitations in representing bacterial remineralization (Hohn et al., 2009) or ammonium and silicate co-limitations (Le Quéré et al., 2016). In contrast to Antarctica, DIN is the primary limiting nutrient for phytoplankton growth, while iron is not limiting in most Arctic systems (Tremblay and Gagnon, 2009; Moore et al., 2013).

While simple lab experiments cannot represent all nutrient dynamics found in the environment (e.g., N excretion by zooplankton), they can focus on the quantitatively most important dynamics, to facilitate the development of simple multi-nutrient models, which are scalable to larger ecosystem models. The present study investigated the impact of silicate, ammonium–nitrate co-limitation and bacterial nutrient regeneration on photosynthesis, nitrogen assimilation, and cellular quotas based on data from a culture-based Arctic spring bloom system. The culture consisted of an axenic isolate of Chaetoceros socialis, dominating a phytoplankton net haul of a Svalbard fjord. The cultivation experiment was conducted either under axenic conditions or after inoculation with mostly free-living bacteria cultures, isolated beforehand from the non-axenic culture. Parametrization and insights from these incubations were then used to develop and parameterize a simple carbon-quota-based dynamic model (based on Geider et al., 1998), aiming to keep the number of parameters and computational costs as low as possible to allow for its implementation in large-scale ecosystem models.

The aims of the study were (i) to study the bloom dynamics of a simplified Arctic coastal pelagic system in a culture experiment consisting of one Arctic diatom species and co-cultured bacteria, (ii) to develop a simple dynamic model representing the observed interactions, and (iii) to discuss the importance of more complex bloom dynamics for their accurate representation in ecosystem models.

We hypothesize that (i) bacterial regeneration extends a phytoplankton growth period and gross carbon fixation, (ii) diatoms continue photosynthesis under silicate limitation at a reduced rate if DIN is available, and (iii) cultivation experiments are powerful for understanding the major spring bloom dynamics.

2 Methods

2.1 Cultivation experiment

The most abundant phytoplankton species from a net haul (20 µm mesh size) in April 2017 in van Mijenfjorden (Svalbard) Chaetoceros socialis was isolated via the dilution isolation method (Andersen et al., 2005) on F/2 medium (Guillard, 1975). Bacteria were isolated on LB-medium (evaluated by Bertani, 2004) agar plates using the algae culture as inoculum and sequenced at GENEWIZ LLC using the Sanger method and standard 16S rRNA primers targeting the V1-V9 region (forward 5-AGAGTTTGATCCTGGCTCAG-3, reverse 5-ACGGCTACCTTGTTACGACTT-3) provided by GENEWIZ LLC for identification via blastn (Altschul et al., 1990). Two strains of Pseudoalteromonas elyakovii, a taxon previously isolated from the Arctic (Al Khudary et al., 2008) and known to degrade algae polysaccharides (Ma et al., 2008) and to excrete polymeric substances (Kim et al., 2016), were successfully isolated and used for the experiments. Before the start of the experiment, all bacteria in the algae culture were killed using a mixture of the antibiotics penicillin and streptomycin. The success was confirmed via incubation of the cultures on LB agar plates and bacterial counts after DAPI staining (Porter et al., 1980). The axenic cultures were diluted in fresh F/2 medium lacking nitrate addition (Guillard, 1975) using sterile filtered seawater of Tromsø sound (Norway) as basis. The algae cultures were transferred into 96 200 mL sterile cultivation bottles with three replicates for each treatment. Half of the incubations were inoculated with bacteria cultures (BAC+), while the other half was kept axenic (BAC). The cultures were incubated at 4 C and 100 µE m−2 s−1 of continuous light and mixed two to three times a day to keep the algae and bacteria in suspension. We ensured sterile conditions during the experiment by keeping the cultivation bottles closed until sampling. However, endospores may survive the antibiotic treatment in low numbers and start growing especially towards the end of the experiment. Over 16 d three axenic and three BAC+ bottles were sacrificed daily for measurements of chlorophyll a (Chl), particulate organic carbon (POC) and nitrogen (PON), bacterial and algal abundances, nutrients (nitrate, nitrite, ammonium, phosphate, silicate), dissolved organic carbon (DOC), and the maximum quantum yield (QY) of photosystem II (PSII) (Fv/Fm) as a measure of healthy photosystems. Due to technical problems not all replicates could be measured on all days, and an overview of replication is given in Supplement Table S2.

Chlorophyll a was extracted from a GF/F (50 mL filtered at 200 mbar) filter at 4 C for 12–24 h in 98 % methanol in the dark before measurement in a Turner Trilogy™ fluorometer (evaluated by Jacobsen and Rai, 1990). POC and PON were measured after filtration onto precombusted (4 h at 450 C) GF/F (Whatman) filters (50 mL filtered at 200 mbar), using a Flash 2000 elemental analyzer (Thermo Fisher Scientific, Waltham, MA, USA) and Euro elemental analyzer (HEKAtech) following the protocol by Pella and Colombo (1973) after removing inorganic carbon by fuming with saturated HCl in a desiccator. Bacteria were counted after fixation of a water sample for 3–4 h with 2 % formaldehyde (final concentration), filtration of 25 mL on a 0.2 µm pore size polycarbonate filter, washing with filtered seawater and ethanol, DAPI staining for 7 min following Porter et al. (1980), and embedding in Citifluor–Vectashield (3:1). Bacteria were counted in at least 20 grids under an epifluorescence microscope (Leica DM LB2, Leica Microsystems, Germany) at 10×100 magnification. In the same sample the average diameter of diatom cells at the start and end of the experiment was measured. Algae were counted in 2 mL wells under an inverted microscope (Zeiss Primovert, Carl Zeiss AG, Germany) at 20×10 magnification after gentle mixing of the cultivation bottle. Algae cells incorporated in biofilms after day 9 in the BAC+ cultures were counted after sonication in a sonication bath until all cells were in suspension. Nutrient and DOC samples were sterile filtered (0.2 µm) and stored at 20 C before measurements. Nutrients were measured in triplicates after using standard colorimetric methods on a nutrient analyzer (QuAAtro 39, SEAL Analytical, Germany) using the protocols no. Q-068-05 rev. 12 for nitrate (detection limit = 0.02 µmol L−1), no. Q-068-05 rev. 12 for nitrite (detection limit = 0.02 µmol L−1), no. Q-066-05 rev. 5 for silicate (detection limit = 0.07 µmol L−1), and no. Q-064-05 rev. 8 for phosphate (detection limit = 0.01 µmol L−1). The data were analyzed using the software AACE. The nutrient analyzer was calibrated with reference seawater (Ocean Scientific International Ltd., United Kingdom). Ammonium was measured manually using the colorimetric method following McCarthy et al. (1977) on a spectrophotometer (Shimadzu UV-1201, detection limit = 0.01 µmol L−1). Ammonium values > 100 µmol L−1 were removed as outliers caused by too high a filtration pressure. DOC was measured by high-temperature catalytic oxidation (HTCO) using a Shimadzu TOC-5000 total C analyzer following methods for seawater samples (Burdige and Homstead, 1994). The photosynthetic quantum yield was determined using an AquaPen-C AP 100 (Photon Systems Instruments, Czech Republic).

Certain factors, such as grazing, settling out of the euphotic zone, and bacterial and algae succession, were not included into the experimental setup to reduce complexity and focus on nutrient dynamics. Trace metals, phosphate, and vitamin B12 in coastal systems are assumed to be not limiting in Arctic coastal systems and were supplied in excess to the culture medium. Realistic pre-bloom DOC concentrations were present in the medium as it was prepared with sterilized seawater from the Fjord outside Tromsø before the onset of the spring bloom (March 2018).

The f ratio as indication for the importance of regenerated production (Eppley, 1981) was calculated based on the average PON fixation in the last 3 d of the experiment (Eq. C1). Here, nitrogen assimilation in the BAC culture was assumed to be based on new (nitrate-based) production, while fixation in the BAC+ experiment was assumed to also be based on regenerated (ammonium-based) production.

2.2 Model structure

This section outlines the overall model structure followed by a description of the chosen parametrization approach for each relevant process. Details regarding model equations are provided in the Appendix (Tables A6 and A7), and a schematic representation of the models is given in Fig. 1. We used a dynamic cell quota model by Geider et al. (1998) to describe the BAC experiment (G98). We then extended the G98 model to represent the role of silicate limitation, bacterial regeneration of ammonium, and different kinetics for ammonium and nitrate uptake (EXT) and fitted it to the BAC+ experiment while retaining the parameter values previously estimated for G98.

Figure 1Schematic representation of the state variables and connections and controls in the G98 model (blue) and EXT model (purple). The EXT model has the same formulations as G98, with the additions shown in purple.


The Geider et al. (1998) model (G98) describes the response of phytoplankton to different nitrogen and light conditions and is based on both intracellular quotas and extracellular dissolved inorganic nitrogen (DIN) concentrations, allowing decoupled C and N growth (Fig. 1). Within this model, light is a controlling factor on photosynthesis and chlorophyll synthesis. C : N ratios and DIN concentrations control nitrogen assimilation, which is coupled to chlorophyll synthesis and photosynthesis. Chl : N ratios are controlling photosynthesis and chlorophyll synthesis. G98 has been used in a variety of large-scale ecosystem models, with some extensions representing the actual conditions in the environment or mesocosms (e.g., Moore et al., 2004; Schartau et al., 2007; Hauck et al., 2013). Photoacclimation dynamics in Geider-type models have been evaluated as quick and robust (Flynn et al., 2001), while the N assimilation component has some shortcomings in regard to ammonium–nitrate interactions. The original model of Geider et al. (1998) for C and N was corrected for minor typographical errors (see Ross and Geider, 2009; Appendix Tables A6, A7).

One aim of the study was to develop a model (EXT) with simplified dynamics of nutrient co-limitation, which is suitable for future implementation in coupled biogeochemistry–circulation models. The EXT model keeps all formulations of the G98 and adds dynamics and interactions of silicate, nitrate and ammonium uptake, carbon and nitrogen excretion, and bacterial remineralization (Fig. 1). The aim of the model was to describe the response in photosynthesis, chlorophyll synthesis, and nitrogen assimilation with a minimal number of parameters. Hence, dynamics in silicate cycling and bacterial physiology were highly simplified. The limitations of these simplifications and the potential need for more complex models are discussed later.

Silicate uptake was modeled using Monod kinetics following Spilling et al. (2010). The effect of silicate limitation on photosynthesis and chlorophyll synthesis was implemented after findings by Werner (1978), Martin-Jézéquel et al. (2000), and Claquin et al. (2002). Werner (1978) found that silicate limitation can lead to a 80 % reduction in photosynthesis and a stop of chlorophyll synthesis in diatoms within a few hours. Hence, we added a parameter for the reduction of photosynthesis under silicate limitation (SiPS) and formulated a stop of chlorophyll synthesis under silicate limitations.

N and Si metabolism have different controls and intracellular dynamics, with N uptake driven by photosynthesis (as PCref in G98) and Si mainly linked to algal respiration (Martin-Jezequel et al., 2000). Besides earlier cultivation studies, the reduction of photosynthesis after Si limitation has been shown via photophysiological (inhibited PSII reaction center, Lippemeier et al., 1999) and molecular (down-regulated photosynthetic proteins, Thangaraj et al., 2019) approaches. In general, we assume that nitrogen metabolism is not directly affected by silicate limitation (Hildebrand, 2002; Claquin et al., 2002), but we expect cellular ratios to be affected by reduced photosynthesis and chlorophyll synthesis under Si limitation (Hildebrand, 2002; Gilpin, 2004).

The algal respiration term included both respiration and excretion of dissolved organic nitrogen and carbon as a fraction of the carbon and nitrogen assimilated. For testing the importance of DON excretion, we also ran the EXT model without DON excretion (EXT-excr). Dissolved organic nitrogen (DON) was recycled into ammonium via bacterial remineralization. It was assumed that this process is faster for freshly excreted DON compared to DON already present in the medium. Thus, we implemented a labile DON pool (DONl) for freshly excreted DON and a refractory (DONr) DON pool with the respective remineralization rates rem and remd. We also assumed that excreted DON and DOC do not coagulate as extracellular polymeric substances (EPS) during the course of the experiment. Following Tezuka (1989), net bacterial regeneration of ammonium occurs at DOM C : N molar ratios below 10 and is proportional to bacterial abundances. Higher thresholds up to 29 have been found (e.g., Kirchman, 2000), but we selected a lower number to stay conservative. DOM C : N ratios are assumed to be proportional to algae C : N ratios (van der Meersche et al., 2004), with algal molar C : N ratios below 10 representing substrate (DOM) molar C : N ratios below 10.5. Hence, we assumed net bacterial ammonium regeneration to occur at molar POC / PON ratios below 10, while higher ratios lead to bacteria retaining more N for growth than they release. Bacteria abundance change was estimated using a simple logistic growth curve as a function of DOM since the number of parameters is low (2) and the fit sufficient for the purpose of modeling algae physiology. Michaelis–Menten kinetics based on bacteria growth on DOM with different labilities could give a more accurate representation of bacterial growth but would not change the fit of the other model parameters aiming for the best fit of the model output to algal PON, POC, Chl, and DIN. Algal nitrate uptake was modeled after the original model by Geider et al. (1998), and ammonium assimilation was based on the simplified SHANIM model by Flynn and Fasham (1997), excluding the internal nutrient and glutamine concentrations. Ammonium uptake is preferred over nitrate (lower half-saturation constant) and reduces nitrate assimilation if available above a certain threshold concentration of ammonium (Dortch, 1990; Flynn, 1999). Ammonium is the primary product of bacterial regeneration after remineralization of DON. Nitrification was assumed to be absent, since the bacteria in our experiment are not known to be capable of nitrification.

2.3 Model fitting

The model was based on a set of ordinary differential equations (ODEs) and was written in R. All model equations are provided in the Appendix (Tables A6 and A7), the newest version of the R code is available on GitHub (, last access: 29 January 2021), and the version used in this paper is archived at zenodo ( The ODEs were solved using the ode function of the deSolve package (Soetaert et al., 2010) with the second- to third-order Runge–Kutta method with automated step size control.

The parameters of the G98 model were fitted to the BAC experiment data, and those of the EXT model were fitted to the BAC+ experiment data. The model fitting started with data from day 1 in order to avoid artifacts during acclimation of the cultures after transfer to a new medium. Si and NOx were not measured at day 1, and the mean of day 0 and day 2 was used. The G98 parameter values were fitted first and retained without changes for the EXT model fitting. The maximum Chl : N ratio (θmaxN), minimum and maximum N : C ratios (Qmin, Qmax), and irradiance (I) were available as experimental data and needed no further fitting (Appendix Table A2). The start values and constraints for the remaining six variables (ζ, RC, αChl, n, KNO3, PrefC, Table A3) were based on model fits of G98 to other diatom culture in previous studies (Geider et al., 1998; Ross and Geider, 2009). The parameters were first fitted manually via graphical comparisons with the experimental data (POC, PON, Chl, DIN, Figs. 5 and 6) and via minimizing the model cost calculated as the root of the sum of squares normalized by dividing the squares with the variance (RMSE Eq. C2, Stow et al., 2009). The initial manual tuning approach allowed control of the model dynamics, considering potential problems with known limitations of the G98 model (e.g., lag phase not modeled; Pahlow, 2005). The manual tuning also allowed obtaining good start parameters for the automated tuning approach and sensitivity/collinearity analyses, which are sensitive to the start parameters.

After the manual tuning, an automated tuning approach was used to optimize the fits. The automated tuning was done using the FME package (Soetaert and Petzoldt, 2010), a package commonly used for fitting dynamic and inverse models based on differential equations (i.e., deSolve) to measured data. The automated analyses were based on minimizing the model cost calculated as the sum of squares of the residuals (SSR, fitted vs. measured data). The experimental data were normalized so that all normalized data were in a similar absolute range of values. This involved increasing Chl and PON values by an order of magnitude while decreasing DIN (NH4+ NO3) data by 1 order of magnitude. The data were not weighted, assuming equal data quality and importance. Prior to the automated fitting, parameters were tested for local sensitivity (sensFun) and collinearity or parameter identifiability (collin; e.g., Wu et al., 2014). sensFun tests for changes in output variables at each time point based on local perturbations of the model parameter. The sensitivity is calculated as L1 and L2 norms (Soetaert et al., 2009, 2010). The sensFun output is further used as input for the collinearity or parameter identifiability analyses. Parameters were considered collinear and not identifiable in combination with a collinearity index higher than 20 (Brun et al., 2001). In this case, only the more sensitive parameter was used for further tuning. Eventually, RC, KNO3, n, and αChl were subject to the automated tuning approach using the modFit function, based on minimizing the SSR within the given constraints. Parameters were first fitted using a pseudorandom search algorithm (Price, 1977) to ensure a global optimum. The resulting parameters were then fine-tuned using the Nelder–Mead algorithm (Soetaert et al., 2010) for finding a local optimum. A model run with the new parameters was then compared to the initial model via graphical comparisons of the model fit to the experimental data and via the RMSE value.

The parameter values obtained for the G98 fit to the BAC experiment were retained without changes or further fitting in the EXT model. The additional parameters of the EXT model were then fitted to the BAC+ experimental data (POC, Chl, PON, DIN). The model was only fitted to total DIN due to the potential uncertainties related to ammonium immobilization in the biofilm, which could be released during filtration and be part of the measured data. In fact, in a test run, fitting the EXT model to NO3 and NH4 separately led to a substantially worse overall fit (RMSE = 3.49). Otherwise, the data were not weighted. Since the aim of the study was to model the effects of silicate and bacteria on algae growth and not to develop an accurate model for bacteria biomass and silicate concentrations, the parameters μbact, bactmax, KSi, and Vmax were only fitted to the corresponding data (bacteria, silicate) prior to fitting the other parameters of the EXT model. Bacterial growth parameters (μbact, bactmax) were fitted to the bacterial growth curve using common bacterial carbon conversion units (20 fg C per cell as described by Lee and Fuhrman, 1987). Silicate-related parameters (KSi, Vmax) were constrained by the study of Werner (1978) and fitted to the measured silicate concentrations. The remaining parameters were subject to the tuning approach described for G98. Ammonium-related parameters (KNH4, NH4thres) were constrained by measured ammonium concentrations and constants available for other diatom taxa described by Eppley et al. (1969). Remineralization parameters for excreted (rem) and background (remd) DOM were constrained by the data with the limitation of rem > remd, assuming that the excreted DOM is more labile. The parameters related to the effect of silicate limitation on photosynthesis and chlorophyll production (smin, SiPS) were constrained by the study of Werner (1978) and fitted as described for G98. None of the added parameters were collinear/unidentifiable or given by the measured data and thus retained for the automated tuning approach. Eventually, the 15 parameters (Table A3) were fitted against 160 data points (Table A1).

Due to the biofilm formation in the stationary phase of the BAC+ experiment, we tested two additional modeling approaches representing different dynamics in biofilms: (i) DOC coagulation to EPS as part of the POC pool, which was assumed absent in our EXT model (Schartau et al., 2007), and (ii) increased DOM excretion in the stationary phase (e.g., Christie-Oleza et al., 2017). However, we suggest that the photosynthesis reduction term SiPS can give very similar model outputs, while being similarly or more sensitive. Thus, we tested the sensitivity of the added parameters of the two extended biofilm models in comparison to SiPS by testing the magnitude of perturbations of SiPS needed to reverse the effects of the added biofilm parameter (Figs. S1–3). The effects could be mostly reversed with similar or lower perturbations of SiPS. However, DOC coagulation to EPS can yield in a better overall better fit than SiPS alone. The main effect of the biofilm that we could not model with the available data appears to be ammonium immobilization in the biofilm, potentially by adsorption, accumulation in pockets, or conversion to ammonia due to the locally reduced pH caused by increased bacterial respiration. Model stability was estimated by extending the model run for 120 d to test for unrealistic model dynamics (Fig. S3).

3 Results

3.1 Cultivation experiment

The concentrations of nitrate and silicate declined rapidly over the course of the experiment (Fig. 2). After 8 d, silicate decreased to concentrations below 2 µmol L−1 – a threshold known to limit diatom dominance in phytoplankton (Egge and Aksnes, 1992), while inorganic nitrogen (nitrate, nitrite, and ammonium) became limiting (< 0.5 µmol L−1, POC : PON > 8–9 DIN : DIP < 16) only in the BAC culture. DIN : DIP ratios far below 16, or DIN concentrations below 2 µmol L−1, have been described as indication for DIN limitation (Pedersen and Borum, 1996), as well as POC : PON ratios > 9 (Geider and La Roche, 2002). Phosphate was not potentially growth limiting, with molar DIN-to-PO4 ratios consistently far below 16 (Redfield, 1934) and concentrations around 15 µmol L−1. Typically, phosphate concentrations below 0.3 µmol L−1 are considered limiting (e.g., Haecky and Andersson, 1999). Regeneration of ammonium and phosphate were important after 8 d as seen by increasing concentrations of both nutrients and showed higher concentrations in the BAC+ experiments compared to the BAC cultures (Fig. 2a, b). Ammonium concentrations were consistently higher, and nitrate was removed more slowly in the presence of bacteria, especially during the exponential phase. With the onset of the stationary phase in the BAC+ experiment, PO4 and NH4 concentrations doubled within 2 to 4 d and stayed high with variations in phosphate concentrations, while they stayed low in BAC. With depletion of NO3 in BAC+, NH4 concentrations remained high, while PO4 concentrations dropped. While not all ammonium measured is also available for algae growth, discussion of the dynamics (decrease in the start, increase with the onset of the stationary phase), especially if also shown in the EXT model, is still useful to understand multi-nutrient dynamics (e.g., regeneration). Considering the overall higher concentrations of NO3, compared to NH4, discussions of total DIN dynamics, DIN : DIP ratios, and limitations are also meaningful. DOC values were very high from the start (approx. 2–4 mmol L−1) and remained largely constant throughout the experiment (Table A8).

Figure 2Nutrient measurements over the experimental incubations of (a) NOx, (NO3-+ NO2-), (b) NH4+, (c) PO42- (the asterisk indicates a presumed measurement error), (d) silicate; blue circles are BAC cultures and red symbols are BAC+ cultures. Circles show median values (blue = BAC, red = BAC+) and the colored polygons show maximum and minimum of measured data (n=1–3, Table S2). The grey line shows the beginning of the stationary growth phase of Chaetoceros socialis.


The diatom Chaetoceros socialis grew exponentially in both treatments until day 8 before reaching a stationary phase with declining cell numbers (Fig. 3). The growth rate of the BAC culture (0.36 d−1) was slightly lower than in the treatment with bacteria present (0.42 d−1) during the exponential phase. Algal cellular abundance was higher in the BAC+ cultures. Towards the end of the exponential phase, the diatom started to form noticeable aggregates in cultures with bacteria present, but only to a limited extent in the BAC cultures. Such aggregate formation with associated EPS production is typical for C. socialis. With the onset of the stationary phase in the BAC+ cultures about 30 % of the cells formed biofilms on the walls of the cultivation bottles (estimated after sonication treatment). Bacteria (Fig. 3) continued to grow throughout the entire experiment, but growth rates slowed down from 0.9 to 0.6 after day 8. In the BAC cultures, bacterial numbers increased after 8 d, but abundances remained 2 orders of magnitude below the BAC+ cultures and effectively axenic over the experimental incubation period. The maximum photosynthetic quantum yield (Fv/Fm) is commonly used as a proxy of photosynthetic fitness (high QY), indicating the efficiency of energy transfer after adsorption in photosystem II. Low values are typically related to stress, including for example nitrogen (Cleveland and Perry, 1987) or silicate (Lippemeier et al., 1999) limitation. We found an increase in QY from approx. 0.62 to 0.67 in the exponential phase and a decrease to approx. 0.62 in the BAC+ treatment after 8 d and to approx. 0.58 in the BAC treatment (Table A8).

Figure 3Abundances of (a) Chaetoceros socialis and (b) bacteria over the 14 d experimental period. Blue data are from BAC cultures and red from BAC+ cultures. The asterisks at day 14 indicate potential outliers as based on only one replicate. Circles represent median values (blue = BAC, red = BAC+) and the colored polygons show maximum and minimum of measured data (n=1–3, Table S2, not visible for bacteria counts in BAC cultures due to very small range). The maximum values of the BAC+ experiment include algae cells in the biofilm (after day 9).


During algal exponential growth, POC and PON concentrations followed changes in algal abundances, increasing 4-, 7-, and 19-fold, respectively, within 8 d (Figs. 3a, 4). Interestingly, with the beginning of the stationary phase, POC and PON continued to increase in the BAC+ cultures, while their concentrations stayed constant (POC) or decreased due to maintenance respiration (PON) in BAC cultures. POC and PON concentrations were consistently higher (1.2 × POC, 1.4 × PON) in BAC+ cultures during the exponential phase. POC : PON ratios decreased in both cultures, but increased again after 11 d in the BAC culture. Chlorophyll a concentrations also increased exponentially over the first 8 d in both treatments, and thereafter they decreased within the stationary phase in the BAC cultures. In contrast, cell numbers remained nearly constant in the BAC+ cultures, before declining at the last sampling day.

Figure 4Total particulate organic (a) carbon (POC) (b) nitrogen (PON), (c) C : N ratios, and (d) chlorophyll a concentration in experimental cultures (the asterisk indicates a presumed measurement error). Blue symbols are BAC cultures and red show BAC+ cultures. Circles show median values (blue = BAC, red = BAC+) and the colored polygons show the maximum and minimum of measured data (n=2–3, Table S2). The grey line indicates the start of the stationary phase.


Overall, both experimental cultures showed similar growth dynamics until day 8, with silicate becoming limiting for both treatments and nitrogen only limiting in BAC cultures. Algal growth with bacteria present was slightly but consistently higher during this phase. After 8 d, algae growth stopped in both treatments, but nitrogen and carbon were continuously assimilated in BAC+ cultures. BAC cultures started to degrade chlorophyll, while it stayed the same in BAC+ cultures. Algal abundances in the BAC+ treatment at the end of the experiment were ca. 30 % higher due to biofilm formation, and considerably more carbon (2× total POC, or 10 %–20 % per cell) and nitrogen (3× total PON) per cell was assimilated and considerably more chlorophyll (2–3× total chlorophyll) was produced at day 16. Cell size differences were not detectable (ca. 4 µm diameter, Table A8). POC-to-PON ratios increased after 11 d in BAC cultures but showed no change in BAC+ cultures. POC-to-Chl ratios were comparable in both treatments (Fig. 5). Assuming BAC N fixation was mostly based on new production (nitrate as N source), while the algal N fixation in bacterial enriched treatments was based on new and regenerated (ammonium as N source) production, two-thirds of the production was based on regenerated production (f ratio = 0.31).

3.2 Modeling

A comparison of the traditional G98 model with the EXT model allowed an estimate of importance of bacterial DIN regeneration and Si co-limitations for describing the experimental growth dynamics. The EXT model led to a slightly worse fit to the BAC experiment (RMSEG98=2.79 RMSEEXT=3.38, Figs. 5 and 6). The real strength of the EXT model was in representing growth dynamics with bacteria present (Figs. 5 and 6). Here, the RMSE was reduced by 47 % from RMSEG98=4.31 down to RMSEEXT=2.31.

Figure 5Model fit of the EXT model to the BAC (blue) and BAC+ (red) experiment. Circles show median values, and the colored polygons show the maximum and minimum of measured data (n=2–3, Table S2). Solid lines show the model outputs of (a) POC, (b) PON, (c) Chl (the asterisk indicates a presumed measurement error) , (d) C : N, (e) C : Chl, and (f) N : Chl. Dotted lines show the model fit without additional carbon excretion term xf. At day 8 the threshold for silicate limitation is reached, leading to reduced photosynthesis (by the factor given by SiPS) and inhibited Chl synthesis, which is visible as sharp transitions in POC and Chl.


Figure 6Model fit of the EXT model to the BAC (blue) and BAC+ (red) experiment. Circles show median values, and the colored polygons show the maximum and minimum of measured data (n=1–3, Table S2). Solid lines show the model outputs of (a) DIN (NOX and NH4), (b) NOx, (c) NH4, and (d) Si(OH)4 (all model fits overlap).


Both the G98 and EXT model fits of the BAC experiment were similarly good for POC and PON with a slightly lower modeled growth rate. PON in the BAC+ experiment was poorly modeled without consideration of silicate limitation or regenerated production specifically towards the end of the exponential phase and during the stationary phase. Maximum PON values were about 3 times lower using the G98 model (Fig. B3). In addition, the start of the stationary phase in the BAC+ experiment was estimated 3 d too late via G98, even though modeled DIN was depleted 3 d too soon (Fig. B3). Under BAC conditions, where silicate limitation does not play a major role the G98 model appears sufficient.

The EXT model allowed representing detailed dynamics in a bacteria-influenced system such as the responses to silicate limitation with a decrease in POC production, continued PON production, and the stagnation of Chl synthesis (Fig. 5). Apart from the lag phase, the mass ratios of gC : gN and gC : gChl were represented accurately (Fig. 5). The model fits of POC, PON and Chl without the separate carbon excretion term (xf) were lower compared to the model with excretion, indicating the importance of excreted dissolved organic matter (DOM) concentrations (RMSEEXT-exr of 5.73).

DIN dynamics caused by ammonium–nitrate interactions were represented well (Fig. 6a). However, at the onset of the stationary phase, ammonium concentrations of the model were 1 order of magnitude lower than in the experiment, showing a major weakness (Fig. 6c). Increased weighting of ammonium during the model fitting led to a slightly better fit to ammonium but a substantially worse fit of the model to POC, PON, and Chl (RMSEEXT=3.49). This indicates that the problem lies with the ammonium data, which include immobilized ammonium in the biofilm released during the filtration but unavailable for diatom growth. Other potential differences in biofilms were tested by means of different model extensions (DOC and DON aggregation to EPS, increase in DOM excretion). After including the DOC and DOM aggregation to the model, the overall fit was improved (RMSEEXT=2.31, RMSEEXT+eps=2.19, Fig. S1). However, in the absence of EPS data, we used the EXT model for the main discussion. Increased DOM excretion could be explained by the SiPS term of the EXT model (Fig. S2). The silicate uptake estimation was highly simplified using simple Monod kinetics, leading to modeled values that are too high in the stationary phase and a depletion that is too quick at the start (Fig. 6d). Carbon excretion (xf) led to NOx depletion after 8 d.

The sensitivity analysis (Table A4) revealed that the sensitivity of the added parameters in EXT is overall comparable to the sensitivity of the original parameters in G98. The model outputs were most sensitive to PCref (L1 = 1.27, L2 = 2.05), which is a parameter in both G98 and EXT. The most sensitive added parameters in EXT were the DON excretion rate (xf, L1 = 0.17, L2 = 0.24), the bacterial growth rate (μbact, L1 = 0.18, L2 = 0.36), the remineralization rate of refractory DON (remd, L1 = 0.09, L2 = 0.21), and the inhibition of photosynthesis under Si limitation (SiPS, L1 = 0.07, L2 = 0.18), which was overall comparable to other sensitive parameters of the G98 model (Qmax, RC, αChl, ζ, n, I, ΘNmax, Table A1). Small perturbations of the parameters only indirectly related to the fitted output variables did not lead to changes in POC, PON, Chl, or DIN.

4 Discussion

The experimental incubations showed that in the presence of bacteria both the growth period and gross carbon fixation were extended (Hypothesis I). The diatoms were able to continue photosynthesis under silicate limitation at a reduced rate as long as inorganic nitrogen was present (Hypothesis II). Overall, the experimental incubations represented typical spring bloom dynamics for coastal Arctic systems, including an initial exponential growth phase terminated by N and Si limitation (Hypothesis III) and the potential for an extended growth period via regenerated production. Our model incorporating these results was able to reflect these dynamics by adding NH4–NO3–Si(OH)4 co-limitations and bacterial NH4 regeneration to the widely used G98 model. In addition, bacteria–algae interactions, DOC, and biofilm dynamics were important in the experiment, but those were not crucial for quantitatively modeling algal C : N : Chl quotas. While C. socialis may not be the dominant species in all coastal Arctic phytoplankton blooms, we argue that it is representative for chain-forming diatoms typically dominating these systems due to their shared needs and responses to nutrient limitations (e.g., Eilertsen et al., 1989; von Quillfeldt, 2005).

4.1 Silicon–nitrogen regeneration

Spring phytoplankton dynamics in Arctic and sub-Arctic coastal areas is typically characterized by an initial exponential growth of diatoms, followed by peaks of other taxa (like Phaeocystis pouchetii) soon after the onset of silicate limitation (Eilertsen et al., 1989). Thus, a shift in species composition for the secondary bloom is linked to silicate limitation prior to final bloom termination caused by inorganic nitrogen limitation. As suggested by our second hypothesis, photosynthesis was reduced by approx. 70 % after silicate became limiting, which is comparable to earlier experimental studies (Tezuka, 1989). However, as suggested by our first hypothesis, the secondary bloom was extended in time by bacterial regeneration of ammonium, allowing regenerated production to contribute about 69 % of the total production (f ratio = 0.31) even during the diatom-dominated scenario in our experimental incubation. With the start of the stationary phase, NH4 and PO4 concentrations doubled, presumably due to decreased assimilation by the silicate-starved diatoms and increased regeneration by bacteria, supplied with increasing labile DOM (doubled remineralization rate in EXT) excreted by the stressed algae. However, NH4 concentrations doubled within 4 d, while PO4 concentrations doubled in only 1 d, indicating some unexplained internal dynamics, potentially via different bacterial uptake and release of N and P. After NO3 depletion at day 15, the PO4 concentrations also dropped, indicating a coupling of N : P metabolism but not of NH4 : P metabolism. Thus, the sudden drop may also indicate dynamics of bottle experiments, which are not accounted for, showing potential limitations of these experiments. The presence of bacteria and thus regenerated production allowed diatom growth to continue 8 d after silicate became limiting (Figs. 2, 3 and 4), nearly doubling the growth period similar to observations in the field (e.g., Legendre and Rassoulzadegan, 1995; Johnson et al., 2007), which supports our hypotheses I and III.

The G98 model has its most severe limitation in the modeling of PON, simply due to the lack of the ammonium pool, supplied via bacterial regeneration. The substantially better fit of PON in the EXT model therefore clearly shows that bacterial remineralization is crucial to successfully modeling spring bloom dynamics, especially near bloom termination. Many biogeochemical models used in the Arctic include remineralization but rely on fixed or temperature-dependent rates and do not consider them bacteria-dependent (MEDUSA, LANL, NEMURO, NPZD; see Table 1). While this simplification allows modeling regenerated production, using bacteria-independent remineralization rates does have limitations under spring bloom scenarios, where bacteria biomass can vary over orders of magnitudes (e.g., Sturluson et al., 2008) as also seen in our experimental study.

Table 1A comparison of major components contributing to the complexity of different models discussed. “No. of param” is the number of parameters. In the case of ecosystem models (SINMOD, BFM, MEDUSA, LANL, NEMURO, NPZD) only the model formulations representing the components of the current model (phytoplankton growth, remineralization, nutrient dynamics) are considered. For the full ecosystem-scale models we give the original reference to the biogeochemical compartment of the ecosystem-scale models and examples for more recent versions with updated formulations of other model compartments (e.g., physical drivers). REM designates those models that include remineralization (Rem); marked with “V” is present and “–” is absent. Ratios shows if the stoichiometry in the model considers variable or fixed ratios of intracellular elements (C : N : Si : P : Fe). The nutrients considered are given under the nutrients column. If DIN is considered as both NH4 and NO3, N is shown as N2. MEDUSA has Fe-dependent Si : N ratios, which makes them fixed in the Arctic (entry fixed* in the Table).

Degrees of freedom after constraints by the measured data are ∗114 and ∗26.

Download Print Version | Download XLSX

While we do not expect the f ratio in our bottle experiment to be directly comparable to an open-ocean system due to the higher biodiversity of field communities, a comparison can aid to identify limitations in our experiment and model. The f ratio can be used as a measure to check how representative the cultivation study was for typical spring bloom dynamics (Hypothesis III). Regenerated production is significant in polar systems, and our estimated experimental f value of 0.31 is slightly below the average for polar systems (Harrison and Cota, 1991, mean f ratio = 0.54). Nitrification is a process supplying about 50 % of the NO3 used for primary production in the oceans, which may lead to a substantial underestimation of regenerated production (Yool et al., 2007), inflating the f ratio interpreted as an estimate for new production, potentially also in the study by Harrison and Cota (1991). The absence of vertical PON export in our experiment may be another explanation for the above-average fraction of regenerated production. In the ocean environment, regenerated production is also affected by vertical export (sedimentation) and grazing, which are not represented in the experimental incubations. Via sedimentation, a fraction of the bloom either in the form of direct algal sinking of fecal pellets is typically exported to deeper water layers, reducing the potential for N regeneration within the euphotic zone (e.g., Keck and Wassmann, 1996). Larger zooplankton grazing can lead to increased export of PON via fecal pellet aggregation, or diel vertical migration (Banse, 1995), but may also release ammonium and urea (Conover and Gustavson, 1999; Saiz et al., 2013).

In contrast, bacterial death by microflagellate grazing and viral lysis may supply additional nutrients, or DON available for N regeneration in the euphotic zone (e.g., Goldman and Caron, 1985), which potentially leads to an overestimation of regenerated production. Another potentially important N source for regenerated production may be urea (Harrison et al., 1985), which would lead to an even higher importance of regenerated production as suggested by our study. Hence, ecosystem-scale models will need to consider these dynamics regarding bacterial abundances, microbial networks and particle export in addition to bacterial remineralization in order to model realistic ammonium regeneration in the euphotic zone. Overall, our cultivation experiment was powerful enough to represent major aspects of spring bloom dynamics but has its limitations, thereby confirming our third hypothesis only to some extent.

Bacteria-mediated silicate regeneration is absent from the modeling approach, as indicated by the identical silicate concentrations in both treatments and models (Fig. 2d). In the environment silicate dissolution is, in fact, mostly described as an abiotic process with temperature as the main control and a minor contribution by bacterial remineralization (Bidle and Azam, 1999). Our experiment indicates that silicate dissolution for Chaetoceros socialis was negligible at cold temperatures and the timescale of the incubations and typical for bloom durations and residence times of algae cells in the euphotic zone (Eilertsen et al., 1989; Keck and Wassmann, 1996). We conclude that silicate dissolution in coastal Arctic systems happens most likely in the sediment or deeper water layers and is only supplied via mixing in winter. In Antarctica substantial silicate dissolution has been observed but not in the upper 100 m, which has been related to the low temperatures (Nelson and Gordon, 1981) in agreement with our conclusion. Hence, modeling silicate regeneration in the euphotic zone is not necessary in these systems.

4.2 Algal growth response to Si and N limitation

The response of diatoms to Si or N limitation is based on different dynamics and different roles of N and Si in diatom growth. N is needed for proteins and nucleic acids, and its uptake is mainly driven by phototrophic reactions (Martin-Jézéquel et al., 2000). Si is only needed for frustule formation and cell division, mostly during a specific time in the cell cycle (G2 and M phase, Hildebrand, 2002), and the assimilation is mostly driven by heterotrophic reactions (Martin-Jézéquel et al., 2000). Once N is limiting, growth rapidly stops (Geider et al., 1998). In the case of Si limitation, however, growth can continue with a slower rate if N is still available (Werner, 1978; Gilpin et al., 2004). Several studies found a reduced growth rate with weaker silicified cell walls (Hildebrand, 2002; Gilpin, 2004) but unaffected nitrogen assimilation under silicate limitation (Hildebrand, 2002; Claquin et al., 2002) in accordance with our study. Claquin et al. (2002) found variable Si : C and Si : N ratios and highly silicified cells under nitrogen limitation, indicating uncoupled Si and N : C metabolism.

Nitrogen is a crucial element as part of amino acids and nucleic acids, which are necessary for cell activity and growth. If N becomes limiting, major cellular processes are affected and growth or chlorophyll synthesis is not possible. Photosynthesis can continue for a while, leading to carbon overconsumption (Schartau et al., 2007), which is well modeled by G98 for both BAC+ and BAC. A part of the excess carbon can be stored as intracellular reserves, and a part is excreted as DOC, which may aggregate as EPS, contributing to the total POC pool. The excess carbon can potentially be toxic for the algae and excretion, and extracellular degradation by bacteria may be crucial for algal survival (Christie-Oleza et al., 2017). Quantitatively, N limitation is well modeled by G98 under BAC conditions, if only one nitrogen source plays a role. However, under longer nitrogen starvation times or higher light intensities, alternative models that include carbon excretion and aggregation (Schartau et al., 2007) or intracellular storage in reserve pools (Ross and Geider, 2009) might be needed. Our growth experiment clearly shows that C : N ratios are not fixed, and variable quotas are needed. Vichi et al. (2007) estimated that carbon-based models that do not consider variable C : N ratios may underestimate net primary production (NPP) by 50 %, arguing for the importance of quota-based models (Fransner et al., 2017). However, most ecosystem-scale models are simplified by using fixed C : N ratios (Table 1). The next step towards quota-based models is the consideration of more detailed cell-based characteristics, such as transporter density, cell size, and mobility, including sedimentation (Aksnes and Egge, 1991). Flynn et al. (2018) discuss a model with detailed uptake kinetics showing that large cells are overall disfavored over small cells due to higher half-saturation constants but that they may still have competitive advantages due to lower investment in transporter production. Increased sedimentation in larger cells also increases the mobility and may offset the disadvantage of a larger size. While this extension is too complex for our aim of a simplified model, the dynamics may become important when modeling different algae taxa.

The type of inorganic nitrogen available also affects nitrogen uptake. Due to the metabolic costs related to intracellular nitrate reduction to ammonium, ammonium uptake is preferred over nitrate, potentially leaving more energy for other processes (Lachmann et al., 2019). Ammonium can even inhibit or reduce nitrate uptake over certain concentrations (Morris, 1974). The dynamics are mostly controlled by intracellular processes, such as glutamate feedbacks on nitrogen assimilation, cost for nitrate conversion to ammonium, or lower half-saturation constants of ammonium transporters (Flynn et al., 1997). The most accurate representation of these dynamics is given in the ANIM model by Flynn et al. (1997), but the model is too complex for implementations in larger ecosystem models. The number of parameters is difficult to tune with the typically limited availability of measured data, and its complexity makes it also computationally costly to scale up the models. Typically, modeling ammonium–nitrate interactions by different half-saturation constants and inhibition of nitrate uptake by ammonium appears sufficient (e.g., BFM, LANL, NEMURO, Table 1) and has been adapted in our model.

Studies on the coupling of silicate limitation on C, N, and Chl show inconclusive patterns, including a complete decoupling (Claquin et al., 2002), a relation of N to Si (Gilpin et al., 2004), and reduction of photosynthesis without new chlorophyll production (Werner, 1978; Gilpin et al., 2004). Cell size is limited by the frustules, but cells may become more nutritious (higher N : C ratio), or simply excrete more DOM, which may aggregate and contribute to the PON and POC pools. A detailed cell-cycle-based model has been suggested by Flynn (2001), but the number of parameters (30) makes the model too complex for ecosystem-scale models. In ecosystem-scale models Si limitation is modeled in various simplifications, such as thresholds triggering a stop (MEDUSA), reduction (e.g., BFM, MEDUSA, SINMOD) of the Si-dependent production (Table 1), or Si : N ratio scaled production (NEMURO, Table 1).

Our cultivation study shows (i) that a threshold value in the model, leading to a stop or solely Si-dependent photosynthesis, has its limitations, since DIN controlled photosynthesis continues at lower rates, and (ii) that coupling of Si : N : C : Chl is present. We do not expect a direct Si : N coupling, due to different controls of Si and N metabolism (Martin-Jézéquel et al., 2000) but suggest indirect coupling via reduced photosynthesis. In fact, detailed photophysiological and molecular approaches under Si limitation found inhibited PSII reaction centers (Lippemeier et al., 1999) similar to the decreased QY in our experiment and down-regulated photosynthetic proteins (Thangaraj et al., 2019) under Si limitation. Thus, we modeled the response of diatom growth to silicate limitation by reducing photosynthesis through a parameterized fraction (SiPS) and a stop of chlorophyll synthesis below a certain threshold, based on experimental studies (Werner, 1978; Lippemeier et al., 1999; Gilpin et al., 2004; Thangaraj et al., 2019) and in accordance to other ecosystem-scale approaches. We suggest that this extension is more accurate than the typical threshold-based dynamics, with one limiting nutrient controlling the growth equally for POC and Chl production (e.g., SINMOD by Wassmann et al., 2006; BFM by Vichi et al., 2007) while still keeping the number of parameters low compared to very detailed cell-cycle-based models (e.g., Flynn, 2001; Flynn et al., 2018).

4.3 Importance of algae–bacteria interactions and DOC excretion

As described above, N or Si limitation can lead to excretion of DON and DOC, which can aggregate as EPS and be available for bacterial regeneration of ammonium. For accurately including EPS dynamics in the model, additional data would be needed. However, the importance of EPS formation is evident in the end of the BAC+ experiment. Firstly, a biofilm was clearly visible, containing about 30 % of the algae cells. While we would not expect biofilms in the open ocean, aggregation of algae cells, facilitated by EPS, is common towards the end of spring blooms, increasing vertical export fluxes (e.g., Thornton, 2002). Chaetoceros socialis is in fact a colony-forming diatom building EPS-rich aggregates in nature (Booth et al., 2002). Secondly, POC and PON concentrations increased, while cell numbers and sizes stayed constant, showing that the additional POC and PON was most likely part of an extracellular pool. Silicate limitation could be one trigger for enhanced exudation. In fact, the biofilm dynamics evaluated (DOC aggregation, increased excretion) showed similar dynamics as the SiPS term. Since the biofilm formation corresponds with silicate limitation, it is difficult to untangle the direct effects of the biofilm or the indirect effects of silicate limitation without additional data or experiments (e.g., EPS measurements, DOM characterization). However, our model run, including EPS aggregation, allowed an improved overall fit while reducing the SiPS term to values more similar to earlier studies (Werner, 1978), pointing to the importance of DOC and DON aggregation. Hence, we suggest that Si limitation and EPS aggregation are both important.

Interestingly, algae–bacteria interactions can be species-specific, with specific organic molecules excreted by the algae to attract specific beneficial bacteria (Mühlenbruch et al., 2018). Thereby bacteria are crucial for recycling ammonium but also to degrade potentially toxic exudates (Christie-Oleza et al., 2017).

In the BAC experiment, carbon excretion after carbon overconsumption could be expected following Schartau et al. (2007), but no indications, such as biofilm formation or increased POC per cell, were found. In fact, a model considering EPS aggregation led to a substantially worse fit to BAC. This indicates that carbon overconsumption has been of minor importance, likely due to the low light levels. An alternative explanation is that bacteria and potentially chemotaxis are important controls on algal carbon excretion (Mühlenbruch et al., 2018). Overall, DOM excretion and EPS dynamics appear to play a major role in quantitatively modeling C : N : Chl quotas in our experiment, with higher errors (RMSEEXT-excr=5.73, RMSEEXT=2.31) for a model run without than with the excretion term xf.

4.4 Considerations in a changing climate

Due to a rapid changing climate, especially in Arctic coastal systems, the dynamics addressed in this study will change (Tremblay and Gagnon, 2009). With warmer temperatures, heterotrophic activities, and thereby bacterial recycling, will increase (Kirchman et al., 2009). Our study showed that regenerated production is crucial for an extended spring bloom. Hence, higher heterotrophic activities may lead to extended blooms (increased bacterial regeneration). At the same time, higher temperatures and increased precipitation will lead to stronger and earlier stratified water columns, which will lead to less nutrients reaching the surface by winter mixing, reducing new production (decreased bacterial regeneration) (Tremblay and Gagnon, 2009; Fu et al., 2016). Consequently, the phenology of Arctic coastal primary production in a warmer climate will likely be increasingly driven by bacterial remineralization, showing the necessity to include this process into biogeochemical models. An earlier temperature-driven water column stratification may also lead to an earlier bloom. However, due to increasing river and lake brownification and sediment resuspension, the spring bloom may also be delayed (Opdal et al., 2019). With decreased light, carbon overconsumption as described by Schartau et al. (2007) may become less important due to decreased photosynthesis. An earlier or later phytoplankton bloom can lead to a mismatch with zooplankton grazers (Durant et al., 2007; Sommer et al., 2007). Reduced zooplankton production would decrease the fecal-pellet-driven vertical export and thereby increase the residence time of particulate organic matter (POM) in the euphotic zone and the potential for ammonium regeneration. Thus, the incorporation of bacterial recycling into ecosystem models may be even more important under this scenario. In fact, global climate change models agree that vertical carbon export is decreasing overall (Fu et al., 2016). Silicate regeneration is thought to be mostly controlled abiotically by temperature (Bidle and Azam, 1999). Thus, increasing temperature and a stronger stratification will allow recycling of silicate in the euphotic zone before sinking out and thus could cause a shift in the algal succession observed during spring with prolonged contributions of diatoms (Kamatani, 1982). Thus, a temperature-dependent silica dissolution may need to be included for models in a substantially warmer climate in further model developments. Increased precipitation will also lead to increased runoff and allochthonous DOM inputs, increasing the importance of terrestrial DOM degradation and decreasing the relative importance of algal exudate regeneration (Jansson et al., 2008). The low light levels and the absence of grazing and export fluxes are simplifications of our study, which are, however, expected to be realistic scenarios under climate change. Hence, we suggest that our experiment and model are well suited as a baseline for predictive ecosystem models investigating the impacts of climate change on coastal Arctic spring blooms. However, climate change may lead to shifts in algae communities with non-silicifying algae dominating over diatoms (e.g., Falkowski and Oliver, 2007), reducing the importance of silicate limitation. Thus, conducting similar experiments and modeling exercises with a wider range of algal taxa and different temperature and nutrient regimes is suggested.

Appendix A

Table A1State variables of the G98 model and the EXT model (marked with “V” if present and “–” if absent) with units and designation if these state variables had been measured in the experiment.

Download Print Version | Download XLSX

Table A2Parameters of the original G98 model and the model extension with associated units.

Download Print Version | Download XLSX

Table A3Parameters and constants of the original G98 model and the EXT model with initial values used in the model and the lower and upper value constraints used for model fitting, unless the parameter was already defined by the data (measured). The constraints are either based on G98 fits to other diatom species, to present experimental data, or to typical values found in the literature.

Download Print Version | Download XLSX

Table A4Output of the sensitivity analysis (sensFun of the FME package in R, EXT fit on BACT+) with the value for each parameter and different sensitivity indices obtained after quantifying the effects of small perturbations of the parameters on the output variables (POC, PON, Chl, DIN). The L1 and L2 norms are normalized sensitivity indices defined as L1=Si,jn and L2=Si,j2n, with Si,j being the sensitivity of parameter i for model output j.

Download Print Version | Download XLSX

Table A5Other parameters calculated and used in the model equations.

Download Print Version | Download XLSX

Table A6Model equations from G98 (Geider et al., 1998) corrected for typographical errors by Ross and Geider (2009) with extensions. DIC stands for dissolved inorganic carbon.

Download Print Version | Download XLSX

Table A7Model equations of the EXT model based on G98.

Download XLSX

Table A8Output of the collinearity or parameter identifiability analysis using the collin function (G98 fit on BACT-data) of the FME R package (Soetaert et al., 2010). A subset of any combination of two parameters with a collinearity above 20 indicating non-identifiable parameter combinations is given (Brun et al., 2001). Parameter combinations tested are marked with a “V” in the left part of the table; the collinearity output is given on the right column.

Download Print Version | Download XLSX

Appendix B

Figure B1Model fit of the G98 model to the BAC (blue) and BAC+ (red) experiment. Circles show median values, and the colored polygons show the minimum and maximum of the measured data (n=2–3, Table S2). Solid lines show the model outputs of (a) POC, (b) PON, (c) Chl (the asterisk indicates a presumed measurement error), (d) C : N, (e) C : Chl, and (f) N : Chl.


Appendix C

The F ratio estimation in the cultivation experiments with the average PON concentrations at day 13 to 15 (PONd13−15) for the BAC and BAC+ treatments is calculated as follows.

(C1) f  ratio = PON BAC - d 13 - 15 PON BAC - d 13 - 15 + PON BAC + d 13 - 15

The normalized RMSE, with i being the different variables (POC, PON, Chl, DIN), j being the different values of each state variable, P being the predicted values, and O being observed values, is calculated as follows.

(C2) RMSE = i = 1 , j = 1 n , p P i , j - O i , j 2 Var O i
Data availability

The experimental data are archived at DataverseNO under the DOI number (, Vonnahme et al., 2021). The R scripts for the model used in this publication are archived at zenodo (, Vonnahme, 2021a) and available at GitHub under (Vonnahme, 2021b).


The supplement related to this article is available online at:

Author contributions

TRV designed the experiment with contributions by RG and ML. TRV isolated and identified the cultures. ML performed the experiment with contributions of TRV and UD. HRH measured DOC and SK measured the nutrients. The other parameters where measured by ML and TRV. TRV programmed the model with contributions from CV, ST, and DvO. TRV wrote the manuscript with contributions from all co-authors.

Competing interests

The authors declare that they have no conflict of interest.


The project was supported by Arctic SIZE – a research group on the productive marginal ice zone at UiT (project number 01vm/h15). We want to thank Paul Dubourg and Elzbieta Anna Petelenz-Kurdziel for the help with nutrient and POC : PON analyses. DOC analyses were supported through a Fulbright Distinguished Scholar Award to Rodger Harvey.

Financial support

This research has been jointly funded by UiT – The Arctic university of Norway and the Tromsø Research Foundation under the project Arctic SIZE – a research group on the productive Marginal Ice Zone, project number 01vm/h15.

The article processing charges for this open-access
publication were covered by a Research
Centre of the Helmholtz Association.

Review statement

This paper was edited by Kenneth Rose and reviewed by three anonymous referees.


Aksnes, D. L. and Egge, J. K.: A theoretical model for nutrient uptake in phytoplankton, Mar. Ecol. Prog. Ser., 70, 65–72, 1991. 

Alcaraz, M., Almeda, R., Calbet, A., Saiz, E., Duarte, C. M., Lasternas, S., Agusti, S., Santiago, R., Movilla, J., and Alonso, A.: The role of arctic zooplankton in biogeochemical cycles: respiration and excretion of ammonia and phosphate during summer, Polar Biol., 33, 1719–1731, 2010. 

Al Khudary, R., Stößer, N. I., Qoura, F., and Antranikian, G.: Pseudoalteromonas arctica sp. nov., an aerobic, psychrotolerant, marine bacterium isolated from Spitzbergen, Int. J. Syst. Evol. Microbiol., 58, 2018–2024, 2008. 

Altschul, S. F., Gish, W., Miller, W., Myers, E. W., and Lipman, D. J.: Basic local alignment search tool, J. Mol. Biol., 215, 403–410, 1990. 

Alver, M. O., Broch, O. J., Melle, W., Bagøien, E., and Slagstad, D.: Validation of an Eulerian population model for the marine copepod Calanus finmarchicus in the Norwegian Sea, J. Mar. Syst., 160, 81–93, 2016. 

Amin, S. A., Parker, M. S., and Armbrust, E. V.: Interactions between diatoms and bacteria, Microbiol. Mol. Biol. Rev., 76, 667–684, 2012. 

Amin, S. A., Hmelo, L. R., Van Tol, H. M., Durham, B. P., Carlson, L. T., Heal, K. R., Morales, R. L., Berthiaume, C. T., Parker, M. S., Djunaedi, B., Ingalls, A. E., Parsek, M. R., Moran, M. A., and Armbrust, E. V.: Interaction and signalling between a cosmopolitan phytoplankton and associated bacteria, Nature, 522, 98–101, 2015. 

Andersen, R. A. and Kawachi, M.: Microalgae isolation techniques, in: Algal culturing techniques, edited by: Andersen, R. A., Elsevier, San Diego, London, 83–100, 2005. 

Anju, M., Sreeush, M. G., Valsala, V., Smitha, B. R., Hamza, F., Bharathi, G., and Naidu, C. V.: Understanding the role of nutrient limitation on plankton biomass over Arabian Sea via 1-D coupled biogeochemical model and Bio-Argo observations, J. Geophys. Res.-Oceans, 125, 1–28,, 2020. 

Banse, K.: Zooplankton: pivotal role in the control of ocean production: I. Biomass and production, ICES J. Mar. Sci., 52, 265–277, 1995. 

Bertani, G.: Lysogeny at mid-twentieth century: P1, P2, and other experimental systems, J. Bacteriol., 186, 595–600, 2004. 

Bidle, K. D. and Azam, F.: Accelerated dissolution of diatom silica by marine bacterial assemblages, Nature, 397, 508–512, 1999. 

Booth, B. C., Larouche, P., Bélanger, S., Klein, B., Amiel, D., and Mei, Z. P.: Dynamics of Chaetoceros socialis blooms in the North Water, Deep Sea Res., 49, 5003–5025, 2002. 

Brun, R., Reichert, P., and Kunsch, H. R.: Practical identifiability analysis of large environmental simulation models, Water Resour. Res. 37, 1015–1030, 2001. 

Burdige, D. J. and Homstead, J.: Fluxes of dissolved organic carbon from Chesapeake Bay sediments, Geochim. Cosmochim. Ac., 58, 3407–3424, 1994. 

Christie-Oleza, J. A., Sousoni, D., Lloyd, M., Armengaud, J., and Scanlan, D. J.: Nutrient recycling facilitates long-term stability of marine microbial phototroph–heterotroph interactions, Nat. Microbiol., 2, 17100,, 2017. 

Claquin, P., Martin-Jézéquel, V., Kromkamp, J. C., Veldhuis, M. J. W., and Kraay, G. W.: Uncoupling of silicon compared with carbon and nitrogen metabolisms and the role of the cell cycle in continuous cultures of Thalassiosira pseudonana (Bacillariophyceae) under light, nitrogen, and phosphorus control, J. Phycol., 38, 922–930, 2002. 

Cleveland, J. S. and Perry, M. J.: Quantum yield, relative specific absorption and fluorescence in nitrogen-limited Chaetoceros gracilis, Mar. Biol., 94, 489–497, 1987. 

Conover, R. J. and Gustavson, K. R.: Sources of urea in arctic seas: zooplankton metabolism, Mar. Ecol. Prog. Ser., 179, 41–54, 1999. 

Degerlund, M. and Eilertsen, H. C.: Main species characteristics of phytoplankton spring blooms in NE Atlantic and Arctic waters (68–80 N), Estuar. Coast, 33, 242–269, 2010. 

Dortch, Q.: The interaction between ammonium and nitrate uptake in phytoplankton, Mar. Ecol. Prog. Ser., 61, 183–201, 1990. 

Durant, J. M., Hjermann, D. Ø., Ottersen, G., and Stenseth, N. C.: Climate and the match or mismatch between predator requirements and resource availability, Clim. Res., 33, 271–283, 2007. 

Egge, J. K. and Aksnes, D. L.: Silicate as regulating nutrient in phytoplankton competition, Mar. Ecol. Prog. ser., 83, 281–289, 1992. 

Eilertsen, H. C. and Frantzen, S.: Phytoplankton from two sub-Arctic fjords in northern Norway 2002–2004: I. Seasonal variations in chlorophyll a and bloom dynamics, Mar. Biol. Res., 3, 319–332, 2007. 

Eilertsen, H. C., Taasen, J. P., and WesIawski, J. M.: Phytoplankton studies in the fjords of West Spitzbergen: physical environment and production in spring and summer, J. Plankton Res., 11, 1245–1260, 1989. 

Eppley, R. W.: Autotrophic production of particulate matter, in: Analysis of marine ecosystems, edited by: Longhurst, A. R., Academic Press, New York, 343–361, 1981. 

Eppley, R. W., Rogers, J. N., and McCarthy, J. J.: Half-saturation constants for uptake of nitrate and ammonium by marine phytoplankton, Limnol. Oceanogr., 14, 912–920, 1969. 

Falkowski, P. G. and Oliver, M. J.: Mix and match: how climate selects phytoplankton, Nat. Rev. Microbiol., 5, 813–819, 2007. 

Field, C. B., Behrenfeld, M. J., Randerson, J. T., and Falkowski, P.: Primary production of the biosphere: integrating terrestrial and oceanic components, Science, 281, 237–240, 1998. 

Firme, G. F., Rue, E. L., Weeks, D. A., Bruland, K. W., and Hutchins, D. A.: Spatial and temporal variability in phytoplankton iron limitation along the California coast and consequences for Si, N, and C biogeochemistry, Global Biogeochem. Cy., 17, 1–13,, 2003. 

Flynn, K. J.: Nitrate transport and ammonium-nitrate interactions at high nitrate concentrations and low temperature, Mar. Ecol. Prog. Ser., 187, 283–287, 1999. 

Flynn, K. J.: A mechanistic model for describing dynamic multi-nutrient, light, temperature interactions in phytoplankton, J. Plankton Res., 23, 977–997, 2001. 

Flynn, K. J.: Modelling multi-nutrient interactions in phytoplankton; balancing simplicity and realism, Prog. Oceanogr., 56, 249–279, 2003. 

Flynn, K. J. and Fasham, M. J.: A short version of the ammonium-nitrate interaction model, J. Plankton Res., 19, 1881–1897, 1997. 

Flynn, K. J., Fasham, M. J., and Hipkin, C. R.: Modelling the interactions between ammonium and nitrate uptake in marine phytoplankton, Philos. T. R. Soc. A, 352, 1625–1645, 1997. 

Flynn, K. J., Marshall, H., and Geider, R. J.: A comparison of two N-irradiance interaction models of phytoplankton growth, Limnol. Oceanogr., 46, 1794–1802, 2001. 

Flynn, K. J., Skibinski, D. O., and Lindemann, C.: Effects of growth rate, cell size, motion, and elemental stoichiometry on nutrient transport kinetics, PLoS Comput. Biol., 14, e1006118,, 2018. 

Fransner, F., Gustafsson, E., Tedesco, L., Vichi, M., Hordoir, R., Roquet, F., Spilling, K., Kuznetsov, I., Eilola, K., Mörth, C., Humborg, C., and Nycander, J.: Non-Redfieldian dynamics explain seasonal pCO2 drawdown in the Gulf of Bothnia, J. Geophys. Res.-Oceans, 123, 166–188, 2017. 

Fritz, M., Vonk, J. E., and Lantuit, H.: Collapsing arctic coastlines, Nat. Clim. Chang, 7, 6–7,, 2017. 

Fu, W., Randerson, J. T., and Moore, J. K.: Climate change impacts on net primary production (NPP) and export production (EP) regulated by increasing stratification and phytoplankton community structure in the CMIP5 models, Biogeosciences, 13, 5151–5170,, 2016. 

Geider, R. and La Roche, J.: Redfield revisited: variability of C: N: P in marine microalgae and its biochemical basis, Eur. J. Phycol., 37, 1–17, 2002. 

Geider, R. J., Maclntyre, H. L., and Kana, T. M.: A dynamic regulatory model of phytoplanktonic acclimation to light, nutrients, and temperature, Limnol. Oceanogr., 43, 679–694, 1998. 

Gilpin, L.: The influence of changes in nitrogen: silicon ratios on diatom growth dynamics, J. Sea Res., 51, 21–35, 2004. 

Goldman, J. C. and Caron, D. A.: Experimental studies on an omnivorous microflagellate: implications for grazing and nutrient regeneration in the marine microbial food chain, Deep Sea Res., 32, 899–915, 1985. 

Gruber, N., Frenzel, H., Doney, S. C., Marchesiello, P., McWilliams, J. C., Moisan, J. R., Oram, J. J., Plattner, G., and Stolzenbach, K. D.: Eddy-resolving simulation of plankton ecosystem dynamics in the California Current System, Deep Sea Res., 53, 1483–1516, 2006. 

Guillard, R. L. L.: Culture of phytoplankton for feeding marine invertebrates, in: Culture of marine invertebrates animals, edited by: Smith, W. L. and Chanley, M. H., Plenum Press, New York, USA, 29–60, 1975. 

Haecky, P. and Andersson, A.: Primary and bacterial production in sea ice in the northern Baltic Sea, Aquat. Microb. Ecol., 20, 107–118, 1999. 

Harrison, W. G. and Cota, G. F.: Primary production in polar waters: relation to nutrient availability, Polar Res., 10, 87–104, 1991. 

Harrison, W. G., Head, E. J. H., Conover, R. J., Longhurst, A. R., and Sameoto, D. D.: The distribution and metabolism of urea in the eastern Canadian Arctic, Deep Sea Res., 32, 23–42, 1985. 

Hauck, J., Völker, C., Wang, T., Hoppema, M., Losch, M., and Wolf-Gladrow, D. A.: Seasonally different carbon flux changes in the Southern Ocean in response to the southern annular mode, Global Biogeochem. Cy., 27, 1–10, 2013. 

Henson, S. A., Cole, H. S., Hopkins, J., Martin, A. P., and Yool, A.: Detection of climate change-driven trends in phytoplankton phenology, Glob. Change Biol., 24, 101–111, 2018. 

Hildebrand, M.: Lack of coupling between silicon and other elemental metabolisms in diatoms, J. Phycol., 38, 841–843, 2002. 

Hohn, S.: A model of the carbon:nitrogen:silicon stoichiometry of diatoms based on metabolic processes, PhD thesis, Universität Bremen, Bremen, Germany, 135 pp., 2009. 

Hünken, M., Harder, J., and Kirst, G. O.: Epiphytic bacteria on the Antarctic ice diatom Amphiprora kufferathii Manguin cleave hydrogen peroxide produced during algal photosynthesis, Plant Biol., 10, 519–526, 2008. 

Iversen, K. R. and Seuthe, L.: Seasonal microbial processes in a high-latitude fjord (Kongsfjorden, Svalbard): I. Heterotrophic bacteria, picoplankton and nanoflagellates, Polar Biol., 34, 731–749, 2011. 

Jacobsen, T. R. and Rai, H.: Comparison of spectrophotometric, fluorometric and high performance liquid chromatography methods for determination of chlorophyll a in aquatic samples: effects of solvent and extraction procedures, Int. Rev. ges. Hydrobio., 75, 207–217, 1990. 

Johnson, M., Sanders, R., Avgoustidi, V., Lucas, M., Brown, L., Hansell, D., Moore, M., Gibb, S., Liss, P., and Jickells, T.: Ammonium accumulation during a silicate-limited diatom bloom indicates the potential for ammonia emission events, Mar. Chem., 106, 63–75, 2007. 

Jansson, M., Hickler, T., Jonsson, A., and Karlsson, J.: Links between terrestrial primary production and bacterial production and respiration in lakes in a climate gradient in subarctic Sweden, Ecosystems, 11, 367–376, 2008. 

Kamatani, A.: Dissolution rates of silica from diatoms decomposing at various temperatures, Mar. Biol., 68, 91–96, 1982. 

Keck, A. and Wassmann, P.: Temporal and spatial patterns of sedimentation in the subarctic fjord Malangen, northern Norway, Sarsia, 80, 259–276, 1996. 

Kim, S. J., Kim, B. G., Park, H. J., and Yim, J. H.: Cryoprotective properties and preliminary characterization of exopolysaccharide (P-Arcpo 15) produced by the Arctic bacterium Pseudoalteromonas elyakovii Arcpo 15, Prep. Biochem. Biotechnol., 46, 261–266, 2016. 

Kirchman, D. L.: Uptake and regeneration of inorganic nutrients by marine heterotrophic bacteria, in: Microbial ecology of the oceans, edited by: Kirchmann, D. L., John Wiley & Sons, 

Kirchman, D. L., Morán, X. A. G., and Ducklow, H.: Microbial growth in the polar oceans – role of temperature and potential impact of climate change, Nat. Rev. Microbiol., 7, 451–459, 2009. 

Kishi, M. J., Kashiwai, M., Ware, D. M., Megrey, B. A., Eslinger, D. L., Werner, F. E., Noguchi-Aita, M., Azumaya, T., Fujii, M., Hashimoto, S., Huang, D., Iizumi, H., Ishida, Y., Kang, S., Kantakov, G. A., Kim, H., Komatsu, K., Navrotsky, V. V., Smith, S. L., Tadokoro, K., Tsuda, A., Yamamura, O., Yamanaka, Y., Yokouchi, K., Yoshie, N., Zhang, J., Zuenko, Y. I., and Zvalinsky, V. I.: NEMURO – a lower trophic level model for the North Pacific marine ecosystem, Ecol. Model., 202, 12–25, 2007. 

Krause, J. W., Schulz, I. K., Rowe, K. A., Dobbins, W., Winding, M. H., Sejr, M. K., Duarte, C. M., and Agustí, S.: Silicic acid limitation drives bloom termination and potential carbon sequestration in an Arctic bloom, Sci. Rep., 9, 1–11, 2019. 

Lachmann, S. C., Mettler-Altmann, T., Wacker, A., and Spijkerman, E.: Nitrate or ammonium: Influences of nitrogen source on the physiology of a green algae, Ecol. Evol., 9, 1070–1082, 2019. 

Lannuzel, D., Tedesco, L., Van Leeuwe, M., Campbell, K., Flores, H., Delille, B., Miller, L., Stefels, J., Assmy, P., Bowman, J., Brown, K., Castellani, G., Chierici, M., Crabeck, O., Damm, E., Else, B., Fransson, A., Fripiat, F., Geilfus, N., Jacques, C., Jones, E., Kaartokallio, H., Kotovitch, M., Meiners, K., Moreau, S., Nomura, D., Peeken, I., Rintala, J., Steiner, N., Tison, J., Vancoppenolle, M., van der Linden, F., Vichi, M., and Wongpan, P.: The future of Arctic sea-ice biogeochemistry and ice-associated ecosystems, Nat. Clim. Change, 10, 1–10, 2020. 

Lee, S. and Fuhrman, J. A.: Relationships between biovolume and biomass of naturally derived marine bacterioplankton, Appl. Environ. Microbiol., 53, 1298–1303, 1987. 

Legendre, L. and Rassoulzadegan, F.: Plankton and nutrient dynamics in marine waters, Ophelia, 41, 153–172, 1995. 

Le Quéré, C., Andrew, R. M., Canadell, J. G., Sitch, S., Korsbakken, J. I., Peters, G. P., Manning, A. C., Boden, T. A., Tans, P. P., Houghton, R. A., Keeling, R. F., Alin, S., Andrews, O. D., Anthoni, P., Barbero, L., Bopp, L., Chevallier, F., Chini, L. P., Ciais, P., Currie, K., Delire, C., Doney, S. C., Friedlingstein, P., Gkritzalis, T., Harris, I., Hauck, J., Haverd, V., Hoppema, M., Klein Goldewijk, K., Jain, A. K., Kato, E., Körtzinger, A., Landschützer, P., Lefèvre, N., Lenton, A., Lienert, S., Lombardozzi, D., Melton, J. R., Metzl, N., Millero, F., Monteiro, P. M. S., Munro, D. R., Nabel, J. E. M. S., Nakaoka, S., O'Brien, K., Olsen, A., Omar, A. M., Ono, T., Pierrot, D., Poulter, B., Rödenbeck, C., Salisbury, J., Schuster, U., Schwinger, J., Séférian, R., Skjelvan, I., Stocker, B. D., Sutton, A. J., Takahashi, T., Tian, H., Tilbrook, B., van der Laan-Luijkx, I. T., van der Werf, G. R., Viovy, N., Walker, A. P., Wiltshire, A. J., and Zaehle, S.: Global Carbon Budget 2016, Earth Syst. Sci. Data, 8, 605–649,, 2016. 

Lippemeier, S., Hartig, P., and Colijn, F.: Direct impact of silicate on the photosynthetic performance of the diatom Thalassiosira weissflogii assessed by on-and off-line PAM fluorescence measurements, J. Plankton Res., 21, 269–283, 1999. 

Loebl, M., Colijn, F., van Beusekom, J. E., Baretta-Bekker, J. G., Lancelot, C., Philippart, C. J., Rousseau, V., and Wiltshire, K. H.: Recent patterns in potential phytoplankton limitation along the Northwest European continental coast, J. Sea Res., 61, 34–43, 2009. 

Ma, L. Y., Chi, Z. M., Li, J., and Wu, L. F.: Overexpression of alginate lyase of Pseudoalteromonas elyakovii in Escherichia coli, purification, and characterization of the recombinant alginate lyase, World J. Microbiol. Biotechnol., 24, 89–96, 2008. 

Martin-Jézéquel, V., Hildebrand, M., and Brzezinski, M. A.: Silicon Metabolism in Diatoms: Implications for Growth, J. Phycol., 36, 821–840, 2000. 

McCarthy, J. J. and Taylor, W. R.: and Taft, J. L., Nitrogenous nutrition of the plankton in the Chesapeake Bay, 1. Nutrient availability and phytoplankton preferences, Limnol. Oceanogr., 22, 996–1011,, 1977. 

Mills, M. M., Brown, Z. W., Laney, S. R., Ortega-Retuerta, E., Lowry, K. E., Van Dijken, G. L., and Arrigo, K. R.: Nitrogen limitation of the summer phytoplankton and heterotrophic prokaryote communities in the Chukchi Sea, Front. Mar. Sci., 5, 362,, 2018. 

Moore, C. M., Mills, M. M., Arrigo, K. R., Berman-Frank, I., Bopp,L., Boyd, P. W., Galbraith, E. D., Geider, R. J., Guieu, C., Jaccard, S. L., Jickells, T. D., La Roche, J., Lenton, T. M., Ma-howald, N. M., Maranon, E., Marinov, I., Moore, J. K., Nakat-suka, T., Oschlies, A., Saito, M. A., Thingstad, T. F., Tsuda, A., and Ulloa, O.: Processes and patterns of oceanic nutrient limitation, Nat. Geosci., 6, 701–710, 2013. 

Moore, J. K., Doney, S. C., and Lindsay, K.: Upper ocean ecosystem dynamics and iron cycling in a global three-dimensional model, Global Biogeochem. Cy., 18, GB4028,, 2004. 

Morris, I.: Nitrogen assimilation and protein synthesis, in: Algal physiology and biochemistry, edited by: Stewart, W. D. P., University of California press, 583–609, 1974. 

Mühlenbruch, M., Grossart, H. P., Eigemann, F., and Voss, M.: Mini-review: Phytoplankton-derived polysaccharides in the marine environment and their interactions with heterotrophic bacteria, Environ. Microbiol., 20, 2671–2685, 2018. 

Nelson, D. M. and Gordon, L. I.: Production and pelagic dissolution of biogenic silica in the Southern Ocean, Geochim. Cosmochim. Ac., 46, 491–501, 1982. 

Nelson, D. M., Treguer, P., Brzezinski, M. A., Leynaert, A., and Queguiner, B.: Production and dissolution of biogenic silica in the ocean: revised global estimates, comparison with regional data and relationship to biogenic sedimentation, Global Biogeochem. Cy., 9, 359–372, 1995. 

Opdal, A. F., Lindemann, C., and Aksnes, D. L.: Centennial decline in North Sea water clarity causes strong delay in phytoplankton bloom timing, Glob. Change Biol., 25, 3946–3953, 2019. 

Pahlow, M.: Linking chlorophyll-nutrient dynamics to the Redfield N: C ratio with a model of optimal phytoplankton growth, Mar. Ecol. Prog. Ser., 287, 33–43, 2005. 

Pedersen, M. F. and Borum, J.: Nutrient control of algal growth in estuarine waters. Nutrient limitation and the importance of nitrogen requirements and nitrogen storage among phytoplankton and species of macroalgae, Mar. Ecol. Prog. Ser., 142, 261–272, 1996 

Pella, E. and Colombo, B.: Study of carbon, hydrogen and nitrogen determination by combustion-gas chromatography, Microchim. Ac., 61, 697–719, 1973. 

Porter, K. G. and Feig, Y. S.: The use of DAPI for identifying and counting aquatic microflora, Limnol. Oceanogr., 25, 943–948,, 1980. 

Price, W. L.: A controlled random search procedure for global optimisation, Computer J., 20, 367–370, 1977. 

Ratkova, T. N. and Wassmann, P.: Seasonal variation and spatial distribution of phyto- and protozooplankton in the central Barents Sea, J. Mar. Syst., 38, 47–75, 2002. 

Redfield, A. C.: On the proportions of organic derivatives in sea water and their relation to the composition of plankton, James Johnstone memorial volume, 176–192, 1934. 

Rey, F. and Skjoldal, H. R.: Consumption of silicic acid below the euphotic zone by sedimenting diatom blooms in the Barents Sea, MEPS, 36, 307–312, 1987. 

Ross, O. N. and Geider, R. J.: New cell-based model of photosynthesis and photo-acclimation: accumulation and mobilisation of energy reserves in phytoplankton, Mar. Ecol. Prog. Ser., 383, 53–71, 2009. 

Saiz, E., Calbet, A., Isari, S., Anto, M., Velasco, E. M., Almeda, R., Movilla, J., and Alcaraz, M.: Zooplankton distribution and feeding in the Arctic Ocean during a Phaeocystis pouchetii bloom, Deep Sea Res., 72, 17–33, 2013. 

Schartau, M., Engel, A., Schröter, J., Thoms, S., Völker, C., and Wolf-Gladrow, D.: Modelling carbon overconsumption and the formation of extracellular particulate organic carbon, Biogeosciences, 4, 433–454,, 2007. 

Schourup-Kristensen, V., Wekerle, C., Wolf-Gladrow, D., and Völker, C.: Arctic Ocean biogeochemistry in the high resolution FESOM 1.4-REcoM2 model, Prog. Oceanogr., 168, 65–81,, 2018. 

Slagstad, D., Wassmann, P. F., and Ellingsen, I.: Physical constrains and productivity in the future Arctic Ocean, Front. Mar. Sci., 2, 1–23,, 2015. 

Smith, K. M., Kern, S., Hamlington, P. E., Zavatarelli, M., Pinardi, N., Klee, E. F., and Niemeyer, K. E.: BFM17 v1.0: Reduced-Order Biogeochemical Flux Model for Upper Ocean Biophysical Simulations, Geosci. Model Dev. Discuss. [preprint],, in review, 2020. 

Soetaert, K. and Herman, P. M. J.: A Practical Guide to Ecological Modelling – Using R as a Simulation Platform, Springer, Springer-Verlag, New York, 2009. 

Soetaert, K. and Petzoldt, T.: Inverse Modelling, Sensitivity and Monte Carlo Analysis in R Using Package FME, J. Stat. Softw., 33, 1–28,, 2010. 

Soetaert, K., Petzoldt, T., and Setzer, R. W.: Solving Differential Equations in R: Package deSolve, J. Stat. Softw., 33, 1548–7660,, 2010. 

Sommer, U., Aberle, N., Engel, A., Hansen, T., Lengfellner, K., Sandow, M., Wohlers, J., Zollner, E., and Riebesell, U.: An indoor mesocosm system to study the effect of climate change on the late winter and spring succession of Baltic Sea phyto-and zooplankton, Oecologia, 150, 655–667, 2007. 

Spilling, K., Tamminen, T., Andersen, T., and Kremp, A.: Nutrient kinetics modeled from time series of substrate depletion and growth: dissolved silicate uptake of Baltic Sea spring diatoms, Mar. Biol., 157, 427–436, 2010. 

Stow, C. A., Jolliff, J., McGillicuddy Jr., D. J., Doney, S. C., Allen, J. I., Friedrichs, M. A., Kenneth, A. R., and Wallhead, P.: Skill assessment for coupled biological/physical models of marine systems, J. Mar. Syst., 76, 4–15, 2009. 

Sturluson, M., Nielsen, T. G., and Wassmann, P.: Bacterial abundance, biomass and production during spring blooms in the northern Barents Sea, Deep Sea Res., 55, 2186–2198, 2008. 

Sverdrup, H. U.: On conditions for the vernal blooming of phytoplankton, Cons. Perm. Int. Expl. Mer., 18, 287–295, 1953. 

Teeling, H., Fuchs, B. M., Becher, D., Klockow, C., Gardebrecht, A., Bennke, C. M., Kassabgy, M., Huang, S., Mann, A. J., Waldmann, J., Weber, M., Klindworth, A., Otto, A., Lange, J., Bernhardt, J., Reinsch, C., Hecker, M., Peplies, J., Bockelmann, F. D., Callies, U., Gerdts, G., Wichels, A., Wiltshire, K. H., Glöckner, F. O., Schweder, T., and Amann, R.: Substrate-controlled succession of marine bacterioplankton populations induced by a phytoplankton bloom, Science, 336, 608–611, 2012. 

Teeling, H., Fuchs, B. M., Bennke, C. M., Krueger, K., Chafee, M., Kappelmann, L., Reintjes, G., Waldmann, J., Quast, C., Glöckner, F. O., Lucas, J., Wichels, A., Gerdts, G., Wiltshire, K. H., and Amann, R.: Recurring patterns in bacterioplankton dynamics during coastal spring algae blooms, Elife, 5, e11888,, 2016. 

Tezuka, Y.: The C: N: P ratio of phytoplankton determines the relative amounts of dissolved inorganic nitrogen and phosphorus released during aerobic decomposition, Hydrobiologia, 173, 55–62, 1989. 

Thangaraj, S., Shang, X., Sun, J., and Liu, H.: Quantitative proteomic analysis reveals novel insights into intracellular silicate stress-responsive mechanisms in the diatom Skeletonema dohrnii, Int. J. Mol. Sci., 20, 2540,, 2019. 

Thornton, D.: Diatom aggregation in the sea: mechanisms and ecological implications, Eur. J. Phycol., 37, 149–161, 2002. 

Tremblay, J. É. and Gagnon, J.: The effects of irradiance and nutrient supply on the productivity of Arctic waters: a perspective on climate change, in: Influence of climate change on the changing arctic and sub-arctic conditions, edited by: Nihoul, J. C., J. and Kostianoy, A. G., Springer, Dordrecht, The Netherlands, 73–93, 2009. 

Uitz, J., Claustre, H., Gentili, B., and Stramski, D.: Phytoplankton class-specific primary production in the world's oceans: Seasonal and interannual variability from satellite observations, Global Biogeochem. Cy., 24, GB3016,, 2010. 

Van den Meersche, K., Middelburg, J. J., Soetaert, K., Van Rijswijk, P., Boschker, H. T., and Heip, C. H.: Carbon-nitrogen coupling and algal-bacterial interactions during an experimental bloom: Modeling a 13C tracer experiment, Limnol. Oceanogr., 49, 862–878, 2004. 

Vichi, M., Pinardi, N., and Masina, S.: A generalized model of pelagic biogeochemistry for the global ocean ecosystem. Part I: Theory, J. Mar. Syst., 64, 89–109, 2007. 

Vonnahme, T. R.: Model code for: Modelling Silicate – Nitrate – Ammonium co-limitation of algal growth and the importance of bacterial remineralisation based on an experimental Arctic coastal spring bloom culture study (Version 6.0),, 2021a. 

Vonnahme, T. R.: Dynamic-Algae-Bacteria-model, available at: (last access: 29 January 2021), 2021b. 

Vonnahme, T. R., Leroy, M., Thoms, S., van Oevelen, D., Harvey, R., Kristiansen, S., Dietrich, U., Gradinger, R., and Völker, C.: Replication data for: Modelling Silicate – Nitrate – Ammonium co-limitation of algal growth and the importance of bacterial remineralisation based on an experimental Arctic coastal spring bloom culture study,, 2021. 

von Quillfeldt, C. H.: Common diatom species in Arctic spring blooms: Their distribution and abundance, Bot. Mar., 43, 499–516,, 2005.  

Wassmann, P., Slagstad, D., Riser, C. W., and Reigstad, M.: Modelling the ecosystem dynamics of the Barents Sea including the marginal ice zone: II. Carbon flux and interannual variability, J. Mar. Syst., 59, 1–24, 2006. 

Weitz, J. S., Stock, C. A., Wilhelm, S. W., Bourouiba, L., Coleman, M. L., Buchan, A., Follows, M. J., Fuhrman, J. A., Jover, L., Lennon, J. T., Middelboe, M., Sonderegger, D. L., Suttle, C. A., Taylor, B. P., Thingstad, T. F., Wilson, W., and Wommack, K. E.: A multitrophic model to quantify the effects of marine viruses on microbial food webs and ecosystem processes, ISME J., 9, 1352–1364, 2015. 

Werner, D.: Silicate metabolism, in: The biology of diatoms, edited by: Werner, D., Blackwell Scientific Publications, Berkeley, 111–149, 1977. 

Werner, D.: Regulation of metabolism by silicate in diatoms, in: Biochemistry of silicon and related problems, edited by: Bendz, G. and Lindqvist, I., Springer, Boston, USA, 149–176, 1978. 

Westberry, T. K., Behrenfeld, M. J., Siegel, D. A., and Boss, E.: Carbon-based primary productivity modeling with vertically resolved photoacclimation, Global Biogeochem. Cy., 22, 1–18, 2008. 

Wu, Y., Liu, S., Huang, Z., and Yan, W.: Parameter optimization, sensitivity, and uncertainty analysis of an ecosystem model at a forest flux tower site in the United States, J. Adv. Model. Earth Sy., 6, 405–419, 2014. 

Yool, A., Martin, A. P., Fernández, C., and Clark, D. R.: The significance of nitrification for oceanic new production, Nature, 447, 999–1002, 2007. 

Yool, A., Popova, E. E., and Anderson, T. R.: Medusa-1.0: a new intermediate complexity plankton ecosystem model for the global domain, Geosci. Model Dev., 4, 381–417,, 2011. 

Zambrano, J., Krustok, I., Nehrenheim, E., and Carlsson, B.: A simple model for algae-bacteria interaction in photo-bioreactors, Algal Res., 19, 155–161, 2016. 

Short summary
Diatoms are crucial for Arctic coastal spring blooms, and their growth is controlled by nutrients and light. At the end of the bloom, inorganic nitrogen or silicon can be limiting, but nitrogen can be regenerated by bacteria, extending the algal growth phase. Modeling these multi-nutrient dynamics and the role of bacteria is challenging yet crucial for accurate modeling. We recreated spring bloom dynamics in a cultivation experiment and developed a representative dynamic model.
Final-revised paper