Variabilidad estacional de la estructura comunitaria y abundancia de copepodos, modificada por la transicion El Nino-La Nina (2010) en el Pacifico, 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.
MATERIALS AND METHODS
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 (http://www.aviso.ahimetry.fr/en/data.html) were processed to identify sea level anomalies and geostrophic currents (with 1/3[degrees] resolution); both parameters were plotted using [R]MATLAB 184.108.40.2064 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. ncep.noaa.gov/products/analysis_monitoring/ 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; https://www.pfeg.noaa.gov/ 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 http://copepodes.obs-banyuls.fr/en
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; firstname.lastname@example.org, email@example.com, firstname.lastname@example.org
(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; email@example.com
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).
TABLE 1 Environmental variable mean for three periods in 2010 Variable  Jan (SE ) May (SE ) 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  Oct (SE ) 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 applicable 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 MEI [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. TABLE 2 Mean densities by sampling period of copepod groups and adult and copepodite taxa with abundance > 0.5 ind./1 000 [m.sup.3] Jan  SE  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  SE  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  SE  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. TABLE 3 SIMPER analysis % contribution of species/families contributing to 90 % abundance of the groups defined by the cluster analysis Group Average # spp Similarity 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)  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.57)  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 vulgaris. TABLE 4 Species and copepodite taxa with significant indicator values for each cluster group Cluster / Group Taxa IV  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 ) 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.