ZHANG Xia,ZHANG Jingping,HUANG Xiaoping,and HUANG Liangmin

Key Laboratory of Tropical Marine Bio-resources and Ecology, South China Sea Institute of Oceanology,Chinese Academy of Sciences,Guangzhou 510301,P. R.China

1 Introduction

Estuaries are transition zones linking freshwater and marine systems,and are therefore characterized by gradients of chemical,physical and biological components in the water column.These gradients strongly influence the spatial and temporal distribution,species composition and abundance of phytoplankton in estuaries (Quinlan and Phlips,2007).Hence,the phytoplankton community will be composed of a series of assemblages along the freshwater to marine continuum as a result of physical tolerances and resource competition (Attrill and Rundle,2002).

Phytoplankton dynamics are mainly forced by bottom-up and top-down controls,and these interactions are more complicate in estuaries,because of the seasonal variations of the river discharge besides wind and tidal energy (Pennock and Sharp,1994; Cloern,1996).In many estuaries,salinity appears to be the main factor determining the shifts in many phytoplankton taxa,and other factors (i.e.,temperature,light and nutrient availability) have coupling effects with salinity to explain the spatial and temporal alteration of the phytoplankton composition (Quinlan and Philips,2007).Phytoplankton in the Suwannee Rive estuary progressively changes from a freshwater to marine dominated community along a salinity gradient (Quinlan and Philips,2007).

It is common accepted that salinity and nutrient concentrations are not necessarily independent in estuaries,as both are a function of the riverine discharge.It has been approved that most of dissolved inorganic nitrogen in the PRE is from the freshwater discharges and its concentration declines along the salinity axis resulting from dilution and biological uptake (Linet al.,2004).Hence nitrate shows the conservative behaviour (Huanget al.,2003; Yinet al.,2000).However,there is no linear relationship between the phosphate and salinity,suggesting that the main source of phosphate not depends on the freshwater input.Huanget al.(2003) indicated that land-based sources outside the estuary brought by coastal current and flood tide current should be the main sources of the phosphate in the PRE.

Phosphorus is considered as a frequent limiting nutrient to phytoplankton growth in the PRE and other estuaries (Smith,2006; Yin and Harrison,2008).In addition,low light availability has been proposed as another main reason for lower yields of phytoplankton (Cloern,1987;Fisheret al.,1988; Kocumet al.,2002) and much of the temporal variability of phytoplankton biomass in estuaries is related to variations in light availability.Consequently,some phytoplankton species serve as seed populations downstream and start to bloom when light penetration or other ambient conditions are improved in the PRE (Yinet al.,1995).Huanget al.(2004) observed that distribution of the dominant phytoplankton species varies with salinity of seawater,and their abundances correlate negatively with nutrients and zooplankton.In addition,Yinet al.(2004) sketched a typical spatial variability of phytoplankton biomass in the PRE: low biomass and productivity due to rapid dilution and light limitation in turbid estuarine waters; a regional maximum of biomass and productivity under intermediate salinities in the coastal plume and low biomass due to nutrient limitation in oceanic waters.

Phytoplankton diversity and distribution,seasonal and spatial variations of biomass and bloom events,and their linear correlations with nutrient have been widely investigated in the PRE (Huanget al.,2004; Yinet al.,2000,2001,2004).However,most of these studies only included one or two distinct seasons and little explored the relationship between dominant phytoplankton species and key environmental factors.In this study,the spatial and temporal composition and abundance of phytoplankton were examined in relation to physico-chemical and zooplankton variables in three seasons under distinct hydrodynamic conditions.To achieve these goals a series of exploratory statistical methods have been employed in the analysis of the combined data sets.

2 Materials and Methods

2.1 Study Area

The Pearl River Estuary (PRE),being located in southern Guangdong Province,is formed by inflows of the Pearl River to the South China Sea (SCS) through 8 entrances,4 of which enter the Ling Ding Yang estuary,a sub-estuary of the Pearl River Estuary.Geographically,the Pearl River has a subtropical climate with a long summer and a short winter; the annual average temperature is 21–22℃ and total rainfall is 1600–2000 mm.The annually averaged river discharge is 10524 m3s?1,with 20% of it occurring during the dry season from October to March and 80% during the wet season from April to September (Zhao,1990).In recent decades,the Pearl River Delta area has been under rapid socio-economic changes,including the population growth,industrialization,and urbanization.A large amount of domestic,industrial,and agricultural effluent is discharged into the river system,causing deterioration of the aquatic environment in this estuary (Huanget al.,2003; Qiuet al.,2010).Then the harmful algal blooms (HAB) happened more frequently in recent years in the PRE (Yinet al.,2001).

2.2 Sampling and Analysis

Three surveys were conducted at 16 stations respecttively in July 2006 (high river flow season,summer),March 2007 (low river flow season,spring),and September 2007 (normal river flow season,autumn) in the PRE(Fig.1).The whole study area is within the sub-tidal zone with strong fluvial runoffs and marine water inter-reactions.Temperature and salinity were measured using YSI6600.Integrated water samples were collected with 5 L Niskin bottles through the full depth of the water column.

Fig.1 Sampling stations of the Pearl River Estuary.

Suspended solid content (SSC) samples were acquired through pre-weighed Whatman GF/F fiber filters (25 mm).The filters were dried and weighed to determine the amount in mg L?1of sample.Dissolved oxygen was determined according to Dickson (1995).Total phosphorus(TP) concentration was measured by colorimetry after digestion of the unfiltered samples with K2S2O8to orthophosphate (Ebinaet al.,1983).Total nitrogen (TN) was digested simultaneously with TP.After digestion,TN was measured as nitrate and absorbance was measured at 543 nm.Nitrate (plus nitrite) and ammonium were determined manually following the procedures of Woodet al.(1967)and Slawyk and MacIsaac (1972),respectively.DIN was the sum of nitrate,nitrite and ammonia.DIP concentration was measured based on the methods by Armstronget al.(1967).

A plankton net (77 μm mesh size) was employed in collecting phytoplankton and samples were stored with Lugol’s solution and identified according to Jinet al.(1965).Zooplankton was sampled using a modified WP-2 net (mouth size,0.5 m2; mesh size,505 μm) by towing vertically from 1 m above the bottom to the surface and the filtered water volume was determined from the rope length multiplied by mouth size (Zhang,1993).The samples collected were preserved immediately in 5% formaldehyde for further identification of species under microscope (Chen and Zhang,1974; Chen and Shen,1974).

2.3 Data Analysis

Non-metric multidimensional scaling (NMDS) was designed to distinguish distinct regions within the estuary based on the phytoplankton assemblages.These analyses were undertaken using PRIMER (v.6) software (Plymouth Routines in Multivariate Ecological Research,Plymouth Marine Laboratory,Plymouth,UK; Clarke and Gorley,2006).On the basis of Bray–Curtis similarities,the similarity percentages analysis (SIMPER) was applied to phytoplankton species abundance,in order to allow the separation of every two groups of stations according to different phytoplankton community structure.We intended to analyze the dominant phytoplankton density as dependent factors in a linear regression analysis with water temperature,salinity,SSC,NO3-N,NH4-N,DIN,TN,DIP,TP concentrations and copepod density as the independent factors.Owing to the dominance of the copepods in zooplankton community,we considered copepods as main efficient consumers.Since several of the variables obviously were correlated,we performed a principal component regression (PCR),i.e.,principal component scores from a principal component analysis (PCA) instead of the actual x-variables in the regression analysis (Quinn and Keough,2002).Statistical analysis was performed using SPSS 17.0 for Windows.

3 Results

3.1 Environmental Variables

In our surveys,the average of water temperature was higher in summer (29.3℃) and autumn (27.6℃) than that in spring (19.9℃).The average of salinity was higher in spring (20.4) than those in autumn (10.9) and summer(9.4).In summer,a stratified layer and salinity gradient were observed in the middle reaches of the estuary resulting from large river discharge and salt wedge intrusion (Fig.2a).A well developed salt wedge could extend to station S8.In spring,salinity difference existing between the surface and bottom layer was much smaller.Salinity was roughly less than 25 in the upper reaches(S1?S5) or on the western side (S9,S10 and S12) of the PRE,while it varied from 30 to 35 in the middle and lower reaches (S8,S11 and S13?S16).About 95% of the sediment load is delivered during the wet season (Harrisonet al.,2008).SSC was higher in high flow season(41.65 mg L?1) than in spring (25.34 mg L?1) and exhibited the peak values near the river mouth in both seasons(Fig.2b).

Fig.2 Variations of salinity in surface and bottom layers (a)and SSC (b) in the Pearl River Estuary.

Generally,the distribution patterns of nutrient were similar to that of SSC and the concentration decreased as the riverine water mixed with the shelf water (Fig.3).DIN mainly came from the runoff of Pearl River and the concentration in the northern area was generally higher than that in the southern part.The maximum DIP concentration was observed at coastal stations S4 and S7 (near Shenzhen Bay) in two seasons of 2007 (Fig.3) which were under intense anthropogenic pressure and implicated that land based pollutants near Shenzhen Bay contributed to phosphate load greatly (Huanget al.,2003).Meanwhile,the lowest DIP values in summer,spring and autumn were 0.28 μmol L?1,0.43 μmol L?1and 0.06 μmol L?1,respectively,which were then observed in the outer estuary.

Fig.3 Spatial variation of phosphorus and nitrogen concentrations in three seasons.

3.2 Abundance and Community Structure of Phytoplankton and Zooplankton

Average density of phytoplankton in summer,spring and autumn was 4.0×105,1.8×105and 2.3×105cells L?1,respectively (Table 1).In these three seasons,the higher abundance occurred in the lower section of the PRE,and the maximum value appeared at the station 15,which was situated in the southern waters of Hong Kong.In the upstream and middle PRE (S1 to S10),phytoplankton abundance was several orders of magnitude lower than those detected in the downstream.Phytoplankton abundance of the three seasons was negatively related with SSC,DIN and DIP while positively with salinity (Fig.4).151 species were identified in the three surveys,belonging to 42 genera of Bacillariophyta (105 species),10 genera of Pyrrhophyta (31 species),4 genera of Chrysophyta(8 species) and 5 genera of Cyanobacteria (7 species).Diatoms almost dominated the community in the three cruises except at some stations in the autumn cruise(Fig.5). The most important diatom species (average abundance >5% of the total) wereCoscinodiscus jonesianus,Skeletonema costatum,Thalassiosira subtilis,Biddulphia sinensis,Pseudonitzschia pungens,P.delicatissima,Aulacoseira granulata,A.granulata v.angustissima,Chaetocerospseudocurvisetus,C.curvisetusandC.lorenzianus.

The NMDS ordination of stations generated by phytoplankton abundance is illustrated in Fig.6.Stress values associated with this two-dimensional plot was 0.04,0.14 and 0.07 for the three seasons,respectively,revealing that this representation of stations was sound (Clarke and Gorley,2006).The successive position of the stations revealed algal community differences along the estuary.

In the summer,stations located in the upper estuary (S1 to S7) were classified as group A characterized by the prevalence of freshwater phytoplankton species (A.granulataandA.granulatav.angustissima).S16 located in the outer estuary,formed group C.Group B (S8 to S15)individualized dominance of coastal and estuarine species(S.costatumandP.delicatissima) (Fig.5,Fig.6a).

Table 1 Mean and ranges of the physico-chemical parameters,phytoplankton and zooplankton abundance of three surveys in the Pearl River Estuary

Fig.4 Linear regression of phytoplankton abundance (log10-transformed) and SSC (a),salinity (b),DIN (c),DIP (d)in the PRE.

Fig.5 Spatial variation of dominant phytoplankton species in three studied seasons.

Fig.6 NMDS analyses of the phytoplankton community structure in three studied seasons.

In the spring,Stations S1 to S14 and S16 were categorized into one group (group A),while S15 constituted another one (group B) (Fig.6b).Components of group A were characterized by dominance of marine species (e.g.Thalassiosira subtilisandBiddulphia sinensis) and freshwater species only made a minor contribution in the community: group B was dominated byC.pseudocurvisetusand total phytoplankton abundance was much larger than that in group A (Fig.5).

In the early autumn cruise,16 stations could be separated into three groups according to phytoplankton structure (Fig.6c).Group A consisted of S1 to S10 with an overwhelming speciesMicrocystisspp.Stations in the inner estuary (S11 to S14) were clustered as group B.Outer stations S15 and S16 were labeled as group C withP.pungens,SkeletonemaandChaetocerosbeing dominant species (Fig.5).

The average of zooplankton abundance was highest in summer (778 ind.m?3),followed by those in autumn (734 ind.m?3) and spring (351 ind.m?3).Average copepod densities were 594,294 and 407 ind.m?3in summer,spring and autumn,respectively.The copepod assemblages were dominated byPseudodiaptomus poplesia,Sinocalanus solstitialis,S.laevidactylus,Acarfia spinicaudaandPavocalanus crassirostrisin the three seasons (Fig.7).

Fig.7 Spatial variation of total zooplankton,dominant copepods abundance in three studied seasons.

3.3 Environmental Variables Related with Phytoplankton Community Explored by PCR

The ten abiotic and biotic variables that potentially explained the horizontal heterogeneity of phytoplankton community structure in the summer and spring were interpreted in terms of principal components (Table 2).In both seasons,SSC and nutrient species were positively associated with each other and both of them were negatively associated with salinity.In the summer of 2006,the first principal component (PC1),which explained 60.07%of the variation of the environmental data,was positively correlated with nutrient species (NO3-N,NH4-N,DIN,TN,DIP and TP),SSC and water temperature were negatively correlated with salinity.The second component(PC2),which explained 11.49% of the variation of the environmental data,was correlated with copepod density.In the spring,PC1 explained 57.89% of the variation of the environmental data,was positively correlated with SSC,salinity,nutrient species and copepods.The second component (PC2) was positively related with temperature and explained 19.95% of the variation of the environmental data.The linear regression models including principal components and dominant phytoplankton species are shown in Table 3.In the summer,PC1 was negatively related withS.costatum,P.delicatissimaand positively related withA.granulataandA.granulatav.angustissima.In the spring,T.subtilisandB.sinensiswere negatively related with PC1.

Table 2 Results of principal component analysis (PCA) describing the major environmental gradients of variation in two seasons

Table 3 Regression analysis data of the principal components (PC) as dependent factors for the dominant phytoplankton species in two seasons

4 Discussion

4.1 Key Environmental Factors Regulating Phytoplankton Abundance

In turbid estuaries like the Schelde Estuary,phytoplankton biomass is mainly regulated by concentrations of suspended matter (Kromkamp and Peene,1995;Kocumet al.,2002).Although suspended sediment concentration in the PRE is lower than the Yangtze River estuary and Yellow River estuary,Yin and Harrison(2008) pointed out low chlorophylladuring summer in the PRE mainly resulted from high turbidity due to high freshwater discharge.The 10 mg L?1has been regarded as a threshold concentration of SSC for phytoplankton inhibition (DeMasteret al.,1996; Ragueneauet al.,2002).In our summer and spring surveys,suspended solid content was much higher than this threshold value except at the outer stations S15 and S16.According to Cloern (1987),average euphotic depth in the inner estuary (S1 to S14) in summer was only 1.43 m,while at the outer sites S15 and S16 such average depths were 3.18 and 3.90 m,respectively.Hence,phytoplankton growth was light-limited and phytoplankton could not accumulate in the inner estuary in spite of the high nutrients (Yinet al.,2000; Yinet al.,2004).In the Chesapeake and Delaware Bays,a chlorophyll maximum occurred downstream of the turbidity maximum (Fisheret al.,1988).Turbidity was considered as the main limiting factor for phytoplankton abundance in the PRE and exhibited an inversed relationship.

It has been widely accepted that the Ks (half-saturation constant) is associated with nutrient limitation for algal uptake (Aksnes and Egge,1991).Typical Ks values of nitrate and phosphate for natural phytoplankton populations are ca.1?2 μmol L?1and 0.1?0.5 μmolL?1,respectively (Fisheret al.,1988).As for phosphate,its concentration was always below 0.5 μmol L?1in outer stations,suggesting a potential P-limitation.In the present study,NO3?concentrations at the outer stations S15 and S16 were both higher than 2 μmol L?1in the three seasons,indicating that N was abundant.Yin and Harrison (2008)figured in the region of the chl.amaximum,phosphorus was rapidly depleted by phytoplankton and excess of N was left in the water.In the present study,phytoplankton abundance was negatively related with DIP in the three surveys.The maxima of phytoplankton density occurred in the outer estuary (S15 and S16) where phosphate was reduced to extremely low level.Therefore,phosphate may be another factor regulating phytoplankton abundance in the PRE,which was also suggested by Huanget al.(2003).

4.2 Environmental Factors Shaping Horizontal Phytoplankton Community Structure

Along the estuary axis,phytoplankton composition shifted from dominant freshwater species to salinitytolerant species,and finally to the increasing dominance of estuarine forms (Wong and Townsend,1999).It is obvious that the different dominating species of diatoms at different sites are expected to use distinct strategies in order to adapt to the diverse trophic conditions and mixing levels (Reynoldset al.,2002; Alves-de-souzaet al.,2008).In our case,configuration explored by NMDS gave a good representation of the distance matrix,allowing an easier identification of the dominant species assemblage at the salinity gradient.A.granulataandA.granulatavar.angustissimahad the highest abundance near Humen mouth (S1) with a surface water salinity of 0.9 in the summer,and then decreased with increasing salinity along the estuarine axis.This distribution pattern was very similar to those reported by Huanget al.(2004).A.granulataandA.granulatavar.angustissimaare typical riverine phytoplankton in large rivers of the world(Rojoet al.,1994; Lewiset al.,1995).With a limited salinity tolerance,these two freshwater species were restricted to the head of the estuary in the summer and exhibited a positive relationship with PC1,which reflected the freshwater inflow effects.

On the contrary,other dominant diatoms such asS.costatum,P.delicatissima,T.subtilisandB.sinensiswere negatively associated with PC1,indicating that their development was restricted by turbid freshwater inflow.S.costatumis a pollution-tolerant and indicator species of eutrophication (Nassar,2000).In the competition experiments with two marine diatoms,S.costatumsupplantedPhaeodactylum tricornutumat high N: P ratios with the medium (Pauw and Naessens-Foucquaert,1991).S.costatumis well adapted to higher N: P conditions andhigher N: P ratio is one of a characteristic of the nutrient regime in the PRE.The optimum temperature and salinity ranges forS.costatumwere 20–26℃ and 25–30,respectively (Jin,1965).Chenet al. (2005) found that the growth ofS.costatumwas little affected by a salinity shift from 18 to 35.7.Another cosmopolitan genus,Pseudonitzschiais well known for its adaption to great ranges in salinity,15–40 for optimal growth in some species(Thessenet al.,2005).In spring,eurythermal speciesT.subtilisandB.sinensisalways appear in coastal waters and estuaries which are generally characterized by high nutrient concentrations (Qasimet al.,1973).Sudden changes in salinity (15?30) and temperature (12?20℃)almost did not influence growth ofB.sinensis(Jahnke and Baumann,1983).In the present study,two independent variables,salinity and turbidity,played important roles in the distribution of the major phytoplankton species.

The hydrographic conditions of the PRE and other estuaries (e.g.,Swan,Yura and Pearl Estuary) drastically change between wet and dry periods (Thompson,1998;Harrisonet al.,2008).In the spring survey,probably because of the decrease of the freshwater flushes from the Pearl River,the inner estuary was relatively homogenous.Hence,sampling stations of the inner estuary were classified as one group.Some estuarine and marine species were detected at the head of the estuary.

In the autumn,an enhanced proliferation of potentially toxicMicrocystiswas detected at the lower section of the PRE.Microcystisblooms have often been associated with varying concentrations of the microcystin in the surrounding water.However,during our study,Microcystisabundance was below a blooming level (105cells L?1),which probably imposed limited impairment to environment.The higher temperatures along with water column stability may have supported cyanobacteria,as have been observed in several estuaries and coastal waters,including the Chesapeake Bay (Marshall and Nesius,1996),the Florida Bay (Phlipset al.,1999),the Neuse River estuary(Valdes-Weaveret al.,2006),the San Francisco Bay(Ninget al.,2000),and the York River estuary (Sinet al.,2000).Cyanobacteria abundance in these ecosystems varies considerably,but in general warmer estuaries tend to be associated with higher abundance of cyanobacteria(Murrel and Lores,2004).During the dominance ofM.aeruginosa,water temperatures were from 24.7 to 33.9℃in the eutrophic Furuike Pond,Japan (Imaiet al.,2009).Furthermore,Christianet al.(1986) found the optimum water temperature for this species to be about 28℃.Our investigation is consistent with these results.

The growth ofM.aeruginosawas restricted in the brackish region and negatively impacted by the salinity in Patos Lagoon (Yuneset al.,1997).In addition,M.aeruginosagrew optimally at salinities up to 4,above which growth rate declined to zero at salinity about 25 in laboratory experiments (Robson and Hamilton,2003).In the present study,the salinity was beyond its tolerance threshold in the middle of the estuary.It is probable that colonies ofMicrocystisdeveloped in this area have been seeded from upstream and transported by river discharge and the wind-induced advection,which has also been observed in San Francisco Bay Estuary (Lehmanet al.,2005) and a Western Australia estuary (Robson and Hamilton,2003).

The absence of significant relationship between phytoplankton and zooplankton density was observed in the present study.Low levels of phytoplankton in the inner estuary was probably not related with top-down force as zooplankton (copepods) abundance was very low and restricted only at given regions of the estuary (Liet al.,2004).Aggregations of copepods and other zooplankton tended to occur at boundaries or fronts in the vertical or horizontal planes (Petipa,1985).Tanet al.(2004) observed that a high value of Chl.ausually coincided with a low ingestion rate in the PRE.Hence,phytoplankton abundance was more likely controlled by abiotic factors than by biotic mechanisms in this estuary.

Acknowledgements

This study was supported by the National Natural Sciences Foundation of China (Nos.31000185,41076069),the Strategic Priority Research Program of the Chinese Academy of Sciences (Grant No.XDA11020205).

Aksnes,D.L.,and Egge,J.K.,1991.A theoretical model for nutrient uptake in phytoplankton.Marine Ecology Progress Series,70: 65-72.

Alves-de-Souza,C.,Gonzalez,M.T.,and Iriarte,J.L.,2008. Functional groups in marine phytoplankton assemblages dominated by diatoms in fjords of southern Chile.Journal of Plankton Research,30: 1233-1243.

Armstrong,F.A.,Stearns,C.R.,and Strickland,J.D.,1967.The measurement of upwelling and subsequent biological processes by means of the Technicon AutoAnalyzer and associated equipment.Deep-Sea Research,14: 381-389.

Attrill,M.J.,and Rundle,S.D.,2002.Ecotone or ecocline:Ecological boundaries in estuaries.Estuarine,Coastal,and Shelf Science,55: 929-936.

Chen,Q.C.,and Zhang,S.Z.,1974.The pelagic copepods of the South China Sea.I.Studia Marina Sinica,9: 101-116 (in Chinese with English abstract).

Chen,Q.C.,and Shen,J.R.,1974.The pelagic copepods of the South China Sea.II.Studia Marina Sinica,9: 125-137 (in Chinese with English abstract).

Chen,Y.T.,1995.Sedimentation divisions of the Pearl River mouth.Acta Scientiarum Naturalium Universitatis Sunyatseni,34 (3): 109-114 (in Chinese with English abstract).

Chen,B.Z.,Wang,Z.L.,Zhu,M.Y.,and Li,R.X.,2005.Effects of temperature and salinity on growth ofProrocentrum dentatumand comparisons between growths ofProrocentrum dentatumandSkeletonema costatum.Advances in Marine Science,23: 60-64 (in Chinese with English abstract).

Christian,R.R.,Bryant,W.L.,and Jr.Stanley,D.W.,1986.The relationship between river flow andMicrocystis aeruginosablooms in the Neuse River,North Carolina.Water Resources Research Institute Report223.North Carolina State University.

Clarke,K.R.,and Gorley,R.N.,2006.Primer v6: User Manual/Tutorial,Primer E: Plymouth.Plymouth Marine Laboratory,Plymouth,UK.

Cloern,J.E.,1987.Turbidity as a control on phytoplankton biomass and productivity in estuaries.Continental Shelf Research,7: 1367-1381.

Cloern,J.E.,1996.Phytoplankton bloom dynamics in coastal ecosystems: A review with some general lessons from sustained investigation of San Francisco Bay,California.Review of Geophysics,34: 127-168.

DeMaster,D.J.,Smith,W.O.,Nelson,D.M.,and Aller,J.Y.,1996.Biogeochemical processes in Amazon shelf waters:Chemical distributions and uptake rates of silicon,carbon and nitrogen.Continental Shelf Research,16: 617-643.

Dickson,A.G.,1995.Determination of dissolved oxygen in sea water by Winkler titration.WOCE Operations Manual.Part 3.1.3 Operations & Methods,WHP Office Report WHPO 91-1.

Ebina,J.,Tsuyoshi,T.,and Shirai,T.,1983.Simultaneous determination of total nitrogen and total phosphorus in water using peroxodisulfate oxidation.Water Research,17: 1721-1726.

Fisher,T.R.,Harding,L.W.,Stanley,D.W.,and Ward,L.G.,1988.Phytoplankton,nutrients,and turbidity in the Chesapeake,Delaware,and Hudson Estuaries.Estuarine,Coastal and Shelf Science,27: 61-93.

Harrison,P.J.,Yin,K.,Lee,J.H.W.,Gan,J.,and Liu,H.,2008.Physical-biological coupling in the Pearl River Estuary.Continental Shelf Research,28: 1405-1415.

Huang,L.M.,Jian,W.J.,Song,X.Y.,Huang,X.P.,Liu,S.,Qian,P.Y.,Yin,K.D.,and Wu,M.,2004.Species diversity and distribution for phytoplankton of the Pearl River Estuary during rainy and dry seasons.Marine Pollution Bulletin,49:588-596.

Huang,X.P.,Huang,L.M.,and Yue,W.Z.,2003.The characteristics and sources of eutrophication in Pearl River Estuary,South China.Marine Pollution Bulletin,47: 30-36.

Imai,H.,Chang,K-H.,Kusaba,M.,and Nakano,S.,2009.Temperature-dependent dominance ofMicrocystis(Cyanophyceae)species:M.aeruginosaandM.wesenbergii.Journal of Plankton Research,31: 171-178.

Jahnke,J.,and Baumann,M.,1983.Chemical and physical effects on the shape and growth of the diatomBiddulphia sinensisGreville in batch cultures: A contribution to bioindication in plankton ecology.Aquatic Ecology,17: 5-20.

Jin,D.,1965.The effects of temperature and salinity on reproduction of three phyto-diatoms.Oceanologia et Limnologia Sinica,7: 273-284.

Jin,D.X.,Chen,J.H.,and Huang,K.G.,1965.Marine Planktonic Diatoms in China Sea.Shanghai Scientific and Technical Publishers,Shanghai,China,1-44.

Kocum,E.,Underwood,G.J.,and Nedwell,D.B.,2002.Simultaneous measurement of phytoplanktonic primary production,nutrient and light availability along a turbid,eutrophic UK east coast estuary (the Colne Estuary).Marine Ecology Progress Series,231: 1-12.

Kromkamp,J.,and Peene,J.,1995.On the possibility of net primary production in the turbid Schelde estuary (SW Netherlands).Marine Ecology Progress Series,121: 249-259.

Lehman,P.W.,Boyer,G.,Hall,C.,Waller,S.,and Gehrts,K.,2005.Distribution and toxicity of a new colonial Microcystis aeruginosa bloom in the San Francisco Bay Estuary,California.Hydrobiologia,541: 87-99.

Lewis,W.M.,Hamilton,S.K.,and Saunders,J.F.,1995.Rivers of Northern South America.In:River and Stream Ecosystems.Cushing,C.E.,Cummins,K.,W.,and Minshall G.W.,ed., Elsevier,Amsterdam,219-256.

Li,K.Z.,Yin,J.Q.,Huang,L.M.,and Tan,Y.H.,2004.Spatial and temporal variations of mesozooplankton in the Pearl River Estuary,China.Estuarine,Coastal and Shelf Science,67: 543-552.

Lin.Y.,Su,J.,Hu,C.,Zhang,M.,Li,Y.,Guan,W.,and Chen,J.,2004.N and P in waters of the Zhujiang River Estuary in summer.Acta Oceanologica Sinica,26: 63-73 (in Chinese with English abstract).

Marshall,H.G.,and Nesius,K.K.,1996.Phytoplankton composition in relation to primary production in Chesapeake Bay.Marine Biology,125: 611-617.

Murrell,M.C.,and Lores,E.M.,2004.Phytoplankton and zooplankton seasonal dynamics in a subtropical estuary: Importance of cyanobacteria.Journal of Plankton Research,26:371-382.

Nassar,M.Z.,2000.Ecophysiological studies on phytoplankton along the Western coast of Suez Gulf.PhD thesis,Faculty of Science,Tanta University,1-22.

Ning,X.,Cloern,J.E.,and Cole,B.E.,2000.Spatial and temporal variability of picocyanobacteriaSynechococcussp.in San Francisco Bay.Limnology and Oceanography, 45: 695-702.

Pauw,N.D.,and Naessens-Foucquaert,E.,1991.Nutrient- induced competition between two species of marine diatoms.Aquatic Ecology,25: 23-27.

Pennock,J.R.,and Sharp,J.H.,1994.Temporal alternation between light-and nutrient-limitation of phytoplankton production in a coastal plain estuary.Marine Ecology Progress Series,111: 275-288.

Petipa,T.S.,1985.Production and concentration of plankton at boundary water masses: Perspectives of investigations.In:Proceedings of the 19th European Marine Biology Symposium.Gibbs,P.E.,ed.,Plymouth,Devon,UK,61-71.

Phlips,E.J.,Badylak,S.,and Lynch,T.C.,1999.Blooms of the picoplanktonic cyanobacteriumSynechococcusin Florida Bay,a subtropical inner-shelf lagoon.Limnology and Oceanography,44: 1166-1175.

Qasim,S.Z.,Bhattathiri,P.M.,and Devassy,V.P.,1973.Growth kinetics and nutrients requirements of two tropical marine phytoplankters.Marine Biology,21: 299-304.

Qiu,D.J.,Huang,L.M.,Zhang,J.L.,and Lin,S.J.,2010.Phytoplankton dynamics in and near the highly eutrophic Pearl River Estuary,South China Sea.Continental Shelf Research,30: 177-186.

Quinn,G.P.,and Keough,M.J.,2002.Experimental Design and Data Analysis for Biologists.Cambridge University Press,Cambridge,1-38.

Quinlan,E.L.,and Philips,E.J.,2007.Phytoplankton assemblages across the marine to lowsalinity transition zone in a blackwater dominated estuary.Journal of Plankton Research,29: 401-416.

Ragueneau,O.,Lancelot,C.,Egorov,V.,Vervlimmmeren,J.,Cociasu,A.,Déliat,G.,Krastev,A.,Daoud,N.,Rousseau,V.,Popovitchev,V.,Brion,N.,Popa,L.,and Cauwet,G.,2002.Biogeochemical transformations of inorganic nutrients in the mixing zone between the Danube River and the North-western Black Sea.Estuarine,Coastal and Shelf Science,54:321-336.

Reynolds,C.S.,Huszar,V.,Kruk,C.,Naselli-Flores,L.,and Melo,S.,2002.Towards a functional classification of the freshwater phytoplankton.Journal of Plankton Research,24:417-428.

Robson,B.J.,and Hamilton,D.P.,2003.Summer flow event induces a cyanobacterial bloom in a seasonal Western Australia Estuary.Marine and Freshwater Research,54: 139-151.

Rojo,C.,Alvarez-Cobelas,M.,and Arauzo,M.,1994.An elementary,structural analysis of river phytoplankton.Hydrobiologia,289: 43-55.

Sin,Y.,Wetzel,R.L.,and Anderson,I.C.,2000.Seasonal variations of size-fractionated phytoplankton along the salinity gradient in the York River Estuary,Virginia (USA).Journal of Plankton Research,22: 1945-1960.

Slawyk,G.,and Macisaa,J.J.,1972.Comparison of two automated ammonium methods in a region of coastal upwelling.Deep-Sea Research,19: 521-524.

Smith,V.H.,2006.Responses of estuarine and coastal marine phytoplankton to nitrogen and phosphorus enrichment.Lim-nology and Oceanography,51 (1): 377-384.

Tan,Y.,Huang,L.,Chen,Q.,and Huang,X.,2004.Seasonal variation in zooplankton composition and grazing impact on phytoplankton standing stock in the Pearl River Estuary,China.Continental Shelf Research,24 (16): 1949-1968.

Thessen,A.E.,Dortch,Q.,Parsons,M.L.,and Morrison,W.,2005.Effect of salinity onPseudonitzschiaspecies (Bacillariophyceae) growth and distribution.Journal of Phycology,41: 21-29.

Thompson,P.A.,1998.Spatial and temporal patterns of factors influencing phytoplankton in a salt wedge estuary,the Swan River,Western Australia.Estuaries,21: 801-817.

Valdes-Weaver,L.M.,Piehler,M.F.,Pinckney,J.L.,Howe,K.E.,Rossignol,K.,and Paerl,H.W.,2006.Long-term temporal and spatial trends in phytoplankton biomass and classlevel composition in the hydrologically variable Neuse-Pimlico estuarine continuum,North Carolina,USA.Limnology and Oceanography, 51 (3): 1410-1420.

Wong,M.,and Townsend,D.,1999.Phytoplankton and hydrography of the Kennebec Estuary,Maine.Marine Ecology Progress Series,178: 133-144.

Wood,E.D.,Armstrong,F.A.J.,and Richards,F.A.,1967.Determination of nitrate in seawater by cadmium-copper reduction to nitrate.Journal of the Marine Biological Association of the United Kingdom,47: 23-31.

Yin,K.,and Harrison,P.J.,2008.Nitrogen over enrichment in subtropical Pearl River Estuarine coastal waters: Possible causes and consequences.Continental Shelf Research,28:1435-1442.

Yin,K.,Harrison,P.J.,Pond,S.,and Beamish,R.J.,1995.Entrainment of nitrate in the Pearl River plume and its biological implications.III.Effects of winds.Estuarine,Coastal and Shelf Science,40: 545-558.

Yin,K.,Qian,P.Y.,Chen,J.C.,Hsieh,D.P.,and Harrison,P.J.,2000.Dynamics of nutrients and phytoplankton biomass in the Pearl River Estuary and adjacent waters of Hong Kong during summer: Preliminary evidence for phosphorus and silicon limitation.Marine Ecology Progress Series,194: 295-305.

Yin,K.,Qian,P.Y.,Wu,M.C.S.,Chen,J.C.,Huang,L.M.,Song,X.Y.,and Jian,W.J.,2001.Shift from P to N limitation of phytoplankton growth across the Pearl River Estuarine plume during summer.Marine Ecology Progress Series,221:17-28.

Yin,K.,Qian,P.Y.,Chen,J.F.,Huang,L.,Zhang,J.,and Wu,M.,2004.Effect of wind events on phytoplankton blooms in the Pearl River Estuary during summer.Continental Shelf Research,24 (16): 1909-1923.

Yunes,J.S.,Mathinsen,A.,Parise,M.,Salomon,P.S.,Ragget,S.L.,Beattie,K.A.,and Codd,G.A.,1997.Microcystis aeruginosagrowth stages and the occurrence of microcystins in Patos Lagoon,Southern Brazil.In:Harmful Algae.VIII.International Conference.Reguera,B.,Blanco,J.,Fernandez,M.L.,and Wyatt T.,eds., Xunta De Galicia,Vigo,18- 21.

Zhang,S.Z.,1993.Distribution characteristics and resources exploitation of zooplankton in the Pearl River Estuary.South China Sea Fishery Research,4: 2-10 (in Chinese with English abstract).

Zhao,H.,1990.Evolution of the Pearl River Estuary.Ocean Press,Beijing,1-27.