DISTRIBUCIÓN ESTACIONAL DE CETÁCEOS EN LAS ISLAS CANARIAS CETACEAN SEASONAL DISTRIBUTION IN CANARY ISLANDS Máster en Biología Marina: Biodiversidad y Conservación Talía Morales Herrera

NATACHA AGUILAR DE SOTO, investigadora postdoctoral del Departamento de Biología Animal, Edafología y Geología de la Universidad de La Laguna, y FERNANDO LOZANO, profesor titular del Departamento de Biología Animal, Edafología y Geología de la Universidad de La Laguna, CERTIFICAN: Que la memoria presentada por la Licenciada en Biología Talía Morales Herrera, titulada “Análisis de la distribución estacional de cetáceos en las islas Canarias”, ha sido realizada bajo nuestra dirección y consideramos que reúne todas las condiciones de calidad y rigor científico requeridas para optar a su presentación como Trabajo de Fin de Máster, en el Máster de Biología Marina: Biodiversidad y Conservación de la Universidad de La Laguna, curso 2013-2015. Y para que así conste y surta los efectos oportunos, firmamos el presente certificado en San Cristóbal de La Laguna, a 1 de junio de 2015

Fdo.: Dra. Natacha Aguilar de Soto

Fdo.: Dr. Fernando Lozano

ii

ÍNDICE

INTRODUCCIÓN JUSTIFICACIÓN ....................................................................................................................................................... 1 SITUACIÓN Y CONTEXTUALIZACIÓN. LOS CETÁCEOS EN CANARIAS .......................................................... 2

METODOLOGÍA ......................................................................................................................................... 6 METODOLOGÍA EN LA ISLA DE EL HIERRO ....................................................................................................... 8 METODOLOGIA EN LAS ISLAS DE LA PALMA, GRAN CANARIA, TENERIFE Y LA GOMERA ....................... 9 METODOLOGÍA USADA EN CETAVIST: RED DE AVISTAMIENTO ARCHIPELÁGICA .................................. 10 METODOLOGÍA UTILIZADA EN LA CAMPAÑA ZIPHYSETER ......................................................................... 11 ANÁLISIS ESTADÍSTICOS DESCRIPTIVOS ......................................................................................................... 12 MODELOS ESTADÍSTICOS ARIMA: SERIES TEMPORALES DE DATOS........................................................... 13

RESULTADOS RESULTADOS GENERALES ................................................................................................................................. 15 DISTRIBUCIÓN ESTACIONAL DE LA PRESENCIA DE ESPECIES .................................................................... 16 ANÁLISIS COMPARATIVO: MODELOS ESTADÍSTICOS ARIMA ...................................................................... 20

DISCUSIÓN FACTORES QUE DETERMINAN LA DISTRIBUCIÓN DE CETÁCEOS EN CANARIAS...................................... 39 DISTRIBUCIÓN ESTACIONAL DE CETÁCEOS ODONTOCETOS ....................................................................... 41 DISTRIBUCIÓN ESTACIONAL DE CETÁCEOS MISTICETOS ............................................................................ 48

CONCLUSIONES ...................................................................................................................................... 52 ANEXO I ............................................................................................................................................................ 53 AGRADECIMIENTOS .......................................................................................................................... 54

iii

Resumen Las características geomorfológicas de las Islas Canarias, así como la influencia de la corriente de Canarias y del up-welling sahariano, hacen del archipiélago el lugar con mayor biodiversidad de cetáceos de las áreas vecinas (n=30), incluyendo especies residentes, transeúntes y migratorias. Usan el área con fines de alimentación, reproductivos, de descanso o como corredor migratorio, sin embargo, la escasez de datos sobre su fenología y distribución dificulta el diseño de medidas eficaces centradas en su conservación. Además, la necesidad de medios costosos que permitan la aproximación a sus áreas de distribución dificulta el seguimiento y control de las poblaciones. Datos procedentes de campañas de investigación dirigida, así como de plataformas de oportunidad, como embarcaciones de observación comercial y ferris de pasajeros, han sido recolectados con el fin de analizar el patrón estacional que presentan las distintas especies de cetáceos en las islas. Modelos estadísticos tipo ARIMA, basados en series temporales de datos, han sido aplicados en un análisis comparativo entre dos especies con diferente ocurrencia temporal. Nuevos datos de especies intensamente estudiadas han sido descubiertos, así como la preferencia de determinadas especies por áreas concretas del archipiélago. Junto a esto, fenómenos oceanográficos ligados al aumento de temperatura del agua han podido ser vinculados con cambios en el comportamiento temporal de grandes cetáceos. Estudios que combinen los muestreos científicos sistemáticos con la potente herramienta de obtención de información adicional que suponen las plataformas de oportunidad, deben ser prioritarios para conocer el efecto que el cambio climático, y otros impactos, tienen sobre la distribución espacio-temporal de los cetáceos. Palabras clave: diversidad cetológica, plataformas de oportunidad, patrón estacional, cambio climático.

Abstract The geomorphological features of the Canary Islands, as well as the influence of the Canary current and saharan up-welling, make the archipelago the place with de most cetacean diversity of neighboring areas (n = 30), including resident, vagrant and migratory species. They use the area for feeding, breeding, resting or like a migration corridor, however, insufficient data on phenology and distribution difficult to design effective measures focusing on their conservation. In addition, the need for costly means to approach to their ranges difficults to monitor and control populations. Data from research directed campaigns and from opportunistic platforms as commercial vessels and passenger ferries, have been collected in order to analyze the seasonal pattern that various species of cetaceans present in the islands. Statistical models based on time series data have been applied in a comparative analysis of two species with different temporal occurrence. New data from intensively studied species have been discovered, as well as the preference of certain species for specific areas of the archipelago. Alongside this, oceanographic phenomena linked to increasing water temperatures have been related to changes in the temporal behavior of large cetaceans. Studies combining scientific systematic sampling with the powerful tool like opportunistic platforms for obtaining further information, should be a priority to know the effect that climate change and other impacts have on the spatio-temporal distribution of cetaceans. Keywords: cetacean diversity, opportunity platforms, seasonal pattern, climate change.

iv

INTRODUCCIÓN Justificación La región macaronésica se erige como un enclave extraordinario para la fauna vertebrada, con una enorme heterogeneidad de hábitats (Brito et al. 2001). Esta diversidad de ambientes se traduce en una elevada riqueza de especies entre las que destacan los mamíferos marinos. Éstos, y concretamente los cetáceos, constituyen uno de los grupos animales que por su desarrollada inteligencia, su estructura y hábitos sociales y sus complejas estrategias de supervivencia,

entre otras razones,

despiertan una gran expectación en la comunidad

científica y en la sociedad en general. Sin embargo, existe un déficit de información en cuanto a su distribución espacio temporal, el uso del hábitat, la existencia y ubicación de corredores migratorios o su ecología trófica. La interacción entre humanos y cetáceos ha experimentado grandes cambios, desde el desarrollo de la industria ballenera, que perdura en la actualidad en algunas regiones, hasta su mantenimiento en cautividad con cuestionados fines comerciales y terapéuticos. Sin embargo, en las últimas décadas estamos asistiendo a un cambio en la forma de aproximación a estos animales gracias al esfuerzo conjunto de las administraciones e investigadores y a la labor de entidades ligadas a movimientos conservacionistas. A esto, se une la implantación y consolidación de actividades eco-turísticas como la observación de cetáceos en libertad desde embarcaciones

comerciales,

que

ha

experimentado un espectacular crecimiento. Se

estima

que

actualmente

aproximadamente diez millones de personas participan

cada

año

en

viajes

de

avistamiento de delfines y ballenas en todo el mundo (Day, 2006). En este contexto, Canarias se establece como destino europeo preferente por sus bondades climáticas y la gran diversidad de especies de cetáceos que alberga,

siendo Tenerife la isla que

Fig.1 Embarcación de observación comercial de cetáceos "Nashira Uno" durante un avistamiento de rorcual tropical (Balaenoptera edeni brydei) en Los Gigantes, Tenerife (Mayo, 2014). Foto: Talía Morales.

concentra el 65% de la actividad (Elejabeitia & Urquiola, 2009). No es de extrañar que el sector científico haya puesto sus ojos en estas y otras plataformas oportunistas como forma de 1

obtención de datos adicionales, que contrarresta el elevado coste que supone el estudio en el mar y la aproximación a sus áreas de distribución, a menudo, lejos de costa. Sin embargo, es este mismo sector científico el que alerta de severos problemas de conservación y de la existencia de graves amenazas que ponen en peligro determinadas poblaciones. La presión que generan ésta y otras actividades turísticas, sobre el hábitat y los individuos (Aguilar et al. 2002), el ruido marino (Aguilar & Brito, 2002), la escasez de alimento, la captura accidental durante maniobras pesqueras (Félix et al.2006), la elevada contaminación marina, las colisiones con embarcaciones (Stephanis & Urquiola, 2006), etc. suponen fuertes impactos que han de ser estudiados en profundidad para conocer el estado de conservación en el que se encuentran. Conocer cuándo, dónde y por cuánto tiempo se dan las distintas especies en Canarias se establece como una herramienta prioritaria de gestión, fundamental para diseñar medidas centradas en su conservación. Es conocido que las distintas especies de cetáceos muestran diferentes patrones de ocurrencia temporal dependiendo de sus necesidades y uno de los factores que en mayor medida estimula estos movimientos es la disponibilidad de alimento (Lusseau et al., 2004; Keiper et al., 2005; Mercuri, 2007). Este estudio aborda la distribución temporal de las distintas especies de cetáceos presentes en el archipiélago, haciendo uso de los datos generados tanto por campañas científicas dirigidas como por plataformas de oportunidad, que engloban tanto a los barcos de avistamiento comercial de cetáceos como a los ferris de pasajeros que conectan vía marítima las Islas Canarias.

Situación y contextualización. Los cetáceos en Canarias. El océano Atlántico y concretamente sus archipiélagos nororientales, albergan una enorme diversidad de cetáceos, sin embargo, son las Islas Canarias las que cuentan con mayor riqueza de las regiones circundantes (Brito et al., 2001) con 30 citadas. Este fenómeno puede ser explicado por varias razones, entre ellas, las características climáticas y oceanográficas. El archipiélago es bañado por la corriente fría de Canarias, rama oriental de la corriente del Golfo, que en su recorrido hacia el suroeste confiere a las aguas insulares una temperatura más baja de las que les corresponde por latitud (Tejedor, 2014). Este fenómeno hace posible que Canarias sea el punto de encuentro de especies típicas de aguas templadas que tendrían aquí su límite sur de distribución y especies con apetencias más tropicales, siendo Canarias su 2

límite norte (Aguilar, 2006). A este fenómeno térmico se une el efecto de los vientos alisios, que tienen un importante papel en la llegada de filamentos de agua fría y ricos en nutrientes desde la costa sahariana (Hernández-Guerra et al., 1993; Arístegui et al., 1994; Arístegui et al., 1997; Hernández-Guerra y Nykjaer, 1997; Pacheco y Hernández-Guerra, 1999; Hernández-Guerra et al., 2001). Debido a este evento, conocido como afloramiento o upwelling, se establece un gradiente térmico, con temperaturas más frías en el sector oriental, próximo al continente africano y más cálidas en el sector occidental, con diferencias de hasta 2° C entre ambos extremos. Este fenómeno tiene importantes repercusiones en la distribución longitudinal de las especies marinas en general y de cetáceos en particular. Por su parte, las características geomorfológicas son especialmente relevantes en las islas con mayor altitud ya que suponen un obstáculo en el recorrido de vientos y corrientes. Es el llamado efecto masa-isla, que crea zonas de turbulencias y zonas de calma o remanso (Gómez, 1991) y como consecuencia, se dan fenómenos oceanográficos mesoescalares, como remolinos ciclónicos y anticiclónicos localizados (Hernández-Guerra et al., 1993) o pequeños afloramientos de agua profunda en las desembocaduras de los barrancos, lo que condiciona la distribución parcheada de nutrientes en un medio eminentemente oligotrófico y que a su vez repercute en la distribución espacial y temporal de la especies marinas pelágicas entre los que se encuentran los cetáceos. La batimetría de las islas también hace posible la presencia conjunta de especies neríticas y especies oceánicas. La práctica inexistencia de plataforma permite que se alcancen elevadas profundidades muy próximas a la costa, lo que determina la presencia de un elevado número de especies llamadas de “buceo profundo”, que destacan por su abundancia (Carrillo & Tejedor, 2004) y que en otras regiones del mundo se encuentran bastante alejadas de tierra. Esto es especialmente reseñable en las islas más jóvenes, donde los procesos eruptivos son recientes, y los de erosión y relleno apenas se han producido, generando un acusado desnivel en distancias muy cortas. Las especies presentes en aguas canarias (tabla 1) incluyen las residentes, como es el caso de las poblaciones de calderón de aleta corta (Heimlich-Boran, 1993) bien estudiada en el suroeste de Tenerife así como la de delfín mular (Carrillo & Peña, 2004), el calderón gris o delfín de Risso con poblaciones en el suroeste y noreste de Gran Canaria, Lanzarote y Fuerteventura así como norte de La Gomera y Tenerife (Martín et al., 2011) o el caso de los zifios de Cuvier y de Blainville (Aguilar de Soto et al., 2005; Aparicio et al., 2006) en la isla de El Hierro. A éstas se une el cachalote, presente todo el año (André, 1998). Existen además poblaciones estables de las que se desconoce su distribución temporal como el delfín de dientes rugosos y el delfín moteado del Atlántico. De otras especies, se sospecha que tienen 3

en las islas un lugar de estancia temporal como el delfín común. Poco se sabe de la orca, con avistamientos puntuales y cuya distribución no está bien conocida especialmente en aguas al sur del Estrecho de Gibraltar (Martín, 2009). Miembros de la familia Balaenopteridae están bien documentados en Canarias. El rorcual tropical con un patrón de visitas regular está entre los más frecuentes de los misticetos (Aguilar, 2006) mientras que el rorcual común, el azul, el boreal, el aliblanco o la yubarta, ocasionales en las islas, transitan el archipiélago en sus rutas con fines reproductivos o de alimentación. Tabla 1. Especies de cetáceos citadas para Canarias. De derecha a izquierda: suborden, familia, nombre científico y nombre común. En negrita, las especies registradas en este estudio.*Especies descritas como residentes en Canarias. **Especies con única cita por varamiento. ? Especies registradas pero de dudosa identificación. SUBORDEN

FAMILIA

SUBORDEN ODONTOCETOS (CETÁCEOS DENTADOS)

Delphinidae

Physeteridae Kogiidae

Ziphiidae

SUBORDEN MISTICETOS (CETÁCEOS CON BARBAS

Phocoenidae Balaenopteridae

Balaenidae

NOMBRE CIENTÍFICO Tursiops truncatus* Steno bredanensis Stenella frontalis Stenella coeruleoalba Stenella longirostris Stenella attenuata? Delphinus delphis Lagenodelphis hosei Grampus griseus* Globicephala macrorhynchus* Globicephala melas? Orcinus orca Pseudorca crassidens Physeter macrocephalus Kogia breviceps Kogia sima Mesoplodon densirostris* Mesoplodon europaeus Mesoplodon bidens** Mesoplodon mirus** Ziphius cavirostris* Hyperoodon ampullatus Phocaena phocaena** Balaenoptera musculus Balaenoptera physalus Balaenoptera borealis Balaenoptera edeni brydei Balaenoptera acutorostrata Megaptera novaengliae Eubalaena glacialis**

NOMBRE COMÚN Delfín mular Delfín de dientes rugosos Delfín moteado del atlántico Delfín listado Delfín acróbata Delfín tropical Delfín común Delfín de Fraser Delfín de Risso Calderón de aleta corta Calderón de aleta larga Orca Falsa orca Cachalote Cachalote pigmeo Cachalote enano Zifio de Blainville Zifio de Gervais Zifio de Sowerby Zifio de True Zifio de Cuvier Zifio calderón boreal Marsopa común Rorcual azul Rorcual común Rorcual norteño Rorcual tropical Rorcual aliblanco Yubarta Ballena franca

4

Estos patrones de presencia sin embargo, podrían estar sufriendo un proceso de reajuste presumiblemente potenciado por los cambios térmicos registrados en las últimas décadas. Las aguas del archipiélago canario han experimentado un aumento progresivo en su temperatura que alcanza los 0,28°C/década (Vélez et al.,com.pers.) Este proceso de calentamiento podría estar eliminando barreras naturales que hasta entonces suponían fronteras infranqueables para determinadas especies, en el caso de los peces, (Brito et al., 2005) o determinando la distribución y disponibilidad de alimento. Las series temporales de datos persiguen explicar la evolución y el comportamiento de determinados valores a lo largo del tiempo así como realizar predicciones futuras. El uso de esta herramienta estadística es de enorme utilidad en este estudio, pues se pretende en primer lugar, dilucidar si existe un fenómeno de variación anual repetitiva en la ocurrencia temporal de las distintas especies y en segundo lugar, encontrar un modelo que recoja las variables elegidas y que permita describir y predecir matemáticamente el patrón de presencia estacional de dichas especies.

5

METODOLOGÍA Los datos disponibles para el presente estudio fueron recolectados desde diversas fuentes que incluyen campañas científicas dirigidas, plataformas de observación comercial de cetáceos y una red de avistamiento desde ferris de cetáceos y aves marinas operativa en todo el archipiélago (www.aviste.me). Los datos recopilados fueron digitalizados y sometidos a diferentes análisis, dependiendo de la procedencia de los mismos. Para conocer la ocurrencia de cada especie de cetáceo en cada estación del año, los datos fueron agrupados y les fue asignada la estacionalidad, atendiendo al mes del año en el que fueron recogidos y teniendo en cuenta que la estación marítima se encuentra desplazada aproximadamente un mes con respecto a la estación terrestre, debido al retardo en el calentamiento y enfriamiento del agua con respecto a la tierra.

La

correspondencia de los meses utilizados en este análisis para definir las estaciones del año se muestra en la tabla 2. Tabla 2. Correspondencia de las estaciones terrestres y marinas. ESTACIÓN MARÍTIMA DEL AÑO

MESES

PRIMAVERA

Abril, Mayo, Junio

VERANO

Julio, Agosto, Septiembre

OTOÑO

Octubre, Noviembre, Diciembre

INVIERNO

Enero, Febrero, Marzo

Fig. 2a Gráfico que muestra la variación anual de la SST en Canarias (1983-2013). Nótese las máximas temperaturas registradas en los meses de septiembre; las mínimas en los marzos.

6

s

s

M

M

s s s

s

s

M

M

M

M

M

s

s s s

M

s

s

M

M

M

M

M

M

s

s

M

M

Fig. 2b Evolución de las temperaturas para la región canaria en los últimos 15 años (2000-2015) según los datos obtenidos del satélite NOAA. Se indican con la inicial los meses de septiembre (más cálidos) y los de marzo (más fríos).

Tabla 3. Tabla resumen que muestra la procedencia de los datos, el periodo en el que fueron recopilados y los años disponibles para el estudio. FUENTE DE DATOS CAM PAÑA CIENTÍFICA

PROYECTO/ COM PAÑÍA

PARÁM ETROS ZIFIOS FANCY II

ESFUERZO: Nº DÍAS O TRAYECTOS / ESTACIÓN

ISLA

AÑOS DISPONIBLES

PRIM AVERA

VERANO

EL HIERRO

2008-2014* *

90

60

71

38

2005-2014* *

478

589

349

267

2004-2014* *

233

276

216

150

1996-2014* * 1997-2013* *

1016 573

1095 918

812 452

779 403

2012-2014* *

118

82

127

89

2009-2010* *

-

4

21

4

LA PALM A PLATAFORM AS DE GRAN SPIRIT OF THE SEA OPORTUNIDAD: CANARIA OBSERVACIÓN ASOCIACIÓN M EER LA GOM ERA COM ERCIAL VARIOS * TENERIFE PLATAFORM AS DE OPORTUNIDAD: CETAVIST TODAS FERRIS DE PASAJEROS CAM PAÑA ZIPHYSETER TODAS CIENTÍFICA

OTOÑO

INVIERNO

* Embarcaciones desde donde se recopilaron los datos en la isla de Tenerife: Kathrine, Lady Shelley, Black Jack, Bonadea II, Montecarlo, Tropical delfín, Aforte, Bahía de Santander, Bahriyeli, Eden Discovery, Free Bird one for you, Free Bird Maxicat, El Teide II, Flipper Uno, Joly Roger, Kosamui, Whale Adventure, Ocean Warrior, Top Class Tour, Giants dolphin, Gladiator, Glass Botton, Muscat I, Muscat II, Nashira Uno, Ocean Exporer, Neptuno, Peter Pam, Paola Rib, Salmon II. Seaquest, Royal Dolphin, Taxi Whale, Shot gun, Sun Caliu. ** No están disponibles para este estudio todas las estaciones de todos los años, debido principalmente a la pérdida del material recopilado a bordo por las empresas. En el caso de la isla de El Hierro, el verano de 2011 no está disponible por encontrarse la isla afectada por los temblores de tierra producidos por la erupción del volcán submarino en la localidad de La Restinga.

Debido a las diferencias entre las metodologías usadas en cada fuente de datos, a continuación se especifica el tratamiento de los mismos según su procedencia.

7

Metodología - Isla de El Hierro Se obtuvieron los datos de las cuatro campañas anuales que se desarrollan en la isla de El Hierro, zonas de Las Playas y Mar de Las Calmas, durante el periodo especificado en la tabla 3, donde los avistamientos desde el barco (estación de mar) fueron dirigidos desde tierra para optimizar el esfuerzo (estación de tierra). Se recopilaron los siguientes parámetros: fecha, hora, especie de mar, especie de tierra, posición geográfica, nº de individuos y

Fig. 3 Mapa del área de estudio (interior de la línea curva), con las dos estaciones, en la isla de El Hierro

comportamiento y fueron agrupados según la estación. Dado que un mismo avistamiento es registrado dos veces (estación de mar y estación de tierra), se eliminaron duplicados escogiendo el nivel de identificación más preciso - desde tierra, existe la limitación de la distancia por lo que muchas observaciones no se identifican a nivel de especie. El elevado número de observaciones hacen que la embarcación solo se aproxime a alguno de los avistamientos-. Los días en los que se realizó esfuerzo sin éxito no fueron contabilizados, dado que el área de estudio se encuentra en zonas de calma con buenas condiciones de mar y viento, por lo que la detectabilidad de los animales es alta. No se realiza esfuerzo cuando las condiciones no son óptimas y el número de días sin avistamientos es despreciable.

Fig. 4a Vista aérea de dos ejemplares de zifio de Blainville, M. densirostris en la Isla de El Hierro. Al fondo la embarcación "Punta Ballena", usada en las campañas "Parámetros Zifios" (Julio, 2014). Foto: "Oceana".

Fig. 4b Estación de Tierra dirigiendo a la embarcación y en esfuerzo de búsqueda de animales. Foto: Equipo Investigacion ULL.

8

Metodología - Islas de La Palma, Gran Canaria, Tenerife y La Gomera Los datos fueron recolectados desde plataformas oportunistas de observación comercial por empresas de avistamiento de cetáceos. En todos los casos, se trata de excursiones que recorren el sector suroeste de las islas, zonas de conocida agregación de muchas especies de cetáceos y de aguas normalmente abrigadas que permiten la actividad. El área de estudio se muestra en la figura 5.

Fig. 5 Mapa que muestra las áreas recorridas por las embarcaciones de observación comercial de cetáceos y los puertos desde donde operan. (Ejebeitia & Urquiola, 2009). No se pudo identificar y por tanto obtener datos de ninguna empresa del sector en Lanzarote y Fuerteventura durante el periodo de estudio.

Los datos proporcionados abarcan los avistamientos registrados durante todos los meses del periodo especificado en la tabla 3. Los casos en los que no se realizó esfuerzo se debieron a malas condiciones meteorológicas que impidieron la navegación. Se estima que las salidas tienen un 90% de probabilidad de avistar al menos una especie de cetáceo, en todas las estaciones, debido a las especies residentes en Canarias. Esto permite utilizar como unidad de esfuerzo el número de días con avistamiento de cualquier especie de cetáceo y referenciar a este esfuerzo la presencia relativa de las distintas especies a lo largo del periodo de estudio. Los datos fueron recogidos por los barcos en fichas de observación en los que se apuntaron los siguientes parámetros: fecha, hora, especie, posición geográfica, rumbo y comportamiento. Fueron transformados a formato digital en este trabajo, y agrupados según la estación del año. Debido a la diferencia en el número de días disponibles de cada año entre las distintas estaciones, los análisis fueron realizados escogiendo únicamente los años donde los días con 9

avistamiento en cada estación superaba un valor mínimo, que asegurara que la especies con poca presencia fueran detectadas y que permitiera usar la mayor cantidad de años posibles. Este valor fue calculado realizando un test de supervivencia (técnica inferencial que tiene como objetivo esencial modelizar el tiempo que se tarda en que ocurra un determinado suceso). (Collect, 2003), aplicado a los datos procedentes de la isla de la Palma y se aplicó el mismo criterio (umbral del número mínimo de días muestreados por estación) al resto de los datos procedentes de barcos de observación comercial. TEST SUPERVIVENCIA - ESTACIONES Nº AÑOS DISPONIBLES

12 10 8

INVIERNO

6

PRIM AVERA

4

OTOÑO

2

VERANO

0 10

20

30

40

50

60

70

80

Nº DÍAS CON AVISTAM IENTO

Fig.6 Test de supervivencia que muestra el número mínimo de días de esfuerzo requeridos en el análisis, n= 20 (cambio de pendiente), para la isla de la Palma.

Los días en las que se realizó esfuerzo sin éxito, no fueron contabilizados puesto que las áreas donde operan estos barcos son siempre zonas de calma, y a que únicamente se navega cuando hay óptimas condiciones meteorológicas y por lo tanto la detectabilidad de los animales es alta.

Metodología utilizada en "Cetavist": red de avistamiento archipelágica Los muestreos se realizan mediante transectos lineales llevados a cabo por voluntarios formados en la recogida científica de datos a bordo de ferris de pasajeros. Éstos, operan en los canales de todas las islas del archipiélago exceptuando el canal Lanzarote - La Graciosa.

Fig.7 Mapa del área de estudio mostrando los transectos lineales realizados por los ferris.

10

La Red de Avistamiento de Cetáceos y Aves Marinas de Canarias (CetAvist) fue creada por la ULL (BIOECOMAC) en 2012 y está financiado (2014 y 2015) por la Fundación Biodiversiad-MAGRAMA. Se incluye en la web "www.aviste.me", albergado en el servidor de la ULL. Un mínimo de 1 observador (situado en posición

alta

y

en

las

mejores

condiciones de luz y viento posibles), o un máximo de 4, (en los trayectos que superan las 4 horas de navegación para mantener un esfuerzo constante) registran todos los cetáceos que se avistan, anotando

los

siguientes

parámetros:

fecha, hora, coordenadas del barco en el momento de la observación, ángulo y distancia al avistamiento, especie, así

Fig.8 Recogida de datos a bordo del ferri "Volcán de Tamadaba" (Naviera Armas) durante un trayecto Cetavist.

como datos sobre el tamaño y comportamiento del grupo avistado. Cada 20 minutos se anotan las condiciones ambientales existentes o siempre que éstas varíen en alguno de sus parámetros. Los datos recopilados durante el periodo especificado en la tabla 3 fueron digitalizados con los mismos campos para los parámetros recogidos a bordo.

Metodología utilizada en la campaña "Ziphyseter" El muestreo recorrió las áreas oceánicas del archipiélago de forma homogénea mediante transectos lineales y se realizó un registro tanto visual como acústico durante el otoño e invierno del periodo especificado en la tabla 3. La metodología detallada desarrollada durante la campaña se especifica en el informe elaborado por el equipo

Fig.9 Mapa del área de estudio mostrando los transectos lineales realizados durante la campaña.

científico que participó en la misma (Aguilar et al. 2010).

11

Los datos proporcionados por esta campaña solo han sido utilizados como fuente de información adicional para el análisis de la distribución invernal de rorcuales en el archipiélago.

Análisis estadísticos descriptivos Para representar la tasa de avistamiento de cada especie a lo largo del tiempo, primero fue necesario ponderar las observaciones de campo por el número de días de esfuerzo. Esto permite eliminar el error que introduce el hecho de la gran diferencia de días de esfuerzo en las distintas estaciones. Se aplicó la siguiente corrección de esfuerzo: (

) =

× 100

nk = nº de días con avistamiento para una especie determinada "k" N= nº de días totales con avistamiento Los barcos de avistamiento de cetáceos no operan en los días con mala mar, por lo que se asume que aunque se realizara esfuerzo, los animales no serían detectados. Por operar en zonas de calma, navegan durante prácticamente todo el año. Para los datos obtenidos en Cetavist, debido a que los barcos navegan siempre bajo casi cualquier condición meteorológica y no exclusivamente en zonas de calma, la corrección del esfuerzo en este caso se hizo ponderando el número de avistamientos de una especie determinada con respecto al número de trayectos en buenas condiciones de mar realizados. Por lo tanto, la tasa de avistamiento, se calculó aplicando la siguiente corrección: (

)=

× 100

nk = nº de avistamientos para una especie determinada "k" NT= nº de trayectos realizados en el que el valor promedio del estado del mar (Coeficiente de Douglas) fue ≤ 4 Para conocer la distribución estacional de las distintas especies, se representó el número de días con avistamiento de cada especie con respecto al número de días con avistamiento totales por estación e isla. Las tendencias de la presencia de las distintas especies a lo largo de los años también fueron representadas.

12

Modelos estadísticos ARIMA- Series Temporales de datos Definimos una serie temporal como un conjunto de “N” observaciones cuantitativas realizadas sobre una variable unidimensional ordenadas en el tiempo. En este caso, la serie temporal está formada por las observaciones de la variable “número de días con avistamiento de la especie x” en cada mes del año. El hecho de que las observaciones estén ordenadas en el tiempo hace que manifiesten una fuerte dependencia o auto correlación que hace que no puedan ser tratados como muestra aleatoria. Cada observación en un momento dado, es modelada en función de los valores anteriores. Los modelos “ARIMA” (Autoregresive Integrated Movieng Average) derivan de sus dos componentes: AR (Auto Regresivo Integrado), MA (Medias Móviles). Estos modelos requieren de series largas de datos. Si se trata de series mensuales, como en este caso, se requiere de entre 6 y 10 años, por lo que no ha sido posible aplicarlo a todas las especies ni a todas las islas. Se usó el software R (www.r-project.org/) con el entorno de programación R studio para llevar a cabo los siguientes análisis. En primer lugar, fue necesaria la representación gráfica de los datos (análisis descriptivo preliminar), de manera que en este punto se observó si la serie era estacionaria en varianza (si no lo es se pueden aplicar transformaciones, por ejemplo la transformación logarítmica) y si existía una tendencia (movimiento a largo plazo que permite observar la evolución general del conjunto de datos), que tendría que ser eliminada, a través de procesos de diferenciación para poder estudiar los patrones estacionales. En el caso de los datos procedentes de las embarcaciones comerciales, existen vacíos probablemente debidos a la pérdida de datos de la embarcación, pero no se puede descartar que sea debido a que no se realizó ningún avistamiento ese día, o bien a que el barco directamente no navegó. Estos valores ausentes pueden ser imputados con métodos de ajuste para asegurar que la serie sea equi-espaciada en el tiempo. Cuanto más bajo sea el nivel en el que se imputan los datos, más fino es el análisis. La función Auto-Arima permite que la serie temporal tenga huecos de datos. La propia serie permite el ajuste completo con un filtro interno tipo Kalman para auto imputar datos a la hora de hacer el modelo, aunque estos datos imputados no los plasma gráficamente. Para comprobar que el modelo tenía un buen ajuste, se aplicaron también otras técnicas que sí requirieron de imputación de datos de forma independiente: - Mediante interpolación de los valores más cercanos por medio de la aplicación de un filtro de Kalman (pondera entre observaciones cercanas y tiene un proceso de retroalimentación en 13

el que no se permiten saltos en los datos estimados que superen las variaciones observadas en los datos previos). - Suavizado exponencial, aplicando en este caso, el método de Holt Winter. Es una técnica de predicción no casual, que permite obtener predicciones basadas en la información histórica contenida en la serie observada. Un análisis visual permite observar si el modelo tiene o no un buen ajuste. Tras aplicar cada metodología, se obtuvieron las funciones de correlación (ACF) y de correlación parcial (PACF) para cada método, lo que permitió obtener las propuestas de los coeficientes del componente (AR), auto regresivo y del estacional (MA), medias móviles del modelo. - Para series sin estacionalidad: ARIMA (p,d,q) - Para series con componente estacional: ARIMA (p,d,q) (P,D,Q) Parte cíclica de la serie temporal (estacionalidad)

Además, se analizaron los residuales que recogen todo aquello que no es explicado por la tendencia o la variación cíclica estacional. Explica el componente aleatorio que pueda existir en el comportamiento de los datos. Para discriminar entre modelos, se consideró el AIC (Akaike Informatio Criterium) descrito por Akaike en 1974 y más recientemente por Burnham y Anderson (2002). El AIC es un valor relativo de la bondad de ajuste de modelos, de forma que el modelo con el valor AIC menor, es el que mejor se ajusta. ARIMA, se aplicó en un análisis comparativo entre dos especies: el delfín mular T. truncatus en la isla de Tenerife cuya presencia en Canarias es permanente durante todo el año, y el delfín común (D. delphis) en la isla de La Gomera, cuya presencia se intuye que presenta estacionalidad. Fueron escogidas las islas de Tenerife y La Gomera por ser de las que mayor cantidad de datos se disponía. La aplicación de estos modelos a las restantes especies e islas no está contemplado en este trabajo pero se realizará una continuación del mismo. Ver ANEXO I.

14

RESULTADOS GENERALES Análisis del esfuerzo realizado El esfuerzo de muestreo se define como el trabajo invertido en la búsqueda de animales y fue calculado para cada isla y cada estación (figuras 10a-10f). DISTRIBUCIÓN ESTACIONAL DEL ESFUERZO (DÍAS) EL HIERRO

14%

DISTRIBUCIÓN ESTACIONAL DEL ESFUERZO (DÍAS) LA PALM A

16%

37%

28% PRIM AVERA

PRIM AVERA

VERANO

VERANO

OTOÑO

OTOÑO INVIERNO

27%

21%

22%

Fig.10a Gráfica que representa el porcentaje del total de días con avistamientos en cada estación del año (I=38; V=60; P=99; O=71 días) en la isla de El Hierro (2008-2014)

DISTRIBUCIÓN ESTACIONAL DEL ESFUERZO (DÍAS) LA GOM ERA

21%

22%

27%

PRIM AVERA

Fig.10b Gráfica que representa el porcentaje del total de días con avistamientos en cada estación del año (I=267; V=589; P=478; O=349 días) en la isla de La Palma (2005-2014)

DISTRIBUCIÓN ESTACIONAL DEL ESFUERZO (DÍAS) GRAN CANARIA

18%

26%

OTOÑO

OTOÑO

30%

INVIERNO

25% 31%

DISTRIBUCIÓN ESTACIONAL DE ESFUERZO TENERIFE

Fig.10d Gráfica que representa el porcentaje del total de días con avistamientos en cada estación del año (I=163; V=277; P=233; O=227 días) en la isla de Gran Canaria (2004-2014)

DISTRIBUCIÓN ESTACIONAL DEL ESFUERZO CETAVIST

18%

29%

PRIM AVERA

PRIM AVERA

VERANO

VERANO OTOÑO

OTOÑO

19%

39%

PRIM AVERA VERANO

INVIERNO

25%

35%

VERANO

Fig.10c Gráfica que representa el porcentaje del total de días con avistamientos en cada estación del año (I=779; V=1095; P=1016; O=812 días) en la isla de La Gomera (1998-2014)

17%

INVIERNO

INVIERNO

Fig.10e Gráfica que representa el porcentaje del total de días con avistamientos en cada estación del año (I=403; V=918; P=573; O=452 días) en la isla de Tenerife (1997-2011)

31%

22%

INVIERNO

Fig.10f Gráfica que representa el porcentaje de trayectos en buenas condiciones de mar (valor promedio del coeficiente Douglas ≤ 4) en cada estación del año (I=3; V=51; P=58; O=62 trayectos) en todo el archipiélago (2012-2014) 15

El esfuerzo de muestreo no es homogéneo en todas las estaciones del año ni en todas las islas. En las plataformas comerciales de observación se intensifica en verano, coincidiendo con la máxima afluencia turística y las buenas condiciones para la navegación. No ocurre lo mismo en la isla del Hierro, en donde a pesar de que parte del estudio se realiza en las zonas de calma, sector suroeste (Mar de Las Calmas), el verano es época de intensas calimas que reducen la visibilidad de manera que se reajusta el reparto de los días de campaña (40/año) hacia estaciones con mejores condiciones como otoño y primavera. El invierno es, para todos los casos, la estación en la que menos esfuerzo se invierte. La distribución estacional del esfuerzo observado en los datos recopilados por Cetavist muestran el otoño como la estación con mayor número de trayectos, lo que obedece no solo a las condiciones climáticas ni marítimas, sino a la disponibilidad de observadores, ya que esta red se surte principalmente de estudiantes cuya tasa de muestreo va en función de su actividad académica. Es esta misma circunstancia la que explica el menor número de trayectos realizados en verano, cuando el curso académico finaliza y desciende el número de observadores disponibles.

Riqueza específica La riqueza específica estacional también fue calculada, entendiéndose como el número de especies observadas en cada estación, corregida por el número de días totales de avistamiento (en cada estación). La evolución de este parámetro a lo largo del tiempo se obtuvo para cada isla y para todos los años, (Fig. 11a-11e). Los datos del eje "y" no son por tanto comparables

Nº ESPECIES CORREGIDO POR UNIDAD DE ESFUERZO

entre sí, solo indican los datos relativos de riqueza estacional para cada isla. RIQUEZA ESPECÍFICA ESTACIONAL EL HIERRO

1 0,8 0,6 0,4 0,2 0

PRIM AVERA VERANO OTOÑO INVIERNO

2008 2009 2010 2011 2012 2013 2014

Fig. 11a Gráfica que representa la riqueza (nº de especies identificadas corregidas por nº de días con avistamiento) en cada estación del año desglosado por año en la isla de El Hierro (I=38; V=60; P=99; O=71).

16

Nº ESPECIES CORREGIDO POR UNIDAD DE ESFUERZO

RIQUEZA ESPECÍFICA ESTACIONAL LA PALM A

0,4 0,3

PRIM AVERA

0,2

VERANO

0,1

OTOÑO

0

INVIERNO

2005 2006 2007 2008 2009 2010 2011 2012 2013 2014

RIQUEZA ESPECÍFICA ESTACIONAL LA GOM ERA

0,4

PRIM AVERA

0,3

VERANO

0,2

OTOÑO

0,1

INVIERNO

2014

2013

2012

2011

2010

2009

2008

2007

2006

2005

2004

2003

2002

2001

2000

1999

1998

1997

0 1996

Nº ESPECIES CORREGIDO POR UNIDAD DE ESFUERZO

Fig.11b Gráfica que representa la riqueza (número de especies identificadas, corregido por nº de días con avistamiento) desglosada por año en la isla de La Palma (I=267; V=589; P=478; O=349).

Nº ESPECIES CORREGIDO POR UNIDAD DE ESFUERZO

Fig. 11c Gráfica que representa la riqueza (número de especies identificadas, corregido por nº de días con avistamiento) desglosada por año en la isla de La Gomera (I=779; V=1095; P=1016; O=812).

RIQUEZA ESPECÍFICA ESTACIONAL GRAN CANARIA 0,4 0,3

PRIM AVERA

0,2

VERANO

0,1

OTOÑO INVIERNO

0

2004

2007

2008

2011

2012

2013

2014

Nº ESPECIES CPRREGIDO POR UNIDAD DE ESFUERZO

Fig. 11d Gráfica que representa la riqueza (número de especies identificadas, corregido por nº de días con avistamiento) desglosado por año en la isla de Gran Canaria (I=163; V=277; P=233; O=227).

RIQUEZA ESPECÍFICA ESTACIONAL CETAVIST

0,8

PRIM AVERA

0,6

VERANO

0,4

OTOÑO

0,2

INVIERNO

0 2012

2013

2014

Fig. 11e Gráfica que representa la riqueza (número de especies identificadas, corregido por nº de trayectos realizados en buenas condiciones de mar (Valor promedio coeficiente Douglas ≤ 4)) en todo el archipiélago (I=31; V=51; P=58;

17

Se han observado un total de 26 especies de las 30 que hay citadas para Canarias (marcadas en negrita en la tabla 1). Tabla 4. Relación del número total de especies, observados en cada isla durante el periodo de estudio. Estos datos no incluyen corrección por número de días de esfuerzo, que es mucho menor en El Hierro que en otras islas. ISLA

Nº ESPECIES ODONTOCETOS

Nº ESPECIES MISTICETOS

TOTAL

EL HIERRO LA PALM A LA GOM ERA GRAN CANARIA TENERIFE

13 17 16 14 2*

1 4 4 5 -

14 21 20 19 2*

* En la isla de Tenerife, solo se recogieron datos para el delfín mular y el calderón de aleta corta.

Tasa de presencia de especies en las islas Para conocer qué especies son las que con más frecuencia se observan en cada isla, se obtuvo su tasa de presencia. Para esto, se calcularon los días en los que estuvo presente cada especie durante todos los años disponibles y se dividió por los días totales de avistamiento en ese mismo periodo (Fig.12a -12e). La distribución no es homogénea y determinadas especies son más frecuentes en unas islas que en otras.

TASA DE AVISTAM IENTO (%)

TASA DE PRESENCIA DE ESPECIES - EL HIERRO 100 80 60 40 20

77,5% 67,1%

ESPECIES QUE HAN SIDO

50,7%

PUNTUALM ENTE AVISTADAS (< 2 DÍAS CON

17,4%

AVISTAM IENTO)

8,9%

7,8%

4,2%

3,3%

1,9%

1,4%

0,94%

0

Fig. 12a Gráfico que representa la tasa de presencia (H= [nk/N] x 100; nk= nº días con avistamiento de la especie "k"; N= nº de días totales con avistamiento) de cada especie en la isla de El Hierro (2008-2014)

18

TASA DE AVISTAM IENTO (%)

TASA DE PRESENCIA DE ESPECIES - LA PALM A 100 80 60 40 20

58,7% ESPECIES QUE HAN SIDO PUNTUALM ENTE AVISTADAS (< 2

28,5% 18%

DÍAS CON AVISTAM IENTO)

12,8%

8,5% 6,1%5,2%

0

3,6% 2,9% 1,7% 0,9% 0,4% 0,4% 0,2% 0,2%

Fig. 12b Gráfico que representa la tasa de presencia (H= [nk /N] x 100; nk= nº días con avistamiento de la especie "k"; N= nº de días totales con avistamiento) de cada especie en la isla de La Palma (2005-2014)

TASA DE AVISTAM IENTO (%)

TASA DE PRESENCIA DE ESPECIES - LA GOM ERA 100 80 60 40 20

ESPECIES QUE HAN SIDO PUNTUALM ENTE AVISTADAS

26,3% 21,7%

(< 2 DÍAS CON AVISTAM IENTO)

16%

9%

4,2% 3,8% 2,2% 1,6% 1,3% 0,9% 0,5% 0,3% 0,2% 0,04% 0,04% 0,04%

0

Fig. 12c Gráfico que representa la tasa de presencia (H= [nk /N] x 100; nk= nº días con avistamiento de la especie "k"; N= nº de días totales con avistamiento) de cada especie en la isla de La Gomera (1998-2014)

TASA DE AVISTAM IENTO (%)

TASA DE PRESENCIA DE ESPECIES - GRAN CANARIA 100 80 60 40 20

41,6% 35,8%

ESPECIES PUNTUALM ENTE AVISTADAS (< 2 DÍAS CON

18,1% 18,1% AVISTAM IENTO) 9,5% 7,3% 6,9% 6,3% 3,3% 0,8% 0,8% 0,8% 0,5% 0,5% 0,3%

0

Fig. 12d Gráfico que representa la tasa de presencia (H= [nk /N] x 100; nk= nº días con avistamiento de la especie "k"; N= nº de días totales con avistamiento) de cada especie en la isla de La Gomera (1998-2014)

19

TASA DE PRESENCIA DE ESPECIES - CETAVIST 100 80

62,5%

60

ESPECIES PUNTUALM ENTE

40

(< 2 AVISTAM IENTOS)

AVISTADAS

15,4%

20

8,7% 7,7% 6,75% 2,6%

1,35% 1,3% 0,6% 0,3% 0,3% 0,3%

0

Fig. 12e Gráfico que representa la tasa de presencia (H= [nk /N] x 100; nk= nº de avistamientos de la especie "k"; N= nº trayectos totales de cada especie en la isla de Gran Canaria (2004 - 2014) en buenas condiciones de mar (Coeficiente de Douglas ≤4)

Las 3 especies más frecuentes son el delfín mular (T. truncatus), el delfín moteado del Atlántico (S. frontalis) y el calderón de aleta corta (G. macrorhynchus). El Hierro constituye una excepción y son los zifios de Cuvier (Z. cavirostris) y el de Blaniville (M. densirostris) las especies más observadas, seguidas del delfín mular. Es interesante observar que el delfín listado (S. coeruleoalba) es común en la isla más oriental estudiada, Gran Canaria y solo ocasionalmente en las más occidentales.

Distribución estacional de la presencia de especies A continuación se muestran los resultados del análisis descriptivo de la presencia estacional de las especies analizadas.

Familia Delphinidae

100 50 0 Fig. 13a

TASA AVISTAM IENTO (%)

Nk=108

PRIM AVERA

VERANO

OTOÑO

T. truncatus - GRAN CANARIA

INVIERNO

50 0 Fig. 13c

VERANO

OTOÑO

T. truncatus - LA GOM ERA

INVIERNO

Nk=988

100 50 0

Fig. 13b

Nk=313

100

PRIM AVERA

TASA AVISTAM IENTO (%)

T. truncatus - EL HIERRO

TASA AVISTAM IENTO (%)

TASA AVISTAM IENTO (%)

El delfín mular - Tursiops truncatus

PRIM AVERA

VERANO

OTOÑO

T. truncatus - LA PALM A

INVIERNO

Nk=1193

100 50 0 Fig. 13d

PRIM AVERA

VERANO

OTOÑO

INVIERNO

20

Nk=1649

100 50 0 PRIM AVERA

Fig. 13e

VERANO

OTOÑO

INVIERNO

TíASA AVSTAM IENTO (%)

TASA AVISTAM IENTO (%)

T. truncatus - TENERIFE

T. truncatus - CETAVIST

Nk=34

100 50 0 PRIM AVERA

VERANO

OTOÑO

INVIERNO

Fig. 13f

Fig. 13a-13d Distribución estacional de T. truncatus. En rojo el valor promedio de la TA en cada estación. En gris, los valores máximos y mínimos registrados. Nk= Nº de días con avistamiento de T. truncatus o Nº avistamientos para los datos extraídos de Cetavist.

Los resultados para el delfín mular confirman el patrón de residencia de esta especie en Canarias con un pico de presencia que varía para cada isla, de forma que en La Palma y La Gomera, la tasa de avistamiento aumenta en verano, mientras que en Tenerife y en El Hierro, es la primavera la que registra mayor porcentaje de días con esta especie. Los datos recabados por Cetavist muestran un pico otoñal.

Tursiops t runcatus, Chloe Yzoard

.

21

Nk=20

100 50 0 PRIM AVERA

VERANO

OTOÑO

INVIERNO

TASA AVISTAM IENTO (%)

Fig. 14a

S. bredanensis - LA PALM A

50 0

Nk=688

100 50 0 PRIM AVERA

VERANO

OTOÑO

INVIERNO

Nk=303

100

Fig. 14b

S. bredanensis - LA GOM ERA

Fig. 14c

TASA AVISTAM IENTO (%)

S. bredanensis - EL HIERRO

PRIM AVERA

VERANO

OTOÑO

S. bredanensis - GRAN CANARIA TASA AVISTAM IENTO (%)

TASA AVISTAM IENTO (%)

Delfín de dientes rugosos - Steno bredanensis

INVIERNO

Nk=60

100 50 0

Fig. 14d

PRIM AVERA

VERANO

OTOÑO

INVIERNO

Fig. 14a-14d Distribución estacional de S. bredanensis. En rojo el valor promedio de la TA en cada estación. En gris, los valores máximos y mínimos registrados. Nk= Nº de días con avistamiento de S. bredanensis o Nº de avistamientos para los datos extraídos de Cetavist.

Aunque con menor tasa de avistamiento que la especie anterior, el delfín de dientes rugosos se encuentra en las aguas del archipiélago durante todo el año, con especial presencia en primavera-verano. Destaca un descenso otoñal, muy evidente en la isla de El Hierro.

St eno bredanensis, Chloe Yzoard

22

50 0

TASA AVISTAM IENTO (%)

Fig. 15a

PRIM AVERA

VERANO

OTOÑO

S. frontalis - LA GOM ERA

INVIERNO

50 0 PRIM AVERA

VERANO

OTOÑO

S. frontalis - CETAVIST

S. frontalis - LA PALM A

INVIERNO

Nk=479

100 50 0

Fig. 15b

Nk=1270

100

Fig. 15c TASA AVISTAM IENTO (%)

Nk=37

TASA AVISTAM IENTO (%)

S. frontalis - EL HIERRO

100

PRIM AVERA

VERANO

OTOÑO

S. frontalis - GRAN CANARIA TASA AVISTAM IENTO (%)

TASA AVISTAM IENTO (%)

Delfín moteado del Atlántico - Stenella frontalis

INVIERNO Nk=364

100 50 0 Fig. 15d

PRIM AVERA

VERANO

OTOÑO

INVIERNO

Nk=49

100 0 Fig. 15e

PRIM AVERA

VERANO

OTOÑO

INVIERNO

Fig. 15a-15e Distribución estacional de S. frontalis. En rojo el valor promedio de la TA en cada estación. En gris, los valores máximos y mínimos registrados. Nk= Nº de días con avistamiento de S. frontalis o Nº de avistamientos para los datos extraídos de Cetavist.

El delfín moteado del Atlántico permanece todo el año en aguas canarias, siendo sus máximos de presencia en primavera e invierno y sus mínimos en otoño. Destaca la elevada tasa de avistamiento registrada en la isla de Gran Canaria, donde el pico primaveral se extiende hasta el verano.

23

Delfín listado -Stenella coeruleoalba

30

0

TASA AVISTAM IENTO (%)

Fig. 16a

PRIM AVERA

VERANO

OTOÑO

S. coeruleoalba - LA GOM ERA

INVIERNO

Nk=166

30 0 PRIM AVERA

VERANO

OTOÑO

INVIERNO

TASA AVISTAM IENTO (%)

Fig. 16c

TASA AVISTAM IENTO (%)

Nk=2

Nk=61

S. coeruleoalba - LA PALM A

30 0 PRIM AVERA

VERANO

OTOÑO

INVIERNO

Fig. 16b TASA AVISTAM IENTO (%)

TASA AVISTAM IENTO (%)

S. coeruleoalba - EL HIERRO

S. coeruleoalba - GRAN CANARIA

Nk=158

30 0 PRIM AVERA

VERANO

OTOÑO

INVIERNO

Fig. 16d

S. coeruleoalba - CETAVIST

Nk=1

30 0 PRIM AVERA

VERANO

OTOÑO

INVIERNO

Fig. 16e

Fig. 16a-16e Distribución estacional de S. coeruleoalba. En rojo el valor promedio de la TA en cada estación. En gris, los valores máximos y mínimos registrados. Nk= Nº de días con avistamiento de S. coeruleoalba o Nº de avistamientos para los datos extraídos de Cetavist.

El delfín listado es observado ocasionalmente en todas las estaciones del año. Su tasa de avistamiento se ve incrementada especialmente en la isla de Gran Canaria, llegando a detectarse hasta en un 30% de los días con un pico presencial en otoño.

St enella coeruleoalba, Chloe Yzoard

24

Nk=4

30

0 PRIM AVERA

VERANO

OTOÑO

INVIERNO

TASA AVISTAM IENTO (%)

Fig. 17a

D. delphis - LA GOM ERA

30

0 Fig. 17c

TASA AVISTAM IENTO (%)

Nk=314

PRIM AVERA

VERANO

OTOÑO

D. delphis - CETAVIST

INVIERNO

TASA AVISTAM IENTO (%)

D. delphis - EL HIERRO

D. delphis - LA PALM A

Nk=143

30

0

Fig. 17b TASA AVISTAM IENTO (%)

TASA AVISTAM IENTO (%)

Delfín común - Delphinus delphis

PRIM AVERA

VERANO

OTOÑO

D. delphis - GRAN CANARIA

INVIERNO

Nk=55

30 0 PRIM AVERA

VERANO

OTOÑO

INVIERNO

Fig. 17d

Nk=25

30 0 PRIM AVERA

VERANO

OTOÑO

INVIERNO

Fig. 17e

Fig. 17a-17e Distribución estacional de D. delphis. En rojo el valor promedio de la TA en cada estación. En gris, los valores máximos y mínimos registrados. Nk= Nº de días con avistamiento de D. delphis o Nº de avistamientos para los datos extraídos de Cetavist.

El delfín común muestra una marcada estacionalidad visitando las aguas insulares en primavera e invierno tal y como muestran los datos obtenidos para todas las islas, abandonando la región en verano y otoño. La Palma y La Gomera parecen ser zonas preferentes para esta especie que alcanza tasas de avistamiento superiores al 20% en esas estaciones.

Delphinus delphis, Chloe Yzoard

25

Nk=1

0 PRIM AVERA

VERANO

OTOÑO

INVIERNO

Fig. 18a TASA AVISTAM IENTO (%)

G. griseus - LA PALM A

Nk=7

10 0 PRIM AVERA

VERANO

OTOÑO

INVIERNO

Fig. 18b

G. griseus - LA GOM ERA

Nk=17

10

0 PRIM AVERA

VERANO

OTOÑO

INVIERNO

Fig. 18c TASA AVISTAM IENTO (%)

TASA AVISTAM IENTO (%9

G. griseus - EL HIERRO

20

TASA AVISTAM IENTO (%)

TASA AVISTAM IENTO (%)

Calderón gris - Grampus griseus

G. griseus - GRAN CANARIA

Nk=83

50 0 PRIM AVERA

VERANO

OTOÑO

INVIERNO

Fig. 18d

G. griseus - CETAVIST

Nk=4

10

0 Fig. 18e

PRIM AVERA

VERANO

OTOÑO

INVIERNO

Fig. 18a-18e Distribución estacional de G. griseus. En rojo el valor promedio de la TA en cada estación. En gris, los valores máximos y mínimos registrados. Nk= Nº de días con avistamiento de G. griseus. o Nº de avistamientos para los datos extraídos de Cetavist.

El calderón gris muestra una distribución temporal desigual según el área y a pesar de estar presente en todas las estaciones del año, sus tasas de avistamiento están por debajo del 10% en todas las islas. La excepción de nuevo es Gran Canaria, donde esta especie llega a ser frecuente, intensificando su presencia en invierno.

Grampus griseus, Chloe Yzoard

26

Calderones - Globicephala sp. Debido a la gran dificultad que supone distinguir a las dos especies de este género desde las embarcaciones (G. melas, calderón de aleta larga o común y G. macrorhynchus, calderón de aleta corta o tropical), se han analizado conjuntamente tratándolo a nivel de género, aunque

50 0

TASA AVISTAM IENTO (%)

PRIM AVERA

VERANO

OTOÑO

Globicephala sp. - LA GOM ERA

INVIERNO

Nk=1722

100 50 0

Fig. 19c

PRIM AVERA

VERANO

OTOÑO

Globicephala sp. - TENERIFE

INVIERNO

50 0 PRIM AVERA

VERANO

OTOÑO

Globicephala sp. - LA PALM A

INVIERNO

Nk=215

100 50 0 PRIM AVERA

Fig. 19b

VERANO

OTOÑO

Globicephala sp. - GRAN CANARIA

INVIERNO

Nk=64

100 50 0

Fig. 19d

Nk=2075

100

Fig. 19e

TASA AVISTAM IENTO (%)

100

Fig. 19a

TASA DE AVISTAM IENTO (%)

Nk=17

TASA AVISTAM IENTO (%)

Globicephala sp. - EL HIERRO

TASA AVISTAM IENTO (%)

TASA AVISTAM IENTO (%)

existe una gran probabilidad de que se trate en todos los casos del calderón de aleta corta.

PRIM AVERA

VERANO

OTOÑO

Globicephala sp. - CETAVIST

INVIERNO

Nk=217

150 0 PRIM AVERA

VERANO

OTOÑO

INVIERNO

Fig. 19f

Fig. 19a-19f Distribución estacional de Globicephala sp. En rojo el valor promedio de la TA en cada estación. En gris, los valores máximos y mínimos registrados. Nk= Nº de días con avistamiento de Globicephala sp. o Nº de avistamientos para los datos extraídos de Cetavist.

El calderón se encuentra en Canarias durante todo el año. Los datos obtenidos para la especie dejan ver un pico otoñal para todos los casos de estudio y unas tasas de avistamiento excepcionalmente altas en las islas de Tenerife y La Gomera. Cetavist también refleja una elevada tasa de encuentro que hace que sea la especie más frecuente registrada en esta red.

Globicephala macrorhynchus, Chloe Yzoard

27

0 PRIM AVERA

VERANO

OTOÑO

O. orca - LA GOM ERA

INVIERNO

VERANO

OTOÑO

Fig. 20c

0 PRIM AVERA

Nk=3

0 PRIM AVERA

3

Fig. 20b

3

INVIERNO

Nk=1

O.orca - LA PALM A

TASA AVISTAM IENTO (%)

30

Fig. 20a TASA AVISTAM IENTO (%)

Nk=3

O. orca - EL HIERRO

VERANO

OTOÑO

O. orca - CETAVIST TASA AVISTAM IENTO (%)

TASA AVISTAM IENTO (%)

Orca - Orcinus orca

INVIERNO

Nk=1

3 0 PRIM AVERA

VERANO

OTOÑO

INVIERNO

Fig. 20d

Fig. 20a-20d Distribución estacional de O. orca. En rojo el valor promedio de la TA en cada estación. En gris, los valores máximos y mínimos registrados. Nk= Nº de días con avistamiento de O. orca. o Nº de avist am ient os para los dat os ext raídos de Cet avist .

La orca ha sido observada de manera esporádica en las islas, sin embargo es de destacar que el 60% de los días con avistamientos suceden en verano.

Orcinus orca, Chloe Yzoard

28

Familia Ziphiidae Zifio de Blainville - Mesoplodon densirostris

100 0

TASA AVISTAM IENTO (%)

Fig. 21a

PRIM AVERA

VERANO

OTOÑO

M . densirostris - LA GOM ERA

0 PRIM AVERA

VERANO

OTOÑO

INVIERNO

Nk=88

M . densirostris - LA PALM A

30 0 Fig. 21b

Nk=87

30

Fig. 21c TASA AVISTAM IENTO (%)

INVIERNO

TASA AVISTAM IENTO (%)

Nk=143

TASA AVISTAM IENTO (%)

TASA AVISTAM IENTO (%)

M . densirostris - EL HIERRO

PRIM AVERA

VERANO

OTOÑO

INVIERNO

M . densirostris - GRAN CANARIA

Nk=7

10

0

Fig. 21d

PRIM AVERA

VERANO

OTOÑO

INVIERNO

Nk=3

M . densirostris - CETAVIST

10

0

Fig. 21e

PRIM AVERA

VERANO

OTOÑO

INVIERNO

Fig. 21a-21e Distribución estacional de M. densirostris. En rojo el valor promedio de la TA en cada estación. En gris, los valores máximos y mínimos registrados. Nk= Nº de días con avistamiento de M. densirostris o Nº de avist amient os para los dat os extraídos de Cet avist.

El zifio de Blainville se encuentra presente todo el año en Canarias como demuestran los resultados obtenidos. Es especialmente frecuente en la isla de El Hierro y también es observado ocasionalmente en La Palma y La Gomera. No se observan diferencias reseñables en su distribución estacional.

M esoplodon densirost ris, Chloe Yzoard

29

Nk=165

100 0 PRIM AVERA

VERANO

OTOÑO

INVIERNO

TASA AVISTAM IENTO (%) TASA AVISTAM IENTO (%)

Z. cavirostris - LA PALM A

Nk=28

30

0 PRIM AVERA

VERANO

OTOÑO

INVIERNO

Fig. 22b

Fig. 22a

Z. cavirostris - LA GOM ERA

Nk=77

30

0 Fig. 22c

TASA AVISTAM IENTO (%)

Z.cavirostris - EL HIERRO

PRIM AVERA

VERANO

OTOÑO

Z. cavirostris - CETAVIST

INVIERNO

TASA AVISTAM IENTO (%)

TASA AVISTAM IENTO (%)

Zifio de Cuvier - Ziphius cavirostris

Z. cavisrostris - GRAN CANARIA

Nk=4

10 0 PRIM AVERA

VERANO

OTOÑO

INVIERNO

Fig. 22d

Nk=7

10 0 PRIM AVERA

VERANO

OTOÑO

INVIERNO

Fig. 22e

Fig. 22a-22e Distribución estacional de Z. cavirostris. En rojo el valor promedio de la TA en cada estación. En gris, los valores máximos y mínimos registrados. Nk= Nº de días con avistamiento de Z. cavirostris o Nº de avist amient os para los dat os extraídos de Cet avist.

El zifio de Cuvier es observado con más frecuencia en el otoño para casi todos los casos, aunque los pocos avistamientos registrados en la isla de Gran Canaria ocurren en verano. El Hierro es de nuevo la isla en la que más zifios de esta especie se detectan. Se puede observar en las islas todo el año.

Ziphius cavirost ris

30

Familia Physeteridae

30

0

TASA DE AVISTAM IENTO (%)

Fig. 23a

PRIM AVERA

VERANO

OTOÑO

P. macrocephalus- LA GOM ERA

INVIERNO

0 PRIM AVERA

VERANO

OTOÑO

P. macrocephalus - CETAVIST

INVIERNO

Nk=48

P. macrocephalus - LA PALM A

30

0

Fig. 23b

Nk=39

10

Fig. 23c TASA AVISTAM IENTO (%)

Nk=9

TASA AVISTAM IENTO (%)

P. macrocephalus - EL HIERRO

TASA AVISTAM IENTO (%)

TASA AVISTAM IENTO (%)

Cachalote - Physeter macrocephalus

PRIM AVERA

VERANO

OTOÑO

INVIERNO

P. macrocephalus - GRAN CANARIA

Nk=7

10 0

Fig. 23d

PRIM AVERA

VERANO

OTOÑO

INVIERNO

Nk=24

30 0 Fig. 23f

PRIM AVERA

VERANO

OTOÑO

INVIERNO

Fig. 23a-23f Distribución estacional de P. macrocephalus. En rojo el valor promedio de la TA en cada estación. En gris, los valores máximos y mínimos registrados. Nk= Nº de días con avistamiento de P. macroceus o Nº de avist amient os para los dat os extraídos de Cet avist.

Los escasos datos disponibles para el cachalote muestran como resultado una gran disparidad en la distribución temporal de esta especie en las distintas islas, observándose en algunos casos una ausencia total de avistamientos en otoño como en el Hierro y en los datos aportados por Cetavist y a la vez picos de presencia en esa misma estación en otras zonas como La Palma o Gran Canaria. Los resultados demuestran que el cachalote está presente en Canarias durante todo el año.

Physet er macrocephalus, Chloe Yzoard

31

Familia Balaenopteridae La familia Balaenopteridae engloba a las 6 especies de rorcuales existentes: el rorcual azul (Balaenoptera musculus), el rorcual común (Balaenoptera physalus), el rorcual boreal o norteño (Balaenoptera borealis), el rorcual tropical (Balaenoptera edeni brydei), el rorcual aliblanco (Balaenoptera acutorostrata) y la yubarta (Megaptera novaengliae). La dificultad para identificar a nivel de especie a los rorcuales requiere que los datos sean tratados a nivel de

30 0 PRIM AVERA

VERANO

OTOÑO

INVIERNO

TASA AVISTAM IENTO (%)

FAM . BALAENOPTERIDAE - LA GOM ERA

TASA AVISTAM IENTO (%)

Nk=258

70 0 PRIM AVERA

VERANO

OTOÑO

FAM .BALAENOPTERIDAE - CETAVIST

INVIERNO

Nk=161

FAM . BALAENOPTERIDAE - LA PALM A

70

0

Fig. 24b

Fig. 24a

Fig. 24c

TASA AVISTAM IENTO (%)

Nk=9

FAM . BALAENOPTERIDAE - EL HIERRO

TASA AVISTAM IENTO (%)

TASA AVISTAM IENTO (%)

familia, aunque el grado de aportación de cada una de ellas sea diferente.

PRIM AVERA

VERANO

OTOÑO

INVIERNO

FAM .BALAENOPTERIDAE - GRAN CANARIA

Nk=411

100 0

Fig. 24d

PRIM AVERA

VERANO

OTOÑO

INVIERNO

Nk=44

30 0 PRIM AVERA

VERANO

OTOÑO

INVIERNO

Fig. 24e

Fig. 24a-24b Distribución estacional de la Familia Balaenopteridae. En rojo el valor promedio de la TA en cada estación. En gris, los valores máximos y mínimos registrados. Nk= Nº de días con avistamiento de la Fam. Balaenopteridae o Nº de avistamientos para los datos extraídos de Cetavist.

La distribución temporal de los rorcuales en Canarias abarca todo el año aunque su presencia disminuye en invierno. Pese a las diferencias que se observan en cada caso, la primavera y el verano son las estaciones en la que están más presentes. Es reseñable la elevada frecuencia de avistamiento de estas especies en la isla de Gran Canaria durante el verano, donde se dan tasas máximas próximas al 80%.

Balaenopt era edeni brydei, Chloe Yzoard

32

Análisis comparativo-modelos estadísticos ARIMA En este trabajo se presenta el ejemplo de este análisis en dos especies de estudio, una con carácter residente en Canarias y otra con presencia estacional. A continuación se muestran las gráficas en bruto para cada una de las especies en los años de estudio utilizados para este análisis, T. truncatus (2002-2012), D.delphis (2000-2013):

Fig. 25a Gráfico que representa la variación de la tasa de avistamiento de T. truncatus (20022012). Las líneas rojas señalan los valores donde se concentran la mayoría de los datos, observándose que en la mayor parte, las TA están entre el 50-90%.

Fig 25b. Gráfico que representa la variación de la tasa de avistamiento de D. delphis (20002013). Las líneas rojas señalan los valores donde se concentran la mayoría de los datos en los meses en los que está presente la especie: TA= 3-30%.

El delfín mular presenta tasas de avistamiento más altas que el delfín común, y las grandes oscilaciones son debidas a que los datos han sido recopilados desde muchas embarcaciones que operan en diferentes puertos de Tenerife. Además, esta especie efectúa movimientos 33

entre islas (Tobeña et al. 2014)y la caída de la TA podría atribuirse a los desplazamientos que realiza esta especie en busca de recursos tróficos. Por su parte, las oscilaciones del delfín común parecen obedecer a un patrón estacional, con máximos en la TA que se repiten en los meses de abril. No hay vacíos en los datos para esta especie analizados aquí, lo que significa que TA=0 representa una verdadera ausencia. Nótese que ésta sucede en las estaciones de verano y otoño. - Primer método de ajuste, Filtro de Kalman: imputación por valores cercanos

Fig. 26a Función de auto correlación y función de auto correlación parcial de T. truncatus usando un Filtro de Kalman. Órdenes de la componente: AR (2,0,0).

Fig. 26b Función de auto correlación y función de auto correlación parcial de D. delphis usando un Filtro de Kalman. Órdenes de la componente: AR (3,0,0) (2,0,0)

Las funciones de auto correlación ACF y las funciones de auto correlación parcial PACF indican que el delfín mular es una especie en la que no se detecta componente estacional. La máxima correlación se da cuando se compara el mes 1 con él mismo, y a partir de ahí comienza un decrecimiento exponencial. La PCAF muestra una sola correlación significativa, síntoma de que la especie no presenta estacionalidad. Por el contrario, el 34

delfín común muestra una ACF con la primera correlación significativa, y otra que se repite a los 12 meses (indica un patrón de repetición anual). Ocurre de forma similar para la PACF lo que indica la presencia de un componente estacional. Las gráficas a continuación muestran los datos observados en negro y las predicciones de los modelos con una línea punteada en rojo, así como una predicción del modelo a partir del fin del periodo de estudio.

Fig.27a Propuesta de modelo usando un Filtro de Kalman para T. truncatus. La línea azul representa la predicción y el sombreado azul los intervalos de confianza. Los coeficientes propuestos por este método = ARIMA (2,0,0)

Fig. 27b Propuesta de modelo usando un Filtro de Kalman para D. delphis. La línea azul representa la predicción y el sombreado azul los intervalos de confianza. Los coeficientes propuestos por este método = ARIMA (3,0,0)(2,0,0). No deben considerarse valores