Printer Friendly

Variabilidad estacional de la estructura comunitaria y abundancia de copepodos, modificada por la transicion El Nino-La Nina (2010) en el Pacifico, Mexico.

Seasonal variability of copepod community structure and abundance modified by the El Nino-La Nina transition (2010), Pacific, Mexico

Copepods are an important planktonic group, and account for most of the total biomass and species diversity in pelagic marine ecosystems. In coastal and oceanic regions of the Eastern Tropical Pacific off Mexico (ETPM), they have been reported to comprise around 66-76% of total zooplankton abundance (Franco-Gordo et al., 2015; Pelayo-Martinez et al., 2017). Subsequently, they greatly contribute to zooplankton production and energy transfer to higher trophic levels. In the ETPM, copepods are exposed to a complex mix of environmental conditions, such as seasonal changes in the water column, principally regulated by the summer/fall poleward Mexican Coastal Current (MCC) (Lavin, Beier, Gomez-Valdes, Godinez, & Garcia, 2006; Salas, Gomis, Olivos-Ortiz, & Garcia-Uribe, 2006; Gomez-Valdivia, Pares-Sierra, & Flores-Morales, 2015), and a winter/spring tropical extension of the California Current (Salas et al., 2006; Godinez, Beier, Lavin, & Kurczyn, 2010). Other important hydrological characteristics include mesoscale eddies, which are usually present year-round with seasonal variations (Kurczyn, Beier, Lavin, & Chaigneau, 2012), winter/spring coastal upwelling (Torres-Orozco, Trasvina, Muhlia-Melo, & Ortega-Garcia, 2005; Lopez-Sandoval, Lara-Lara, Lavin, Alvarez-Borrego, & Gaxiola-Castro, 2009) and an extremely shallow upper layer of the oxygen minimum zone (OMZ, 9.0 [micro]mol L-1) which reaches depths of ~50 m along the coast (Zamudio, Hurlburt, Metzger, & Tilburg, 2007; Cepeda-Morales, Beier, Gaxiola-Castro, Lavin, & Godinez, 2009). These environmental conditions make this area a dynamic transitional environment that influences patterns of temperature, dissolved oxygen, and primary and secondary productivity (Ambriz-Arreola, Gomez-Gutierrez, Franco-Gordo, & Kozak, 2015; Kozak, Franco-Gordo, Palomares-Garcia, Gomez-Gutierrez, & Suarez-Morales, 2017; Pelayo-Martinez et al., 2017).

Distribution, feeding and reproduction of zooplankton are influenced by abiotic factors like bathymetry, thermocline depth, temperature, primary productivity and food availability. These factors are in turn modulated by dynamic structures such as eddies (cyclonic and anticyclonic), upwelling events and filaments, which due to their 1-100 km size are referred to as dynamic mesoscale structures (DMS). These DMS are usually caused by gravity, density gradients and wind forcing, and provoke a resuspension or sinking of particles and organisms, which generates patches of organisms according to hydrological conditions, reflected in the abundance of distinct zooplankton groups (Pelayo-Martinez et al., 2017). Within the zooplankton, copepods form important prey for the most common fish larvae found in the ETP off Mexico (ETPM), Bregmaceros bathymaster (Siordia-Cermeno, Sanchez-Velasco, Sanchez-Ramirez, & Franco-Gordo, 2006; Davies et al., 2015). This permits the establishment of complex food webs with a seasonal component and ecological importance for a region considered to be oligotrophic, but which supports a variety of artisanal fisheries (Rojo-Vazquez et al., 2008).

Previous studies about the variability of copepod communities in the ETPM were either carried out exclusively in coastal sites (Kozak, Franco-Gordo, Suarez-Morales, & Palomares-Garcia, 2014; Jimenez-Perez, 2016) or over a greater spatial range but restricted to a single period without seasonal replicates (Chen, 1986; Lopez-Ibarra, Hernandez-Trujillo, Bode, & Zetina-Rejon, 2014). We hypothesize that variability in the abundance and diversity of the copepod community was principally driven by seasonal mesoscale activity (upwelling and eddies) and the moderate El Nino-La Nina event, which impacted the tropical Pacific during this period. The goal of this study was to provide the first analysis of seasonal changes in the copepod community structure of the ETPM through the El Nino-La Nina transition of 2010, and to examine the relationship between the distribution of copepods in the ETPM and environmental variables.


Sampling: Three oceanographic cruises were carried out from 15-27 January (Jan), 25 May- 4 June (May), and 18-29 October (Oct) 2010. Twenty-one stations were sampled in each cruise distributed along three transects parallel to the coast at a distance of 10 nm, 50 nm and 100 nm from Cabo Corrientes, Jalisco to Maruata, Michoacan, Mexico (Fig. 1). Six samples were unavailable for analysis (Jan 3, 4, EX2; Mar 9, 13, 14; Oct 11, 15, EX1), leaving 18 stations for each of the three sampled months. Sampling methods are as in Pelayo-Martinez et al. (2017). Briefly, oblique zooplankton tows were performed during daytime and nighttime at each station, from 200 m depth or as deep as possible, using a bongo net with a mesh size of 505 pm, and with a mouth opening of 0.6 m in diameter (Smith & Richardson, 1977). A calibrated flowmeter was attached to the mouth of the net to estimate the amount of water filtered. Samples were fixed in a 4% formalin solution buffered with sodium borate (Griffiths, Fleminger, Kimor, & Vannucci, 1976). Before zooplankton sampling, temperature and salinity profiles were recorded at each station using a CTD (SBE-19 plus), with the depth of each cast ranging from 75 to 500 m depending on the bathymetry of each station. Water samples were collected using Niskin bottles at depths of 200, 150, 100, 75, 50, 25, 10 and 0 m and analyzed for inorganic nutrients (N[O.sub.3.sup.-] + N[O.sub.2.sup.-], N[H.sub.4.sup.+], P[O.sub.4.sup.-3], Si[(OH).sub.4]) using an auto-analyzer Skalar San Plus II and chlorophyll a (Chl a) with a Perking Elmer's UV/Vis spectrophotometer following the techniques proposed by Grasshoff, Kremling, and Ehrhardt (1983) and Strickland and Parsons (1972). In accordance with Talley, Pickard, Emery, and Swift (2011), altimetry data from the AVISO program ( were processed to identify sea level anomalies and geostrophic currents (with 1/3[degrees] resolution); both parameters were plotted using [R]MATLAB software.

Species identification: In the laboratory, samples were fractioned using a Folsom plankton splitter to obtain final sub-samples of 300500 specimens (Harris, Wiebe, Lenz, Skjoldal, & Huntley, 2000). Adult copepods were identified to species level and copepodites were identified to family level. Although probably under-sampled as a group, cyclopoid copepods were still identified to species level (except in the case of Oncaea spp., which due to uncertainty were identified at the genus level only). We included these organisms as they are part of the total biomass of the larger size fraction of copepods. Identifications followed the keys of Palomares, Suarez, and Hernandez-Trujillo (1998) and Razouls, de Bovee, Kouwenberg, and Desreumaux (2005-2017); when necessary other specialized keys were used.

Data analysis: Standardized copepod abundance data were explored to detect outliers, then [Log.sub.10] [x+1] transformed before performing statistical or multivariate analyses to decrease the variance of the data set and to avoid violating assumptions of normality, except for non-parametric analyses which do not require such transformations. Because zooplankton samples were taken during both night and day, an analysis of variance (ANOVA) was first used to check for significant differences between total copepod abundance and time of sampling. Daytime and nighttime classifications were based on the time of sunrise and sunset during each sampling period. Daytime = 1 hr before sunrise to 1 hr after sunset. Nighttime = 1 hr after sunset to 1 hr before sunrise. No significant differences were found, so samples were defined a priori to statistical analysis by month and inshore (0, 4, 8, 12, 17, 19, EX1) or offshore (remaining stations). Mean abundance and abiotic data of the different sampling periods Jan, May and Oct) were compared using an ANOVA to determine overall changes in the region. If significant results were detected, a Tukey's test was used post hoc to define the differences. The ANOVA and post hoc analyses were calculated using Statistica 10.0 software (Stat Soft).

A cluster analysis (CA) was performed to define copepod assemblages using the Sorensen distance measure and [beta]-flexible linkage method ([beta] = 0.25), followed by a multi-response permutation procedure (MRPP) to test the significance of the resulting groups. The similarity percentage (SIMPER) was used to identify the species that account for 90% of abundance of the species/families that made up the groups defined by the CA. An indicator value analysis (IVA) was also applied to determine which species could serve as indicators of certain hydrological conditions based on the CA groups (Dufrene & Legendre, 1997). The IVA provides the percentage of indication by combining values for relative abundance and relative frequency, with results ranging from zero (no indication) to 100 (perfect indication). Canonical correspondence analysis (CCA) defined the relationships between species assemblages and environmental characteristics in multivariate space. The environmental (explanatory) matrix in the CCA consisted of surface (0-50 m) and subsurface (75-200 m) means of temperature, salinity, Chl a (mg [m.sup.-3]), dissolved oxygen (DO; [micro]mol [kg.sup.-1]), nitrate+nitrite ([micro]M), ammonium ([micro]M); phosphate (P[O.sup.3-.sub.4], [micro]M), and silica ([micro]M). The biological (response) matrix consisted of 76 species and copepodite taxa. To test for the significance of the first ordination axis and all canonical axes together, a Monte Carlo permutation test (499 permutations) was applied. Rare species, defined as those appearing in < 5% of samples, were not included in either CA or CCA analyses. Multivariate analyses were carried out using PC-ORD v.4 software (McCune & Mefford, 1999).

Temperature/Salinity profiles of the water masses found in the study area during the three oceanographic cruises were created by integrated plots using [R]MATLAB 7.11 software according to the limits proposed by Portela et al. (2016). Based on the environmental variables that explained most of the variability seen in the sample distribution, contour maps representing the horizontal distribution of Chl a (0-50 m), salinity (0-50 m), N[O.sub.3.sup.-] + N[O.sub.2.sup.-] (0-50 m), and N[H.sub.4.sup.+] (75-200 m) were created using Surfer 11 software.


Environmental conditions: The survey period of the present study included the second half of the 2009-2010 El Nino (EN, Jan 2010), the neutral transition period (May 2010), and the first half of the 2010-2011 La Nina (LN, Oct 2010). [The entire ETP experienced an ENSO cycle from 2009-2011 (http://www.cpc. ensostuff/ensoyears.shtml)]. The water masses in the upper water column during the three periods were primarily characterized by Tropical Surface Water (TSW), Gulf of California Water (GCW) and a high proportion of transitional water as a result of the confluence of those water masses in the region. In the intermediate layer, the water column was dominated by Subtropical Subsurface Water (StSsW). Jan samples were mostly TSW, probably due to the effect of EN. May samples were dominated by transitional waters and a slight influence of GCW. Oct samples were more varied, with a stronger influence of GCW and a slight influx of California Current Water (CCW) (Fig. 2). The mean temperature of the surface layer (0-50 m) was warmest during Jan (26.8 [degrees]C), which coincided with the expanded mixed layer depth (MLD 45.8 m), the highest value of the Multivariate Enso Index (MEI 1.52) and the lowest value of the Coastal Upwelling Index (CUI; products/PFEL/modeled/indices/upwelling/ NA/data_download.html).

Eddy processes: Geostrophic currents showed seasonal variations in direction and intensity, along with a subsequent increase-decrease of temperatures (Fig. 3A, Fig. 3B and Fig. 3C). Jan was dominated by a geostrophic flow from the south. A very weak cyclonic eddy was apparent near the center of the surveyed area (Fig. 3A), but it had little effect on surface temperatures, which were relatively high and uniform across the region in the upper water column (0-50 m). A cyclonic eddy in May showed a much stronger flow (Fig. 3B) and resulted in cooler temperatures at the core. During this period a geostrophic current parallel to the continent was apparent; it is probable that the effect of movement near the coast along with the topography of the platform generated internal waves. This could have generated a point of lifting of the water column, resulting in lower temperatures and brief episodes of primary productivity in the surface waters near the coast. In Oct a large anticyclonic (warm core) eddy was apparent in the central study area, with a much smaller cold core eddy directly below it (Fig. 3C).

Upwelling and Chl a: Over the entire study area, nutrient concentrations at surface (0-50 m) and sub-surface (75-200 m) depths were averaged and analyzed for differences between surveyed periods (Table 1), and several of the most relevant variables were expressed as contour maps (Fig. 3). Both 0-50 m and 75-200 m ammonium values were highest in Jan. Values of 0-50 m nitrates+nitrites (mineralized nitrogen found in deeper waters), silicates and phosphates were significantly higher in May, representative of upwelling processes along the coast as well as in the large cyclonic eddy (Fig. 3B and Fig. 3K). Jan showed very low nitrate+nitrite values over the entire region, except for a slight increase near the weak cyclonic eddy (Fig. 3A and Fig. 3J). Downwelling conditions dominated during Jan, as shown by significantly higher temperatures, ammonium and DO concentrations, and a deeper MLD than the other two periods (Table 1). The monthly mean of daily CUI values corroborated that coastal upwelling processes were highest in May, were present during Oct, and absent during Jan (Table 1). Despite the lack of active upwelling conditions, very high 0-50 m Chl a concentrations were observed at an inshore station during Jan (Fig. 3D), resulting in significantly higher concentrations than in May and Oct (Table 1).

Copepod abundance and diversity:

A total of 78 copepod species belonging to 19 families and 2 orders were identified (Appendix). Total copepod abundance per station varied from 3 798 to 52 841 ind./1 000 [m.sup.3] across the three sampling periods. The dominant order was Calanoida, whereas the relative abundance of Cyclopoida was much lower. Among calanoid copepods, Eucalanidae, Euchaetidae, and Calanidae were the dominant families. The dominant families of Cyclopoida were Corycaeidae, and Oithonidae. There was no significant difference in the total abundance of copepods between periods (Table 2). The Cyclopoida showed significant variations between months, with its highest abundances during Jan. Species richness and diversity showed the same seasonal and spatial patterns. Richness ranged from 11 to 47 species per station, with seasonal averages from 25 species (May) to 35 species (Jan). Average richness and diversity was significantly lower in May than in the other two seasons (Table 2). Contrary to what was observed in Jan (species richness and diversity were high and evenly distributed) and Oct (species richness and diversity showed latitudinal differences), species richness and diversity showed an inshore-offshore gradient in May (ANOVA, P > 0.05). (Fig. 4D, Fig. 4E and Fig. 4F).

Most adult and copepodite taxa with mean densities of >0.5 ind./1 000 [m.sup.3] across all stations presented significant differences between surveyed periods. Post hoc Tukey tests showed that the two most abundant groups had reverse seasonal patterns; copepodites of the family Euchaetidae were significantly more abundant in Jan than in the other two periods, while copepodites of the family Eucalanidae were significantly more abundant in May and Oct, explaining the lack of significant differences in total abundance between periods (Table 2). After copepodites of the families Eucalanidae and Euchaetidae, the most abundant species were Subeucalanus subcrassus and Euchaeta indica. Distribution plots of the total abundance showed slightly higher abundances along the coast during Jan (Fig. 4A). In May the most coastal inshore station (station 17) had extremely high abundances, but overall abundances were lower along the coast and in the region of the cyclonic eddy than in the Northern oceanic region (Fig. 4B). Abundance was slightly more evenly distributed in Oct, but did have higher values in the same inshore region as observed in May (Fig. 4A, Fig. 4B and Fig. 4C). Distribution plots of S. subcrassus and E. indica visually demonstrate their opposite abundances between periods. Moreover, the highest abundances of E. indica in Jan were in the same region as the extremely high Chl a levels (Fig. 4G, Fig. 4H and Fig. 4I).

Copepod distribution and species groups: The CA indicated that there were four principal assemblages present across the surveyed periods (Fig. 5A). The Jan EN cluster (G1) was comprised of nearly all the stations sampled in Jan. There was a small sub-cluster, which included the coastal stations in this group, but otherwise there was no relevant grouping pattern. The May offshore cluster (G2) was totally comprised of stations from the May survey period located 100 nm from the coast, but did not include the stations closest to the cyclonic eddy. The Oct LN cluster (G3) was principally comprised of stations from the Oct survey period. The May coastal and eddy upwelling cluster (G4) was composed of May stations affected by coastal upwelling or associated with the cold-core cyclonic eddy, but also included stations from Oct which were located in areas of coastal upwelling and the small cyclonic eddy in the Southern region during that period. The MRPP confirmed that the groups defined by the CA were significantly different (P < 0.05). ANOVA indicated that there were no significant differences in total abundance across the cluster groups (P > 0.05), but species richness and diversity were significantly lower in May upwelling than all other groups (P < 0.05). SIMPER analysis showed that Eucalanidae and Euchaetidae were the most abundant families, with copepodites contributing the most abundance all four of the cluster groups (Table 3). However, the relative abundance varied between periods. Adults and copepodites of Euchaetidae contributed the highest percentage of abundance to the relatively more diverse Jan EN (26.76%) and May offshore (21.21%) clusters, while adults and copepodites of Eucalanidae contributed the highest abundance to the Oct LN (25.91%) and May upwelling (41.83%) clusters. Copepodites of Euchaetidae dropped to their lowest relative contribution during May upwelling. The IVA resulted in 9 significant indicators for Jan EN, 8 for May offshore, 3 for Oct LN, and only 1 for May upwelling (Table 4). Pontellina sobrina had the highest percentage IVA of all species, with 65.3 for Jan EN, followed by Paraucalanus sewelli with 59.3 for May offshore.

Association with environmental variables: The seasonal and interannual distribution of environmental variables and their effect on the copepod community were analyzed using CCA, with stations labeled a priori as inshore/offshore and by month, with CA groups shown as dotted lines connecting the stations (Fig. 5B). The MRPP confirmed that the inshore/offshore groups were significantly different (P < 0.05). The first two axes explained 21.5% of the variance seen in species data, with a Pearson Species-Environment Correlation of 0.930. Axis 1 explained the greatest amount of variation, primarily through surface (0-50 m) temperature, salinity, phosphorus, nitrates+nitrites, and sub-surface (75-200 m) ammonium. Axis 2 was most closely associated with dissolved oxygen. The CCA biplot showed that the groups had a well-defined distribution associated with environmental conditions, although with the presence of some slight overlap (Fig. 4A). Stations from the G4 upwelling cluster were distributed along the upper right side of the graph, showing correlations with lower surface temperature, higher surface nitrites+nitrates and salinity, indicative of upwelling conditions. G3 Oct LN stations were primarily grouped in the center of the biplot, showing the modification of what is usually the warmest and least productive period in the region to slightly cooler, more productive conditions in the water column. The G2 May offshore stations did not show a high association with the environmental variables, indicating that other processes affected the species composition of those groups. G1 EN cluster stations were primarily associated with high temperatures and subsurface ammonium.


Tropical copepods form a large part of zooplankton biomass and are therefore important drivers of zooplankton productivity, but the impact of physical and biological variables on copepod community structure is poorly known in the ETPM. The analysis of the effect of hydrological changes such as increased temperature on species distribution and abundance can provide a baseline to understand the effects of short-term variability such as ENSO and long-term hydrological modifications driven by climate change on copepod community structure. The present study includes three seasonally distinct survey periods -Jan (transitional), May (tropical CC and coastal upwelling), and Oct (MCC and coastal downwelling) that were also distinguished by EN (Jan), ENSO neutral (May) and LN (Oct) conditions- over the course of one year. We found that during the EN sampling period, the structure and abundance of the copepod community was affected by the modification of regional hydrodynamics, enhancing the intrusion of carnivorous species and causing a decrease in total abundance which was maintained throughout the ENSO neutral and Oct LN periods.

The contrasting effects of EN and LN on temperature and water column stability appeared to modify the established pattern of currents and coastal upwelling during the study period, processes which have been reported to have strong effects on zooplankton and ichthyoplankton communities in both oceanic and coastal regions of the ETPM (Ambriz-Arreola, Gomez-Gutierrez, Franco-Gordo, Lavaniegos, & Godinez-Dominguez 2012; Leon-Chavez, Beier, Sanchez-Velasco, Barton, & Godinez, 2015). Portela et al. (2016) discussed water masses and circulation in the central Mexican Pacific; the geostrophic flows reported in their study varied from those seen in 2010, probably due to the El Nino (Jan) bringing TSW across the entire region and La Nina (Oct) disrupting the northward flowing MCC and inducing an anomalous intrusion of sub-surface CC water -the oceanic, tropical extension of this current is usually strongest in spring (Kessler, 2006). Two specimens of Calanus pacificus, a species of temperate affinity (Palomares et al., 1998), were recorded in the Northern offshore transect during Oct 2010, corroborating our reported intrusion of CCW. To the best of our knowledge, this is the first report of this CC species indicator, which has been reported along the Baja California Peninsula and the Gulf of California, in the central Mexican Pacific (Brinton, Fleminger, & Siegel-Causey, 1986; Jimenez-Perez & Lavaniegos, 2004). Although Lopez-Ibarra et al. (2014) only included species with abundances > 80% in their study along the entire Mexican Pacific, and it is possible that some few unreported specimens were also present further South. Based on historical CUI values and on Lopez-Sandoval et al. (2009), coastal upwelling and primary productivity, coastal upwelling is usually low in Jan (transitional), relatively high in May, and relatively low to absent in Oct. This pattern was apparently modified by both EN and LN, as Chl a concentrations were very high in Jan (EN), and the CUI was positive for Oct (82.4), with patches of relatively cool surface waters due to extensive coastal upwelling and a small cyclonic eddy to the south, probably resulting from the thermocline-lifting effect of the LN. Kurczyn et al. (2012) characterized mesoscale activity in this region (Cabo Corrientes in that study, corresponding to the ETPM) and found relatively high levels of both anticyclonic and cyclonic eddies, which were associated with seasonal peaks in currents and CUI and no statistically significant association with EN events. Due to the short duration of the present study, it is not possible to confirm a definite effect of EN on eddy production; however, considering that both coastal upwelling events and current processes were modified, it is very possible that an effect occurred on eddy generation, subsequently affecting the ecology of the study area.

There were no significant differences in copepod abundances between the survey periods of the present study. Kozak et al. (2014) studied adult calanoid copepods from the inshore region of the present study area (19[degrees] N, 104[degrees] W), using the same methodology (oblique hauls using a 505 pm mesh net), and reported abundances between 50-200 ind./[m.sup.3] per month, with the lowest values in the EN and post EN periods. In the present study, the inshore abundances of adult calanoid copepods per month averaged 7.2 ind./[m.sup.3] (Jan), 7.7 ind./[m.sup.3] (May) and 4.4 ind./[m.sup.3] (Oct), with a maximum value of 34.3 ind./[m.sup.3] during May (station 17). Therefore, we hypothesize that the lack of significant differences in abundance between the three periods resulted from the EN event at the start of our study period, which depressed this larger size fraction of copepod abundance through at least Oct. Jimenez-Perez (2016) also reported significantly lower copepod abundances during EN in Bahia de Banderas (Nayarit, Mexico), demonstrating that this is a region-wide response to EN events. There is still a much uncertainty regarding the impact of climate change on ENSO (Collins et al., 2010). However, this study indicates that if EN anomalies increase in frequency and/or strength due to global climate change, the rapid succession of disturbances without a sufficient recovery period could have marked repercussions in the overall zooplankton abundance and community resilience, impacting marine ecosystem productivity.

The species composition in this study is very similar to that reported in Kozak et al. (2014); there were no significant alterations in species composition between surveys and most of the same species were present across the different hydrological periods. However, the relative abundances of species did vary between surveys in both studies. During the Jan EN period of this survey, the influx of TSW across the entire studied region induced high species diversity and no inshore-offshore gradients of abundance, S or H'. Oct (LN period), which usually presents uniformly warm temperatures across the region due to the influence of the MCC and relatively high inshore species diversity (Kozak et al. 2014), showed comparatively low S and H'. During May, which was ENSO neutral, a cyclonic mesoscale eddy, coastal upwelling, and an inshore-offshore temperature gradient were re-established, along with lower values of S and H'. Recent studies in the ETPM have shown the important effects of mesoscale processes (upwelling, eddies) on zooplankton distribution and structure (Leon-Chavez, et al., 2010; Davies et al., 2015). Cruz-Hernandez et al. (2018) reported a significant vertical stratification of calanoid copepod species in a mature cyclonic eddy in the Gulf of California, and found that herbivorous/ omnivorous species were dominant. We did not find significant correlations with total abundance, however, species richness in the May upwelling cluster (cyclonic eddy and coastal stations) was significantly lower than all other cluster-defined groups (P < 0.05) with larger, generally herbivorous copepodites of Eucalanidae and Subeucalanus subcrassus dominating the assemblage. Even though the abundances of adults and copepodites of Euchaetidae were low in the upwelling assemblage, in the May offshore assemblage they maintained high levels. Eucalanidae (primarily Subeucalanus subcrassus) may dominate during the more productive periods, but Euchaetidae appears to maintain high abundance values of adults and copepodites in the warm offshore regions that are not affected by cyclonic eddies and which can invade the coast rapidly when conditions permit. This appears to have occurred in Jan, when there was no upwelling grouping despite what appeared to be a mature cyclonic eddy in the oceanic region; the EN appeared to homogenize the offshore/inshore/eddy species distribution and the relatively high temperatures precipitated dominance of Euchaetidae.

Despite upwelling events in May and Oct with high nutrient values and low temperatures that should correlate with high primary productivity, 0-50 m Chl a values were highest in Jan, especially in one coastal station and in the weak cyclonic eddy. High ammonium values could be from increased zooplankton metabolism (due to high temperatures) or microbial loop activity, both of which primarily excrete nitrogen as ammonium (Corner & Newell, 1967). This could have provided some of the nutrients required for the observed high phytoplankton production (Alcaraz, Saiz, & Estrada, 1994). Pelayo-Martinez et al. (2017) hypothesize that the values seen in this same study period could be due to removal of the top down control by herbivorous zooplankters on phytoplankton growth by carnivorous zooplankters. Our study aligns with these hypotheses, as we found that oceanic, carnivorous copepod species (primarily Euchaeta adults and Euchaetidae copepodites) were most abundant during the Jan EN period across the entire region. This also has important implications for the trophic structure of the region. Due to the relatively small size of these organisms in comparison to other carnivorous zooplankters (e.g. chaetognaths), they prey on much smaller organisms, probably eggs and early larval stages of planktonic organisms. Therefore, the transfer of carbon biomass could be reduced, as well as the overall productivity of the ecosystem. In addition, as warming trends continue, not only can tropical copepod size be reduced due to increased metabolism, but there could be a regime shift to smaller species. A "tropicalization" of the copepod community during the 1997/98 EN was reported off the coast of Baja California by Jimenez-Perez and Lavaniegos (2004), and it is generally accepted that continued warming events will contribute to that anomaly, with smaller tropical copepod species expanding their distribution into more temperate regions. However, the indicator species assemblage (highest percentage IVA) for the Jan EN period was primarily composed of relatively smaller (< 1 100 [micro]m) species than most of the May and October surveys significant indicators (> 1 500 [micro]m). This suggests that even what has been a fairly stable tropical assemblage is not immune to the biogeographic redistribution of species that has been observed in higher latitudes (Beaugrand, Reid, Ibanez, Lindley, & Edwards, 2002).

This study was limited by the short sampling duration (< 1 year) that did not include non-anomalous seasons for comparison. In addition, there were certainly species (especially cyclopoid copepods) that were under-sampled due to the relatively large-sized mesh (505 [micro]m) of the sampling net. A smaller mesh size would have been ideal for a more complete analysis of the copepod community. Regardless, the results of this study are still valuable due to the proportionally high abundance of these larger copepods in this marine ecosystem; Pelayo-Martinez et al. (2017) reported that the copepods from this study accounted for 71% of total zooplankton abundance. As the dominant organisms in this larger size fraction, they play an important role in overall zooplankton productivity and trophic dynamics, and variations in their abundance could have a large impact on the ecology of the region.

In summary, the species distribution of copepods in the central Mexican Pacific appears to be principally defined by seasonal modifications of the water column as well as mesoscale eddy processes in the oceanic region and by upwelling episodes in the inshore region. The mesoscale processes are affected by ENSO events, which in the case of EN seem to disrupt the inshore -offshore gradient, import carnivorous, offshore species into the inshore area, cause marked reductions in copepod abundance across the entire ETPM, and increase dominance of smaller sized species. Future studies in the region should incorporate smaller mesh sized nets, stratified samples, and a longer time series to study three-dimensional variability of copepods in a transitional oceanographic region that is markedly affected by interannual anomalies.


The Comision Federal de Electricidad (CFE) and Universidad de Colima supported this study through the project "Marine mammals and turtles monitoring in Mexican Central Pacific" during 2010. We thank the crew, students and researchers who assisted during the oceanographic surveys on board the BIPXII vessel. E.R.K., A.O.O. and C.F.G. are SNI fellows and G.P.M. is a doctoral student (UdG) supported by a CONACyT grant #329882. We also thank several anonymous reviewers for their comments that greatly improved an earlier version of this manuscript.


Alcaraz, M., Saiz, E., & Estrada, M., (1994). Excretion of ammonia by zooplankton and its potential contribution to nitrogen requirements for primary production in the Catalan Sea (NW Mediterranean). Marine Biology, 119(1), 69-76.

Ambriz-Arreola, I., Gomez-Gutierrez, J., del Carmen Franco-Gordo, M., Lavaniegos, B. E., & Godinez-Dominguez, E. (2012). Influence of coastal upwelling-downwelling variability on tropical euphausiid abundance and community structure in the inshore Mexican central Pacific. Marine Ecology Progress Series, 451, 119-136.

Ambriz-Arreola, I., Gomez-Gutierrez, J., Franco-Gordo, C., & Kozak, E. R. (2015). Reproductive biology, embryo and early larval morphology, and development rates of krill (Euphausia lamelligera and Euphausia distinguenda), endemic to the Eastern Tropical Pacific. Sexuality and Early Development in Aquatic Organisms, 1(2), 143-161.

Beaugrand, G., Reid, P. C., Ibanez, F., Lindley, J. A., & Edwards, M., (2002). Reorganization of North Atlantic marine copepod biodiversity and climate. Science, 296(5573), 1692-1694.

Brinton, E., Fleminger, A., & Siegel-Causey, D. S. (1986). The temperate and tropical planktonic biotas of the Gulf of California. California Cooperative Oceanic Fisheries Investigations Report, 27, 228-266.

Cepeda-Morales, J., Beier, E., Gaxiola-Castro, G., Lavin, M. F., & Godinez, V. M. (2009). Effect of the oxygen minimum zone on the second chlorophyll maximum in the Eastern Tropical Pacific off Mexico. Ciencias Marinas, 35(4), 389-403.

Chen, Y. Q. (1986). The vertical distribution of some pelagic copepods in the eastern tropical Pacific. California Cooperative Oceanic Fisheries Investigations Report, 27, 205-227.

Collins, M., An, S. I., Cai, W., Ganachaud, A., Guilyardi, E., Jin, F. F., ... Vecchi, G. (2010). The impact of global warming on the tropical Pacific Ocean and El Nino. Nature Geoscience, 3(6), 391.

Corner, E. D. S., & Newell, B. S. (1967). On the nutrition and metabolism of zooplankton IV. The forms of nitrogen excreted by Calanus. Journal of the Marine Biological Association of the United Kingdom, 47(1), 113-120.

Cruz-Hernandez, J., Sanchez-Velasco, L., Godinez, V M., Beier, E., Palomares-Garcia, J. R., Barton, E. D., & Santamaria-Del-Angel, E. (2018). Vertical distribution of calanoid copepods in a mature cyclonic eddy in the Gulf of California. Crustaceana, 91(1), 63-84.

Davies, S. M., Sanchez-Velasco, L., Beier, E., Godinez, V. M., Barton, E. D., & Tamayo, A. (2015). Three-dimensional distribution of larval fish habitats in the shallow oxygen minimum zone in the eastern tropical Pacific Ocean off Mexico. Deep Sea Research Part I: Oceanographic Research Papers, 101, 118-129.

Dufrene, M., & Legendre, P. (1997). Species assemblages and indicator species: the need for a flexible asymmetrical approach. Ecological Monographs, 67, 345-366.

Franco-Gordo, C., Ambriz-Arreola, I., Kozak, E. R., Gomez-Gutierrez, J., Plascencia-Palomera, V., Godinez-Dominguez, E., & Hinojosa-Larios, A. (2015). Seasonal succession of zooplankton taxonomic group assemblages in surface waters of Bahia de Navidad, Mexico (November 2010-December 2011). Hidrobiologica, 25(3), 335-345.

Godinez, V. M., Beier, E., Lavin, M. F., & Kurczyn, J. A. (2010). Circulation at the entrance of the Gulf of California from satellite altimeter and hydrographic observations. Journal of Geophysical Research: Oceans, 115(C4).

Gomez-Valdivia, F., Pares-Sierra, A., & Flores-Morales, A. L. (2015). The Mexican Coastal Current: A subsurface seasonal bridge that connects the tropical and subtropical Northeastern Pacific. Continental Shelf Research, 110, 100-107.

Grasshoff, K., Kremling, K., & Ehrhardt, M. (1983). Methods of seawater analysis. Weinheim, Germany: Verlag Chemie.

Griffiths, F. B., Fleminger, A., Kimor, B., & Vannucci, M. (1976). Shipboard and curating techiniques. In H. F. Steedman (Ed.), Zooplankton fixation and preservation--monographs on oceanographics methodology (pp. 17-19). Paris, France: UNESCO.

Harris, R., Wiebe, P, Lenz, J., Skjoldal, H. R., & Huntley, M. (Eds.). (2000). ICES zooplankton methodology manual. Jamestown Road, London: Academic Press.

Jimenez-Perez, L. C. (2016). Copepod community structure in Bahia de Banderas during the 2008-2009 La Nina and their transition to the 2009-2010 El Nino. Revista Bio Ciencias, 4(2), 82-103.

Jimenez-Perez, L. C., & Lavaniegos, B. E. (2004). Changes in dominance of copepods off Baja California during the 1997-1999 El Nino and La Nina. Marine Ecology Progress Series, 277, 147-165.

Kessler, W. S., (2006). The circulation of the eastern tropical Pacific: A review. Progress in Oceanography, 69(2), 181-217.

Kozak, E. R., Franco-Gordo, C., Palomares-Garcia, R., Gomez-Gutierrez, J., & Suarez-Morales, E. (2017). Annual egg production rates of calanoid copepod species on the continental shelf of the Eastern Tropical Pacific off Mexico. Estuarine, Coastal and Shelf Science, 184, 138-150.

Kozak, E. R., Franco-Gordo, C., Suarez-Morales, E., & Palomares-Garcia, R. (2014). Seasonal and interannual variability of the calanoid copepod community structure in shelf waters of the Eastern Tropical Pacific. Marine Ecology Progress Series, 507, 95-110.

Kurczyn, J. A., Beier, E., Lavin, M. F., & Chaigneau, A. (2012). Mesoscale eddies in the northeastern Pacific tropical-subtropical transition zone: Statistical characterization from satellite altimetry. Journal of Geophysical Research: Oceans, 117(C10).

Lavin, M. F., Beier, E., Gomez-Valdes, J., Godinez, V. M., & Garcia, J. (2006). On the summer poleward coastal current off SW Mexico. Geophysical Research Letters, 33(2), 1-4.

Leon-Chavez, C. A., Beier, E., Sanchez-Velasco, L., Barton, E. D., & Godinez, V. M. (2015). Role of circulation scales and water mass distributions on larval fish habitats in the Eastern Tropical Pacific off Mexico. Journal of Geophysical Research: Oceans, 120(6), 3987-4002.

Leon-Chavez, C. A., Sanchez-Velasco, L., Beier, E., Lavin, M. F., Godinez, V. M., & Farber-Lorda, J. (2010). Larval fish assemblages and circulation in the Eastern Tropical Pacific in Autumn and Winter. Journal of Plankton Research, 32(4), 397-410.

Lopez-Ibarra, G. A., Hernandez-Trujillo, S., Bode, A., & Zetina-Rejon, M. J. (2014). Community structure of pelagic copepods in the eastern tropical Pacific Ocean during summer and autumn. Cahiers de Biologie Marine, 55, 453-462.

Lopez-Sandoval, D. C., Lara-Lara, J. R., Lavin, M. F., Alvarez-Borrego, S., & Gaxiola-Castro, G. (2009). Primary productivity in the eastern tropical Pacific off Cabo Corrientes, Mexico. Ciencias Marinas, 35(2), 169-182.

McCune, B., & Mefford, M. J. (1999). PC-Ord. Multivariate Analysis of Ecological Data (Version 5.0). Gleneden Beach, Oregon, USA: MjM Software Design.

Palomares, R., Suarez, E., & Hernandez-Trujillo, S. (1998). Catalogo de los copepodos (Crustacea) pelagicos del PacificoMexicano. Mexico: CICIMAR/ECOSUR.

Pelayo-Martinez, G., Olivos-Ortiz, A., Franco-Gordo, C., Quijano-Scheggia, S., Gavino-Rodriguez, J., Kono-Martinez, T., & Castro-Ochoa, F. (2017). Physical, chemical and zooplankton biomass variability (inshore-offshore) of Mexican Central Pacific during El Nino-La Nina 2010. Latin American Journal of Aquatic Research, 45(1), 67.

Portela, E., Beier, E., Barton, E. D., Castro, R., Godinez, V., Palacios-Hernandez, E., ... Trasvina, A. (2016). Water Masses and Circulation in the Tropical Pacific off Central Mexico and Surrounding Areas. Journal of Physical Oceanography, 46(10), 3069-3081.

Razouls, C., De Bovee, F., Kouwenberg, J., & Desreumaux, N. (2005-2017). Diversity and geographic distribution of marine planktonic copepods. Retrieved from

Rojo-Vazquez, J. A., Quinonez-Velazquez, C., Echavarria-Heras, H. A., Lucano-Ramirez, G., Godinez-Dominguez, E., Ruiz-Ramirez, S., ... Sosa-Nishizaki, O. (2008). The fish species composition and variation of catch from the small-scale gillnet fishery before, during and after the 1997-1998 ENSO event, central Mexican Pacific. Revista de Biologia Tropical, 56(1), 133-152.

Salas, J., Gomis, D., Olivos-Ortiz, A., & Garcia-Uribe, G. (2006). Seasonal hydrodinamical features on the continental shelf of Colima (west coast of Mexico). Scientia Marina, 70(4), 719-26.

Siordia-Cermeno, M. P., Sanchez-Velasco, L., Sanchez-Ramirez, M., & Franco-Gordo, M. C. (2006). Variacion temporal de la dieta de larvas de Bregmaceros bathymaster (Pisces: Bregmacerotidae) en las costas de Jalisco y Colima, Mexico, durante un ciclo anual (1996). Ciencias Marinas, 32(1a), 13-21.

Smith, P., & Richardson, S. (1977). Standard techniques for pelagic fish egg and larva surveys. FAO Fisheries Technical Paper, 175, 1-100.

Strickland, J. D., & Parsons, T. R. (1972). A practical handbook for the sea water analysis. Ottawa, Canada: Bulletin Fisheries Research Board of Canada.

Talley, L. D., Pickard, G. L., Emery, W. J., & Swift, J. H. (2011). Descriptive Physical Oceanography: An Introduction. Boston, Massachusetts: Academic Press.

Torres-Orozco, E., Trasvina, A., Muhlia-Melo, A., & Ortega-Garcia, S. (2005). Mesoscale dynamics and yellowfin tuna catches in the Mexican Pacific. Ciencias Marinas, 31(4), 671-683.

Zamudio, L., Hurlburt, H. E., Metzger, E. J., & Tilburg, C. E. (2007). Tropical wave-induced oceanic eddies at Cabo Corrientes and the Maria Islands, Mexico. Journal of Geophysical Research: Oceans, 112(C5), 1-17.

Order CALANOIDA Sars 1902

Family Acartiidae Sars 1903
Acartia danae Giesbrecht, 1889
Acartia negligens Dana, 1849

Family Aetideidae Giesbrecht, 1892
Aetideus bradyi Scott, 1909
Euchirella lisettae Klein, 1989

Family Calanidae Dana, 1846
Calanus pacficus Brodsky, 1948
Canthocalanus pauper (Giesbrecht, 1888)
Cosmocalanus darwini (Lubbock, 1860)
Mesocalanus tenuicornis (Dana, 1849)
Nannocalanus minor (Claus, 1863)
Undinula vulgaris (Dana, 1849)

Family Candaciidae Giesbrecht, 1892
Candacia catula (Giesbrecht, 1889)
Candacia curta (Dana, 1849)
Candacia pachydactyla (Dana, 1849)
Candacia tenuimana (Giesbrecht, 1889)
Candacia truncata (Dana, 1849)
Candacia varicans (Giesbrecht, 1892)

Family Centropagidae Giesbrecht 1892
Centropages calaninus (Dana, 1849)
Centropages furcatus (Dana, 1849)
Centropages gracilis (Dana, 1849)

Family Clausocalanidae Giesbrecht, 1892
Clausocalanus furcatus (Brady, 1883)
Clausocalanus jobei Frost & Fleminger, 1968
Clausocalanus minor Sewell, 1929

Family Eucalanidae Giesbrecht, 1892
Eucalanus inermis Giesbrecht, 1892
Pareucalanus sewelli (Fleminger, 1973)
Subeucalanus crassus (Giesbrecht, 1888)
Subeucalanus subcrassus (Giesbrecht, 1888)
Subeucalanus subtenuis (Giesbrecht, 1888)

Family Euchaetidae Giesbrecht, 1892
Euchaeta indica Wolfenden, 1905
Euchaeta longicornis Giesbrecht, 1888
Euchaeta magniloba Park, 1978
Euchaeta rimana Bradford, 1974

Family Lucicutiidae Sars, 1902
Lucicutia gaussae Grice, 1963
Lucicutia gemina Farran, 1926

Family Metrinidae Sars, 1902
Pleuromamma gracilis (Claus, 1863)
Pleuromamma johnsoni Ferrari & Saltzman, 1998
Pleuromamma piseki Farran, 1929
Pleuromamma robusta (Dahl, 1893)

Family Paracalanidae Giesbrecht, 1892
Acrocalanus andersoni Bowman, 1958
Acrocalanus gracilis Giesbrecht, 1888
Acrocalanus longicornis Giesbrecht, 1888
Calocalanus pavo (Dana, 1849)
Calocalanus plumulosus (Claus, 1863)
Paracalanus aculeatus Giesbrecht, 1888
Paracalanus parvus (Claus, 1863)

Family Pontellidae Dana, 1853
Calanopia minor Scott, 1902
Labidocera acuta (Dana, 1849)
Labidocera detruncata (Dana, 1849)
Pontella agassizi Giesbrecht, 1895
Pontellina sobrina Fleminger & Hulsemann, 1974

Family Rhincalanidae Geletin, 1776
Rhincalanus nasutus Giesbrecht, 1888
Rhincalanus rostrifrons (Dana, 1849)

Family Scolecitrichidae Giesbrecht, 1892
Scolecithricella abyssalis (Giesbrecht, 1888)
Scolecithricella marginata (Giesbrecht, 1888)
Scolecithrix bradyi Giesbrecht, 1888
Scolecithrix danae (Lubbock, 1856)
Scolecitrichopsis ctenopus (Giesbrecht, 1888)

Family Temoridae Giesbrecht, 1893
Temora discaudata Giesbrecht, 1889

Order CYCLOPOIDA Burmeister, 1834

Family Corycaeidae Dana, 1852
Corycaeus crassiusculus Dana, 1849
Corycaeus flaccus Giesbrecht, 1891
Corycaeus furcifer Claus, 1863
Corycaeus pacficus Dahl, 1894
Corycaeus robustus Giesbrecht, 1891
Corycaeus speciosus Dana, 1849

Family Oithonidae Dana, 1853
Oithona plumifera Baird, 1843
Oithona robusta Giesbrecht, 1891
Oithona setigera (Dana, 1849)
Oithona similis- Group Claus, 1866

Family Oncaeidae Giesbrecht, 1893
Family Sapphirinidae Thorell, 1859
Copilia mirabilis Dana, 1849
Sapphirina gemma Dana, 1849
Sapphirina metallina Dana, 1849
Sapphirina nigromaculata Claus, 1863
Sapphirina opalina Dana, 1849
Sapphirina sinuicauda Brady, 1883
Sapphirina stellata Giesbrecht, 1891

Eva R. Kozak (1), Aramis Olivos-Ortiz (2), Carmen Franco-Gordo (1) & Gloria Pelayo-Martinez (1)

(1.) Departamento de Estudios para el Desarrollo Sustentable de Zonas Costeras, Universidad de Guadalajara, Gomez Farias 82, San Patricio Melaque, Jalisco 48980, Mexico;,,

(2.) Centro de Universitario de Investigaciones Oceanologicas, Universidad de Colima, Carretera Manzanillo-Barrade Navidad km 20, Col. El Naranjo, C.P.28860 Manzanillo, Colima, Mexico;

Received 23-V-2018.

Corrected 26-VII-2018.

Accepted 30-VIII-2018.

Caption: Fig. 1. Study area in the central Mexican Pacific during January (Jan), May-June (May), and October (Oct) 2010.

Caption: Fig. 2. Temperature/Salinity profiles of the water masses found in the study area during the three oceanographic cruises of 2010. TSW = Tropical Surface Water; CCW = California Current Water; StSsW = Subtropical Subsurface Water; GCW = Gulf of California Water; PIW = Pacific Intermediate Water.

Caption: Fig. 3. Contour maps of important environmental variables in the region during the three study periods. Months are January (Jan) (A, D, G, J, M); May-June (May) (B, E, H, K, N) and October (Oct) (C, F, I, L, O). Note the different color scales for each variable.

Caption: Fig. 4. Monthly total abundance (A-C), species richness (D-F), Subeucalanus subcrassus abundance (G-I) and Euchaeta indica abundance (J-L). Months are January (Jan) (A, D, G, J); May-June (May) (B, E, H, K) and October (Oct) (C, F, I, L). Note the different color scales for each variable.

Caption: Fig. 5. A. Cluster analysis dendrogram of sample sites resulting in 4 groups: G1= Jan EN cluster; G2 = May offshore cluster; G3 = Oct LN cluster; G4 = May coastal and eddy upwelling cluster. Filled symbols are inshore sites, open symbols are offshore sites. B. Canonical Correspondence Analysis biplot of sample sites and environmental variables; sample symbols were defined a priori as inshore and offshore sites by month. Sample groupings (dotted lines) are results of cluster analysis. Environmental variables (arrows) in multivariate space with > 0.400 [r.sup.2] value are shown: T50 = 0-50 mean temperature ([degrees]C); NH200 = mean 75-200 m ammonium ([micro]M); S50 = 0-50 m mean salinity; N50 = 0-50 m mean nitrate+nitrite ([micro]M); P50 = 0-50 mean phosphorus ([micro]M).
Environmental variable mean for three periods in 2010

Variable [1]                   Jan (SE [2])        May (SE [2])

T 0-50                       26.8     (0.21)     23.9     (0.34)
T 75-200                     15.24    (0.48)     13.83    (0.26)
DO 0-50                      174.4    (2.60)     150.7    (7.03)
DO 75-200                    17.65    (4.84)     7.76     (1.28)
S 0-50                       34.15    (0.05)     34.66    (0.02)
S 75-200                     34.73    (0.03)     34.76    (0.01)
Chl a 0-50                   1.51     (0.43)     0.89     (0.16)
Chl a 75-200                 0.22     (0.04)     0.26     (0.02)
CUI                          -9.62    (10.9)     153.5    (4.37)
NO 0-50                      4.14     (0.58)     12.07    (0.80)
NO 75-200                    17.09    (0.38)     18.69    (0.31)
N[H.sub.4.sup.+] 0-50        3.33     (0.09)     2.85     (0.51)
N[H.sub.4.sup.+] 75-200      3.41     (0.10)     2.11     (0.08)
P[O.sup.3-.sub.4] 0-50       1.00     (0.18)     1.74     (0.14)
P[O.sup.3-.sub.4] 75-200     2.75     (0.35)     2.94     (0.13)
Si[O.sub.2] 0-50              7.8     (0.66)     18.35    (1.87)
Si[O.sub.2] 75-200           22.01    (1.17)     29.87    (1.78)
MLD                          45.8     (2.28)     26.3     (1.86)
MEI                          1.52                -0.33

Variable [1]                   Oct (SE [2])        p          F

T 0-50                       23.5     (0.54)    < 0.001     21.36
T 75-200                     13.09    (0.10)    < 0.001     11.66
DO 0-50                      107.1    (6.34)    < 0.001     36.29
DO 75-200                    1.91     (0.55)     0.001       7.48
S 0-50                       34.39    (0.02)    < 0.001     57.11
S 75-200                     34.77    (0.01)     0.340       1.10
Chl a 0-50                   0.58     (0.08)     0.006       5.56
Chl a 75-200                 0.15     (0.03)     0.046       3.23
CUI                          82.4     (6.32)    < 0.001     81.32
NO 0-50                      9.68     (0.90)    < 0.001     27.61
NO 75-200                    19.15    (0.23)    < 0.001     11.84
N[H.sub.4.sup.+] 0-50        2.22     (0.11)     0.042       3.33
N[H.sub.4.sup.+] 75-200      2.39     (0.14)    < 0.001     40.05
P[O.sup.3-.sub.4] 0-50       1.10     (0.10)     0.001       7.97
P[O.sup.3-.sub.4] 75-200     2.37     (0.12)     0.198       1.67
Si[O.sub.2] 0-50             9.36     (1.02)    < 0.001     19.54
Si[O.sub.2] 75-200           22.31    (0.55)    < 0.001     12.28
MLD                          23.5     (1.16)    < 0.001     44.41
MEI                          -1.52                           Not

Variable (1)                post hoc

T 0-50                       J > All
T 75-200                      J >All
DO 0-50                     J > M > O
DO 75-200                    J > Oct
S 0-50                      M > O > J
S 75-200
Chl a 0-50                  Jan > All
Chl a 75-200                  M > O
CUI                         M > O > J
NO 0-50                      All > J
NO 75-200                    All > J
N[H.sub.4.sup.+] 0-50         J > O
N[H.sub.4.sup.+] 75-200     Jan > All
P[O.sup.3-.sub.4] 0-50      May > Oct
P[O.sup.3-.sub.4] 75-200
Si[O.sub.2] 0-50             M > All
Si[O.sub.2] 75-200           M > All
MLD                          J > All

[1.] 0-50 indicates mean values from 0-50 m in water column. 75-200
indicates mean values from 75-200 m in water column. T =
temperature ([degrees]C); DO = dissolved oxygen ([micro]mol
[L.sup.-1]); S= salinity; CUI = daily coastal upwelling index; Chl
a = Chlorophyll a (mg [m.sup.-2]); NO = nitrate+nitrite ([micro]M);
N[H.sub.4.sup.+] = ammonium ([micro]M); P[O.sup.3-.sup.4] =
phosphate ([micro]M); Si[O.sub.2] = silica ([micro]M); MLD = mixed
layer depth (m); MEI = multivariate ENSO index (MEI was not tested
for significance as it is a monthly value).

[2.] SE = standard error.

Mean densities by sampling period of copepod groups and adult and
copepodite taxa with abundance > 0.5 ind./1 000 [m.sup.3]

                                    Jan [1] SE [2]

Total abundance                  9 728.6     (1 200.49)
Calanoida                        9 115.8     (1 139.38)
Cyclopoida                        612.7       (88.98)
Richness                           35.2        (0.98)
Diversity                          2.9         (0.05)
Eucalanidae copepodites           733.6       (161.86)
Euchaetidae copepodites          2 394.3      (364.00)
Subeucalanus subcrassus           284.5       (90.44)
Euchaeta indica                   977.8       (179.79)
Canthocalanus pauper              671.2       (289.20)
Nannocalanus minor                 34.9       (13.31)
Subeucalanus subtenuis            290.0       (80.05)
Centropages furcatus               23.1       (12.83)
Metridinidae copepodites          223.3       (62.32)
Temora discaudata                  61.0       (27.75)
Calanidae copepodites             337.4       (74.40)
Oncaea spp                        296.4       (50.69)
Pleuromamma johnsoni              257.4       (88.46)
Clausocalanus minor               248.4       (53.88)
Euchaeta longicornis              116.5       (14.78)
Cosmocalanus darwini              100.6       (22.47)
Acrocalanus longicornis           221.7       (54.67)
Scolecithrix bradyi                73.8       (12.70)
Euchaeta rimana                   204.3       (37.21)
Mesocalanus tenuicornis            12.8        (6.43)
Temoridae copepodites              13.9        (5.60)
Scolecitrichopsis ctenopus        145.4       (18.24)
Scolecithricella marginata        135.1       (28.35)
Calanopia minor                   134.4       (38.35)
Acrocalanus gracilis              133.6       (36.18)
Aetideus bradyi                    35.2        (7.54)
Undinula vulgaris                 111.8       (35.64)
Candacia truncata                 103.9       (31.00)
Subeucalanus crassus               64.3       (22.83)
Pleuromamma piseki                 30.9        (8.84)
Scolecitrichidae copepodites       77.4        (8.79)
Sapphirina metallina               39.9       (18.02)
Scolecithricella abyssalis         44.8       (12.54)
Aetideidae copepodites             66.5       (15.01)
Acrocalanus andersoni              58.1       (30.66)
Pontellidae copepodites            53.6       (13.89)

                                    May [1] SE [2]

Total abundance                  1 1223.3    (2 695.72)
Calanoida                        1 1003.2    (2 677.32)
Cyclopoida                        220.1       (60.47)
Richness                           25.1        (1.77)
Diversity                          2.3         (0.14)
Eucalanidae copepodites          2 838.5      (924.95)
Euchaetidae copepodites           942.4       (360.79)
Subeucalanus subcrassus          2 353.2     (1 323.05)
Euchaeta indica                   154.9       (52.88)
Canthocalanus pauper              170.2       (51.04)
Nannocalanus minor                536.2       (74.18)
Subeucalanus subtenuis            412.3       (131.25)
Centropages furcatus              382.9       (166.56)
Metridinidae copepodites          399.1       (81.70)
Temora discaudata                 362.4       (157.59)
Calanidae copepodites             206.7       (67.17)
Oncaea spp                         12.8        (6.01)
Pleuromamma johnsoni              266.4       (82.07)
Clausocalanus minor                9.8         (5.03)
Euchaeta longicornis              233.2       (78.40)
Cosmocalanus darwini              230.4       (118.99)
Acrocalanus longicornis            26.9       (18.85)
Scolecithrix bradyi               186.1       (41.36)
Euchaeta rimana                    97.4       (35.40)
Mesocalanus tenuicornis           195.5       (89.37)
Temoridae copepodites             184.8       (66.56)
Scolecitrichopsis ctenopus        144.2       (48.85)
Scolecithricella marginata         3.8         (3.87)
Calanopia minor                    3.1         (2.42)
Acrocalanus gracilis               6.2         (3.87)
Aetideus bradyi                    81.8       (16.83)
Undinula vulgaris                  20.8        (9.46)
Candacia truncata                  60.0       (32.22)
Subeucalanus crassus              103.8       (37.52)
Pleuromamma piseki                 26.8       (12.25)
Scolecitrichidae copepodites       88.1       (24.85)
Sapphirina metallina               87.7       (36.22)
Scolecithricella abyssalis         69.0       (16.17)
Aetideidae copepodites             66.5       (10.87)
Acrocalanus andersoni              60.6       (27.57)
Pontellidae copepodites            67.4       (38.44)

                                     Oct [1] SE [2]

Total abundance                  9 401.9     (1 367.06)
Calanoida                        9 138.5     (1 348.38)
Cyclopoida                        263.4       (55.85)
Richness                           32.1        (1.85)
Diversity                          2.9         (0.10)
Eucalanidae copepodites          3 926.3     (1 039.81)
Euchaetidae copepodites           675.9       (130.47)
Subeucalanus subcrassus           512.2       (138.72)
Euchaeta indica                   149.5       (44.86)
Canthocalanus pauper              144.4       (21.42)
Nannocalanus minor                377.6       (45.88)
Subeucalanus subtenuis            233.5       (44.25)
Centropages furcatus              408.7       (237.85)
Metridinidae copepodites          240.5       (58.60)
Temora discaudata                 216.8       (43.01)
Calanidae copepodites             184.8       (42.04)
Oncaea spp                         40.1       (13.95)
Pleuromamma johnsoni              180.5       (73.67)
Clausocalanus minor                30.7        (9.79)
Euchaeta longicornis              111.7        (9.57)
Cosmocalanus darwini               27.0        (8.15)
Acrocalanus longicornis           150.1       (33.10)
Scolecithrix bradyi               211.3       (28.51)
Euchaeta rimana                    32.7        (5.68)
Mesocalanus tenuicornis            40.4        (9.57)
Temoridae copepodites             118.4       (58.60)
Scolecitrichopsis ctenopus         127        (32.14)
Scolecithricella marginata         79.2       (35.45)
Calanopia minor                    68.9       (27.79)
Acrocalanus gracilis               45.4       (20.17)
Aetideus bradyi                   131.7       (25.14)
Undinula vulgaris                  59.3       (16.43)
Candacia truncata                  24.4        (6.60)
Subeucalanus crassus              108.6       (18.96)
Pleuromamma piseki                 98.4       (29.13)
Scolecitrichidae copepodites       81.6       (16.01)
Sapphirina metallina               37.1       (14.10)
Scolecithricella abyssalis         73.4       (16.99)
Aetideidae copepodites             73.4       (14.34)
Acrocalanus andersoni              70.1       (14.48)
Pontellidae copepodites            33.2       (10.56)

                                    p            F         post hoc

Total abundance                   0.874         0.13
Calanoida                         0.712        0.342
Cyclopoida                        0.001        9.461       Jan > All
Richness                         < 0.001       10.58        All > M
Diversity                        < 0.001       12.19        All > M
Eucalanidae copepodites           0.001         8.55       All > Jan
Euchaetidae copepodites           0.001         8.14       Jan > All
Subeucalanus subcrassus           0.010         5.11        M > Jan
Euchaeta indica                  < 0.001       17.10       Jan > All
Canthocalanus pauper              0.021         4.18        Jan > M
Nannocalanus minor                0.002         7.32       All > Jan
Subeucalanus subtenuis            0.902
Centropages furcatus             < 0.001       10.92       All > Jan
Metridinidae copepodites          0.530
Temora discaudata                < 0.001       21.90       All > Jan
Calanidae copepodites             0.041         3.41        Jan > M
Oncaea spp                       < 0.001       33.18      Jan > Oct > M
Pleuromamma johnsoni              0.466
Clausocalanus minor              < 0.001       25.82       Jan > All
Euchaeta longicornis              0.947
Cosmocalanus darwini              0.189
Acrocalanus longicornis          < 0.001       11.42        All > M
Scolecithrix bradyi              < 0.001       10.16       All > Jan
Euchaeta rimana                  < 0.001       11.07       Jan > All
Mesocalanus tenuicornis           0.006         5.66        M > Jan
Temoridae copepodites             0.002         7.21       All > Jan
Scolecitrichopsis ctenopus        0.531
Scolecithricella marginata         0.53         0.64
Calanopia minor                  < 0.001       31.02      Jan > Oct > M
Acrocalanus gracilis             < 0.001       14.41        All > M
Aetideus bradyi                   0.011         4.88       Oct > Jan
Undinula vulgaris                 0.002         7.39        All > M
Candacia truncata                 0.001         8.25       Jan > All
Subeucalanus crassus              0.006         5.75       All > Jan
Pleuromamma piseki                0.014         4.66        Oct > M
Scolecitrichidae copepodites      0.198
Sapphirina metallina              0.131
Scolecithricella abyssalis        0.358
Aetideidae copepodites            0.663
Acrocalanus andersoni             0.092
Pontellidae copepodites           0.274

[1.] Abundance values are expressed as ind./1 000 [m.sup.3].

[2.] SE = Standard error.

SIMPER analysis % contribution of species/families contributing to 90
% abundance of the groups defined by the cluster analysis

Group                      Average     # spp

G1: Jan (El Nino)           62.47        32
G2: May Offshore            67.41        25
G3: October (La Nina)       61.53        25
G4: May Eddy and            60.48        18
Coast upwelling

Group                       Species (% Contribution abundance) [1]

G1: Jan (El Nino)        EhJu (12.55), Ein (8.14), EucJu (7.00), Onc
                         (4.09), CalJu (3.98), Cpau (3.73), Ssut
                         (3.44), Er (3.42), Sc (3.36), Clm(3.31), MeJu
                         (3.17), Pj (2.83), Ssu (2.70), El (2.65), Al
                         (2.30), Sma (2.10), Cm (2.07), Ctr (1.84), Sb
                         (1.80), ScJu (1.76), Pso (1.66), ClJ (1.63),
                         Ag (1.49), AtJu (1.48), Cd (1.33), Uv (1.33),
                         CanJu (1.16), Os (1.08), Ab (0.97), Scr
                         (0.82), Ppi (0.80), Td (0.74)

G2: May Offshore         EhJu (11.33), Ssut (7.98), EucJu (6.93), Cd
                         (5.95), Ein (5.81), CalJu (4.67), Mt (4.58),
                         El (4.07), (3.97), Sb(3.63), (3.45), Sd
                         (3.40), Pj (3.00), MeJu (2.81), Ps (2.68),
                         Ctr (2.39), Scu (2.30), Sm (2.03), Ssu
                         (1.64), Nm (1.61), Scr (1.43), Ab (1.37), Td
                         (1.33), AtJu (1.28), Pg (1.16)

G3: October (La Nina)    EucJu (14.12), EhJu (7.22), Ssu (5.32), CalJu
                         (5.19), Sb (4.91), Ssut (4.59), Cpau (4.44),
                         MeJu (4.07), Td (3.92), Nm (3.44), El (3.11),
                         Al (2.90), Ab (2.46), ScJu (2.25), Sa (2.23),
                         Acr (2.19), Cf (2.13), Ein (1.88), AtJu
                         (1.68), Pj (1.52), Sc (1.48), Uv (1.48), Mt
                         (1.39), Coc (1.33), Aa (1.29), TeJu (1.27),
                         PoJu (1.23), Ppi (114)

G4: May Eddy and         EucJu (21.12), SSu (17.57), Cf (7.01), Td
Coast upwelling          (6.97), MeJu (5.27), Sb (5.02), TeJu (4.60),
                         Ssut (3.14), Ab (3.10), Sc (2.89), Scr
                         (2.51), Pj (1.79), Cpau (1.78), Sa (1.78),
                         AtJu (1.75), EhJu (1.68), Nm (1.66), El

[1] Aa = Acrocalanus andersoni, Ag = A. gracilis, Al = A.
longicornis, Ab = Aetideus bradyi, AtJu = Aetideidae copepodites,
CalJu = Calanidae copepodites, CanJu = Candaciidae copepodites, Cd
= Cosmocalanus darwini, Cf = Centropages furcatus, Clj =
Clausocalanus jobei, Clm = C. minor, Cm = Calanopia minor, Coc =
Corycaeus crassiusculus, Cpau = Canthocalanus pauper, Ctr =
Candacia truncata, Ein = Euchaeta indica, El = E. longicornis, Er =
E. rimana, EhJu = Euchaetidae copepodites, EucJu = Eucalanidae
copepodites, MeJu = Metridinidae copepodites, Mt = Mesocalanus
tenuicornis, Nm = Nannocalanus minor, Onc = Oncaeidae spp, Os =
Oithona setigera, Pg = Pleuromamma gracilis, Pj = P. johnsoni, Ppi
= P. piseki, PoJu = Pontellidae copepodites, Ps= Pareucalanus
sewelli, Pso = Pontellina sobrina, Sa = Scolecithricella abyssalis,
Sb = Scolecithrix bradyi, Sc = Scolecitrichopsis ctenopus, ScJu =
Scolecitrichidae copepodites, Sd = Scolecithrix danae, Sm =
Sapphirina metallina, Sma = Scolecithricella marginata, Scr =
Subeucalanus crassus, Ssu = S. subcrassus, Ssut = S. subtenuis, Td
= Temora discaudata, TeJu = Temoridae copepodites, Uv = Undinula

Species and copepodite taxa with significant indicator values for each
cluster group

Cluster / Group      Taxa                          IV [1]      Mean

G1: Jan (El Nino)    Pontellina sobrina             65.3       21.2
                     Clausocalanus jobei            58.2       21.8
                     Clausocalanus furcatus         57.9       16.5
                     Paracalanus aculeatus          49.5       17.5
                     Acrocalanus gracilis           46.8       22.7
                     Candaciidae copepodites        40.1       24.1
                     Lucicutia gaussae              38.8       19.6
                     Paracalanidae copepodites      36.9       17.8
                     Acartia negligens              28.6       11.4
G2: May Offshore     Pareucalanus sewelli           59.3       20.5
                     Pleuromamma gracilis           56.6       14.9
                     Scolecithrix danae             55.9       22.4
                     Centropages gracilis           47.5       16.8
                     Corycaeus robustus             38.7       14.8
                     Corycaeus speciosus            38.4       20.7
                     Aetideus armatus               35.9       19.8
                     Lucicutia flavicornis          29.2       16.8
G3: October          Eucalanus inermis              43.4       14.8
(La Nina)            Corycaeus crassiusculus        38.8       22.2
                     Rhincalanus nasutus            35.7       10.9
G4: May Eddy and     Rhincalanus rostrifrons        21.4       9.4
Coast upwelling

Cluster / Group      Taxa                         (SD [2])      p

G1: Jan (El Nino)    Pontellina sobrina            (6.33)    < 0.001
                     Clausocalanus jobei           (5.96)    < 0.001
                     Clausocalanus furcatus        (6.79)     0.001
                     Paracalanus aculeatus         (6.80)     0.002
                     Acrocalanus gracilis          (5.91)     0.003
                     Candaciidae copepodites       (5.37)     0.012
                     Lucicutia gaussae             (6.24)     0.009
                     Paracalanidae copepodites     (6.66)     0.027
                     Acartia negligens             (6.29)     0.017
G2: May Offshore     Pareucalanus sewelli          (6.62)    < 0.001
                     Pleuromamma gracilis          (6.90)    < 0.001
                     Scolecithrix danae            (5.97)    < 0.001
                     Centropages gracilis          (6.63)     0.002
                     Corycaeus robustus            (6.95)     0.005
                     Corycaeus speciosus           (6.30)     0.015
                     Aetideus armatus              (6.56)     0.030
                     Lucicutia flavicornis         (6.64)     0.047
G3: October          Eucalanus inermis             (6.86)     0.006
(La Nina)            Corycaeus crassiusculus       (5.91)     0.025
                     Rhincalanus nasutus           (6.13)     0.006
G4: May Eddy and     Rhincalanus rostrifrons       (5.49)     0.047
Coast upwelling

[1.] IV = Indicator value.

[2.] SD= Standard deviation.
COPYRIGHT 2018 Universidad de Costa Rica
No portion of this article can be reproduced without the express written permission from the copyright holder.
Copyright 2018 Gale, Cengage Learning. All rights reserved.

Article Details
Printer friendly Cite/link Email Feedback
Author:Kozak, Eva R.; Olivos-Ortiz, Aramis; Franco-Gordo, Carmen; Pelayo-Martinez, Gloria
Publication:Revista de Biologia Tropical
Date:Dec 1, 2018
Previous Article:Humedales del Parque Nacional Chirripo, Costa Rica: caracteristicas, relaciones geomorfologicas y escenarios de cambio climatico.
Next Article:Perfiles estructurales, histoquimicos y fotosinteticos de las agallas inducidas por Eugeniamyia dispar (Diptera: Cecidomyiidae) en las hojas de...

Terms of use | Privacy policy | Copyright © 2020 Farlex, Inc. | Feedback | For webmasters