• Keine Ergebnisse gefunden

Microzooplankton stoichiometric plasticity inferred frommodelling mesocosm experiments in the Peruvian Upwellingregion

N/A
N/A
Protected

Academic year: 2022

Aktie "Microzooplankton stoichiometric plasticity inferred frommodelling mesocosm experiments in the Peruvian Upwellingregion"

Copied!
10
0
0

Wird geladen.... (Jetzt Volltext ansehen)

Volltext

(1)

Supplementary Material

Microzooplankton stoichiometric plasticity inferred from modelling mesocosm experiments in the Peruvian Upwelling

region

Marki, A.* and Pahlow, M.

* Correspondence: Marki, Alexandra: amarki@geomar.de

1 Supplementary Data

Table S1 shows parameter estimates for dinoflagellates (representing Gymnodinium sp.) and thus supplements Table 1 in the main text.

Figure S1 shows all our model configurations and supplements Figure 2 of the main text.

Figures S2 and S3 show our results of the specialist NNPZ-s and NNPZ-s-zooQP configurations.

Please consult the main text for further details.

2 Model configurations

We set up several model configurations (Figure S1), which differ in model complexity in terms of the number of trophic levels and the trophic strategies of the microzooplankton grazers. We simulate two nutrients, dissolved inorganic nitrogen and phosphorus (DIN and DIP, respectively), phytoplankton, and one to three trophic levels, representing phytoplankton, P, and two microzooplankton grazers, Z1 and Z2, to experiment with ecosystem models of different complexity (Figure S1).

We simulate "bottom-up" control with the optimality based chain model (OCM) representing primary production , and "top-down" control with the optimal current feeding model (OCF) (Pahlow and Prowe, 2010). Although the OCF was originally developed for current feeders, it can also describe other foraging strategies or feeding modes (Pahlow and Prowe, 2010). For example, ciliates (and copepods) create a feeding current (Jørgensen, 1983;Stoecker, 1984), whereas dinoflagellates are cruise feeders, which increases the foraging efficiency by increasing prey encounter rate, but also the risk of becoming a prey by increasing encounters with higher predators (Visser et al., 2008;Pahlow and Prowe, 2010). Since the microzooplankton community in the mesocosms was identified as comprising ciliate and dinoflagellate species (Hauss et al., 2012), the foraging strategy in our model is defined by the dinoflagellate and/or ciliate parameter set (Table 1, Figure 2, Table S1 and Figure S1). Furthermore, we investigate the role of specialised (strict herbivores and/or carnivores) vs. omnivorous feeding behaviour and stoichiometric plasticity of the zooplankton as possible physiological response mechanisms to changes in food quality (in terms of C:N:P ratios).

3 Model complexity 1

2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30

(2)

Supplementary Material

The simplest (NNP) configuration contains only the nutrient and phytoplankton compartments and has 6 state variables (DIN, DIP, phyto POC, phyto PON, phyto POP and Chl; see Figures 3 and 4). In the intermediate (NNPZ) configuration we introduce a second trophic level, the zooplankton guild as primary predators, by employing one additional state variable: zoo POC (Z) (Figure 2), representing the total microzooplankton biomass (BMtot) (Figure 2). Our most complex model configuration (NNPZZ) comprises three trophic levels. Here we include two microzooplankton compartments, the primary and secondary predators, with the corresponding state variables Z1, representing the dinoflagellate biomass, and Z2, representing the ciliate biomass (Figure S1). In all our model configurations we employ the same phytoplankton parameters as given in Tables 1 and S1. The microzooplankton compartment(s) in the NNPZ and NNPZZ model configurations are parameterised with the pre-calibrated parameter sets of Pahlow and Prowe (2010) employing either the parameters of the ciliate Strobilidium spiralis or the parameters of the dinoflagellate Gymnodinium spec.

(indicated by the suffix "-D-" in the configuration name; Table S1 and Figure S1).

Our three main specialist configurations (NNPZ-s, NNPZ-Ds and NNPZZ–s) represent the simplest food web structure, a linear food chain, where the organisms are treated as strict food specialists feeding on only one type of food. In the specialist configurations we assume that the only (Z) or first (Z1) predator is grazing on a strict phytoplankton diet. In the specialist NNPZZ-s configuration the second predator (Z2) is preying on the first predator (Z1) and thus is assumed to be a strict carnivore. The omnivore NNPZ-oD configuration is parameterized with the dinoflagellate parameters (supplementary Table S1 and Figure S1). Since both the specialist NNPZ-sD and the omnivore NNPZ-oD simulations produce very similar results and underestimate zooplankton growth (not shown), our model evaluation in the main text focuses mainly on the specialist NNPZ-s and omnivore NNPZ-o model configurations (Figures 2-8, Figures S1-S3).

The NNPZZ-v model configuration represents partial omnivory, where the primary predator (Z1) is solely grazing on phytoplankton, and the secondary predator (Z2) has a mixed diet consisting of his preferred food (Z1=dinoflagellates) and (less preferred) phytoplankton (Kiørboe et al., 1996) (Figure S1). The extended omnivory (NNPZZ-o) configuration is obtained by also allowing for intraguild predation within the primary and secondary predator compartments (Figure S1). In the NNPZZ-o simulation we lower 𝜙 for Z2-intraguild predation to one-fourth, because we assume that two other food sources (Z1 and phytoplankton) were available, and feeding on the own guild (Z2) would increase the risk of becoming a prey (Figure S1). Due to the observed size differences we assume that dinoflagellates (Z1) do not prey on ciliates (Z2). However, the addition of a third trophic level with different feeding strategies in our most complex configuration (NNPZZ, Figure S1) cannot explain the behaviour of the high and low DIN:DIP treatments at the same time; indeed, the NNPZZ configurations only show marginal differences compared to the NNPZ-o configurations with no significant improvement in model performance for both experiments (not shown).

3.1 Specialist configurations (NNPZ-s/NNPZ-s-zooQP):

Figures S2 and S3 show modelling results of the PU1 and PU2 mesocosm experiments, respectively, obtained by implementing a simple specialist food chain, i.e. nutrients are taken up by phytoplankton, which is grazed by microzooplankton. We use a pre-calibrated parameter-set for ciliates according to Pahlow and Prowe (2010) (Table 1). These NNPZ-s (low microzooplankton P:C quota, Table1) and NNPZ-s-zooQP (high P:C quota, Table 1) configurations clearly show the rather strong top-down control by the microzooplankton compartment, which here is the top of the food-web. To control the

2 2

31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 65 66

67 68 69 70 71 72 73 3

(3)

described and shown in the main text of the manuscript in the NNPZ-o and NNPZ-o-zooQP configurations (Figures 5 to 8). This prevents excessive grazing of phytoplankton and controls the zooplankton compartment.

74 75 76 77

(4)

Supplementary Material 4 Supplementary Figures and Tables

Table S1: Symbol definitions, units and parameter estimates of the optimality-based chain model for phytoplankton and the optimal current feeding model for zooplankton;

Symbol Units Estimates Definition

phytoplankton parameters

�0 m3 mmol−1 d−1 0.15 nutrient affinity

α mol m2 E−1 (g Chl)−1 0.9 light absorption coefficient

Q0N molN molC−1 0.07 N subsistence quota

Q0P molP molC−1 0.0019 P subsistence quota

ζ Chl molC (g Chl)−1 0.5 cost of photosynthesis

ζ N molN molC−1 0.6 cost of DIN uptake

V0 mol molC−1 5 maximum rate parameter

microzooplankton parameters for dinoflagellates, representing Gymnodinium sp.

-- 0.33 cost of assimilation coefficient

-- 0.25 cost of foraging coefficient

𝑚�� d−1 2.9 max. specific ingestion rate

𝜙 m3 mmolC−1 2.64 prey capture coefficient

N molN molC−1 0.2 N:C ratio (N quota)

P molP molC−1 0.013—0.026* P:C ratio (P quota)

d−1 0.05 specific maintenance respiration

* range of different constant microzooplankton P:C ratios, tested for identifying the most suitable microzooplankton phosphorus quota;

4 5

78

6

(5)
(6)

Supplementary Material

Figure S1: Model configurations with prey capture coefficients and the main compartments NN=Nutrients, P=Phytoplankton and Z=Zooplankton; the suffix -D indicates that the single Z compartment is represented by dinoflagellates, assumed to behave as specialists (“-sD”) or omnivores (“-oD”), respectively; the suffix “-zooQP” represents a higher zooplankton phosphorus quota (QP = 0.0195 molP molC-1), compared to the other NNPZ configurations; the NNPZZ-s configuration simulates strict food specialists; the NNPZZ-v configuration simulates partial omnivory and the NNPZZ-o configuration simulates extended omnivory. Numbers are prey capture coefficients in m-

3mmolC-1; solid lines represent prey capture coefficients of dinoflagellates or ciliates for phytoplankton and/or microzooplankton - set to 100 %, either representing the preferential food source or food of equal quality for the predator; blue dash-dotted lines represent intraguild prey capture coefficients - either set to 50 % assuming that the microzooplankton community is split into 50% intraguild prey and 50% intraguild predators or set to 25% assuming that two equal food sources are available and that the risk of becoming a prey is higher (Z2); in the omnivore NNPZZo configuration, intraguild predation occurs for Z1 and Z2 but Z1 does not prey on Z2. aDinoflagellates are represented by Gymnodinium sp. Parameters (Table S1); bCiliates are represented by Strobilidium spiralis parameters (see main text, Table 1).

6 8

9

(7)

0 1 0 2 0 3 0 4 0

0 5 0 1 0 0 1 5 0

0 1 0 2 0 3 0

0 1 0 2 0 3 0 4 0

0 5 0 1 0 0 1 5 0

0 2 4 6

0 1 0 2 0 3 0

0 2 4 6 0 2 4 6

0 0 . 5 1 1 . 5 2

5 1 0 1 5 2 0

0 0 . 0 5 0 . 1 0 . 1 5 0 . 2

0 0 . 5 1 1 . 5 2

5 1 0 1 5 2 0

0 0 . 0 5

0 . 1

0 . 1 5

0 . 2

(8)

Supplementary Material

Figure S2: Effect of different microzooplankton elemental phosphorus quotas on specialist model feeding behaviour of the PU1-NNPZ-s configuration (A:C) with lower microzooplankton P:C quota ( QZP = 0.013 molP molC-1) and the NNPZ-s-zooQP configuration (D:E) with higher microzooplankton P:C quota ( QZP = 0.0195 molP molC-1) (Table S1); Microzooplankton biomass is initialised with the total initial zooplankton biomass (BMtot) of ciliates and dinoflagellates (Figure S1). The microzooplankton compartment is initialised with the parameters of the ciliates (Strobilidium spiralis). Left y-axis: (A,D) DIN, (B,E) phytoplankton POC and (C,F) phytoplankton PON:POP ratio; right y-axis: (A,D) DIP, (B,E) (micro)zooplankton POC and (C,F) phytoplankton Chl:C ratio; units of DIN, DIP, phytoplankton POC and (micro)zooplankton POC are mmol m-3; phytoplankton PON:POP ratio is given in mol mol-1, phytoplankton Chl:C ratio is given in g mol-1 and time is given in days (d); model discontinuities are due to dilutions; observations (marks) are daily averages, where vertical bars indicate the range between the lowest and the highest measurement of all mesocosms within one treatment.

8 11

12

(9)

0 5 10 15

0 50 100 150

0 10 20 30

0 5 10 15

0 50 100 150

0 2 4 6

0 10 20 30

0 2 4 6 0 2 4 6 0 2 4 6

0 0.5 1 1.5 2

5 10 15 20

0 0.05 0.1 0.15 0.2

0 0.5 1 1.5 2

5 10 15 20

0 0.05 0.1 0.15 0.2 NNPZs

NNPZs-zooQP

A1 A2 A3 A4

B1 B2 B3 B4

C1 C 2 C 3 C 4

D1 D2 D3 D4

E1 E2 E3 E4

F1 F2 F3 F4

DINPhytoplankton POCPhytoplankton PON:POPDINPhytoplankton POCPhytoplankton PON:POP DIPZooplankton POCChl:CDIPZooplankton POCChl:C

DIN DIP POC N:P ratio C hl:C ratio POC Observations:

Model:

time (days)

HIGH HIGH LOW (ambient) LOW

A

B

C

D

E

F

Figure S3: Same as Figure S2, but showing data and ensemble model simulations for PU2 for the specialist NNPZ-s and NNPZ-s-zooQP configurations.

(10)

Supplementary Material 5 References

Hauss, H., Franz, J.M.S., and Sommer, U. (2012). Changes in N:P stoichiometry influence taxonomic composition and nutritional quality of phytoplankton in the Peruvian upwelling. Journal of Sea Research 73, 74-85. doi: 10.1016/j.seares.2012.06.010.

Jørgensen, C.B. (1983). Fluid mechanical aspects of suspension feeding. Marine Ecology Progress Series 11, 89-103. doi: 10.3354/meps011089.

Kiørboe, T., Hansen, J.L.S., Alldredge, A.L., Jackson, G.A., Passow, U., Dam, H.G., Drapeau, D.T., Waite, A., and Garcia, C.M. (1996). Sedimentation of phytoplankton during a diatom bloom: Rates and mechanisms. Journal of Marine Research 54, 1123-1148. doi: 10.1357/0022240963213754.

Pahlow, M., and Prowe, A.E.F. (2010). Model of optimal current feeding in zooplankton. Marine Ecology Progress Series 403, 129-144. doi: 10.3354/meps08466.

Stoecker, D.K. (1984). Particle production by planktonic ciliates. Limnol. Oceanogr., 29, 930--940.

Visser, A.W., Mariani, P., and Pigolotti, S. (2008). Swimming in turbulence: zooplankton fitness in terms of foraging efficiency and predation risk. Journal of Plankton Research 31, 121-133. doi:

10.1093/plankt/fbn109.

10 14

80 81 82 83 84 85 86 87 88 89 90 91 92 93 94

15

Abbildung

Table S1: Symbol definitions, units and parameter estimates of the optimality-based chain model for phytoplankton and the optimal current feeding model for zooplankton;
Figure S3: Same as Figure S2, but showing data and ensemble model simulations for PU2 for the specialist NNPZ-s and NNPZ-s-zooQP configurations.

Referenzen

ÄHNLICHE DOKUMENTE

The second part of this chapter describes the models and algorithms that are used to obtain proven optimal solutions. The column generation algorithm depends on initial variables

b The verification of the tri-lineage differentiation ability of adipogenesis, osteogenesis, as well as chondrogenesis of hBM-MSC by Oli Red O staining, Alizarin Red Staining and

In December 2011, the programme was decommissioned due to the following reasons: (i) a clash of competences between local and regional governments, (ii) that the process did not

We conclude that branching of lower order roots can be considered a leading root trait of the plant economics spectrum of temperate trees, since it relates to the

Hyphens indicate no estimation due to too small a number of observed events in a subgroup to develop a reliable regression model.. CVD, cardiovascular disease;

(a-d) MiR-100 mice showed no altered glucose tolerance in the GTT (a+b) and insulin sensitivity in the ITT (c+d) compared to wildtype animals fed normal chow diet before the

The model has been proposed for and used with K-12 to university teachers with varying teaching experiences (pre-service, in- service) and is also not dedicated to any

Kaplan-Meier estimates of PFS according to RECIST 1.1 and IRR evaluation (A) and OS (B) among subjects in the ITT population with CD24 - /CD44 + CSC in metastatic tissue..