MARINE ECOLOGY PROGRESS SERIES Vol. 228: 103–117, 2002 Published March 6 Mar Ecol Prog Ser

Photoacclimation and nutrient-based model of light-saturated photosynthesis for quantifying oceanic primary production

Michael J. Behrenfeld1,*, Emilio Marañón2, David A. Siegel3, Stanford B. Hooker1

1National Aeronautics and Space Administration, Goddard Space Flight Center, Code 971, Building 33, Greenbelt, Maryland 20771, USA 2Departamento de Ecología y Biología Animal, Universidad de Vigo, 36200 Vigo, Spain 3Institute for Computational Earth System Science, University of California Santa Barbara, Santa Barbara, California 93106-3060, USA

ABSTRACT: Availability of remotely sensed biomass fields has greatly advanced pri- mary production modeling efforts. However, conversion of near-surface chlorophyll concentrations to carbon fixation rates has been hindered by uncertainties in modeling light-saturated photosynthesis b b (P max). Here, we introduce a physiologically-based model for P max that focuses on the effects of photoacclimation and nutrient limitation on relative changes in cellular chlorophyll and CO2 fixation b capacities. This ‘PhotoAcc’ model describes P max as a function of light level at the bottom of the mixed layer or at the depth of interest below the mixed layer. Nutrient status is assessed from the relationship between mixed layer and nutricline depths. Temperature is assumed to have no direct b influence on P max above 5°C. The PhotoAcc model was parameterized using photosynthesis- irradiance observations made from extended transects across the Atlantic Ocean. Model perfor- mance was validated independently using time-series observations from the Sargasso Sea. The PhotoAcc model accounted for 70 to 80% of the variance in light-saturated photosynthesis. Previously described temperature-dependent models did not account for a significant fraction of the variance in b P max for our test data sets.

KEY WORDS: Photosynthesis · Modeling · Primary production

Resale or republication not permitted without written consent of the publisher

INTRODUCTION sis has remained problematic. Whereas terrestrial pro- ductivity models suffer from a lack of observational Photosynthesis is a fundamental process of nearly all data for parameterization and testing (Field et al. known ecosystems, such that the level of photoau- 1998), high-sensitivity measurements of net primary totrophic carbon fixation supported by a given envi- production in aquatic systems have been routine since ronment broadly dictates the local biomass of subse- the introduction of the 14C method by Steemann quent trophic levels and the biogeochemical exchange Nielsen (1952). of elements between systems. Models of biospheric Analyses of vertical profiles of phytoplankton photo- primary production have been greatly aided by global- synthesis revealed early on that, when normalized scale satellite observations (Field et al. 1998), but con- to depth-specific chlorophyll concentrations, primary version of measured plant biomass to net photosynthe- production can be modeled to first order simply as a function of subsurface irradiance (Ryther 1956, Ryther & Yentsch 1957, Talling 1957). A variety of analytical *E-mail: [email protected] expressions have consequently been developed de-

© Inter-Research 2002 · www.int-res.com 104 Mar Ecol Prog Ser 228: 103–117, 2002

scribing this relationship between vertical light attenu- comparable at this point-source scale of field measure- ation and chlorophyll-normalized carbon fixation ments, but a successful model of this genre has not yet (reviewed by Platt & Sathyendranath 1993, Behrenfeld been described. & Falkowski 1997b). Such models can account for most Here we introduce a model, belonging to this sec- b of the observed variance in depth-integrated photo- ond category, that captures P opt variability in natural synthesis (∑PP), particularly when measurements phytoplankton assemblages. The model largely fo- b encompass a wide phytoplankton biomass range, pro- cuses on changes in P opt resulting from photoaccli- vided model input includes measured values for: mation and is thus referred to as the ‘PhotoAcc (1) the vertical distribution of chlorophyll, (2) the Model’, although a nutrient-dependence is also pre- downwelling attenuation coefficient (Kd) for photosyn- scribed. In the following section, we describe the- thetically active radiation (PAR), and (3) the maximum conceptual basis and underlying equations of the b carbon fixation rate per unit of chlorophyll (P opt) PhotoAcc model, the field data used for model para- (Behrenfeld & Falkowski 1997a,b). meterization and testing, and our approach to assess- For over 40 yr, developments in phytoplankton pri- ing nutrient status and photoacclimation irradiances. mary production models have focused on refining char- Model limitations, directions for expansion, and po- acterizations of the above 3 critical water-column fea- tential avenues for global implementation are ad- tures, with clearly the greatest achievements realized dressed in the section ‘Discussion’. in the description of the underwater light field (e.g. Platt & Sathyendranath 1988, Morel 1991, Antoine et al. 1996). Progress has also been made in predicting verti- METHODS cal profiles of chlorophyll (Platt & Sathyendranath 1988, b Morel & Berthon 1989), but models of P opt have re- The PhotoAcc model is intended to provide esti- b b mained rudimentary and inconsistent (Behrenfeld & mates of both P opt and P max that are robust to varia- b Falkowski 1997b). The importance of accurate P opt tions in environmental conditions and taxonomic b estimates cannot be overstated, especially when model composition. P opt is measured under ambient light performance is evaluated by comparison with point- conditions during ~6 to 24 h in situ or simulated in source field observations. For globally representative situ incubations, and is defined as the maximum data sets, phytoplankton biomass alone accounts for chlorophyll-normalized carbon fixation in a water col- <40% of ∑PP variability, while inclusion of measured umn (the superscript ‘b’ denotes normalization to b P opt values can account for >80% of the variance in chlorophyll) divided by the incubation light period ∑PP (Balch & Byrne 1994, Behrenfeld & Falkowski (Wright 1959, Behrenfeld & Falkowski 1997a). Nat- 1997a). In oligotrophic regions where the range in ural fluctuations in sunlight cause photosynthesis to chlorophyll concentration is further constrained, accu- vary during such incubations, particularly in light- b b rate estimates of P opt are even more critical (Banse & limited samples deep within the water column. P opt, Yong 1990, Siegel et al. 2000). however, is generally observed near the surface b The function of P opt models is to capture spatial and where light levels are sufficient to maintain photo- b temporal changes in assimilation efficiencies (i.e. car- synthesis near the light-saturated rate, P max. Thus, bon fixed per unit of chlorophyll) resulting from phy- unless photoinhibition is excessive or ambient PAR is b siological acclimation to environmental variability. very low, modeling P opt is synonymous with model- b Currently, the 2 principal approaches for estimating ing near-surface P max. b P opt in regional- to global-scale models are: (1) to The longer-term incubations described above treat assign fixed, climatological values to biogeographical the water column as a compound photosynthetic unit b provinces (Longhurst 1995, Longhurst et al. 1995), and and yield a single P opt value for each vertical profile. b (2) to define predictive relationships between P opt and In contrast, short-term (~0.5 to 2 h) photosynthesis- 1 or more environmental variables (e.g. temperature, irradiance (PI) measurements provide unique relation- nutrient concentration) (Megard 1972, Balch et al. ships between light and carbon fixation for each popu- 1992, Antoine et al. 1996, Behrenfeld & Falkowski lation sampled within a water column. This technique 1997a). Both techniques have advantages, but differ involves exposing phytoplankton from a given depth in their intended application. For example, the first to a range of constant light intensities. PI measure- approach furnishes broad-scale average values for ments provide critical information on depth-dependent b designated regions and is not intended to accurately changes in P max. reproduce the much finer scale physiological variabil- PhotoAcc model. The light-saturated rate of photo- ity in 14C-uptake rates corresponding to a given day, synthesis is constrained by the capacity of the Calvin depth and location. As for the second approach, an cycle reactions (Pmax), which fix CO2 into carbohy- b effective model should, ideally, provide P opt estimates drates (Stitt 1986, Sukenik et al. 1987, Orellana & Perry Behrenfeld et al.: Modeling phytoplankton photosynthesis 105

1992). If chlorophyll (chl) concentration is assumed b proportional to light absorption, then P max is the ratio of the Calvin cycle capacity to light-harvesting (i.e. b –1 P max = Pmax × chl ). Light absorption and electron transport by the photosynthetic light reactions function primarily to supply ATP and NADPH to the Calvin cycle. Thus, parallel changes in chlorophyll and the b Calvin cycle can occur without a change in P max. b Variations in P max result when acclimation to a new growth condition requires a shift in the balance between light harvesting and carbon fixation. Thus, b the essential requirement for an effective P max model is a description of how specific environmental factors cause Pmax and chlorophyll to change relative to each other. b The PhotoAcc model estimates P max by describing separate, light-dependent relationships for Pmax and chlorophyll for 3 environmental conditions: (1) nutri- ent-sufficient growth in the mixed layer, (2) nutrient- sufficient growth below the mixed layer, and (3) nutri- ent-depleted growth above the nutricline. The chal- Fig. 1. Compilation of published studies reporting relation- ships between chlorophyll concentration and growth irradi- lenging aspect of this approach is that the underlying ance (Ig). Chlorophyll concentrations were divided by a scalar relationships are not resolved in the field and are to yield unitless relative values. Scalar values varied primarily rarely measured in the laboratory. In the field, Pmax and as a function of cell size, but were also influenced by factors chlorophyll are typically measured per unit volume of such as growth temperature, taxonomic division, spectral characteristics of growth irradiance, and units in which water, so changes resulting from physiological accli- chlorophyll concentrations were expressed. These data were mation cannot be distinguished from changes due to compiled to illustrate the general relationship between variations in cell abundance. However, since variabil- chlorophyll and Ig, but do not represent an exhaustive review b of all published results. Inset shows the data without connect- ity in P max is dependent on relative changes in Pmax and chlorophyll, parameterization of the PhotoAcc ing lines, the full scale for Ig, and the fit (line) of Eq. (1). Irra- diance is expressed in units of mol quanta m–2 h–1 (3.6 mol model can be achieved using a qualitative light- quanta m–2 h–1 = 1000 µmol quanta m–2 s–1). Data sources and chlorophyll relationship from laboratory photoacclima- species tested were as follows (outliers excluded from figure b tion studies and field P max data. are identified in parentheses by incubation temperature [°C], –2 –1 Laboratory photoacclimation experiments suggest photoperiod [h], and growth irradiance [mol quanta m h ]) — Myers (1946): Chlorella pyrenoidosa; Eppley & Sloan that phytoplankton exhibit a reasonably conserved (1966): Dunaliella tertiolecta (16 h, 0.88 mol quanta m–2 h–1); response to changes in growth irradiance (Ig). To illus- Paasche (1967): Coccolithus huxleyi; Paasche (1968): Ditylum trate, we assembled chlorophyll-irradiance relation- brightwellii (10 h, 0.88 mol quanta m–2 h–1), Nitzschia ships for 23 phytoplankton species from 23 studies turgidula; Beale & Appleman (1971): Chlorella vulgaris; Durbin (1974): Thalassiosira nordenskioldii (5°C, 15 h, 0.70 described in the literature (Fig. 1). Reported irradiance –2 –1 –2 –1 –2 –1 mol quanta m h ; 10°C, 15 h, 0.11 mol quanta m h ; 15°C, values were converted to units of mol quanta m h 9 h, 0.70 mol quanta m–2 h–1); Beardall & Morris (1976): following Richardson et al. (1983). Data from a given Phaeodactylum tricornutum; Chan (1978): Chaetoceros sp., study were binned according to species, photoperiod, Skeletonema costatum, Cylindrotheca fusiformis, Thalas- –2 –1 and growth temperature, and then each bin was siosira floridana (0.03 mol quanta m h ), Gymnodinium sim- plex, Amphidinium carterae; Yoder (1979): Skeletonema divided by a scalar to yield normalized chlorophyll costatum (22°C, 0.25 mol quanta m–2 h–1); Falkowski & Owens concentrations (chlnorm) (see Fig. 1 legend). For the 342 (1980): Skeletonema costatum, Dunaliella tertiolecta; Fal- observations in our data set, the relationship between kowski (1980): Dunaliella tertiolecta, Skeletonema costatum; Falkowski et al. (1981): Dunaliella tertiolecta; Verity (1981): chlnorm and Ig was described by: Leptocylindrus danicus (15°C, 0.23 mol quanta m–2 h–1); × –1.1 Ig chlnorm = 0.036 + 0.2 × e (1) Cosper (1982a): Skeletonema costatum; Cosper (1982b): Ske- etonema costatum; Faust et al. (1982): Prorocentrum mariae- Eq. (1) expresses cellular chlorophyll as a function of lebouriae; Raps et al. (1983): Microcystis aeruginosa; Terry light (Fig. 1). If Calvin cycle capacities (cell–1) had been et al. (1983): Phaeodactylum tricornutum; Geider et al. simultaneously measured, then the parameters in (1985): Phaeodactylum tricornutum; Post et al. (1985): Oscilla- toria agardhii; Dubinski et al. (1986): Thalassiosira weisflogii, Eq. (1) could be adjusted to describe light-dependent Isochrysis galbana, Prorocentrum micans; Sukenik et al. changes in chlorophyll relative to Pmax (i.e. chlrel). (1987): Dunaliella tertiolecta; T. Fisher (unpubl. data): Tetra- Unfortunately, Pmax measurements were not made. edron minimum, Nannochloropsis sp. 106 Mar Ecol Prog Ser 228: 103–117, 2002

b However, if Pmax exhibits a monotonic dependence on (Fig. 2b: inset). Below the nutricline, P max was esti- Ig, then chlrel should still follow the relationship: mated using the nutrient-sufficient equations (Fig. 2a). × Parameterizing the PhotoAcc model with field data chl = a + b × e–c Ig (2) rel required an environmental index of nutrient status and For the PhotoAcc model, we assumed that Eq. (2) pro- an assessment of surface mixing depths. Our approach vides an adequate description of chlrel. We then to these issues and the field data employed are de- derived values for the parameters (a, b, c) using field tailed in the 3 subsections that follow the brief model b measurements of P max. summary below. For mixed layer, nutrient-sufficient phytoplankton To recapitulate: the PhotoAcc model describes light-

(Condition 1), the PhotoAcc model attributes variabil- dependent changes in chlorophyll and Pmax for 3 envi- b ity in P max solely to changes in Ig. Since chlorophyll is ronmental conditions (Fig. 2). For each condition, expressed relative to the Calvin cycle capacity, Pmax changes in Pmax and chlrel are expressed relative to –3 –1 was simply assigned an irradiance-independent value the Pmax value of 1 mg C m h assigned to nutrient- –3 –1 of 1 mgC m h (Fig. 2a). Light-dependent changes sufficient, mixed-layer phytoplankton. Chlrel is assumed in chlrel were then determined by fitting Eq. (2) to to vary with light in a similar manner as chlorophyll per b the inverse of measured P max as a function of Ig (i.e. cell, which has been repeatedly measured in the labo- b –1 b chlrel = Pmax × P max = 1/P max), where Ig was taken as ratory (Fig. 1). Calvin cycle capacities are not directly the average daily PAR at the bottom of the active measured in the field, so the model is parameterized in b mixed layer. a stepwise manner using measurements of P max. The For nutrient-sufficient phytoplankton below the first step is to fit the light-chlrel relationship (Eq. 2) to b mixed layer (Condition 2), Ig was calculated as the inverse P max values observed in the mixed layer under average daily PAR at the depth of interest and chlrel nutrient-sufficient conditions. Applying this same was assigned the same relationship with Ig as deter- light-chlrel relationship to the other 2 conditions allows mined for the mixed layer. Relative to Condition 1, Pmax the remaining Pmax equations to be solved. Specifically, was assumed to decrease at very low light below the the light-dependent change in Pmax below the mixed mixed layer (Fig. 2a: inset) because a high Calvin cycle layer (Fig. 2: insets) is determined by fitting Eq. (3) to b capacity is unnecessary in chronically light-limited the product of chlrel and measured P max values from phytoplankton deep within a stratified water column below the mixed layer. The lower Pmax value for nutrient- (Geider et al. 1986, Orellana & Perry 1992). This light- depleted conditions (Fig. 2b) is determined as the b b dependent change in Pmax was described by: mean of the products, chlrel × P max (where P max values × are measured in nutrient-depleted waters). P = d + f × (1 – hn Ig)(3) max Field data. The PhotoAcc model was parameterized where n has inverse units of Ig to yield a dimensionless using Atlantic Meridional Transect (AMT) data. Each exponent. Parameter values (d, f, h, n) were derived by AMT field campaign involves a suite of complimentary fitting Eq. (3) to the product of chlrel and measured physical, chemical, and biological observations from b P max as a function of Ig. a diversity of oceanographic regimes. We employed Few laboratory studies have specifically investigated results from 2 of these program transects, AMT-2 (22 nutrient-dependent changes in the relative abundance April to 22 May 1996) and AMT-3 (16 September to of light harvesting and Calvin cycle components. We 25 October 1996) each spanning an 11 000 km region b assumed that P max would be lower under nutrient-de- between the United Kingdom and the Falkland pleted conditions because associated increases in the Islands. In addition to encompassing 5 of the major energetic cost of extracting nutrients from the environ- biogeochemical provinces defined by Longhurst et al. ment would favor the use of light reaction products for (1995), these AMT studies were also ideal for model ATP generation rather than carbon fixation. This shift in parameterization because they exhibited significantly b the balance between light harvesting and carbon fixa- divergent latitudinal patterns in P max despite being tion was effectuated in the PhotoAcc model by assum- nearly identical geographical transects (Marañón & ing the same light-chlrel relationship as under nutrient- Holligan 1999, Marañón et al. 2000). sufficient conditions, but a lower, light-independent Protocols employed for optical and biological mea- value for Pmax in the mixed layer (Fig. 2b). This nutrient- surements during AMT-2 and AMT-3 have been previ- depleted Pmax value was determined by multiplying ously reported (Robins et al. 1996, Marañón & Holligan b measured P max and calculated chlrel values and taking 1999, Marañón et al. 2000). Optical measurements the median of these products. Between the mixed layer were conducted using the SeaWiFs Optical Profiling and nutricline, Pmax was further reduced as a function of System (SeaOPS) with a set of 7-channel light sensors light following the light-dependent relationship de- (Robins et al. 1996). Optical measurements with the rived for nutrient-sufficient cells below the mixed layer SeaOPS have a characteristic uncertainty of 3.4%, Behrenfeld et al.: Modeling phytoplankton photosynthesis 107

Fig. 2. PhotoAcc model relationships between growth irradiance (Ig) and chlorophyll concentration, Calvin cycle capacity (Pmax), b and chlorophyll-normalized light-saturated photosynthesis (P max). (a) Nutrient-sufficient conditions; (b) nutrient-depleted con- ditions. Inset in (a): light-dependent decrease in Pmax for nutrient-sufficient growth below the mixed layer; inset in (b): light- dependent decrease in Pmax for nutrient-depleted phytoplankton between nitracline and bottom of mixed layer. The light- –2 –1 chlorophyll relationship is the same for all growth conditions and is normalized here to a value of 1 at Ig = 0 molquanta m h . –3 –1 b Pmax is expressed relative to a value of 1 mgC m h for mixed-layer, nutrient-sufficient conditions; P max is expressed relative to nutrient-sufficient, high-light conditions with 1.5% of this from ship-shadow contamination 20 m generally decreased due to longer periods of light- (Hooker & Maritorena 2000). Photosynthesis-irradiance limited photosynthesis. On a few occasions, the P b value measurements were conducted on samples collected at at 20 m greatly exceeded values immediately above and each station from 7 m, near the base of the mixed layer, below this depth, in which case the second highest P b b and near the chlorophyll maximum (Marañón & Holli- value was taken as P opt. Mixed-layer primary produc- b gan 1999). This sampling strategy provided P max tion is phosphate-limited at the BATS site during the values from within and below the mixed layer. Data summer months (Michaels et al. 1994, 1996). collected between 40 and 49° N during AMT-2 (3 sta- Nutrient limitation index. An environmental index tions) were not used in our analysis because of in- was required to discriminate between nutrient-sufficient consistencies in parallel chlorophyll measurements and nutrient-depleted water masses. The dissolved (specifically, the 2 records had similar patterns, but concentration of nutrients does not provide an effec- were juxtaposed by 1 station), nor did we use the tive index because, over a majority of the open ocean, AMT-3 data from 60 m at 38° S because of a large dis- the daily allocation of resources available for growth is crepancy between chlorophyll concentration and spec- largely determined by food-web recycling rates. We tral attenuation coefficients (Kdλ). therefore used the relationship between the physical Performance of the PhotoAcc model was additionally mixed layer depth (ZM) and the nitracline depth (ZN) as tested using a 6 yr record of light and primary production a first-order index of nutrient status. The basis for this from the US-JGOFS Bermuda Atlantic Time Series index is that when physical forcing (e.g. a storm) deep- (BATS) program and Bermuda BioOptics Program ens the mixed layer such that mixing penetrates the (BBOP). BATS and BBOP measurement protocols have nutricline, the upper water column becomes ‘charged’ been described previously (Knap et al. 1993, Michaels with nutrients and either a close correspondence exists

& Knap 1996, Siegel et al. 1995a,b, 2000). Chlorophyll- between ZM and ZN or nutrients are elevated through- normalized rates of primary production (P b) were based out the water column. When stratification follows, the 14 on light-dark in situ C-uptake measurements con- mixed layer shoals, and ZM separates from ZN. Simulta- ducted from sunrise to sunset and HPLC determinations neously, recycling inefficiencies export nutrients out of of phytoplankton pigment. Each vertical profile of P b was the active biological pool and into particulate refrac- visually inspected to extract the appropriate value for tory and deep-ocean pools (Eppley & Peterson 1979). If b b b P opt. With few exceptions, P opt was taken as the P the time scale of these 2 events is similar, then the sep- value measured at 20 m. Above this depth, a slight pho- aration of ZN and ZM should function as a correlate to toinhibitory effect was common, while P b values below increasing nutrient stress. 108 Mar Ecol Prog Ser 228: 103–117, 2002

– Clearly, conditions will exist when ZN > ZM, yet recy- depth at which NO3 was first detected. For nearly the cling is sufficient to transiently maintain phytoplankton entire AMT-3 transect, either a close correspondence ex- – in a nutrient-sufficient physiological state. The relation- isted between ZM and ZN, or NO3 was detectable to the ship between ZN and ZM was thus treated as a first-order surface (Fig. 3a). The AMT-3 data were thus used to pa- index to separate nutrient-sufficient from nutrient-de- rameterize the nutrient-sufficient PhotoAcc equations pleted data. For this analysis, ZM was evaluated from (Fig. 2a). During AMT-2, a general separation between vertical profiles of density (σt) and ZN was assigned the ZM and ZN was observed from 30° N to 42° S (Fig. 3b).

Fig. 3. Results for the Atlantic Meridional Transect studies: (a, c, e) AMT-3 and (b, d, f) AMT-2. (a, b) Physical mixed layer (ZM: r, r) and nitracline (ZN: y) depths; gray symbols identify central section of AMT-2 transect that was treated as nutrient-depleted; – – ZN was taken as the depth at which NO3 was first detectable or assigned a value of 0 m when NO3 was elevated throughout the b water column. (c–f) Measured (h) and modeled (D, D) chlorophyll-normalized light-saturated photosynthesis (P max: mgC mg –1 –1 b b chl h ); (D) P max estimated using the nutrient-sufficient PhotoAcc model (Fig. 2a); (D) P max estimated using the nutrient- depleted PhotoAcc model (Fig. 2b). (c, d) Mixed-layer samples collected at 7 m depth. (e, f) Mid-depth samples collected near ZM Behrenfeld et al.: Modeling phytoplankton photosynthesis 109

– Fig. 4. Exemplary vertical profiles of photosynthetically active radiation (PAR) (——), chlorophyll concentration (F), NO3 (s), 3– PO4 (h), and sigma-t (σt: D) used to identify photoacclimation mixing depths (ZACCL: ). (a) Nutrient-depleted water column – 3– where NO3 and PO4 concentrations are exhausted well below the mixing depth, indicated by σt; in this case, ZACCL was taken as the inflection depth in the σt profile. (b) Typical nutrient ‘charged’ water column where all profiles indicate a similar depth for – ZACCL. (c) Example of infrequent condition where σt was vertically uniform but NO3 and chlorophyll profiles indicated shallower – mixing; in this case, ZACCL was set to the depth indicated by the NO3 profile. These exemplary profiles were extracted from the BATS data set (31° 50’ N, 64° 10’ W) on (a) 14 October 1992, (b) 17 May 1994, and (c) 12 December 1996

Data from this region were used to parameterize the nu- RESULTS trient-depleted model (Fig. 2b). Photoacclimation depths. The PhotoAcc model Parameterized PhotoAcc model requires an assessment of the mixing depth to which surface phytoplankton are photoacclimated (ZACCL). The first step in model parameterization was to de- We determined ZACCL by visually inspecting vertical fine the light-chlrel relationship for nutrient-sufficient, – 3– profiles of σt, NO3 , PO4 , and chlorophyll for each mixed-layer phytoplankton (Fig. 2a). This was accom- –2 –1 sampling date and location, and assigned ZACCL the plished by calculating Ig (mol quanta m h ) at ZACCL shallowest depth at which a vertical gradient was (m) using measured surface PAR (mol quanta m–2 h–1) observed in any of these 4 water-column properties. and spectral downwelling attenuation coefficients (Kdλ: –1 Invariably, ZACCL corresponded to the depth of ZM m ) and then fitting Eq. (2) to the inverse of mixed b 2 under nutrient-depleted conditions (Fig. 4a). Typically, layer P max values measured during AMT-3 (r = 0.62). – 3– σt, NO3 , PO4 , and chlorophyll all indicated a similar The resultant relationship was: depth for Z under nutrient-sufficient conditions ACCL –3 × Ig chlrel = 0.036 + 0.3 × e (4) (Fig. 4b). However, occasionally when ZM was large, nutrient and chlorophyll profiles indicated a shallower where the units (from Eq. 2) are: a = mgchl m–3, b = –3 2 –1 depth for ZACCL than σt (Fig. 4c). This infrequent condi- mgchl m , and c = m h molquanta . tion probably corresponds to a period shortly after a Eq. (4) was used as the light-chlrel relationship for deep mixing event when physical forcing on the sys- nutrient-sufficient conditions both above and below tem has relinquished and the upper water column is the mixed layer. The light-dependent relationship for beginning to stratify. Pmax below the mixed layer (Fig. 2a: inset) was para- 110 Mar Ecol Prog Ser 228: 103–117, 2002

b b meterized by multiplying P max values measured be- of P opt (~20 m) and minimal photoinhibition ensured b low ZACCL during AMT-3 by calculated chlrel values that light-saturated photosynthesis (P max) was main- and then fitting Eq. (3) to these products as a function tained during each measurement over most of the 2 b of Ig (r = 0.72), where Ig is the average daily PAR at photoperiod. P max was calculated for each month each sampling depth. The resultant relationship was: using both the nutrient-sufficient and nutrient- depleted models. The nutrient-sufficient model pro- × × –9 Ig Pmax = 0.1 + 0.9 [1 – (5 10 ) ](5) b vided the best estimates of P opt during and shortly where the units (from Eq. 3) are: d and f = mgC m–3 h–1, following the seasonal deep mixing events (black h = dimensionless, and n = 1 m2 h molquanta–1. symbols in Fig. 5a), which were evidenced in the pro-

For nutrient-depleted conditions, Eq. (4) was again files of ZM and ZN (Fig. 5b). These winter mixing used to describe light-dependent changes in chlrel. events were followed each year by a period of b Mixed-layer P max values measured during AMT-2 increasing stratification, whereby ZM gradually sepa- between 30° N and 42° S (Fig. 3b) were then multiplied rated from ZN (Fig. 5b). Accordingly, the nutrient- by corresponding values of chlrel. The median of these depleted model generally provided the best estimates b products gave a nutrient-depleted value for Pmax of of P opt during the stratified periods (gray symbols in 0.40 mgC m–3 h–1 (SD = 0.15 mgC m–3 h–1). The addi- Fig. 5a). Overall, the PhotoAcc model accounted for b tional, light-dependent decrease in Pmax between the 83% of the variance in P opt for the BATS data set mixed layer and nutricline (Fig. 2b: inset) was calcu- (n = 79). lated by multiplying 0.40 mgC m–3 h–1 by Eq. (5). As discussed above in the ‘Methods’ (subsection In Eqs. (4) & (5), parameter values were chosen that ‘Nutrient limitation index’), the relationship between had the least number of significant digits yet provided ZM and ZN simply provides a first-order index of nutri- a fit to the observational data that was not significantly ent status. Our model results for the BATS time series b different from least-squares-derived values. P max is (Fig. 5a) clearly indicate that conditions can exist b calculated with the PhotoAcc model as the inverse of where ZM and ZN are separated, yet P max exhibits Eq. (4) for mixed-layer, nutrient-sufficient conditions. values consistent with the nutrient-sufficient model. b Below the mixed layer, P max is calculated as Eq. (5) Alternative, physiologically-based indices of nutrient divided by Eq. (4). Under nutrient-depleted conditions, status should therefore be investigated. For example, b –3 P max in the mixed layer is calculated as 0.40 mgC m photochemical quantum efficiencies (Fv/Fm) could pro- –1 b h divided by Eq. (4), and P max between the mixed vide a useful index and can now be routinely mea- layer and nutricline is calculated by multiplying 0.40 sured using stimulated fluorescence techniques (Kol- mgC m–3 h–1 by Eq. (5), then dividing by Eq. (4). ber & Falkowski 1993, Falkowski & Kolber 1995).

b Model performance Comparison with temperature-dependent P opt models b Comparison of modeled and measured P max values b for nutrient-sufficient regions of the AMT transects Most regional- and global-scale models for P opt are (black symbols in Fig. 3) indicated that an effective simple functions of sea surface temperature (excep- b 2 representation of spatial variability in P max (r = 0.72, tions by: Balch & Byrne 1994, Longhurst 1995), with n= 79) was achieved both above (e.g. Fig. 3c,d) and considerable variability between models (Behrenfeld & below (e.g. Fig. 3e,f) ZACCL. Likewise, the nutrient- Falkowski 1997b). The AMT and BATS studies encom- b depleted model captured spatial variability in P max passed a wide range of sea surface temperatures (5° to observed above the nitracline during AMT-2 between 28°C) and thus provided an opportunity to compare 30° N and 42° S (gray symbols in Fig. 3d,f) (r2 = 0.71, the PhotoAcc model with the temperature-dependent n = 41). Overall, the PhotoAcc model accounted for models of Megard (1972), Balch et al. (1992), Antoine b 76% of the variance in measured P max for AMT-2 and et al. (1996) and Behrenfeld & Falkowski (1997a) AMT-3 (n = 120). (Fig. 6). None of these earlier models accounted for b b A level of agreement between PhotoAcc estimates more than 9% of the variance in P max and P opt b of P max and values measured during the AMT stud- observed during the AMT and BATS studies (n = 199; ies was anticipated, since the AMT data were used Fig. 6a–d). In contrast, the PhotoAcc model accounted for model parameterization. Performance of the Pho- for 70% of the variance when switching between the toAcc model was therefore further tested using the nutrient-sufficient and nutrient-depleted equations independent, 6 yr BATS/BBOP data set (Fig. 5). For was strictly based on the correspondence between ZM b these data, a single P opt value was extracted from and ZN (Fig. 6e). Performance of the PhotoAcc model each monthly profile of Pz. The near-surface location was improved by 10% when the best-fit estimates Behrenfeld et al.: Modeling phytoplankton photosynthesis 111

Fig. 5. PhotoAcc results for the 6 yr Bermuda Atlantic Time-Series (BATS) record. (a) Measured (h) and modeled (D, D) optimal b –1 –1 2 b b chlorophyll-normalized carbon fixation (P opt: mgC mgchl h ) (r = 0.83); P opt was equated to P max value estimated with the b b PhotoAcc model; (D) P max estimated using nutrient-sufficient PhotoAcc model (Fig. 2a); (D) P max estimated using nutrient- depleted PhotoAcc model (Fig. 2b). (b) Relationship between physical mixed-layer depth (ZM: f, f) and depth of the nutricline – – (ZN: y); ZN was taken as the depth at which NO3 was first detected or assigned a value of 0 m when NO3 was elevated through- – 3– out the water column; NO3 was chosen as an index of nutrient-sufficient or nutrient-depleted conditions rather than PO4 because the later was often undetectable even when ZM was very large; (f) sampling dates when nutrient-sufficient PhotoAcc model yielded best results; (f) sampling dates when nutrient-depleted PhotoAcc model yielded best results were chosen (Fig. 6f). We conclude from this compari- related with causative environmental forcing factors. In son that: (1) the PhotoAcc model addresses significant contrast, the PhotoAcc model attempts to explicitly de- b b sources of variability in P max and P opt that are not ade- scribe primary causative relationships at the level of quately captured by simple temperature-dependent chlorophyll synthesis and changes in the Calvin cycle functions, and (2) optimal performance of the PhotoAcc capacity. With this approach, the model effectively re- model requires replacing our first-order index of nutri- produced spatial and temporal variability in light-satu- ent stress with a more quantitative physiological index. rated photosynthesis for the AMT and BATS studies. These field programs do not, however, fully encompass global nutrient, light, and temperature conditions. In DISCUSSION the following sections, we evaluate the predictions and hypotheses of the PhotoAcc model, discuss potential al- The importance of effectively modeling light- terations to account for additional growth constraints, saturated photosynthesis at the difficult local scale of and propose avenues for global implementation. daily primary production measurements has been rec- ognized for over 40 yr (Ryther 1956, Ryther & Yentsch 1957). The most common approach to this problem has Maximum light-saturated photosynthesis b been to describe P opt as a function of temperature. Such models capture only a fraction of the observed The PhotoAcc model predicts a maximum photosyn- b b –3 variability in P opt, because temperature is weakly cor- thetic rate of P max = 1/0.036 = 27.5 mgchl m (Eq. 4), 112 Mar Ecol Prog Ser 228: 103–117, 2002

corresponding to nutrient-sufficient, high light condi- suggest that PSUO2 can range from 300 to 5000 (Fal- tions. Falkowski (1981) proposed a maximum value of kowski et al. 1981), while τ may vary from ≥5 ms to 25 mgC mgchl–1 h–1 by assuming a minimum photo- slightly less than 2 ms (Falkowski et al. 1981, Behren- synthetic unit (PSUO2) size of 2000 chlorophyll mole- feld et al. 1998) depending on the excess capacity of τ cules per O2 evolved and a PSUO2 turnover time ( ) of the photosystems (Kok 1956, Weinbaum et al. 1979, τ 1 ms. Choosing a lower PSUO2 or value would yield a Heber et al. 1988, Leverenz et al. 1990, Behrenfeld et b b higher maximum for P max. Laboratory measurements al. 1998). Thus, a potential exists for P max to exceed 25 mgC mgchl–1 h–1, although values >20 mgC mg chl–1 h–1 are rarely ob- served in the laboratory (e.g. Glover 1980 reported values between 1.6 and 21 mgC mgchl–1 h–1 for 11 species]. In b b the field, P max (or P opt) values in excess of 20 mgC mgchl–1 h–1 are occasionally reported that exhibit well-defined di- urnal patterns (e.g. Malone et al. 1980, Hood et al. 1991). Thus, the PhotoAcc b –1 maximum for P max of 27.5 mgC mgchl h–1 is (1) close to the theoretical maxi- mum value proposed by Falkowski (1981), (2) well within the range of vari- ability suggested by laboratory mea- τ surements of PSUO2 and , and (3) consis- b tent with the maximum P max values reported from previous field studies.

Temperature revisited

We consider temperature to have a b negligible direct influence on P max above 5°C for the following reasons: (1) Laboratory studies with phyto- plankton monocultures indicate a clear temperature-dependence for maximum algal growth rates (Eppley 1972), but a variable, species-dependent influence on b P max (e.g. Steemann Nielsen & Hansen 1959, Steemann Nielsen & Jørgensen 1968, Harris 1978, Morris 1981, Post et al. 1985). Enzymatic temperature optima differ between species according to the ambient temperature from which they were isolated. Thus, in field samples with taxonomically diverse phytoplank- ton assemblages acclimated to their b b Fig. 6. Measured P max and P opt values (n = 199) versus modeled values from ambient conditions, temperature-depen- 2 temperature-dependent functions of (a) Megard (1972) (r = 0.08), (b) Balch et dence of P b is likely to be minimal. al. (1992) (r2 = 0.04), (c) Antoine et al. (1996) (r2 = 0.09), and (d) Behrenfeld & max Falkowski (1997a) (r2 = 0.01) and the PhotoAcc model when (e) switching (2) Although the activity of Calvin between the nutrient-depleted and nutrient-sufficient models was objec- cycle enzymes, such as ribulose-1-5- tively based on the relationship between mixed-layer and nutricline depths bisphosphate carboxylase (RuBisCO), 2 2 (r = 0.70) and (f) the best fit from either model was chosen (r = 0.80). Observa- exhibits classic Arrhenius temperature- tional data are from the Atlantic Meridional Transect studies (AMT-2, AMT-3) and Bermuda Atlantic Time-Series study (BATS). Dashed line indicates 1:1 rela- dependence, any decreases in activity at tionship. Equations for each of the temperature-dependent functions are given lower temperatures may be offset by in Behrenfeld & Falkowski (1997b) increasing enzyme concentrations (Stee- Behrenfeld et al.: Modeling phytoplankton photosynthesis 113

mann Nielsen & Hansen 1959, Steemann Nielsen & high light result from (1) dilution by coincident cell Jørgensen 1968, Geider et al. 1985, Geider 1987). In division and cessation of chlorophyll synthesis, (2) b other words, sensitivity of P max to temperature may be photooxidation, and (3) enzymatic breakdown. If these diminished if increases in enzyme concentration keep 3 processes combine to yield a kinetic response of sim- pace with decreasing activities at lower temperatures. ilar magnitude and opposite sign as a shift to low light (3) Temperature-dependent changes in the Calvin (Lewis et al. 1984), then phytoplankton acclimate to cycle capacity may be paralleled by changes in chloro- the average light in the mixed layer by continuously phyll (e.g. Durbin 1974, Yoder 1979, Verity 1981, synthesizing and degrading chlorophyll at time scales

Lapointe et al. 1984), such that normalization of Pmax of <1 h as they transit between the surface and ZACCL. to chlorophyll masks temperature dependence. Light-shift experiments, however, do not mimic nat- Below 5°C, the physiological adjustments described ural conditions well, since mixed-layer phytoplankton above may not be sufficient to eliminate temperature are preconditioned to periodic high-light exposures b effects on P max. Currently, the minimum mixed-layer and thus are physiologically poised to effectively dissi- b P max value predicted by the PhotoAcc model is 3 mgC pate excess excitation energy (Behrenfeld et al. 1998). –1 –1 b mgchl h . This value is near the maximum P opt Goericke & Welschmeyer (1992) demonstrated that value (<0.5 to ~3 mgC mgchl–1 h–1) reported by chlorophyll turnover is independent of light level in Dierssen et al. (2000) for Antarctic phytoplankton sam- preconditioned cells, and concluded that photooxida- b pled from –2 to 2°C waters. Low P opt values measured tion and enzymatic degradation of chlorophyll are in the Southern Ocean partly result from prolonged artifacts of light-shift experiments. Chlorophyll syn- periods of subsaturated photosynthesis causing a thesis at subsaturating light (Escoubas et al. 1995) b b weaker relationship between P max and P opt. Never- should thus dominate over the dilution effect of cell theless, a direct temperature effect is likely, and could division, causing photoacclimation to scale to a lower be accounted for in the PhotoAcc model by treating than average light level. This conclusion implies that

Pmax as a temperature-dependent variable in the nutri- cellular chlorophyll changes at the longer time scale ent-sufficient model. For example, the following rela- of mixed-layer deepening and shoaling. Additional tion (where T = temperature): analyses with mixed-layer phytoplankton are clearly needed to better resolve photoacclimation strategies P = 0.4 + 0.6 × [1 – exp–0.6 × (T + 6)]50 (6) max in nature. would provide a temperature-dependent decrease in b P max below 5°C while essentially leaving modeled val- ues above 5°C unaltered. Nutrient limitation

The PhotoAcc model simply switches between nutri-

Acclimation irradiance (Ig) ent-sufficient and nutrient-depleted equations rather than employing a more gradual transition as a function Modeling photoacclimation as a function of light at of nutrient stress. Falkowski et al. (1989) reported a the bottom of the mixed layer is the simplest formula- similar switch in the ratio of photosynthetic electron tion mathematically, but its physiological representa- transport (PET) components and RuBisCO for Isochry- tion is dependent on the kinetic balance between sis galbana grown under 9 different levels of nitrogen responses to low and high light. Specifically, if photo- limitation in chemostats. Increasing nutrient stress acclimation follows first-order kinetics with similar rate (i.e. decreasing growth rates) decreased cellular con- b constants for both high and low light, P max will vary as centrations of all measured components, with a con- a function of the average light in the mixed layer stant stoichiometry for growth rates between 0.59 and (Falkowski & Wirick 1981, Lewis et al. 1984, Cullen & 0.96 d–1. At growth rates <0.59 d–1, the relationship b Lewis 1988). Likewise, P max will scale to an irradiance between PET components and RuBisCO abruptly less than the average if responses to low light are changed by 60% to a new equilibrium (Fig. 7). The stronger (Cullen & Lewis 1988) and to higher than magnitude of this shift is the same as the decrease in average light if high-light effects dominate (Vincent et Pmax used in the nutrient-depleted PhotoAcc model. al. 1994, Geider et al. 1996). We are unaware of any physiological explanation for Photoacclimation rate constants are typically derived why such an abrupt change occurs, and few laboratory from laboratory ‘light-shift’ experiments, where cells studies have specifically addressed this problem. acclimated for multiple generations to a given irradi- Nutrient conditions represented by the AMT and ance are transferred to much higher or lower light (e.g. BATS studies included nutrient-sufficient, phosphate- Lewis et al. 1984). Riper et al. (1979) suggested that depleted (BATS) (Michaels et al. 1994, 1996), and pos- decreases in cellular chlorophyll following a shift to sibly nitrogen-depleted (AMT-2). We anticipate that 114 Mar Ecol Prog Ser 228: 103–117, 2002

sets (e.g. US National Oceanographic Data Center, Silver Spring, MD), but eventually data coincident with satellite ocean-color measurements will be pre- ferred. Coupled ocean models that assimilate satellite observational data (e.g. surface wind) could be used to generate mixed-layer depth fields. Alternatively, remote-sensing technologies may soon permit global surveys of surface mixing depths, and research is now underway investigating potential lidar techniques for such measurements. Combining remotely sensed sea- surface temperature fields with nutrient-depletion/ temperature relationships may prove useful for assess- ing temporal changes in nutrient status. Again, cou- pled ocean models may represent an alternative tech- nique. A variety of approaches may therefore be ex- plored for implementing the PhotoAcc model, and Fig. 7. Isochrysis galbana. Relative changes in photosynthetic each presents an exciting challenge for future re- electron transport (PET) components and RuBisCO for I. gal- search. bana maintained in nitrogen-limited chemostats at growth rates (µ) ranging from 0.18 to 0.96 d–1. Data are from Table 2 in Falkowski et al. (1989). PET components are: (y) Reaction Center II (RCII), (e) cytochrome-f, and (n) Reaction Center I; Perspective (D) ratio of RCII:RuBisCO. Data are normalized to average value for the 4 highest growth rates. If RuBisCO concentra- Deciphering decadal-scale changes in oceanic pri- tion is assumed proportional to light-saturated carbon fixa- mary production from the much larger amplitude sig- tion (Pmax) (Sukenik et al. 1987), then the decrease in RCII:RuBisCO between µ = 0.59 and 0.48 divisions d–1 is of the nals of seasonal and interannual variability requires same magnitude (60%) as the decrease in Pmax assigned to foremost a long-term commitment to remote-sensing nutrient-depleted conditions in the PhotoAcc model (Fig. 2b) programs monitoring changes in plant biomass, as well as process-oriented models of depth-dependent changes in phytoplankton biomass and assimilation an additional set of PhotoAcc equations will be neces- efficiencies. An from empirical to analytical b sary for modeling P max in iron-limited regions. Iron- models of phytoplankton photosynthesis has been limitation causes unique stoichiometric changes in realized primarily in the spectral description of the the composition of the photosystems, most notably a underwater light field (e.g. Platt & Sathyendranath decrease in cytochrome b6 f and a higher Photosystem 1988, Morel 1991, Platt et al. 1991, Antoine et al. 1996). II to Photosystem I ratio (PSII:PSI) (Guikema & Sher- Demonstrable improvements in model performance man 1983, Sandmann 1985, Greene et al.1992, Straus from these achievements have, unfortunately, been 1994, Vassiliev et al. 1995, Behrenfeld & Kolber 1999). masked by the overwhelming influence of ineffective A high PSII:PSI ratio will diminish chlorophyll-specific local-scale models of physiological variability (Balch efficiencies for converting photochemical energy to et al. 1992, Behrenfeld & Falkowski 1997a,b, Siegel et the reducing equivalents necessary for carbon fixa- al. 2000). The PhotoAcc model presented here pro- b tion, and should thus decrease P max relative to similar vides a framework for future improvements in the ana- levels of growth limitation by nitrogen or phosphate. lytical characterization of physiological variability and, based on comparisons with AMT and BATS data, already incorporates 2 critical forcing factors: light and Global implementation nutrients. The importance of environmental physical-chemical The immediate goal of the research described here variability to physiological traits of algae is well illus- was to develop an effective field-based model for trated in the BATS data by the correspondence b b P max. Our broader objective, however, is to establish a between changes in P opt and the temporal rhythm of b robust P opt algorithm for estimating global oceanic seasonal mixing and stratification (Fig. 5). Nutrient primary production from satellite ocean-color data. availability at this site clearly influences light-saturated Implementing the PhotoAcc model will require tem- photosynthesis at the seasonal scale, while higher fre- b porally resolved mixed-layer depth and nutrient- quency variations in P opt are dominated by photo- depletion fields. As an initial attempt, monthly climato- acclimation responses to changes in the mixed-layer logical fields can be assembled from historical data depth, surface PAR, and phytoplankton biomass (due Behrenfeld et al.: Modeling phytoplankton photosynthesis 115

to associated changes in vertical light attenuation). Beardall J, Morris I (1976) The concept of light intensity adap- The BATS data also suggest a significance in the tation in marine phytoplankton: some experiments with recent history of a water column, as occasional mid- Phaeodactylum tricornutum. Mar Biol 37:377–387 Behrenfeld MJ, Falkowski PG (1997a) Photosynthetic rates summer switches between nutrient-depleted and nutri- derived from satellite-based chlorophyll concentration. b ent-sufficient P opt values conceivably reflect the pass- Limnol Oceanogr 42:1–20 ing of a mesoscale eddy with a significantly different Behrenfeld MJ, Falkowski PG (1997b) A consumer’s guide history than the water masses sampled the month to phytoplankton primary models. Limnol Oceanogr 42:1479–1491 before and the month after. Behrenfeld MJ, Kolber ZS (1999) Widespread iron limitation The conceptual basis of our PhotoAcc model entails a of phytoplankton in the south Pacific Ocean. Science 283: variety of hypotheses regarding sources of variability 840–843 in light-saturated, chlorophyll-normalized photosyn- Behrenfeld MJ, Prasil O, Kolber ZS, Babin M, Falkowski PG thesis. These include a dominant kinetic response (1998) Compensatory changes in photosystem II electron turnover rates protect photosynthesis from photoinhibi- to low-light over high-light exposure, an irradiance- tion. Photosynth Res 58:259–268 dependence of the Calvin cycle capacity under chronic Chan AT (1978) Comparative physiological study of marine low-light conditions, a nutrient-dependent switch in diatoms and dinoflagellates in relation to irradiance and the stoichiometric balance between light-harvesting cell size. I. Growth under continuous light. J Phycol 14: 396–402 and Calvin cycle components, and a temperature- Cosper E (1982a) Effects of variations in light intensity on the b independence of P max above 5°C in natural phyto- efficiency of growth of Skeletonema costatum (Bacillario- plankton assemblages. Additional field and laboratory phyceae) in a cyclostat. J Phycol 18:360–368 studies will be necessary to test these hypotheses and Cosper E (1982b) Influence of light intensity on diel variation to expand the PhotoAcc model to encompass iron- in rates of growth, respiration, and organic release of a marine diatom: comparison of diurnally constant and fluc- limiting conditions, temperature effects at <5°C, and tuating light. J Plankton Res 4:705–724 possibly species-dependent variability in phytoplank- Cullen JJ, Lewis MR (1988) The kinetics of algal photoaccli- ton pigment composition. mation in the context of vertical mixing. J Plankton Res 10: 1039–1063 Dierssen HM, Vernet M, Smith RC (2000) Optimizing models Acknowledgements. M.J.B. thanks John Chapman, Ricardo for remotely estimating primary production in Antarctic Letelier, Paul Falkowski, Richard Geider, and Zbigniew Kol- coastal waters. Antarct Sci 12:20–32 ber for helpful discussions and inspiration, Elizabeth Stanley Dubinsky Z, Falkowski PG, Wyman K (1986) Light harvesting for encouragement, and the US National Aeronautics and and utilization by phytoplankton. Plant Cell Physiol 27: Space Administration (NASA) (Grant UPN161-35-05- 08) for 1335–1349 its continued support. BATS data were assembled by Toby Durbin EG (1974) Studies on the autecology of the marine Westberry and Margaret O’Brien. The US JGOFS BATS pro- diatom Thalassiosira nordenskiöldii Cleve: I. The influ- gram is supported by the US National Science Foundation. ence of daylength, light intensity, and temperature on D.A.S’s contributions were supported by NASA. E.M. thanks growth. J Phycol 10:220–225 Manuel Varela and Beatriz Mouriño for assistance during the Eppley RW (1972) Temperature and phytoplankton growth in AMT studies and Alan Pomroy, Colin Griffiths and Malcolm the sea. Fish Bull 70:1063–1085 Woodward for providing nutrient data. S.B.H. thanks J. Eppley RW, Peterson BJ (1979) Particulate organic matter flux Brown, C. Dempsey, S. Maritorena, and G. Moore for AMT and planktonic new production in the deep ocean. Nature support and J. Aiken for his diligence and commitment to 282:677–680 supporting high-quality AMT optical data. Eppley RW, Sloan PR (1966) Growth rate of marine phyto- plankton: correlations with light absorption by cell chloro- phyll. Physiol Plant 19:47–59 LITERATURE CITED Escoubas JM, Lomas M, LaRoche J, Falkowski PG (1995) Light intensity regulation of cab gene transcription is sig- Antoine D, André JM, Morel A (1996) Oceanic primary pro- naled by the redox state of the plastoquinone pool. Proc duction 2. Estimation at global scale from satellite (coastal Natl Acad Sci USA 92:10237–10241 zone color scanner) chlorophyll. Global Biogeochem Falkowski PG (1980) Light-shade adaptation in marine phyto- Cycles 10:57–69 plankton. In: Falkowski PG (ed) Primary production in the Balch WM, Byrne CF (1994) Factors affecting the estimate of sea. Plenum Publising Corporation, New York, p 99–119 primary production from space. J Geophys Res 99(C4): Falkowski PG (1981) Light-shade adaptation and assimilation 7555–7570 numbers. J Plankton Res 3:203–216 Balch W, Evans R, Brown J, Feldman G, McClain C, Esaias W Falkowski PG, Kolber Z (1995) Variations in chlorophyll (1992) The remote sensing of ocean primary productivity: fluorescence yields in phytoplankton in the world oceans. use of new data compilation to test satellite models. J Geo- Aust J Plant Physiol 22:341–355 phys Res 97(C2):2279–2293 Falkowski PG, Owens TG (1980) Light-shade adaptation. Banse K, Yong M (1990) Sources of variability in satellite- Plant Physiol 66:592–595 derived estimates of phytoplankton production in the Falkowski PG, Wirick CD (1981) A simulation model of the Eastern Tropical Pacific. J Geophys Res 95(C5):7201–7215 effects of vertical mixing on primary production. Mar Biol Beale SI, Appleman D (1971) Chlorophyll synthesis in Chlo- 65:69–75 rella: regulation by degree of light limitation of growth. Falkowski PG, Owens TG, Ley AC, Mauzerall DC (1981) Plant Physiol (Wash DC) 47:230–235 Effects of growth irradiance levels on the ratio of reaction 116 Mar Ecol Prog Ser 228: 103–117, 2002

centers in two species of marine phytoplankton. Plant response curve of green plant cells (Chlamydomonas rein- Physiol 68:969–973 hardtii). Planta 182:161–168 Falkowski PG, Sukenik A, Herzig R (1989) Nitrogen limita- Lewis MR, Cullen JJ, Platt T (1984) Relationship between tion in Isochrysis galbana (Haptophyceae). II. Relative vertical mixing and photoadaptation of phytoplankton: abundance of chloroplast proteins. J Phycol 25:471–478 similarity criteria. Mar Ecol Prog Ser 15:141–149 Faust MA, Sager JC, Meeson BW (1982) Response of Proro- Longhurst A (1995) Seasonal cycles of pelagic production and centrum mariae-lebouriae (Dinophyceae) to light of differ- consumption. Prog Oceanogr 36:77–167 ent spectral qualities and irradiances: growth and pig- Longhurst A, Sathyendranath S, Platt T, Caverhill C (1995) An mentation. J Phycol 18:349–356 estimate of global primary production in the ocean from Field CB, Behrenfeld MJ, Randerson JT, Falkowski PG (1998) satellite radiometer data. J Plankton Res 17:1245–1271 Primary production of the biosphere: integrating terres- Malone TC, Garside C, Neale PJ (1980) Effects of silicate trial and oceanic components. Science 281:237–240 depletion on photosynthesis by diatoms in the plume of Geider RJ (1987) Light and temperature dependence of the the Hudson River. Mar Biol 58:197–204 carbon to chlorophyll ratio in microalgae and cyanobacte- Marañón E, Holligan PM (1999) Photosynthetic parameters of ria: implications for physiology and growth of phytoplank- phytoplankton from 50° N to 50° S in the Atlantic Ocean. ton. New Phytol 106:1–34 Mar Ecol Prog Ser 176:191–203 Geider RJ, Osborne BA, Raven JA (1985) Light dependence of Marañón E, Holligan PM, Varela M, Mouriño B, Bale A (2000) growth and photosynthesis in Phaeodactylum tricornutum Basin-scale variability of phytoplankton biomass, produc- (Bacillariophyceae). J Phycol 21:609–619 tion, and growth in the Atlantic Ocean. Deep-Sea Res Part I Geider RJ, Osborne BA, Raven JA (1986) Growth, photosyn- 47:825–857 thesis and maintenance metabolic cost in the diatom Phaeo- Megard RO (1972) Phytoplankton, photosynthesis, and phos- dactylum tricornutum at very low light levels. J Phycol 22: phorus in Lake Minnetonka, Minnesota. Limnol Oceanogr 39–48 17:68–87 Geider RJ, MacIntyre HL, Kana TM (1996) A dynamic model Michaels AF, Knap AH (1996) Overview of the U.S. JGOFS of photoadaptation in phytoplankton. Limnol Oceanogr Bermuda Atlantic time-series study and the Hydrostation 41:1–15 S program. Deep-Sea Res Part II 43:157–198 Glover HE (1980) Assimilation numbers in cultures of marine Michaels AF, Knap AH, Dow RL, Gundersen K and 8 others phytoplankton. J Plankton Res 2:69–79 (1994) Seasonal patterns of ocean biogeochemistry at the Goericke R, Welschmeyer NA (1992) Pigment turnover in JGOFS Bermuda Atlantic time-series study 14 the marine diatom Thalassiosira weisflogii. I. The CO2- site. Deep-Sea Res Part I 41:1013–1038 labeling kinetics of chlorophyll. J Phycol 28:498–507 Michaels AF, Olson D, Sarmiento JL, Ammerman JW and 5 Greene RM, Geider RJ, Kolber Z, Falkowski PG (1992) Iron- others (1996) Inputs, losses, and transformations of nitro- induced changes in light harvesting and photochemical gen and phosphorous in the pelagic north Atlantic Ocean. energy conversion processes in eukaryotic marine algae. Biogeochemistry 35:181–226 Plant Physiol 100:565–575 Morel A (1991) Light and marine photosynthesis: a spectral Guikema JA, Sherman LA (1983) Organization and function model with geochemical and climatological implications. of chlorophyll in membranes of cyanobacteria during iron Prog Oceanogr 26:263–306 starvation. Plant Physiol 73:250–256 Morel A, Berthon JF (1989) Surface pigments, algal biomass Harris GP (1978) Photosynthesis, productivity, and growth: profiles, and potential production of the euphotic layer: the physiological ecology of phytoplankton. Arch Hydro- relationships reinvestigated in view of remote-sensing biol Spec Issues (Ergeb Limnol) 10:1–171 applications. Limnol Oceanogr 34:1545–1562 Heber U, Neimanis S, Dietz KJ (1988) Fractional control of Morris I (1981) Photosynthetic products, physiological state, photosynthesis by the QB protein, the cytochrome f/b6 and phytoplankton growth. Can Bull Fish Aquat Sci 210: complex and other components of the photosynthetic 83–102 apparatus. Planta (Berlin) 173:267–274 Myers J (1946) Culture conditions and the development of the Hood RR, Abbott MR, Huyer A (1991) Phytoplankton and photosynthetic mechanism. J Gen Physiol 29:419–427 photosynthetic light response in the coastal transition zone Orellana MV, Perry MJ (1992) An immunoprobe to measure of northern California in June 1987. J Geophys Res 96: Rubisco concentrations and maximal photosynthetic rates 14769–14780 of individual phytoplankton cells. Limnol Oceanogr 37: Hooker SB, Maritorena S (2000) An evaluation of oceano- 478–490 graphic optical instruments and deployment methodolo- Paasche E (1967) Marine plankton algae grown with light- gies. J Atmos Oceanic Techol 17:811–830 dark cycles. I. Coccolithus huxleyi. Physiol Plant 20: Knap AH and 11 others (1993) BATS methods — March 1993, 946–956 BATS method manual, version 3. US JGOFS Planning and Paasche E (1968) Marine plankton algae grown with light- Coordination Office, Woods Hole, MA dark cycles. II. Ditylum brightwellii and Nitzschia Kok B (1956) On the inhibition of photosynthesis by intense turgidula. Physiol Plant 21:66–77 light. Biochim Biophys Acta 21:234–244 Platt T, Sathyendranath S (1988) Oceanic primary production: Kolber Z, Falkowski PG (1993) Use of active fluorescence estimation by remote sensing at local and regional scales. to estimate phytoplankton photosynthesis in situ. Limnol Science 241:1613–1620 Oceanogr 38:1646–1665 Platt T, Sathyendranath S (1993) Estimators of primary pro- Lapointe BE, Dawes CJ, Tenore KR (1984) Interactions be- duction for interpretation of remotely sensed data on tween light and temperature on the physiological eco- ocean color. J Geophys Res 98:14561–14567 logy of Gracilaria tikvahiae (Gigartinales: Rhodophyta). II. Platt T, Caverhill C, Sathyendranath S (1991) Basin-scale esti- Nitrate uptake and levels of pigments and chemical con- mates of oceanic primary production by remote sensing: stituents. Mar Biol 80:171–178 the North Atlantic. J Geophys Res 96:15147–15159 Leverenz JW, Falk S, Pilström CM, Samuelsson G (1990) Post AF, deWit R, Mur LR (1985) Interaction between temper- The effects of photoinhibition on the photosynthetic light- ature and light intensity on growth and photosynthesis of Behrenfeld et al.: Modeling phytoplankton photosynthesis 117

the cyanobacterium Oscillatoria agardhii. J Plankton Res Steemann Nielsen E, Jørgensen EG (1968) The adaptation of 7:487–495 plankton algae. I. General part. Physiol Plant 21:401–413 Raps S, Wyman K, Siegelman HW, Falkowski PG (1983) Stitt M (1986) Limitation of photosynthesis by carbon metabo- Adaptation of the cyanobacterium Microcystis aeruginosa lism, I. Evidence for excess electron transport capacity in to light intensity. Plant Physiol 72:829–832 leaves carrying out photosynthesis in saturating light and Richardson K, Beardall J, Raven JA (1983) Adaptation of uni- CO2. Plant Physiol 81:1115–1122 cellular algae to irradiance: an analysis of strategies. New Straus NA (1994) Iron deprivation: physiology and gene Phytol 93:157–191 regulation. In: Bryant DA (ed) The molecular biology of Riper DM, Owens TG, Falkowski PG (1979) Chlorophyll turn- cyanobacteria. Kluwer Academic, New York, p 731–750 over in Skeletonema costatum, a marine plankton diatom. Sukenik A, Bennett J, Falkowski PG (1987) Light-saturated Plant Physiol 64:49–54 photosynthesis-limitation by electron transport or carbon Robins DB, Bale AJ, Moore GF, Rees NW and 6 others (1996) fixation? Biochim Biophys Acta 891:205–215 AMT-1 cruise report and preliminary results. In: Hooker Talling JF (1957) The phytoplankton population as a com- SB, Firestone ER (eds) NASA Tech Memo 104566, Vol 35. pound photosynthetic system. New Phytol 56:133–149 NASA Goddard Space Flight Center, Greenbelt, MD Terry KL, Hirata J, Laws EA (1983) Light-limited growth of Ryther JH (1956) Photosynthesis in the ocean as a function of two strains of the marine diatom Phaeodactylum tricornu- light intensity. Limnol Oceanogr 1:61–70 tum Bohlin: chemical composition, carbon partitioning, Ryther JH, Yentsch CS (1957) The estimation of phytoplank- and the diel periodicity of physiological processes. J Exp ton production in the ocean from chlorophyll and light Mar Biol Ecol 68:209–227 data. Limnol Oceanogr 2:281–286 Vassiliev IR, Kolber Z, Wyman KD, Mauzerall D, Shukla VK, Sandmann G (1985) Consequences of iron deficiency on Falkowski PG (1995) Effects of iron limitation on photo- photosynthetic and respiratory electron transport in blue- system II composition and light utilization in Dunaliella green algae. Photosynth Res 6:261–271 tertiolecta. Plant Physiol 109:963–972 Siegel DA, Michaels AF, Sørensen J, O’Brien MC, Hammer Verity PG (1981) Effects of temperature, irradiance, and day- MA (1995a) Seasonal variability of light availability and length on the marine diatom Leptocylindrus danicus its utilization in the Sargasso Sea. J Geophys Res 100: Cleve: I. Photosynthesis and cellular composition. J Exp 8695–8713 Mar Biol Ecol 55:79–91 Siegel DA, O’Brien MC, Sørensen JC, Konnoff D, Fields E Vincent WF, Bertrand N, Frenette JF (1994) Photoadaptation (1995b) BBOP data processing and sampling procedures. in intermittent light across the St. Lawrence estuary fresh- US JGOFS Planning Rep No 19, US JGOFS Planning and water-saltwater transition zone. Mar Ecol Prog Ser 110: Coordination Office, Woods Hole, MA 283–292 Siegel DA, Westberry TK, O’Brien MC, Nelson NB and 10 Weinbaum SA, Gressel J, Reisfeld A, Edelman M (1979) others (2000) Bio-optical modeling of primary production Characterization of the 32,000 Dalton chloroplast mem- on regional scales: the Bermuda biooptics project. Deep- brane protein: probing its biological function in Spirodela. Sea Res Part II 48:1865–1896 Plant Physiol 64:828–832 Steemann Nielsen E (1952) The use of radio-active carbon Wright JC (1959) Limnology of Canyon Ferry Reservoir: phyto- (C14) for measuring organic production in the sea. J Cons plankton standing crop and primary production. Limnol Int Explor Mer 18:117–140 Oceanogr 4:235–245 Steemann Nielsen E, Hansen VK (1959) Light adaptation in Yoder JA (1979) Effect of temperature on light-limited growth marine phytoplankton populations and its interrelation and chemical composition of Skeletonema costatum (Bacil- with temperature. Physiol Plant 12:353–370 lariophyceae). J Phycol 15:362–370

Editorial responsibility: Otto Kinne (Editor), Submitted: October 24, 2000; Accepted: May 10, 2001 Oldendorf/Luhe, Germany Proofs received from author(s): February 12, 2002