Introducción
Las variaciones climáticas tienen efectos importantes y diversos sobre los ecosistemas marinos y terrestres (Stenseth et al., 2002; Champagnon, Lebreton, Drummond, & Anderson, 2018). Este es el caso de la temperatura del ambiente pelágico, que ha sido utilizada para comprender algunos procesos biológicos, debido a que su variación puede influir en la estructura trófica, al repercutir en la energía producida por los productores primarios, influyendo, por ejemplo, en la distribución de especies, tamaño poblacional, supervivencia de adultos y la proporción de individuos reproductivos (Jenouvrier, Barbraud, & Weimerskirch, 2005; Frederiksen, Edwardst, Richardson, Halliday, & Wanless, 2006). Por esto la variación en la temperatura del mar podría estar asociada con la abundancia de los consumidores tope del ambiente pelágico como las aves marinas.
El efecto de la variación en las condiciones de la temperatura del mar sobre la dinámica de las poblaciones de aves marinas ha sido reportado frecuentemente en la literatura científica (Barbraud et al., 2011; Grémillet & Boulinier, 2009; Sydeman, Thompson, & Kitaysky, 2012). Por ejemplo, bajo condiciones El Niño, se ha registrado una disminución en el éxito reproductivo, aumento de mortalidad, migraciones en busca de alimento y alteraciones en la diversidad de dieta y régimen reproductivo de algunas poblaciones de aves marinas en el Pacífico Oriental Tropical (Wang & Fiedler, 2006; Ancona, Calixto-Albarrán, & Drummond, 2012; Quillfeldt & Masello, 2013). Estos eventos podrían relacionarse con el aumento en la temperatura del mar, y consecuentemente con el reclutamiento de algunos peces, como Engraulis ringens (Jenyns, 1842) (Barber & Chávez, 1986; Sánchez, Calienes, & Zuta, 2000; Wang & Fiedler, 2006; Bakun et al., 2015), la principal presa de Sula nebouxii (Milne-Edwards, 1882)y Sula variegata (von Tschudi, 1843) en la costa peruana (Zavalaga, Benvenuti, Dall’antonia, & Emslie, 2008; Zavalaga, Halls, Mori, Taylor, & Dell’omo, 2010).
Isla Gorgona es una localidad importante para las aves marinas del Pacifico colombiano, por tener colonias permanentes de anidación y descanso de cuatro especies (Fregata magnificens Mathews, 1914, Pelecanus occidentalis Linnaeus, 1766, Sula leucogasteretesiaca Thayer & Bangs, 1905 y S. nebouxii), cuyas poblaciones han sido monitoreadas mensualmente desde el año 2002. Este monitoreo representa la mayor serie de abundancias mensuales de aves marinas registradas en Colombia. También se ha registrado la presencia de S. variegata, que al igual que S. nebouxii, es considerada como especie migratoria proveniente de la costa peruana. Gracias a este monitoreo se tiene un conocimiento preliminar de la variación poblacional de estas especies y que la isla puede actuar como refugio temporal para S. nebouxii y S. variegata cuando la disponibilidad de alimento en las costas peruanas disminuye (Cadena-López & Naranjo, 2010; Payan, Ortega, & Cuéllar, 2017).
La isla Gorgona fue declarada en 2005 como Área de Importancia para la Conservación de las Aves (AICA Co120) (Franco, Devenish, Barrero, & Romero, 2009; Zamudio, 2014). Además, forma parte del Complejo Marino-Costero Iscuandé-Sanquianga-Gorgona, que es el área más importante de Colombia para aves playeras y marinas tanto migratorias como locales (Ruiz-Guerra, 2011). Sin embargo, esta región se encuentra recurrentemente bajo el efecto de variaciones oceanográficas regionales como El Niño, que puede ocasionar alteraciones en la estructura trófica del ecosistema pelágico circundante (Crick, 2004; Santisteban, Benkman, Fetz, & Smith, 2012), generando probablemente efectos indirectos sobre el número de individuos o éxito reproductivo de la colonia reproductiva de S. leucogaster etesiaca, subespecie endémica y vulnerable del Pacífico tropical húmedo americano (Ospina-Álvarez, 2004; Renjifo, Amaya-Villarreal, Burbano-Girón, & Velásquez-Tibatá, 2016). En este contexto, esta investigación pretende describir la dinámica temporal de la abundancia de las poblaciones de S. leucogaster, S nebouxii y S. variegata en isla Gorgona y evaluar su posible relación con la variabilidad temporal de la temperatura superficial del mar (TSM), en la isla, la región Niño 1+2 y San José en la costa peruana.
Materiales y métodos
Isla Gorgona y el islote de Gorgonilla (2˚55’45”-3˚00’55” N & 78˚09’00’’-78˚14’30” W) se ubican a 30 km del continente, y conforman el área insular emergida más extensa de la plataforma continental del Pacífico colombiano con una extensión de 13.7 km2 (Giraldo, Diaz-Granados, & Gutiérez-Landázuri, 2014). En la isla se encuentran playas coralinas, arenosas y rocosas, acantilados y promontorios rocosos emergidos de origen volcánico, en los cuales predomina la roca desnuda, aunque puede haber presencia de arbustos y herbáceas. Estas formaciones rocosas se encuentran representadas en el sur por las rocas alrededor de Gorgonilla, y en el norte y noroccidente de Gorgona (Giraldo & Valencia, 2012; Zamudio, 2014).
La Zona de Convergencia Intertropical condiciona el régimen climático de la isla, caracterizado por una alta pluviosidad (4 164-8 176 mm anuales). La temperatura media anual del aire es 26 °C y la humedad relativa promedio es superior al 90 %. La época lluviosa transcurre entre mayo y octubre, y es seguida por otra de menos precipitaciones entre diciembre y febrero (Diaz, Pinzón, Perdomo, Barrios, & López-Victoria, 2001; Blanco, 2009). En esta localidad han sido descritos dos periodos oceanográficos opuestos, uno cálido de baja salinidad superficial y con una profundidad promedio de la termoclina (profundidad de la isoterma de 22 °C) de 47 m, entre mayo y diciembre, y uno frío de alta salinidad y con una profundidad promedio de la termoclina de 7.5 m, entre enero y abril (Giraldo, Rodríguez-Rubio, & Zapata, 2008). Este enfriamiento ocurre simultáneamente en el sector oriental y occidental de la isla, y consiste en una característica general del Pacífico Oriental Tropical colombiano, originado por la intensificación del chorro de viento de Panamá durante los primeros meses del año, que genera una lengua fría que desciende desde el norte, sin embargo, la temperatura del mar a 1 m de profundidad varía a lo largo del año entre 27.2 y 28.0 °C (Giraldo, 2008; Giraldo, Valencia, Acevedo, & Rivera, 2014).
Se utilizaron los registros de abundancia recopilados durante el desarrollo del programa de monitoreo de aves marinas que implementa el Parque Nacional Natural Gorgona (PNNG). En este programa se realiza un muestreo visual de aves marinas en isla Gorgona, durante dos días consecutivos de la primera semana de cada mes, desde el 2002 hasta la fecha, en el cual se abarca en cada día un sector de la isla, norte o sur. Cada recorrido se inicia a las 0545 h, a bordo de una embarcación con motor fuera de borda y es realizado por dos observadores entrenados. Durante el muestreo se registra, según la especie, número de adultos, juveniles, machos y hembras. Además de las dos especies de piqueros objeto de monitoreo (Sula leucogaster y S. nebouxii), también se ha registrado la presencia de S. variegata en algunas temporadas. En cada registro se incluye la actividad de los individuos (vuelo o descanso), y en caso de vuelo, solo aquellos que se desplazan en dirección contraria al recorrido de la embarcación. El muestreo de las tres especies de piqueros en isla Gorgona abarca un recorrido de 24.04 km, distribuidos en tres trayectos lineales (Poblado-Horno, Horno-Palmasola y Palmasola-Palmeras) y uno circular (Horno) en la zona norte de la isla, además de un trayecto lineal (Poblado-Tasca) y dos circulares (Gorgonilla y El Viudo) en la zona sur (Payan, 2016).
Los datos fueron estandarizados dividiendo la abundancia (sumatoria de registros en los dos días) por la distancia recorrida, ya que no siempre se visitaron todos los sectores de la isla debido a condiciones climáticas adversas. Los datos recolectados en esfuerzos equivalentes a un solo día de muestreo por mes (un sector de la isla) no fueron considerados para el desarrollo de los análisis numéricos. Para evaluar la variación poblacional de S. leucogaster y S. nebouxii, se utilizaron los registros obtenidos desde el 2002. Para S. variegata se emplearon las abundancias registradas a partir del 2012, debido a que durante los primeros años del monitoreo no se registró la abundancia de esta especie. Además, la población de S. leucogaster fue analizada excluyendo los pollos.
Se seleccionaron tres localidades para construir las series de tiempo de TSM: 1. Isla Gorgona, Colombia, como condición local, 2. San José, Perú, debido a la presencia en isla Gorgona de especies de piqueros migratorias probablemente provenientes de la costa del Perú y ser la localidad más cercana a la mayor colonia de estas aves (Islas Lobos de Afuera y Lobos de Tierra), y 3. Región Niño 1+2, Ecuador-Norte Perú, territorio asociado al evento ENOS. Los registros de la TSM de la región Niño 1+2 fueron obtenidos de la página web del Climate Prediction Center (CPC, 2018), los de San José de la página web del Instituto del Mar de Perú (Grados & Vásquez, 2018) y los de isla Gorgona a partir de las imágenes mensuales satelitales registradas por el sensor MODIS-Aqua del Ocean Biology Processing Group (OBPG), centrada en la posición 3°0’0’’ N & 78°0’0’’ W (OBPG, 2015). Las anomalías térmicas fueron calculadas al substraer el promedio histórico mensual de los periodos 1981-2010, 1991-2012 y 2002-2008, correspondientemente, de la serie temporal mensual de la TSM. Así mismo, se tuvo acceso al Índice de Oscilación Peruano (POI), que explica el 72 % de la variación de las anomalías mensuales de la TSM a lo largo de la costa peruana (Purca-Cuicapusa, 2005).
Para determinar el desfase temporal de la posible influencia de la temperatura sobre las abundancias de los piqueros, se realizó una correlación cruzada de los registros de la abundancia mensual de cada especie con la TSM y las anomalías térmicas de las tres localidades, así como con el POI. En total se realizaron 21 correlaciones cruzadas, 7 para cada especie. El desfase temporal (Lag) de cada correlación fue seleccionado entre 0 y 12 meses, ya que es el mínimo intervalo que abarca completamente la variación de la abundancia de las tres especies de piqueros. La fuerza de la correlación es representada por la función de autocorrelación (ACF), cuyo valor varía entre -1 y 1, siendo -1 la mayor correlación negativa y 1 la mayor correlación positiva. A continuación, se desarrolló un modelo lineal generalizado (GLM) para cada especie de manera independiente, con la finalidad de explicar la posible asociación entre la abundancia de cada piquero y la condición de temperatura que presentó la mayor correlación, considerando el desfase temporal encontrado. Cada modelo fue desarrollado con un nivel de significancia de 0.05. Como se analizaron datos de conteo se utilizó la distribución quasipoisson, y en caso de sobredispersión se implementó la distribución binomial negativa. Todos los análisis y figuras fueron desarrolladas en el programa R (R Core Team, 2018), utilizando los paquetes ggplot2 para las figuras (Wickham, 2016), astsa para la correlación cruzada (Stoffer, 2017), MASS para los modelos GLM (Venables & Ripley, 2002), y funciones de uso general para la organización de los datos (Wickham, 2007; Fox & Weisberg, 2011; Trapletti & Hornik, 2018).
Resultados
Durante el muestreo de aves marinas en isla Gorgona se realizaron registros de 184 abundancias mensuales de Sula leucogaster y S. nebouxii y 81 de S. variegata. La abundancia de S. leucogaster presentó la menor variación intra-anual, con una mediana de 6.68 ind/km por mes, y una abundancia máxima y mínima de 16.35 y 1.33 ind/km, respectivamente. Por otro lado, S. nebouxii fue la especie más abundante, con una mediana de 79.64 ind/km por mes. Esta especie presentó la mayor variación, siendo más abundante a mitad de año, con un máximo de 550.97 ind/km y un mínimo de 0 ind/km en octubre 2008, que representó el único mes en el que no hubo ningún registro de esta especie. A diferencia de las especies anteriores, S. variegata no fue observada durante todos los meses del año, siendo menos abundante en el primer semestre. La mediana mensual fue 0.15 ind/km con un máximo y mínimo de 74.05 y 0 ind/km, respectivamente. Asimismo, se destacó que los cinco mayores registros fueron obtenidos en 2014.
Los registros de la TSM de la región Niño 1+2 (mediana = 23.01 °C) generalmente fueron menores que los de isla Gorgona (mediana = 27.26 °C), aunque presentaron mayor dispersión, fluctuando entre 19.5 y 28.10 °C, mientras que en isla Gorgona las temperaturas oscilaron entre 25.47 y 28.05 °C (Fig. 1). Sula nebouxii y S. variegata presentaron la mayor correlación con la TSM de la región Niño 1+2. Además, el desfase temporal encontrado sugiere que la abundancia de S. nebouxii y S.variegata están asociadas con la temperatura registrada cuatro y cinco meses antes, respectivamente (Tabla 1). Por su parte, S. leucogaster estuvo correlacionada con la TSM de isla Gorgona, sin ningún desfase (Tabla 1).
Temperatura | Sula leucogaster | Sula nebouxii | Sula variegata | |||
ACF | Lag | ACF | Lag | ACF | Lag | |
ANOM Niño 1+2 | -0.173 | -1 | 0.181 | -2 | - | - |
TSM Niño 1+2 | -0.252 | -2 | 0.568 | -4 | 0.309 | -5 |
ANOM Gor | -0.153 | 0 | - | - | - | - |
TSM Gor | -0.252 | 0 | 0.421 | -1 | 0.299 | -2 |
ANOM SJ | - | - | -0.175 | 7 | - | - |
TSM SJ | -0.253 | -2 | 0.427 | -4 | - | - |
POI | -0.284 | -1 | 0.198 | 0 | - | - |
ACF = fuerza de la correlación, ANOM = anomalía térmica, Gor = isla Gorgona, Lag: desfase temporal de la correlación medido en meses, Niño 1+2 = región Niño 1+2, POI = Índice de Oscilación Peruano, SJ = San José, TSM = temperatura superficial del mar.
El modelo lineal generalizado, desarrollado entre la TSM y la abundancia con desfase temporal, sugirió que el número de individuos de S. nebouxii (coeficiente TSM = 1.36, P < 0.001) y de S. variegata (coeficiente TSM = 2.77, P < 0.001) en isla Gorgona tendieron a aumentar con el incremento de la TSM en la región Niño 1+2 (Fig. 2A). En cambio, la abundancia de S. leucogaster (coeficiente TSM = 0.82, P ≤ 0.001), al parecer, tiende a disminuir al aumentar la TSM de isla Gorgona (Fig. 2B).
Discusión
Generalmente, los adultos de piqueros permanecen cerca de la colonia reproductiva mientras hay abundante alimento, lo cual se observa principalmente en aquellas especies que habitan aguas productivas, como el caso de Sula nebouxii y S. variegata, que forrajean en las aguas de la corriente de Humboldt (Ludynia, Garthe, & Luna-Jorquera, 2010; Ancona et al., 2012). Recientemente, Payan et al. (2017) identificaron una correlación positiva entre la abundancia de S. variegata y las anomalías de la TSM de la región Niño 1+2, y sugieren que esta especie pudo haber migrado desde la costa peruana hacia el Pacífico colombiano, como respuesta al efecto de El Niño sobre la productividad de esta región, cómo ha sido propuesto por Sánchez et al. (2000), Wang y Fiedler (2006) y Bakun et al. (2015). Aun así, este estudio no evaluó un posible desfase en la influencia de la TSM sobre la abundancia de los piqueros. Al observar que las temperaturas más altas de la región Niño 1+2 fueron registradas en marzo, mientras que las mayores abundancias de S. nebouxii y S. variegata fueron obtenidas entre julio y agosto, se sugiere que el desfase encontrado en la correlación cruzada podría considerarse como el tiempo que demora la migración de estos individuos desde las costas peruanas hasta isla Gorgona, cuando ocurre un aumento en la TSM en la región Niño 1+2.
Es posible que otras variables estén influyendo en la migración de S. variegata, ya que las mayores abundancias (año 2014) no estuvieron asociadas con los mayores registros de TSM (aproximadamente entre 25 y 26 °C), e incluso, a mayores temperaturas (28 °C) se observaron menos individuos en la isla. Asimismo, por presentar menor abundancia que S. nebouxii, y asumir que ambas especies provienen de la costa peruana, se podría llegar a considerar que S. variegata no sería tan sensible como S. nebouxii ante la variabilidad o el incremento de la TSM, y probablemente se desplaza cuando las condiciones de variación térmica son más drásticas. Cabe mencionar que esta especie se caracteriza por forrajear más cerca de la colonia en comparación con S. nebouxii, por lo que podría realizar migraciones más cortas (Zavalaga et al., 2008, Zavalaga et al., 2010), condición que es sustentada por los reportes de esta especie en Ecuador durante El Niño de 1997 (Haase, 1997).
Sula nebouxii se puede reproducir en cualquier mes del año (Nelson, 1978), por lo que la variación poblacional intra-anual de esta especie en isla Gorgona podría estar más asociada con la fluctuación de la TSM de la región Niño 1+2, que con su ciclo reproductivo. Además, es conocido que la población peruana ha migrado hacia el sur de Perú y norte de Chile durante el fenómeno de El Niño, asociado con la escasez de recursos presa (Cushman, 1926; Schulenberg, Stotz, Lane, O’Neill, & Parker III, 2007), menor acceso a ellas debido a la profundización de la termoclina, o alteraciones en la distribución horizontal y vertical de algunas especies de peces pelágicos (Gilman, Chaloupka, Peschon, & Ellgen, 2016). Se podría asumir que este piquero también migra hacia el norte debido al calentamiento del mar, como ocurrió durante El Niño de 1982-83, en el cual hubo registros de más 7 000 individuos de S. variegata y 10 000 de S. nebouxii en el Golfo de Panamá, (Jaksic & Fariña, 2010) mientras que en islas Galápagos se reportaron pérdidas de nidadas y partida de individuos de S. nebouxii(Jaksic & Fariña, 2010). La permanencia de estos piqueros en el Golfo de Panamá hasta abril de 1985, su condición física, y el desarrollo de actividades reproductivas, sugieren que esta localidad presenta las condiciones alimentarias requeridas para abastecer sus necesidades (Smith, 1990; D’Croz & O’Dea, 2007). Finalmente, cabe mencionar que la distancia recorrida al borde de la costa desde Lobos de Afuera hasta Gorgona (aproximadamente 1 200 km) es menor que la recorrida hasta el norte de Chile (aproximadamente 1 800 km). Esto podría favorecer la migración hacia el norte, al encontrarse un refugio a menor distancia, a pesar de que E. ringens no se encuentre en isla Gorgona, y si esté presente en Chile (Mujica, Nava, Matsuda, & Vargas, 2016).
Con respecto a S. leucogaster, en el modelo se observó que la abundancia disminuye cuando aumenta la TSM mensual, relacionándose posiblemente con la disminución en la productividad, y consecuentemente en la disponibilidad de alimento, cuando aumenta la temperatura del mar (Wang & Fiedler, 2006). Aun así, se desconoce por qué las mayores abundancias no fueron obtenidas durante los primeros meses del año, cuando la productividad es mayor en isla Gorgona (Giraldo, Valencia, Acevedo, & Rivera, 2014), por lo que otras variables podrían influir en la disponibilidad de presas.
En conclusión, se sugiere que la TSM, tanto de isla Gorgona como de la región Niño 1+2, está asociada, posiblemente a través de su efecto sobre la disponibilidad y distribución de peces pelágicos, con las variaciones poblacionales de los tres piqueros presentes en isla Gorgona. Además, se destaca la importancia de esta localidad para la conservación de los piqueros, al albergar una colonia reproductiva de una subespecie endémica de S. leucogaster, considerada vulnerable y susceptible a posibles variaciones ambientales. Así mismo, isla Gorgona recibe individuos de dos especies de piqueros migratorias, S. nebouxii y S. variegata, que al parecer encuentran en esta localidad, alimento y sitio de descanso cuando las condiciones ambientales en la costa peruana no son adecuadas por el incremento de la TSM del mar y la consecuente reducción de la disponibilidad de presas. Sin embargo, considerando la menor abundancia de S. variegata en isla Gorgona, se sugiere que esta especie estaría migrando cuando las condiciones adversas son más intensas.
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.