Dynamics of Mesozooplankton Assemblage in Relation to Environmental Factors in the Maryland Coastal Bays

: The mesozooplankton composition and dynamics in coastal lagoons of Maryland, mid-Atlantic region, USA have received little scientiﬁc attention despite the fact that the lagoons have undergone changes in water quality in the past two decades. We compared mesozooplankton abundance and community structure among sites and seasons, and between 2012, a year of higher than average salinity (33.4), and 2013 with lower than average salinity (26.6). It was observed that the composition, diversity, and abundance of mesozooplankton in 2012 di ﬀ ered from those of 2013. Barnacle nauplii were abundant in 2012 contributing 31% of the non-copepod mesozooplankton abundance, whereas hydromedusae were more dominant in 2013 and contributed up to 83% of non-copepod zooplankton abundance. Gastropod veliger larvae were more abundant in 2013 than in 2012 while larvae of bivalves, polychaetes, and decapods, in addition to cladocerans and ostracods had higher abundances in 2012. The abundance and diversity of mesozooplankton were explained by variations in environmental factors particularly salinity, and by the abundance of predators such as bay anchovy ( Anchoa mitchelli ). Diversity was higher in spring and summer 2012 (dry year) than in 2013 (wet year). The reduction of salinity in fall 2012, due to high freshwater discharge associated with Hurricane Sandy, was accompanied by a decrease in mesozooplankton diversity. Spatially, diversity was higher at sites with high salinity near the Ocean City Inlet than at sites near the mouth of tributaries with lower salinity, higher nutrient levels and higher phytoplankton biomass. Perhaps, the relatively low salinity and high temperature in 2013 resulted in an increase in the abundance of hydromedusae, which through predation contributed to the reduction in the abundance of bivalve larvae and other taxa. This study a joint e ﬀ Conceptualization, P.C. and Investigation,


Introduction
Planktonic organisms ranging in size 0.2-20 mm are referred to as mesozooplankton [1]; examples are copepods, decapod zoea, rotifers, gastropod veliger, barnacle nauplii, larvaceans, cladocerans and bivalve larvae. Mesozooplankton are major grazers of primary producers, and important predators of microzooplankton, including ciliates, thereby acting as a pathway for transfer of materials and energy from microbial food web to higher trophic levels [2,3]. Copepods are the most abundant of, and contribute the most biomass to, the mesozooplankton group in marine and estuarine zooplankton [4].
Although non-copepods do not constitute the bulk of mesozooplankton in United States mid-Atlantic estuaries, they occasionally become abundant [5][6][7] and play an important role in the diet of zooplanktivores [8,9]. Mesozooplankton composition, abundance and diversity in estuaries and coastal lagoons vary seasonally, spatially and between years in response to abiotic and biotic factors [10][11][12][13][14][15][16], and have been well documented in large river-dominated estuaries of the mid-Atlantic such as the Chesapeake Bay [12,17] and Delaware Bay [18,19]. Nevertheless, we know little about mesozooplankton composition and dynamics in Maryland Coastal Bays (MCBs), which are poorly flushed, eutrophic, largely polyhaline lagoonal systems that receive relatively small amount of freshwater directly from tributary rivers and creeks. Because of differences in the hydrography, the composition, abundance and distribution of mesozooplankton in MCBs may be different than those of river-dominated estuaries, and are expected to vary between years due to variations in climatic factors.
The coastal bays support commercially and recreationally important species of shellfish such as blue crabs (Callinectes sapidus) and hard clams (Mercenaria mercenaria) which rely on the successful recruitment of their meroplanktonic larvae to sustain the adult populations. The blue crab is the most abundant shellfish in the MCBs and requires high salinity waters for spawning and larval development [20,21]. In the Chesapeake Bay, C. sapidus spawn near the mouth of the estuary, or in the nearby coastal ocean [6,22]. The larvae (zoeae) are then transported into adjacent coastal waters on the continental shelf [23][24][25] where they develop through seven or eight zoeal stages before subsequently returning as post-larvae and juveniles to the estuary [25], while the adults remain resident in the system. Due to the high salinity (>26) levels in a large portion of the MCBs, it has been hypothesized that resident adult populations of the C. sapidus spawn in the high saline waters of the coastal bays and the larvae complete their development within the Bays [26]. In the 2014 updated Coastal Bays Blue Crab Fishery Management Plan Implementation document, the need for studies to determine the "level of localized reproduction and entrapment of larvae" of blue crabs in the MCBs was emphasized.
Furthermore, salinity in the MCBs increased after the opening of the Ocean City Inlet in 1933, which favored the development of hard clams (M. mercenaria) and, by the 1960s, the fishery superseded the commercial landings and value of oysters in the bay [27]. The population of hard clams is currently dominated by large and older adults; the recruitment is low and intermittent except in areas close to the Ocean City Inlet, and densities are lower than historical levels [27] although increasing. Efforts to restore the population including stoppage of the hydraulic clam dredging activities have not resulted in significant increases in clam density, the cause of which is unknown, although intense predation by species such as blue crab, and the frequent occurrence of brown tide, Aureococcus anophagefferens, blooms that negatively impact bivalve larvae have been speculated to be the cause [27]. Information on the distribution and abundance of bivalve larvae, and the environmental factors influencing their abundance would be useful for designing effective management strategies to enhance bivalve populations in the MCBs.
Within the embayments in the MCBs, nutrient levels, phytoplankton biomass and fish abundance vary spatially, seasonally and inter-annually [28][29][30], which in addition to variations in the abundance of gelatinous zooplankton predators can influence the distribution and abundance of mesozooplankton. For example, mortality of bivalve larvae due to predation by gelatinous zooplankton can potentially slow recovery of hard clams despite the implementation of effective management plans [31,32].
The objectives of the study were: 1.
To compare mesozooplankton, particularly non-copepod taxa, abundance and community structure among the MCBs embayments; 2.
To determine seasonal patterns in mesozooplankton abundance and community structure; 3.
To assess the relationships between mesozooplankton taxa abundance and abiotic and biotic factors; and 5.
To use information on decapod larval occurrence to determine whether blue crabs spawn within the MCBs.
The study addressed the following questions: Are mesozooplankton abundance and diversity higher at sites closer to the Ocean City Inlet with higher salinity, than at sites farther from the Inlet with lower salinity? Are mesozooplankton abundance and diversity higher in spring and summer 2012, a year with relatively higher salinity and lower temperature, than in 2013 when salinity was lower, and temperature higher? Do blue crabs spawn within the MCBs?

Study Area
The MCBs located in the Mid-Atlantic region on the east coast of the United States is separated from the Atlantic Ocean by barrier islands. This lagoonal system consists of five embayments ( Figure 1) that are well mixed mostly by wave and tidal actions that facilitate the exchange of water and materials between the MCBs and the nearby coastal ocean [29,33]. The system, however, is poorly flushed, has low freshwater input from groundwater, rivers and tributaries, and therefore high salinity that is close to that of the ocean except in upstream areas of rivers and creeks [28,29]. The MCBs are relatively shallow (average depth about 1.2 m) with a maximum depth of 9.3 m at the Ocean City Inlet [29]. Some characteristics of the northern bays watershed include low forest covers, non-point sources of pollution, dead-end canals, high stream nitrate, and nutrients from agricultural activities and development [34]. Increased nutrient input has resulted in eutrophic conditions in the northern bays and tributaries, although a similar condition exists in the Newport Bay located in the southern MCBs.

Sample Collection
Water quality parameters such as temperature, salinity, dissolved oxygen (DO) concentration, and pH were measured in situ using YSI 6600 sonde. Zooplankton samples were collected monthly from February 2012 to December 2013 from 13 sites (Figure 1). Sites 1-5 are located in Chincoteague Bay (CB), Site 6 in Newport Bay (NB), Sites 7 and 8 in Sinepuxent Bay (SB), Site 9 in Isle of Wight Bay (IWB), Site 10 at the mouth of St. Martin River (SMR) and Sites 11-13 in Assawoman Bay (AB). Zooplankton sampling was conducted using a plankton net, with 60-cm diameter mouth opening, 240-cm-long conical section, a mesh size of 200 µm, and a cod-end made of PVC with a volume of 1000 mL and diameter of 11 cm. Because of the shallow nature of the lagoons (average depth < 2 m), the net was towed horizontally in the upper meter of the water column similar to how zooplankton samples were collected in the shallow Barnegat Bay, NJ [35]. Sampling was for a duration of 2 min. The mouth opening of the net was equipped with a General Oceanic flow meter (Model 2030R) attached to the frame to measure the volume of water filtered by the net. Zooplankton samples were preserved in 5% buffered formalin solution immediately after collection. Biovolumes of gelatinous zooplankton (Mnemiopsis leidyi) were determined with a measuring cylinder.
Preserved plankton samples were processed in the lab using a dissecting microscope. A Folsom splitter was used to split samples when necessary. Samples were diluted to a known volume, and a minimum of 5 mL subsample was extracted with a Stemple pipette and 200-500 zooplankton individuals were counted. Cladocerans and crab zoeae were identified to species level when possible using a compound microscope and by reference to various publications [1,36,37]. Data on bay anchovy (Anchoa mitchelli) catch per unit effort (CPUE) were obtained from the Maryland Department of Natural Resources (MD DNR). Water samples for analyses of chlorophyll a and 19'-Butanoyloxy-Fucoxanthin (But-fuco), a diagnostic pigment for pelagophytes such as the Brown tide, Aureococcus anophagefferens, were collected from 0.5 m below the surface using a peristaltic pump, and placed in 2-L acid-rinsed polyethylene amber bottles. The samples were kept in ice until arrival at the lab. Water samples were then filtered using a Whatman GF/F 47 mm filter (GE Healthcare, Fisher Scientific, Hanover Park, IL, USA) and a vacuum pump for 1-4 min. A minimum of 1 L by volume of water was filtered except in situations when the filter was clogged. Chlorophyll a and 19'-Butanoyloxy-Fucoxanthin (But-fuco) concentrations were determined, as part of the integrated monitoring of the MCBs, using a High Performance Liquid Chromatography [38]. Freshwater discharge data for Birch branch at Showell, Maryland, USA were downloaded from the United States Geological Survey (USGS) website (http://waterdata.usgs.gov/nwis/monthly).

Data Analyses
The Mann-Whitney U test was used to compare mean abundance of mesozooplankton between years for each season. One-way ANOVA was used to compare differences in mean abundance of mesozooplankton among sites and embayments after log(x+1) transformation of data. The relationships between non-copepod taxa distribution and environmental factors were examined by canonical correspondence analysis (CCA), using CANOCO version 4.5 (Microcomputer Power, Ithaca, New York, USA.) [39]. Non-copepod mesozooplankton abundance was log(x+1) transformed; rare species were down-weighted when necessary and environmental data were standardized. Non-copepod mesozooplankton abundance matrices included only taxa that appeared in >10% of the sampling sites, representing >1% in 2012 and >0.1% in 2013 of the non-copepod community. Shannon diversity index [40] was computed using PRIMER 6 (PRIMER-e, Auckland, North Island, New Zealand. [41]. Spearman's rank correlation was applied to identify the environmental factors (water quality variables and biota) that were associated with diversity.

Environmental Factors
Mean monthly water temperature in 2012 varied from 5 • C in February to 25 • C in July and then decreased to about 5 • C in December (Figure 2a). In 2013, mean temperature was 30.7 • C in July and 7 • C in December. Water temperature was on average lower in 2012 than in 2013. Lower water temperatures were observed at sites close to the inlets, while St. Martin River had the highest mean water temperature (Table 1).  (Table 1). Salinity was also lower at sites in St. Martin River and Assawoman Bay relative to sites in Chincoteague and Isle of Wight Bays.
Monthly dissolved oxygen (DO) levels on the average were high in fall and winter, but relatively lower in spring and summer, although dissolved oxygen levels were > 4.6 mg L −1 at all sites throughout the duration of this study (Figure 2d). The lowest DO level was recorded in July (4.6 mg L −1 ), and the highest in March (12.7 mg L −1 ).
The seasonal patterns of chlorophyll a concentrations are similar in 2012 and 2013. Monthly levels of chlorophyll a averaged between 0.2 and 3.4 µg L −1 (Figure 2e). Chlorophyll a levels were relatively high in winter but decreased in spring, particularly in April. Chlorophyll a peaked in summer (July) in both years. Site 6 (Newport Bay) had the highest mean level of chlorophyll a in the MCBs during this study (Table 1). Monthly mean pH levels in the MCBs ranged from 7.8 to 8.2 and did not vary much within the system (Figure 2f).  Seasonal and spatial patterns of abundance of mesozooplankton are presented in Figure 4. For each of the seasons examined, there were no significant differences among sites in mesozooplankton abundance (ANOVA, p > 0.05) except in summer 2012 (ANOVA, p = 0.05) when the highest density was recorded ( Figure 4 and Table 3) at Site 7 (Tukey test, p < 0.05). When mesozooplankton abundance was compared among bays, higher abundance was recorded in Sinepuxent Bay than in any other bay in winter and summer 2012 (ANOVA, p < 0.05), as shown in Table 3.

Community Composition, Relative Abundance and Distribution of Non-Copepod Mesozooplankton in MCBs
Mesozooplankton community in the MCBs was represented by 14 major taxa dominated by copepods, but only groups that are relatively abundant are presented in Figure 5. Non-copepods contributed about 5% of the total mesozooplankton abundance. A detailed description of copepod species composition and abundance in the MCBs has been made [42]. The non-copepod community on the average was dominated by hydromedusae (3%), followed by gastropod veliger (0.5%), and barnacle nauplii (0.5%) ( Figure 5). Cladocerans (0.1%) and larval stages of decapods (0.2%), polychaetes (0.2%), bivalves (0.1%), and fish (0.1%) were also relatively abundant within the non-copepod community. Larvaceans, nematodes, amphipods, and isopods made up <1% of the zooplankton community.

Relationships Between Non-Copepod Mesozooplankton and Environmental Factors
In 2012 (Figure 9a), the first two CCA axes explained 47.1% of the cumulative percentage variance of taxa-environment relationship and the correlations were 0.89 and 0.89 for the first and second axes, respectively. The biplot in 2012 showed that the main variability (29.6% of variance) was due to the positive association of larval shrimp and bivalves with salinity (p = 0.046). In contrast, barnacle nauplii and ostracods were negatively related to salinity (p < 0.05). Gastropod and polychaete larvae showed a positive relationship with pH (p = 0.022), while crab zoeae, Evadne spp. and Pleopis sp. related inversely to pH. Hydromedusae neither associated closely with salinity nor pH in 2012. In 2013 (Figure 9b), the first two axes of the CCA explained 70% of the cumulative variation in taxa-environment relationship, and the correlation was 0.98 for the first axis and 0.91 for the second axis. The taxa-environment biplot (Figure 9b.) showed that the variability (63.7% of variance) was accounted for by the positive association of gastropod veliger with salinity (p = 0.026), and inverse relationship of hydromedusae to salinity (p < 0.05). Barnacle larvae were negatively aligned to Chl. a (p = 0.034); bay anchovy larvae were positively correlated to DO (p = 0.04), which was also plotted inversely with cladocerans, decapod zoae and polychaete larvae. Temperature, pH, and ctenophores correlated strongly with DO and were not selected by the model.

Mesozooplankton Diversity
The diversity of mesozooplankton varied temporally and spatially and was affected by salinity such that diversity was higher at higher salinity than at lower salinity ( Figure 10). The Shannon-Wiener diversity index was high (≥2.5) in 2012 at all sites except at Site 6 (H = 2.2) in Newport Bay with the lowest diversity ( Figure 10

Discussion
This is the first study to describe mesozooplankton assemblage dynamics in the MCBs. The mesozooplankton composition, abundance and diversity varied spatially and temporally in the MCBs, driven in part by fluctuations in environmental factors, particularly salinity. Lower salinity due to increased freshwater input coupled with higher ambient temperature in 2013 ( Figure 2) relative to that of 2012, perhaps promoted high densities of hydromedusae in winter and early spring. Gelatinous zooplankton can exert high predatory pressure on crustacean zooplankton species as well as on gastropod veligers, barnacle nauplii, polychaete larvae, and ichthyoplankton (fish larvae) [32,[43][44][45][46][47]. The feeding rate of ctenophores, in particular, increases with water temperature [48], possibly contributing to the decline in annual peak abundance observed for some non-copepod taxa (polychaete larvae, ostracods, bivalve and fish larvae). Nevertheless, the high freshwater inflow associated with Hurricane Sandy in fall 2012 accompanied by the relatively low salinity that occurred through spring 2013 might have contributed to the observed low abundances of the taxa in the MCBs.

Mesozooplankton Community Composition and Taxa-Environment Relationships
The percent contribution of non-copepods (e.g., decapod larvae, cladocerans) to mesozooplankton abundance in this study was highest in late spring and early summer when copepod density was relatively low [42]. Similar observations were made in the lower Narragansett Bay [16] and in other mid-Atlantic estuaries [6,19,35,49]. The increase in non-copepod taxa relative abundance that occurred when copepod density was low is important since they serve as alternate prey for zooplanktivorous fishes. Bay anchovy is abundant in summer in MCBs [50], and the juveniles usually make up the majority of estuarine bay anchovy population biomass in summer and fall [51]. Although copepods contribute the bulk of bay anchovy's diet, non-copepod taxa, especially larvae of barnacles, bivalves, and polychaetes, as well as amphipods, crab zoeae, mysids, ostracods, cladocerans, and ichthyoplankton seasonally dominate the diet in mid-Atlantic estuaries such as Long Island Sound, Cheseapeake Bay and MCBs [8,9,[52][53][54][55].
Temporal variations in freshwater discharge and salinity influenced the relative abundance of mesozooplankton taxa. For example, Barnacle nauplii dominated the non-copepod community, especially in Newport Bay from January to August 2012. During this period, freshwater discharge into MCBs was relatively low and salinity was high. Naupliar stages of barnacles were absent in fall of 2012, which might have been due to the dramatic reduction in salinity caused by excessive rainfall that was associated with Hurricane Sandy that increased freshwater inputs into the system possibly increasing mortality of the larvae and perhaps the adult stages. Mortality of barnacle larvae is higher at low salinity (<6) than at higher (>14) salinity [56], and high rainfall with associated high freshwater discharge has been reported to depress the abundance of some species of barnacles. Higher densities of barnacles were collected during dry season than wet season in coastal lagoons in Panama [56] and in Lagos Harbor, Nigeria [57]. The salinity observed in the MCBs in spring of 2013 (wet year) was much higher (23)(24)(25)(26) than the levels observed during the rainy season in these other lagoons, though the salinity in MCBs during the period Hurricane Sandy occurred might have been much lower than the values we recorded after the event. Spring and summer peaks of barnacle naupliar observed in the MCBs were also reported for the Delaware River estuary [19]. In 2013, a year of above average freshwater discharge into MCBs and lower salinity, hydromedusae dominated non-copepod community particularly in Assawoman Bay and at the mouth of St. Martin River. Larvae of polychaetes and bivalves, in addition to ostracods, and larvaceans, were better represented in the plankton in 2012 when salinity was higher, than in 2013 when salinity was lower. In estuaries that experience much horizontal gradient in salinity, seawater penetration can bring in zooplankton species from the coastal ocean during periods of low freshwater inputs [11].
There were seasonal changes in the relative abundance of mesozooplankton taxa in the MCBs. Hydromedusae were present from winter through summer, but scarce in fall (October-December). Hydromedusae were abundant in Delaware Bay River estuary during winter [18] and in Long Island estuary in April [49]. Polychaete and bivalve larvae were also abundant in winter to spring, with a shift in seasonal dominance from polychaetes to barnacles, gastropod larvae, crab zoea, barnacles, shrimps, and finally polychaete larvae in 2012. This seasonal pattern was different in 2013 when a shift in dominance of the major taxa occurred from hydromedusae to gastropod larvae, crab zoea, barnacles, and finally polychaetes.
The spatial distributions of mesozooplankton taxa were related to variations in salinity. Larvaceans were positively correlated to salinity, and found only in Sinepuxent Bay close to the Ocean City inlet in 2013. Bivalve larvae also showed a similar preference for high salinity bay areas in 2013, and gastropod larval densities decreased at sites near the mouths of MCBs tributaries. The relatively high abundance of barnacle nauplii in Newport Bay and St. Martin River in both years as well as the occurrence of polychaete and bivalve larvae in high abundance in St. Martin River in the dry year (2012) and within the open waters of Sinepuxent Bay in the wet year (2013), suggests that salinity is an important factor that influenced their distributions and abundance. As salinity increases towards the ocean, clams, crustaceans and polychaete worms dominate the benthos of the coastal bays [27], thus accounting for the high abundance of meroplanktonic larvae of these organisms around Sinepuxent Bay. The occurrence of gastropod and polychaete larvae in high abundance in embayments closest to the coastal ocean (downstream) is contrary to the findings in Pagan River, a sub-estuary of Chesapeake Bay [17] and in Mondego estuary, Western Portugal [58] where the larvae were abundant upstream, although this might be due to differences in the species reported in the systems. In the Delaware River estuary, polychaete larvae occurred from the entrance of the river to the middle of the bay and peaked in abundance from winter to spring [19]. The lowest salinity at which they were observed was 24.1. In this study, annual peaks in the abundance of polychaete larvae were recorded at salinity levels >30.

Occurrence of Blue Crab Larvae in the MCBs
Blue crab, C. sapidus zoeae, were common in July in this study which is similar to the month during which they were observed by other investigators in nearby estuaries [18,37]. The study confirms what was previously suspected based on the capture of gravid females by the MD DNR that blue crabs spawn in the MCBs. Blue crab zoeae were most prevalent in Sinepuxent Bay and the northern bays and rarely occurred in southern bay areas with adequate salinity for larval development, although no sampling occurred near the Chincoteague Inlet in Virginia. This supports results of previous studies that indicated that blue crabs hatch their eggs near the mouth of estuaries [6,22]. The occurrence of zoeae at sites located in Assawoman Bay close to tributary creeks, however, suggests that spawning occurs not only close to the Ocean City Inlet, but also in areas that are some distance away from the inlet. It is also possible that crab larvae that hatched near the inlets were dispersed into Assawoman Bay by waves and tidal action. The proximity of spawning locations to the inlet enhances easy transportation of zoeae from the MCBs into the nearby continental shelf waters where further development takes place [22,25]. Additional studies are needed, however, to determine if blue crab larvae complete their development from zoeae to megalopae within the MCBs.

Diversity of Mesozooplankton
Mesozooplankton assemblage was most diverse in spring and summer of 2012 when salinity was relatively high, especially in areas close to the Ocean City Inlet, connecting the bays to the Atlantic Ocean. Intrusion of seawater in estuaries can act as a major driver for zooplankton diversity [59]. In Pearl River estuary China, zooplankton diversity and species richness were enhanced in the middle and lower portions with high diversity at salinity level >25 [11]. The higher zooplankton diversity in spring and summer months in 2012 might have been contributed by the increase in the abundance of cladocerans, decapod larvae, and other meroplanktonic forms in the system during that period, some of which were likely transported into the bays from coastal ocean by tidal action.
The most diverse assemblage was observed in areas with high salinity, but low chlorophyll a. The observed negative correlation between chlorophyll a and diversity was because chlorophyll a levels were higher in areas close to MCBs tributaries (Newport Bay and St. Martin River) with lower salinity and high nutrient levels due to freshwater inflow [38,60], which also had lower mesozooplankton diversity than in areas close to the inlets. Howson et al. [35] observed higher zooplankton diversity in the southern bay of the Barnegat Bay, NJ, with more oceanic influence and less anthropogenic impact, than in the more nutrient enriched northern bay.

Implications of Temporal Variability in Environmental Factors on Mesozooplankton
Variability in freshwater discharge altered salinity and influenced the abundance and distribution of mesozooplankton in the MCBs. High freshwater discharge that began in fall 2012 and continued through spring 2013 lowered salinity that likely favored hydromedusae whose abundance was highest in March and April, but perhaps disfavored bivalve and barnacle larvae. During this period, salinity was about 23 (March) and 26 (April) compared with about 30.5 (March) and 34 (April) in 2012. Salinity has been reported to influence the abundance and distribution of hydromedusae in other estuaries [61,62]. In Mississippi Sound, the incidence of hydromedusae was highest at 25.1 to 30 below and above which the abundance decreased [61]. With the occurrence and high abundance of hydromedusae, ctenophores [42] and scyphozoan jellyfish in the MCBs, and their documented negative impacts on the standing stock of planktonic communities [31,45,47,48], the slow rate of recovery of the hard clam population abundance in the MCBs is not surprising. This is because, peak abundance of bivalve larvae (March) coincided with the increase in the abundance of hydromedusae (March and April), which was followed by ctenophores (April-June) especially in 2013 when salinity was relatively low [42]. Hydromedusae are known predators of zooplankton [44] and when numerically abundant, are able to exert their potential maximum clearance impact [63]. Ctenophores as predators of bivalve larvae [31,43,64] cause high mortality, with a potential of consuming up to 94.1% of the larvae during spawning periods [32].
Non-copepod zooplankton are important prey items for carnivorous zooplankton, larval fish, and adult planktivorous and forage fish such as bay anchovy. These organisms are relatively abundant when copepod abundance is low in the MCBs. During years with high freshwater input, the MCBs salinity levels in some areas may not be favorable for the development of some decapod crab zoeae and bivalve larvae (Mercenaria mercenaria) that prefer high salinity. In addition, the temperature range (~17-30 • C) tolerated by M. mercenaria larvae is reduced when salinity decreases from the optimum range of 26-27 [65]. Thus, increased temperature, reduced salinity as well as high densities of gelatinous zooplankton (hydromedusae and ctenophores) can affect the breeding success of finfish and shellfish species in the MCBs. Gelatinous zooplankton can change food web dynamics and decrease recruitment of fish and shellfish populations by preying on the larval forms or via competition for zooplankton prey [45,46].

Conclusions
Results of this study indicate that variations in freshwater discharge affect the abundance, distribution and diversity of mesozooplankton in MCBs. These changes in mesozooplankton structure likely influence trophic relationships and recruitment of finfish and shellfish larvae to the adult populations. The percent contribution of non-copepods to mesozooplankton abundance in MCBs was highest in late spring and early summer when copepod density was relatively low. The composition and abundance of mesozooplankton differed between 2012 and 2013 such that barnacle larvae were the most abundant of non-copepod mesozooplankton, whereas hydromedusae were the most abundant group in 2013. Mesozooplankton diversity varied temporally and spatially, and was higher at higher salinity than at lower salinity. Seasonally, diversity was higher in spring and summer 2012 than in fall. Additionally, in 2012 when salinity was comparatively high because of low precipitation, diversity was higher than in 2013 when salinity was relatively low. Blue crab zoeae were common in the MCBs in the summer, especially in Sinepuxent Bay and the northern bays confirming that the adults spawn in the MCBs, although further studies are needed to determine if they complete their life cycle in the bays. Additional studies are also needed to investigate the predation rates of ctenophores, jellyfish and zooplanktivores such as bay anchovy on various zooplankton taxa. This will be useful for understanding the role of non-copepod community as alternate prey for fish in the MCBs during spring and summer minimum of copepods, and will also provide more information on the impact of predation on the survival of bivalve larvae and their recruitment to the adults.