Mesozooplankton biomass and copepod estimated production in a temperate estuary (Mondego estuary): effects of processes operating at different timescales

A fundamental step in characterizing marine planktonic communities is the knowledge of the variability of their biomass and production, in order to assess the flux of matter through the food web and estimate the fate of this energy. The zooplankton assemblage’s dynamics of the Mondego estuary (Portugal) have been extensively characterized, but none have focused on biomass and secondary production. In the present study, zooplankton biomass and secondary production were investigated with a seasonal basis in a fixed station, located at the mouth of the Mondego estuary (Portugal), with emphasis on the copepod community. In order to study the influence of physical processes operating at different temporal scales, sub-surface and near-bottom samples were collected every 2 h in diel and tidal cycles, over neap tides. Copepods dominated the biomass during spring (0.585 ± 0.671 mg C m−3), summer (0.287 ± 0.383 mg C m−3), and winter (0.221 ± 0.128 mg C m−3). Lower observed biomass values were probably due to the higher mesh size used in comparison to previous studies in other estuarine systems. No differences were detected between seasons in total zooplankton biomass. Concerning diel cycle, significant variation occurred mainly in spring and winter, with substantially higher biomass values being registered at night for both the surface and bottom. Moreover, copepod biomass was consistently higher at the bottom during winter. In this season, bottom samples presented higher copepod biomass both at ebb and flood tides. From the GLM analyses run for each season, chlorophyll a, diel phase, and depth were significant factors. The considerable complexity in the processes structuring the dynamics and composition of biomass and production rates has been shown. Our results demonstrated that along spring and winter, secondary production varied according to diel phase and depth. In winter, night and bottom samples revealed consistently higher daily secondary production rates. From this, it is possible to conclude that local factors (river flow, predation, food limitation) regulate the seasonal variations in zooplankton biomass and production in the Mondego estuary.


Background
The dynamic nature of estuarine ecosystems is defined by varying physical and chemical characteristics, including fluctuations in salinity, temperature, turbidity, and oxygen in the water column (Elliott and McLusky 2002). Water circulation and land influence also induce high variability on the distribution and structure of planktonic populations, which tend to be spatially and temporally heterogeneous when compared to other aquatic systems (Marques et al. 2007). Although these habitats are considered stressful and harsh, few organisms are able to survive and adapt to these conditions and manage to escape from competition in the sea or surrounding freshwater, and often thrive in large numbers. In fact, all estuaries show a gradual reduction in diversity but not in abundance or productivity of species when compared to adjacent aquatic habitats (Elliott and McLusky 2002).
Due to their significant abundance and biomass in estuarine and nearshore areas, zooplankton play a key role as major grazers in aquatic food webs, being critical intermediaries in the flow of energy and matter through marine food chains, from primary producers to consumers at higher trophic levels, such as fish, marine mammals, and turtles (Leandro et al. 2014, Richardson 2008. The composition and abundance of estuarine zooplankton are distinct and display a non-homogeneous, variable spatial distribution (Ré et al. 2005). Zooplankton distribution is influenced by factors acting on many scales, with multiple interactions in an unstable aquatic environment. For instance, estuaries are strong advective environments as a direct consequence of the tidal currents (ebb and flood) and river flow, imposing an additional stress factor for zooplankton populations (McLusky and Elliott 2004). Moreover, zooplankton exhibit a variety of daily cycles including vertical and horizontal migrations, changes in feeding behavior, and alternative reproductive states (Hays 2003). This suggests that there is an important link between hydrological cycles and zooplankton community ecology in these habitats. In most estuarine ecosystems, copepods comprise the bulk of the holoplankton (Day et al. 1989), even outnumbering insects by possibly three orders of magnitude. The relevance of copepods is such that the reproductive and recruitment success of several pelagic fish and shellfish species of high economic value are extremely dependent on the dynamics of their populations (Beaugrand et al. 2003, Conover et al. 1995. The copepod community of the Mondego estuary is dominated by the Acartia genus (Azeiteiro et al. 1999, Vieira et al. 2003. These organisms reproduce throughout the year, and several annual growth generations can be found in temperate systems like this one because of favorable temperatures and high food availability (Kleppel 1992, Leandro et al. 2007).
Productivity has become a central and extensively studied aspect in marine plankton research over the last few decades (Runge and Roff 2000). It is of interest as a measure of energy flow through a population and as an indicator of its physiological or nutritional state (Kimmerer 1987). The zooplankton communities of the Mondego estuary have been extensively characterized and studied on various levels by a number of researchers (e.g., Marques et al. 2006, Primo et al. 2009, Vieira et al. 2003, but none have focused on the exchange of biomass and secondary production with the estuary and coastal waters. In order to better comprehend the role of zooplankton in the material flow of this ecosystem, the principle aims of this study were as follows: (1) to quantify zooplankton biomass categorized by taxonomic groups; (2) to estimate copepod production by using the global model of Huntley and Lopez (1992); (3) to evaluate the effect of processes operating at different timescales on biomass transport at the interface between the estuary and the adjacent coast: seasons, tidal, and diel cycles.

Study site
The Mondego estuary is located on the Atlantic west coast of Portugal (40°08′ N, 8°50′ W) and has a warm temperate climate. With an area of about 8.6 km 2 , it comprises two arms (north and south) that are separated by a central island (Fig. 1). The two arms have different hydrographic characteristics, being the north arm deeper (5 to 10 m during high tide, tidal range of 1 to 3 m) and presenting lower residence times (<1 day). The southern arm is almost silted up in upstream areas (2 to 4 m during high tide, tidal range of 1 to 3 m), causing freshwater to flow mainly through the north arm. Water circulation in this arm is mostly due to the tides and a small input of freshwater from a connecting tributary. Due to the differences in depth, tide penetration is faster in the north arm and causes greater daily changes in salinity, whereas temperature shows more daily variations in the south arm (Marques et al. 2002).

Field sampling and laboratory procedures
Zooplankton samples were collected at a fixed station located at the mouth of the Mondego estuary ( Fig. 1), from the summer of 2005 to the spring of 2006. The sampling station was characterized by depths of 6 to 13 m. In this area, there is a strong influence of both river and coastal waters. Seasonal sampling took place in June, October, and December of 2005 and in March of 2006, during neap tides. Collection was performed every 2 h over diel cycles, with one sample being gathered from each of two depth ranges: sub-surface and 1 m above the bottom. Concerning the bottom samples, the net was attached to a winch cable and a weight was added to the lower end of the rope. A bridle was adjusted so that the frame of the net (not the mesh netting) supported the weight. The zooplankton net was then lowered very slowly and carefully to the desired depth and dragged for 3 min. The net was raised out of the water in one motion until the plankton bucket was just above the surface. Samples were classified as day or night, with the day phase being taken from sunrise to sunset and the night phase from the corresponding period. Zooplankton was gathered by horizontal tows at low speed, using a bongo net of 335 μm mesh size (0.5 m diameter) equipped with a Hydro-Bios flowmeter to calculate the filtered volume (average 20 m 3 ± 10.3 m 3 ). Samples were fixed and preserved in a 5 % borax-buffered formaldehyde seawater solution. In parallel with the tows, water temperature (T) and salinity (S) were recorded at each collection point using a WTW Cond 330i meter at both depths. Water samples of 500 ml were also filtered for determination of chlorophyll a (Chl a) and suspended particulate matter (SPM). The bottom water samples were obtained with a Van Dorn sampler. For determination of Chlorophyll a, 500-1000 ml water samples were filtered onto Whatman GF/C glass-fiber filters followed by extraction following the protocol of Parsons et al. (Parsons et al. 1985). SPM was estimated by filtering 700 ml water through pre-combusted and pre-weighted Whatman GF/C filters and dried at 60°C for 72 h and combusted at 450°C for 8 h (APHA 1995).

Biomass determination
Organisms were separated by taxonomic group into small previously weighed aluminum capsules under Leica M80 stereomicroscopes. Samples with very high numbers of individuals were sub-sampled as necessary by using a Folsom plankton splitter. Subsequently, the capsules were placed in heat resistant acrylic multiwell trays and put in a drying oven set to 60°C for at least 24 h, as recommended (Lovegrove 1962(Lovegrove , 1966. This procedure assured water removal without affecting organic content (Ré et al. 2005). Zooplankton biomass was expressed as mg m −3 . Dry weight were converted according to carbon weight (μg C) assuming that it corresponded to 40 % of the dry weight, according to the methods employed by Omori and Ikeda (1984) and Båmstedt (1986). These values represent a realistic estimate of biomass and the nutritive value of zooplankton (Ré et al. 2005).

Estimation of copepod production
Daily secondary production rate was determined with the following equation: where P is the production (mg C m −3 d −1 ), B is the copepod biomass (mg C m −3 ), and g is the weight-specific growth rate (d −1 ). Estimates of copepod growth rates were obtained from the temperature-dependent predictive method of Huntley and Lopez (1992) with the expression: being T equal to temperature (°C). These authors have shown that for marine copepods, the temperature dependence of growth rates surpasses species differences. As a result, the exponential growth model presented should be applicable to production estimates of the entire communities, regardless of species composition or size frequency. In this study, growth was directly estimated from the weight of adult organisms for each taxonomic group and recorded temperatures.

Statistical analysis
Permutational univariate analysis of variance (PERMA-NOVA, Anderson et al. 2008) was used to test for significant differences in both environmental and biological variables. Tests were applied to non-transformed data, and based on Euclidian distances between samples, considering all the factors as fixed and with unrestricted permutation of raw data.
Environmental variables were tested with a two-way design (season/depth) and seasonal variations in zooplankton biomass by a one-way design (season). Copepod vertical distribution of biomass and production were tested for each season separately by two-way design (diel/ depth and tide/depth). PERMANOVA tests were applied with PERMANOVA+ for PRIMER software (PRIMER v6 and PERMANOVA +v1, PRIMER-E Ltd.). To identify which environmental variables influenced copepod production, Generalized Linear Models (R software; R Development Core Team 2008) were applied. Again, each season was analyzed separately and environmental variables included were the following: salinity, temperature, chlorophyll a, and suspended particulate matter. Tide, diel cycle, and depth were included as nominal variables, as well as interaction tide: depth and diel cycle: depth. Significant results in interaction terms indicate that the relationships between the explanatory variables are not the same for each condition. Variables were transformed whenever necessary (presence of extreme values): log transformation of Chl a and SPM in summer, winter, and spring and square root transformation of secondary production in spring. The Spearman correlation coefficient was used to test for collinearity between continuous variables, and variables were excluded whenever the values were superior to 0.7. A Gaussian error distribution was used for the response variables (copepod production). The selected variables included in the final model were obtained using an automatic stepwise selection (drop 1), and the Akaike information criterion (Akaike 1974) was used to select the best model.

Environmental characterization
The hydrological conditions recorded during the study period are shown in Fig. 2 Salinity and water temperature revealed a marked seasonal variation typical of temperate estuaries (Fig. 3). Salinity ranged between 31.4-35.2 in summer and 6.7-33.0 in autumn. Water temperature varied from 18.6-20.1°C in spring and 11.3-13.1°C in autumn. There were significant differences in the water temperature values between all seasons (pseudo-F = 293.90, p < 0.05), while salinity displayed considerably higher values in spring and summer (pseudo-F = 31.86, p < 0.05). A similar pattern of variation was found for SPM, which presented significantly higher values in warmer seasons (0.044 ± 0.025 mg m −3 spring, 0.033 ± 0.017 mg m −3 summer, pseudo-F = 25.01, p < 0.05). Chlorophyll a concentration increased noticeably in spring (4.8 ± 3.3 mg m −3 ) and autumn (4.6 ± 2.7 mg m −3 ) (pseudo-F = 9.41, p < 0.05).
A vertical variation of environmental parameters in the water column was also detected (Fig. 3). Mean salinity values were generally higher at the bottom during summer and autumn (post hoc t test, p < 0.05), showing no differences for the remaining seasons. At the surface, a sharp decrease in salinity was noted in autumn (6.7 ± 3.8), indicating the occurrence of freshwater intrusion. Concerning temperature, autumn presented significantly higher values for the bottom, while values increased at  the surface during spring (post hoc t test, p < 0.05). Chlorophyll a concentration was always higher at the surface (pseudo-F = 4.85, p < 0.05) in opposition to the SPM load, which had higher values at the bottom (pseudo-F = 32.69, p < 0.05).
The diel and tidal vertical variations of copepod biomass are shown in Figs. 5 and 6, respectively. Diel variation occurred mainly in spring (pseudo-F = 10.70, p < 0.05) and winter (pseudo-F = 8.31, p < 0.05), with substantially higher biomass values being registered at night for both the surface and bottom, respectively. Moreover, copepod biomass was consistently higher at the bottom during winter (pseudo-F = 6.27, p < 0.05) (Fig. 5).
Regarding tidal phases, no significant differences were observed between flood and ebb tides, even when considering vertical distribution, bottom, and surface values. In winter, bottom samples presented higher copepod biomass, both at ebb and flood tides (pseudo-F = 4.99, p < 0.05) (Fig. 6).
Secondary production values for copepods were much lower during autumn (0.04 ± 0.03 mg C m −3 d −1 , pseudo-F = 2.67, p < 0.05) (Fig. 7). Daily secondary production rates exhibited significant differences in spring and winter, for both diel and depth factors. Throughout winter, the highest values were found at night (0.13 ± 0.12 mg C m −3 d −1 , pseudo-F = 5.25, p < 0.05) and in the bottom samples (0.13 ± 0.13 mg C m −3 d −1 , pseudo-F = 2.96, p < 0.05). In spring, diel variation changed according to depth (pseudo-F = 4.68, p < 0.05). The surface presented considerably higher values during the night (post hoc t test, p < 0.05), while at the bottom, there were no noteworthy differences. Copepod production showed no variation considering tidal phases.
Relationship between copepod production and environmental factors GLM analyses identified, for each season, which environmental factor had higher influence on copepod secondary production variation, resulting in the best model that fits observed values (Table 1, Fig. 7). From the GLM analyses ran for each season, chlorophyll a, diel phase, and depth were significant factors for the models (Table 1, Fig. 7). Spring and winter displayed higher coefficients of  (Table 1).
During spring and winter, diel phase and depth influenced copepod secondary production, and thus, data were presented according to day/night periods and bottom/surface depth (Fig. 7). For spring, copepod secondary production displayed lower values during the day. Furthermore, night values showed higher variability. Also, vertical distribution seemed to vary according to diel phase: higher production in bottom samples during the day and in surface samples during the night (Fig. 7).
In the winter, copepod production was also higher during the night, both in surface and bottom samples. Additionally, copepod secondary production was generally higher in bottom samples, both during the day and night (Fig. 7). The highest variation in copepod production was observed at night in bottom samples.
The GLM analysis pointed out Chl a as the best predictor for secondary production during summer and autumn; however, the relationship varied according with the season. During summer, higher Chl a resulted in lower secondary production (negative relationship) while in autumn; increased secondary production was observed with higher Chl a (positive relationship) (Fig. 8).

Discussion
Zooplankton is composed of animals from several taxonomic groups, but crustaceans are generally dominant. In terms of the overall contribution to zooplankton biomass, the most important group found during this study was Copepoda. This is in accordance with previous studies of the zooplankton community in the Mondego estuary (Marques et al. 2006). In fact, total zooplankton biomass reflected the seasonal variation of the copepod population quite well. The results are also parallel to findings in other geographical areas, which revealed that copepods usually contribute with the majority of zooplankton biomass and species diversity in estuaries (Arashkevich et al. 2002, Leandro et al. 2007, Uye and Liang 1998. Other organisms encountered included gelatinous species (consisting primarily of medusae, chaetognaths, and appendicularians). Due to low dry weight, their contribution to biomass was minimal. Nevertheless, their total predatory impact on the zooplankton community is possibly important since these carnivores have high grazing rates (Marshalonis andPinckney 2008, Purcell et al. 1994). When numerous, they significantly affect plankton species occurrence and size composition.
Although estuarine ecosystems are generally characterized by high zooplanktonic biomass and secondary production values (David et al. 2006), the results obtained for the Mondego estuary are lower than those encountered in the literature for other regions, such as the Ria de Aveiro in Portugal (Leandro et al. 2007), Westerschelde in the Netherlands (Escaravage and Soetaert 1995), Kattegat in Denmark (Kiørboe and Nielsen 1994), and the Inland Sea of Japan (Uye and Liang 1998). These authors found numbers ranging from a minimum of 1.18 mg C m −3 d −1 (Kattegat) to a maximum of 6.85 mg C m −3 d −1 (Inland Sea of Japan) ( Table 2). Maximum production for the Mondego estuary was only 0.125 mg C m −3 d −1 when considering the contribution of copepods alone (main taxonomic group).
Methods for zooplankton biomass estimation are relatively standard and calculated values have been reported for different areas worldwide (Fernández de Puelles et al. 2003, Hays et al. 2001, Irigoien and Castel 1995, Kimmerer and McKinnon 1987, Melo Júnior et al. 2007, Rawlinson et al. 2005, Roman et al. 2002. However, data are generally obtained from samples collected with nets of mesh diameter ranging from 64 to 300 μm. It is important to note that net selectivity could have a significant influence on the generated data since it can differentially represent the dominant species and/or the proportion of developmental stages present in the environment. Most coastal mesozooplankton assemblages are composed of small-sized organisms, as well as early developmental stages of larger species, which generally dominate in terms of abundance and biomass (Turner 2004). Consequently, the differences observed in this study were probably due to the higher mesh size used in comparison to other studies 64 and 125 μm-Ria de Aveiro, 55 μm-Westerschelde, 50 μm-Kattegat and 62 μm-Inland Sea of Japan). Due to their small body size, microzooplankton, particularly ciliated protozoans, have higher weight-specific physiological rates such as feeding, respiration, excretion, and growth than large metazoans. Therefore, the frequent undersampling of small copepod species may lead to a limited view of the planktonic systems and must be considered when comparing results. In this respect, the importance of small copepods is well recognized as a fundamental link in marine food webs, serving as major grazers of phytoplankton and prey for ichthyoplankton and other larger pelagic carnivores (Turner 2004, Verity et al. 2002. While it is now recognized that many factors contribute to fish recruitment (Houde 2008), the amount of zooplankton prey available to larval fish remains without a doubt vital. Concerning the estimation of secondary production rate, the current lack of a completely accepted method for its determination is a consequence of the variety of factors affecting secondary production, being difficult to determine the efficacy of each one. Moreover, in pelagic research, efforts on zooplankton growth rates are mainly focused on copepods (Hirst and McKinnon 2001) and Table 2 Daily secondary production values encountered in different coastal regions for zooplankton Location Production (mg C m −3 d −1 ) Empirical model or method Reference Kattegat, Denmark 1.2 Weight-specific egg production rates Kiørboe and Nielsen 1994 Westerschelde, Netherlands 1.9 Growth rates (temperature-dependent growth model) Escaravage and Soetaert 1995 Inland Sea of Japan 6.9 Length-weight relationships, instantaneous growth rates Uye and Liang 1998 Ria de Aveiro, Portugal 3.7 Huntley and Lopez (1992) Leandro et al. 2007 2.9 Hirst and Bunker (2003) Mondego estuary, Portugal 0.13 Huntley and Lopez (1992) Present study Fig. 8 Generalized linear models fitted to secondary production at each season. Only significant variables at p < 0.05 are shown. Dashed lines are approximate 95 % confidence intervals not on communities of different species. As mentioned earlier, several methods can be used to determine secondary production for zooplankton, but in addition to biomass values, all of them require the determination of growth rates (Leandro et al. 2007), which cannot always be done directly. Due to the nature of the preserved samples in this study, available field parameters, and the lack of more detailed information, it was decided that the mathematical growth model described by Huntley and Lopez (1992) would be the most appropriate one to conduct the intended analysis. However, this model depends only on temperature to make predictions and does not assume that copepod/zooplankton growth may be limited by food at some points during the year (Burkill andKendall 1982, Peterson et al. 1991). Its general nature does not take into account particular characteristics of the ecosystem either. Hence, caution must be used in its application, as it tends to overestimate growth rates (Kleppel et al. 1996). Despite this, it is important to understand that this type of model can be a reliable tool for estimating secondary production because it considers the ecologically relevant parameters of biomass and growth rates, which have been extensively studied. During the study period, both copepod biomass and production rate displayed a weak seasonal effect. This lack of cyclic variation may be related to species replacement when they are seasonally substituted by each other yet maintaining average biomass (e.g., freshwater species replaced by marine ones). Changes in copepod production were only detected in autumn, which had the lowest values. During autumn, the increased prevalence of predators (e.g., mysids and chaetognaths) most likely influenced copepod community structure by favoring species with lower growth rates. Salinity values detected at the surface during this season indicated a strong input of freshwater that may also affect the planktonic community. According to Kimmerer (2002), variations on estuarine organisms' abundance and survival may occur through attributes of physical habitat that vary with river flow, and negative responses can be expected mainly from marine species due to the seaward displacement of their habitat, and therefore, population centers.
During autumn, the low production rates were positively related with chlorophyll a, showing no diel, tidal or vertical patterns. The same occurred during summer, although presenting an inverse (negative) relationship with primary productivity. These results further suggest a seasonal shift in the copepods' diet in the estuary. Contrasting with autumn, which exhibited high values of chlorophyll a, the amount of primary production in summer in the Mondego estuary was quite certainly well below the carbon requirements of the copepod population, not being sufficient enough to cover their nutritional needs. Phytoplankton appears to be the most important carbon source over the course of autumn, while in summer, alternative sources are utilized. Many researchers have noted the paradox between high copepod biomass in estuarine areas with a large quantity of particles and low phytoplankton production (David et al. 2006). This suggests the possible use of detritus (essentially originating from terrestrial plants) as a food source for estuarine copepods. For instance, Diodato and Hoffmeyer (2008) found that Acartia tonsa, the most abundant species in the Mondego estuary, consumed detritus in the Bahía Blanca Estuary (Argentina), and a study of the cycling of organic material in the Kariega Estuary (South Africa) revealed that mesozooplankton used both phytoplankton and detritus as food (Richoux and Froneman 2008).
In spring and winter, GLM analysis identified diel and vertical patterns as having a key role on production rate variations, as opposed to summer and autumn. In fact, during spring and winter, both biomass and production rates were significantly higher at nighttime than at daylight hours. Additionally, spring copepod production rates seem to show a vertical variation associated with depth (changes in depth according to diel phase); a strong evidence of diel vertical migrations. It is generally accepted that the daily sunlight cycle plays an important role in the diel vertical migration (DVM) behavior of zooplankton, which has been observed in coastal regions (e.g., Zhu et al. 2000, Devreker et al. 2008, Primo et al. 2012. Nocturnal vertical movements in the water column have commonly been reported for different taxa, including copepods (Rawlinson et al. 2005). Synergistic benefits might have shaped the evolution of this adaptive complex behavior in response to variations in the environment. Explanations for the vertical migration mechanism include environmental factors (e.g., light, temperature, salinity, oxygen, and diet), predator-avoidance, and energy and resource utilization (for a detailed review see Hays 2003). This behavior causes active transport of dissolved organic and inorganic carbon and nitrogen and plays a welldocumented role in biogeochemical cycling (Steinberg et al. 2012). Furthermore, there are associated consequences for higher trophic levels of this behavior. For example, Hays (2003) noted that some predators at higher trophic levels modify their activity to exploit the vertical movement of the food source. The food chain of the Mondego estuary supports an important fish community (Martinho et al. 2007). Studies performed by Martinho et al. (2008) and Dolbeth et al. (2008), who analyzed the feeding ecology of the main fish community, concluded that copepods were an important component of the fish diet. As a result, DVM behavior could lead to important trophodynamic effects in the fish community.
Vertical patterns were stronger during winter, when bottom samples showed increased biomass and production rates compared to surface samples. Despite autumn samples presenting higher salinity stratification, during the winter, salinity at the surface in the Mondego estuary is highly variable due to changes in freshwater flow ). The instability of the water column influences planktonic organisms to gather near the bottom in order to avoid being flushed out of the estuary. Contrary to diel and vertical distributions, variations in biomass and secondary production showed no relation with tides. This can indicate that the amount of organisms imported (flood) was equivalent to those exported (ebb) from the estuary. However, results from vertical salinity variation pointed out low stratification during spring and summer. At this time a persistent landward current seemed to occur, increasing the import of biomass and transportation to upstream areas. The higher freshwater flow detected in autumn and its strong associated advective effects (export) are mainly seen at the surface, but the same amount of biomass and production rate appear to be inputted near the bottom, hence, the absence of vertical differences. During the winter, copepods tend to agglomerate close to the bottom, which results in a lower transport of organisms.

Conclusions
Of the eight taxonomic groups identified, Copepoda, Mysidacea, and Cladocera were the most significant. Copepods represented a large percentage of the total biomass, but production calculated with the Huntley and Lopez (1992) model was lower than expected, possibly due to the large mesh size used (335 μm). Production rates revealed the influence of diel vertical patterns, with higher values being registered in depth at nocturnal periods. The Mondego estuary did not display the typical seasonal patterns of high zooplankton biomass in summer and low in winter, common to similar sites (Fernández de Puelles et al. 2003, Rawlinson et al. 2005. Nonetheless, in the western Mediterranean basin, the range of zooplankton biomass and occurrence is not always well defined, and comparisons between data of other systems should be done carefully due to the great variety of sampling methods. Tidal exchange also revealed a non-significant variation, with little difference between imported and exported biomass and production. This is not uncommon, and studies in other locations have yielded equivocal results with different explanations, as discussed by Melo Júnior et al. (2007). In a shallow coastal ecosystem like the Mondego estuary, the influence of river discharge and varying salinity levels are significant, and it is possible to conclude that local factors such as freshwater flow, predation, and food limitation may regulate the seasonal variations in biomass and production. This study represents the first attempt at estimating biomass and secondary production of the zooplankton community of the Mondego estuary. It would be important to conduct future studies with smaller mesh size nets in order to estimate the contribution of smaller sized zooplankton, as well as to define a copepod growth model specific to this ecosystem for achieving more accurate production estimates.