- MITgcm - https://mitgcm.org -

Ecological Control of Subtropical Nutrient Concentrations

 

Ecological Control of Subtropical Nutrient Concentrations

Helen Hill [1] and Stephanie Dutkiewicz.

Figure 1: Multiple-Resource Experiment. (top) Emergent biogeographical provinces, defined by most dominant species, reminiscent of Longhurst (1995). (bottom) Biogeography of four major functional groups: (i) Diatom-analogs (red), (ii) other large phytoplankton (orange), (iii) <i>Prochlorococcus</i>-analogs (green), and (iv) other small phytoplankton (yellow-green). [2]

Figure 1: Multiple-Resource Experiment. (top) Emergent biogeographical provinces, defined by most dominant species, reminiscent of Longhurst (1995). (bottom) Biogeography of four major functional groups: (i) Diatom-analogs (red), (ii) other large phytoplankton (orange), (iii)Prochlorococcus-analogs (green), and (iv) other small phytoplankton (yellow-green).

In this article we spotlight recent work by Darwin Project [3] team members Stephanie Dutkiewicz [4], Mick Follows [5] and Jason Bragg [6], who have been examining the utility of resource control theory to interpret the relationships between organisms and resources in a global coupled physical-biogeochemistry-ecosystem model built around MITgcm [7].

The team find that in regions of low seasonality, resource competition theory (Tilman, ’77)  not only anticipates the competitive outcome amongst organisms but also provides a quantitative diagnostic of ecological control of nutrient concentrations. DFB’s sensitivity experiments clearly indicate control on the ambient nutrient by phytoplankton physiology as predicted by the theory.

Ecological Framework

Consider a simple system where many phytoplankton (Pj ) are nourished by nutrient (N ) resupplied at rate S. Phytoplankton are lost at rate m, and have physiology controlled by growth rate (μj ) and nutrient half-saturation (KNj). The time rate of change of nutrient concentration and phytoplankton can be simply related through equations 1 and 2.

eqn1 [8] (1)
eqn2 [9] (2)

1) In steady state (“resource competition theory”):

eqn3 [10] (3)
eqn4 [11] (4)

where R* is the ambient concentration of limiting nutrient with which a population of phytoplankton would come to equilibrium.

According to this model:

  • Nutrients are drawn down to a concentration in equilibrium with organisms which have the lowest R* thus preventing organisms with higher R* from being able to compete (K-selection); (Coexistence can occur if phytoplankton have the same R* in a perfectly steady environment).
  • Nutrient concentration is determined by phytoplankton physiology and the tight coupling between growth and loss
  • Phytoplankton concentrations are controlled by the balace between their losses and resource re-supply.

2) At the beginning of spring, in a highly seasonal region:

Growth rate matters most (r-selection). The organism with the fastest per capita growth under ideal resource conditions will exclude all others over several seasonal cycles.

eqn5 [12] (5)

Numerical Simulations

The team initialized the model simulations with 78 phytoplankton types (exhibiting an assortment of growth parameters) and 2 grazers (Follows et al, 2007, with minor modifications). This biogeochemical-ecosystem was then embedded in circulation state estimates, provided by the ECCO-GODAE [13] consortium (Wunsch and Heimbach, 2007).

Figure 2: Growth μ = μmaxγT γI γN is function of temperature (γT ), light (γI ), and nutrients (γN ). (top) Single Resource Case, (bottom) Multiple Resource Case. Green and yellow: K-strategy phytoplankton types with low half-saturation constants and low μmax. Red: r-strategy phytoplankton types with high μmax and high Wunsch, half-saturation constants. Green: small phytoplankton that cannot use nitrate (Prochlorococcus analogs). [14]

Figure 2: Growth μ = μmaxγT γI γN is function of temperature (γT), light (γI ), and nutrients (γN). (top) Single Resource Case, (bottom) Multiple Resource Case. Green and yellow: K-strategy phytoplankton types with low half-saturation constants and low μmax. Red: r-strategy phytoplankton types with high μmax and high half-saturation constants. Green: small phytoplankton that cannot use nitrate (Prochlorococcus analogs).

The team considered 2 experiments:

One nutrient with phytoplankton divided into two groups one group having low R* (low KN, high μ) and one having a higher R* ( high μ, high KN). (K-strategy)

N, P, Si and Fe cycling (ie multiple resources) with phytoplankton growth traits randomly assigned from within reasonable ranges (Fig. 2).

Single Resource Case

A look at behaviour in the single resource case helps illustrate the mechanisms at work in the more complex multi-resource case to follow:

Figure 3: Single Resource case. Fraction of biomass in r-strategy types relative to total biomass. [15]

Figure 3: Single Resource case. Fraction of biomass in r-strategy types relative to total biomass.

What type lives where?

Ecological Nutrient Control

Figure 4: Single Resource Case. Ratio of macro-nutrients of sensitivity exper- iment to control case. (top) KN doubled, (bottom) KN halved. Contours at .5 (sensitivity study has half the macro-nutrient) and 2 (sensitivity study has double the macro-nutrients). [16]

Figure 4: Single Resource Case. Ratio of macro-nutrients of sensitivity experiiment to control case. (top) KN doubled, (bottom) KN halved. Contours at .5 (sensitivity study has half the macro-nutrient) and 2 (sensitivity study has double the macro-nutrients).

  • Eq.3 suggests that nutrient concentrations will be affected by changes to the phytoplankton physiology.
  • In sensitivity studies, the team repeated the simulation doubling and halving KN, the nutrient saturation concentration from the control run value.
  • In most of the region dominated by K-strategy types, nutrient concentrations could be controlled by the phytoplankton physiology (Eq. 3); Phytoplankton concentrations however remained almost identical to the control run (Eq. 4).

Biogeography

  • Nutrients are ecologically controlled in relatively
    stable regimes (mixed layer depth annual range < 250m)

Figure 5: Annual range of mixed layer depth. Dashes contour indicates where macro-nutrient concentration can be predicted by resource control theory (see Fig. 4). Solid line indicates the region where the Prochlorococcus-analogs dominate in the multiple resource case. [17]

Figure 5: Annual range of mixed layer depth. Dashes contour indicates where macro-nutrient concentration can be predicted by resource control theory (see Fig. 4). Solid line indicates the region where the Prochlorococcus-analogs dominate in the multiple resource case.

Multiple Resource Case

How well does the theory hold for a more realistic simulation? The team
consider an experiment with with the cycling of N, P, Si and Fe and more strategy types of phytoplankton (e.g. additional nutrient limitation, differing light requirement, differing grazing potential)

Which types live where?

Figure 6: Multiple Resource case. Fraction of biomass in “large” (faster growing) relative to total biomass (compare to Fig. 3). Dashed line indicates the 0.5 contour. Solid line indicates the region where the Prochlorococcus-analogs dominate. [18]

Figure 6: Multiple Resource case. Fraction of biomass in “large” (faster growing) relative to total biomass (compare to Fig. 3). Dashed line indicates the 0.5 contour. Solid line indicates the region where the Prochlorococcus-analogs dominate.

Ecological Nutrient Control

Figure 7: Multiple Resource Case: Ratio of (left) Fe, (right) Inorganic nitrogen (N O3+N O2+N H4) of sensitivity experiment to control case. (top) KN doubled, (bottom) KN halved. [18]

Figure 7: Multiple Resource Case: Ratio of (left) Fe, (right) Inorganic nitrogen (N O3+N O2+N H4) of sensitivity experiment to control case. (top) KN doubled, (bottom) KN halved.

  • Sensitivity studies involved doubling and halving K for all nutrients.
  • Most limiting nutrient in the low and mid-latitudes is controlled by the phytoplankton assemblage
  • Iron in the Equatorial Pacific and inorganic nitrogen elsewhere.

Biogeography

In summary the team found that:

Want to find out more? Contact Steph… [19]

Stephanie Dutkiewicz [20]

Steph has been using MITgcm since 1997. When she’s not doing research she enjoys gardening and running.

Mick Follows [21]

Mick has also been using MITgcm since 1997.

References

Longhurst, A. Seasonal cycles of pelagic production and consumption, 1995: Progress in Oceanography,36, 77, doi:10.1016/0079-6611(95)00015-1 [22]

Tilman, D. Resource competition between planktonic algae: an experimental and theoretical approach, 1977, Ecology 58:338-348.

Wunsch and Heimbach, Practical global oceanic state estimation, 2006: Physica D: Nonlinear Phenomena, 230, 197-208, doi:10.1016/j.physd.2006.09.040 [23]

Publications

Dutkiewicz, S., M.J. Follows and J. Bragg. Modeling the coupling of ocean ecology and biogeochemistry, 2009: Global Biogeochemical Cycles, 23, GB4017, doi:10.1029/2008GB003405 [24]. PDF [25].

Follows, M.J., S. Dutkiewicz, S. Grant and S.W. Chisholm, 2007: Emergent biogeography of microbial communities in a model ocean. Science, 315, 1843-1846, doi: 10.1126/science.1138544 [26]. Read [27].