Introducción
Conocer los atributos ecológicos de los ensamblajes de larvas de peces ha tomado mayor relevancia recientemente, toda vez que son una etapa crítica en la historia natural de un recurso alimentario fundamental para la humanidad (Flores-Coto et al., 2009; Hamilton et al., 2021; Llopiz et al., 2014; Schilling et al., 2020). Particularmente, el conocimiento de las fases iniciales del ciclo de vida de los peces posee implicaciones directas para la gestión, administración y control de los recursos pesqueros, adquiriendo un valor significativo en la protección y explotación sostenible de estos recursos (Harrison et al., 2012; Hsieh et al., 2016; Kough et al., 2013; Saville & Schnack, 1981), ya que puede utilizarse para evaluar el tamaño de la población (Crowder et al., 2018; Van der Lingen & Huggett, 2003), así como para identificar las zonas (Kurogi et al., 2012; Kuroki et al., 2020; Tsukamoto, 2006) y temporadas de desove (Pritt et al., 2015; Sassa & Hirota, 2013; Somarakis et al., 2019). Además, conocer la distribución de las larvas de peces en áreas y periodos de tiempo, permite comprender la movilidad espacio-temporal que presentan algunas especies, actividad que se encuentra íntimamente ligada con el hábito reproductivo y puede variar en escalas temporales que van desde mensuales, semestrales o anuales, e inclusive con el ciclo de marea o el ciclo lunar (Gallego & Giraldo, 2018; Miller & Kendall, 2009).
Diversas investigaciones han sugerido que los ensamblajes de larvas de peces son el resultado de estrategias de desove convergentes de múltiples especies que aprovechan condiciones favorables en el entorno que permitan incrementar la probabilidad de supervivencia de las larvas (Doyle et al., 1993; Sherman et al., 1983), llegando a estar la variabilidad espacial y temporal de estas agrupaciones relacionadas con el comportamiento de las larvas (Gray & Miskiewicz, 2000), su distribución en la columna de agua (Olivar & Sabatés, 1997), la disponibilidad de presas (Olivar et al., 2010), la estrategia de desove (Moser & Smith, 1993) y el comportamiento migratorio de peces adultos (Hare et al., 2001). En estas etapas iniciales de desarrollo, las larvas de peces son muy sensibles a los cambios en las condiciones hidrológicas, debido a su poca capacidad natatoria, siendo también vulnerables a la inanición y a la depredación (China & Holzman, 2014; Downie et al., 2020), por lo tanto, los factores ambientales pueden llegar a tener una fuerte influencia en la estructura y composición de los ensamblajes de larvas de peces (de Macedo-Soares et al., 2014; Houde, 2008; Keskin & Pauly, 2019; León-Chávez et al., 2010; León-Chávez et al., 2015; Whitfield & Pattrick, 2015). En este sentido, procesos oceanográficos como el transporte advectivo o los procesos de mezcla en la columna de agua (Auth, 2008; Muhling et al., 2007), así como los patrones de circulación, los frentes termales o los procesos de surgencia (Govoni & Pietrafesa, 1994; Sabatés et al., 2007), se han considerado como mecanismos que pueden llegar a promover la retención, concentración o dispersión de los individuos, aumentando o disminuyendo la supervivencia de las larvas (Kiorboe et al., 1988) y por ende modulando la estructura de los ensamblajes y la dinámica de sus poblaciones.
En un ciclo anual han sido ampliamente descritos los cambios en la composición y abundancia de los ensamblajes de larvas de peces (Calle-Bonilla et al., 2017; Flores-Coto et al., 2000; Gallego et al., 2018; Jiménez-Rosenberg et al., 2007; Leis & McCormick, 2002; León-Chávez et al., 2015; Mavruk et al., 2018; Ribeiro et al., 2015). Aunque la magnitud de estas variaciones intra-anuales dependen de la distancia al ecuador, en el Pacífico Oriental Tropical (POT), la variación intra-anual en los patrones climáticos afecta de manera drástica las condiciones oceanográficas regionales y locales (Alongi, 1990; Fernández-Álamo & Färber-Lorda, 2006; Pennington et al., 2006), lo que estaría modulando la dinámica temporal de las comunidades biológicas que allí habitan. En esta región, el desarrollo de procesos de surgencias estacionales influenciadas por el desplazamiento anual de la Zona de Convergencia Intertropical (ZCIT) (Amador et al., 2006), se convierte en un factor determinante en la variabilidad de las condiciones ambientales del ambiente pelágico (Fernández-Álamo & Färber-Lorda, 2006; Pennington et al., 2006). Particularmente, en el área colombiana del POT, se han descrito dos periodos oceanográficos contrastantes como consecuencia del desarrollo de procesos de surgencia costero en la Ensenada de Panamá: un periodo de aguas frías, entre enero y abril, en donde ocurre un afloramiento de aguas frías y saladas provenientes de capas subsuperficiales debido a giros ciclónicos impulsados por corrientes de viento, que conlleva a una disminución de la temperatura de la superficie, un aumento en la salinidad y nutrientes, y una mayor abundancia de plancton, y un periodo de aguas cálidas, que va de mayo a diciembre, en donde la columna de agua está más cálidas y menos productivas, y con baja salinidad superficial debido a las intensas precipitaciones en el área (Corredor-Acosta et al., 2020; Giraldo et al., 2008; Torres et al., 2023; Valencia et al., 2019).
En esta investigación se describe la estructura y composición del ensamblaje de larvas de peces del Golfo de Tribugá, el cual se localiza en la zona norte del Pacífico Colombiano, y es una región de gran importancia para el aprovisionamiento de recursos biológicos a través de actividad pesquera de subsistencia y comercial (Neira et al., 2016; Tobón-López et al., 2008). Específicamente, se evaluó la variación temporal del ensamblaje de larvas en dos periodos oceanográficos contrastantes que han sido descritos para esta zona (Valencia et al., 2013), a partir de las siguientes preguntas: (1) ¿Cuál es la identidad taxonómica de las larvas de peces presentes en el golfo de Tribugá? (2) ¿Cuál es la abundancia de las larvas presentes en cada periodo oceanográfico descrito para la región? (3) ¿Varía la abundancia de larvas de peces en el Golfo de Tribugá entre los periodos oceanográficos contrastantes de la región?, y (4) ¿Cambia la estructura de los ensamblajes de larvas de peces del Golfo de Tribugá entre los periodos oceanográficos de región? Se planteó como hipótesis que la estructura, composición y abundancia del ensamblaje de larvas en cada periodo oceanográfico es significativamente diferentes, de tal manera que la abundancia será mayor durante la temporada de influencia de las aguas de surgencia y se presentará un recambio de especies superior al 50 %.
Materiales y métodos
Área de Estudio: El Golfo de Tribugá es una amplia escotadura en la costa norte del Pacífico colombiano que corresponde al sector más meridional de la Ecorregión Marino-Costera colombiana “Pacífico Norte” (Díaz & Acero, 2003). Abarca un área aproximada de 103 110 ha y se extiende desde Cabo Corrientes (5º53´ N & 77º09´ W), en el sur, hasta la desembocadura del río Valle (6º11´ N & 77º09´ W), en jurisdicción de los municipios de Nuquí y Bahía Solano (Zuluaga, 2010) (Fig. 1). Esta área ofrece una zona litoral y nerítica con una alta variedad de hábitats, entre los que se destacan las zonas de bosque de manglar, planos lodos, fondos arenosos, sustratos rocosos, acantilados y formaciones coralinas (Prahl et al., 1990).

Fig. 1 Ubicación geográfica de las estaciones de muestreo en el Golfo de Tribugá, Pacífico norte de Colombia. / Fig. 1. Geographic location of the sampling stations in the Gulf of Tribugá, northern Pacific coast of Colombia.
Recolección de muestras: Para la recolección de muestras, se estableció una malla de muestreo compuesta por 18 estaciones, distribuidas sistemáticamente en el golfo. Cada línea de estaciones comenzaba aproximadamente a 1 km de la costa, con una distancia de 2 km entre cada estación (Fig. 1). En total fueron recolectadas 72 muestras de zooplancton en dos campañas oceanográficas realizadas durante el periodo de agua fría (marzo) y agua cálida (septiembre) de 2022. En cada estación, se realizaron registros verticales de temperatura y salinidad utilizando un perfilador CastAway®. Además, se tomaron muestras discretas de agua a cuatro profundidades (1, 10, 30 y 50 m) utilizando una botella Niskin en las que se registró el pH y oxígeno disuelto con una sonda multiparamétrica Hanna® HI9829 y se estableció la concentración de clorofila-a in vivo y la turbidez con un fluorómetro portátil AquaFlour®, previamente calibrado. La recolección de muestras de ictioplancton se realizó mediante arrastres oblicuos de 20 m a superficie durante 5 min, utilizando una red cilíndrica-cónica de plancton de 71 cm de diámetro con malla de 500 µm, equipada con medidor de flujo Hydrobios® para estimar el volumen de agua filtrada. Las muestras colectadas fueron conservadas en formaldehído buferizado al 4 % y transportadas al laboratorio de Ecología Animal de la Universidad del Valle.
Procesamiento de muestras e identificación taxonómica: Se separaron e identificaron las larvas de peces de cada muestra hasta el nivel taxonómico más bajo posible siguiendo la nomenclatura taxonómica establecida por Nelson (2006). Para la confirmación taxonómica se utilizó el atlas 33 de Moser (1996), el atlas de Richards (2005) y bibliografía especializada (Jiménez-Rosenberg et al., 2003; Jiménez-Rosenberg et al., 2006; Matsumoto, 1959; Miyashita et al., 2001; Ortíz-Galindo et al., 2008). Además, cuando fue necesario se realizaron comparaciones con especímenes de referencia depositados en las colecciones científicas de larvas de peces del Centro de Investigación en Ciencias del Mar (La Paz, México CICIMAR IPN) y de la Colección de Prácticas Zoológicas de la Universidad del Valle (Cali, Colombia).
Análisis de datos: Se describieron las condiciones oceanográficas registradas durante cada periodo y se utilizó una prueba no paramétrica de Mann-Whitney (MW-U) para comparar estas condiciones entre los dos periodos oceanográficos establecidos. Se calculó la abundancia relativa por familias y por especie de larvas de peces para cada uno de los periodos de estudio y se evaluaron diferencias entre periodos utilizando una prueba MW-U. Se utilizaron los índices de atributos comunitarios de Shannon y Margalef para describir la diversidad y riqueza de los ensamblajes, y se estableció el recambio de especies a partir de la composición de los ensamblajes de larvas de peces durante cada periodo oceanográfico utilizando el índice de Whittaker (Bw) como:
Bw = (S/a)/N
Donde:
S es el número total de especies.
a es el número promedio de especies
N es el número de muestras, considerando que el esfuerzo de muestro fue igual durante los dos periodos de estudio.
Para evaluar las diferencias en la estructura comunitaria de larvas de peces entre periodos, se empleó el análisis multidimensional no paramétrico (nMDS) y se determinó el aporte de las especies estimando el porcentaje de similaridad (Simper). Finalmente, se realizó un análisis de correspondencia canónica (ACC) utilizando el software libre Rstudio v. 4.1.2 (R Core Team, 2021) para evaluar el grado de relación entre la abundancia de larvas de peces y los parámetros fisicoquímicos de la columna de agua en el Golfo de Tribugá.
Resultados
Condiciones ambientales: El promedio de la temperatura superficial del mar de 1 a 10 m, durante el periodo cálido, fue de 28.2 °C, mientras que la salinidad promedio fue de 24.8 PSU. Para la profundidad de 10 a 30 m, la temperatura promedio se mantuvo similar a la superficial, siendo de 28.0 °C, mientras que la salinidad promedio aumentó, siendo de 29.2 PSU. Por otro lado, en el periodo frío, la temperatura promedio de 1 a 10 m fue de 26.4 °C, con una salinidad promedio de 31.6 PSU. Mientras que, de 10 a 30 m, la temperatura promedio fue de 23.3°C y la salinidad promedio fue de 33.6 PSU. La variabilidad en las variables oceanográficas permitió observar diferencias significativas entre los periodos para la temperatura superficial del mar (MW-U: 44, p < 0.05), la salinidad (MW-U: 324, p < 0.05), el oxígeno disuelto (MW-U: 311, p < 0.05), el pH (MW-U: 0, p < 0.05) y la clorofila-a (MW-U: 308, p < 0.05) (Tabla 1).
Tabla 1 Parámetros físico-químicos registrados en el Golfo de Tribugá, Pacífico norte de Colombia, durante los periodos Frío y Cálido de 2022. / Table 1. Physicochemical parameters recorded in the Gulf of Tribugá, northern Pacific coast of Colombia, during the Cold and Warm periods of 2022.
| Frío | Cálido | W | p | |||||||
|---|---|---|---|---|---|---|---|---|---|---|
| D.s. | Min | Max | D.s. | Min | Max | |||||
| Temperatura (°C) | 27.54 | 0.54 | 27.18 | 29.37 | 28.14 | 0.32 | 27.67 | 28.79 | 44 | < 0.05 |
| Salinidad (UPS) | 30.48 | 0.78 | 29.42 | 33.05 | 24.20 | 0.46 | 23.22 | 25.12 | 324 | < 0.05 |
| Oxígeno disuelto (ml/l) | 4.73 | 0.29 | 3.97 | 5.15 | 3.86 | 0.5 | 2.58 | 4.44 | 311 | < 0.05 |
| pH | 8.15 | 0.023 | 8.07 | 8.16 | 8.19 | 0.009 | 8.18 | 8.21 | 0 | < 0.05 |
| Clorofila-a (mg/m3) | 0.86 | 0.33 | 0.38 | 1.62 | 0.36 | 0.14 | 0.27 | 0.81 | 308 | < 0.05 |
μ: promedio, D.s.: desviación estándar, Min: valor mínimo registrado. Max: valor máximo registrado. W: valor de Wilcoxon, p: valor de probabilidad. / μ: mean, SD: standard deviation, Min: minimum recorded value, Max: maximum recorded value, W: Wilcoxon value, p: probability value.
Abundancia larval: Se capturaron e identificaron 602 larvas de peces durante el periodo de estudio (199 larvas durante el periodo frío y 403 larvas durante el periodo cálido). En la mayoría de las estaciones de muestreo se capturaron larvas de peces en ambos periodos de estudio. Durante el periodo de agua fría la densidad de larvas de peces estuvo entre 12.49 larvas/1 000 m3 y 464.05 larvas/1 000 m3, mientras que durante el periodo de agua cálida la densidad de larvas de peces estuvo entre 53.49 larvas/1 000 m3 y 725.03 larvas/1 000 m3, estableciéndose diferencias significativas en la densidad de larvas peces entre los periodos de muestreo (MW-U: 74, p < 0.05) (Fig. 2A).

Fig. 2 Variación temporal de: A. Densidad de larvas de peces (larvas/1 000 m3), B. Diversidad (H’), y C. Riqueza del ictioplancton, durante los periodos Frío y Cálido del 2022, en el Golfo de Tribugá, Pacífico norte de Colombia. / Fig. 2. Temporal variation of: A. Fish larvae density (larvae/1 000 m³), B. Diversity (H’), and C. Ichthyoplankton richness, during the Cold and Warm periods of 2022, in the Gulf of Tribugá, northern Pacific of Colombia.
Composición taxonómica: Durante el estudio, del total de larvas recolectadas y separadas, el 2.5 % no pudieron ser identificadas debido a su deterioro. Las larvas taxonómicamente identificadas, fueron clasificadas en 29 familias y 89 especies (Tabla 2). Durante el periodo oceanográfico frío, se registraron 43 especies, mientras que en el periodo cálido se identificaron 64 especies, con 18 especies comunes entre ambos periodos. A lo largo de todo el estudio, la familia con el mayor número de especies fue Carangidae, con un total de 12 especies identificadas, seguida por Sciaenidae con 11, Gerreidae con 9 y Engraulidae con 8 especies. De las 89 especies registradas, 15 fueron identificadas solo hasta el nivel de familia, 41 hasta género y 33 hasta especie.
Tabla 2 Composición taxonómica del ensamblaje de larvas de peces del Golfo de Tribugá los periodos Frío y Cálido de 2022. / Table 2. Taxonomic composition of the fish larvae assemblage in the Gulf of Tribugá during the Cold and Warm periods of 2022.
| Familia | Especie | Densidad promedio (Larvas/1 0003) | AR (%) | ||
| Frío | Cálido | Frío | Cálido | ||
| Nemichthyidae | Nemichthys sp. | 0 | 10.51 | 0 | 0.22 |
| Engraulidae | Anchoa spp. | 292.11 | 12.31 | 11.77 | 0.26 |
| Anchoa sp.1 | 111.94 | 36.19 | 4.51 | 0.77 | |
| Anchoa sp.2 | 33.58 | 124.07 | 1.35 | 2.65 | |
| Anchoa sp.3 | 40.48 | 93.51 | 1.63 | 2.00 | |
| Anchoa sp.4 | 0 | 21.02 | 0 | 0.45 | |
| Anchoa sp.5 | 0 | 42.20 | 0 | 0.90 | |
| Anchoa sp.6 | 16.57 | 31.11 | 0.67 | 0.66 | |
| Anchoa sp.7 | 0 | 86.89 | 0 | 1.86 | |
| Phosichthyidae | Vinciguerria lucetia (Garman, 1899) | 12.29 | 0 | 0.50 | 0 |
| Synodontidae | Synodus sp.1 | 12.29 | 0 | 0.50 | 0 |
| Myctophidae | Benthosema panamense (Tåning, 1932) | 61.12 | 0 | 2.46 | 0 |
| Lampanyctus parvicauda Parr, 1931 | 10.09 | 0 | 0.41 | 0 | |
| Bythitidae | Bythitidae sp.1 | 0 | 23.16 | 0 | 0.50 |
| Mugilidae | Mugil curema Valenciennes, 1836 | 11.81 | 0 | 0.48 | 0 |
| Mugil sp.1 | 0 | 12.49 | 0 | 0.27 | |
| Exocoetidae | Prognichthys sp. | 11.73 | 0 | 0.47 | 0 |
| Carangidae | Carangidae sp. | 33.15 | 0 | 1.34 | 0 |
| Caranx sp.1 | 11.74 | 0 | 0.47 | 0 | |
| Caranx sp.2 | 0 | 21.25 | 0 | 0.45 | |
| Chloroscombrus orqueta Jordan & Gilbert, 1883 | 10.09 | 0 | 0.41 | 0 | |
| Decapterus sp.1 | 0 | 10.46 | 0 | 0.22 | |
| Decapterus sp.2 | 0 | 173.15 | 0 | 3.70 | |
| Decapterus sp.3 | 0 | 10.51 | 0 | 0.22 | |
| Hemicaranx sp.1 | 0 | 10.46 | 0 | 0.22 | |
| Oligoplites saurus inornatus (Bloch & Schneider, 1801) | 33.15 | 134.95 | 1.34 | 2.88 | |
| Oligoplites refulgens Gilbert & Starks, 1904 | 0 | 168.53 | 0 | 3.60 | |
| Selar crumenophthalmus (Bloch, 1793) | 11.74 | 0 | 0.47 | 0 | |
| Selene peruviana (Guichenot, 1866) | 13.47 | 10.46 | 0.54 | 0.22 | |
| Lutjanidae | Lutjanus sp.1 | 9.70 | 10.51 | 0.39 | 0.22 |
| Lutjanus sp.2 | 0 | 35.42 | 0 | 0.76 | |
| Lutjanus sp.3 | 10.22 | 30.62 | 0.41 | 0.65 | |
| Lutjanus sp.4 | 20.19 | 108.60 | 0.81 | 2.32 | |
| Lobotidae | Lobotes pacifica Gilbert, 1898 | 0 | 77.89 | 0 | 1.66 |
| Gerreidae | Diapterus peruvianus (Cuvier, 1830) | 0 | 9.70 | 0 | 0.21 |
| Diapterus sp.1 | 0 | 34.23 | 0 | 0.73 | |
| Diapterus sp.2 | 0 | 65.72 | 0 | 1.40 | |
| Eucinostomus gracilis (Gill, 1862) | 0 | 254.14 | 0 | 5.43 | |
| Eucinostomus entomelas Zahuranec, 1980 | 10.22 | 0 | 0.41 | 0 | |
| Eugerres lineatus (Humboldt, 1821) | 0 | 1265.36 | 0 | 27.04 | |
| Gerreidae sp.1 | 0 | 115.34 | 0 | 2.47 | |
| Gerreidae sp.2 | 0 | 11.91 | 0 | 0.25 | |
| Gerreidae sp.3 | 0 | 21.91 | 0 | 0.47 | |
| Haemulidae | Haemulidae sp.1 | 43.93 | 59.09 | 1.77 | 1.26 |
| Haemulidae sp.2 | 0 | 10.46 | 0 | 0.22 | |
| Haemulidae sp.3 | 13.47 | 14.80 | 0.54 | 0.32 | |
| Haemulidae sp.4 | 0 | 14.80 | 0 | 0.32 | |
| Haemulidae sp.5 | 0 | 23.00 | 0 | 0.49 | |
| Pomadasys sp.1 | 0 | 10.85 | 0 | 0.23 | |
| Sparidae | Calamus brachysomus (Lockington, 1880) | 22.42 | 0 | 0.90 | 0 |
| Sciaenidae | Cynoscion sp.1 | 0 | 40.05 | 0 | 0.86 |
| Cynoscion sp.2 | 0 | 44.39 | 0 | 0.95 | |
| Cynoscion sp.3 | 0 | 13.28 | 0 | 0.28 | |
| Isopisthus remifer Jordan & Gilbert, 1882 | 0 | 20.00 | 0 | 0.43 | |
| Larimus sp.2 | 19.22 | 87.01 | 0.77 | 1.86 | |
| Macrodon mordax (Gilbert & Starks, 1904) | 0 | 12.49 | 0 | 0.27 | |
| Micropogonias altipinnis (Günther, 1864) | 0 | 35.00 | 0 | 0.75 | |
| Stellifer illecebrosus Gilbert, 1898 | 27.07 | 20.00 | 1.09 | 0.43 | |
| Stellifer sp.3 | 0 | 12.49 | 0 | 0.27 | |
| Sciaenidae sp.1 | 0 | 12.49 | 0 | 0.27 | |
| Umbrina sp.1 | 16.57 | 0 | 0.67 | 0 | |
| Pomacentridae | Abudefduf troschelii (Gill, 1862) | 0 | 176.33 | 0 | 3.77 |
| Azurina atrilobata (Gill, 1862) | 46.53 | 0 | 1.88 | 0 | |
| Chromis sp. | 0 | 24.62 | 0 | 0.53 | |
| Labridae | Halichoeres sp. | 33.15 | 12.31 | 1.34 | 0.26 |
| Blenniidae | Hypsoblennius sp.1 | 0 | 11.58 | 0 | 0.25 |
| Ophioblennius steindachneri Jordan & Evermann, 1898 | 0 | 23.82 | 0 | 0.51 | |
| Labrisomidae | Paraclinus sp.1 | 0 | 11.81 | 0 | 0.25 |
| Gobiesocidae | Gobiesocidae sp. | 90.20 | 188.93 | 3.64 | 4.04 |
| Gobiesox sp.1 | 0 | 160.06 | 0 | 3.42 | |
| Gobiidae | Bathygobius andrei (Sauvage, 1880) | 13.54 | 0 | 0.55 | 0 |
| Bathygobius lineatus (Jenyns, 1841) | 23.50 | 0 | 0.95 | 0 | |
| Ctenogobius manglicola (Jordan & Starks, 1895) | 0 | 11.26 | 0 | 0.24 | |
| Lythrypnus sp.1 | 13.47 | 0 | 0.54 | 0 | |
| Microdesmidae | Microdesmus sp.1 | 10.22 | 0 | 0.41 | 0 |
| Microdesmus sp.2 | 11.21 | 0 | 0.45 | 0 | |
| Microdesmus sp.3 | 13.54 | 0 | 0.55 | 0 | |
| Ephippidae | Chaetodipterus zonatus (Girard, 1858) | 18.70 | 0 | 0.75 | 0 |
| Sphyraenidae | Sphyraena ensis Jordan & Gilbert, 1882 | 26.18 | 22.32 | 1.05 | 0.48 |
| Sphyraena sp.1 | 11.74 | 23.82 | 0.47 | 0.51 | |
| Sphyraena sp.2 | 0 | 22.09 | 0 | 0.47 | |
| Scombridae | Auxis rochei (Risso, 1810) | 901.11 | 0 | 36.32 | 0 |
| Auxis thazard (Lacepède, 1800) | 0 | 131.93 | 0 | 2.82 | |
| Thunnus sp.1 | 0 | 86.40 | 0 | 1.85 | |
| Nomeidae | Cubiceps pauciradiatus Günther, 1872 | 201.00 | 0 | 8.10 | 0 |
| Paralichthyidae | Etropus ectenes Jordan, 1889 | 13.47 | 0 | 0.54 | 0 |
| Etropus peruvianus Hildebrand, 1946 | 13.47 | 0 | 0.54 | 0 | |
| Achiridae | Achirus scutum (Günther, 1862) | 0 | 14.80 | 0 | 0.32 |
| Tetraodontidae | Sphoeroides sp.1 | 0 | 192.81 | 0 | 4.12 |
Densidad: abundancia estandarizada por el volumen de agua filtrada. AR: abundancia relativa. / Density: abundance standardized by the volume of filtered water. AR: relative abundance.
En el periodo frío, las familias con mayor abundancia relativa fueron: Scombridae (41.92 %), Engraulidae (16.67 %) y Nomeidae (8.08 %), mientras que en el periodo cálido destacaron Gerreidae (37.72 %), Carangidae (10.92 %) y Engraulidae (9.18 %). Respecto a las especies más abundantes, durante el periodo frío se observó una prevalencia significativa de Auxis rochei (36.32 %), Anchoa spp. (11.77 %) y Cubiceps pauciradiatus (8.10 %), mientras que en el periodo cálido predominaron Eugerres lineatus (27.04 %), Eucinostomus gracilis (5.43 %) y Sphoeroides sp.1 (4.12 %) (Tabla 2).
Según el índice de diversidad de Shannon-Wiener, el periodo de aguas cálidas presentó un mayor valor medio de diversidad (1.90 bits/ind) en comparación con el periodo de aguas frías (1.12 bits/ind), evidenciando una variación significativa en la diversidad entre los periodos de estudio (MW-U: 21.5, p < 0.05) (Fig. 2B). La riqueza, medida mediante el índice de Margalef, fue de 0.76 en el periodo frío y de 1.64 en el cálido, estableciéndose diferencias significativas para la riqueza entre los periodos (MW-U: 23.5, p < 0.05), consistente con la tendencia observada en la diversidad (Fig. 2C).
Se establecieron diferencias significativas en la estructura de los ensamblajes de larvas de peces (ANOSIM, R = 0.6695, p = < 0.05) entre los dos periodos de estudio, registrando una mayor variabilidad en la composición y estructura durante el periodo frio que durante el periodo cálido (Fig. 3). El recambio de especies (BW) entre los dos periodos de estudio fue del 66.33 %, siendo Eugerres lineatus, Auxis rochei y Cubiceps pauciradiatus las especies que más contribuyeron en esta diferenciación temporal (Tabla 3).
Tabla 3 Resumen del análisis Simper de las especies de ictioplancton presentes en el Golfo de Tribugá. / Table 3. Summary of the Simper analysis of ichthyoplankton species present in the Gulf of Tribugá.
| Especies | Frío | Cálido | Promedio de contribución (%) | Contribución Acumulada (%) |
| Densidad promedio | Densidad promedio | |||
| Eugerres lineatus | 0 | 70.30 | 10.95 | 10.95 |
| Auxis rochei | 50.06 | 0 | 7.78 | 18.74 |
| Cubiceps pauciradiatus | 11.17 | 0 | 3.63 | 22.36 |
| Oligoplites refulgens | 0 | 9.36 | 3.43 | 25.79 |
| Eucinostomus gracilis | 0 | 14.12 | 3.37 | 29.16 |
| Decapterus sp.2 | 0 | 9.62 | 2.82 | 31.99 |
| Oligoplites saurus inornatus | 1.84 | 7.50 | 2.54 | 34.53 |
| Anchoa sp.2 | 1.87 | 6.89 | 2.52 | 37.06 |
| Sphoeroides sp.1 | 0 | 10.71 | 2.52 | 39.57 |
| Larimus sp.2 | 1.07 | 4.83 | 2.47 | 42.04 |
| Auxis thazard | 0 | 7.33 | 2.40 | 44.44 |
| Anchoa sp.1 | 6.22 | 2.01 | 2.40 | 46.84 |
| Abudefduf troschelii | 0 | 9.80 | 2.27 | 49.11 |
| Gerreidae sp.1 | 0 | 6.41 | 2.23 | 51.34 |
| Gobiesox sp.1 | 0 | 8.89 | 2.06 | 53.40 |
| Anchoa spp. | 16.23 | 0.68 | 1.99 | 55.39 |
| Gobiesocidae sp.1 | 5.01 | 10.50 | 1.84 | 57.23 |
| Haemulidae sp.1 | 2.44 | 3.28 | 1.81 | 59.04 |
| Thunnus sp.1 | 0 | 4.80 | 1.77 | 60.81 |
| Lutjanus sp.4 | 1.12 | 6.03 | 1.74 | 62.55 |
| Anchoa sp.3 | 2.25 | 5.19 | 1.71 | 64.26 |
| Benthosema panamense | 3.40 | 0 | 1.64 | 65.90 |
| Lobotes pacifica | 0 | 4.33 | 1.63 | 67.53 |
| Diapterus sp.2 | 0 | 3.65 | 1.38 | 68.90 |
| Anchoa sp.7 | 0 | 4.83 | 1.28 | 70.19 |

Fig. 3 Diagrama de Análisis de escalamiento multidimensional (MDS) de la estructura comunitaria de larvas de peces, durante los periodos Frío y Cálido de 2022, en el Golfo de Tribugá. / Fig. 3. Multidimensional Scaling (MDS) diagram of the fish larvae community structure during the Cold and Warm periods of 2022 in the Gulf of Tribugá.
Los dos primeros ejes de ordenación del análisis de correspondencia canónica (ACC) explicaron el 28.75 % de la varianza total del ensamblaje de larvas de peces en el Golfo de Tribugá. Se realizó una prueba de permutación para evaluar la significancia global del modelo de CCA, y los resultados mostraron que el modelo global es altamente significativo (p < 0.0001). Al evaluar la contribución individual de cada variable ambiental al modelo de CCA se encontró que las variables temperatura (Temp) y salinidad (Sal) fueron significativas (p < 0.0001), mientras que Chl-a mostró una tendencia a la significancia (p = 0.05), lo que indica que la variabilidad en la composición de especies está influenciada por variaciones en la temperatura y salinidad, y que podría haber una influencia de la clorofila-a en la estructura comunitaria evaluada (Fig. 4).

Fig. 4 Análisis de correspondencia canónica utilizando la abundancia de las especies encontradas, durante los periodos Frío y Cálido de 2022, en el Golfo de Tribugá. / Fig. 4. Canonical Correspondence Analysis using the abundance of species found during the Cold and Warm periods of 2022 in the Gulf of Tribugá.
Discusión
Las variaciones temporales en las condiciones ambientales presentaron características típicas de aguas del Pacífico Oriental Tropical (POT), con aguas superficiales cálidas (> 25 °C) y de baja salinidad (< 31 PSU). Además, debido a que el Pacífico colombiano se encuentra dentro de la Ensenada de Panamá, las salinidades en esta región son aún más bajas ocasionadas por el exceso de precipitación (Fiedler & Lavín, 2017; Fiedler & Talley, 2006), condición que fue particularmente evidente en el área de estudio durante el periodo de aguas cálidas, donde las salinidades superficiales y la transparencia de la columna de agua fueron más bajas debido a la descarga de los ríos Jurubidá, Baudó, Docampadó, Catripe, Pilizá y Purricha, cobrando este dato vital importancia, puesto que la salinidad del agua también es un factor ambiental significativo que puede afectar la distribución y fisiología del ictioplancton (Houde, 2008).
Se registró una variabilidad estacional en la temperatura del agua en el golfo de Tribugá, similar a la reportada en otras localidades del Pacífico colombiano (Giraldo, 2008; Giraldo et al., 2008; Jerez-Guerrero et al., 2017; Valencia et al., 2019), condición que puede llegar a tener importantes implicaciones para la reproducción, el crecimiento y la supervivencia de las larvas de peces, ya que la temperatura influye en la tasa de desarrollo de los huevos y larvas, así como en la disponibilidad de alimento (Houde, 2008). Además, los resultados obtenidos también sugieren que la productividad de la columna de agua en el Golfo de Tribugá parece estar más influenciada por los procesos de surgencia estacional que afectan a la Ensenada de Panamá durante la estación seca (enero-marzo). En este sitio, la Chl-a fue más alta en el periodo frío, debido al afloramiento estacional provocado por la formación de la surgencia de Panamá, similar a lo reportado por D’Croz y O’Dea (2007) para las aguas de la plataforma en el Golfo de Panamá y por Valencia et al. (2019) para la Ensenada de Utría, localizada en el límite norte del Golfo de Tribugá.
En términos generales, la abundancia y composición de especies de larvas de peces en el golfo de Tribugá exhibió una variación estacional asociada con la variabilidad oceanográfica registrada en la zona de estudio, con un recambio temporal de especies del 66 %. Esta tendencia de variación estacional del ictioplancton es similar a la reportada para otros ambientes marino-costeros tropicales, donde la composición de los ensamblajes de larvas de peces responde de manera muy estrecha a la variabilidad de las condiciones ambientales locales, con taxones y grupos distintos que muestran máximos de abundancia interanuales (Aceves-Medina et al., 2003; Aceves-Medina et al., 2008; de Macedo-Soares et al., 2014; Franco-Gordo et al., 2002; Funes et al., 2002; Giraldo, 2020; Houde, 2008; Keskin & Pauly, 2019; León-Chávez et al., 2010; León-Chávez et al., 2015; Valencia et al., 2024; Whitfield & Pattrick, 2005).
Es relevante destacar que la mayor diversidad taxonómica y abundancia de larvas de peces se registró durante el periodo de aguas cálidas en el golfo de Tribugá. Esto podría estar relacionado con que las características oceanográficas durante este periodo fueron más favorables para el desarrollo larval. Aunque se esperaba que la surgencia y el afloramiento de aguas frías incrementaran la productividad primaria, y, por ende, la disponibilidad de alimentos para las larvas, factores como la estabilidad de la columna de agua, la temperatura y la disponibilidad de nutrientes en la capa superficial pueden haber favorecido una mayor abundancia y diversidad durante el periodo cálido (Bakun, 1996; Pattrick et al., 2021). En ambientes tropicales, es común observar que las aguas cálidas, aunque menos ricas en nutrientes, pueden ofrecer condiciones estables y temperaturas óptimas para el desarrollo de huevos y larvas de peces (Longhurst & Pauly, 1987; Pauly et al., 2023). Se ha documentado que un incremento en la temperatura del agua aumenta la probabilidad de supervivencia de las larvas, ya que el agua presenta una menor viscosidad, lo que mejora la locomoción de las larvas, incrementando su habilidad para escapar de depredadores y facilitando su alimentación (Miller & Kendall, 2009). Además, la combinación de temperaturas moderadas y altas crea ventanas ambientales que favorecen el crecimiento de las larvas de peces. Así, es indiscutible que la variación temporal de las condiciones ambientales modula la distribución, abundancia y, en general, la estructura de los ensambles de larvas en el ambiente pelágico (Aceves-Medina et al., 2003; Aceves-Medina et al., 2008; Franco-Gordo et al., 2002; Funes et al., 2002).
Otro factor relevante que podría explicar estos resultados es la fenología reproductiva de las especies. Es posible que un mayor número de especies se reproduzcan y desoven durante el periodo cálido, lo que explicaría la mayor riqueza y abundancia observadas (Sundby, 2000; Zapata-Londoño et al., 2020). Algunas especies podrían estar adaptadas a condiciones específicas de temperatura y estratificación de la columna de agua, favoreciendo su desarrollo larval en estos periodos (Longhurst & Pauly, 1987).
En cuanto al número de taxa encontrados en el Golfo de Tribugá (29 familias y 89 especies) y las familias con mayor abundancia relativa (Scombridae, Engraulidae y Nomeidae durante el periodo frío y Gerreidae, Carangidae y Engraulidae durante el periodo cálido) fueron consistentes con los resultados reportados por otros estudios de ictioplancton realizados en diferentes localidades del Pacífico colombiano. Por ejemplo, Giraldo (2020) señaló la presencia de 43 familias y 183 especies en Cabo Manglares, destacando la importancia de familias como Engraulidae, Gerreidae, Scombridae y Carangidae. Gallego-Zerrato et al. (2023) reportaron 46 familias y 88 especies en Isla Gorgona, siendo las familias más abundantes Engraulidae, Bregmacerotidae, Haemulidae, Myctophidae y Carangidae, estando el ensamblaje dominado por especies pelágico-costeras (Bregmaceros bathymaster), pelágica (Cetengraulis mysticetus) y batipelágica (Diaphus pacificus). Martínez-Aguilar et al. (2010) reportaron un total de 68 especies pertenecientes a 37 familias para la zona nerítica del Pacífico colombiano, destacando por su abundancia familias como Gobiidae y Engraulidae, cada una el 5 % del total de individuos. Además, familias como Scombridae, Gerreidae y Nomeidae también mostraron importancia numérica. La dominancia en número de estas familias podría estar asociada a patrones ontogénicos de migración vertical.
El Golfo de Tribugá ha sido ampliamente reconocido como una localidad con una alta diversidad de fauna íctica (Tobón-López et al., 2008) y amplia oferta de servicios ecosistémicos (Rojas et al., 2019), entre los que se destaca el aprovisionamiento a través de la pesca (Estupiñán-Cardona et al., 1990; Gómez, 1992; Tobón, 2001; Vieira, 2003). Sin embargo, los diferentes esfuerzos de ordenamiento y manejo de recursos ícticos desarrollados en esta región adolecen de información sobre las fases larvales de las especies de peces presentes, incluidas aquellas especies que son de interés para la actividad pesquera artesanal, la cual es una de las principales actividades económicas para los pobladores de la región (Neira et al., 2016). Poder integrar la información sobre los ensamblajes de larvas de peces en una localidad con los estudios de dinámica poblacional de los estadios adultos y la caracterización de las condiciones oceanográficas propias de cada sitio de interés es requerido para comprender la dinámica espacio-temporal de la diversidad íctica en una localidad, y permitiría incluso fortalecer las acciones de manejo, conservación y ordenamiento de uso de recursos marinos de interés para las comunidades humanas.
Declaración de ética: los autores declaran que todos están de acuerdo con esta publicación y que han hecho aportes que justifican su autoría; que no hay conflicto de interés de ningún tipo; y que han cumplido con todos los requisitos y procedimientos éticos y legales pertinentes. Todas las fuentes de financiamiento se detallan plena y claramente en la sección de agradecimientos. El respectivo documento legal firmado se encuentra en los archivos de la revista.











uBio 


