A Review with Some General Lessons from

1 downloads 0 Views 5MB Size Report
J A S ON D. Oj FM A M J J A S ON D. Years 1990-1994. Years 1990-1994 ...... Striven and Fyne, Proc. R. Soc. Edinburgh, Sect. B, 90,. 223-238, 1986. Therriault ...
PHYTOPLANKTON ECOSYSTEMS' FROM

BLOOM A REVIEW

SUSTAINED

DYNAMICS

WITH

SOME

INVESTIGATION

IN COASTAL GENERAL

OF SAN

LESSONS

FRANCISCO

BAY, CALIFORNIA James E. Cloern

U.S. GeologicalSurvey,Menlo Park, California

Abstract. Phytoplanktonbloomsare prominent features of biologicalvariability in shallowcoastalecosystems such as estuaries,lagoons,bays, and tidal rivers. Long-term observationand researchin San Francisco Bay illustratessome patterns of phytoplanktonspatial and temporalvariabilityand the underlyingmechanisms of thisvariability.Bloomsare eventsof rapid production and accumulationof phytoplanktonbiomassthat are usuallyresponsesto changingphysicalforcingsoriginating in the coastalocean (e.g., tides), the atmosphere (wind), or on the land surface(precipitationand river runoff).Thesephysicalforcingshavedifferenttimescales of variability,so algalbloomscanbe short-termepisodic events,recurrent seasonalphenomena,or rare events associatedwith exceptionalclimatic or hydrologicconditions.The biogeochemicalrole of phytoplanktonprimary productionis to transformand incorporatereactive inorganic elementsinto organic forms, and these transformationsare rapid and lead to measurablegeochemicalchangeduring blooms.Examplesinclude the depletionof inorganicnutrients(N, P, Si), supersaturation of oxygenand removal of carbondioxide,shiftsin the isotopiccompositionof reactiveelements(C, N), production of climaticallyactive trace gases (methyl

bromide,dimethylsulfide), changesin the chemicalform and toxicityof trace metals(As, Cd, Ni, Zn), changesin the biochemicalcompositionand reactivityof the suspendedparticulatematter, and synthesis of organicmatter requiredfor the reproductionand growthof heterotrophs, including bacteria, zooplankton, and benthic consumeranimals.Some classesof phytoplanktonplay specialroles in the cyclingof elementsor synthesisof specificorganicmolecules,but we haveonlyrudimentary understandingof the forcesthat selectfor and promote blooms of these species.Mounting evidence suggests that the natural cyclesof bloom variability are being altered on a global scaleby human activitiesincluding the input of toxic contaminantsand nutrients,manipulation of river flows, and translocationof species.This hypothesiswill be a key component of our effort to understandglobal changeat the land-seainterface. Pursuit of this hypothesiswill require creative approaches for distinguishing natural and anthropogenicsourcesof phytoplanktonpopulationvariability,aswell as recogni-

INTRODUCTION

where natural changeoriginatesfrom processesin the coastalocean,in the atmosphere,and on the land surface, and where human disturbanceis highly focused. Seventy-fivepercent of the U.S. population will live within 75 km of a coastby the year 2010 [Williamset al., 1991].As a resultof this densehuman settlementalong the continental margins,coastal ecosystemsare influencedby diversehuman activities,includingagricultural practices[Fleischeret al., 1987;Nixon, 1995]; the damming of riversand manipulationof their flows[Dynesius and Nilsson,1994];inputsof wastes,includingnutrients [Justif et al., 1995] and toxic contaminants[Windom, 1992];changingland use [Cooper,1995;Hopkinsonand Vallino, 1995] and habitatloss[Nicholset al., 1986];and disturbancesof biologicalcommunitiesthrough harvest [Rothschild et al., 1994]or introductionsof exoticspecies

Earth

science of the 1980s and 1990s has been moti-

vated partly by the challenge to understand global change,which in its broadestsenseincludesthe myriad impactsof the human populationon the Earth system. This problempersistsasa difficultchallengebecausethe componentelementsof global change,such as the climate system,hydrologiccycle, and biological populations,all haveinherentlargenaturalfluctuations.Superimposedonto thesenatural changesare thosecausedby humandisturbance.This scientificchallengepersistsbecause human disturbance is often subtle, indirect, and

sometimesconfoundedby natural changesthat themselvesare not well understood.This problem applies particularly to ecosystemsat the continental margins,

This paper is not subjectto U.S. copyright.

tion that

the modes

of human

disturbance

of coastal

bloom cyclesoperateinteractivelyand cannotbe studied as isolatedprocesses.

Reviewsof Geophysics,34, 2 / May 1996 pages 127-168

Publishedin 1996 by the AmericanGeophysicalUnion.

Paper number 96RG00986 e127e

128 ß Cloern: COASTAL PHYTOPLANKTON BLOOMS Table 1.

34, 2 / REVIEWS OF GEOPHYSICS

Examples of Recent Ecological Change in Global Coastal Waters Phenomenon

Location(s)

Episodesof anoxiaand associatedmortalities Baltic Sea of fish and shellfish

Adriatic Black

Sea

Sea

Harmful algal outbreaks global coastalwaters Oysterpopulationdeclinesand disappearance ChesapeakeBay of vascularplants Changesin the communitycompositionof Bay of Aratu Salvador, Brazil phytoplankton Dutch Wadden Sea Doubling of the biomassand shiftsin the communities

Reference (s)

Rosenberg et al. [1990] Justi•et al. [1987] Mee [1992] Smayda[1989],Hallegraeff[1993] Smithet al. [1992],Orth and Moore [1983] Cowgill[1987] Beukema[1991]

of invertebrates

Invasionby exoticinvertebratesand restructuringof biologicalcommunities

north San FranciscoBay

Interannual

Peel-Harveyestuary,Australia Laveryet al. [1991]

fluctuations in abundance and

Alpineand Cloem[1992],Kimmereret al. [1994]

speciesof macroalgae Mass mortalities

of diverse marine biota

Seabirdmortality Episodesof thickfoam accumulation on

Scandinavian coastal waters

Monterey Bay southern North

Sea

Underdalet al. [1989] Walzet al. [1994] BtitjeandMichaelis[1986]

beaches

Persistent closures of commercial harvest

shellfish

New Zealandcoastalwaters

Increasedfrequencyof red tidesand associated Tolo Harbor, Hong Kong fish kills SetoInland Sea,Japan

[Carlton, 1985; Carltonand Geller, 1993]. Nearly a fifth of the total marine

fish catch comes from these zones of

intensehuman activity[Paulyand Christensen, 1995]. In

recent

decades

we

have

observed

remarkable

changesin coastalwatersof all the continents;examples are listed in Table 1. Many of thesechangesare related either directly or indirectly to changesin the species composition,abundance,or productionrate of the phytoplankton, so this one biologicalcommunityis central to the problemof environmentalchangeat the land-sea interface. This theme is prominent in the recent initiatives to understandmechanismsof change in coastal ecosystems (Table 2). One fundamentalfeature of phytoplankton dynamicsis the episodicrapid population increase as events that we traditionally refer to as "blooms," presumablyin reference to Sch•itt's[1892] descriptionof the seasonalplanktoncyclesin Kiel Bight, (translationfrom Mills [1989,p. 125]):

Mackenzie et al. [1995] Hodgkiss and Yim [1995] Prakash[1987]

"and this play repeats itself year after year with the same regularityas every springthe trees turn green and in autumn losetheir leaves;withjust suchabsolutecertaintyasthe cherries bloom before the sunflowers, so Skeletonema arrives at their

yearly peak earlier than Ceratium."

My purpose here is to review some principles of phytoplanktonbloom dynamicsin their contextas features of changein estuarineand nearshorecoastalwaters.This reviewis organizedto addressbasicquestions: What are phytoplanktonblooms?What are their underlying mechanisms?How are they related to changesin ecosystemprocessesand the geochemistryof estuarine and coastal waters? The review is based around

obser-

vations,models, and insightsthat have come from the U.S. GeologicalSurvey(USGS) programof research focusedon San FranciscoBay. I begin with a general description and conceptual model of phytoplankton bloom dynamicsin shallowcoastalecosystems. Next is a

Table 2. Some Contemporary Research Initiatives to Describe and Explain Processesof Change in Coastal Ecosystemsof North America and Europe Program

Sponsor(s) *

NSF LMER (Land Margin EcosystemResearch) NSF, NOAA, ONR CoOP (CoastalOceanProcesses) NSF, NOAA GLOBEC (GlobalOceanEcosystems Dynamics) NECOP (NutrientEnhancedCoastalOceanProductivity) NOAA ECOHAB (Ecologyand Oceanography of HarmfulAlgal NSF, NOAA Blooms) ICSU LOICZ (Land-OceanInteractions in the CoastalZone) CEC EROS (EuropeanRiver OceanSystem)2000

Reference (s)

Boyntonet al. [1992] Brinket al. [1992] U.S. GLOBEC [1995] Wenzeland Scavia[1993],Atwood et al. [1994] Anderson[1995] Holliganand de Boois[1993] MartinandBarth[1989]

*Abbreviationsare NSF, U.S. National ScienceFoundation;NOAA, U.S. National Oceanicand AtmosphericAdministration;ONR, U.S. Office of Naval Research;ICSU, InternationalCouncilfor ScientificUnions;CEC, Commissionof the EuropeanCommunities.

34, 2 / REVIEWSOF GEOPHYSICS

Cloern: COASTAL PHYTOPLANKTON BLOOMS ß 129

ShallowCoastalEcosystem Local watershed ß..........

.. -..:.....ß

Shallows "" -.

Processesof phytoplankton population variability

ß.. '•....

•Zoca?•?us

Za/era/ •ran•?orZc '.i" ./

Longitudinal"• i7•'::. :.

Primary.production

.tr.•._..::•a..nsport:•.:..• -':-:?_:-_.-

'•:' -i)/:::: ':"'•:';• -,i-"2i:'?•:i":' -::" ............ eel'2gic• "'"• ':":'" Turbulent '"::':: ' •nputs [ ß"...Resusp,e.i•Sion •i•' • ,&_grazing mixing ... ß . ::'.•!•:..i,: Sedir•entation '....•.•"i•.'-'::.. -'-":I:'--')'" ":! -': '!i • -'-•-_d• Benthic grazing

Plate 1. An idealized shallowcoastal ecosystem(SCE) showingthe key processesof phytoplankton populationvariabilityat the land-seainterface.Localprocesses includeprimaryproduction(the sourceof new phytoplanktonbiomass),consumptionby benthicand pelagicgrazers,exchanges betweenthe sediment surfaceand water columnby sedimentationand resuspension, and turbulentmixing.Horizontal transport processes includeexchanges with the coastalocean,riverineinputsof phytoplankton, and longitudinaland lateral transportsgeneratedby density-,tide-, and wind-drivencirculationand mixing.The colortransition fromgreento bluerepresents the salinitygradientalongthe river-oceancontinuum.Brightness represents the bathymetry, andwhiteshadings indicatethelateralshallowdomainscommonto coastalplainestuaries andlagoons.

detailedreviewof the approachesusedin our analysisof the long term observationsin southSan FranciscoBay, with emphasison the searchfor patterns and mechanismsof phytoplanktonpopulationvariability.Then I describethe significanceof bloomsin the coastalzone, emphasizingthe role of phytoplanktonproductionas a mechanismof geochemicaland ecologicalchange.I conclude with some thoughtson the compellingquestions and hypothesesthat will guide researchat the land-sea interface into the next century. The USGS researchprogram in San FranciscoBay beganin 1968and continuesasone of the few long-term efforts of combined observation

and research in the U.S.

coastal zone. Other aspectsof this program are describedby Conomos[1979b],Kockelmanet al. [1982], Cloernand Nichols[1985],Nicholset al. [1986],Nichols and Pamatmat [1988], Peterson[1989], Cheng [1990], Luoma et al. [1990], and Petersonet al. [1995].

THE

GENERAL

PROBLEM

Domains of Interest

San FranciscoBay has been the focus of sustained investigationbecauseit has features common to many shallowcoastalecosystems (SCEs)that are influencedby natural and anthropogenicsourcesof variability.These SCEs include tidal rivers, estuaries, embayments,lagoons,and coastalriver plumes,which are distinctly different ecosystems from the open ocean.The differencesarise partly from the physicalfeaturesillustrated in Plate 1, and they includethe following. 1.

SCEs are transition

zones at the land-sea

inter-

face,sothey are influencedboth by inputsfrom the land surfaceand exchangeswith the coastalocean. Connections to the land surfaceare made through rivers that carryrunofffrom large-scalecatchmentsor from smaller tributariesthat carry runoff from the local watershed

130 ß Cloern: COASTAL PHYTOPLANKTON BLOOMS

34, 2 / REVIEWSOF GEOPHYSICS

tions.The phytoplankton(Table 3) include5000marine species[Hallegraeff,1993] of unicellularalgae havinga broad diversityof cell sizes(mostlyin the rangeof 1 to 100 txm), morphologies,physiologies, and biochemical compositions (Margalef[1978],Soumia[1982],andFogg 2. The domains considered here are shallow. The [1991]giveexcellentreviewsof the form and functionof mean depth of San FranciscoBay is less than 10 rn the phytoplankton).All phytoplanktonspeciesare capaand many have the capacityfor [Conomoset al., 1985], and maximum depth of the ble of photosynthesis, central channelsis of the order of 20-30 m. Exchanges rapid cell division and population growth, up to four of materialsbetweenthe pelagicdomain (open water doublingsper day [e.g.,Fahnenstielet al., 1995]. Popucolumn)and the benthos(bottomsediments)are rapid, lation dynamicsof the phytoplanktoncanbe interpreted so strongbenthic-pelagiccouplingis a definingfeature as responsesto changesin individual processesthat of SCEs. regulatethe biomass(total quantity,in measuressuchas 3. River-influencedSCEs are very differentphysical carbon,nitrogen,or chlorophyllconcentration),species environmentsfrom the open ocean.For example,turbu- composition,and spatialdistributionof the phytoplanklent mixingis a key physicalprocessthat determinesthe ton population.These processes includein situ (local) vertical fluxes of heat, salt, nutrients, and plankton. processesthat causepopulation changewithin a water Vertical mixing in the open ocean is regulatedby the parcel, and horizontal transportsthat displaceor mix seasonalcycle of heat input and thermal stratification water parcelsand their phytoplankton(Plate 1). that retards mixing between surfaceand deep waters. In situprocesses of populationchangeinclude(1) the However, in estuaries and other marine "regions of productionof new biomass,which is controlledby the freshwaterinfluence"[Simpsonet al., 1991],verticalmix- availabilityof visiblelight energyrequiredfor photosyning is regulatedby a larger and more variable sourceof thesisand nutrient resourcesrequiredfor the biosynthebuoyancy,the riverine input of fresh water that acts to sisof newalgalcells,(2) mortality,includingthat caused stabilizethe water columnthroughsalinitystratification. by parasitesor viruses,(3) grazinglossesto pelagicand 4. SCE's are particle-richrelativeto the open ocean. benthicconsumeranimals,(4) turbulentmixingby tideIn San FranciscoBay the near-surfaceconcentrationsof and wind-inducedmotions in the water column, (5) suspended particulatematter (SPM) rangefrom 300mgL-• [Conomos etal., 1979];at highest concen- bottom sediments,and (6) resuspension of bottom-detrations the SPM is dominatedby mineral particlesde- posited microalgaeby tidal currents and wind waves. livered by river flow or resuspendedoff the bottom by Horizontal transportsfollow water circulationsthat are tidal and wind wave currents.Suspendedparticles ab- driven by tidal currents, wind stresseson the water sorb and scatter light, so SCEs are turbid habitats in surface,and horizontal gradientsof water density [Fiwhichphytoplanktongrowthcan be limited by the avail- scheret al., 1979].Thesetransportsdisplacephytoplankabilityof sunlightto sustainphotosynthesis [Wofsy,1983; ton biomasslongitudinallyalong the river-oceancontinuum and laterally between shallowand deep domains, Cloem, 1987]. 5. Many SCEs are nutrient-rich becauseof inputs which are very different habitats for phytoplankton from the land surface[Maloneet al., 1988;Nixon, 1995; growth [Cloem and Cheng,1981; Malone et al., 1986]. Justif et al., 1995] and geochemicaland biologicalpro- Although techniquesexist to measureor estimateeach cesses that act as "filters" to retain nutrients within processshown in Plate 1, comprehensiveprocess-speestuaries[Sharpet al., 1984].In southSanFranciscoBay, cific measurementprogramsare expensive,logistically summerphosphateconcentrationsoften exceed10 gM challenging,and rarely done.Rather, phytoplanktondycomparedwith concentrationsof 10 m), but larger than in the micface is the tide, whichpropagatesinto San FranciscoBay rotidal Gulf of Mexico, Baltic Sea, and Mediterranean throughthe Golden Gate. In this regionof the northeast Sea (tidal amplitude< 0.5 m).

136 ß Cloern' COASTAL PHYTOPLANKTON BLOOMS

34, 2 / REVIEWS OF GEOPHYSICS The Issue of Scale

150OO a

•10000 5000

z

o 300

200 lOO

-• o "*•i '• 15

Each of the physicalforcingsillustratedin Figure 2 contributesto phytoplanktonpopulationvariabilityby influencingthe ratesof verticalmixing,horizontaltransport, production,or grazing.Each forcinghascharacteristictimescalesof variability,suchas the 12.42-hourand •-14-daytidal periods;the diel (24 hours)light cycle;2to 5-daystormeventsof enhancedstreamflowandwind stress;seasonalcyclesof irradiance and temperature; and the pronouncedinterannualvariabilityof river flow. Therefore phytoplanktonpopulationsin SCEs are exposedto physicalforcingsthat have timescalesof variability ranging from hours to years: "a hierarchy of forcingfunctionswhich drive the variousbiologicalresponsemechanismsat differenttime and length scales" [Mackasetal., 1985,p. 653].Each of thesetimescales can be identified in the populationfluctuationsof the phytoplankton. Biomass fluctuations

• -o

10

measured

5

•, 6o •

40

c:

20

0.7

0.5

0.4 1983

1986

1989

1993

Figure 2. Daily fluctuations of five physical forcings that influencepopulationdynamicsof phytoplanktonin San Francisco Bay. Series are for four contrastinghydrologicyears

(1983,1986,1989,and 1993),showing(a) net freshwaterinflow from the Sacramentoand San JoaquinRivers (a calculated quantity,the Delta Outflow Index, from the CaliforniaDepartmentof Water Resources), (b) dailydischarge of a localstream (PattersonCreek) that dischargesinto south San Francisco R•,x,(from 11R •oc•lc•glcml Suv:ey (L•, ni•_ trict), (c) meandailywind speedat the SanFranciscoairport (from the National Climatic Data Center, National Oceanic and AtmosphericAdministration),(d) daily irradiance,measured as quantumflux of photosynthetically activeradiation, and (e) m•imum daily tidal currentspeedin the channelof south San FranciscoBay, calculatedwith harmonicconstants from long-term current measurements[Chengand Ganner, 1985].

with

moored

at the short timescales fluorometers

that

can be

detect

and

record chlorophyllfluorescenceat a fixed location.A samplerecord in Figure 3a showsa 2-week seriesof hourly measurementsin south San FranciscoBay. The first week of the record showsperiodic fluctuationsin chlorophyllfluorescenceover the semidiurnaltidal period, with chlorophyllpeaksat the two low slacktides each day. This high-frequencyvariabilityis causedby tidal advectionas chlorophyllspatialgradientsoscillate over the sensorwith the tide [Cloemet al., 1989]. Spectral analysesof suchchlorophyllseriesconfirma high varianceat the tidal frequencies[Litaker et al., 1993]. Short-termfluctuationscan result from other processes suchasdiel cyclesof chlorophyllsynthesis andzooplankton grazing[Litakeret al., 1993] and wind-waveresuspensionof algalcellsoff the bottom[Demerset al., 1987; de Jongeand van Beusekom,1995]. The secondhalf of the record in Figure 3a showsa 6-day period of exponentialchlorophyllincrease.Here the semidiurnaltidal variabilityis overwhelmedby the rapid growthof phytoplanktonbiomassalongthe entire channel;this period of the record illustratesdaily-scalevariability during a bloom event. Day-to-dayfluctuationsin phytoplankton biomassor physiologicalconditionare commonlyassociatedwith hydrologic-meteorologic events,suchasrainfall pulses,wind events,or periodsof abrupt warming [C6t• and Platt, 1983], and with fluctuationsin tidal mixingoverthe neap-springperiod [Sinclairet al., 1981; Cloem, 199lb ]. Longer-termfluctuationsof phytoplanktonbiomass are illustratedin Figure 3b, a seriesof monthlychlorophyll a measurements in north San FranciscoBay from 1974 to 1995. Three scalesof variability are evident in this record: a seasonalcyclewith peak biomassduring the low-flow summers,interannual variability with a dampedsummerbloom duringyearsof extremelyhigh river flow and shortresidencetime (e.g., 1983), and an apparentpermanentchangein the natureof the record, with the virtual disappearanceof summerbloomsafter

34, 2 / REVIEWSOF GEOPHYSICS

Cloern' COASTAL PHYTOPLANKTON BLOOMS ß 137

40

•ø• o

i

Semi-diurnal

30 _

variability

'

Y



(D

oo

20

'•

10 0

1

2

4

5

7

i

i

i

i

i

•--"•----'-----':•;% daily variabilityF

' - ofa ,bo•t•,)•9

8

10

11

12

13

14

15

March,1995

b

'• 4o,-• /



30

o

20

e

10



Seasonal Interannual , ""'"':•:•:•• •"•'....'•• .:.• variabili I

variabili•

-

persistent .

-

o

1976

1978

198o

1982

1984

1986

1988

199o

1992

1994

1996

Figure3. Sometimescales of phytoplankton biomassvariabilityin SanFranciscoBay. (a) A 14-dayrecord of hourlymeasurements with an in situfluorometerplaced2 m abovethe bottomnear station28 in southSan FranciscoBay (Figure5); the fluorometerwascalibratedwith discretechlorophylla measurements made at the beginningandendof the record(fluorometerdatafrom D. A. Cacchione,personalcommunication, 1995). (b) A 21-yearrecordof monthly(or semimonthly) measurements of surfacechlorophylla concentration near station5 in north SanFranciscoBay (Figure 1); this recordincludesdata from the CaliforniaDepartmentof Water Resources(1974-1986) and USGS (1986-1995).

the North Bay was colonizedin large numbersby the clam Potamocorbulaamurensisin 1987.Jassbyand Powell [1994] examineda period of this chlorophyllrecord and concludedthat most of the interannualvariability is driven by two hydrodynamicprocesses,one associated with fluctuationsin river flow (a natural sourceof variability), and one associated with manageddiversionsof fresh water from the upper estuary(an anthropogenic sourceof variability). Spatialvariability(patchiness) is alsoobservedacross a spectrumof scales[Mackaset al., 1985].Phytoplankton patchinesscan be measuredwith continuousprofilesof chlorophyllfluorescence alonghorizontaltransects[Wilson and Okubo, 1980; Childerset al., 1994;A.D. Jassby et al., Towards the design of sampling networks for characterizingwater quality changesin estuaries,submitted to Estuadne,Coastal,and Shelf Science,1996]. Horizontal variability in south San Francisco Bay is illustratedin Figure 4, which showsthe spatialstructure of chlorophyllfluorescencealonga longitudinalchannel transect (Figure 4a) and along a transversetransect betweenthe channeland adjacentshallows(Figure 4b). The longitudinal profiles show measurementsspaced every 25 m, and variabilityis evidenteven at this small spatialscale.This small-scalevariabilityis superimposed onto larger-scalepatterns (mesoscalevariability) that

often showtrendsof decreasingchlorophyllfluorescence in the seawarddirection(February28, 1995;Figure 4a) and increasingchlorophyllacrossthe shallows(March 21, 1985; Figure 4b). Both small-scaleand mesoscale patchinessare influencedby hydrodynamicprocesses becausethe plankton are transportedby their fluid environment [Mackaset al., 1985]. Therefore the spatial variability patterns,both of chlorophyllbiomassand of individualspecies[e.g.,Kononenet al., 1992;Zingoneet al., 1995], are shapedby the turbulent advectivetransportsof phytoplanktonwithin a spatiallyvariablegrowth environment.These featuresof spatial structureare not stableor persistent.The mesoscaletrendschangefrom seasonto season[Glibert et al., 1995] and day to day [Wilsonand Okubo, 1980],while the small-scalevariability changesover the tidal period [Dustanand Pinckney, 1989]. Figures3 and 4 showthat phytoplanktonbiomassin shallow coastal ecosystemsvaries at timescalesfrom hours to decades and at spatial scalesfrom tens of meters to tens of kilometers.Levin [1992] suggeststhat the centralproblemin ecologyis the problemof pattern and scale,where pattern is the descriptionof spatial or temporalvariabilityand the mechanismsof pattern formation are scale-dependent. Therefore studiesof population fluctuationssuchasphytoplanktonbloomsrequire

138 ß Cloern' COASTAL PHYTOPLANKTON BLOOMS

34 2 / REVIEWS OF GEOPHYSICS b

• 50 A• ,T

"•B•-:•:•.•__• • • •

]

40

50







• D

40

'• 30

30

Fa 20

20-



10

•0

-

c

05

10

15

20

25

30

35

40

o

o



2



4



6

i

8



lO

12

Distance (kin)

Figure4. Horizontalvariabilityof phytoplankton biomass from continuous measures of chlorophyll fluorescence along(a) longitudinal transects and(b) a lateraltransectof southSanFrancisco Bay.Near-surface waterwaspumpedto a shipboardfluorometerthat wascalibratedwith 6-10 discretemeasurements of chlorophyll a concentration takenalongeachtransect.Insetmapsshowlocations of the transects.

explicitchoicesabout the scalesat whichpopulation variabilitycanbe observedandexplained.The long-term observationalprogramin San FranciscoBay was designedto characterize mesoscale spatialvariabilityalong the longitudinalaxis, at timescalesof weeks to years. Even though(large)populationvarianceexistsat other scales,the data set describedbelow is not appropriate 122ø15; '

122030 '

for analysisof the small-scalefluctuationsillustratedin Figures3 and 4. The SouthSan FranciscoBay Data Set

Patternsof Variability. South San FranciscoBay hasbeen the site of focusedresearchon phytoplankton 122000 '

.:.:.-......:•

Central

Bay...._

37045 '

Figure 5. Map of southSanFranciscoBayshowing locationsof fixedsamplingstationsalongthe longitudinalchannel.

37o30 '

34, 2 / REVIEWSOF GEOPHYSICS Table 4. Phytoplankton SpeciesCommonly Observed During Spring Blooms in South San Francisco Bay

Cloern: COASTAL PHYTOPLANKTON BLOOMS ß 139

the seawardand landwardestuary.Local streamscarry runoff to the lower (landward)estuary.The open connection to the central bay allows fresh water from the Class Species Sacramentoand San JoaquinRivers to intrude into the southbay during periodsof high river flow (Figure 1) Diatoms and allowstidal exchangebetween the South Bay and Diatomophyceae Chaetoceros debile Chaetoceros decipiens coastalocean. The basin'sresidual (tidally averaged) Chaetoceros didymus circulationis slow,with mean seawardflow along the Chaetoceros gracilis eastern shallowsand landward flow along the channel Chaetoceros socialis [Chengand Gartner, 1985] and a hydraulic residence Chaetoceros vistulae time of severalmonths[Walterset al., 1985]. This weak Chaetoceros wighami Coscinodiscus curvatulus(Actinocyclus mean circulationcanbe disruptedby strongwind events curvatulus) or freshwaterinputsthat alter the strengthand direction Coscinodiscus lineatus(Thalassiosira of the residualflows[Huzzeyet al., 1990].Water density leptopus) is often verticallyuniform, indicatingrapid verticalmixCoscinodiscus radiatus ing of the water column.The channelcan becomesalinCyclotellameneghiniana Cyclotellasp. ity stratified after periods of runoff that deliver lowCyclotellastriata density fresh water as a source of buoyancy, but Ditylum brightwelii stratificationoccursonly during neap tides when the Eucampiazoodiacus tidal stirringis weak [Cloern,1984]. PhytoplanktonpriLeptocylindrus minimus Nitzschiaseriata(Pseudo-nitzschia seriata) mary productionis the largestsourceof organiccarbon Paralia sulcata to south San FranciscoBay [Jassbyet al., 1993], and Rhizosoleniasetigera much of the total annual primary production occurs Skeletonema costatum during the spring[Coleet al., 1986]. Thalassiosira decipiens(Thalassiosira The southSan FranciscoBay springbloom has been angulata) Thalassiosira rotula followed every year since 1978, and the patterns of variability within this series can be used to illustrate Nondiatoms somegenerallessonsof phytoplanktonbloom dynamics Chlorella marina Chlorophyceae Chlorella salina in SCEs. The core measurementprogram.in south San Monoraphidiumconvolutum FranciscoBay includesvertical profilesof salinity,temNannochloris atomus perature, chlorophyll(by calibratedfluorometry),and Chrysophyceae Chromulinasp. turbidity at samplingstations(Figure 5) spacedabout Kephyrionsp. every3-4 km along the channel(detailed methodsare Ochromonassp. givenin annualdatareports[e.g.,Edmundsetal., 1995]). Cryptophyceae Chroomonas acuta ( Teleaulaxacuta) Chroomonas amphioxeia(Teleaulax Water samplesare also taken for microscopicexaminaamphioxeia ) tion to determinethe abundanceof phytoplanktonspeChroomonas salina (Rhodomonas salina) cies present during the springblooms (Table 4). The Dinophyceae Gonyaulaxtamarensis(Alexandrium sampling frequency in 1978-1979 was once monthly; ostenfeldii) after 1979the samplingfrequencywasincreasedto once Heterocapsatriquetra Katodinium rotundatum or twice weekly during the period of the springblooms. Prorocentrum minimum Bloom dynamicsare characterizedhere with the chloProtoperidiniumclaudicans rophyll measurementsmade along the channeltransect Prasinophyceae Pyramimonas micron(or P. orientalis) (Figure 5). Near-surfacechlorophyllis used as a meaTetraselmis gracilis Mesodinium rubrum sureof biomassin the photiczone (i.e., the quantitythat Photosynthetic ciliates contributesto primaryproduction),andverticalvariabilFrom samplestaken during 1992-1995 and microscopicenumera- ity is not consideredhere. Lateral variability is also not tions/identifications of R. G. Dufford (personal communications, considered,although there is coherence between the 1992-1995). Names in parenthesesare revisionsbasedon the taxon- seasonal dynamics of phytoplankton biomass in the omy of Tomas[1993, 1996]. channeland adjacentshallows[Cloernand Cheng,1981; Cloernet al., 1985; Vidergaret al., 1993]. Plate 2 depicts the chlorophyllfluctuationsin southSan FranciscoBay bloom dynamicsbecauseit has a recurrent, somewhat with color, and severalprominent patternsof variability predictableperiod of rapid biomassincreaseduringthe stand out. First, the solid backgroundshowsthat biospringmonths. This lagoon-likebasin has an irregular mass isusually lessthan5 IxgL- 1chlorophyll a. Second, bottom topographywith sharpbathymetrictransitions a period of rapid biomassgrowth (a bloom) occurs between the 10- to 25-m-deep axial channel and the duringthe springof everyyear. Third, bloom intensityis lateral shallows(Figure 5), and a large transverseshoal often greatest in the landward estuary,consistentwith (San Bruno) that slowshorizontalexchangesbetween the trends in the high-resolutionchlorophylltransects.

140 ß Cloern' COASTAL PHYTOPLANKTON

BLOOMS

34 2 / REVIEWS OF GEOPHYSICS

San Mateo Bri%ge Dumbarton Bridge

Ba/• Bridge 78

8O

82

84



86

88

9o

92

94

0

5

10

15

20

25

30

35

Distance along channel transect, in kilometers

....... I ....

5

10

15

I

I

I

I

I

20

25

30

35

40

Chlorophylla Plate 2. Color representationof spatial-temporalvariability of phytoplanktonbiomass (chlorophyll a concentration)in southSanFranciscoBay.Bloomsare shownasdeparturesfrom the low backgroundbiomass (green),with bloomintensityproportionalto color brightness. The verticalaxisrepresentstime from January 1978 to July 1995, and the horizontal axis representsspatialvariability along the channeltransectfrom the seawardestuary(station21) to the landwardestuary(station32). Color contourswere producedfrom a matrix of interpolatedvaluesbasedon 4449 measurementsof surfacechlorophylla concentrationon 417 dates.

34, 2 / REVIEWSOF GEOPHYSICS

Cloern' COASTAL PHYTOPLANKTON BLOOMS ß 141

1.0

PC1 = 52% of Variance 0.8

0.6

Figure 6. Coefficients(loadings)of the first two principalcomponents,which accountfor 87% of the chlorophyllvariance along the south San Francisco Bay transect, for the period January1978 to July 1995. Principal

0.4

0.2

component analysis wasdoneon •hecorrela-

1.0

tion matrix, and results are varimax-rotated

solutions[Jassby and Powell,1990].This analysissuggests two independentspatialmodes of phytoplankton biomass variability: one (PC1) in the landwardestuarinebasinsouth of station 25, and a second(PC2) in the seawardestuarynorth of station 25.

0.8

PC2 = 35% of Variance 0.6

0.4

0.2

21

22

23

24

25

26

27

28

29

30

31

32

Station

Finally, the nature of the springbloom changesfrom in south San FranciscoBay can be explainedwith only year to year. For example,the 1986 springbloomwas a two spatial modes (patterns): the first is expressed monthlongevent of very high chlorophyllconcentration stronglyin that region of the estuarylandwardof station (maximum >70 !xgL-•) alongtheentiretransect. The 25, and the secondis expressedstronglyin the region of 1991 springbloom was a short-livedevent confinedto the estuaryseawardof station25 (Figure 6). This sepathe landwardestuary,where maximumchlorophyllcon- ration of the two spatialmodesat station25 corresponds centration wasonly12!xgL-•. Thepatterns revealed in to the location of the San Bruno Shoal (Figure 5), a this image illustratetwo lessonsthat apply generallyto topographiccontrol of horizontal mixing. The tidally shallowcoastalecosystems: (1) phytoplankfon blooms averagedcirculationof southSan FranciscoBay is charare spatiallyheterogeneous,and (2) the timing and acterizedby a slowlyrotatinggyrethat actsto retain fluid magnitudeof the seasonalbloomschangefrom year to within the central basin of the estuary,below the San Bruno Shoal [Chengand Casulli,1982].The PCA result year. The patterns of population variability in time and is consistent with this mean circulation that maintains spacecan be formalized with multivariate techniques two spatialdomainsseparatedby a topographiccontrol such as principal componentsanalysis(PCA). If we of mixing [Powellet al., 1986]. This lessonapplies to considermeasurements at samplingstationsasvariables, other coastal ecosystemswhere the mesoscalespatial and samplingdates as individual cases,then PCA can variabilityof planktonreflectsthe water circulationpatreveal patterns of spatial coherencein the temporal ten [Jouffreet al., 1991] and where distinct spatial dofluctuationsof biomassalonga transect.This PCA of the mains can have different temporal patterns of bloom southSanFranciscoBay data setidentifiestwo principal evolution [Therriaultand Levasseur,1986; Kahru and componentsthat togetheraccountfor 87% of the total Nbmmann, 1990; Glibertet al., 1995]. A secondresult of PCA is the seriesof scalaramplichlorophyll variance. The first principal component (PC1), whichaccountsfor 52% of the chlorophyllvari- tudes(or scores)that expressthe relativeimportanceof and Powell, ance,hassmallcoefficients (or loadings)at the seaward eachprincipalcomponentover time [Jassby stationsand progressively largercoefficientsat the land- 1990].Figure 7 showsthe time seriesof amplitudesfor ward stations(Figure6). The secondcomponent(PC2), PC1. This series is a representationof pattern in the accountingfor an additional 35% of the variance,has temporal variability of phytoplanktonbiomassin the largest coefficientsat the seawardstations.This PCA landwardbasin.This temporal pattern is dominatedby showsthat mostof the phytoplanktonbiomassvariability episodicspikesthat occurin the springof eachyear,with

142 ß C!oern' COASTAL PHYTOPLANKTON BLOOMS

'•'

34, 2/REVIEWS OF GEOPHYSICS

30

o•.•,:20 lO I

I

I

I

I

I

I

I

I

I

I

I

6

Figure7. (bottom)Time seriesof amplitudes (scores)for the first principalcomponentof chlorophyllvariability in south San Francisco Bay. Large amplitudescorrespondto eventsof high phytoplanktonbiomass(blooms)in the landwardestuary.(top) Parallelseriesof nearsurfacesalinityin the landwardestuary(mean of measurementsat stations26-32), illustrating that bloomsoccurduringthe wet seasonof high river flow and annualsalinityminima.

m

m

E

1-0 m

_1 m

_2 /

78

79

80

81

82

83

84

85

86

87

88

89

90

91

92

93

94

95

Year

prominentyear-to-yearfluctuationin the timing, duration, and magnitudeof these spikes. Mechanismsof variability. The seriesof principal componentamplitudesin Figure7 is consistent with the notion of phytoplanktonblooms as departures from population quasi-equilibrium.Blooms are population outburststhat can appear and dissipatewithin a period of weeks.What determinesthe timing of theseoutbursts, their duration, and their magnitude,and why do these features change from year to year? The two series in Figure 7 showthat bloomsoccurin southSan Francisco Bay during the wet season,when surface salinity is dilutedby freshwater.This seasonal coherencebetween bloom developmentand low salinity results, in part, from the stabilizationof the water columnby the salinity stratification induced by freshwater inflows [Cloern, 1991a,hi. The absenceof bloomsin summer-autumnis

that thisbalanceis sensitiveto the rate of verticalmixing in the water column[Cloern,1991b].Vertical mixingis broughtaboutby tidal stresses appliedat the bottomand wind stresses on the water surface[Simpsonet al., 1991], and in mesotidal

SCEs the tidal stresses often dominate.

Weak tidal stirringpromotesbloomsby (1) slowingthe turbulent diffusivelossof phytoplanktonbiomassfrom the photiczone [Koseffet al., 1993], (2) deepeningthe photic zone (growth habitat) as tidal resuspension of bottomsediments weakens[Schoellhamer, 1996],and(3) reducingthe vertical flux of phytoplanktonto benthic consumeranimals[Cloern,1991b].Bloomsdevelopduring periods of weak tidal energy, and they dissipate during periods of strong tidal energy. This lessonis illustratedin Figure 8, which showsthe evolutionof the 1985 springbloom and the simultaneouschangesin the dailytidal currentamplitudeU. Rapid biomass(chloroalso a responseto the seasonalcYCles of the benthic phyll) growth occurredduring the neap tide in late suspensiOn-feeding animals,whosebiomassand grazing March; this bloom dissipatedsoonafter the subsequent rate are highest in summer. Therefore the seasonal springtide in earlyApril. distribution of bloomsin thisparticularestuaryappears The associationbetween weekly-scalechlorophyll vanatomyand the uuaxre -me•s a prominent........ of to be a responseto seasonatnuctuatlonsin nyurmogy (riverflowandits influenceon densitystratification)and phytoplanktonpopulationdynamicsin southSan Frantrophic interactions(grazinglossesto benthicconsum- cisco Bay. This associationis illustrated in Figure 9, which showsa linear relationshipbetween the rate of ers). At shorter timescales,too, the timing of blooms is biomasschangeand the antecedenttidal current speeds regulatedby physicalprocessesthat influencethe bal- for all the large springbloomsobservedbetween 1978 ance between phytoplankton production and losses. and 1995.The rate of biomasschange,R, wascalculated Simulation experimentswith numerical models show as ........

1



_•

_ •._

....

34, 2 / REVIEWSOF GEOPHYSICS

Cloern: COASTAL PHYTOPLANKTON BLOOMS ß 143

0.7 0.6 0.5 0.4 I

I

I

I

I

I

I

Figure 8. The 1985 springbloom in south San FranciscoBay: (bottom) weeklyfluctuations of phytoplanktonbiomassin the landward estuary(meanchlorophylla concentration at stations 26-32) and (top) daily fluctuationsin the tidal currentamplitudeU. Current amplitudes were computed from tidal harmonicconstantsderived from long-

30 /

20•

term

current

measurements

in the central

southbay channel[Chengand Gartner,1985] (predictedU from J. W. Gartner (personal communication,1995)). Biomassincreased rapidly after March 20, when U waslessthan

10 /

0.5 ms -•. .e'

I

'

I

I

'

I

I

I

79-,%b 10-/l•a r 20-1•ar30-graf

19-4,0 r

1985

0.2

0.1



Figure 9. Tidal forcingand phytoplankton bloomsin southSanFranciscoBay. The net rate of biomasschange,R, was calculated from the week-to-week changesin mean chlorophylla concentrationin the landward estuary(stations26-32). The tidal regimeis representedby the 7-day mean current amplitude, U7. Pointsrepresentall the weekly-

/

scale events

of

1978-1995

in which

the

meanchlorophylla concentrationfell below

ß

or increased above10IxgL- • (updated and

'e-

89

revisedfrom Figure 4 of Cloern [1991b]). The dashedline is the least squaresregressionof R againstU7 (r = 0.66). Deviations around this regressionin 1983, 1989, 1990, and 1995 are examinedin Figures10-13 to illustrate

4

ability.

95

0.4

0.5

0.6

MeanCurrent Amplitude U7 (m/s)

0.7

other mechanisms

of bloom vari-

144 ß Cloern' COASTAL PHYTOPLANKTON BLOOMS

34, 2 / REVIEWS OF GEOPHYSICS

Distance (km) 0

10

Distance (km)

20

30

40

O0

10

20

30

g'-.2j:•-•i!::iii•:'"----:-'---•::-.... "-'•il;-":g•i:•iiii?::::1111iii'#.. .";':'::!iii!::i::iii::i::i::::::i:" ".:i...:'•i•;;--':.-i•.. :ii!::i•i::i:;::i•!• ':•

40

!•;•i•.g•.:•i•:;..i.•.....•.•`..i•;•`..•.•.....?.`:.:..•i•i• ':":•½i•i•.' •:"•-"•:•i[?-::.-'i-'•.'.'• :':!•i•i![!•i•i!::i•i.:-::!..:-•:i•::::: !:•.• ....... -...........Mar-MayRiverFlow(m3/s) 5000

10"•:•' ..........

.•i?.'ii•!,?,•i,•. ....

15-

3000

i

i

1000

ii

.:: ........... :?.....-:-

0 1978

1980 20

20>26

24

22

20

18

16

14

1982