The marine zooplankton is composed by many organisms that together help to modulate the life of the planet . In this community are found the holoplanktonic polychaetes; organisms that come from their benthic congeners . They have evolved to fully live their entire life cycle in the water column  . The holoplanktonic polychaetes are taxonomically represented by nine families that group about 100 species distributed in all oceans  , particularly in subtropical and tropical regions  .
Although holoplanktonic polychaetes are not abundant or diverse in the different oceanic planktonic assemblages, they are frequent  . Studies have found that the protein content of holoplanktonic polychaetes is of excellent quality  . In addition, it has been shown that they have many feeding strategies, as well as different types of food  . These characteristics allow them to be considered as an ecological group of great importance in the planktonic community because they help in the mineralization processes of organic matter and in the transfer of energy to higher trophic levels such as fish larvae  . They are also considered by some authors as indicators of water masses, being found in greater abundance with waters from upwelling  .
Despite these considerations, studies of holoplanktonic polychaetes are scarce compared to other zooplankton groups such as copepod, chaetognaths, and euphausiids    . These studies have been carried out based on transoceanic campaigns of the last century e.g.  , mainly describing species and their geographical distributions to macroscale . In the Pacific Ocean, the most information of this group has been obtained in the Eastern Tropical Pacific Ocean—ETP from the work by Fernández-Álamo , who identified 44 species from five families, with Alciopidae and Lopadorrhynchidae as the most frequent and diverse.
In the Colombian Pacific basin (CPB), which is part of the ETP, few studies have been carried out on holoplanktonic polychaetes. In addition to the three taxonomic works carried out by Fernández-Álamo   , only one ecological study has been found in the neritic province of the CPB that describes the relationship of the group with the concentration of chlorophyll-α, temperature and salinity during 1996 and 1997 . However, the analyses were not performed at the species level. With these four studies,  species of holoplanktonic polychaetes have been recorded in the CPB for the Colombian Pacific basin, all extracted from stations located in the Neritic province.
It has become imperative to carry out ecological studies on various biological groups to know the different population responses to climate variability and thus complement and improve the understanding of life on the planet, providing new knowledge, as well as data for the predictive models created to understand the effect of climate change    . Studies have shown that organisms can have different responses to similar environmental conditions when contrasted between regions . An example of this can be evidenced by the responses of population densities. When the climate variability increases the temperature, salinity or quantity and quality of the available food from the water, population increases of some species have been found . The same response, in the same species, has also been observed when the values of the environmental variables are inversed or have combinations  .
In the ETP region, one of the main modulators of climate variability is the El Niño Southern Oscillation—ENSO   . This event changes environmental conditions by increasing or decreasing sea temperatures, as well as stratifying or mixing the water column  . Marine communities have been found to adapt to changes caused by ENSO   . In the case of the Colombian Pacific basin, as in the ETP, the studies have focused on the effect and response of warming waters (e.g. El Niño). These studies have found that populations in the Neritic province and coastal region decrease their densities due to the decrease in the quality and quantity of food  .
However, no works have been found in the literature that attempt to study the response at the population level of the different holoplanktonic polychaetes to climate variability. This has caused gaps in knowledge and a lack of necessary supplies to better understand the dynamics of our ecosystems. Therefore, and considering that described above, the aims of this study were 1) described the influence of climate variability on the populations of some species of holoplanktonic polychaetes present in the oceanic province of the Pacific Basin of Colombia; and 2) described the relation between the presence of the species and their population densities with the environmental variables of temperature and salinity, as well as chlorophyll-α concentration.
2.1. Study Area
The Region I of the central ETP is located within the Colombian Pacific Basin—CPB . It has an approximate area of 49,000 km2 (14% of the CPB) (Figure 1) and depths greater than 3000 m. The water column is stable and has no direct influence on coastal environmental factors, making it a sector of purely oceanic water mass  . The annual pattern of water temperature and salinity, as well as rainfalls, are modelled by the Tropical Northeast Pacific warm pool, the Inter Tropical Convergence Zone—ITCZ and the ENSO  .
2.2. Environmental Variables and Climate Variability
Monthly time series of temperature, salinity and chlorophyll-α concentration of
Figure 1. Geographic location of sampling stations in the Colombian Pacific Basin.
the sea surface were used from the works carried out by Cabarcas-Mier  and O’Brien et al. , who obtained the information from satellite images. The Southern Oscillation Index—SOI  was used to identify climate variability caused by El Niño or La Niña events (Appendix—Table A1).
2.3. Field Sampling
The ERFEN program established eight stations in the study area (Figure 1). Samples between 2004 and 2012 were annual between September and October. Additional samples were also obtained between the months of February and March of the 2007-2009 period (Appendix—Table A2). To simplify the analyzes, the samples obtained between February and March are considered as the first half of the year (I); while the samples obtained between September and October are the second half of the year (II). Sampling time depended on the arrival of the vessel at each station. All samples could be used because no significant differences were found between population density values obtained at day and night (Mann-Whitney test p-value > 0.05).
The samples were obtained with oblique trawls in the first 200 m depth using a bongo system with 60 cm diameter rings, 250 µm mesh pore nets and General Oceanic mechanical flowmeter to measure the filtered volume. The samples were fixed and preserved in formalin at 5% final concentration, following the basic suggestions for handling and storage . More detailed procedures are recorded in Uribe-Palomino et al.  and CPPS .
A total of 85 samples were used. The analyzes were made with a fraction that corresponded to 25% of each sample. To do this, each sample had to be divided twice with a Folsom splitter. Polychaetes were removed entirely from each fraction using a stereoscope. Preliminarily, polychaetes were separated at the family level. Only the organisms of the families Alciopidae, Iospilidae and Lopadorrhynchidae were identified and quantified to species level. The Tomoteridae and Typhloscolecidae families, as well as the larvae, were not identified due to deterioration of the specimens or the difficulty of identification. The taxonomic identifications were based on the keys of Dales , Orensanz and Ramírez , Støp-Bowitz , Suárez-Moralez et al.  and de León-González et al. . Identities and the validation of species were verified in the World Registry of Marine Species . The identified material will be deposited in the biological collections of the Museum of Marine Natural History of Colombia (MHNMC)—Makuriwa. The population density was expressed in individuals per 100 cubic meters (ind/100m3).
2.5. Analysis of Data
An environmental description of the study area was made before determining the influence given by both, spatial and temporal environmental conditions on the population densities of holoplanktonic polychaetes, as well as the relation with environmental variables.
To establish if the study area was homogeneous, we carry out a hierarchical agglomerative analysis of weighted average similarity—CLUSTER and non-metric multidimensional scaling—NMDS. These analyses were performed from the databases of temperature, salinity, and chlorophyll-α, using a triangular similarity matrix from the Euclidean distance algorithm. The analyses were realized using the statistical program Primer-e V.7 .
Climate variability was described by plotting monthly Southern Oscillation Index—SOI values between 2004 and 2012, as well as the average SOI values corresponding to the sampling times (i.e., three months before the sampling month, and the sampling month). This sampling time is due to the ecological time of renewal of the zooplankton community  . The neutral period was considered when the SOI values were between −0.25 and 0.25, the presence of El Niño when the values were <−0.25, and the presence of La Niña with values > 0.25. Additionally, anomalies of the environmental variables were calculated for each sampling moment using the analysis of transform Z (E.1) . The general average per month and its standard deviation were calculated from each data set until 2003.
· = Transform Z
· x = Individual value
· = Mean
· S.D. = Standard deviation
To determine the influence by both spatial and temporal environmental conditions on the population densities of holoplanktonic polychaetes, graphic descriptions of the distribution and population density of the species were made considering the environmentally formed sectors and the climate variability of the study period. The population density is expressed in individuals per 100 cubic meters (ind/100m3). To obtain the population density of each species, the following formula was used (E.2):
· Pd = Population density
· Oc = Organism count
· P = Proportion of the revised sample
· F.V. = Filtered volume
Additionally, non-parametric Kruskal-Wallis tests were carried out with the six most frequent species since they met the minimum statistical requirements of the test. The spatial factors were equivalent to the sectors formed and the temporal factors to the ENSO conditions at the time of sampling. The significant tests were complemented with the Dunn test, which establishes the factor that generated the difference.
Finally, to describe the relation between the presence of the species and their population densities with the environmental variables of temperature, salinity, and chlorophyll-α, the environmental ranges of all species were calculated from the presence data. Additionally, environmental variables were correlated with population densities using the Pearson and Spearman tests. Both, the Kruskall-Wallis and Dunn tests and the Spearman and Pearson correlation tests were performed with the statistical package for Excel Real-Statistics v5.4 .
3. Results and Discussion
The Cluster and NMDS justified divide the study area into three sectors (Figure 2). Figure 3 shows the spatial differences. The northern sector was the warmest and less saline ( : 27.5˚C ± 1.0˚C and 32.4˚C ± 0.7˚C), while the southern sector was the less warm and most saline ( : 26.1˚C ± 0.8˚C and 33.3˚C ± 0.7˚C). The central sector presented intermediate values ( : 26.8 ± 0.9˚C and 32.8˚C ± 0.7˚C). The chlorophyll-α concentration, showed similar values among the three sectors ( : north 228.4 ± 84.4 µg/m3; central 270.7 ± 135.3 µg/m3; south 270.1 ± 102.1 µg/m3). However, the central sector presented the most extreme peaks and valleys, while the northern sector had the less time variations.
The condition of La Niña predominated between 2004 and 2012, increasing durations and intensities annually. The period of 2010-2011 was the La Niña strongest. The condition of El Niño appeared at the beginning of the period study and at the end of 2009. The neutral periods presented were short, the most prominent being that of 2012 (Figure 4(a)). Sampling time reflected the general annual pattern of ENSO (Figure 4(b)).
The environmental variables allowed sectorizing the study area latitudinally. The highest values were measured in the northern sector and gradually decreased towards the southern sector, which allowed inferring that it was the dominant variable in this spatial delimitation. Additionally, it coincides with the temperature
Figure 2. Graphical representation of multivariate non-parametric analyses based on data from environmental variables. (a) Cluster; (b) nMDS.
Figure 3. Time series between 2004 and 2012 of environmental values by sector formed in the study area. Vertical lines: sampling month.
Figure 4. Southern Oscillation Index (SOI) during the 2004-2012 period. (a) SOI monthly. (b) SOI during the sampling time (three months before and the sampling month). Vertical lines: sampling month; horizontal lines: limit chosen for the ENSO influence. I: First semester; II: Second semester (see Methods).
patterns modulated by the Tropical Northeast Pacific warm pool , as well as the circulation at the mesoscale level of the region . Previous studies had classified the study area as homogeneous  . However, these previous studies classified the area from a purely physical perspective and not consider the biological communities. It has been shown that the change in a Celsius degree is enough to modify the behaviour of many of the marine species . Therefore, from an ecological perspective, it was necessary to consider the oceanic province of the Colombian Pacific Basin as a heterogeneous area.
The distribution of surface salinity was inverse to temperature, but with its main axis coming from the north coast of Colombia. Coastal rainfall is what modulates the distribution of salinity in the Colombian Pacific Basin—CPB  and the North Pacific coast of Colombia is the place of greatest precipitation in the Tropic . This could cause the E079 station to end up belonging more to the central sector, while the E109 to the north sector.
The distribution of the chlorophyll-α concentration is according to expectations. The Central sector has factors that generate higher productivity and more changes during the year. In this sector is Malpelo Island , which generates local conditions such as upwelling areas , benefiting phytoplankton. Furthermore, due to the movement of nutation of the earth, the ITCZ in the Central sector generates contrasting seasons during the year. In the Northern sector the ITCZ influence is more constant, so rainfall occurs throughout most of the year; different from the Southern sector, where the influence of the ITCZ is less, causing a prolonged drought during the year . Productivity in the North and South sectors reflects this condition, with a lower concentration of chlorophyll-α and fewer changes during the year with respect to the Central sector.
In general, the behaviour of the environmental variables during the sampling time was influenced by climate variability similarly between the sectors and response patterns were evidenced according to ENSO conditions (Figure 5). Temperature, salinity and chlorophyll-α were more influenced in periods of La Niña. For the first variable there was a decrease in the values, while for the other two, the values increased.
The climate variability moderately changed the normal conditions during the study period, increasing or decreasing the values of the environmental variables analysed as described in other studies . It was also determined that for the sampling periods no differentiation was observed in the impact according to the
Figure 5. Anomalies of the environmental variables during the sampling moments (three months before and the sampling month). Vertical lines: northern sector; Dots: central sector; Horizontal lines: southern sector.
sector. ENSO progressively impacts the different coastal regions of the Pacific, with the equatorial sector being the first to receive it. For the study area, the southern sector receives the impact two months after it occurs, while for the northern sector it can take up to four months . The short time of involvement between sectors is reflected in the anomalies since only four months were considered.
The 14 identified species of holoplanktonic polychaetes belonged mostly to the Alciopidae and Lopadorrhynchidae families; only one species belonging to the family Iospilidae appeared. The frequency of occurrence of the species did not exceed 50% of the samples reviewed. Lopadorrhynchus cf. henseni, L. cf. brevis and Pelagobia cf. longicirrata (Lopadorrhynchidae); and Plotohelmis cf. capitata, Rhynchonereella cf. petersii and R. cf. gracilis (Alciopidae) were the species that had occurrence frequencies greater than 10% and represented 85% of the total population density (Appendix—Figure A1).
The low frequency of occurrence and population density of holoplanktonic polychaetes was consistent with that registered in various studies conducted in open water of the Pacific Ocean    , contrary to what some authors have mentioned for studies done in coastal areas in other oceans  . This difference in the ecological system can cause patterns to change : low productivity areas similar to ocean province, have greater competition for the resource, while areas with external contributions (i.e. river plume or upwelling) as coastal sectors, allow all organisms to increase their densities more easily .
Regardless of the technical reasons related to capture and analysis methods , it was hypothesized that low frequencies of occurrence and population densities due to aspects such as the great mobility of holoplanktonic polychaetes , the rapid reaction of movement to water disturbances , the capacity that allows them to fractionate and regenerate again , and natural migrations to mesopelagic depths . These factors would facilitate evasion to be captured. It could also be considered that, compared to copepods, holoplanktonic polychaetes are less adapted to the environment  and their presence in the zooplanktonic community is more recent than most other groups , which would give them disadvantages of colonization of niches.
The distribution and population density by sector and sampled period were different for all species (Figure 6) (Appendix—Table A3). The species of Alciopidae appeared in more sampling periods and were better distributed in the three sectors contrasting for to Lopadorrhynchidae species. In general, the species increased population densities at the end of the period studied, with a decrease in densities during 2009, as observed more clearly with Plotohelmis cf. capitata and Pelagobia cf. longicirrata (Figure 6(g), Figure 6(m)).
Of the Alciopidae, Rhynchonereella cf. moebii and Krohnia cf. lepidota appeared in two periods with averages of population density below 6 ind/100m3. The first species appeared during 2007 in the northern and central sectors, while the second was between the second semester of 2007 and the first of 2008 in the
Figure 6. Mean population density and its standard deviations by sector of quantified holoplankton polychaete species during the 2004-2012 period. Square: northern sector; Circle: central sector; Triangle: southern sector. Red: El Niño condition; Blue: La Niña condition; No color: Neutral condition.
central and southern sectors (Figure 6(a), Figure 6(b)). Plotohelmis cf. alata and Alciopina cf. parasitica were collected in four and five sampling periods respectively, with 2008 being the best spatially represented for the second species. The average population densities were less than 5 ind/100m3. Both species were collected in the three sectors, although for A.cf. parasitica it was more predominant in the northern sector, representing more than 50% of appearance and density population (Figure 6(c), Figure 6(d)). Rhynchonereella cf. gracilis was collected in five sampled periods and, like A. cf. parasitica, 2008 was the year in which it was best represented spatially, but with average population densities less than 2 ind/100m3. This species appeared in the three sectors, with the south being the least present (Figure 6(e)).
Rhynconereella cf. petersii and Plotohelmis cf. capitata, were those that occurred in the greatest number of sampled periods, with eight and ten times, respectively. Both species were not collected in 2006, and the case of R. cf. petersii, they did not appear in 2010 and 2012. The average population densities of R. cf. petersii were mostly below 5 ind/100m3; while for P. cf. capitata the average population densities were divided equally below and above the 5 ind/100m3. Like Alciopina cf. parasitica, for R. cf. petersii and P. cf. capitata there was a decrease in population density during 2008. Both R. cf. petersii and P. cf. capitata were presented in the three sectors, the central one being the one with the highest frequency of appearance and the south with the highest population densities (Figure 6(f), Figure 6(g)).
As for the species of the Lopadorrhynchidae, specimens of Maupasia cf. gracilis were collected only during 2004 with an average density above 5 ind/100m3 (Figure 6(h)). Pedinosoma cf. curtum and Lopadorrhynchus cf. krohnii were presented during four sampled periods, the first species in the middle and end of the sampling period, while the second was distributed throughout the study period. The average population density of P. cf. curtum was close to 4 ind/100m3, while for L. cf. krohnii it started above 10 ind/100m3 and was decreasing over time to a value less than 5 ind/100m3. Pedisoma cf. curtum was collected in all three sectors, while L. cf. krohnii only in the north (Figure 6(i), Figure 6(j)).
Lopadorrhynchus cf. brevis was collected at the beginning and end of the study period for a total of six of the sampled periods. It presented average population densities variable in time, being the highest at the beginning and end of the study (2004 and 2012), showing a similar pattern to some of the species of the Alciopidae with a decrease in densities towards the middle of the period. This species was present in all three sectors, being predominant in the south in the first sampling periods, and in the north and central the last sampled periods (Figure 6(k)).
The last two species of the Lopadorrhynchidae, Lopadorrhynchus cf. henseni and Pelagobia cf. longicirrata were collected in ten of the eleven sampling periods. They differed in that L. cf. henseni specimens were not collected in 2006 and P. cf. longicirrata in 2004. The population densities of both species varied over time, showing a pattern of peaks and valleys. Both species were found in all three sectors, however, the highest densities for L. cf. henseni were in the northern sector, while for P. cf. longicirrata in the central and southern areas (Figure 6(l), Figure 6(m)). Finally, Phalacrophorus cf. pictus, of the Iospilidae family, appeared once during the second semester of 2007 in the central sector (Figure 6(n)).
Regarding the spatio-temporal distribution, Phalacrophorus cf. pictus, Alciopina cf. parasitica, Rhynchonereella cf. moebii, Kronhia cf. lepidota, Plotohelmis cf. alata, Maupasia cf. gracilis, Pedisoma cf. curtum, and Lopadorrhynchus cf. krohnii were considered rare species for CPB, as has been recorded in other regions such a North Pacific, ETP or Caribbean Sea   . Not knowing the environmental requirements of these species, it is difficult to find an answer to this result. However, it seems that in some species its appearance coincided with some environmental characteristics. For example, L. cf. krohnii only appeared in the northern sector, when climate variability allowed increases in water temperature; P. cf curtum appeared when the intensity of La Niña began to decrease, and A. cf. parasitica had its greatest presence during La Niña from 2008-2009.
The other six species identified in the present study, Rhynchonereella cf. gracilis, R. cf. petersii, Plotohelmis cf. capitata, Lopadorrhynchus cf. brevis, L. cf. henseni and Pelagobia cf. longicirrata presented a wider geographical and temporal distribution. This same pattern has been described in other studies from Pacific Ocean   , Atlantic Ocean and Caribbean Sea   . In all of them, P. cf. longicirrata was the most common species.
No statistical differences in population density between sectors were determined from the three most frequent species of both families, however, climate variability did influence the Lopadorrhynchus cf. henseni and L. cf. brevis (Kruskall-Wallis p-value < 0.05). In both cases, the population densities during the El Niño periods were significantly higher compared to the periods with La Niña and neutral for the first of the species, and with La Niña for the second (Dunn p-value < 0.02) (Appendix—Table A4, Figure A2).
The different sectors found from the environmental variables did not influence population densities in the six most common species. This could be due to the fact that the environmental differences were not marked enough to generate a natural barrier , and that the species have sufficient tolerance ranges to be able to live in other regions   . The latter fact could be verified with the environmental ranges of the species in the present study.
The holoplanktonic polychaetes presented different ranges of the environmental variables conditioned with the frequency of appearance (Figure 7). In the case of the species of the Alciopidae they tended to warmer temperatures, and higher salinities and chlorophyll-α concentrations, while the species of the Lopadorrhynchidae tended to colder temperatures and less salinity waters but with high chlorophyll-α concentrations. The most frequent species with the highest population densities (Plotohelmis cf. capitata and Pelagobia cf. longicirrata) were present in the entire range of temperature, salinity and chlorophyll-α concentration. The species with intermediate frequency and population density (Rhynchonereella cf. petersii, R. cf. gracilis, Lopadorrhynchus cf. henseni, and L. cf. brevis)
Figure 7. Ranges of environmental variables where holoplankton polychaete species appeared in the oceanic province of Colombian Pacific Basin during 2004-2012. (n): number of samples; (A): Alciopidae; (L): Lopadorrhynchidae; (I): Iospilidae.
coincided with the decrease in the range of at least one of the variables.
Less frequent species with low population densities (Alciopina cf. parasitica, Plotohelmis cf. alata, Lopadorrhynchus cf. krohnii and Pedisoma cf. curtum) coincided with a greater decrease in environmental ranges in general, with tendencies to warmer waters and medium salinities. They differ in chlorophyll-α concentrations, in which species of the Lopadorrhynchidae tended to appear with low concentrations. Finally, the species that were collected at two or fewer sampling times (Rhynchonereella cf. moebii, Krohnia cf. lepidota, Maupasia cf. gracilis, and Phalacrophorus cf. pictus) tended to appear with average values of each of the variables (Figure 7).
Regarding the correlations between population density of species and environmental variables, significant and positive statistical differences with temperature were found in Lopadorrhynchus cf. henseni (Pearson 0.32; p-value < 0.1) and L. brevis (Pearson 0.60; p-value < 0.05), with salinity in Pelagobia cf. longicirrata (Pearson 0.30; p-value < 0.1) and with chlorophyll-α concentration in Rhynchonereella cf. gracilis (Pearson 0.80; p-value < 0.05). Rhynchonereella cf. petersii was the only species to having a negative statistical significance with temperature (Pearson −0.58; p-value < 0.1) (Appendix—Table A5, Figure A3).
The population density of Lopadorrhynchus cf. henseni and L. cf. brevis showed a positive correlation with temperature and increased the population densities with El Niño compared to La Niña and Neutral condition. There are no data that support an explanation, but it was hypothesized that their densities are related to the fact that these species are considered omnivorous , added to the evidence that during El Niño the presence and abundance of omnivorous species increases  .
No evidence why Rhynchonereella cf. petersii showed negative correlation between population density and temperature. According to Jumars et al. , this species is omnivorous with a tendency to herbivore. Continuing with the approach to the increase in omnivores due to the influence of El Niño, it can be hypothesized that this species is herbivorous in the study area; therefore it is affected with food decrease.
The decrease of population density by El Niño has also been described for the entire assembly of holoplanktonic polychaetes during the period with the El Niño of 1997, one of the strongest in history . El Niño condition in the present study it can be considered weak due to its low intensity and duration. This reason allows us to infer that R. cf. petersii may be the only one that is at the limit of environmental tolerance with respect to the other holoplanktonic polychaetes.
Rhynchonereella cf. gracilis showed a contrary pattern to that found in coastal regions with respect to correlation between population density and chlorophyll-α concentration  . It is a species considered carnivorous , so possibly its positive correlation is due to the presence of its prey. However, there is no evidence to support it.
Finally, the population density of Pelagobia cf. longicirrata showed a positive correlation with salinity. As has happened with the species analyzed above, there are no previous studies that have related these variables and attributed any reason. It is hypothesized for this case that its morphology, a depressed body, different from the other identified species  may allow it an advantage to move when water density increases due to salinity. This morphological advantage, in turn, would explain why it is the most frequent species of holoplankton polychaetes. It would have an advantage in capturing prey from the pico- and nanoplankton since these organisms, being smaller, would find it difficult to move to avoid being predated .
This study provides new ecological information on 14 species of holoplanktonic polychaetes with respect to environmental variables of temperature, salinity, and chlorophyll-α concentration, as well as climate variability in a tropical region.
In ecological terms, the Colombian Pacific basin during the 2004-2012 climate variability presented a spatial differentiation from north to south, with the greatest differences being found between the North and South sectors, while the central sector was the most fluctuating annually. Climate variability during the study period was dominated by La Niña and did not differentially influence the sectors.
The present study confirms that the planktonic assemblages of the oceanic province present a high frequency and low population densities of holoplanktonic polychaetes. The 14 species identified did not show a spatially or temporally differentiated distribution. Only the species Lopadorrhynchus cf. henseni and L. cf. brevis were influenced by climate variability, increasing their populations with the influence of El Niño.
Finally, five species showed correlations with the environmental variables. L. cf. henseni and L. cf. brevis had a positive correlation with temperature; Rhynchonereella cf. petersii had a negative correlation with temperature, and R. cf. gracilis and Pelagobia. cf. longicirrata had a positive correlation with chlorophyll-α concentration and salinity, respectively.
J.M.G.-S. carried out the separation, quantification and taxonomic identification of organisms, data analysis and wrote the manuscript as part of his Ph.D. thesis. A.C.-O. carried out the taxonomic identification and validation of the organisms. N.H.C. guided the structuring, edited, and made the final revision of the manuscript. The authors declare that they have no conflict of interest
We express our gratitude to 1) CECIMAR of the Universidad Nacional de Colombia (UNAL) Caribbean headquarters for allowing us to review the samples and have all the logistic support of the laboratory, 2) the Universidad del Valle, administrators of the analysed samples and allowed their unrestricted use, 3) to the MAKURIWA museum of INVEMAR for allowing work in their laboratories and guarding the biological samples throughout the process and 4) COLCIENCIAS for financial support to JMG-S. to devote full time to doctoral work and completion of the internship. Personally, J.M.G.-S. Thanks to his parents Alicia and José Ángel for the support always and in every way to complete the doctoral stage. This manuscript is equivalent to contribution No. 513 of CECIMAR and No. CTRB-1295 of INVEMAR.
This work was supported by the Departamento Administrativo de Ciencia, Tecnología e Innovación—COLCIENCIAS with the grant of national doctorates [Grant number 727 to J.M.G.-S.]; the Universidad Nacional de Colombia—UNAL, Caribbean headquarters [39195 to J.M.G.-S.]; and the Museo de Historia Natural Marina de Colombia—MAKURIWA of the Instituto de Investigaciones Marinas y Costeras—INVEMAR contributes with the loan of the equipment, and the laboratory for the research in question.
Supplementary data from sampling stations and environmental data; results of population analyses; and data on population density by species.
Table A1. Oceanographic and climatic data sources used for the present study. SST: sea surface temperature; SSS: sea surface salinity; SSC: sea surface chlorophyll-α concentration; SOI: Southern Oscillation Index. (): no units.
Table A2. Information stations sampled. The sample identifier (ID.) is the combination of the station number, the semester, and the sampled year.
Table A3. Population density (ind/100m3) of holoplankton polychaetes found in the Pacific Ocean Basin of Colombia during the climate variability from 2004 to 2012.
Figure A1. Proportion (%) of the biological attributes of quantified holoplanktonic polychaetes during the 2004-2012 period. A. Frequency of occurrence; B. Population density. A.p.: Alciopina cf. parasitica; K.l.: Krohnia cf. lepidota; P.a.: Plotohelmis cf. alata; P.ca.: Plotohelmis cf. capitata; R.g.: Rhynchonereella cf. gracilis; R.m.: Rhynchonereella cf. moebii; R.p.: Rhynchonereella cf. petersii; L.b.: Lopadorrhynchus cf. brevis; L.h.: Lopadorrhynchus cf. henseni; L.k.: Lopadorrhynchus cf. krohnii; M.g.: Maupasia cf. gracilis; P.cu.: Pedinosoma cf. curtum; P.l.: Pelagobia cf. longicirrata; P.p.: Phalacrophorus cf. pictus.
Table A4. P-values of the Kruskall-Wallis test to determine if there are differences in the population densities of holoplankton polychaete species between sectors (North, Central and South) and by climate variability (El Niño—No, La Niña—Na and Neutral—N). *p-value < 0.05. **p-value < 0.017.
Figure A2. Graphical representation of the population density of two species of holoplanktonic polychaetes with statistical significance from the Kruskal-Wallis tests.
Table A5. Correlations between population density and environmental variables. n: number of data used; L: data transformed with logarithm base 10 to obtain normality; D: population density; T: temperature; S: salinity; C: chlorophyll-α concentration; P: Pearson test; S: Spearman test; *p < 0.1; **p < 0.05.
Figure A3. Relations between environmental variables and population density of holoplanktonic polychaete species with statistical significance based on correlation tests.
 Jumars, P.A., Dorgan, K.M. and Lindsay, S.M. (2015) Diet of Worms Emended: An Update of Polychaete Feeding Guilds. Annual Review of Marine Science, 7, 497-520.
 Guglielmo, R., Gambi, M.C., Granata, A., Guglielmo, L. and Minutoli, R. (2014) Composition, Abundance, and Distribution of Holoplanktonic Polychaetes within the Strait of Magellan (Southern America) in Austral Summer. Polar Biology, 37, 999-1015.
 Jiménez-Cueto, S., Suárez-Morales, E. and Morales-Ramírez, A. (2012) Algunos poliquetos holoplanctónicos (Annelida: Polychaeta) del Parque Nacional Isla del Coco, Costa Rica. Revista de Biología Tropical, 60, 207-222.
 Guglielmo, L., Minutoli, R., Bergamasco, A., Granata, A., Zagami, G. and Antezana, T. (2011) Short-Term Changes in Zooplankton Community in Paso Ancho Basin (Strait of Magellan): Functional Trophic Structure and Diel Vertical Migration. Polar Biology, 34, 1301-1317.
 Suárez-Morales, E., Jiménez-Cueto, S. and Salazar-Vallejo, S.I. (2005) Catálogo de los poliquetos pelágicos (Polychaeta) del Golfo de México y mar Caribe mexicano. El Colegio de la Frontera Sur (ECOSUR), Campeche, México.
 Salazar-Vallejo, S.I. and Londoño-Mesa, M.H. (2004) Lista de especies y bibliografía de poliquetos (Polychaeta) del Pacífico Oriental Tropical. Anales del Instituto de Biología, Universidad Nacional Autónoma de México, Serie Zoología, 75, 9-97.
 Fernández-álamo, M.A. (2006) Composition, Abundance, and Distribution of Holoplanktonic polychaetes from the Expedition “El Golfo 6311-12” of Scripps Institution of Oceanography. Scientia Marina, 70, 209-215.
 Støp-Bowitz, C. (1992) Polychêtes pélagiques des campagnes de ‘L’Ombago’ dans les eaux équatoriales et tropicales ouest-africaines. Institut Francais de Recherche Scientifique pour le Développement en Coopération, France.
 Buzhinskaja, G.N. (2017) Descriptions of Pedinosoma curtum Reibisch from North-West Pacific and a Late Larva of Pedinosoma polaris sp. nov. from the Arctic Basin (Polychaeta: Phyllodocida: Lopadorhynchidae). Zoosystematica Rossica, 26, 3-10.
 Fernández-álamo, M.A. (1983) Los poliquetos pelágicos (Annelida-Polychaeta) del Pacífico Tropical Oriental: sistemática y zoogeografía. PhD Thesis, Universidad Nacional Autónoma de México, Mexico City.
 Fernández-álamo, M.A. (2004) Distribution of Holoplanktonic Typhloscolecids (Annelida-Polychaeta) in the Eastern Tropical Pacific Ocean. Journal of Plankton Research, 26, 647-657.
 Carvajal-Pinilla, L.A., Vergara-Castillo, C.F. and López-Peralta, R.H. (2009) Chaetognatha, Thaliacea, Euphausiacea and Pelagic Polychaeta in the Colombian Pacific Ocean during Two Periods in 1996 (La Niña) and Two Periods in 1997 (El Niño). Revista Facultad de Ciencias Básicas, 5, 172-185.
 Gajbhiye, S.N. (2002) Zooplankton-Study Methods, Importance and Significant Observations. Proceedings of the National Seminar on Creeks, Estuaries and Mangroves Pollution and Conservation, Mumbai, India, 28-30 November 2002, 21-27.
 Mair, L., Hill, J.K., Fox, R., Botham, M., Brereton, T. and Thomas, C.D. (2014) Abundance Changes and Habitat Availability Drive Species’ Responses to Climate Change. Nature Climate Change, 4, 127-131.
 McQuatters-Gollop, A., Johns, D.G., Bresnan, E., Skinner, J., Rombouts, I., Stern, R., Aubert, A., Johansen, M., Bedford, J. and Knights, A. (2017) From Microscope to Management: The Critical Value of Plankton Taxonomy to Marine Policy and Biodiversity Conservation. Marine Policy, 83, 1-10.
 Uribe-Palomino, J., López, R., Gibbons, M.J., Gusmão, F. and Richardson, A.J. (2018) Siphonophores from Surface Waters of the Colombian Pacific Ocean. Journal of the Marine Biological Association of the United Kingdom, 99, 67-80.
 Molina-González, O., Lavaniegos, B.E., Gómez-Valdés, J. and de la Cruz-Orozco, M. (2018) Holoplanktonic mollusks off Western Baja California during the weak El Niño 2006-07 and further transition to La Niña. American Malacological Bulletin, 36, 79-95.
 Velasco, J., Gutiérrez-Cánovas, C., Botella-Cruz, M., Sánchez-Fernández, D., Arribas, P., Carbonell, J.A., Millán, A. and Pallarés, S. (2019) Effects of Salinity Changes on Aquatic Organisms in a Multiple Stressor Context. Philosophical Transactions of the Royal Society B: Biological Sciences, 374, Article ID: 20180011.
 Amador, J.A., Rivera, E.R., Durán-Quesada, A.M., Mora, G., Sáenz, F., Calderón, B. and Mora, N. (2016) The Easternmost Tropical Pacific. Part I: A Climate Review. Revista de Biología Tropical, 64, S1-S22.
 Amador, J.A., Durán-Quesada, A.M., Rivera, E.R., Mora, G., Sáenz, F., Calderón, B. and Mora, N. (2016) The Easternmost Tropical Pacific. Part II: Seasonal and Intraseasonal Modes of Atmospheric Variability. Revista de Biología Tropical, 64, S23-S57.
 Akbari, E., Alavipanah, S., Jeihouni, M., Hajeb, M., Haase, D. and Alavipanah, S. (2017) A Review of Ocean/Sea Subsurface Water Temperature Studies from Remote Sensing and Non-Remote Sensing Methods. Water, 9, Article No. 936.
 Rebstock, G.A. (2002) Climatic Regime Shifts and Decadal-Scale Variability in Calanoid Copepod Populations off Southern California. Global Change Biology, 8, 71-89.
 Eiane, K., Espinasse, M. and Espinasse, B. (2018) Environmental Effects on Zooplankton Abundance on a Sub-Arctic Shelf off Northern Norway. Aquatic Biology, 27, 75-86.
 Gómez, F.A., Lee, S.-K., Hernández, F.J., Chiaverano, L.M., Muller-Karger, F.E., Liu, Y. and Lamkin, J.T. (2019) ENSO-Induced Co-Variability of Salinity, Plankton Biomass and Coastal Currents in the Northern Gulf of Mexico. Scientific Reports, 9, Article No. 178.
 Cabarcas-Mier, A.M. (2017) Anomalías de temperatura superficial del mar y su relación con eventos atmosféricos en la Cuenca Pacífica colombiana. MSc Thesis, Escuela Naval de Cadetes Almirante Padilla, Cartagena, Colombia.
 O’Brien, T.D., Lorenzoni, L., Isensee, K. and Valdés, L. (2017) What Are Marine Ecological Time Series Telling Us about the Ocean? A Status Report. IOC Technical Series, International Olympic Committee, Paris, France, 297 p.
 de León-González, J.A., García-Garza, M.E., Peña-Rivera, A., Carreña-Parra, L.F., Salazar-Vallejo, S.I., Solís-Weiss, V. and Bastida-Zavala, J.R. (2009) Poliquetos (Annelida: Polychaeta) de México y América Tropical. Universidad Autónoma de Nuevo León, Monterrey.
 Hidalgo, P., Escribano, R. and Morales, C.E. (2005) Annual Life Cycle of the Copepod Eucalanus inermis at a Coastal Upwelling Site off Mejillones (23˚S), Northern Chile. Marine Biology, 146, 995-1003.
 Shimode, S., Takahashi, K., Shimizu, Y., Nonomura, T. and Tsuda, A. (2012) Distribution and Life History of the Planktonic Copepod, Eucalanus californicus, in the Northwestern Pacific: Mechanisms for Population Maintenance within a High Primary Production Area. Progress in Oceanography, 96, 1-13.
 Marcaillou, B., Charvis, P. and Collot, J.-Y. (2006) Structure of the Malpelo Ridge (Colombia) from Seismic and Gravity Modelling. Marine Geophysical Researches, 27, 289-300.
 Quintanilla, E., Madurell, T., Wilke, T. and Sánchez, J.A. (2019) Dynamic Interplay of ENSO Events and Local Hydrodynamic Parameters Drives Demography and Health Status of Gorgonian Sea Fan Populations on a Remote Tropical Eastern Pacific Island. Frontiers in Marine Science, 6, Article No. 694.
 Christensen, J.H., Krishna-Kumar, K., Aldrian, E., An, S.-I., Cavalcanti, I.F.A., de Castro, M., Dong, W., Goswami, P., Hall, A., Kanyanga, J.K., Kitoh, A., Kossin, J., Lau, N.-C., Renwick, J., Stephenson, D.B., Xie, S.-P. and Zhou, T. (2013) Climate Phenomena and Their Relevance for Future Regional Climate Change. In: Stocker, T.F., Qin, D., Plattner, G.-K., Tignor, M., Allen, S.K., Boschung, J., Nauels, A., Xia, Y., Bex, V. and Midgley, P.M., Eds., Climate Change 2013: The Physical Science Basis. Contribution of Working Group I to the 5th Assessment Report of the Intergovernmental Panel on Climate Change, Cambridge University Press, Cambridge and New York, 1217-1308.
 Fernández-álamo, M.A. (1996) Holoplanktonic polychaetes off the Southwestern Coast of Baja California, México, in March, 1997. Anales del Instituto de Biología, Universidad Nacional Autónoma de México, Serie Zoología, 67, 51-66.
 Griffiths, J.R., Hajdu, S., Downing, A.S., Hjerne, O., Larsson, U. and Winder, M. (2016) Phytoplankton Community Interactions and Environmental Sensitivity in Coastal and Offshore Habitats. Oikos, 125, 1134-1143.
 Pelayo-Martínez, G., Olivos-Ortiz, A., Franco-Gordo, C., Quijano-Scheggia, S., Gavino-Rodríguez, J., Kono-Martínez, T. and Castro-Ochoa, F. (2017) Physical, Chemical and Zooplankton Biomass Variability (Inshore Offshore) of Mexican Central Pacific during El Niño-La Niña 2010. Latin American Journal of Aquatic Research, 45, 67-78.
 Jiménez-Cueto, S. and Suárez-Morales, E. (2008) An Account of Alciopina, Torrea, and Rhynconereella (Polychaeta: Alciopidae) of the Western Caribbean Sea. The Belgian Journal of Zoology, 138, 70-80.
 Cárdenas-Oliva, A.V., Díaz, O.F. and Márquez, B. (2010) Caracterización taxonómica de los poliquetos holoplanctónicos (Annelida: Polychaeta) de la plataforma norte de la Península de Paria y Golfo de Paria, Venezuela. Boletín del Instituto Oceanográfico de Venezuela, 49, 53-63.
 Collazo, N., Hernández, F., Soldevilla, F., Vera, A., Núñez, J. and Fraile-Nuez, E. (2017) Poliquetos planctónicos relacionados con enclaves de vulcanismo reciente en Canarias. Vieraea, 45, 89-118.
 Castillo, F. and Vizcaino, Z. (1992) Los indicadores biológicos del fitoplancton y su relación con el fenómeno de El Niño 1991-92 en el Pacífico colombiano. Boletín Científico CIOH, No. 12, 13-22.