Changes in phytoplankton composition in response to tides , wind-induced mixing conditions , and freshwater outflows in an urbanised estuarine complex

Recent reports have shown an increase in potentially harmful phytoplankton in Santos bay (Southeastern Brazilian Coast), located in a highly urbanised estuarine complex. Prediction of blooms is, thus, essential but the phytoplankton community structure in very dynamic regions is difficult to determine. In the present work, we discriminate bloom forming microphytoplankton dominance and their relationship to physical and meteorological variables to look for patterns observed in different tides and seasons. Comparing 8 distinct situations, we found five scenarios of dominance that could be related to winds, tides and rainfall: i) Surfers, diatoms occurring during high surf zone energies; ii) Sinkers, represented by larger celled diatoms during spring tide, after periods of high precipitation rates; iii) Opportunistic mixers, composed of chain forming diatoms with small or elongate cells occurring during neap tides; iv) Local mixers, microplanktonic diatoms and dinoflagellates which occurred throughout the 298 sampling stations; and v) Mixotrophic dinoflagellates, after intense estuarine discharges. Results suggest alterations in the temporal patterns for some bloom-forming species, while others appeared in abundances above safe limits for public health. This approach can also illustrate possible impacts of changes in freshwater discharge in highly urbanised estuaries.


Introduction
Understanding how changes in environmental conditions are capable of driving phytoplankton dominance is key for comprehending coastal dynamics and to develop environmental monitoring programmes.Changes in phytoplankton taxonomic composition, in a given time and place, are well documented and are manifestations of the adaptive strategies and modes of nutrition of occurring species (e.g.Margalef, 1997).In shallow estuaries, short-term and transient gradients of river outflows (Acha et al., 2008), tide regimes, vertical mixing and turbidity (Baek et al. 2008), and nutrient inputs (Heil et al. 2005;Glibert et al., 2006), especially in areas with intense human occupation, hamper the attempts for structure determination of phytoplankton communities, either taxonomic or as functional groups.Human induced eutrophisation can trigger harmful algal blooms, so the development of conceptual models for phytoplankton dominance is imperative (Cloern, 1996).Furthermore, environmental policies demand quantitative assessments of variables altered by increasing urbanisation and their potential relationships to bloom-forming phytoplankton growth conditions (Whitea et al., 2004;Burkholder et al., 2006).
The relationships between increase in both input and residence times of anthropogenic nutrients in estuarine systems and enhancements of phytoplankton biomass are well documented (e.g., Pinckney et al., 2001).Studies in costal regions, influenced by estuaries, show that inputs of nutrients drive changes on the frequency of local algal blooms (Glibert et al., 2006), with occasional changes in taxa that may even include potential harmful species (harmful algal blooms-HABs - Hallegraeff, 2004;Zingone and Enevoldsen, 2000;Heil et al., 2005), but in many regions of the coastal oceans, the available information is not enough to link the occurrence of HABs and anthropogenic effects (Anderson, 2008).
Highly urbanised coastal areas are influenced by domestic and industrial sewage, often mixed with pollutants and a range of dissolved and particulate materials and the effects of these mixtures on phytoplankton composition are still little understood (Thomas et al., 1980;Walter et al., 2002).Anthropogenic activities in coastal areas can alter the overall freshwater inputs to estuaries (Paerl et al., 1998) as well, which may drive additional changes of phytoplankton community structure (Kimmerer, 2002;Hallegraeff, 2004).
The Estuarine complex of Santos (SE Brazil) (Figure 1) receives pollutants from petrochemical and fertiliser plants (Braga et al., 2000) and high concentrations of nutrients and organic matter (Moser et al., 2004;Aguiar and Braga, 2007).The literature shows high ammonium, phosphate and chlorophyll concentrations in the Santos estuarine complex year round (Gianesella et al., 2000;Aguiar and Braga, 2007), and the occurrence of local summer Skeletonema cf.costatum blooms have been associated with increases in light availability resulting from enhanced vertical stability (Moser et al., 2005).During the winter months, passages of atmospheric cold frontal systems increase wave energy in Santos bay and adjacent sand beaches, leading to blooms of Asterionellopsis glacialis to occur (Villac and Noronha, 2008).To date, there are no conceptual relationships developed among freshwater inputs, tides and wave action on phytoplankton species composition for the Santos estuarine complex area.
Important changes, driven by human impacts in the past decade, are documented for Santos bay and adjacent estuaries, which include a growing volume of domestic sewage disposal and continuous dredging of bottom sediments from the Santos estuarine channel (Schmiegelow et al., 2008).Recent phytoplankton taxonomic surveys report increases in both occurrence and abundance of potentially harmful and mixotrophs dinoflagellates such as Dinophysis acuminata and a number of species from the Prorocentrum genera (Villac et al., 2008).
In this study we discuss the relationships of bloomforming microphytoplankton abundances in an urbanised estuarine complex (Santos bay, Brazil), systematically observed during distinct tides, wind-induced mixing conditions, and freshwater outflow episodes, in the summer and winter months of 2006 and 2007.Our goals were to develop a preliminary conceptual model of phytoplankton dominance in a highly impacted subtropical estuarine complex and to identify scenarios for potential growth of bloom forming microphytoplankton species.

Material and Methods
Santos bay has an average depth of 9.3 m and an area of 44.45 km 2 , and is semi-open to the Atlantic Ocean.Gradients of water temperature between winter and summer are in the order of 4.5 °C.Precipitation rates are significantly higher in the summer, and wind velocities are higher during winter (Harari et al., 2008).Tides are semi-diurnal and are of less than a metre in amplitude, but together with remote (mesoscale) winds, control the dynamics of local water masses and vertical mixing (Alfredini et al., 2008).
The sampling design used in the present study consisted of eight surveys in Santos bay, performed during March 2006and March 2007, and during August 2006and September 2007, to represent the end of summer and winter, respectively.In each selected month, two 10-hour surveys were conducted during spring (ST) and neap tides (NT) (as shown in Table 1), sampling a grid of 16 stations (Figure 1).Two vessels were used simultaneously to minimise sampling time.
Mesoscale wind speeds and current velocities were modelled hourly for each month of observation, using the National Centre for Environmental Prediction/National Centre for Atmospheric Research model (NCEP/NCAR, Kalnay et al., 1996).In situ precipitation rates were obtained from a gauge station, located at 23°57' S and 46°18' W. Wind and precipitation rates model outputs and in situ data, respectively, were averaged for a period of 5 days (four prior the samplings) to minimise the influence of short-term events.Wind direction is presented as N-S and E-W components.
Vertical profiles of temperature and salinity were obtained for stations located above the 10 m isobath with a CTD (Seabird, Seacat).Water transparency was estimated with a Secchi disk.Water samples were collected with clean plastic buckets at the surface, and at 1-2 m above the bottom using horizontal van Dorn bottles.Samples were analysed for salinity, total and organic particulate matter (PM and POM) concentrations, chlorophyll-a concentration, and microphytoplankton composition and enumeration.Discrete surface and near-bottom salinity samples were measured with an AUTOSAL salinometer (Beckman RS10).The instantaneous degree of vertical mixing was estimated by differences in salinity measured near bottom and at the surface (∆S).
From 250 to 500 mL were filtered, immediately after sampling, onto pre-dried and weighted Whatman GF/F filters, for quantification of particulate matter (mg.L -3 ).Filters were kept in the dark in tightly sealed containers containing silica gel pellets until laboratory analysis for total (PM) and organic suspended particulate matter (POM) by the gravimetric method (APHA, 1985).
A volume of 250 mL of sample was concentrated on GF/F filters, and stored in liquid nitrogen immediately after filtration for quantification of chlorophyll-a concentration (chl, mg.m -3 ).The filters were extracted at -10 °C, for at least 24 hours, in pre-cooled 90% acetone: DMSO solution (6:4 by volume) (Shoaf and Lium, 1976).Chlorophyll was measured in a routinely calibrated Turner Designs Trilogy fluorometer, following the method proposed by Welschmeyer (1994).
Phytoplankton composition were analysed only in surface samples, from aliquots of 100 ml, preserved with buffered formaldehyde (0.4 % final concentration) and stored at room temperature in the dark.Sub samples of 10 to 25 mL were placed in sedimentation chambers (Utermöhl, 1958) and phytoplankton was identified and enumerated Chlorophyll-a CHL (mg.m -3 ); particulate matter PM (mg.L -1 ); particulate organic matter POM (% of total PM); coloured disolved organic matter CDOM (absorption); Secchi disk depth SECCHI (m); winds WIND (m.s -1 ), precipitation rates PPT (mm.day -1 ); and densities (cell.L -1 ) of abundant species: Skeletonema cf.under an inverted light microscope (Zeiss, Axiovert) at 400× magnification.At least 200 cells were counted for microphytoplankton (Lund et al., 1958), identified according to Tomas (1997) and Tenenbaum et al. (2004).Species were considered abundant when their observed numbers were above the mean of a particular cruise in at least 55% of the stations (Lobo and Leighton, 1986), in this case in 9 of the 16 stations.
For the statistical analyses, density of abundant phytoplanktonic organisms were grouped according to i) tidal conditions: neap and spring tides; ii) time of the year: summer and winter seasons; and iii) year of survey: 2006 and 2007.One-way ANOVA tests were applied to log-transformed cell density data (log 10 x + 1), after a priori homogeneity test.Results from ANOVA were tested a posteriori with a Tukey test (Zar, 1999).Spearman correlation tests were applied to the entire data set (Table 1) to detect significant correlations (p < 0.05) between abundant species and environmental data.A Canonical Correspondence Analysis (CCA) was applied to verify the relationships between the density of abundant microphytoplanktonic species and the complementary observed variables, after log-transforming the data (Table 1).A Monte Carlo permutation test, followed by a stepwise forward selection was carried out to determine the significance of the correlations between the environmental and species data.The analyses were carried out with the CANOCO software package, version 4.0.

Results
In most surveys, weak to moderate northern winds (around 2.5 m.s -1 ) prevailed during the 5-day window preceding the samplings.Southern winds were observed in March 2006 NT and in August 2006 ST surveys (Figure 2), and both September 2007 surveys, moderate easterly winds were unusually persistent.Hereafter, the individual surveys will be identified as their respective months and years, followed by NT or ST to discriminate between neap and ST.
Precipitation rates were low during the March surveys (mean values bellow 6 mm.day -1 , see Table 2), except for the ST survey of 2007, when values were closer to the summer climatology (about 30 mm.day -1 ) (DAEE, 2009).In all winter surveys (August 2006 and September 2007), average precipitation rates were also lower than climatological values, 15 mm.day -1 (DAEE, 2009).
Salinities were generally higher than 30 except during September 2007, when values were below 28 (Figure 3); salinities were also low in near bottom waters during September 2007-NT.Higher salinity was observed in March 2006-ST, and a clear vertical salinity gradient (∆S) was observed during both March 2007 surveys.Salinity was not directly related to precipitation rates, with high values observed following rainy events, in March 2006 and 2007 (as shown in Table 2 and Figure 3), but were negatively correlated to winds from EW quadrant, and positively correlated to the vertical salinity gradient (∆S) (Table 3).
Secchi disk depths were around 2 m for all surveys, with slightly more transparent waters observed during March 2007-NT, when salinities varied from 31 to 32 (Figure 3).The highest total particulate matter (PM) of 109.3 mg.L -1 was detected in August 2006 ST.Low PM values were observed in all NT surveys (Figure 3).In both September 2007 observations, PM were about one third of the remaining surveys, followed by an increase of the relative contribution of POM, notably during that of the NT.Percentage of POM were positively correlated to southern winds, chlorophyll-a and suspended matter concentrations.Chlorophyll concentrations varied from 1.1 mg.m -3 (March 2006 NT) and 37.9 mg.m -3 (August 2006-NT), and were positively correlated to salinity and PM.
One-way ANOVA tests applied to the densities of abundant microphytoplankton and using as discrimant factors: tide (Neap × Spring); seasons (summer × winter), and year (2006 × 2007); showed significant differences (p < 0.05) between years (number of samples = 126; F = 11.8;p = 0.0008) and to a lesser extent between tides (number of samples = 126; F = 13.79;p = 0.0003), showing no seasonal significant variability.These results were confirmed with a posteriori Tukey tests.
Both spearman correlations (Table 3) together with CCA analysis (Figure 4) suggested that five groups with distinct microphytoplanktonic dominance, as well as chlorophyll-a and PM and POM values, could be related to salinity, tides, precipitation rates, and wind direction and speed (Tables 4 and 5; Figure 5).Group 4, composed by mixotrophic dinoglagellates (Dinophysis acuminata and Gymnodinium species), a flagellate (Chattonella cf.antiqua) and a pennate diatom (Navicula sp.), was abundant in 2 surveys, March 2006, NT and ST, probably due to ressuspension and associated with a higher POM contribution to PM, since the higher densities were observed in stations 1, 4, 5 and 8, near the estuarine channel mouth.Dinophysis acuminata, treated individually, occurred in 3 situations, March 2006 NT, and September 2007, spring and NTs, with higher densities in the last surveys, associated with E winds and estuarine discharge (lower salinity values).Dinophysis acuminata, Chattonella cf.antiqua and Gymnodinium species were positively correlated to East and West winds.The lowest salinities in both September 2007 surveys, accompanied an overall decrease in microphytoplankton densities and chlorophyll-a concentrations, as well (as shown in Tables 2  and 4 and also in Figure 3), while POM contribution to PM were the highest of all surveys.Precipitation rates previous to September 2007 observations were low, but eastern winds favour estuarine discharge.Group 5, Skeletonema cf.costatum, Prorocentrum minimum, Prorocentrum gracile and Scrippsiella trochoidea, local mixers, was abundant in all surveys, and higher densities were attained in March 2006, NT conditions, high salinity and strong North winds.59) 11 ( 16) 38 ( 29) 41 ( 22) 29 ( 25) 30 ( 35

Patterns for species composition
Santos bay has an average depth of only 9.5 m, thus wind and tide-driven currents force local transport of dissolved and particulate materials (Geyer, 1997;Durand et al., 2002).Due to the south orientation of Santos bay opening, shifts in wind direction will influence turbulence levels, circulation patterns, and modify estuarine outflows.Our observations showed that changes in chlorophyll-a concentration followed closely the abundances of microphytoplankton species, especially during NTs at the end of summer, always with dominance of diatoms.Previously described mechanisms linking vertical mixing and consequent light availability to local phytoplankton abundance were not always repeated (e.g.Moser et al., 2004;Moser et al., 2005).Nonetheless, the observed contribution of diatoms, dinoflagellates and other phytoplanktonic groups (silicoflagellates, ciliophora, microphytoplanktonic flagellates, and cyanobacteria) to the overall phytoplankton community richness were similar to world coastal seas (Sournia et al., 1991), and systems with large inputs of nutrients (Muylaert et al., 2009).Responses of plankton communities to short-term physical fluctuations allow many species to coexist, but even in a rapidly changing environment, a specific phytoplankton group can dominate.Diatoms have a number of physiological adaptations for local vertical mixing (Cullen and McIntyre, 1998).The dominant diatoms found in our study generally support high metabolism, function of favourable surface-volume ratios (Nogueira and Figueiras, 2005) resulting from small-sized cells and chain-forming (Skeletonema cf.costatum) or elongated cells (Guinardia and Pseudo-nitzschia species).The elongated penate chain forming diatom species observed in Santos Bay may belong to the Pseudo-nitzschia seriata complex (Hasle, 1965).Higher densities were observed during 2006 surveys, varying from 10 3 cell.L -1 to 10 5 cell.L -1 , close to the upper limit considered a risk (Mafra-Júnior et al., 2006).
A larger contribution of dinoflagellates to total community (29%) than previous studies in our study region (less than 10%) raises concerns towards possible effects of local eutrophication (Aguiar and Braga, 2007).In addition, routine-dredging activities in the Santos estuarine channel can promote dinoflagellates cists suspension (Pitcher and Joyce, 2009).The Santos estuarine channel has been dredged continuously since 2002, and in average, 300000 m 3 of sediments are removed every month, and disposed in the adjacent inner continental shelf (ANTAQ, 2009).Local dinoflagellate abundances have coincidentally increased since 2002 (Villac et al., 2008).Prorocentrum minimum is a neritic bloom-forming dinoflagellate, responsible for harmful episodes in many estuarine and coastal environments that have been subjected to significant geographical expansion (Taylor et al., 2003;Heil et al., 2005).Prorocentrum minimum occurred in all surveys in densities above 10 4 cell.L -1 , similar to other eutrophised sites, such as Paranaguá bay (Mafra-Júnior et al., 2006) and Santo André Lagoon (Macedo et al., 2001), where densities reached 2.5 10 7 cell.L -1 during blooms.
As in other eutrophic costal regions, Prorocentrum minimum high abundances seem to co-occur with, or develop after, high abundances of Skeletonema cf.costatum (Heil et al., 2005).The processes governing the temporal relationships between these two species are still unknown (Heil et al., 2005).Skeletonema cf.costatum is in fact a composite of distinct species (Sarno et al., 2005;Jinfeng et al., 2008).Thus further studies regarding their co-occurrence are still required, as well as any link to the presence of these species and changes in local eutrophication (Macedo et al., 2001).The lack of seasonal trends observed, in part result from our sample design, was contrary to what is reported in earlier studies.Skeletonema cf.costatum blooms are rather responses to episodic events than characteristic summer-bloom species.Skeletonema cf.costatum was negatively correlated to water transparency, following the expected adaptation to low light conditions (e.g.Sakshaug and Andersen, 1986).Schmiegelow et al. (2008) reported nano-sized flagellates as dominant phytoplankton groups in Santos bay year-round and that shifts to microphytoplanktonic-dominance occur during decreases in mixing and increases in precipitation rates.Neither Macedo et al. (2001) nor Heil et al. (2005) have considered the role of wind intensity and direction that will probably be more decisive establishing light level conditions in a time scale of days to weeks.

Microphytoplankton assemblages and the role of weather-induced variability
Our results show at least five distinct and episodic scenarios of microphytoplankton dominance as responses to winds, tides and rainfall.Northern winds prevailed during the observations, but a number of atmospheric cold frontal systems previous to the August 2006 ST survey, changed winds to the southern quadrant.Increases in Asterionellopsis glacialis abundance, a diatom commonly observed in high-energy surf zones (McLachlan and Brown, 2006), followed the increase in strength of the south wind component.High abundances of A. glacialis in the area, is perhaps a result of auxospores suspension or physical accumulation of vegetative cells (Villac and Noronha, 2008).In addition to watercolour changes, the large amount of mucilage released by these blooms can become vehicles for pollutants (Koukal et al., 2007), thus their occurrence and consequences must be better understood in urbanised areas.
During September 2007 eastern winds atypically persisted for many days prior to both sampling dates, during which we observed the lowest chlorophyll-a concentrations, salinity and microphytoplankton abundances comparing all surveys (Figure 3).Considering the unsignificant local precipitation rates, probably high estuarine outflows were forced by the wind in association with remote precipitation, which occurred in São Paulo, in the vicinity of Congonhas airport (IPMET, 2007).During both September 2007 surveys, mixotrophic dinoflagellates, such as Dinophysis acuminata were abundant.Gymnodinium species and Chattonella cf.antiqua, also potentially harmful bloom-forming species (Hallegraeff, 2004), were also abundant during September 2007.Dinophysis acuminata showed densities of 10 4 cell.L -1 in September 2007 ST survey, above the limits (10 3 cell.L -1 ) as a risk for public health due to potential DSP toxin production (Mafra-Júnior et al., 2006).Blooms of Dinophysis acuminata and Gymnodinium species, as well as other dinoflagellates, have been observed after strong rainfall events and in periods of low turbulent mixing in estuaries and bays (Hallegraeff, 2004).In September 2007, these dinoflagellates were able to compete with Skeletonema cf.costatum.
Fluxes of organic material through rivers and estuaries affect coastal phytoplankton growth and composition (Klug and Cottingham, 2001), especially through mixotrophic nutrition where inorganic nutrients (Arenovski et al., 1995) or light (Bockstahler and Coats, 1993) are low.The magnitudes of CDOM light absorption, proxy for DOM concentration, showed highest mean value in September 2007 NT survey (Ciotti, A.M. unpublished data).Humic substances must undergo chemical transformations to become available to bacteria that in turn support mixotrophic organisms.DOM degradation, mainly photochemical, can be an additional source of ammonia and phosphate (Rochelle-Newal and Fisher, 2002).A number of anthropogenic substances (e.g., chlorine) in domestic sewage may perhaps alter DOM characteristics and its availability to microorganisms.Nonetheless, the effects of sewage inputs on chemical and optical characteristics of DOM have been little studied (Staehr et al., 2009), and its consequences to changes in phytoplankton community structures must be addressed.
Tides and precipitation rates were also related to some phytoplanktonic assemblages.Chlorophyll-a concentration and microphytoplanktonic cell densities, mainly Chaetoceros, Pseudo-nitzschia species, Thalassiosira, Guinardia striata and Leptocylindrus minimus were higher during neap than STs, and also after periods of rainfall.Only Hemiaulus hauckii and Guinardia flaccida showed larger densities during STs, (H.hauckii during both 2006 surveys and September 2007; and G. flaccida during summer surveys), probably a result of their larger cell size and hence higher sinking rates (Scharek et al., 1999).On the other hand, Skeletonema cf.costatum, Prorocentrum minimum, Prorocentrum gracile and Scrippsiella trochoidea occurred in all surveys, with enhanced densities during winds from the northern quadrant, high PM and low Secchi disk depths.One possible explanation for the high abundance of Skeletonema cf.costatum observed in August 2006 could be ressuspension of auxospores and resting stages (Harris et al., 1998) from bottom sediments.As S. cf.costatum is the most cited bloom-forming organism in the region, it is reasonable to assume the presence of resting stages in Santos bay sediments that could eventually nourish bloom formation.However, as for the majority of diatoms (Davidovich and Bates, 2002), the life cycle mechanisms of Skeletonema cf.costatum are not properly understood.

Conclusions
Changes in microphytoplankton communities observed in Santos bay coastal waters are episodic, mainly related to wind speed and direction and tides, and to a lesser extent, to precipitation rates.Chain forming diatoms with small or elongate cells ("mixers") occurred throughout the surveys, increasing during NTs.Densities of larger celled diatoms, such as Hemiaulus hauckii ("sinkers") increased during ST, after periods of high precipitation rates.
Skeletonema cf.costatum, also a "mixer", was observed in high densities during August 2006 (austral winter) linked to strong N-NE winds, contrary to what has been previously reported, and suggesting an alteration in the temporal patterns of these species blooms, or the existence of local mechanisms for resuspension of resting stages and auxospores.
Strong southern winds increases surf zone energy in Santos bay, common after atmospheric frontal systems.In these situations, Asterionellopsis glacialis blooms are expected on a scale of a few hours to days.The enhancement of atmospheric instabilities will most likely alter the frequency of these blooms in the area, as well as in the remaining dissipative coastal in southern Brazil.
Mixotrophic dinoflagellates, such as Dinophysis acuminata, appear to be able to compete with mixers (bloom-forming diatoms) after intense estuarine discharges, raising concerns about uncontrolled disposals of freshwater into eutrophised systems, especially those surrounding highly urbanised areas.
Densities of Dinophysis acuminata and species of Pseudo-nitzschia seriata complex were above published safe limits, emphasizing the importance of urgent and constant monitoring programmes in the area.

Figure 1 .
Figure 1.Santos bay location and sampling stations.
Group 1 (surfers), composed by Asterionellosis glacialis (surf zone diatom-Group 1) was abundant in 2 of our 8 observations, August 2006 NT and ST surveys.Abundances increased when southern winds prevailed (notably in August 2006-ST).Wind speeds were above 5 m.s -1 in several occasions resulting from the passage of an intense and persistent atmospheric cold front.PM also increased, probably from bottom sediments resuspension.Group 2, composed of larger diatoms, Hemiaulus hauckii and Guinardia flaccida, called here sinkers, was abundant in 4 of our 8 observations March 2006-NT, March 2006-ST, March 2007-ST and August 2006-NT.Moderate north winds and rainfall, as observed in March 2006 and March 2007, were associated with high abundances of Hemiaulus hauckii and Guinardia flaccida, notably during STs.Both chain forming diatoms with large cell-sizes (diameter > 10 µm).Guinardia flaccida highest densities were observed in August 2006, NT, associated with winds from North and ressuspension, as observed for Asterionellopsis glacialis.

Figure 3 .
Figure 3. Variation among surveys of: a) Salinity (psu); b) Secchi disk depth (m); c) chlorophyll (mg m -3 ); d) microphytoplankton density (x 10 6 cell L -1 ); e) ∆S; f) particulate matter PM (mg L -1 ); g) particulate organic matter POM (as % of total PM).Surveys codes: March 2006 neap tide (M06_N), March 2006 spring tide (M06_S), August 2006 neap tide (A06_N), August 2006 spring tide (A06_S), March 2007 neap tide (M07_N), March 2007 spring tide (M07_S), September 2007 neap tide (S07_N), and September 2007 spring tide (S07_S).box plots-top of the tail is the highest score for the variable and the bottom of the tail is the lowest score; the box signifies the upper and lower quartiles and the median is represented by a short black line within the box.

Figure 4 .
Figure 4. Factorial plane of environmental variables, surveys (a) and species (b) resulting from Canonical Correlation Analysis (CCA).

Figure 5 .
Figure 5. Abundant species grouped by Spearman correlations and Canonical Correspondence Analysis observed in 2006 and 2007 surveys.Densities in x10 5 cell L -1 .a) Group 1 (G1), composed by Asterionellopsis glaciallis; b) Group 2 (G2), composed by Hemiaulus hauckii and Guinardia flaccida; c) Group 3 (G3), composed by Chaetoceros species, Thalassiosira sp1, Pseudo-nitzschia species, Leptocylindrus minimus and Guinardia striata; d) G4-showing only Dinophysis acuminata distribution; e) Group 4 (G4), composed by Navicula, Chatonella cf.antiqua, Gymnodinium sp1 and Dinophysis acuminata;; f) Group 5 (G5) Skeletonema cf.costatum, Prorocentrum spp and Scripsiella trochoidea.box plots-top of the tail is the highest score for density of the group, and the bottom of the tail is the lowest score; the box signifies the upper and lower quartiles and the median is represented by a short black line within the box.

Table 1 . Sampling surveys, total of sampling points, sampling days for neap (NT) and spring tides (ST) surveys and variables abbreviators utilized in statistical analysis. Months/ Sampling days Surveys Total of sampling points Variables
Salinity SAL (p.s.u.); ∆S;

Table 2 .
Average and standard deviation (between brackets) of WIND (speed, m.s -1 ) and PPT (precipitation rates mm.day -1 ) for surveys held in 2006 and 2007, spring and neap tides.

Table 3 .
Spearman correlations for environmental data and species.Only significant data (p < 0.05) are presented.Symbols follow Table1.

Table 4 .
Average and standard deviation (between brackets) of abundant species (cell.L -1) for surveys held in 2006 and 2007, spring and neap tides.

Table 5 .
Abundant species observed in different wind, salinity and tidal conditions, in 2006 and 2007 surveys.