Biogeosciences www.biogeosciences.net 1726-4170 1726-4189 3 4 2006 10.5194/bg-3-585-2006 http://www.biogeosciences.net/3/585/2006/ http://www.biogeosciences.net/3/585/2006/bg-3-585-2006.html http://www.biogeosciences.net/3/585/2006/bg-3-585-2006.pdf 585 606 2006-11-28 Multi-nutrient, multi-group model of present and future oceanic phytoplankton communities E. Litchman litchman@msu.edu C. A. Klausmeier J. R. Miller O. M. Schofield P. G. Falkowski Institute of Marine and Coastal Sciences, Rutgers University, New Brunswick, NJ 08901, USA Michigan State University, Kellogg Biological Station, MI 49060, USA Department of Ecology and Evolutionary Biology, Princeton University, Princeton, NJ 08544, USA Phytoplankton community composition profoundly affects patterns of nutrient cycling and the dynamics of marine food webs; therefore predicting present and future phytoplankton community structure is crucial to understand how ocean ecosystems respond to physical forcing and nutrient limitations. We develop a mechanistic model of phytoplankton communities that includes multiple taxonomic groups (diatoms, coccolithophores and prasinophytes), nutrients (nitrate, ammonium, phosphate, silicate and iron), light, and a generalist zooplankton grazer. Each taxonomic group was parameterized based on an extensive literature survey. We test the model at two contrasting sites in the modern ocean, the North Atlantic (North Atlantic Bloom Experiment, NABE) and subarctic North Pacific (ocean station Papa, OSP). The model successfully predicts general patterns of community composition and succession at both sites: In the North Atlantic, the model predicts a spring diatom bloom, followed by coccolithophore and prasinophyte blooms later in the season. In the North Pacific, the model reproduces the low chlorophyll community dominated by prasinophytes and coccolithophores, with low total biomass variability and high nutrient concentrations throughout the year. Sensitivity analysis revealed that the identity of the most sensitive parameters and the range of acceptable parameters differed between the two sites. We then use the model to predict community reorganization under different global change scenarios: a later onset and extended duration of stratification, with shallower mixed layer depths due to increased greenhouse gas concentrations; increase in deep water nitrogen; decrease in deep water phosphorus and increase or decrease in iron concentration. To estimate uncertainty in our predictions, we used a Monte Carlo sampling of the parameter space where future scenarios were run using parameter combinations that produced acceptable modern day outcomes and the robustness of the predictions was determined. Change in the onset and duration of stratification altered the timing and the magnitude of the spring diatom bloom in the North Atlantic and increased total phytoplankton and zooplankton biomass in the North Pacific. Changes in nutrient concentrations in some cases changed dominance patterns of major groups, as well as total chlorophyll and zooplankton biomass. Based on these scenarios, our model suggests that global environmental change will inevitably alter phytoplankton community structure and potentially impact global biogeochemical cycles. Ammerman, J. W., Hood, R. R., Case, D. A., and Cotner, J. B.: Phosphorus deficiency in the Atlantic: an emerging paradigm in oceanography, EOS, 84, 165, 2003. Balch, W. M., Holligan, P. M., and Kilpatrick, K. A.: Calcification, photosynthesis and growth of the bloom-forming coccolithophore, \textitEmiliania huxleyi, Continental Shelf Res., 12, 1353–1374, 1992. Barlow, R. G., Mantoura, R. F. C., Peinert, R. D., Miller, A. E. J., and Fileman, T. W.: Distribution, sedimentation and fate of pigment biomarkers following thermal stratification in the Western Alboran Sea, Marine Ecology Progress Series, 125, 279–291, 1995. Bissett, W. P., Walsh, J. J., Dieterle, D. A., and Carder, K. L.: Carbon cycling in the upper waters of the Sargasso Sea. I. Numerical simulation of differential carbon and nitrogen fluxes, Deep-Sea Res. I, 46, 205–269, 1999. Brand, L. E. and Guillard, R. R. L.: The effects of continuous light and light intensity on the reproduction rates of twenty-two species of marine phytoplankton, J. Exp. Mar. Biol. Ecol., 50, 119–132, 1981. Berlow, E. L.: Strong effects of weak interactions in ecological communities, Nature, 398, 330–334, 1999. Béthoux, J. P., Morin, P., and Ruiz-Pino, D. P.: Temporal trends in nutrient ratios: chemical evidence of Mediterranean ecosystem changes driven by human activity, Deep-Sea Res. II, 49, 2007–2016, 2002. Boyd, P. W. and Harrison, P. J.: Phytoplankton dynamics in the NE subarctic Pacific, Deep-Sea Res. II, 46, 2405–2432, 1999. Boyd, P. W., Sherry, N. D., Berges, J. A., Bishop, J. K. B., Calvert, S. E., Charette, M. A., Giovannoni, S. J., Goldblatt, R., Harrison, P. J., Moran, S. B., Roy, S., Soon, M., Strom, S., Thibault, D., Vergin, K. L., Whitney, F. A., and Wong, C. S.: Transformations of biogenic particulates from the pelagic to the deep ocean realm, Deep-Sea Res. II, 46, 2761–2792, 1999. Boyd, P. W., Watson, A. J., Law, C. S., et al.: A mesoscale phytoplankton bloom in the polar Southern Ocean stimulated by iron fertilization, Nature, 407, 695–702, 2000. Bowie, A. R., Whitworth, D. J., Achterberg, E. P., Mantoura, R. F. C., and Worsfold, P. J.: Biogeochemistry of Fe and other trace elements (Al, Co, Ni) in the upper Atlantic Ocean, Deep-Sea Res. I, 49, 605–636, 2002. Brock, T. D.: Calculating solar radiation for ecological studies, Ecol. Modell., 14, 1–19, 1981. Broerse, A. T. C., Ziveri, P., van Hinte, J. E., and Honjo, S.: Coccolithophore export production, species composition, and coccolith-CaCO₃ fluxes in the NE Atlantic (34 degrees N 21 degrees W and 48 degrees N 21 degrees W), Deep Sea Res. II, 47, 1877–1905, 2000. Bury, S. J., Boyd, P. W., Preston, T., Savidge, G., and Owens, N. J. P.: Size-fractionated primary production and nitrogen uptake during a North Atlantic phytoplankton bloom: implications for carbon export estimates, Deep-Sea Res. I, 48, 689–720, 2001. Chavez, F. P., Ryan, J., Lluch-Cota, S. E., and Niquen, C. M.: From anchovies to sardines and back: multidecadal change in the Pacific Ocean, Science, 299, 217–221, 2003. Denman, K. L. and Peña, M. A.: A coupled 1-D biological/physical model of the northeast subarctic Pacific Ocean with iron limitation, Deep-Sea Res. II, 46, 2877–2908, 1999. Doney, S. C., Wallace, D. W. R., and Ducklow, H. W.: The North Atlantic carbon cycle: new perspectives from JGOFS and WOCE, in: The changing ocean carbon cycle, edited by: Hanson, R. B., Ducklow, H. W., and Field, J. G., 375–391, 2000. Doney, S. C., Kleypas, J. A., Sarmiento, J. L., and Falkowski, P. G.: The US JGOFS Synthesis and Modeling Project – an introduction, Deep-Sea Res. II, 49, 1–20, 2002. Dortch, Q.: The interaction between ammonium and nitrate uptake in phytoplankton, Mar. Ecol. Prog. Series, 61, 183–201, 1990. Droop, M. R.: Some thoughts on nutrient limitation in algae, J. Phycol., 9, 264–272, 1973. Dugdale, R. and Wilkerson, F.: Silicate regulation of new production in the equatorial Pacific upwelling, Nature, 391, 270–273, 1998. Egge, J. K.: Are diatoms poor competitors at low phosphate concentrations?, J. Mar. Syst., 16, 191–198, 1998. Engel, A., Zondervan, I., Aerts, K., et al.: Testing the direct effect of CO₂ concentration on a bloom of the coccolithophorid \textitEmiliania huxleyi in mesocosm experiments, Limnol. Oceanogr., 50, 493–507, 2005. Evans, G. T.: A framework for discussing seasonal succession and coexistence of phytoplankton species, Limnol. Oceanogr., 33, 1027–1036, 1988. Evans, G. T. and Parslow, J. S.: A model of annual plankton cycles, Biol. Oceanogr., 3, 327–347, 1985. Falkowski, P. G., Laws, E. A., Barber, R. T., and Murray, J. W.: Phytoplankton and their role in primary, new, and export production, in: Ocean biogeochemistry: the role of the ocean carbon cycle in global change, edited by: Fasham, M. J. R., Springer, 99–121, 2003. Falkowski, P. G., Katz, M. E., Knoll, A. H., Quigg, A., Raven, J. A., Schofield, O., and Taylor, F. J. R.: The evolution of modern eukaryotic phytoplankton, Science, 305, 354–360, 2004. Fasham, M. J. R., Ducklow, H. W., and McKelvie, S. M.: A nitrogen-based model of plankton dynamics in the oceanic mixed layer, J. Mar. Res., 48, 591–639, 1990. Fasham, M. J. R. and Evans, G. T.: Advances in ecosystem modeling within JGOFS, in: The changing ocean carbon cycle, edited by: Hanson, R. B., Ducklow, H. W., and Field, J. G., 417–446, 2000. Fennel, K., Abbott, M. R., Spitz, Y. H., Richman, J. G., and Nelson, D. M.: Impacts of iron control on phytoplankton production in the modern and glacial Southern Ocean, Deep-Sea Res. II, 50, 833–851, 2003. Field, C. B., Behrenfeld, M. J., Randerson, J. T., and Falkowski, P. G.: Primary production of the biosphere: integrating terrestrial and oceanic components, Science, 281, 237–240, 1998. Freeland, H., Denman, K., Wong, C. S., Whitney, F., and Jacques, R.: Evidence of change in the winter mixed layer in the Northeast Pacific Ocean, Deep-Sea Res. I, 44, 2117–2129, 1997. Frost, B. W. and Kishi, M. J.: Ecosystem dynamics in the eastern and western gyres of the Subarctic Pacific – a review of lower trophic level modelling, Prog. Oceanogr., 43, 317–333, 1999. Fulton, E. A., Parslow, J. S., Smith, A. D. M., and Johnson, C. R.: Biogeochemical marine ecosystem models II: the effect of physiological detail on model performance, Ecological Modelling, 173, 371–406, 2004. Gaul, W. and Antia, A. N.: Taxon-specific growth and selective microzooplankton grazing of phytoplankton in the Northeast Atlantic, J. Mar. Syst., 30, 241–261, 2001. Gregg, W. W., Ginoux, P., Schopf, P. S., and Casey, N. W.: Phytoplankton and iron: validation of a global three-dimensional ocean biogeochemical model, Deep-Sea Res. II, 50, 3143–3169, 2003. Grover, J. P.: Dynamics of competition among microalgae in variable environments: experimental tests of alternative models, Oikos, 62, 231–243, 1991. Gruber, N., Keller, K., and Key, R. M.: What story is told by oceanic tracer concentrations?, Science, 290, 455–456, 2000. Guillard, R. R. L., Kilham, P., and Jackson, T. A.: Kinetics of silicon-limited growth in the marine diatom \textitThalassiosira pseudonana Hassle and Heimdal (\textitCyclotella nana Hustedt), J. Phycol., 9, 233–237, 1973. Hanson, R. B., Ducklow, H. W., and Field, J. G. (Eds.): The changing ocean carbon cycle: a midterm synthesis of the Joint Global Ocean Flux Study, Cambridge University Press, Cambridge, 2000. Harrison, P. J., Boyd, P. W., Varela, D. E., Takeda, S., Shiomoto, A., and Odate, T.: Comparison of the factors controlling phytoplankton productivity in the NE and NW subarctic Pacific gyres, Prog. Oceanogr., 43, 205–234, 1999. Harrison, P. J., Whitney, F. A., Tsuda, A., Saito, H., and Tadokoro, K.: Nutrient and plankton dynamics in the NE and NW gyres of the subarctic Pacific Ocean, J. Oceanogr., 60, 93–117, 2004. Holligan, P. M. and Robertson, J. E.: Significance of carbonate formation for the ocean carbon cycle, Global Change Biol., 2, 85–95, 1996. Huisman, J. and Weissing, F. J.: Light-limited growth and competition for light in well-mixed aquatic environments: an elementary model, Ecology, 75, 507–520, 1994. Hutchins, D. A. and Bruland, K. W.: Iron-limited diatom growth and Si:N uptake ratios in a coastal upwelling regime, Nature, 393, 561–564, 1998. Iglesias-Rodrígues, M. D., Brown, C. W., Doney, S. C., Kleypas, J., Kolber, D., Kolber, Z., Hayes, P. K., and Falkowski, P. G.: Representing key phytoplankton functional groups in ocean carbon cycle models: coccolithophorids, Global Biogeochem. Cycles, 16, 47(1)–47(20), 2002. Iida, T., Saitoh, S. I., Miyamura, T., Toratani, M., Fukushima, H., and Shiga, N.: Temporal and spatial variability of coccolithophore blooms in the eastern Bering Sea, 1998–2001, Prog. Oceanogr., 55, 165–175, 2002. Johnson, K. S., Gordon, R. M., and Coale, K. H.: What controls dissolved iron concentrations in the world ocean?, Mar. Chem., 57, 137–161, 1997. Joint, I. R., Pomroy, A., Savidge, G., and Boyd, P.: Size fractionated primary productivity in the North East Atlantic in Spring 1989, Deep-Sea Res. II, 40, 423–440, 1993. Karl, D. M., Bidigare, R. R., and Letelier, R. M.: Long-term changes in plankton community structure and productivity in the North Pacific Subtropical Gyre: the domain shift hypothesis, Deep-Sea Res. II, 48, 1449–1470, 2001. Kirk, J. T. O.: Light and photosynthesis in aquatic ecosystems, Cambridge University Press, Cambridge, 509 pp., 1994. Kleypas, J. A. and Doney, S. C.: Nutrients, Chlorophyll, Primary Production and Related Biogeochemical Properties in the Ocean Mixed Layer – A Compilation of Data Collected at Nine JGOFS Sites, NCAR Technical Note NCAR/TN-447+STR, 55 pp., 2001. Klausmeier, C. A., Litchman, E., and Levin, S. A.: Phytoplankton growth and stoichiometry under multiple nutrient limitation, Limnol. Oceanogr., 49, 1463–1470, 2004. Lam, P. J., Tortell, P. D., and Morel, F. M. M.: Differential effects of iron additions on organic and inorganic carbon production by phytoplankton, Limnol. Oceanogr., 46, 1199–1202, 2001. Langdon, C.: On the causes of interspecific differences in the growth-irradiance relationship for phytoplankton. II. A general review, J. Plankton Res., 10, 1291–1312, 1988. Lochte, K., Ducklow, H. W., Fasham, M. J. R., and Stienen, C.: Plankton succession and carbon cycling at 47 N 20 W during the JGOFS North Atlantic Bloom experiment, Deep-Sea Res. II, 40, 91–114, 1993. Longhurst, A. R.: Ecological geography of the sea, Academic Press, New York, 398 p., 1998. Louanchi, F. and Najjar, R. G.: A global monthly climatology of phosphate, nitrate, and silicate in the upper ocean: spring-summer export production and shallow remineralization, Global Biogeochem. Cycles, 14, 957–977, 2000. Maldonado, M. T., Boyd, P. W., Strzepek, R., LaRoche, J. L., Waite, A., Croot, P. L., and Price, N. M.: Phytoplankton physiological responses to changing iron chemistry during a mesoscale Southern Ocean iron-enrichment, Limnol. Oceanogr., 46, 1802–1808, 2001. Manabe, S., Stouffer, R. J., Spelman, M. J., and Bryan, K.: Transient responses of a coupled ocean-atmosphere model to gradual changes of atmospheric CO₂. Part I: Annual mean response, J. Climate, 4, 785–818, 1991. Marty, J.-C., Chiavérini, J., Pizay, M.-D., and Avril, B.: Seasonal and interannual dynamics of nutrients and phytoplankton pigments in the western Mediterranean Sea at the DYFAMED time-series station (1991–1999), Deep-Sea Res. II, 49, 1965–1985, 2002. Miller, J. R. and Russell, G. L.: Investigating the interactions among river flow, salinity, and sea ice using a global coupled atmosphere-ocean-ice model, Ann. Glaciol., 25, 121–126, 1997. Montagnes, D. J. S., Berges, J. A., Harrison, P. J., and Taylor, F. J. R.: Estimating carbon, nitrogen, protein, and chlorophyll a from volume in marine phytoplankton, Limnol. Oceanogr., 39, 1044–1060, 1994. Montagnes, D. J. S. and Lessard, E. J.: Population dynamics of the marine planktonic ciliate Strombidinopsis multiauris: its potential to control phytoplankton blooms, Aquat. Microb. Ecol., 20, 167–181, 1999. Moore, J. K., Doney, S. C., Kleypas, J. A., Glover, D. M., and Fung, I. Y.: An intermediate complexity marine ecosystem model for the global domain, Deep-Sea Res. II, 49, 403–462, 2002. Moore, J. K., Doney, S. C., and Lindsay, K.: Upper ocean ecosystem dynamics and iron cycling in a global three-dimensional model, Global Biogeochem. Cycles, 18, GB4028, 2004. Muggli, D. L. and Harrison, P. J.: Effects of nitrogen source on the physiology and metal nutrition of \textitEmiliania huxleyi grown under different iron and light conditions, Mar. Ecol. Prog. Ser., 130, 255–267, 1996. Napp, J. M. and Hunt, G. L.: Anomalous conditions in the south-eastern Bering Sea 1997: linkages among climate, weather, ocean, and biology, Fisheries Oceanogr., 10, 61–68, 2001. Orr, J. C., Fabry, V. J., Aumont, O., et al.: Anthropogenic ocean acidification over the twenty-first century and its impact on calcifying organisms, Nature, 437, 681–686, 2005. Pahlow, M. and Riebesell, U.: Temporal trends in deep ocean Redfield ratios, Science, 287, 831–833, 2000. Peña, M. A.: Modelling the response of the planktonic food web to iron fertilization and warming in the NE subarctic Pacific, Prog. Oceanogr., 57, 453–479, 2003. Platt, T., Fuentes-Yaco, C., and Frank, K. T.: Spring algal bloom and larval fish survival, Nature, 423, 398–399, 2003. Press, W. H., Teukolsky, S. A., Vetterling, W. T., and Flannery, B. P.: Numerical recipes in FORTRAN, 2nd ed., Cambridge University Press, 1992. Quigg, A., Finkel, Z. V., Irwin, A. J., Rosenthal, Y., Ho, T.-Y., Reinfelder, J. R., Schofield, O., Morel, F. M. M., and Falkowski, P. G.: The evolutionary inheritance of elemental stoichiometry in marine phytoplankton, Nature, 425, 291–294, 2003. Rappé, M. S., Suzuki, M. T., Vergin, K. L., and Giovannoni, S. J.: Phylogenetic diversity of ultraplankton plastid small-subunit rRNA genes recovered in environmental nucleic acid samples from the Pacific and Atlantic coasts of the United States, Appl. Environ. Microbiol., 64, 294–303, 1998. Richardson, K., Beardall, J., and Raven, J. A.: Adaptation of unicellular algae to irradiance: an analysis of strategies, New Phytologist, 93, 157–191, 1983. Ridgwell, A. J.: Dust in the Earth system: the biogeochemical linking of land, air and sea, Phil. Trans. R. Soc. Lond. A, 360, 2905–2924, 2002. Riegman, R., Stolte, W., Noordeloos, A. A. M., and Slezak, D.: Nutrient uptake and alkaline phosphatase (EC 3:1:3:1) activity of \textitEmiliania huxleyi (Prymnesiophyceae) during growth under N and P limitation in continuous cultures, J. Phycol., 36, 87–96, 2000. Russell, G. L., Miller, J. R., and Rind, D.: A coupled atmosphere-ocean model for transient climate change studies, Atmos.-Ocean, 33, 683–730, 1995. Sarmiento, J. L., Hughes, T. M. C., Stouffer, R. J., and Manabe, S.: Simulated response of the ocean carbon cycle to anthropogenic climate warming, Nature, 393, 245–249, 1998. Savidge, G., Boyd, P., Pomroy, A., Harbour, D., and Joint, I.: Phytoplankton production and biomass estimates in the northeast Atlantic Ocean, May–June 1990, Deep-Sea Res. I, 42, 599–617, 1995. Smayda, T. J.: Harmful algal blooms: Their ecophysiology and general relevance to phytoplankton blooms in the sea, Limnol. Oceanogr., 42, 1137–1153, 1997. Smetacek, V.: Diatoms and the ocean carbon cycle, Protist, 150, 25–32, 1999. Staal, M., Meysman, F. J. R., and Stal, L. J.: Temperature excludes N₂-fixing heterocystous cyanobacteria in the tropical oceans, Nature, 425, 504–507, 2003. Sunda, W. G. and Huntsman, S. A.: Iron uptake and growth limitation in oceanic and coastal phytoplankton, Mar. Chem., 50, 189–206, 1995. Sverdrup, H. U.: On conditions for the vernal blooming of phytoplankton, J. Cons. Perm. Int. Explor. Mer, 18, 287–295, 1953. Tilman, G. D.: Resource competition and community structure, Princeton University Press, Princeton, NJ, 1982. Tyrrell, T., Holligan, P., and Mobley, C.: Optical impacts of oceanic coccolithophore blooms, J. Geophys. Res., 104, 3223–3242, 1999. Tsuda, A., Takeda, S., Saito, H., et al.: A mesoscale iron enrichment in the Western Subarctic Pacific induces a large centric diatom bloom, Science, 300, 958–961, 2003. Van der Staay, S. Y. M., De Wachter, R., and Vaulot, D.: Oceanic 18S rDNA sequences from picoplankton reveal unsuspected eukaryotic diversity, Nature, 409, 607–610, 2001. Varela, D. E. and Harrison, P. J.: Seasonal variability in nitrogenous nutrition of phytoplankton assemblages in the northeastern subarctic Pacific Ocean, Deep-Sea Res. II, 46, 2505–2538, 1999. Weeks, A., Conte, M. H., Harris, R. P., et al.: The physical and chemical environment and changes in community structure associated with bloom evolution – the joint global flux study north-Atlantic bloom experiment, Deep-Sea Res. II, 40, 347–368, 1993. Woody, C., Shih, E., Miller, J., Royer, T., Atkinson, L., and Moody, R.: Salinity measurements in the coastal ocean, Workshop review, September 1998, CCPO Technical Report, 99–01, 1999. Wroblewski, J. S.: A model of phytoplankton bloom formation during variable Oregon upwelling, J. Mar. Res., 35, 357–394, 1977.