1
Valentín Pérez-Mellado
Departamento de Biología Animal, Ecología
Edafología, Química Agrícola y Parasitología
Facultad de Biología
Edificio de Farmacia, 5ª planta, Campus Miguel de Unamuno
Universidad de Salamanca, Salamanca 37071. España (Spain)
e-mail: [email protected]
INFORME FINAL
Programa de seguimiento de las poblaciones de la lagartija balear,
Podarcis lilfordi (Squamata, Lacertidae) en Menorca
Introducción
En este proyecto se lleva a cabo el seguimiento de todas y cada una de las
poblaciones de la lagartija balear, Podarcis lilfordi existentes en la actualidad en los
islotes costeros de Menorca. Las citadas poblaciones incluyen los islotes de: Aire,
Binicodrell gros, Binicodrell petit, Bledas, Sanitja (o Porros), Sargantana, Rovells,
Porros (de Fornells), s'en Tosqueta, Addaia gran, Addaia petita, ses Aguiles, en Carbó,
en Carbonet, ses Mones, Colom, Mel y Rei.
2
Nuestra evaluación se basa en varios grupos de informaciones: las densidades y
efectivos poblaciones estimados a partir de las mismas, la razón de sexos y la estructura
demográfica básica en edades y sexos, así como la identificación de potenciales
amenazas: presión humana, usos de los islotes, presencia de especies competidoras y
potenciales depredadores. Este proyecto se articula como un Programa de
Seguimiento, a desarrollar anualmente en las citadas poblaciones y que permitirá
identificar con rapidez y eficacia cualquier amenaza o perturbación de las mismas. Se
propone más abajo un protocolo preciso para cada una de las poblaciones conocidas,
estrictamente diferenciado para cada una de ellas y con un primer estudio más
pormenorizado en esta primera temporada, en el cual se establece el protocolo definitivo
para el resto de seguimientos anuales.
Objetivos
Se estiman los siguientes parámetros:
- Densidad de población por medio de las metodologías adecuadas para cada uno de los
enclaves (ver más adelante), en todos los islotes. Se expresa dicha densidad en
individuos/hectárea, con un descriptor de tendencia central (media aritmética) y los
correspondientes descriptores de variabilidad, incluyendo los límites al 95%.
- Efectivos de todas y cada una de las poblaciones insulares de Podarcis lilfordi.
- Análisis comparativo de la densidad de población de cada enclave en relación a los
factores ecológicos potencialmente relevantes como la superficie del islote, su grado de
aislamiento, la disponibilidad de recursos tróficos, la cobertura vegetal, la presencia de
potenciales competidores y la presencia de potenciales depredadores.
- Variaciones sufridas por las poblaciones de las que se disponga de información
anterior sobre sus efectivos poblacionales.
- Estructura demográfica básica de cada una de las poblaciones, incluyendo en la
misma: la proporción de sexos adultos o razón de sexos terciaria, esto es, de individuos
adultos reproductivos, así como la proporción de individuos juveniles o subadultos en
3
cada población.
Equipo investigador
Equipo investigador de la Universidad de Salamanca dirigido por el Dr.
Valentín Pérez Mellado (DNI: 2843962N), Catedrático de Zoología de la Universidad
de Salamanca, que actuará también como investigador principal, junto con tres
colaboradores, la Dra. Ana Pérez Cembranos (DNI: 70870827-F) y los licenciados en
Biología, Alicia León Lobera (DNI: 71948951Y), Becaria FPU de la Universidad de
Salamanca y Gonzalo Manuel Rodríguez Ruíz (DNI: 70818608K), alumno del
Programa de Doctorado de Biología y Conservación de la Biodiversidad de la
Universidad de Salamanca, de modo que en las salidas de campo el trabajo será
efectuado por no más de tres o cuatro personas. Este reducido número de investigadores
garantiza la ausencia de perturbaciones en las visitas a los islotes. En años sucesivos, la
participación de dos investigadores sería suficiente para cubrir los objetivos del
seguimiento.
Área de estudio
El presente trabajo abarca la totalidad de los islotes con poblaciones de lagartija
balear en Menorca (ver listado en la Introducción).
Metodología
Visitas a los enclaves insulares
Se llevaron a cabo entre el mes de abril y el mes de agosto de 2015,
coincidiendo con el período de máxima actividad anual y en el cual se pueden efectuar
estimaciones de densidad fiables y comparables. En años sucesivos, el seguimiento se
llevaría a cabo dentro de estos mismos períodos de tiempo.
Los islotes han sido visitados por un máximo de tres investigadores, en muchos
casos sólo por dos o uno solo, especialmente en el caso de los transectos lineales. Cada
islote ha sido visitado al menos en una ocasión, preferentemente durante el período
4
diario de máxima actividad que, en primavera, va de las 08.00 a las 12.00 hora solar y
en verano comienza y termina una o dos horas después. Los muestreos se realizaron
preferentemente durante dicho período. En todos los casos, los transectos lineales han
sido realizados por el mismo observador (el director del proyecto), con objeto de
minimizar las diferencias en las estimaciones de dos o más observadores y permitir las
comparaciones entre las estimaciones actuales y las llevadas a cabo en años anteriores.
Se procedió a la detección de potenciales depredadores y presencia de
potenciales competidores. Toda esta información ha sido complementada con la ya
disponible en la literatura y en los datos inéditos del director del proyecto (ver más
abajo).
Trabajo de campo
Estimación de los efectivos y densidad de población
El objetivo principal de este proyecto es obtener estimaciones de densidad y
efectivos poblacionales de varios enclaves dentro de un período de tiempo limitado, de
modo que las mismas tengan un valor comparativo entre poblaciones y entre años. En
los enclaves de mayor superficie y densidad media o alta se empleó el transecto lineal
(ver más adelante).
La movilidad topográfica y la vegetación presente en los islotes determinan en
gran medida la linealidad del transecto, pero de todos modos, el método analítico
empleado presupone que los transectos son una sucesión de líneas rectas de variable
longitud que se articulan unas con otras de modo que se eviten la pseudorréplicas, esto
es, la duplicación en la observación de un mismo individuo. Otro aspecto a tener en
cuenta es la dificultad de establecer con precisión la distancia contacto-eje dada la
rápida huida de Podarcis lilfordi una vez realizado el avistamiento. Este hecho depende
también de la población estudiada, ya que el comportamiento ante la presencia del
hombre se ha mostrado diferente en función del islote.
En los islotes, como se ha mencionado más arriba, se emplea como
método de estimación de la densidad el transecto lineal, siguiendo para ello las premisas
básicas de Burnham et al. (1980). Para ello, usamos como estimador estadístico de
densidad el algoritmo basado en las series de Fourier, así como otros algoritmos
5
disponibles en el programa Rdistance, utilizando las distancias perpendiculares a los
individuos como dato base y una anchura de banda del censo fijada a posteriori gracias
a la distancia del individuo más alejado a dicha línea. El transecto lineal permite
estimaciones fiables de la densidad si se cumplen las premisas básicas expuestas por
Burnham et al. (1980, ver una actualización del método en Thompson et al., 1998 y
referencias incluidas), siendo las principales que el número de contactos en el censo sea
superior a los 30, que la probabilidad de detección de los individuos en la línea de
transecto sea 1 y que las distancias perpendiculares de los individuos a la línea de
transecto sean correctamente medidas. En islotes de elevada densidad y orografía
relativamente llana este método rinde resultados excelentes.
Metodología analítica
Existen numerosos programas para el análisis de transectos lineales, pero cada
uno de ellos tiene sus propias particularidades, ventajas y desventajas. Hemos hecho
análisis con varios programas como "Distance" (Miller, 2014) y "unmarked" (Fiske &
Chandler, 2011), así como "Rdistance", todos ellos de libre acceso en el entorno R (R
Developmental Core Team, 2009). El programa que ha rendido unos resultados más
satisfactorios en cuanto al ajuste de los modelos empleados y a la coherencia de las
estimaciones ha sido "Rdistance" (McDonald, 2012).
Todos los datos obtenidos se incluyen en hojas de cálculo de Excel de las cuales
se generan después archivos .csv, esto es, un archivo delimitado por comas, que es el
legible por los paquetes estadísticos del entorno R. Una vez leído el archivo en R y
asignado a un objeto "data.frame", se carga el correspondiente programa para llevar a
cabo los análisis sobre las selecciones de cada transecto. En el caso de programas cortos
de captura-recaptura, se llevará a cabo el siguiente protocolo: Se procederá a tomar
sobre cada individuo capturado los datos de: sexo, edad, peso (por medio de
dinamómetro o bien balanza digital portátil), longitud cabeza-cuerpo (LCC), longitud
del píleo, anchura de la cabeza, altura de la cabeza, anchura de la comisura bucal,
longitud de la pata posterior izquierda, longitud de la cola y estado de regeneración de
la misma. Todos estos datos se obtuvieron en campo. La identificación de los
individuos se llevó a cabo por medio de fotografías digitales tomadas con objetivo
macro y analizadas posteriormente con un software de análisis de imágenes apropiado
como Image Tool ® o Wild-ID. La totalidad de los individuos capturados han sido
6
inmediatamente liberados de nuevo en cada población.
Parra realizar el análisis de densidad en la sintaxis propia de
R se hace usos de la siguiente instrucción general:
“densidad<-F.automated.CDA(datos,area=1000*10,total.trans.len=385,w.hi=2,likelihoods=c("h
alfnorm","uniform","gamma"),series=c("cosine","simple","hermite"),expansions=c(0,1,2),plot=
TRUE)”.
Esta instrucción general un objeto al que llamamos “densidad” y al que se le asigna el
análisis
que realiza la función “F.automated.CDA”, que es un análisis de distancias clásico y
“automatizado”. Esto significa que al contrario de otros métodos, en “Rdistance” es posible
hacer todo el análisis con una sola función, mientras que otros métodos precisan primero de la
creación de un objeto con la función de detectabilidad y después, a partir de dicha función,
estimar la densidad.
Como se puede observar, en la función “F.automated.CDA” se incluyen entre paréntesis
distintos elementos. A continuación se describe cada uno de ellos:
-- “datos”: vector numérico extraído previamente del “data.frame”, que contiene el
listado de distancias perpendiculares a la línea de transecto (en metros).
“area=1000*10”: se da el área o superficie de referencia. Como la unidad de trabajo es
el “metro”, al indicar un área de 1000x10=10000, la densidad quedará expresada como
individuos/hectárea.
“total.trans.len.=385”: indica la longitud total del transecto. En este caso, la longitud
del transecto ascendería a 385 metros.
“wi.hi=2”: con este elementos se señala al modelo que la anchura de banda máxima
será de 2 metros. De este modo, excluimos de la estimación cualquier individuo que se
haya localizado a más de 2 metros de distancia. La razón es que tales individuos
constituyen en casi todos los transectos llevados a cabo un porcentaje minoritario, no
superior al 10% de las observaciones y a tal distancia, la experiencia indica que la
mayoría de los individuos no son detectados, de modo que los observados, por
ejemplo, a 2.5 metros, deben considerarse puntos externos (“outliers”) que afectan de
modo muy significativo a las estimaciones y, por ello, se recomienda su exclusión
(Buckland et al., 2001).
7
“likelihooods=c(“halfnorm”,”uniform”,”gamma”,”negexp”,”hazrate”)”: aquí, se
incluye en el modelo los parámetros de ajuste. El caso de “Rdistance” pueden ser cinco
diferentes: “halfnorm”, ”uniform”, ”gamma”, ”negexp” y ”hazrate”. El más empleado (y
que suele ser también el que mejor se ajusta en la mayoría de los casos) es el
“halfnormal”. De hecho, es el modelo único en muchos paquetes estadísticos de
transectos lineales. Sin embargo, en algunos casos no es un modelo que se ajuste
correctamente a las observaciones disponibles. La estrategia seguida en este estudio ha
sido probar todos los modelos disponibles y escoger el de mejor ajuste a los datos
observados.
“series=c(“cosine”,”simple”,”hermite”)”: el elemento “series” incluye los vectores de
los tipos de series de expansión, que son términos matemáticos asociados a cada
modelo para mejorar su ajuste. Cada una de las series puede combinarse con cada uno
de los modelos. De modo que la combinación del modelo de máxima verosimilitud
halfnorm” y la serie “coseno” es una de las más empleadas, es la conocida como serie
de Fourier, que se ha empleado como “modelo para todo” en numerosos estudios
(Burnham et al., 1980; Buckland et al.,2001).
“expansions=c(0,1,2)”: las denominadas más arriba series de expansión lo pueden ser
empleando ninguno, uno o dos parámetros (de ahí el intervalo, ver más detalles en
Mc Donald, 2011).
“plot=TRUE”: se trata de una de las instrucciones lógicas del programa que le indica al
análisis que se construyan, para cada modelo ajustado, las gráficas con las
probabilidades de detección y las frecuencias de observaciones en función de la
distancia a la línea de transecto. En los resultados se incluyen dichas gráficas para los
modelos mejor ajustados y que finalmente se han empleado para estimar las
densidades.
Los resultados obtenidos con esta instrucción se explican en el primer transecto,
correspondiente a la isla del Aire, del apartado de Resultados. Los resultados del resto de
transectos son todos ellos similares.
8
Resultados
Aire. Realización anual de un transecto general en la zona central de la isla, de
aproximadamente 450 metros de longitud, desde la formación de tamarindos en la base
de la península septentrional, hasta la entrada del recinto del faro. Este transecto
discurre por la zona central de la isla cubierta por Suaeda vera, Carlina corymbosa y
otras plantas vasculares y representa la mayor zona de la isla ocupada por la lagartija
balear. El transecto se ha llevado a cabo en época de máxima actividad de forma
intermitente desde 1989, de modo que su realización anual permite establecer una
comparación plurianual fiable sobre una de las poblaciones más importantes de
Menorca. Además, en este primer año se propuso y se han llevado a cabo transectos
lineales de menor longitud (unos 100 metros cada uno) en toda la superficie central de
la isla, a intervalos de unos 50 metros, desde el recinto del faro hasta el muelle, a ambos
lados del camino. Se completan dichos transectos con cuatro transectos lineales que
abarcan la totalidad de la superficie de la península occidental, así como dos transectos
en la zona más oriental de la isla y otros dos en la península norte. Se posee información
similar de tres años anteriores, lo que permite establecer con elevada fiabilidad la
densidad media de lagartijas en toda la isla, así como los efectivos poblacionales.
Aquí presentamos únicamente los datos correspondientes a la superficie central
de la isla, que es la comparable con los años anteriores.
Coefficients: Sigma 0.7619686 (coeficiente del modelo empleado).
Function: HALFNORM (indica que se ha empleado el algoritmo Halfnorm para
estimar la densidad).
Strip: 0 to 2 (indica que la anchura de banda fue de 0 a 2 metros a cada lado de la línea
de transect).
Effective strip width: 0.9467049 (anchura de banda efectiva del modelo, es la ESW de
las tablas).
Scaling: g(0) = 1 (Indica que el modelo se ajusta a la premisa de una probabilidad igual a
1 de detección de los individuos en la línea de transecto (esta premisa es común a casi
todos los modelos disponibles)). Log likelihood:53.90918 AIC: 55.94019 (“Akaike Information Criterion”, o criterio de
información de Akaike: es un parámetro esencial que indica el ajuste del modelo. Cuanto menor
es AIC, mejor es el modelo. De modo que de todos los modelos posibles, se escoge el que tiene
9
menor
AIC).
Abundance estimate: 2730.361 individuos/hectárea; 95% CI=( 2364.575 to 3155.457 ).
Estimación final de la abundancia que, equivale a la densidad de individuos por hectárea.
A continuación se muestra la función de detección y frecuencias de distancias perpendiculares
para el transecto “Aire” correspondiente al recorrido habitual en la isla del Aire:
0.0
0.5
1.0
1.5
0 0.5 1 1.5 2
Distance
Pro
babi
lity
of d
etec
tion
halfnorm, 0 expansions
ESW = 0.947
BEST FITTING FUNCTION
Figura 1.-Función de detección y frecuencias de distancias perpendiculares en la isla del Aire (para esta y el resto de curvas de detectabilidad, la curva aparece dibujada en rojo y superpuesta al diagrama de barras de las frecuencias de probabilidad de detección a cada distancia perpendicular a la línea de transecto, en abcisas). Nótese que la curva de
detectabilidad comienza siempre en 1.
10
Tabla1. Densidades registradas en la isla del Aire a lo largo de los años (Contactos= número de individuos detectados en los transectos, Longitud= longitud en metros de cada transecto), Densidad= estimación de la densidad en número de individuos por hectárea, SE= error estándar de la estimación,
AIC= función de Akaike de ajuste del modelo, ESW= anchura efectiva de banda, D LCL= límite inferior de confianza al 95%, D UCL= límite superior de confianza al 95%, D CV= coeficiente de variación en
tanto por uno) Año Mes Contactos Longitud Densidad SE AIC ESW D LCL D UCL D CV 1982 4 43 70 766,27 148,29 552,39 400,83 520,39 1128,3 0,1935 1987 8 42 544 541,53 135,98 427,69 71,29 328,4 892,97 0,2511 1988 8 28 273 337,59 112,84 310,25 151,91 172,93 659,04 0,3343 1989 6 157 273 1666 168,2 1817,8 172,6 1365,5 2032,7 0,101 1991 7 162 446,6 3160,5 426,51 1623,1 57,39 2423,9 4120,9 0,1349 1995 8 154 700 1432,8 188,73 1586,6 76,77 1085 1892 0,1317 1996 7 85 367,5 2225,8 354,7 795,8 51,96 1624,4 3049,8 0,1594 1997 8 93 343,7 2105,1 354,22 943,87 64,27 1510,4 2934 0,1683 1998 8 84 345 1962,7 352,13 854 62,03 1377,2 2797,2 0,1794 1999 4 184 360,5 6066,6 781,23 1727,1 42,07 4710,2 7813,6 0,1288 1999 4 170 360,5 7885,1 1101,6 350,7 29,9 5992,4 10375 0,1397 2000 3 233 315 4199,3 492,82 2317 88,45 3335,1 5287,5 0,1174 2001 5 187 322 5527,2 655,53 1797,8 52,54 4377,5 6978,8 0,1186 2002 10 144 329,7 3523,5 466,16 1397,2 61,98 2715,7 4571,6 0,1323 2003 3 175 326,9 3601,7 434,66 1714,33 74,32 2840,8 4566,5 0,1207 2004 6 98 322 2614,8 384,94 934,41 58,2 1955,3 3496,8 14,72 2005 3 154 306,6 4582,9 934,97 -- --- 1808,5 11613 20,4 2006 3 117 319,9 4098,6 586,6 -- -- 3091,5 5433,7 14,3 2008 3 36 177,1 4380,4 686,69 301,99 23,2 3187,9 6019 15,68
La densidad de P. lilfordi lilfordi registrada en 2015 en la isla del Aire es muy
similar a la estimada en 2004 y corresponde a los años de densidades medias, pues en
los años 90 del siglo XX hubo períodos de densidades menores y en otras décadas las
densidades han sido mayores (Tabla 1). La fluctuación de la densidad en Aire,
probablemente gobernada por la propia fluctuación de los recursos tróficos, se mueve
siempre en torno a cifras por encima de los 1500 ejemplares/hectárea. Se trata pues de
una población en excelentes condiciones de conservación, siempre con elevadas
densidades y que se mantienen, con fluctuaciones normales, a lo largo de los años.
No obstante, la presión humana sobre la isla es evidente. No poseemos
información de cómo eran las densidades de lagartijas antes de que la isla fuera habitada
y tampoco cuando la misma albergaba a los fareros y sus familias, así como a una cierta
cantidad de ganado doméstico que incluía cabras y otras especies. también ignoramos la
influencia que puedan ejercer de modo indirecto especies como el conejo, introducido
periódicamente como pieza cinegética.
Addaia gran. En esta isla, los mejores resultados se han obtenidos por medio de tres
largos transectos lineales que recorren la isla en el sentido de su eje mayor.
11
Coefficients:Sigma 0.5171111
Function: HALFNORM
Strip: 0 to 2 Effective strip width: 0.6480315 Scaling: g(0) = 1
Log likelihood: 5.878108 AIC: 7.923563
Abundance estimate: 1176.768 individuos/hectárea; 95% CI=( 1032.961 to 1436.731 )
CI based on 495 of 500 successful bootstrap iterations.
0.0
0.5
1.0
1.5
0 0.5 1 1.5 2
Distance
Pro
babi
lity
of d
etec
tion
halfnorm, 0 expansions
ESW = 0.648
BEST FITTING FUNCTION
Figura 2. Curva de detectabilidad para la totalidad de la isla de Addaia gran (abril 2015)
12
Tabla 1. Densidades en Addaia gran a lo largo de los años
Año Mes Contactos Longitud Densidad SE AIC ESW D LCL D UCL D CV
1982 4 143 177 2355,3 2018 1740,7 171,51 0,18956 29264000 0,8568
1987 8 60 533,9 582,43 96,191 648,06 96,48 405,92 835,7 0,1652
1987 8 60 533,9 622 109,74 -- -- 427,28 905,45 0,1764
1988 8 26 315 571,19 142,64 -- -- 343,81 948,95 0,2497
1989 6 53 420 1452,9 275,81 508,92 43,43 995,25 2121 0,1898
1991 7 103 749 1095,3 118,79 212,82 62,78 811,25 1478,9 0,1085
1997 3 66 490,7 1306,6 619,57 647,25 51,47 4,2746 399360 0,4742
1997 9 59 672,7 937,93 131,84 541,09 46,76 708,57 1241,5 0,1406
1999 3 197 950,6 2925,8 769,9 1783,9 35,42 1426,4 6001,3 0,2631
1999 4 106 770 2297,1 887 924,98 29,97 461,9 11424 0,3861
1999 6 36 245 1207,5 688,08 468,31 44,39 123,31 11824 0,5698
2000 8 19 742 277,42 80,114 172,22 46,15 144,44 532,81 0,2888
2002 4 80 1108,1 1062,4 391,81 673,47 33,98 228,58 4938,2 0,369
2002 7 104 1248,4 1227,3 150,32 928,94 33,93 952,2 1581,8 0,1225
La densidad detectada en 2015 para Addaia gran es muy similar a la encontrada
en 2002, 1999 y otros años anteriores. Existen, sin embargo, algunas fluctuaciones,
como la detectada en 1988 y 1987. La densidad actual se sitúa entre los valores
habituales e la isla, de modo que podemos considerar que la población de lagartijas de
Addaia gran está relativamente estabilizada y que su estado de conservación es bueno.
No se han detectado amenazas recientes sobre la misma, si bien habría de controlarse el
acceso de los visitantes a la parte superior de la isla e incluir, como tarea rutinaria del
Parque Natural de s'Albufera des Grau, al que pertenece la isla, la limpieza anual de la
costa occidental de la isla, en la cual se acumulan desechos inorgánicos como
recipientes plásticos, que son trampas mortales para las lagartijas.
Addaia petita. Al igual que en el caso anterior, se propuso llevar a cabo un total de tres
largos transectos lineales que cubran la totalidad de la superficie insular. Sin embargo,
en este caso hemos observado en años anteriores que los transectos no parecen permitir
un número elevado de contactos, de modo que es posible que hubiera que completarlos
con un estudio de captura-recaptura.
Tabla 2. Densidades detectadas en Addaia petita en 2002 y 2005
Año Mes Contactos Longitud Densidad SE AIC ESW D LCL D UCL D CV 2002 10 21 490 289,14 94,72 212,73 74,11 148,21 564,09 0,3276 2005 9 28 742 412,57 113,35 252,19 45,73 189,19 899,69 27,47
13
En Addaia petita no obtuvimos resultados positivos en los transectos lineales
llevados a cabo y tampoco tuvimos un número suficiente de capturas para poder
establecer un programa de captura-recaptura. Los transectos se llevaron a cabo en junio
de 2015, con temperatura excelente y condiciones óptimas. Pero en un recorrido de 207
metros obtuvimos únicamente dos contactos. Es obvio que la densidad en el islote es
extremadamente baja. Sólo capturamos y liberamos posteriormente dos individuos. En
años sucesivos, el único modo de estimar la densidad de población en este islote será
dedicar varios días consecutivos a un programa intensivo de captura-recaptura. La
población actual de lagartijas es extremadamente reducida e ignoramos porqué. Es
urgente averiguar si dicha población está sufriendo una reducción progresiva y cuales
son las causas de la misma, si es el caso.
ses Aguiles. Se trata de un islote muy peculiar, con una muy escasa cobertura vegetal y
un asombroso polimorfismo de los individuos que lo habitan. Proponemos un estudio
anual de captura-recaptura en 4 sesiones consecutivas. Tenemos una escasa información
de esta población que, sin embargo, resulta muy interesante, dadas las condiciones
extremas en que viven las lagartijas.
Las condiciones del mar no permitieron en este caso más que dos visitas a ses
Aguiles en las cuales hemos podido comprobar que la población de lagartijas sigue
presente en el islote. Sin embargo no hemos podido acumular recapturas que nos
permitieran un estudio de captura-recaptura. La población debe ser extremadamente
escasa, pero parece ocupar la totalidad de la superficie insular, hasta la cumbre del
islote.
en Carbó. En este islote se han probado en años anteriores dos transectos lineales
paralelos en la totalidad de la superficie insular. Se propuso su realización, con el
complemento de un estudio de captura-recaptura de 4 sesiones consecutivas.
Tabla 3. Densidades de población estimadas en el islote de En Carbó en 2002 y 2004
Año Mes Contactos Longitud Densidad SE AIC ESW D LCL D UCL D CV 2002 4 8 194,6 475,33 261,43 73,97 43,24 135,11 1672,2 0,55 2004 8 8 121,8 595,26 308,22 78,72 55,17 180,63 1961,6 51,78
Aunque en este islote hemos podido realizar la captura de algunos ejemplares y
comprobar que la población de lagartijas parece en un buen estado de conservación, ha
14
sido imposible obtener un número de contactos suficiente para establecer la densidad de
las mismas. En 2002 y 2004, con un número muy bajo de contactos, pudimos estimar
densidades bastante coherentes. Es muy probable que la densidad actual sea similar y
que esta población, a pesar de las precarias condiciones del islote, se mantenga estable a
lo largo de los años. El islote acumula grandes cantidades de desechos inorgánicos que
deberían retirarse periódicamente. Si bien es cierto que algunos elementos ajenos al
islote y aportados por el mar, como grandes troncos y otros desechos, actúan como
refugios para las lagartijas y su eliminación quizás no sería deseable.
en Carbonet. Muy recientemente, se ha descubierto una escasa población de lagartijas
en este islote. Ignoramos si proceden de en Carbó, porque en numerosas visitas
anteriores no se habían detectado. Se propone un estudio de captura-recaptura de cuatro
sesiones consecutivas, complementado con la realización de un total de cinco transectos
puntuales desde cinco lugares preestablecidos del islote.
Hemos visitado en dos ocasiones en Carbonet, pero las observaciones y capturas
no permiten una estimación adecuada de la densidad de población. Realmente, estimar
la densidad de población en un islote de tan reducidas dimensiones es probablemente
muy complicado, incluso en años venideros. Sin embargo, se ha podido confirmar la
presencia de lagartijas y creemos que este dato es del mayor interés y sería deseable
obtener dicha confirmación en años subsiguientes.
ses Mones. Este islote reviste una particular importancia, ya que en el mismo se han
detectado, sucesivamente primero y simultáneamente después, dos especies, la lagartija
balear y la lagartija italiana, Podarcis siculus. Este hecho plantea, obviamente, la
necesidad de un seguimiento preciso de la población y el estudio de la potencial
interacción entre ambas especies. En este caso proponemos un exhaustivo estudio de
captura-recaptura de 4 sesiones consecutivas en dos ocasiones diferentes, en primavera
y en verano, unido a recorridos lineales de la totalidad de la superficie insular, así como
un estudio de la selección del hábitat de ambas especies, el uso de los recursos tróficos y
las potenciales interacciones entre ambas. Esto implicaría al menos 8-10 visitas al islote
para completar el primer año dicho estudio y para planificar el protocolo de seguimiento
para años sucesivos.
Hemos llevado a cabo un total de tres visitas al islote en las cuales hemos
capturado varios ejemplares de Podarcis lilfordi que parece relativamente abundante,
15
especialmente en las laderas orientales de la isla. Lo más sorprendente es que, a pesar
del esfuerzo de muestreo con cubos cebados y con recorridos exhaustivos del islote, no
hemos observado ni capturado ningún ejemplar de Podarcis siculus. Es prematuro
afirmar que a lagartija italiana ha desaparecido del islote, pero lo cierto es que su
presencia no es dominante en el mismo. Es obvio que hemos de continuar el
seguimiento de esta población en años venideros, pues reviste el máximo interés.
Hemos solicitado un pequeño proyecto adicional a la Societas Herpetologica Europaea
para estudiar específicamente esta población. Esto permitiría abordar dicho estudio de
un modo mucho más pormenorizado en el año 2016.
Binicodrell gros. Dos transectos lineales paralelos que cubren la totalidad de la
longitud mayor del islote. Los datos obtenidos en años anteriores se obtuvieron con este
mismo diseño.
Coefficients: Sigma 0.6607104
Function: HALFNORM Strip: 0 to 2 Effective strip width: 0.8260323
Scaling: g(0) = 1 Log likelihood: 2.075368
AIC: 4.875368
Abundance estimate: 151.3258 individuos/hectárea; 95% CI=( 106.8811 to 273.4772 ).
16
0.0
0.2
0.4
0.6
0.8
1.0
0 0.5 1 1.5 2
Distance
Pro
babi
lity
of d
etec
tion
halfnorm, 0 expansions
ESW = 0.826
BEST FITTING FUNCTION
Figura 3. Curva de detectabilidad en Binicodrell gros (2015)
Tabla 4. Densidades detectadas en Binicodrells gros a lo largo de los años
Zona Año Mes Contactos Longitud Densidad SE AIC ESW D LCL D UCL D CV
total 1987 8 55 220 1003,9 351,32 606,38 124,52 13,365 75401 0,35
total 1991 7 20 175 489,02 143,07 201,32 116,85 268,43 890,88 0,2926
total 1996 4 45 224 1188 216,99 445,01 84,55 544,89 2590,2 0,1827
total 1997 9 14 226,8 828,89 309,65 120,23 37,24 377,13 1821,8 0,3736
Oeste 1999 7 11 217 330,62 82,552 102,62 76,66 189,55 576,68 0,2497
Este 2002 7 11 272,3 202 78,01 103,31 99,99 88,01 463,61 0,3862
Oeste 2004 8 24 210 571,43 147,45 223,05 100 337,96 966,17 0,258
total 2006 8 34 333,4 1180,5 264,45 313,61 43,13 751,76 1853,8 22,4
En Binicodrell gros las densidades de población parecen fluctuar de un modo
espectacular a lo largo de los años. Obsérvese en la Tabla 4 que los valores actuales de
2015 son los más bajos registrados hasta la fecha, pero que en 2002 también hubo
densidades muy reducidas, mientras que en 1987, 1996 y 2006 las densidades superaban
17
los 1000 individuos por hectárea. Es pues complicado afirmar que la población se halla
en declive o describir la tendencia de la misma. Es necesario continuar este seguimiento
en los años venideros, para establecer si las fluctuaciones extremas son, también en este
caso, una de las características intrínsecas a estos pequeños islotes o bien, asistimos a
una erosión progresiva de los efectivos poblacionales.
Binicodrell petit. Las reducidas dimensiones de este islote y la densa cobertura
arbustiva impiden llevar a cabo transectos lineales. Se propuso un estudio de captura-
recaptura con 3-4 sesiones consecutivas y el empleo de trampas de caída cebadas con
fruta. Nuestros datos preliminares indican que en este pequeño islote hay una población
de lagartijas de al menos 20 individuos.
Hemos llevado a cabo dos visitas al islote en las cuales hemos podido
comprobar la existencia de lagartijas y su densidad relativamente elevada, aunque se
traduce en unos efectivos muy reducidos, dada la exigua extensión del islote. La
población no sufre ninguna amenaza apreciable, si exceptuamos la inherente a su propio
tamaño.
Bledas. La población de lagartijas de este islote tienen un carácter extremadamente
huidizo. Prousimos llevar a cabo dos largos transectos paralelos que abarcaran la
totalidad de la superficie insular en las laderas occidental y oriental de la isla. Se
propuso llevar a cabo un estudio de captura-recaptura con trampas de fruta en la zona
central de la isla, con 3-4 sesiones consecutivas.
Coefficients: Beta 0.3785486
Function: NEGEXP Strip: 0 to 2
Effective strip width: 1.402659 Scaling: g(0) = 1
Log likelihood: 8.035199
AIC: 10.39884
Abundance estimate: 2648.031 ; 95% CI=( 2099.196 to 5009.624 )
CI based on 338 of 500 successful bootstrap iterations
18
0.0
0.2
0.4
0.6
0.8
1.0
1.2
1.4
0 0.5 1 1.5 2
Distance
Pro
babi
lity
of d
etec
tion
negexp, 0 expansions
ESW = 1.403
BEST FITTING FUNCTION
Figura 4. Función de detectabilidad en Bledas (2015)
La visita llevada a cabo en 2015 permitió establecer que la población de
lagartijas de Bledas sigue presente en el islote y que posee unos efectivos
probablemente reducidos. La erosión de las laderas de la isla es intensa en los últimos
años y ha impedido llevar a cabo un transecto lineal adecuado. Sólo hemos podido hacer
un transecto en la ladera occidental de la isla. Por otro lado, se intentó aplicar un
programa de captura-recaptura para estimar las densidades actuales, pero no fue posible
capturar y recapturar un número suficiente de individuos. Es necesario repetir los
esfuerzos de captura-recaptura en años sucesivos, para tratar de obtener una
información más fiable. De cualquier modo, la estimación de densidad arroja una cifra
superior a la de estimaciones anteriores, lo que indicaría que esta población, muy
raramente visitada, posee unos efectivos adecuados.
19
Tabla 5. Densidad detectada en la isla de Bledas en 1997
Año Mes Contactos Longitud Densidad SE AIC ESW D LCL D UCL D CV
1997 3 15 149,8 929,52 362,19 139,04 53,86 412,59 2094,1 0,3897
1997 3 13 100,8 1119,8 461,34 126,93 57,59 468,42 2676,8 0,412
1997 3 28 250,6 1060,8 248,07 265,58 52,66 634,37 1773,8 0,2338
Colom. En esta isla se propuso la realización de 10-15 transectos lineales. Cinco de
ellos en las zonas adyacentes a las dos playas de la isla, las zonas más antropizadas y
con mayor presión humana. El resto de transectos se reparten en toda la zona central de
la isla (camino central y era), así como en todo la zona septentrional y oriental, que
apenas recibe visitas humanas.
Coefficients: Sigma 0.5701582
Function: HALFNORM
Strip: 0 to 2 Effective strip width: 0.7142643 Scaling: g(0) = 1
Log likelihood: 5.618123 AIC: 7.739335
Abundance estimate: 541.8119 individuos/hectárea; 95% CI=( 429.8651 to 672.7353 ).
20
0.0
0.2
0.4
0.6
0.8
1.0
1.2
0 0.5 1 1.5 2
Distance
Pro
babi
lity
of d
etec
tion
halfnorm, 0 expansions
ESW = 0.714
BEST FITTING FUNCTION
Figura 5. Curva de detectabilidad en la totalidad de la isla de Colom
Tabla 6. Densidades de la isla de Colom a lo largo de los años
Año Mes Contactos Longitud Densidad SE AIC ESW D LCL D UCL D CV
1982 4 43 52,9 1266,8 233,59 535,49 320,84 875,86 1832,2 0,1844
1987 8 30 238,6 859,9 200,08 301,82 73,11 537,25 1376,3 0,2327
1988 8 38 342,3 538,83 123,02 410,16 103,01 341,12 851,11 0,2283
1989 6 69 388 465,23 88,305 783,9 191,13 319,55 677,33 0,1898
1991 7 35 435,3 400,39 193,89 375,97 100,41 1,174 136540 0,4843
1996 5 90 700,7 846,97 162,55 923,16 75,83 561,6 1277,4 0,1919
1997 3 76 290,5 1047,2 276,51 802,66 124,91 537,31 2041,1 0,264
1997 8 53 200,2 1970,5 232,42 520,82 67,18 1556,5 2494,7 0,1179
1999 3 78 446,6 694,53 107,44 812,59 125,74 425,73 1133 0,1547
2001 10 70 402,5 1396 231,39 640,52 62,29 982,9 1982,8 0,1657
2002 74 608,3 809,67 187,07 728,38 75,12 483,36 1356,3 0,231
2003 2 46 431,9 410,19 192,2 481,37 129,82 7,46 22540 0,4686
2005 3 166 899,5 1615,5 294,93 1076,9 2423,7 18,26
21
Aparentemente, la isla de Colom tiene una población de lagartijas cuyos
efectivos fluctúan en ciclos plurianuales relativamente largos. En general, las
densidades no superan los 1000 individuos por hectárea, como suele ser habitual en
islas de una cierta extensión en el Mediterráneo occidental (Pérez-Mellado et al., 2008 y
referencias incluidas). De cualquier modo, los valores de densidad actualmente
detectados son comparables a los que hubo en 1988, 1989, 1991 0 2003, en la zona baja
de las fluctuaciones. No parece detectarse ningún factor actuar que afecte de modo
significativo a las lagartijas, cuyo estado de conservación puede considerarse bueno.
Sería deseable continuar este seguimiento y establecer una zonificación más precisa de
toda la isla, tratando en el futuro de relacionar las fluctuaciones poblacionales con
factores como la disponibilidad de alimento, a través de la precipitación anual u otro
indicador similar.
Rey. La isla del Rey ha sido objetivo de seguimiento desde hace casi dos décadas. Se
propone, con periodicidad anual, la realización de 5-10 transectos lineales que abarquen
la totalidad de la superficie insular, desde las zonas periféricas de vegetación natural,
especialmente en la zona norte de la isla, hasta la zona de la basílica paleocristiana,
pasando por el recinto central del patio del hospital.
Coefficients: Sigma 1.004361
Function: HALFNORM
Strip: 0 to 2 Effective strip width: 1.200312 Scaling: g(0) = 1
Log likelihood: 27.25726 AIC: 29.34422
Abundance estimate: 110.0732 individuos/hectárea; 95% CI=( 93.73289 to 136.2791 ).
22
0.0
0.2
0.4
0.6
0.8
1.0
1.2
1.4
0 0.5 1 1.5 2
Distance
Pro
babi
lity
of d
etec
tion
halfnorm, 0 expansions
ESW = 1.2
BEST FITTING FUNCTION
Figura 6. Curva de detectabilidad en la isla del Rey (2015)
Tabla 7. Densidades de población en la isla del Rey a lo largo de los años
Año Mes Hora Contactos Longitud Densidad SE AIC ESW D LCL D UCL D CV 1987 8 12 36 375,9 253,56 47,828 414,96 188,85 173,47 370,62 0,1886 1988 8 10 73 213,2 903,15 93,205 861,27 189,56 735,62 1108,8 0,1032 1989 6 8 96 96,8 5156 630,68 1007,5 96,17 4047,8 6567,4 0,1223 1991 7 7 151 476 2348,4 357,95 1553,7 67,54 1740,6 3168,5 0,1524 2002 10 10 26 368,2 307,3 92,23 260,91 114,89 86,85 1087,3 0,3001 2006 8 8 29 722,4 402,66 75,65 271,12 49,85 271,31 597,6 18,79
La densidad de población de la lagartija balear ha disminuido de modo evidente
en la isla del Rey, especialmente a partir de 2002. Nuestros resultados de 2015
confirman pues la tendencia observada en 2002 y 2006 e indican que en esta población
se ha producio un declive muy claro de la población de P. lilfordi. Dicho declive
coincide con el inicio de actividades humanas en la isla, después de varias décadas de
abandono. En la isla la vegetación autóctona se recuperaba lentamente, en parte gracias
a la acción dispersora de la lagartija balear, particularmente sobre determinadas especies
como Phillyrea media y Rubia peregrina. En la actualidad, las obras de reconstrucción
del hospital, el ajardinamiento de la zona central del patio del mismo y el resto de
23
actuaciones llevadas a cabo ha provocado la concentración de la población en áreas
relativamente subóptimas del oeste y norte de la isla, así como del recinto de la basílica
paleocristiana. El camino perimetral del islote es otro de los lugares en los cuales las
lagartijas sobreviven con concentraciones adecuadas, pero sus densidades han caído
dramáticamente en estos últimos diez años. Es obvio que la población de Rey no se
encuentra al borde la extinción, porque la extensión del islote permite la supervivencia
de la lagartija en determinados enclaves, pero se han roto algunos mutualismos del
máximo interés con varias plantas vasculares y se han eliminado numerosos refugios
que las lagartijas tenían en los antiguos muros del hospital, ahora sellados con cemento.
Es imperativo que las futuras acciones llevadas a cabo en la isla tengan en
consideración la supervivencia de la lagartija balear y preserven al menos toda la zona
norte en su estado natural. Sería deseable situar acúmulos de rocas en zonas del este y
oeste de la isla para incrementar el número de refugios de las lagartijas y suplir así la
ausencia de los anteriores refugios eliminados en los edificios. Por último, ha de
evitarse a toda costa la entrada de animales domésticos en la isla, así como el control de
roedores por medio de venenos.
Sanitja (o Porros). En este islote remoto se propuso llevar a cabo anualmente una
estimación de densidad por medio de un total de 5 transectos lineales paralelos que
cubran, en recorridos SE-NO, la totalidad de la superficie insular.
Coefficients:Sigma 0.4475044
Function: HALFNORM
Strip: 0 to 2 Effective strip width: 0.5608591 Scaling: g(0) = 1
Log likelihood: -3.448174 AIC: -1.352936
Abundance estimate: 353.9253 individuos/hectárea; 95% CI=( 255.8906 to 491.7978 )
CI based on 499 of 500 successful bootstrap iterations.
24
0.0
0.5
1.0
1.5
0 0.5 1 1.5 2
Distance
Pro
babi
lity
of d
etec
tion
halfnorm, 0 expansions
ESW = 0.561
BEST FITTING FUNCTION
Figura 7. Curva de detectabilidad en la isla de Sanitja (2015)
Tabla 8. Variación de la densidad de lagartijas baleares en la isla de Sanitja
Año Mes Contactos Longitud Densidad SE AIC ESW D LCL D UCL D CV 1987 8 80 197 2458,5 327,78 816,42 82,59 1887,6 3202,1 0,1333 1988 8 41 221,1 2490,4 1300,1 360,26 37,23 301,23 20589 0,5221 1991 7 57 203 2126,1 772,51 553,22 66,04 24,229 186560 0,3634 1991 7 90 280 2497,7 227,38 876,67 64,35 2052,6 3039,3 0,091 1992 7 58 280 1583 354,77 542,77 65,43 95,051 26365 0,2241 1992 8 264 563,5 5169,1 682,64 698 45,32 3987,8 6700,3 0,1321 1993 8 63 435,4 1535,2 407,9 587,74 47,12 668,66 3524,9 0,2657 1995 7 85 317,1 2029,3 941,3 880,75 66,05 303,71 13560 0,4639 1996 7 50 396,2 420,7 54,856 503,06 149,99 322,93 548,07 0,1304 1997 7 30 133 1716,3 509,88 292,5 65,71 945,97 3114 0,2971 1999 7 22 162,4 846,73 330,62 220,01 80 385,91 1857,8 0,3905 2000 3 74 326,9 1622,6 312,92 717,67 68,81 1044,3 2521 0,1929 2001 6 21 222,6 898,33 203,52 200,27 52,51 562,44 1434,8 0,2266 2002 5 42 261,1 1608,8 270,66 392,16 49,99 1124,3 2302 16,82 2005 9 72 550,2 2623,3 1577,9 606,05 24,94 9,33 730000 0,6015 2006 7 12 382,2 598,52 250,31 96,25 26,23 244,67 1464,1 41,82
La densidad de P. lilfordi registrada en Sanitja en 2015 es la más reducida de
todas las densidades estimadas desde 1987. Es cierto que en 1996 se produjo también
una drástica disminución de efectivos, recuperados en años subsiguientes. Después, en
25
2006, se produjo otra acusada caída de la densidad de lagartijas que, aparentemente,
persiste actualmente. El hecho de que existan algunas fluctuaciones evidentes en años
anteriores nos impide extraer conclusiones de la reducida densidad observada en 2015,
pero plantea el interrogante de qué está sucediendo en esta población. las lagartijas de
Sanitja fuerom las primeras a nivel mundial en las cuales se demostró
experimentalmente la capacidad de un reptil para polinizar una planta vascular. Después
se han acumulado evidencias empíricas en otras poblaciones, pero en ninguna hemos
observado una interacción tan intensa entre una lagartija y una planta vascular, el hinojo
marino, Crithmum maritimum. Además, la subespecie exclusiva del islote, P. lilfordi
fenni posee rasgos diferenciales muy claros, tanto a nivel morfológico como molecular.
Por todo ello, consideramos esencial la conservación de esta extraordinaria población y
proponemos un seguimiento anual de la misma y un estudio exhaustivo de las posibles
causas de su actual declive.
Sargantana. En esta isla se propone llevar a cabo un total de tres transectos lineales en
el sentido del eje mayor de la isla, que permitirá una estimación de densidad de la
totalidad de la isla, comparable también a la llevada a cabo en años anteriores. Se
propone una periodicidad anual.
Coefficients: Sigma 0.5150385
Function: HALFNORM
Strip: 0 to 2 Effective strip width: 0.6454385 Scaling: g(0) = 1
Log likelihood: 5.337515 AIC: 7.384574
Abundance estimate: 848.2824 individuos/hectárea; 95% CI=( 719.385 to 1012.718 )
CI based on 494 of 500 successful bootstrap iterations.
26
0.0
0.2
0.4
0.6
0.8
1.0
1.2
1.4
0 0.5 1 1.5 2
Distance
Pro
babi
lity
of d
etec
tion
halfnorm, 0 expansions
ESW = 0.645
BEST FITTING FUNCTION
Figura 8. Curva de detectabilidad en la isla de Sargantana (2015)
Tabla 9. Densidades detectadas en Sargantana en años anteriores
Año Mes Contactos Longitud Densidad SE AIC ESW D LCL D UCL D CV
1987 8 66 306 1056,5 347,26 704,51 102,07 381,25 2927,9 0,3287
1991 7 81 411,6 995,13 280,79 820,08 98,88 302,47 3274 0,2822
1993 8 25 233,8 486,35 134,5 250,35 109,93 277,73 851,68 0,2766
1997 3 50 212,8 997,43 228,62 527,08 117,78 632,84 1572,1 0,2292
1997 5 19 193,9 401,87 162,8 102,09 70,58 168,6 957,84 0,4051
1997 5 11 114,1 313,75 142,25 78,57 111,74 112,13 877,92 0,4566
1997 5 8 308 377,85 77,746 181,31 81,63 240,03 594,79 0,2058
2003 10 43 599,7 687,66 155,99 417,91 52,13 426,88 1107,8 0,2268
En Sargantana las densidades, desde un primer transecto en 1987, se
establecieron en menos de 1000 individiduos por hectárea. Exceptuando el primer
transecto, el resto se han llevado a cabo tras la ocupación de la isla por su propietario
actual y la presión que ha supuesto la reconstrucción de la vivienda y la actividad
27
humana en la isla. Aún así, las densidades de lagartijas se mantienen relativamente
elevadas con valores notablemente superiores este año 2015 a los detectados en 2003, a
pesar de que en aquella ocasión los transectos llevados a cabo incluían numerosos
individuos juveniles nacidos durante el verano. Por lo tanto, la densidad actual puede
considerarse buena y el estado de conservación de las lagartijas también. En Sargantana
es imprescindible controlar de modo estricto la ausencia de gatos y otros animales
domésticos. Se ha detectado la presencia de una cabra en la isla que podría afectar a la
vegetación y, por lo tanto, a los recursos tróficos disponibles para las lagartijas. Por lo
demás, el propietario actual parece concienciado de la importancia de las lagartijas y de
su conservación.
s'en Tosqueta. En este islote de superficie relativamente llana y densidad media se
puede establecer la densidad de población adecuadamente por medio de 3-4 transectos
lineales aproximadamente paralelos y que abarquen la totalidad de la superficie insular.
Se propone una periodicidad anual.
Coefficients: Sigma 0.6491517
Function: HALFNORM
Strip: 0 to 2 Effective strip width: 0.8119119 Scaling: g(0) = 1
Log likelihood: 5.89969 AIC: 8.110216
Abundance estimate: 652.8237 individuos/hectárea; 95% CI=( 447.9802 to 967.3437 ).
28
0.0
0.2
0.4
0.6
0.8
1.0
1.2
1.4
0 0.5 1 1.5 2
Distance
Pro
babi
lity
of d
etec
tion
halfnorm, 0 expansions
ESW = 0.812
BEST FITTING FUNCTION
Figura 9. Curva de detectabilidad en la isla de s'en Tosqueta (2015)
Tabla 10. Densidades en s'en Tosqueta desde 1991
Zona Año Mes Contactos Longitud Densidad SE AIC ESW D LCL D UCL D CV
Centro 1991 7 8 96,6 1035,3 556 71,78 39,99 302,1 3548,3 0,537
Centro 1996 7 14 168 483,9 162,38 130,68 86,11 238,93 980,02 0,3356
Centro 1999 7 17 78,4 5067,5 1891,2 159 23,91 2347,1 10941 0,3732
total 2002 6 14 245 285,71 98,65 130,94 100 138,38 589,94 0,3453
total 2002 7 21 197,4 1703 501,41 176,3 31,23 931,36 3113,9 0,2944
S'en Tosqueta se caracteriza por presentar densidades extremadamente variables
a lo largo de los años, quizás porque se trata de un islote de muy reducidas dimensiones
en el cual la disponibilidad de recursos fluctúa de modo notable. Las densidades
actuales con netamente menores a las detectadas en 2002 o 1991, pero la población
parece que se halla en un buen estado de conservación.
Rovells. Debido a la densa vegetación arbustiva del islote, es imposible llevar a cabo
29
transectos lineales fiables, de modo que proponemos, con una periodicidad anual, un
estudio de captura-recaptura de 3-4 sesiones consecutivas por medio de trampas de
fruta, que se colocarían en toda la zona occidental y central del islote.
Coefficients: Sigma 0.3368048
Function: HALFNORM
Strip: 0 to 2 Effective strip width: 0.4221223 Scaling: g(0) = 1
Log likelihood: -11.59873 AIC: -9.465396
Abundance estimate: 1244.785 individuos/hectárea; 95% CI=( 923.9241 to 1747.55 ).
0.0
0.2
0.4
0.6
0.8
1.0
1.2
0 0.5 1 1.5 2
Distance
Pro
babi
lity
of d
etec
tion
halfnorm, 0 expansions
ESW = 0.422
BEST FITTING FUNCTION
Figura 10. Curva de detectabilidad en la isla de Rovells (2015)
Estos son los primeros datos existentes sobre la densidad de P.lilfordi en el islote de
Rovells. En la anterior revisión de densidades (Pérez-Mellado et al., 2008) no se incluía
esta población ya que la densa vegetación arbustiva del islote impedía la realización de
transectos lineales. Sin embargo, toda la periferia de la isla se halla hoy día
30
relativamente despejada, permitiendo un recorrido perimetral de la isla y una evaluación
de la densidad precisamente en las zonas donde las lagartijas son más comunes. Una
densidad de más de 1200 individuos por hectárea es, desde nuestro punto de vista,
adecuada para un islote de las características de Rovells y permite concluir que la
población de lagartija balear se hallan en unas buenas condiciones de conservación.
Porros (de Fornells). Se trata de un islote de muy reducidas dimensiones que
permitiría, sin embargo, llevar a cabo dos transectos lineales cortos y, en jornadas de
trabajo diferentes, un estudio de captura-recaptura por medio de trampas de fruta. La
densidad de lagartijas es asombrosamente elevada, de modo que su seguimiento reviste
un gran interés.
Coefficients: Sigma 0.25
Function: HALFNORM
Strip: 0 to 2 Effective strip width: 0.3133285 Scaling: g(0) = 1
Log likelihood: -2.642012 AIC: 1.357988
Abundance estimate: 2532.967 individuos/hectárea; 95% CI=( 1791.078 to 157936.5).
31
0.0
0.2
0.4
0.6
0.8
1.0
0 0.5 1 1.5 2
Distance
Pro
babi
lity
of d
etec
tion
halfnorm, 0 expansions
ESW = 0.313
BEST FITTING FUNCTION
Figura 11. Función de detectabilidad en Porros de Fornells (2015)
Tabla 11. Densidad de P. lilfordi en el islote de Porros de Fornells en 1987 Año Mes Contactos Longitud Densidad SE AIC ESW D LCL D UCL D CV 1987 8 13 42 1676 599,31 139 92,34 781,1 3596,3 0,3576
Este año las tres visitas efectuadas al islote han permitido estimar la densidad
actual de lagartijas que, dado el escaso número de contactos, tiene una muy limitada
fiabilidad. En todo caso, dicha densidad parece similar a la estimada en 1987, con más
de 1600 individuos por hectárea (ver también Pérez-Mellado et al., 1987), ya que en una
cifra similar se situaría el límite inferior al 95% de la estimación de 2015 (ver más
arriba). De cualquier modo, se detectaron en varias ocasiones serios impactos sobre la
vegetación del islote, como consecuencia del desembarco de pescadores deportivos y
practicantes de vela ligera en el mismo. Al pisoteo de las plantas se une ocasionalmente
el impacto de fuegos de campamento en el islote. Esta población reviste un
extraordinario interés y debería preservarse de cualquier impacto exterior. En este caso,
32
sería deseable la adopción de una normativa que prohibiera el desembarco en el islote.
Mel. Este islote posee una densa cobertura arbustiva y su estudio sólo puede llevarse a
cabo por medio de 4-5 sesiones consecutivas de captura-recaptura con trampas de fruta
situadas en todo el perímetro del islote.
En este caso hemos llevado a cabo una única visita en 2015 durante la cual nos
ha sido imposible establecer la densidad de población. La densa vegetación arbustiva
impide la realización de transectos lineales, de modo que las estimaciones en años
sucesivos deberán llevarse a cabo con programas de captura-recaptura. De cualquier
modo, el islote mantiene una población en buen estado de conservación que muy
raramente recibe visitas y que se halla adecuadamente protegida en el interior del
Parque Natural de s'Albufera des Grau.
Conclusiones
En este primer año de seguimiento hemos visitado la totalidad de las poblaciones
de lagartija balear conocidas en Menorca. El esfuerzo de campo ha rendido, sin
embargo, unos resultados variables. En primer lugar, podemos confirmar que la
totalidad de las poblaciones conocidas en Menorca siguen albergando lagartijas, incluso
las más reducidas como Addaia petita, ses Aguiles, en Carbó, en Carbonet, ses Mones,
Binicodrell petit y Mel. Estas siete poblaciones son tan reducidas que los transectos
lineales y los programas de captura-recaptura se han revelado inaplicables en las
condiciones actuales. Hemos de diseñar programas específicos de transectos puntuales y
de captura-recaptura a largo plazo para estimar la densidad real. Ya hemos comenzado
con la identificación individual en base a fotografías ventrales digitales que permiten la
progresiva construcción de una base de datos sobre los individuos presentes en cada
población.
Las poblaciones de Aire, Addaia gran, Bledas, Colom, Sargantana, Rovells y
Porros de Fornells se mantienen en buenas condiciones de conservación, con efectivos
poblaciones muy elevados en Aire y Colom y poblaciones más reducidas pero densas en
el resto de islotes mencionados.
Los islotes de Binicodrell gros y s'en Tosqueta exhiben fluctuaciones extremas
de sus poblaciones de lagartijas, de modo que sólo ulterior información en los próximos
años permitirá establecer si existe alguna amenaza sobre las mismas, aunque en
principio no lo parece.
33
Finalmente, las lagartijas que habitan en las islas de Sanitja y Rey han sufrido un
fuerte disminución de sus efectivos en las últimas décadas. En el caso de Rey parece
claro que la causa es la actividad humana en el islote, que se ha incrementado
exponencialmente desde hace una década, con trabajos de reconstrucción y visitas
frecuentes. En el caso de Sanitja, las causas de la dramática disminución de su
población de lagartijas son complejas. El islote sufrió la presión ejercida por un rebaño
de cabras durante años, lo que provocó, con toda probabilidad, la drástica reducción de
su cobertura vegetal. La realmente llamativo es que, si tal efecto tuvo lugar hace ya
varios años, la cobertura no se haya recuperado plenamente y tampoco la población de
lagartijas que se vio afectada por esta disminución de recursos tróficos. Es necesario
investigar más a fondo qué causas pueden estar detrás del colapso poblacional
observado en Sanitja.
En definitiva, las 18 poblaciones estudiadas albergan lagartijas baleares, pero
muchas de ellas con efectivos muy reducidos y una situación relativamente precaria
como consecuencia, justamente, de sus bajos efectivos. Otras siete poblaciones tienen
una situación aceptable, dos más presentan fluctuaciones poblaciones que parecen
debidas a causas naturales y, por fin, dos poblaciones han visto seriamente reducidos
sus efectivos en las últimas décadas. Todo ello aconseja la continuidad de este programa
de seguimiento en los próximos años, a ser posible con una periodicidad anual en todas
las poblaciones de reducidos efectivos y con un especial seguimiento de las que han
sufrido reducciones extremas.