Introduction
Copepods are a group of microcrustaceans considered the most abundant metazoans on the planet, particularly the small-sized forms (Boxshall & Halsey, 2004; Huys & Boxshall, 1991; Turner, 2004). They are among the best studied organisms of the marine zooplankton (Campos-Hernández & Suárez-Morales, 1994). This interest is due to their importance within the dynamics of the marine ecosystem where they represent between 60 and 97% of the zooplankton biomass in neritic-coastal and oceanic areas (Björnberg, 1981; Bradford-Grieve et al., 1999). They are an essential part of the biological pump, which represents a key process to understand the carbon flux in the oceans, where the microbial decomposition of their remains contributes to the biogeochemical transformation of organic matter flowing into the deep sea (Glud et al., 2015).
Their trophic role is highly relevant and broad in the pelagic community, as the group includes omnivorous, herbivorous and carnivorous forms. Most are phytoplankton-consuming herbivores, thus playing a pivotal role in transferring this energy to higher trophic levels and showing a variety of adaptations that make them a highly successful group in the pelagic environments (Chen et al., 2018; Kiørboe, 2011). Therefore, fisheries worldwide depend, at least partly, on the biomass, distribution and trophic dynamics of the zooplankton community (Saiz et al., 2007), as copepods are the main food item for larval and juvenile fish (Tilley et al., 2016). Consequently, the diversity and distribution patterns of the planktonic copepods have been extensively studied (Piontkovski & Landry, 2003). These works are a current priority, as the plankton community is sensitive to processes related to global warming (Comeau et al., 2009; Rombouts et al., 2009). Despite its ecological and economic importance, the information on the biology and diversity of pelagic copepods carried in the Southern Caribbean ecoregion is scarce.
The Southern Caribbean (SCA) is a Marine Ecoregion of the World (MEOW) that belongs to the Tropical Northwestern Atlantic province (TNWA) or “Greater Caribbean” (Spalding et al., 2007, 2012). The SCA ecoregion is characterized by its high productivity, a condition resulting from prevalent coastal upwelling processes and nutrient inputs from the Orinoco and Amazon rivers (Castellanos et al., 2002; Muller-Karger & Varela, 1990; Okolodkov, 2003; Rueda-Roa et al., 2018). It includes the Exclusive Economic Zones (EEZs) of several Caribbean countries including: Aruba, Curaçao, Bonaire, Trinidad and Tobago, Venezuela and part of Colombia (La Guajira Peninsula). Venezuela’s EEZ is about 58% of the total area of the ecoregion, which extends to 16º44’49” N, because of Isla de Aves position (Flanders Marine Institute, 2021). The main physical and bioecological characteristics of the ecoregion were described by Villamizar and Cervigón (2017) and Correa-Ramírez et al. (2020).
Since its beginnings, several decades ago, research on the SCA area has focused on the neritic zone off eastern Venezuela (Cervigón, 1963; Legaré, 1961; Zoppi, 1961), and some other studies involved oceanic and mesopelagic environments of this area (i.e., MHI Archive, 1970; Owre & Foyo, 1964; Wilson, 1942), known to host an increased zooplankton biodiversity (Angel, 1993). The first diversity evaluations of the Venezuelan pelagic copepods yielded a total of 115 to 210 species (Miloslavich et al., 2005; Rodríguez & Suárez, 2003), but these data only incorporate general references. In a recent review for eastern Venezuela (Gulf of Cariaco), Márquez-Rojas et al. (2020) reported 136 copepod species, a figure that includes some parasitic forms found in the plankton. Contrastingly, the species richness of other, less abundant groups of planktonic crustaceans like amphipods and euphausiids is better documented in Venezuela (González-Cebrero et al., 2017; Miloslavich et al., 2005). This work aims to offer basic information on the diversity of the pelagic copepods from the SCA in order to expand the knowledge of the zooplankton community in this ecoregion and motivate future studies on this group of pelagic crustaceans.
Materials and methods
In this work we compiled unpublished data obtained in 6 oceanographic cruises carried out on board the M/S “La Salle” between September 1967 and July 1968. Three stations were visited (Fig. 1), station 11 near the Venezuelan coastline (10º28’60” N, 64º26’00” W), at the western part of the Cariaco Basin. Station 63 (11º13’60” N, 65º00’00” W), was located between La Tortuga and La Blanquilla islands, and station 76 (11º45’00” N, 66º56’00” W), south of Los Roques Archipelago.
During the oceanographic cruises, zooplankton samples were taken with a Clark-Bumpus collector (net No. 20, mesh size = 76 µm), performing horizontal hauls at 0, 25, 50, 100, 200, 500, and 800 m depths. In addition, vertical trawls were made from depths of 800, 1,000, or 1,200 m to the surface with a standard plankton net No. 2 (mesh size = 333 µm). All samples were obtained during the day between 09:00 and 18:00 h, preserved with 5% formaldehyde buffered with sodium borate and later processed for their taxonomic identification. The taxonomic keys, descriptions and illustrations by Rose (1933), Farran (1936), Mori (1937), Sewell (1947), Tanaka (1960), Vervoort (1963, 1965) and Owre & Foyo (1967) were used for initial taxonomic identification of the material. Additional sampling information, species composition and vertical distribution is available at https://doi.org/10.15468/jzyssw
In order to compile this checklist of the pelagic copepod species from the SCA (Appendix), the information was supplemented by a review of the available regional literature, from which we selected 28 references including articles, technical reports and dissertations. References containing data already reported in the previous works selected for this review were excluded. The references are arranged chronologically, these studies were conducted under varying conditions and sampling gears (Table 1).
# | Reference | Studies conditions and locations |
---|---|---|
1 | Wilson (1942) | Cruise VII (1928-1929), R/V Carnegie. Venezuela Basin, Sta. Nº 31, 14º46’00” N, 63º26’00” W. Sta. Nº 32, 15º18’00” N, 68º11’00” W. |
2 | Bowman (1957) | Wilson (1942) emendation, only new species description. Venezuela Basin, Sta. Nº 32, 15º18’00” N, 68º11’00” W. |
3 | Legaré (1961) | Cariaco project, Cruise 1960, R/V Guaiquerí. Gulf and Cariaco Basin, 10º27’45” to10º33’34” N, 63º44’54” to 64º23’30” W. |
4 | Zoppi (1961) | Cariaco project, Cruise 1960, R/V Guaiquerí. Gulf and Cariaco Basin, 10º30’00” to 10º30’52” N, 63º58’49” to 64º21’06” W. |
5 | Cervigón (1963) | Cruise 1960, M/S La Salle. Margarita Island (Nueva Esparta State) and Cariaco Basin, 10º15’00” to 10º55’00” N, 63º45’00” to 65º10’00” W. |
6 | Cervigón (1964) | Cruises (1963-1964), M/S La Salle. Only Corycaeidae. Cariaco Basin, north of the Sucre State, Gulf of Paria (Venezuela) and north of Trinidad, 10º20’38” to 12º02’22” N, 63º13’03” to 65º26’02” W. |
7 | Legaré (1964) | Cruises (1960-1961), R/V Guaiquerí. Single Station. Cariaco Basin, 10º30’15” N, 64º20’22” W. |
8 | Owre and Foyo (1964) | Caribbean Sea Expedition “Carib”, R/V Spencer F. Baird. Venezuela Basin, Sta. CT-8, 15º00’00” N, 67º05’00” W, and Sta. CT-17, 14º55’00” N, 64º50’00” W. |
9 | Cervigón and Marcano (1965) | Cruises (1962-1965), M/S La Salle. Cariaco Basin, north of the Sucre State (Venezuela), north of Trinidad, Gulf of Paria and Venezuelan Atlantic, 10º20’38” to 12º02’22” N, 63º13’03” to 65º26’02” W. |
10 | Owre and Foyo (1967) | Owre and Foyo (1964) update. Venezuela Basin, Sta. CT-8, 15º00’00” N, 67º05’00” W. Sta. CT-17, 14º55’00” N, 64º50’00” W. |
11 | MHI Archive (1970) | Cruise AV03, R/V Akademik Vernadskiy. Colombian Guajira, Bonaire, Cariaco and Venezuela Basin, 10º37’48” to 16º00’00” N, 65º35’60” to 73º04’48” W. |
12 | Owre and Foyo (1972) | Cruise P-6602, R/V John Elliott Pillsbury. Venezuela Basin, Sta. Nº 14, 15º00’00” N, 64º01’00” W. Venezuelan Atlantic, Sta. Nº 31, 10º01’00” N, 58º10’00” W. |
13 | Fleminger and Hulsemann (1974) | R/V Thomas Washington. Cumaná (Sucre State), Venezuela, Sta. Nº 8, 10º29’00” N, 64º12’00” W. |
14 | Zoppi de Roa (1977) | Cruises 1968-1970, R/V Gulf of Cariaco. Between Margarita Island, north of the Sucre State and Gulf of Paria, Venezuela, Sta. 6, 10º28’00” N, 62º18’00” W. Sta. M, 10º43’00” N, 63º07’00” W. |
15 | Pineda-Polo (1979) | February 1979, R/V USNS Bartlett. New species description. Cariaco Basin, Venezuela, Sta. 10º30’00” N, 64º40’00” W. |
16 | Ferrari and Bowman (1980) | Cruise 1977, R/V Alpha Helix. Only Oithonidae. Carupano harbor (Sucre State, Venezuela), Sta. PN-13-60, 10º41’12” N, 63º14’48” W. Netherlands Antilles, Bonaire, Sta. PN-14-60, 12º10’12” N, 68º18’12” W. Aruba, Sta. PN-18-60, 12º30’18” N, 70º02’42” W. |
17 | Walter (1989) | National Museum of Natural History collection, Smithsonian Institution (USNM) and others researchers. Only Pseudodiaptomidae. Trinidad and Venezuela, Sta. PP, 10º53’00” N, 64º06’00” W. Sta. TT, 10º36’19” N, 61º28’08” W. Sta. IC, 10º45’50” N, 63º58’00” W. |
18 | Zoppi de Roa and Palacios-Cáceres (2005) | Project “Estudio Geológico - Ambiental de la Fachada Atlántica”. Cruise October 2001, R/V ARBV “Punta Brava”. Venezuelan Atlantic, 08º48’12” to 10º16’36” N, 57º42’18” to 60º33’54” W. |
19 | Zoppi de Roa et al. (2007) | Project “Línea Base Ambiental Plataforma Deltana”. Cruises (2004-2005), R/V Hermano Ginés. Gulf of Paria and Venezuelan Atlantic, 09º14’03” to 10º22’29” N, 59º31’34” to 62º28’21” W. |
20 | Rojas-Márquez (2008) | CARIACO project, La Salle Foundation of Natural Sciences. Cruises 2006-2007, R/V Hermano Ginés. Single station, Cariaco Basin, Sta. 10º30’00” N, 64º40’00” W. |
21 | Morales (2014) | Campaigns 2012-2013. Single station, Guaracayal Depression, Gulf of Cariaco (Sucre State, Venezuela), Sta. 10º28’30” N, 63º58’00” W. |
22 | Casanova (2017) | Cruise Hydro-Oceanografic Isla de Aves I-2009, R/V ARBV “Punta Brava”. Isla de Aves, Venezuela, 15º39’50” to 15º40’29” N, 63º36’46” to 63º37’28” W. |
23 | Márquez-Rojas et al. (2021) | Campaigns 2003-2004. Only Corycaeidae. Gulf of Cariaco (Sucre State, Venezuela), 10º27’34” to 10º33’09” N, 64º04’44” to 64º13’02” W. |
24 | Márquez et al. (2021) | Campaigns 2009-2010. Gulf of Cariaco (Sucre State, Venezuela), 10º28’18” to 10º30’00”N, 63º39’50” to 63º43’42” W. |
25 | Casanova et al. (2021) | Campaign 2000. Los Roques Archipelago, Venezuela, 11º44’51” to 11º56’00” N, 66º35’00” to 67º54’20” W. |
26 | Camisotti and Pérez (2021) | Project “Estudio Geológico - Ambiental de la Fachada Atlántica”. Cruise October 2001, R/V ARBV “Punta Brava”. Venezuelan Atlantic, 08º48’12” to 10º16’36” N, 57º42’18” to 60º33’54” W. |
27 | Orrell and Informatics Office (2022) | May, 04, 1962. NMNH Extant Specimen Records (USNM, US). Single station, North Atlantic Ocean, Caribbean Sea. Sta. 4 BSP, 13º01’48” N, 65º28’12” W. |
28 | Segovia-Kancev et al. (2022) | Campaigns Mach 2016. Mochima Bay (Sucre State, Venezuela), 10º22’40” to 10º23’57” N, 64º20’38” to 64º20’45” W. |
Taxonomic considerations
The copepod nomenclature and classification used were based on Walter and Boxshall (2021). The following considerations were taken into account during the list preparation: 1) 3 taxa with open nomenclature at the genus level were included: Saphirella sp., Oithona sp. 1 and 2, whose taxonomic status could not be established; 2) Ditrichocorycaeus affinis (McMurrich) is designated here as D. cf. affinis according to Márquez-Rojas et al. (2014b), awaiting a detailed revision of the material available in the collection; 3) the SCA records of the genus Saphirella consist of larval stages of the family Clausidiidae (Walter & Boxshall, 2021); their taxonomic position should be clarified in the future; 4) the current position of Paracalanus pygmaeus (Claus) and Oncaea gracilis (Dana) is uncertain, so they are included in the list with reservations; 5) Due to the difficulty in establishing its affinities with other species, Pachos punctatum (Claus) is provisionally included in Cyclopoida incertae sedis (Walter & Boxshall, 2021); 6) the families Clausidiidae, Corycaeidae, Lubbockiidae, Oncaeidae and Sapphirinidae belonging to the poecilostome lineage, are included in the order Cyclopoida following Khodami et al. (2017) and Walter and Boxshall (2021).
Species of the family Caligidae are not considered in this review, because it is a group with parasitic forms of fish that are often found in plankton samples (Ohtsuka et al., 2018). Additional data on the Caligidae diversity of the SCA can be obtained from the works by Ho and Bashirullah (1977), González et al. (1986), Lagarde (1989), Díaz (2000), Suárez-Morales et al. (2012), and Kim et al. (2019).
Other components of copepod diversity excluded from our analysis are the benthic and symbiotic forms. For example, this is the case of the siphonostomatoid Tychidion guyanense Humes, a symbiont associated with deep-water annelids from the Atlantic shelf (Humes, 1973) and the cyclopoid Lobosomatium enigmaticum Stock, a parasite of small polychaetes in Curaçao (Stock, 1995). Reid (1990) presented a detailed checklist that included benthic and demersal copepods present in the SCA (Aruba, Bonaire, Curaçao, Trinidad and Tobago, and coasts of South America). The benthic copepods diversity of the order Harpacticoida recorded in the Caribbean Sea was evaluated by Suárez-Morales et al. (2006).
In order to compare the number of species by family present in the SCA, which largely coincide with those documented in the TNWA (Venezuela, Caribbean Sea, Gulf of Mexico, Florida and Sargasso Sea), we used information from Razouls et al.’s (2021) online database.
Estimation of copepod diversity
Species richness is a widely used parameter to describe communities, thus providing a quick idea of the local alpha diversity in the area (Chiarucci et al., 2011; van der Spoel, 1994). In order to analyze the SCA pelagic copepod diversity, we used the program EstimateS v.9.1.0 (Colwell, 2013). From the selected references, the stations visited during each oceanographic cruise (Fig. 2) were established as the unit of sampling effort (n = 613). In this way, it was possible to prepare a dataset with the species incidence (Supplementary material) to estimate the species rarefaction curve, using 100 randomizations without replacement. We also integrate data from Wilson (1942), Bowman (1957) (Table 1, references 1, 2) and Owre and Foyo (1964, 1967) (references 8, 10), because the records come from the same samples. The Clench model (Soberón & Llorente, 1993) was used to calculate the fitted function of the accumulated curve, as:
where S represents the species richness and n is the sampling effort. Values for a and b were obtained by non-linear estimation using the Simplex & quasi-Newton algorithm based on geometrical procedures to reduce the loss function. The values were calculated in the program STATISTICA v.10 (StatSoft, 2011). For the Clench model, the first derivative represents the slope of the curve: S’ = (a) / [(1 + b * n)] 2 . While the ratio (a/b) is the asymptote of the curve (Soberón & Llorente, 1993), so the species proportion recorded can be represented as:
A complete survey of large marine areas is usually logistically infeasible, so rare species are frequently missed (Branco et al., 2018). One solution is to use species richness estimators (Chiarucci et al., 2011). Expected richness was estimated from the Bootstrap and second-order Jackknife non-parametric models, based on the rare species incidence such as “uniques” and “duplicates” (Colwell et al., 2004; Smith & van Belle, 1984). These estimators have the advantage of making no assumptions about the species abundance distribution, which reduces the bias (number of missing species) and thus allowing inferences about the least possible number of species present in the community (Branco et al., 2018).
Results
We determined that the SCA pelagic copepod community is represented by 5 orders, 42 families, 114 genera, 346 species and 2 subspecies (Appendix). The order Calanoida is the most species-rich (244 spp.), followed by Cyclopoida (88), and Harpacticoida (11). Only 3 species are reported of the orders Siphonostomatoida and Monstrilloida.
At the family level, the largest number of species found is contained in Corycaeidae (30 spp.), followed by Paracalanidae (26), Aetideidae (25), Scolecithricidae (25), and Augaptilidae (24). The most diverse genera in the SCA ecoregion are Oithona (19 spp.), Sapphirina (13), Calocalanus (12), Candacia (11), and Euaugaptilus (11). We recorded 2 subspecies: Rhincalanus cornutus atlanticus Schmaus and Oncaea venusta venella Farran, both reported in the SCA only once in different works (Owre & Foyo, 1967; MHI Archive, 1970).
From the literature data it was possible to assess the frequency in which the species were reported (Appendix). In the current inventory, only 40 species have been mentioned in more than 10 publications and can be considered frequent elements of the SCA pelagic copepod community. Eight of these species belong to the family Corycaeidae; Corycaeus speciosus is the species with the highest number of records, mentioned in 22 of the regional works consulted. It is followed by Clausocalanus arcuicornis and Temora stylifera, each with 20 reports. Up to 141 species (40.8% of our inventory) are linked to single reports.
This work includes 11 new records of copepods for the SCA ecoregion (Table 2). Most of the species not hitherto recorded were found at station 63, between La Tortuga and La Blanquilla islands. Among the new findings, the mesopelagic Mimocalanus cultrifer Farran, Arietellus setosus Giesbrecht, and Pseudoamallothrix ovata (Farran) had the highest number of specimens. The depth at which they were obtained is indicated for each record; most were found at depths greater than 200 m. Only Euchaeta pubera Sars and M. cultrifer were found above 100 m.
Specie | Specimen | Station | Date | Depth (m) |
---|---|---|---|---|
Arietellus setosus Giesbrecht, 1893 | 1 ♂ | 63 | 1968-03-28 | 200 |
1 CV | 63 | 1967-04-10 | 200 | |
1 CV | 63 | 1968-01-31 | 1,000 | |
1 ♂ | 76 | 1967-11-25 | 200 | |
Candacia elongata (Boeck, 1872) | 1 ♀ | 63 | 1967-12-08 | 500 |
Candacia tenuimana (Giesbrecht, 1889) | 1 ♀ | 63 | 1967-08-12 | 500 |
Euchaeta pubera Sars G.O., 1907 | 1 ♀ | 63 | 1967-08-12 | 50 |
Euchirella maxima Wolfenden, 1905 | 1 ♀ | 63 | 1967-08-12 | *1,000-0 |
Mimocalanus cultrifer Farran, 1908 | 7 ♀ | 63 | 1967-08-12 | 100 |
Paraeuchaeta barbata (Brady, 1883) | 1 ♀ | 63 | 1967-08-12 | *1,000-0 |
1 ♀ | 63 | 1968-01-31 | 1,000 | |
Paraeuchaeta bisinuata (Sars, 1907) | 1 ♀ | 63 | 1967-08-12 | *1,000-0 |
Pleuromamma borealis Dahl, 1893 | 1 ♀ | 63 | 1968-01-31 | 1,000 |
Pontoeciella abyssicola (Scott, 1893) | 1 ♀ | 63 | 1967-08-12 | 500 |
Pseudoamallothrix ovata (Farran, 1905) | 1 ♀ | 76 | 1968-01-30 | 200 |
1 ♀ | 76 | 1968-01-30 | 1,000 | |
1 ♀ | 76 | 1968-03-27 | 500 |
The comparison between the number of species by family recorded in the present study that concurs with those documented for the TNWA (i.e., Venezuela, Caribbean Sea, Gulf of Mexico, Florida, Sargasso Sea) is shown in Figure 3. For 27 of the 42 families of copepods recorded in the Southern Caribbean, the number of species occurring in both areas exceeds 50%. Of these, the species of 12 families are represented by all those documented in the TNWA province. In the case of both Corycaeidae and Paracalanidae, the number of species documented in this work is greater than that reported for the TNWA province.
From the information compiled, it was possible to model the species rarefaction curve (Fig. 4A). With the selected data (n = 613), the greatest estimated species richness was S = 346.0 ± 10.5 (mean ± SD) which, as expected, coincides with observed species richness (346). With Clench’s model we obtained the fitted curve (R2 = 0.9593, n = 613) whose estimated asymptote (a/b) reaches 368.7 species. Thus, we determined that the theoretical observed richness (q%) represents 93.9% of the inventory of species reported in the SCA. The value of the calculated slope S’ is 0.06 species by sampling unit.
The most conservative estimator of richness was the Bootstrap, with 387.5 species expected (n = 613), whereas the most extreme value was obtained for Jack 2 (second-order Jackknife) with 511.7 expected species. These values predict that the inventory could increase by 12 to 48% over the currently reported figure (Fig. 4B). The curves of these estimators show a clear increasing tendency, mostly related to rare species incorporation (i.e., “uniques” and “duplicates”) (Fig. 4C). In the case of species that appear only once in the sample (“uniques”), we noticed that their number increases proportionately to the sampling effort (98 spp.), thus representing 28.3% of the biodiversity to be evaluated, whereas the species appearing twice (“duplicates”) stabilize at 30 species and progressively decrease with increasing sampling effort.
Discussion
Diversity of pelagic copepods in the Southern Caribbean (SCA)
The diversity found in this work (346 species and 2 subspecies) represents a substantial increase compared to previous reports made in the SCA ecoregion (Márquez-Rojas et al., 2020; Miloslavich et al., 2005; Rodríguez & Suárez, 2003). This value corresponds to 48% of the 724 copepods species documented by Razouls et al. (2021) in the TNWA province (including Venezuela, western Caribbean Sea, Gulf of Mexico, Florida and Sargasso Sea). The copepod species richness found differs numerically from similar inventories conducted in other ecoregions of the TNWA province. For example, 216 species have been recorded in Florida and adjacent waters (Owre & Foyo, 1967), 201 in the Western Caribbean (Mexican Caribbean) (Suárez-Morales & Gasca, 1998), and 247 in the Southwestern Caribbean (Colombian Caribbean), as reported by Medellín-Mora and Navas (2010), and subsequently updated by Gaviria et al. (2019) and Dorado-Roncancio et al. (2021). For the Gulf of Mexico, the 335 species reported represent a value close to that found here and as expected, the species composition between these ecoregions have similarities (Campos-Hernández & Suárez-Morales, 1994; Suárez-Morales et al., 2009). Following Reid’s (1990) review of the Gulf of Mexico and the Caribbean Sea pelagic copepods, we recognize that these ecoregions share a total of 155 species (44.7%) with the present work. While species reported here in 12 families for the Southern Caribbean coincide with all those documented in the Tropical Northwestern Atlantic (Razouls et al., 2021).
It is recognized that 40 species are often reported, so that, as a whole, they can be considered characteristics of the ecoregion community. Corycaeidae species have been continuously reported in the literature, which includes tropical and subtropical species (Campos-Hernández & Suárez-Morales, 1994; Owre & Foyo, 1967), which are very frequent and sometimes abundant in the eastern region of Venezuela (Cervigón, 1964; Márquez-Rojas et al., 2014a, 2014b). This has been attributed to the great adaptive capacity of these copepods under the changing conditions of coastal areas (Álvarez-Cadena et al., 1998; Björnberg, 1981; Suárez, 1989) and their affinity to tropical water masses (Márquez-Rojas et al., 2014a).
The greatest diversity of marine copepods species in the Atlantic is found in the tropical and subtropical latitudes (Piontkovski & Landry, 2003; Rombouts et al., 2009; Wilson, 1942), a pattern also valid for other pelagic groups like: fish, ostracods, decapods and euphausiids (Angel, 1993). The SCA is considered the area with the highest nutrient input in the TNWA because of the widespread upwelling processes along its coastline (Castellanos et al., 2002; Correa-Ramírez et al., 2020; Okolodkov, 2003; Rueda-Roa et al., 2018; Villamizar & Cervigón, 2017), which adds to the huge nutrient input from the Orinoco and Amazon rivers (Muller-Karger & Varela, 1990; Rueda-Roa et al., 2018). This condition may promote species diversity because of the increased energy flow of these ecosystems (Saiz et al., 2007). While for oceanic environments, Angel (1993) indicates that the high values in the pelagic species numbers are usually associated with oligotrophic conditions. Consistent with this pattern, most of the rare species and new reports in the SCA resulted from studies conducted in oceanic environments and mesopelagic samples (MHI Archive, 1970; Owre & Foyo, 1964; Park, 1970; Wilson, 1942). In these areas, the greatest diversity can be reached at depths close to 1000 m, well below the depth where the ecological processes supporting it take place (Angel, 1993), it may therefore be appropriate to evaluate how these conditions interact and influence the diversity of the pelagic components of the SCA. Unfortunately, studies of the deep-water pelagic fauna of tropical latitudes are still scarce, it is expected that future work in these strata will allow us to expand our knowledge of the marine diversity of the ecoregion.
Quantitative tools application for biodiversity estimation allow to test the accuracy species inventories (Colwell et al., 2004; Soberón & Llorente, 1993; van der Spoel, 1994). They improve the confidence in the methodologies employed and allow proving adequate sample sizes for analysis (Branco et al., 2018). Our analysis is deemed as an adequate approximation of the SCA pelagic copepods diversity, as indicated by the 93.9% of the expected inventory and the low slope value obtained (i.e., species incorporation rate by sampling unit S’ = 0.06). Jiménez-Valverde and Hortal (2003) stated that, in order to consider that the inventory is representative of the species diversity, it is necessary to obtain a slope of less than 0.10 species by sampling unit. It has been hitherto reported that community attributes, such as spatial aggregation or uneven distribution of abundance (both characteristics typical of planktonic communities) tend to weaken the predictive trait of non-parametric estimators (Branco et al., 2018; Jiménez-Valverde & Hortal, 2003; Reese et al., 2014). The values obtained in the present study from non-parametric estimators show that the present copepod inventory for the SCA is underestimated between 12 and 48%. This is shown by the increase in the number of rare species (“uniques” and “duplicates”) that are incorporated to the inventory at increasing sampling efforts. Although another characteristic of pelagic ecosystems is their high proportion of rare species (Angel, 1993), this fact must be weighted when predicting species richness using non-parametric estimators. These rare species can also be considered transient members of the community resulting from local dispersal (McManus & Woodson, 2012; van der Spoel, 1994).
Distribution of copepods in the Southern Caribbean (SCA)
Large Marine Ecosystems (LME) arrangement is based on coastal areas classification and continental shelf (neritic zone) (Spalding et al., 2007, 2012). On the other hand, the biogeographic approaches to the mesopelagic fauna proposed by Sutton et al. (2017) we consider may be more suitable to understand the ecological processes that determine the structure and distribution of the zooplankton community in the SCA. According to this biogeographic classification, the Tropical Northwestern Atlantic (TNWA) province (Greater Caribbean) is part of a larger ecoregion called the Central North Atlantic, bounded between the North Atlantic Drift and the Equatorial Atlantic ecoregions. In contrast to the LME classification, where the SCA ecoregion is geographically well delimited, based on Sutton et al. (2017), it is possible to consider that distribution in some planktonic species present greater geographic amplitude, extending throughout the Central North Atlantic ecoregion, as they are influenced by factors such as the physiology of organisms, and environmental factors such as the water mass transport, temperature and salinity (Angel, 1993; McManus & Woodson, 2012; Rombouts et al., 2009). For example, the South Atlantic (Bradford-Grieve et al., 1999) shares 47.8% of the species reported in this work, despite presenting differentiable environmental characteristics such as oligotrophic waters and complex current circulation patterns (Sutton et al., 2017).
The similarity in species composition with contiguous areas, such as the Equatorial Current, has been previously hypothesized (Björnberg, 1981). This is possible due to the relationship between hydrographic regime patterns between major ocean basins (Angel, 1993). In addition, biogeographic boundaries may be asymmetric/semi-permeable by concerning to specific components of the planktonic community (Sutton et al., 2017). In this sense, the Amallothrix tenuiserrata, Eucalanus elongatus and Pareucalanus attenuatus distribution is documented in different works conducted in the SCA and areas adjacent to the Tropical Northwestern Atlantic (TNWA) province (Cervigón & Marcano, 1965; Legaré, 1964; MHI Archive, 1970; Owre & Foyo, 1964, 1972). For this reason, we consider the inclusion of these species justified as part of the biota of the TNWA. Rendón et al. (2003), suggest that E. elongatus may be an introduced species, indicating that it is present in ship ballast waters in Cartagena Bay, Colombia (Southwestern Caribbean). This statement is in contrast to the records of the species in the province (Casanova et al., 2007; Cervigón & Marcano, 1965; MHI Archive, 1970; Owre & Foyo, 1964, 1967; Wilson, 1942; Zoppi de Roa et al., 2007), which we consider enough to rule out that it is an introduced species and instead consider the species part of pelagic copepods community of the TNWA province.
On the other hand, reports of the following species have been eventual: Calanus propinquus, Calocalanus elongatus, C. gresei, C. kristalli, C. longisetosus, C. neptunus, C. pseudocontractus, Candacia tenuimana, Canthocalanus pauper, Cymbasoma rigidum, Ditrichocorycaeus cf. affinis, D. americanus, D. subtilis, Euchirella formosa, Farranula curta, F. gibbula, Labidocera wollastoni, Lucicutia grandis, Oithona fallax, Oncaea gracilis, Onychocorycaeus pumilus, Paracalanus pygmaeus, Pontella lobiancoi, P. mediterranea, Sapphirina maculosa, Scolecitrichopsis tenuipes (Márquez-Rojas et al., 2014a; MHI Archive, 1970; Owre & Foyo, 1964; Wilson, 1942). Despite these few specific reports, the documented distribution areas for the species are contiguous to the TNWA (Razouls et al., 2021), which may explain their presence in the SCA. The opposite is the case for Ditrichocorycaeus andrewsi (Farran), which is not included in the present list because its known natural distribution is centered in the Western Pacific, with some reports from the Indian and Eastern Pacific (GBIF.org, 2021; Razouls et al., 2021). We do not rule out the species presence in the Gulf of Cariaco (Márquez-Rojas et al., 2020; Morales, 2008), but we await its confirmation in the TNWA or adjacent provinces.
Information gaps and perspectives
The underestimation in the number of species mentioned above is notorious in the documented reports of genera in the SCA that so far have not been associated with species present in the TNWA, as is the case with the benthopelagic and bipolar Bradyidius (Legaré, 1964) and Monstrilla (Cervigón & Marcano, 1965). Legaré (1961) indicates that 8 to 10 undetermined species of the order Harpacticoida have been found in samples, although most of them probably correspond to demersal species.
The families Oithonidae and Paracalanidae stand out due to their species richness and unique reports (Ferrari & Bowman, 1980; MHI Archive, 1970). In particular, identification of these species has represented a difficulty and they correspond to groups in which greater efforts must be made to evaluate their real diversity. Oithona hebes reported by Ferrari and Bowman (1980) has also been recorded in brackish waters of the coastal lagoon of Tacarigua (Zacarías & Zoppi de Roa, 1981), an environment in which 2 or more species of the genus are presumed to be present. Other examples of the taxonomic difficulty associated with Oithona species is evidenced in the work of Cervigón (1963) and more recently by Segovia-Kancev et al. (2022), who reports an aggregate of unidentified species belonging to the genus. Regarding Paracalanus, Reid (1990) indicates that there is a probability that P. indicus has been misidentified in the literature as P. parvus and P. quasimodo. Therefore, we do not rule out the possibility that this may be the case in some of the historical reports and we consider of interest the more detailed evaluation of this genus in the Southern Caribbean. In order to support species presence with eventual reports, we urge the various research groups working in the ecoregion to confirm and incorporate the specimens in biological collections during future activities.
Finally, we consider that knowledge of meso and bathypelagic copepods are in the early stages of development; a study of these communities in the Gulf of Mexico and the Northwestern Caribbean provided 58 new records of planktonic copepods for these ocean basins (Park, 1970), while more recently Dorado-Roncancio et al. (2021) reported 33 new records for the Southwestern Caribbean (Colombian Caribbean). Most research on plankton in the Southern Caribbean (SCA) has focused on the neritic zone east of Venezuela. Therefore, to increase the knowledge of the community and associated ecological processes, it is necessary to extend the research area to the western and oceanic zone of the SCA. This work was intended to cover the deficit in the pelagic copepod information in the SCA, which we consider of major importance to carry out with intention of facilitating additional studies of the group. We hope that this work represents a starting point for documentation of species with few records and that it allows a better understanding of the associated oceanographic processes, due to the great relevance of ecological studies in meso and bathypelagic waters for global changes evaluation.
From the information analyzed, we can conclude that the value we give for pelagic copepods diversity is an adequate reflection of the current documented knowledge. However, as suggested by the richness estimators, it is very likely that the actual diversity of pelagic copepods is underestimated. One element that supports this fact is the high number of rare species, typical of pelagic ecosystems. The most conservative value of expected richness was obtained from the non-parametric Bootstrap estimator, while Jack 2 (second-order Jackknife) expresses the most extreme increase in expected richness. In the pelagic copepod community of the Southern Caribbean, the members of Corycaeidae stand out with a significant number of species that are often reported. We consider that presence of Amallothrix tenuiserrata, Eucalanus elongatus and Pareucalanus attenuatus in the Tropical Northwestern Atlantic is supported by documented reports from the ecoregion.