JOURNAL OF PLANKTON RESEARCH j VOLUME 31 j NUMBER 10 j PAGES 1249 – 1263 j 2009 Evaluation of estuarine mesozooplankton dynamics at a fine temporal scale: the role of seasonal, lunar and diel cycles SÓNIA COTRIM MARQUES1*, ULISSES MIRANDA AZEITEIRO1,2, FILIPE MARTINHO1, IVAN VIEGAS1 AND MIGUEL ÂNGELO PARDAL1 1 INSTITUTE OF MARINE RESEARCH (IMAR), C/O DEPARTMENT OF ZOOLOGY, UNIVERSITY OF COIMBRA, SCIENCE AND TECHNOLOGY, UNIVERSIDADE ABERTA, 3004-517 COIMBRA, PORTUGAL AND 2DEPARTMENT OF 4200-055 PORTO, PORTUGAL *CORRESPONDING AUTHOR: [email protected] Received April 21, 2009; accepted in principle May 28, 2009; accepted for publication July 11, 2009; published online 1 August, 2009 Corresponding editor: Mark J. Gibbons In order to study the influence of physical forcing at different temporal scales, zooplankton was sampled at a fixed station located at the mouth of Mondego estuary (southern Europe). Samples were collected during diel cycles, over neap and spring tides. Zooplankton abundance and diversity were estimated for each sampling period. Holoplankton were dominated by three taxa: Copepoda, with 48% of total zooplankton abundance, Cladocera (16%) and Medusae (12%). Meroplankton occurred mainly as barnacle and decapod larvae. Copepoda were the most diverse group, represented by 26 species followed by Decapoda larvae (21) and Medusae (16). In order to assess significant differences between seasons, a univariate analysis was carried out. Higher abundance and diversity were found in warm months, particularly at neap tides, when water temperature and salinity were higher. Multivariate analysis revealed significant seasonal differences in species composition. The estuarine community was strongly dependent on allochthonous events, such as tidal exchange and river inflow. The results of our study show that the period of higher river flow, coincident with winter, resulted in changes in the zooplankton community. Short-term temporal variations in the species composition and abundance were also attributed to tidal and diel cycles. Zooplankton reached significantly higher densities at night (P , 0.05), suggesting the occurrence of vertical migrations. By emphasizing the importance of different timescale changes in the zooplankton community of the Mondego estuary, this study will be useful for the design of more efficient sampling programs, aiming at documenting changes in the zooplankton at a broad but also at a fine temporal scale. I N T RO D U C T I O N One of the major goals in marine ecological research is to understand the response of planktonic communities to physical forcing. In addition, trophic relationships and environmental parameters (such as temperature, salinity, turbidity and pollution) are known to influence the structure of zooplankton communities (Roman et al., 2001; Uriarte and Villate, 2005). In estuarine and coastal areas, the environmental conditions fluctuate at different time and space scales, subjecting organisms that inhabit the pelagic realm to tidal, diurnal and seasonal environmental changes. Knowledge of the variability of zooplankton communities at different spatial and temporal scales is therefore a prerequisite to understanding their dynamics. Physical processes such as tidal currents and estuarine circulation associated with coastline configurations and bottom topography are well-known phenomena that contribute to horizontal patchiness of zooplankton (Alldredge and Hamner, 1980; Jessopp et al., 2007). The study of these rather complex hydrodynamic processes, however, is often complicated because of their multiple doi:10.1093/plankt/fbp068, available online at www.plankt.oxfordjournals.org # The Author 2009. Published by Oxford University Press. All rights reserved. For permissions, please email: [email protected] JOURNAL OF PLANKTON RESEARCH j VOLUME interactions. Moreover, biological mechanisms may also account for a significant part of the temporal variation in zooplankton community structure. Among the biological mechanisms involved in temporal changes in zooplankton distribution, active vertical migrations are the most consistent (Rawlinson et al., 2004). Diel vertical migrations (DVMs) occur in a wide range of both freshwater and marine zooplankton taxa, probably representing the largest synchronized animal migration in terms of biomass on the planet, and is of a paramount importance for ecosystem functioning and carbon cycling (Hays, 2003). The majority of DVM studies have been focused on freshwater Cladocera and Copepoda (e.g. Dauvin et al., 1998; Armengol and Miracle, 2000; Haupt et al., 2009) since these organisms are good swimmers and show vertical movements of higher amplitude. Despite the fact that in recent years there has been an apparent growing interest in this research topic in marine ecosystems (Hays, 2003), few studies have been made in southern European coastal systems (e.g. Queiroga et al., 1997; Rawlinson et al., 2004; David et al., 2005; Morgado et al., 2006). These studies focused mainly on macrozooplankton populations or on specific species. The mesozooplankton abundance and specific composition are well documented in the Mondego estuary (Azeiteiro et al., 1999; Marques et al., 2006; Primo et al., 2009), but the results obtained in previous studies were based on monthly sampling and did not address shortterm variations. However, to examine stresses on a community, the composition, the structure and dynamics should be reasonably understood. Up to this date, only one study addressed the dynamics of decapod crustacean larval under different lunar phases and during nyctemeral cycles in this estuary (Gonçalves et al., 2003). Thus, the present study focuses on short-term temporal and spatial variations in the mesozooplankton community of the Mondego estuary. This estuarine system covers a marine, brackish and freshwater system, providing useful data for comparison with other estuaries. The specific aims were (i) to describe the taxonomic composition and seasonal patterns of zooplankton at different lunar phases and (ii) to determine the short-term changes in the community in relation to the diel light/ dark cycle and the semidiurnal tidal cycles. METHOD Study site The Mondego River estuary, located on the west coast of Portugal (408 080 N, 88 500 W), has an area of 8.6 km2 31 j NUMBER 10 j PAGES 1249 – 1263 j 2009 and a volume of 0.0075 km3 (Fig. 1). The hydrological basin of the Mondego, with an area of 6 670 km2, provides an average freshwater flow rate of 8.5 109 m3 s21. This system is located in a warm temperate region with a basic Mediterranean temperate climate. At about 5.5 km from the sea, the river branches into two arms (north and south), which converge again near the mouth. Tides in this system are semi-diurnal, and at the inlet the tidal range is 0.35– 3.3 m, with an average residence time of 2 days for the north arm and 9 days for the south arm. Sampling and laboratory analysis The zooplankton community was sampled at one fixed station at the mouth of the estuary (Fig. 1), during 1 year: June, September, October, December 2005 and March, April 2006, over neap and spring tides. In this area, the influence of both the river flow and neritic waters is strong. The sampling station was characterized by depths of 6 – 13 m. Zooplankton samples were collected by horizontal tows, using a bongo net of 335 mm mesh sized equipped with a hydro-bios flowmeter to calculate the filtered volume (average 20 m3). Considering copepods, the mesh size used certainly under-estimates the early life stages and smaller organisms. Yet, only adult copepods were taken into account in this study and smaller organisms are discussed with some caution. Hourly samples were collected from subsurface and 1 m above the bottom, over diel cycles. The samples were classified as day, from sunrise to sunset, and night, periods. Net samples were immediately preserved in a 4% borax-buffered formalin seawater solution for further analysis. A total of 388 samples were examined and sub-sampled using a Folsom-splitter, and a minimum of 500 individuals (all categories) were counted. The zooplankton analysis was carried out under a stereoscopic microscope and identification was made to species level whenever possible. Abundance estimates were standardized as the number of individuals per m3 for each depth level. In parallel with the zooplankton tows, water temperature (T), salinity (S), dissolved oxygen (DO) and pH were recorded in situ with appropriate sensors (WTW) at both depths; turbidity was measured using a Secchi disk. Also, at each sampling event, a 500 – 1000 mL water sample was filtered for determination of chlorophyll a (Chl a) and suspended particulate matter (SPM). All data concerning precipitation and freshwater runoff were acquired from the Portuguese Water Institute (INAG, www.inag.pt) stations Soure 13/01G and Açude Ponte de Coimbra 12G/01A (nearby city of Coimbra, 1250 S. C. MARQUES ET AL. j ZOOPLANKTON DYNAMICS AT A FINE TEMPORAL SCALE Fig. 1. Map of the Mondego estuary and location of sampling site. since no meteorological station was present in the study area), respectively. Statistical analysis The zooplankton community structure was described for each sampling occasion, over the temporal and spatial scales, by determining their abundance (ind.m23) and two diversity indices: Shannon diversity index (H0 ) and species richness (defined as the number of species). The comparison of these indices will indicate whether diversity variations are due to a change in the number of species, or a modification in the relative contribution of taxa. The general linear model (GLM) procedure was used for two analyses of variance (ANOVAs). The first analysis considered spatial and temporal variations in zooplankton abundance and diversity at a longer timescale (season, tide and depth). The second analysis considered variation at shorter timescales (diel, tidal cycle and depth) within seasons. When the differences were significant, Bonferroni adjusted pairwise comparisons amongst groups were computed. Prior to statistical analyses, normality and homogeneity of variance were checked for all data (Kolmogorov –Smirnov test and Levene’s test, respectively). When necessary, log (x þ 1) or square-root transformation was performed. Univariate data analysis was performed in SPSS v16.0 software (SPSS). Additionally, multivariate statistical analyses were performed to examine the variation in community composition through non-metric multidimensional scaling (MDS) plots. These analyses were based on triangular matrices of the Bray –Curtis similarities, by performing a fourth-root transformation of the species abundance data in order to down-weight the influence of highly abundant species, using only those taxa that represent .0.1% of the zooplankton community. One-way analysis of similarity (ANOSIM) was subsequently used to test whether the community structure differed significantly between groups. The contribution of each species to the average Bray – Curtis similarity within each group was analyzed using the similarity percentage procedure (SIMPER), and these results were used to help interpret the faunal change responsible for the pattern observed in the MDS ordination. Finally, the relationships between the hydrological variables were analyzed separately, using the BIOENV procedure developed by Clarke and Ainsworth (Clarke 1251 JOURNAL OF PLANKTON RESEARCH j VOLUME and Ainsworth, 1993), to define the environmental variables that better explain the assemblages’ structure. The BIOENV procedure (using Spearman’s rank correlation method) was used on log-transformed environmental data. All these analyses were carried out by means of the PRIMER-6 software (Clarke and Gorley, 2006). R E S U LT S Hydrological conditions The hydrological conditions recorded during the study period are shown in Fig. 2. A clear seasonal pattern of precipitation and runoff was observed, with the highest values in the winter. However, in the hydrological year of 2005/2006, high values were also observed in March and April. As in most river-influenced estuaries, salinity was greatly driven by the seasonality of the Mondego river freshwater discharge (Fig. 2). Mean salinity values were generally higher in spring and summer months at both depths (surface: 29.3 + 2.65, bottom: 32.37 + 4.09) (Fig. 3). At the surface, a sharp decrease in salinity was observed in autumn and winter (mean 11.6 + 7.24), showing the occurrence of the intrusion of low saline water. At the bottom, this turnover was more pronounced in the winter (average of 23.3 + 8.4). The mean salinity gradient between the high water and low water ranged up to 4.0 at both depths, being more pronounced in autumn/winter. The water temperature was only slightly different between surface and bottom (1.1 + 0.98C on average), suggesting thermal homogenization of water column (Fig. 3). The seasonal cycle was 31 j NUMBER 10 j PAGES 1249 – 1263 j 2009 marked, with high water temperature values recorded in spring and summer (mean 17.7 + 0.688C) and lower in autumn and winter (mean 13.1 + 0.68C). Secchi disc visibility decreased with increasing river discharge, inversely with the load of SPM, which generally increased. Fluctuations of SPM values over the tidal cycle were more pronounced at surface depths than on the bottom. In the winter, exceptionally high values of SPM were measured for both depths (average.140 + 5.1 mg L21, Fig. 3). Chlorophyll a values generally increased in spring and autumn being higher at the bottom, independent of the tidal cycle. Zooplankton abundance and diversity: temporal and spatial distribution Pooling all samples, a total of 112 species of zooplankton belonging to 18 taxonomic groups were identified in this survey (Table I). Generally, holoplanktonic taxa (80 + 20%) predominated over the meroplanktonic organisms (20 + 3%). Concerning holoplanktonic abundance, the highest values were attained by three taxa: Copepoda, with 50 + 12% of total abundance, Cladocera (13 + 12%) and Medusae (10 + 11%). Meroplankton occurred mainly as barnacle and decapod larvae (4 + 10 and 5 + 4%, respectively). Copepoda were the most diverse group, represented by 26 species followed by Decapoda larvae (21) and Medusae (16). Cumacea and Cladocera were composed of nine and seven species, respectively. At longer timescales, total zooplankton abundance showed highly variability within the study period (Fig. 4), ranging between 3 and 2548 ind. m23, with an average of 272 + 406 ind. m23. Zooplankton abundance Fig. 2. Monthly variation of precipitation (mm) and river runoff (dam3) from June 2005 to April 2006. 1252 S. C. MARQUES ET AL. j ZOOPLANKTON DYNAMICS AT A FINE TEMPORAL SCALE Fig. 3. Environmental parameters recorded during each sampling period (mean + SD). T, temperature (8C); Chl a; chlorophyll a (mg m23); SPM, suspended particulate matter (mg L21) (open circle, surface samples, closed triangle, bottom samples). exhibited significant seasonal differences for all biological descriptors (Table II), being generally higher in the spring and summer. Only Mysidacea and Chaetognatha exhibited higher abundance in colder seasons. Regarding the Shannon diversity index and species richness, higher values were also significantly recorded in warm months (Table I). Regarding the lunar phase, total zooplankton abundance was generally higher at neap tide (310 + 411 ind. m23; spring tide 251 + 405 ind. m23), with the exception of the winter period, where abundance was significantly higher at spring tide (Table II). In addition, a significant interaction between 1253 JOURNAL OF PLANKTON RESEARCH j VOLUME 31 j NUMBER 10 j PAGES 1249 – 1263 j 2009 Table I: Mean abundances (ind. m23) for the principal zooplankton taxonomical groups (% contribution to the total zooplankton density, average+SD) and main taxa identified in Mondego estuary for each season and lunar phase (neap and spring) and sampled depth [surface (S) and bottom (B)] Shaded boxes indicate significant differences in the three-way interaction between season, tide and depth (ANOVA, Bonferroni, P , 0.05). Bold numbers indicate the significant higher depth. Values in boxes represent higher abundance tide. Others include Amphipoda, Polychaeta larvae, Echinodermata larvae, Bivalvia and Gastropoda post-larvae, Cumacea and Isopoda groups. –, not found. lunar phase and depth (Bonferroni, Table II) indicates high abundance near bottom. Although species richness values were not significantly different between lunar phases, higher values were recorded at neap tide in warm months and at spring tide in colder months. Shannon diversity index was significantly lower at neap tide, driven by the dominance of a few species. In the spring, the marine copepod Acartia clausi and barnacle larvae constituted the bulk of zooplankton abundance at neap tide, contributing 34 and 29%, respectively (Table I). At spring tide, the estuarine copepod Acartia tonsa showed the highest values (20%), followed by the appendicularian Oikopleura dioica (16%) and the hydrozoan Muggiaea atlantica (14%), being more evenly distributed. In summer, A. clausi composed more than 20% of total abundance. While at neap tide Penillia avirostris was also an important contributor (20%), at spring tide Temora longicornis (17%) and O. dioica (18%) were the most abundant. In autumn, unlike former seasons, the proportion of the copepods increased and species composition was slightly altered. At both tides, the freshwater copepod Copidodiaptomus numidicus exhibited the highest abundance (56 and 31% of total abundance, respectively). Up to this point, a decrease in abundance was observed throughout the seasons. However, the exception noted for winter, at spring tide, was due to A. clausi (38%) and the marine cladoceran Evadne nordmanni (30%). At neap tide, the main contributors were C. numidicus (18%), A. clausi (15%) and Carcinus maenas larvae (14%). Another strong response was observed for depth. With few exceptions, higher abundance was found near the bottom, with an average of 302 + 413 ind. m23, when compared with the surface (258 + 403 ind. m23) (Tables I and II). Considering a short timescale, zooplankton community was strongly affected by the phases of the tidal and 1254 S. C. MARQUES ET AL. j ZOOPLANKTON DYNAMICS AT A FINE TEMPORAL SCALE Fig. 4. Variation in total zooplankton abundance (ind. m23) at each season and lunar phase according to tidal cycle (LT, low tide; HT, high tide) and diel phase at each sampled depth. Arrows delimit each lunar phase. diel cycles. A detailed analysis of the zooplankton abundance and diversity highlighted the differences between the daylight and dark periods (ANOVA, abundance: F ¼ 45.872, P , 0.001; species richness: F ¼ 51.168, P , 0.001; Shannon diversity index: F ¼ 11.880, P ¼ 0.001) with significantly higher values at night, suggesting the occurrence of vertical migration (Figs 5 and 6). Concerning the tidal cycle, an inconsistent pattern was observed for zooplankton abundance and diversity (Figs 4 and 5). Community structure: seasonal patterns The results of multivariate analyses demonstrated that there were significant differences in the zooplankton community structure, defining four major groups (Fig. 6). A clear cycle of species abundance and composition is revealed by the MDS plot, in which the samples are grouped by season. Group I represents the spring samples, which occurred in warmer water. Group II, which clustered closely to Group I, occurred concomitantly in warm water and was associated with higher salinities, representing summer samples. Group III included only samples from autumn. Group IV was associated with higher river flows and represented the winter communities. Despite moderate convergence of shared species between the four groups obtained by the MDS (stress ¼ 0.18, Table III), divergence was statistically significant at all groups (ANOSIM global R ¼ 0.651, P , 0.001 for 1255 JOURNAL OF PLANKTON RESEARCH j VOLUME 31 j NUMBER 10 j PAGES 1249 – 1263 j 2009 Table II: Results of factorial three-way ANOVA to evaluate significant differences in biologic descriptors between season (S), tide (T) and depth (D) and their interactions Season Tide Depth Variables F- value F- value F- value Zooplankton Species number Shannon diversity Holopankton Meroplankton Copepoda Cladocera Medusae Cirripedia Appendicularia Decapoda Mysidacea Chaetognatha 14.65** 45.130** 37.818** 9.682** 138.139** 8.032** 54.045** 129.539** 230.091** 129.756** 28.991** 11.075** 80.939** S SP W A SP S A W S SP A W S SP W A SP S W A S W A SP S SP W A SP S W A SP S W A SP S W A SP W S A A W SP S A S SP W 0.187 0.208 20.438** 0.473 5.314* 4.478* 26.749** 5.838* 43.776** 121.204** 9.936* 15.225** 0.014 NS NS SN SN NS NS SN SN NS SN NS NS NS 9.892* 15.054** 7.849* 10.336* 12.838** 9.373* 1.585 6.806 0.536 38.202** 17.228** 9.095* 9.546* B B B B B B B B B B B B B S S S S S S S S S S S S S S3T F- value S3D F- value T3D F- value S3T3D F- value 21.704** 17.345* 10.395** 24.212** 19.659** 21.329** 76.946** 42.841** 58.416** 5.099* 6.718** 3.776* 1.127 0.992 0.590 1.849 0.730 3.382* 0.864 0.427 1.991 0.087 0.961 1.007 1.423 5.647* 0.211 0.799 2.066 0.288 0.418 0.962 0.633 1.923 0.072 2.451 0.004 0.382 1.120 6.730 1.512 0.085 0.674 1.088 0.710 2.009 2.354* 0.487 6.949** 0.135 0.281 1.025 Pairwise comparisons are listed from highest mean (left) to lowest mean (right) (Bonferroni, P , 0.05). Interactions are highlighted in Table I. Season: SP, spring; S, summer; A, autumn; W, winter. Tide: N, neap; S, spring. Depth: S, Surface; B, bottom. P 0.05, *P 0.05, **P 0.001. all paired comparisons tests). Since most of these species occurred in all seasons, differences between groups appeared to be mostly the result of the variations in its relative abundance. Yet, markedly differences were noted between spring/summer and autumn/winter samples driven by a distinct zooplankton composition. The BIOENV procedure showed that among all possible combinations of the eight environmental variables, water temperature and salinity (that fluctuated with river discharge) were the major variables influencing the faunal pattern showed in the MDS ordination. Table IV shows that SPM, in combination with the above parameters, was also closely correlated with the zooplankton community structure (third best result; rho ¼ 0.501). Community structure: tidal and diel fluctuations As described earlier, zooplankton community was affected by the phases of the tidal and diel cycles. During the study period, the most abundant taxa showed significant diel variations in surface waters. Other taxa were not present in large enough numbers over the study period for the analysis. Concerning tidal cycle, differences in zooplankton community structure were noted between different tidal cycle phases (SIMPER). Such differences were greatest between high water (+1 h) and low water (+1 h) than for other pairwise comparisons. Table V shows the species contributing most to the dissimilarity observed between high and low tide. A general pattern regarding the community structure was observed throughout the sampling period. Although differing somewhat in abundance, the zooplankton composition in spring and particularly in summer were quite similar. Marine species which made up the bulk of the zooplankton community were responsible for the differences observed between the two tidal phases, with higher abundance at high water. Particularly in spring, the estuarine species A. tonsa was observed with high abundance at low water, when the extrusion of estuarine waters is at its maximum. Conversely, the zooplankton community differed considerably in species composition in the autumn and winter, showing a more pronounced effect of the tidal cycle during periods of high river flow. While there was a trend towards a freshwater community during low tide (e.g. copepods Copidodiaptomus numidicus, A. robustus, Daphnia longispina and Ceriodaphnia sp.), a large number of species with marine affinities were found at high tide. Considerably high numbers of the decapod Carcinus maenas larvae (Zoea 1), during ebb and low tide at neap tide were also found during the winter period. DISCUSSION Analysis of the zooplankton community indicated a seasonal pattern of high abundance and species richness in warm months, common to most European estuaries (Dauvin et al., 1998; Rawlinson et al., 2005). A seasonal succession of numerically dominant taxa was evident 1256 S. C. MARQUES ET AL. j ZOOPLANKTON DYNAMICS AT A FINE TEMPORAL SCALE Fig. 5. Species number and Shannon diversity index at each season and lunar phase according to tidal cycle (LT, low tide; F, flood; HT, high tide; E, ebb) and diel phase (D, day; N, night) at each sampled depth. Shaded box represents night period. over the study period, due to pulses of abundance variation in several holoplanktonic and meroplanktonic taxa. As expected, the most abundant and diverse holoplanktonic taxa were the calanoid copepods. This taxonomic group comprised 34 to 72% of total zooplankton, a usual contribution for marine coastal areas all over the world, especially for marine and brackish part of estuaries (Humes, 1994; Kibirige and Perissinotto, 2003; David et al., 2005; Leandro et al., 2007), as already stated for this estuarine system (Azeiteiro et al., 1999; Marques et al., 2006, Primo et al., 2009). Considering meroplankton, the decapod larvae were the most important contributors. This is not surprising, since the local benthic community is composed of important decapod populations (Baeta et al., 2005; Viegas et al., 2007), which in turn have life histories that include a planktonic larval phase. The Mondego estuary food chain supports an important fish community (Martinho et al., 2007). Studies performed by Martinho et al. (Martinho et al., 2008) and Dolbeth et al. (Dolbeth et al., 2008), who analyzed the feeding ecology of the main fish species of the Mondego estuary, concluded that copepods were an important component of fish diet, mainly for the juveniles of the European sea bass Dicentrarchus labrax. Consequently, high zooplankton 1257 JOURNAL OF PLANKTON RESEARCH j VOLUME 31 j NUMBER 10 j PAGES 1249 – 1263 j 2009 Fig. 6. Multidimensional scaling (MDS) plot for zooplankton assemblages sampled in the Mondego estuary. Dashed line shows the division between the four main groups. circle: ebb tides, triangle: neap tides. Table III: SIMPER analysis (cut-off 50%) for zooplankton densities per group/season (average group similarity) determined with the MDS Group and characteristic species Group I. Spring (53.59) Acartia clausi Oikopleura doica Muggiaea atlantica Diphyes sp. Nauplii cirripedia Evadne nordmanni Lizzia blondina Group II. Summer (60.47) Acartia clausi Temora longicornis Oikopleura doica Penilia avirostris Muggiaea atlantica Sagitta friderici Podon polyphemoides Group III. Autumn (52.45) Copidodiaptomus numidicus Temora longicornis Sagitta friderici Acanthocyclops robustus Acartia clausi Group VI. Winter (43.87) Acartia clausi Carcinus maenas Zoea 1 Copidodiaptomus numidicus Temora longicornis Clausocalanus arcuicornis %, percentage explained. % Faunal group Table IV: Results of the BIOENV analysis of the influence of the eight environmental parameters on zooplankton community in Mondego estuary Best singles 9.6 9.4 8.1 7.9 7.4 7.0 6.8 Copepod Appendicularian Siphonophore Siphonophore Cirriped Cladoceran Hydromedusa 10.7 10.3 8.4 7.8 6.0 5.7 5.5 Copepod Copepod Appendicularian Cladoceran Siphonophore Chaetognatha Cladoceran 15.1 10.9 9.4 8.5 6.7 Copepod Copepod Chaetognatha Copepod Copepod 18.6 12.7 10.4 7.4 5.5 Copepod Decapod Copepod Copepod Copepod Best set Variable name Rho Variable name Rho Temp Sal SPM pH Chl a 0.513 0.304 0.265 0.216 0.079 Sal, temp SPM, temp SPM, temp, sal DO, sal SPM, Temp 0.513 0.513 0.501 0.497 0.491 The six single variables most highly correlated with the biota (best single) and the six sets of variables which give the highest correlation (best set) and the Spearman’s correlation coefficient (rho) are reported. densities may lead to important trophodynamic effect on the fish community. Furthermore, pelagic filter feeders, such as gelatinous zooplankton and larger zooplankters, play an important role in this estuarine ecosystem. The gelatinous planktonic species also showed abundance peaks during summer months, which have been reported to have a significant impact on the zooplankton community (Azeiteiro et al., 1999). When jellyfish occur in high numbers, their collective prey-consumption rate can be so high that this predation directly or indirectly controls the population size of other zooplankton organisms, such as fish larvae (Hansson et al., 2005). This may lead to a 1258 S. C. MARQUES ET AL. j ZOOPLANKTON DYNAMICS AT A FINE TEMPORAL SCALE Table V: SIMPER analyses for zooplankton assemblages considering tidal cycle: average Bray– Curtis percentage dissimilarities and ANOSIM results among high tide (HT) and low tide (LT) for each sampling occasion Mean abundance LT HT Spring neap 49.87 (ANOSIM, R ¼ 0.396, P ¼ 0.05) Nauplii cirripedia 425 34 Acartia tonsa 45 0 Acartia clausi 44 439 Oikopleura doica 31 28 Pachygrapsus marmoratus 3 63 Summer neap 43.44 (ANOSIM, R ¼ 0.427, P ¼ 0.03) Diphyes sp. 2 51 Muggiaea atlantica 3 56 Penilia avirostris 25 182 Oithona plumifera 8 97 Lizzia blondina 1 5 Autumn neap 48.89 (ANOSIM, R ¼ 0.344, P ¼ 0.05) Copidodiaptomus numidicus 119 29 Ceriodaphnia sp. 5 0 Temora longicornis 2 15 Schistomysis spiritus 0 6 0 1 Evadne nordmanni Winter neap 57.59 (ANOSIM, R ¼ 0.563, P ¼ 0.03) Copidodiaptomus numidicus 52 0 Acartia tonsa 3 0 Carcinus maenas 26 305 Calanus helgolandicus 3 8 Daphnia longispina 1 0 Mean abundance % LT HT Spring spring 48.46 (ANOSIM, R ¼ 0.875, P ¼ 0.03) Diphyes sp. 3 69.0 Acartia clausi 2 78 Muggiaea atlantica 5 89 Acartia tonsa 16 0 Liriope tetraphylla 1 27 Summer spring 44.59 (ANOSIM, R ¼ 0.552, P ¼ 0.05) Penilia avirostris 2 31 Evadne nordmanni 1 39 Podon leuckarti 0 5 Podon polyphemoides 2 31 Acartia clausi 14 60 Autumn spring 56.81 (ANOSIM, R ¼ 0.979, P ¼ 0.03) Copidodiaptomus numidicus 46 1 Daphnia longispina 5 0 Acanthocyclops robustus 7 0 Penilia avirostris 0 3 Muggiaea atlantica 0 2 Winter spring 58.64 (ANOSIM, R ¼ 0.688, P ¼ 0.03) Muggiaea atlantica 0 18 Copidodiaptomus numidicus 65 0 Diphyes sp. 0 13 Acartia clausi 72 149 Acanthocyclops robustus 5 0 8.2 7.7 5.3 3.1 3.0 6.2 5.9 5.6 5.5 4.3 7.2 6.4 5.0 4.3 4.2 12.2 6.7 5.3 5.2 5.0 % 6.4 5.9 5.7 5.4 5.1 4.3 4.3 3.9 3.4 3.3 9.9 6.1 5.8 4.6 4.2 6.6 6.6 6.3 5.6 5.4 Mean abundance and percentage explained by the five taxa contributing most to dissimilarity. shift in the trophic structure of the pelagic community as a result of a trophic cascade effect. The temporal turnover in species composition, measured by the average similarities between seasons, had a consistent relationship with environmental predictors. Temperature was found to be the most important environmental factor, enough to determine species’ seasonal distribution. This fact agrees with other studies, which demonstrated that metabolic processes in zooplankton are related to temperature (Hirst and Kiørboe, 2002; Leandro et al., 2006). Nevertheless, it is well known that the factors that determine species distribution interact in complex ways. Temperature itself is not necessarily the proximal factor of zooplankton distribution patterns. In fact, a number of studies carried out in the Mondego estuary have emphasized the significant influence that salinity has on zooplankton abundance and composition (Azeiteiro et al., 1999; Marques et al., 2006, 2007; Primo et al., 2009) this environmental parameter being an indicator of different water masses. Hence, the speciesspecific relationship with salinity can indicate that the zooplankton distribution is related to different water masses. Previous studies have shown that winter climate, which favors high river discharge, can affect the mesozooplankton community, especially A. tonsa abundance, an important resident estuarine species in the Mondego estuary (Marques et al., 2007; Primo et al., 2009). It is well known that the diffusive and advective properties of freshwater discharge play a critical role in population distribution patterns and richness, as well as in their temporal variability (Kimmel et al., 2006). While the upper estuary was dominated by a freshwater community at periods of highest flow, a different scenario was observed for dry years, the zooplankton community being dominated by A. tonsa (Marques et al., 2007). It was suggested that the brackish population of this species may have been drifted away during higher freshwater flow. Moreover, during dry years, the decreasing trend of river flow was associated with variations of hydrological parameters such as an increase in salinity which contributed to high abundance of A. tonsa, namely in the upstream areas. While this work has emphasized the effects of physical forcing at short-term and vertical scales, as well as at horizontal spatial scales resulting from the movement of water passing through the fixed station, the brackish waters and associated species were not represented. To take into account spatial variability in population dynamics, at least a second sampling location would be useful. Indeed, this study showed that A. tonsa was not well represented (numerically) in the zooplankton community. Besides, 1259 JOURNAL OF PLANKTON RESEARCH j VOLUME high densities were found mostly in bottom waters. Salinity can be an important factor, limiting the occurrence of A. tonsa population (Marques et al., 2008) and may lead to high mortality outside the optimal salinity range. Our findings suggested that A. tonsa might have a retention mechanism to maintain its population within a range of brackish salinities, acting through the combined effect of the swimming behavior that allows regulation of vertical distribution (Kouassi et al., 2001) probably synchronized with tidal cycle, the production of diapause eggs (Castro-Longoria, 2001) and life history traits such as short generation times and high reproductive potential. Further sampling and experiments are required to fully understand these mechanisms. In dry months, low river flow makes tidal exchange the main process of water movement, which explains the higher importance of marine neritic species in the estuary. Such a composition is similar to what has been reported in other temperate estuaries (e.g. Dauvin et al., 1998). The present results showed that the period of higher river flow, coincidently with autumn and winter, resulted in changes in the zooplankton community. We suggest that this water mass is characterized by its own zooplankton community that is drifted away through upstream estuarine areas. In addition, a replacement in the community composition and abundance was observed, and the taxa responsible for the differences in assemblage structure were identified. Moreover, at times of low freshwater discharge, high species richness was observed. Exchange of water with the adjacent coastal area broadened the salinity gradient, providing a wider variety of marine and brackish habitats. Increasing the likely range of habitats increases the diversity of zooplankton species that will be accommodated, which will also result in a more even species distribution. Conversely, previous studies, conducted in this estuarine system by Marques et al. (Marques et al., 2006), found a lack of seasonality in zooplankton abundance. Besides the fact that fluctuations in abundance may be the result of specific and population patterns, the sampling design conducted in this study did not take into account the lunar cycle (spring and neap tides), which was probably a source of bias in the zooplankton abundance comparisons. As reported by Kennish (Kennish, 1990), the structure of the zooplankton community along an estuary depends not only on the ebbflood cycle, but also on the neap-spring tidal cycle, as has been concluded in the presently study. Indeed, regarding lunar phase, neap tide samples had generally greater total zooplankton abundance, with the exception of winter, where abundance was significantly higher at spring tide. At Lough Hyne (Ireland), a coastal marine system connected to the Atlantic Ocean, 31 j NUMBER 10 j PAGES 1249 – 1263 j 2009 a similar pattern in mesozooplankton abundance was observed by Rawlinson et al. (Rawlinson et al., 2005), but much debate has focused on the underlying reason for such observed pattern. On the Portuguese coast, neap high tides occur early in the evening, and many of the benthic species release their larvae during crepuscular neap high tides (Paula, 1989). These findings suggested that rhythmic cycles of larval export are best explained by the timing of the change in light intensity rather than tidal amplitude and is the main factor promoting synchrony with lunar cycle. Yet, studies concerning the effect of the lunar cycle on the zooplankton community led to the conclusion that minimum amplitude tides may be easier to swim against, and also that nights become progressively darker between the thirdquarter and new moon, potentially reducing losses to predation (Sponaugle and Cowen, 1997). In the present work, spring tide samples were carried out at full moon, when nights were more illuminated, possibly reflecting the mechanisms described above. According to McLusky and Elliot (McLusky and Elliot, 2004), two features potentially limit zooplankton abundance in estuaries: turbidity and, often more important, tidal currents. On one hand, turbidity can limit phytoplankton production and thus the food available for the zooplankton. In the Gironde estuary (France), the high SPM values, and consequently, turbidity were found to affect copepod selective feeding, given that it was responsible for a very reduced primary production, which indirectly controlled the temporal variability of copepods Acartia spp. and mysids (David et al., 2005). In the Mondego estuary, SPM, in combination with other parameters, was also closely correlated to the community structure of zooplankton. However, SPM concentration is very low comparatively with the Gironde estuary (1 g L21 and up to 10 g L21 in the bottom against 0.1 g L21 in Mondego estuary). In Chesapeake Bay, Roman et al. (Roman et al., 2001) found that some copepod species had the ability to ingest suspended sediments and detritus with their associated microfauna, allowing species to prosper on the high particulate concentration. In the Mondego estuary a question arises: is the correlation between SMP and zooplankton a consequence of nutritionally favorable or unfavorable environment? To answer to this question more sampling and experiments are undoubtedly needed. Additionally, currents can also transport zooplankton out of the estuary. McLusky and Elliot (McLusky and Elliot, 2004) stated that for overcoming such stressors and maintaining themselves inside the estuary, zooplankton species stay near the bottom and utilize as far as possible the inflowing marine currents. In our study, we found evidence of such mechanism, 1260 S. C. MARQUES ET AL. j ZOOPLANKTON DYNAMICS AT A FINE TEMPORAL SCALE since a clear higher abundance was observed near the bottom. Nevertheless, currents are tidally intermittent, and the best a species may achieve is to be carried to and fro about one location. Yet, in the Seine estuary, a high-frequency sampling concluded that vertical variability of zooplankton in the estuary should be considered more like a trade-off between flushing avoidance and predation (Devreker et al., 2008). Vertical migrations of estuarine zooplankton may serve as a way to control their horizontal transport when combined with tidal streams (Hill, 1991). Such behavioral mechanisms imply a process of tidally timed vertical migration synchrony, which has been suggested for copepods (Dauvin et al., 1998; Devreker et al., 2008), for meroplankton (Queiroga et al., 1997; Pereira et al., 2000) and for mysids (Hough and Naylor, 1992). In our study, we have evidence to suggest tidally oriented motion. Nonetheless, these results must be interpreted with some caution. Rapid variations in estuarine environmental and biological parameters require a sampling frequency high enough to capture consistent patterns of variation within the tidal scale (Devreker et al., 2008). However, it has been shown that decapod larval stages undergo vertical migrations. The most common situation in first stage zoea of Carcinus maenas is an upward migration at receding tide, in order to maximize seaward transportation (Queiroga et al., 1997). This pattern of migration was also observed for the Carcinus maenas zoea in this study. Macho et al. (Macho et al., 2005) noted that on the NW Iberian Peninsula, barnacle larvae were released during high tide. Here, we found high abundances of barnacle larvae around low tide. We can speculate that this drift could be related to the location of hard substrates along the shore. The distribution pattern could represent tidal flow of larvae, being displaced downstream by a combination of ebb tides and surface currents. However, except in the case of some vertical movements that are coupled with tidal flow in coastal regions, horizontal movements are probably just a consequence of DVM (Hays, 2003). Nocturnal vertical movements in the water column have been commonly reported for different taxa, including copepods (Rawlinson et al., 2005), and may have different causes. Synergistic benefits might have shaped the evolution of this complex behavior. However, predation and physical stress are believed to be the major selective pressure for the development of vertical migrations (for a detailed review, see Hays, 2003). Yet, there are associated consequences for higher trophic levels of this behavior. For example, Hays (Hays, 2003) noted that some predators at higher trophic levels modify their behavior to exploit the vertical movement of the food source, while at the same time possibly minimizing their own risk of predation. In conclusion, our data demonstrate that besides tidal changes, short-term variations in the zooplankton community may also result from diurnal cycles. Zooplankton abundance attained significantly higher densities at night than during the day, independent of tidal current, following the most common pattern of vertical migration: nocturnal migration, characterized by an evening ascent and a morning descent. The amplitude and pattern of migration may differ greatly between species and between ontogenetic stages of the same species (Rawlinson et al., 2004). The pattern recorded in this study may also reflect the use of zooplankton net with the 335 mm mesh size. This coarse mesh net may have under estimated the early developmental stages and smaller specimens, which may have masked some seasonal trends corresponding to the reproductive activities of the invertebrates within the estuary. A more suitable spatial and temporal survey of zooplankton has recently started with 64 and 200 micron mesh size net in Mondego estuary, in which it will be possible to evaluate the extent of this underestimation and its biological implications. The present study, by emphasizing the importance of different timescale changes in the zooplankton structure, will be useful for the design of more efficient sampling programs and to document changes in zooplankton abundance at a broad but also at a fine temporal scale. AC K N OW L E D G E M E N T S A special thanks to all colleagues that helped during this fieldwork. The authors wish to thank Juan Ignacio González-Gordillo for help in the identification of decapod larvae. Statistical advices of the colleague Susana Mendes were very helpful for the final version of the present manuscript. FUNDING The present work was supported by I.I.I. (Instituto de Investigação Interdisciplinar of the University of Coimbra) through a PhD grant awarded to S.C.M. (III/AMB/28/2005). REFERENCES Alldredge, A. L. and Hamner, W. M. (1980) Recurring aggregation of zooplankton by a tidal current. Estuarine Coastal Mar. Sci., 10, 31–37. Armengol, X. and Miracle, M. R. (2000) Diel vertical movements of zooplankton in lake La Cruz (Cuenca, Spain). J. Plankton Res., 22, 1683–1703. 1261 JOURNAL OF PLANKTON RESEARCH j VOLUME 31 j NUMBER 10 j PAGES 1249 – 1263 j 2009 Azeiteiro, U. M. M., Marques, J. C. and Ré, P. (1999) Zooplankton annual cycle in the Mondego river estuary (Portugal). Arq. Mus. Bocage, 3, 239–263. Kibirige, I. and Perissinotto, R. (2003) The zooplankton community of the Mpenjati Estuary, a South African temporarily open/closed system. Estuarine Coastal Shelf Sci., 58, 727– 741. Baeta, A., Cabral, H., Neto, J. N. et al. (2005) Biology, population dynamics and secondary production of the green crab, Carcinus maenas (L.) in a temperate estuary. Estuarine Coastal Shelf Sci., 65, 43–52. Kimmel, D. G., Miller, W. D. and Roman, M. R. (2006) Regional scale climate forcing of mesozooplankton dynamics in Chesapeake Bay. Estuar. Coasts, 29, 375–387. Castro-Longoria, E. (2001) Comparative observations on the external morphology of subitaneous and diapause eggs of Acartia species from Southampton water. Crustaceana, 74, 225– 236. Kouassi, E., Pagano, M., Saint-Jean, L. et al. (2001) Vertical migrations and feeding rhythms of Acartia clausi and Pseudodiaptomus hessei (Copepoda: Calanoida) in a tropical lagoon (Ebrié, Côte d’Ivoire). Estuarine Coastal Shelf Sci., 52, 715 –728. Clarke, K. R. and Ainsworth, M. (1993) A method of linking multivariate community structure to environmental variables. Mar. Ecol. Prog. Ser., 92, 205 –219. Clarke, K. R. and Gorley, R. N. (2006) PRIMER v6: User Manual/ Tutorial. PRIMER-E, Plymouth. Leandro, S. M., Queiroga, H., Rodriguez, L. et al. (2006) Temperature dependent development and somatic growth in two allopatric populations of Acartia clausi (copepoda: calanoida). Mar. Ecol. Prog. Ser., 322, 189– 197. doi:10.3354/meps322189. Dauvin, J. C., Thiébaut, E. and Wang, Z. (1998) Short-term changes in the mesozooplanktonic community in the Seine ROFI (Region of Freshwater Influence) (eastern English Channel). J. Plankton Res., 20, 1145– 1167. Leandro, S. M., Morgado, F., Pereira, F. et al. (2007) Temporal changes of abundance, biomass and production of copepod community in a shallow temperate estuary (Ria de Aveiro, Portugal). Estuarine Coastal Shelf Sci., 74, 215 –222. David, V., Sautour, B., Chardy, P. et al. (2005) Long-term changes of the zooplankton variability in a turbid environment: the Gironde estuary (France). Estuarine Coastal Shelf Sci., 64, 171–184. Macho, G., Molares, J. and Vásquez, E. (2005) Timing of larval release by three barnacles from the NW Iberian Peninsula. Mar. Ecol. Prog. Ser., 298, 251– 260. Devreker, D., Souissi, S. and Molinero, J. C. (2008) High frequency Eulerian sampling of the population of the calanoid copepod Eurytemora affinis in the Seine estuary: role of the tidal regime. J. Plankton Res., 30, 1329– 1342. Marques, S. C., Azeiteiro, U. M., Marques, J. C. et al. (2006) Zooplankton and ichthyoplankton communities in a temperate estuary: spatial and temporal patterns. J. Plankton. Res., 28, 297 –312. doi:10.1093/plankt/fbi126. Dolbeth, M., Martinho, F., Leitão, R. et al. (2008) Feeding patterns of the dominant benthic and demersal fish community in a temperate estuary. J. Fish Biol., 72, 2500–2517. Marques, S. C., Azeiteiro, U. M., Martinho, F. et al. (2007) Climate variability and planktonic communities: the effect of an extreme event (severe drought) in a southern European estuary. Estuarine Coastal Shelf Sci., 73, 725–734. doi:10.1016/j.ecss.2007.03.010. Gonçalves, F., Ribeiro, R. and Soares, A. M. V. M. (2003) Comparison between two lunar situations on emission and larval transport of decapod larvae in Mondego estuary (Portugal). Acta Oecol., 24, 183–190. Hansson, L. J., Moeslund, O., Kiørboe, T. et al. (2005) Clearance rates of jellyfish and their potential predation impact on zooplankton and fish larvae in a neritic ecosystem (Limfjorden, Denmark). Mar. Ecol. Prog. Ser., 304, 117– 131. Haupt, F., Stockenreiter, M., Baumgartner, M. et al. (2009) Daphnia diel vertical migration: implications beyond zooplankton. J. Plankton Res., 31, 515–524. Hays, G. (2003) A review of the adaptive significance and ecosystem consequences of zooplankton diel vertical migration. Hydrobiologia, 503, 163 –170. Hill, A. E. (1991) Vertical migration in tidal currents. Mar. Ecol. Prog. Ser., 75, 39– 54. Hirst, A. G. and Kiørboe, T. (2002) Mortality of marine planktonic copepods: global rates and patterns. Mar. Ecol. Prog. Ser., 230, 195–209. doi: 10.3354/meps230195. Hough, A. R. and Naylor, E. (1992) Distribution and position maintenance behaviour of the estuarine mysid Neomysis integer. J. Mar. Biol. Assoc. UK, 72, 869–876. Humes, A. G. (1994) How many copepods? Hydrobiologia, 292–293, 1– 7. Jessopp, M., Mulholland, O., McAllen, R. et al. (2007) Coastline configuration as a determinant of structure in larval assemblages. Mar. Ecol. Prog. Ser., 352, 67– 75. Kennish, M. J. (ed.) (1990) Ecology of Estuaries: Biological Aspects. Vol. II. CRC Press, Boca Raton, Florida. Marques, S. C., Azeiteiro, U. M., Leandro, S. M. et al. (2008) Predicting zooplankton response to environmental changes in a temperate estuarine system. Mar. Biol., 155, 531–541. Martinho, F., Leitão, R., Viegas, I. et al. (2007) The influence of an extreme drought event in the fish community of a southern Europe temperate estuary. Estuarine Coastal Shelf Sci., 75, 537–546. Martinho, F., Leitão, R., Neto, J. N. et al. (2008) Estuarine colonization, population structure and nursery functioning for 0-group sea bass (Dicentrarchus labrax), flounder (Platichthys flesus) and sole (Solea solea) in a mesotidal temperate estuary. J. Appl. Ichthyol., 24, 229 –237. McLusky, D. S. and Elliot, M. (2004) The Estuarine Ecosystem- Ecology, Threats and Management. Oxford University Press, 214 pp. Morgado, F. M., Pastorinho, M. R., Quintaneiro, C. et al. (2006) Vertical distribution and trophic structure of the macrozooplankton in a shallow temperate estuary (Ria de Aveiro, Portugal). Sci. Mar., 70, 177– 188. Paula, J. (1989) Rhythms of larval release of decapod crustaceans in the Mira Estuary, Portugal. Mar. Biol., 100, 309–312. Pereira, F., Pereira, R. and Queiroga, H. (2000) Flux of decapod larvae and juveniles at a station in the lower Canal de Mira (Ria de Aveiro, Portugal) during one lunar month. Invertebr. Reprod. Dev., 3, 183–206. Primo, A. L., Azeiteiro, U. M., Marques, S. C. et al. (2009) Changes in zooplankton diversity and distribution pattern under varying precipitation regimes in a southern temperate estuary. Estuarine Coastal Shelf Sci., 82, 341–347. Queiroga, H., Costlow, J. D. and Moreira, M. H. (1997) Vertical migration of the crab Carcinus maenas (L.) first zoea in an estuary: implications for tidal stream transport. Mar. Ecol. Prog. Ser., 149, 121–132. 1262 S. C. MARQUES ET AL. j ZOOPLANKTON DYNAMICS AT A FINE TEMPORAL SCALE Rawlinson, K. A., Davenport, J. and Barnes, D. K. A. (2004) Vertical migration strategies with respect to advection and stratification in a semi-enclosed lough: a comparison of mero- and holozooplankton. Mar. Biol., 144, 935 –946. Rawlinson, K. A., Davenport, J. and Barnes, D. K. A. (2005) Tidal exchange of zooplankton between Lough Hyne and the adjacent coast. Estuarine Coastal Shelf Sci., 62, 205– 215. Roman, M. R., Holliday, D. V. and Sanford, L. P. (2001) Temporal and spatial patterns of zooplankton in the Chesapeake Bay turbidity maximum. Mar. Ecol. Prog. Ser., 213, 215– 227. Sponaugle, S. and Cowen, R. K. (1997) Early life history traits and recruitment patterns of Caribbean wrasses (Labridae). Ecol. Monogr., 67, 177– 202. Uriarte, I. and Villate, F. (2005) Differences in the abundance and distribution of copepods in two estuaries of the Basque coast (Bay of Biscay) in relation to pollution. J. Plankton Res., 27, 863 –874. Viegas, I., Martinho, F., Neto, J. et al. (2007) Population dynamics, distribution and secondary production of the brown shrimp Crangon crangon (L.) in a southern European estuary. Latitudinal variations. Sci. Mar., 71, 451–460. 1263
© Copyright 2026 Paperzz