Introducción

La temperatura es la principal variable ambiental que influye en la sobrevivencia de los ectotermos, pues de ésta depende que se realicen correctamente sus procesos ecofisiológicos (^{Dubois et al., 2009}; ^{Huey, 1982}; ^{Huey y Slatkin, 1976}). Para que estos organismos mantengan su temperatura corporal dentro de un intervalo óptimo, deben termorregular empleando diferentes estrategias. En el caso de los reptiles, la termorregulación es principalmente conductual y en menor medida mediante ajustes fisiológicos (^{Blouin-Demers y Nadeau, 2005}; ^{Hertz et al., 1993}; ^{Huey, 1974a}; ^{Huey y Slatkin, 1976}; ^{Labra et al., 2008}). En específico, las lagartijas pueden obtener el calor mediante 2 vías: 1) tigmotermia, obtención de calor mayormente mediante conducción y 2) heliotermia, obtienen calor principalmente por radiación solar directa (^{Castellini, 2009}; ^{Cowles y Bogert, 1944}; ^{Willmer et al., 2005}). Aunado a ésto, para lograr una termorregulación adecuada, las lagartijas dependen de la temperatura disponible en distintos micrositios en su hábitat (^{Hertz et al., 1993}; ^{Sinervo et al., 2010}). La temperatura ambiental determina la cantidad de horas al día en las cuales las lagartijas pueden estar activas en función de sus requerimientos fisiológicos, dicho rendimiento se rige en respuesta directa de las temperaturas corporales y se ha estudiado mediante las curvas de desempeño (^{Huey y Slatkin, 1976}; ^{Labra et al., 2008}).

Las curvas de desempeño tienen un comportamiento general, el cual inicia a partir de la temperatura crítica mínima de los organismos, se incrementa a medida que se eleva la temperatura corporal hasta llegar a la temperatura óptima y posteriormente decrece abruptamente hasta la temperatura crítica máxima (^{Angilletta et al., 2002}). El límite crítico máximo y mínimo, así como la amplitud de cada curva, son generalmente características de cada especie (^{Angilletta et al., 2002}; ^{Currie et al., 1998}). La amplitud del rendimiento térmico es un intervalo de temperaturas en las cuales la temperatura corporal se mantiene en niveles óptimos, mediante estrategias conductuales y en ausencia de restricciones externas (e.g., depredación, competencia). En general, las curvas con gran amplitud indican que los organismos son generalistas térmicos (i.e., euritérmicos) y, por el contrario, las curvas con una amplitud estrecha nos indican organismos especialistas (i.e., estenotérmicos; ^{Hertz et al., 1993}; ^{Rowe et al., 2020}).

Los organismos tienen una amplitud del intervalo de temperaturas en el cual son capaces de desempeñar sus funciones y actividades vitales de manera óptima (^{Grigg y Buckley, 2013}; ^{Huey y Slatkin, 1976}). Cuando las temperaturas ambientales son más altas que el límite crítico máximo de dicho intervalo, pueden incrementar el riesgo de muerte de los organismos por sobrecalentamiento. Para evitar sobrecalentarse, las lagartijas deben resguardarse en refugios con temperaturas favorables (^{Sinervo et al., 2010}). Si los organismos pasan demasiado tiempo en dichos refugios, su aporte energético puede llegar a ser insuficiente, lo cual puede causar efectos poblacionales negativos, especialmente en la época de reproducción, ya que puede causar interrupción en el cortejo y disminución de la energía destinada a la reproducción, lo que podría suponer consecuencias demográficas negativas para las poblaciones en años consecutivos (^{Sinervo et al., 2010}; Vicenzi et al., 2017). Por tanto, dicho aumento en la “inactividad” forzada (también conocido como horas de restricción térmica), principalmente en la época de reproducción, podría causar la extinción local de hasta el 39% de las poblaciones de lagartijas en México para el año 2080, como consecuencia del incremento en la temperatura (^{Sinervo et al., 2010}). Debido a lo anterior, para evaluar la sensibilidad fisiológica y el riesgo de extinción o resiliencia de las especies de lagartijas ante el cambio climático, es necesario estudiar su ecología térmica y el desempeño térmico de diferentes poblaciones (^{O’Connor et al., 2006}).

Así como el cambio climático es una amenaza en contra de la biodiversidad, también altera los patrones meteorológicos, ya sea aumentando la probabilidad o severidad de eventos extremos (e.g., sequías, ondas de calor, heladas) ligados al calor, que a su vez, afectan los patrones de circulación del aire en la atmósfera (^{Stott, 2016}). Particularmente, el cambio climático ha provocado alteraciones a fenómenos oceáno-atmosféricos cíclicos como El Niño-Oscilación del Sur (ENSO, por sus siglas en inglés), provocando que su fase cálida, conocida como “El Niño” sea aún más drástica (^{Snell y Rea, 1999}). El fenómeno de El Niño se asocia con temperaturas altas en la superficie del mar en el área este del Pacífico tropical, entre las costas de Perú y Ecuador, lo que causa una disminución en la intensidad de los vientos alisios y grandes evaporaciones de agua a la atmósfera, lo cual genera lluvias y fuertes tormentas en Sudamérica, además de alterar los efectos de la corriente fría de Humboldt (^{Snell y Rea, 1999}). Como consecuencia de dichas fluctuaciones climáticas, las interacciones ecológicas se ven fuertemente afectadas en muchas especies marinas y terrestres, principalmente de la región costera del Pacífico oriental ecuatorial, incluyendo el archipiélago de Galápagos (^{Snell y Rea, 1999}).

Se ha documentado que el fenómeno de El Niño puede acarrear consecuencias ecológicas en las islas Galápagos (^{Snell y Rea, 1999}). Por ejemplo, poblaciones de aves y mamíferos marinos cuya distribución límite es ecuatorial, han experimentado mortandades alarmantes, las cuales se pueden exacerbar con las actividades humanas o por el impacto de especies introducidas. Particularmente, las islas Galápagos presentan un alto grado de endemismos tanto en su flora como en su fauna, las cuales están altamente especializadas a las condiciones ambientales locales (^{Benavides et al., 2009}; ^{Muñoz-Barriga, 2015}). Tal es el caso de las lagartijas de lava (Microlophus), uno de los grupos mejor representados en el archipiélago. Cada una de las islas principales tiene una especie diferente de lagartija, con poblaciones que se distribuyen exclusivamente en la zona costera y xerofítica (^{Benavides et al., 2009}; ^{Kizirian et al., 2004}). El objetivo de esta investigación fue estudiar la ecología térmica de la lagartija de lava de Santa Cruz, Microlophus indefatigabilis. Particularmente, buscamos evaluar el desempeño locomotor, para así calcular su umbral térmico de actividad y las posibles horas de restricción térmica de la especie ante el aumento de la temperatura ambiental, efectos que podrían presentarse ante la incidencia del fenómeno de El Niño.

Materiales y métodos

La lagartija de lava de Santa Cruz (Microlophus indefatigabilis), es una especie endémica de las islas Galápagos que se distribuye en las islas Santa Cruz, Santa Fe, Baltra y en pequeños islotes asociados. Al igual que el resto de las especies del género Microlophus, esta especie presenta un marcado dimorfismo sexual, los machos son más robustos y de mayor tamaño que las hembras. Los machos presentan una combinación de manchas negras y anaranjadas en la garganta y en los flancos, así como una cresta dorsal que va desde la cabeza hasta la cola. Las hembras exhiben un color naranja en la cabeza, garganta y vientre, el cual es muy conspicuo durante la temporada reproductora (^{Stebbins et al., 1967}). En este estudio consideramos a los organismos como adultos a partir de 59 mm de longitud hocico-cloaca (LHC) en machos y a partir de 57 mm de LHC en hembras (^{Stebbins et al., 1967}).

Realizamos el estudio durante la estación fría, entre el 31 de julio y el 6 de agosto de 2017. Recolectamos organismos de M. indefatigabilis en 3 poblaciones de 2 islas. Una población en la isla Baltra (00º27’00” S, 90º16’00” O) y 2 poblaciones en la isla Santa Cruz: Punta Estrada (00°45’14.0’’ S, 90°18’41.8’’ O) y la Estación Científica Charles Darwin (referida de aquí en adelante como Estación Científica, 00º44’34.81” S, 90º18’12.40” O). Las 3 zonas de estudio presentan montículos de rocas basálticas de distintos tamaños y vegetación característica de zonas áridas de tierras bajas, con especies como monte salado (Cryptocarpus pyriformis), cactus (Opuntia echios), sesuvium (Sesuvium edmondstonei), parkinsonia (Parkinsonia aculeata), palo santo (Bursera graveolens) y el árbol de la muerte (Hippomane mancinella), entre otras (^{Snell et al., 1995}; ^{Stebbins et al., 1967}).

Recolectamos individuos de M. indefatigabilis durante su periodo de actividad (07:00 -19:00 h, GMT-6), mediante cañas de pescar modificadas con un nudo corredizo en el extremo. Registramos su temperatura corporal (cloacal; T_c) mediante un termómetro digital de lectura rápida (Fluke modelo 51-II), con un termopar tipo T (± 0.1°C) y descartamos aquellas lecturas en que demoramos más de 20 segundos en capturar al individuo. También registramos la temperatura del sustrato (Ts) de donde fue recolectado cada ejemplar y la temperatura del aire (T_a) a la sombra, a 5 cm de altura, así como la hora de recolecta (^{Hertz et al., 1993}).

En la población de la Estación Científica registramos la temperatura operativa (T_o) mediante lectores de temperatura (HOBO®Pro v2). Colocamos 2 registradores con 2 sensores, respectivamente. Cada sensor estuvo conectado a un modelo nulo (ver más adelante). También registramos la temperatura del aire (T_air) con un tercer registrador de temperatura y humedad (HOBO®Pro v2 U23-002), con un sensor descubierto a metro y medio de altura. Finalmente, en la población de Baltra registramos la T_o mediante 20 registradores de temperatura (iButton®), cada uno de los cuales fue colocado dentro de un modelo nulo.

El modelo nulo fue diseñado para emular las tasas de ganancia y pérdida de calor de los organismos en ausencia de termorregulación y representan la T_o en diferentes microhábitats (^{Bakken, 1992}; ^{Dzialowski, 2005}). Para calibrar los modelos nulos, sometimos a 2 organismos del género Microlophus, junto con diferentes posibles modelos a intervalos de luz artificial para incrementar su temperatura, y periodos de sombra para disminuirla. Posteriormente, analizamos las temperaturas de los organismos y los modelos nulos por medio de una correlación lineal, y seleccionamos el modelo con el coeficiente de correlación (r) más alto (^{Dzialowski, 2005}; ^{Hertz et al., 1993}). El modelo seleccionado consistió en un tubo de policloruro de vinilo (PVC) relleno de silicón y pintado de gris, mismo que presentó una alta correlación con el organismo de referencia (r = 0.93; p < 0.001). Colocamos los modelos nulos en diferentes microhábitats donde previamente observamos individuos de M. indefatigabilis, tomando en cuenta condiciones de sol y sombra. Los modelos de la Estación Científica registraron la temperatura cada 30 minutos del 4 al 31 agosto de 2017, mientras que los modelos de Baltra registraron la temperatura cada 10 minutos entre el 5 y el 9 de agosto de 2017.

Adicionalmente, obtuvimos las temperaturas ambientales promedio (T_air) y máximas al día (T_max) de la estación meteorológica Estación Científica Charles Darwin, correspondientes a la estación fría de 1997, 2015 y 2017 (www.darwinfoundation.org), y la T_air de cada hora entre el 5 y el 9 de agosto de 2017 de la estación meteorológica Baltra (www.windfinder.com). Cabe destacar que solo se emplearon las temperaturas del intervalo correspondiente al periodo de actividad diario de la especie.

Para conocer las temperaturas seleccionadas de la especie (T_sel; temperaturas corporales que presenta un organismo cuando se le brinda un abanico de temperaturas en ausencia de limitantes ambientales; también llamadas temperaturas preferidas), colocamos a los organismos de la Estación Científica y de Punta Estrada en un gradiente térmico bajo condiciones de laboratorio. Construimos el gradiente térmico con láminas de policarbonato con dimensiones de 1.80 x 1.00 x 0.10 m (largo-ancho-alto) y lo dividimos en 10 carriles de 0.10 m. Generamos un extremo frío de 20 °C por medio de bolsas de hielo, colocadas dentro de los carriles y un extremo caliente de 45 °C con una hilera de focos de 100 W, colocados a 30 cm de altura. Realizamos los experimentos de gradiente térmico durante un día a lo largo del periodo de actividad de la especie observado en campo.

Brindamos a los organismos 1 hora de ambientación, previa a la primera lectura de T_sel en el gradiente, y posteriormente registramos la temperatura corporal de cada organismo en intervalos de 1 hora, con un termómetro digital de lectura rápida (Fluke modelo 51- II) con un termopar tipo T (± 0.1 °C). Al final, con todas las lecturas determinamos el 50% de las T_sel centrales, usando como intervalos de T_sel los cuartiles 25% y 75% (^{Hertz et al., 1993}). Los valores correspondientes fueron empleados como intervalos para calcular los índices de termorregulación y calidad térmica del ambiente.

Calculamos la eficiencia térmica (E) con base en el protocolo de ^{Hertz et al. (1993)}. Con los datos recabados de T_c, T_sel y T_o, obtuvimos los siguientes parámetros necesarios para el cálculo de E. La precisión en la termorregulación (d_b), la cual calculamos mediante el promedio de las desviaciones de las T_c con respecto al intervalo de T_sel, es decir, la diferencia en el número absoluto entre la T_c y T_sel, donde valores altos de d_b implican una baja precisión, mientras que valores de 0 o cercanos a 0 representan una alta precisión en la selección de la T_c del organismo.

También calculamos la calidad térmica del ambiente (de) con el promedio de las desviaciones entre cada una de las T_o correspondientes al intervalo de actividad y el intervalo de T_sel, los valores altos indican baja calidad térmica del ambiente, y valores cercanos a 0 representan que el ambiente es ideal para los organismos. Si el valor de T_c o de T_o están por debajo del intervalo de T_sel, se contrastan con el valor inferior del intervalo T_sel, de lo contrario, se contrastan con el valor superior de T_sel. Por último, obtuvimos la eficiencia de la termorregulación (E) al restar a 1 los cocientes de d_b/d_e. Si el valor final es cercano a 0, esto sugiere que el organismo presenta una tendencia al termoconformismo, mientras que los valores cercanos a 1 supondrían que el organismo termorregula activamente.

Calculamos también la eficiencia en la termorregulación de acuerdo con el índice de ^{Blouin-Demers y Weatherhead (2001)}. Dicho índice estima la magnitud de la desviación del extremo termoconformista, para esto, se resta la calidad térmica del ambiente a la precisión de la termorregulación (d_e - d_b). Valores altos significan que los organismos están termorregulando activamente, mientras valores cercanos a cero indican que los organismos están comportándose como termoconformistas (Blouin-Demers y Weatherhead, 2001).

A fin de no causar efectos de estrés térmico, determinamos las temperaturas medias críticas máximas (CT_Max) y mínimas (CT_Min) de la especie posteriormente a la toma de T_sel. Para la estimación de la CT_Min, colocamos a 10 organismos de la Estación Científica dentro de un recipiente y este sobre hielo, para reducir gradualmente su T_c. En el caso de la CT_Max, colocamos a 5 organismos de la Estación Científica debajo de una lámpara de 100 W para elevar su T_c de forma gradual. En ambos casos, los organismos fueron monitoreados y colocados sobre su dorso constantemente, y en el momento en que dejaron de reincorporarse, registramos su T_c (^{Du et al., 2000}; ^{Lutterschmidt y Hutchison, 1997}). Inmediatamente después de registrar dichos parámetros de tolerancia, llevamos a los organismos al intervalo de T_sel, y en el caso de aquellos sometidos a la CT_Max, les proporcionamos agua para rehidratarlos.

Medimos el desempeño locomotor de las 3 poblaciones a 5 distintos tratamientos de temperatura (^{Angilletta et al., 2002}). Para ésto, resguardamos durante 30 min a cada organismo en una incubadora rústica (construida en una hielera de unicel de 31 x 24 x 22 cm, con un foco de 100 W controlado por una placa de desarrollo de hardware libre Arduino^TM Mega) a la temperatura del tratamiento, realizamos el experimento cuando cada organismo llegó a la temperatura deseada ± 0.5 °C. Colocamos a cada organismo solo en uno de los 5 tratamientos térmicos: 18, 25, 30, 35 o 38 °C. Posteriormente, colocamos a cada ejemplar en el extremo de un carril de polimetilmetacrilato (1 x 0.12 x 0.12 m) y estimulamos a cada uno a correr la longitud del mismo durante 3 ocasiones. El carril tenía un sistema de diodos infrarrojos espaciados a intervalos de 20 cm, conectados a una placa de desarrollo de hardware libre Arduino^TM Mega, los cuales estuvieron programados para registrar la velocidad (cm/s) de los organismos a su paso por el carril en una computadora portátil, de manera que promediamos la velocidad obtenida de los diodos para cada carrera, y usamos la velocidad máxima de cada ejemplar. Elegimos al azar a los organismos en cada uno de los tratamientos, pero consideramos que las 3 categorías estuvieran equitativamente representadas en cada uno de ellos. De cada ejemplar registramos la LHC y el sexo.

Construimos la curva de desempeño locomotor de M. indefatigabilis mediante un modelo aditivo generalizado (GAM) empleando la LHC como covariable, la velocidad máxima de carrera de cada organismo en alguno de los 5 tratamientos de temperaturas, y el promedio de la CT_Max y de la CT_Min para anclar las curvas. A partir de la curva determinamos la temperatura óptima (T_opt), y con ésta calculamos 80% del umbral térmico (B₈₀, amplitud térmica de desempeño 80) para la especie. Las estimaciones de la curva fueron realizadas con el programa R (R Core Team, 2020) usando la paquetería minpack.lm (^{Elzhov et al., 2016}).

Calculamos la “tolerancia al calentamiento” restando el promedio de la T_o del periodo de actividad al promedio de la CT_Max. Dicho valor se aproxima a la cantidad promedio de calentamiento ambiental que una lagartija puede tolerar antes de que el rendimiento descienda a niveles fatales. También calculamos el “margen de seguridad térmica” con la diferencia entre la T_opt y el promedio de la T_o del periodo de actividad. Este valor indica el intervalo térmico en el cual las especies pueden disminuir o aumentar su rendimiento (^{Deutsch et al., 2008}).

Cabe destacar que los organismos no fueron sometidos a todos los experimentos, que las 3 poblaciones no fueron muestreadas el mismo día y que todos los experimentos los realizamos al día siguiente de la captura de los ejemplares. Para la realización de las pruebas, se acondicionó un laboratorio donde mantuvimos a los ejemplares en sacos de tela individuales. Posteriormente a los experimentos, hidratamos a todos los organismos y los liberamos en el sitio donde fueron encontrados inicialmente. Todas las pruebas se realizaron con respaldo del permiso No PC-66-17 otorgado por la Dirección del Parque Nacional Galápagos.

Estimamos los efectos de El Niño en la T_o de la Estación Científica (sitio con abundante vegetación) y de Baltra (sitio con poca cobertura vegetal) durante las fechas de muestreo. Para ello, utilizamos regresiones lineales y no lineales entre T_air y las T_o de los 2 sitios. Seleccionamos el mejor modelo para las regresiones comparando los valores de Akaike (AIC) de un modelo lineal, uno potencial y uno exponencial. Después, determinamos el promedio de la Tair correspondiente a la estación fría de 2017 (TP₂₀₁₇) y el promedio de la Tair para la misma estación de El Niño 2015, el cual corresponde al evento más cercano al muestreo, y de El Niño 1997 (TP _niño) que fue el evento más extremo del que se tiene registro en la zona (^{Snell y Rea, 1999}; ^{Yong, 2019}). Para esto, utilizamos datos climáticos históricos de la estación meteorológica de la Estación Científica y de Baltra. A continuación, calculamos el incremento en la TP (ΔTP = TP_niño - TP₂₀₁₇). Una vez que obtuvimos el ΔTP, sumamos este valor a las temperaturas del aire de la Estación Científica y de Baltra. Por medio de las regresiones, estimamos la T_o para esa misma temporada, pero bajo el supuesto de la presencia del fenómeno de El Niño. Finalmente, obtuvimos el promedio de las horas de restricción térmica (Hr, horas del día en las cuales la actividad de un organismo se ve limitada por altas temperaturas ambientales; ^{Sinervo et al., 2010}) de los sitios de estudio durante las fechas del muestreo (en 2017 y con las T_o reconstruidas durante un hipotético fenómeno de El Niño). Consideramos una Hr cuando la T_o de una hora del día supera el valor máximo de la B₈₀.

Analizamos la normalidad y homocedasticidad de los datos mediante las pruebas Kolmogorov-Smirnov y Levene. Las T_c de los organismos de la Estación Científica fueron analizadas mediante un análisis de varianza (Andeva), mientras que la Tc de los ejemplares de Baltra y el promedio de la Tsel (entre organismos de Punta Estrada y la Estación Científica) se analizaron mediante pruebas U de Mann-Whitney. Usamos la prueba Kruskal-Wallis para comparar la Tc de los ejemplares de Punta Estrada, la T_c de las 3 poblaciones y la LHC entre tratamientos de desempeño locomotor. Hicimos una prueba DHS de Tukey para conocer entre qué poblaciones existían diferencias. Presentamos el promedio de la T_c y su gris. desviación estándar, y el intervalo mínimo y máximo entre paréntesis. Realizamos correlaciones de Pearson entre las temperaturas ambientales (T_a y T_s) y la T_c. Realizamos los análisis en el programa SPSS 21 (SPSS Inc.), las gráficas se construyeron en el programa SigmaPlot 11.0 (Systat Software) y con la paquetería ggplot2 del programa R (^{R Core Team, 2020}).

Resultados

Se recolectó un total de 223 individuos de M. indefatigabilis entre el 31 de julio y el 06 de agosto de 2017. De éstos, 117 ejemplares fueron capturados en la isla Santa Cruz en 2 poblaciones: Punta Estrada (n = 66) y la Estación Científica (n = 51). Asimismo, recolectamos 116 organismos en la isla Baltra. No existieron diferencias significativas en la T_c entre categorías (machos, hembras y jóvenes) en la población de Punta Estrada (H _66,2 = 2.091, p = 0.351), en la Estación Científica (F _50,2 = 1.145, p = 0.327), ni en la isla Baltra (U_39,1 = 217.0, p = 0.408), por lo que agrupamos los datos por población para los siguientes análisis.

La T_c promedio para la especie fue 34 ± 3.79 °C con un intervalo entre 22.3 y 41.1 °C, mientras que la T_sel promedio para la especie fue de 30.5 °C. La Tc promedio para la Estación Científica fue 32.5 ± 3.34 °C, Punta Estrada 33.9 ± 4.76 °C y en Baltra 34.4 ± 3.19 °C. Encontramos que existen diferencias significativas entre la Tc de las 3 poblaciones (H_233,2 = 13.674, p < 0.05), en específico, entre los organismos de Baltra y los de la Estación Científica (F_164,1 = 1.993, p < 0.05). Por el contrario, no observamos diferencias significativas del promedio de T_sel entre la Estación Científica y Punta Estrada (U20,1 = 46.0, p = 0.762). Los organismos de las 3 poblaciones mostraron correlaciones positivas y significativas entre su T_c y las temperaturas ambientales (T_s y T_a ; fig. 1).

Figura 1 Relación de las temperaturas ambientales (temperatura del aire, Ta y temperatura del sustrato, T_s) con la temperatura corporal (T_c) de Microlophus indefatigabilis, de la población Punta Estrada (A), Estación Científica (B) y Baltra (C). Se muestra el valor de correlación de Pearson (r) y la línea de mínimos cuadrados de las rectas de regresión en gris 

En la Estación Científica registramos 2,800 lecturas de T_o, mientras que en Baltra obtuvimos 6,820 tomas entre las 07:00 y las 19:00 h (fig. 2). Con las T_o, la T_c y los cuartiles 25% - 75% de la T_sel (26.4 - 34 °C) de la Estación Científica, calculamos la eficiencia térmica de M. indefatigabilis. De acuerdo con el índice de ^{Hertz et al. (1993)}, la precisión de la termorregulación de los organismos fue alta (d_b = 0.8), en tanto que la calidad térmica del ambiente indica una calidad térmica moderada (d_e = 3.5), en consecuencia, se presenta una alta eficiencia en la termorregulación (E = 0.79). De igual forma, el índice de ^{Blouin-Demers y Weatherhead (2001)} sugiere que los organismos se esfuerzan en aumentar su Tc por lo que tienden a termorregular (d_e- d_b= 2.7).

Figura 2 Promedio de temperaturas operativas y temperatura máxima durante la estación fría en la Estación Científica (Estación) y en Baltra. Temperatura corporal promedio (T_c promedio) y cuartiles 25% y 75% de las temperaturas seleccionadas (T_sel) de Microlophus indefatigabilis. 

El promedio de la CT_Min de M. indefatigabilis fue de 12.5 °C (11 - 14.3 ºC), mientras que el promedio de la CT_Max fue de 40.6 °C (39.9 - 41.4 ºC). No hubo diferencias significativas entre la LHC de los organismos entre los tratamientos (H₄ = 3.837, p = 0.429). Para el análisis de desempeño locomotor, primero probamos diferentes modelos, con los datos de velocidad máxima con distintas particiones de datos, por ejemplo: por sitio o totales, y con o sin el uso de LHC como covariable. El mejor modelo para la especie fue descrito por un modelo de la familia Gaussiana (AIC = 149.9773). La Topt del desempeño fue de 30.4 °C, la velocidad máxima fue a los 1.54 m/s, mientras que la B₈₀ de la especie tuvo una amplitud de 11.6 °C (24.6 - 36.2 °C; fig. 3). La tolerancia al calentamiento fue de 11.3 °C y el margen de seguridad térmica de 1.09 °C.

Figura 3 Curva de desempeño locomotor de Microlophus indefatigabilis, la cual se ancló con la temperatura media crítica mínima (12.5 °C) y la temperatura máxima (40.6 °C). Las líneas punteadas corresponden a la amplitud térmica de desempeño 80 (B₈₀ 24.6 - 36.2 °C), el rectángulo gris son los cuartiles 25% - 75% de las temperaturas seleccionadas (T_sel 26.4 - 34 °C) de la Estación Científica, y la línea amarilla es la temperatura óptima (T_opt 30.4 °C). 

Encontramos un ΔTP de 3.4 °C durante el fenómeno de El Niño en la estación fría, mientras que para la estación cálida fue de 1.8 °C. Por otro lado, utilizamos una regresión potencial para estimar el aumento en las T_o para la Estación Científica (T_o Niño = 0.2207*((T_air + ΔTP) ^ 1.5114)) y una regresión lineal para Baltra (T_oNiño = -38.494527 + (2.835572*(Tair + ΔTP)); tabla 1). Encontramos que M. indefatigabilis no presentó Hr en condiciones normales en la Estación Científica, ni cuando se presenta el fenómeno de El Niño. Sin embargo, en Baltra registramos 1 Hr promedio durante los días de muestreo, y estimamos 8.4 Hr promedio por día, durante la temporada fría de El Niño.

Discusión

En este estudio no encontramos diferencias significativas en la T_c de M. indefatigabilis entre jóvenes, hembras y machos, lo que concuerda con lo observado en Microlophus (Tropidurus) peruvianus y Microlophus bivittatus, las cuales tampoco presentaron diferencias significativas en la T_c entre categorías (^{Catenazzi et al., 2005}; ^{Huey, 1974b}; ^{Rowe et al., 2020}).

El clima característico de las islas Galápagos presenta 2 temporadas, la estación cálida y lluviosa (de enero a mayo) y la estación fría, típicamente nublada y seca en las zonas bajas (de junio a diciembre; ^{Trueman y d’Ozouville, 2010}). Las lagartijas pueden presentar una Tc más baja durante las estaciones secas y frías, en comparación con las estaciones cálidas y húmedas (^{Rowe et al., 2020}). No obstante, nosotros encontramos que el promedio de la Tc de M. indefatigabilis en la estación fría es similar a la Tc promedio observada en estudios previos para esta misma especie en la estación cálida (entre 33.7 y 35 °C; ^{Carpenter, 1970}; ^{Stebbins et al., 1967}). La Tc de M. indefatigabilis también es cercana a la de M. peruvianus (36.3, 30.6 y 33.0 °C en la estación fría y 36.1 °C en la cálida; ^{Catenazzi et al., 2005}; ^{Huey, 1974b}), pero es menor a la registrada para M. bivittatus en ambas estaciones (38.7°C en estación cálida y 38 °C en la estación fría; ^{Rowe et al., 2020}).

Encontramos que en la estación fría la actividad de M. indefatigabilis es unimodal y que el intervalo de actividad sucede entre 22.3 y 41.1 °C de T_c. Resultados similares fueron encontrados en una población de Microlophus delanonis en la isla Española (^{Altamirano, 1996}). Anteriormente, había sido considerado que el intervalo normal de actividad en M. indefatigabilis era bimodal entre 31 y 38 °C (^{Stebbins et al., 1967}). No obstante, en investigaciones del género Microlophus sobresale la importancia de la cobertura de las nubes en la fluctuación de la T_c. En el caso de M. indefatigabilis anteriormente se consideró que las T_c por debajo de 29 °C correspondían a organismos recolectados en la mañana o en días con nubosidad típica en la zona (^{Huey, 1974b}; ^{Rowe et al., 2020}; Stebbins et al., 1967). En M. bivittatus existen variaciones en las Tc de los organismos entre la hora del día, posiblemente ligadas a la presencia de nubes (Rowe et al., 2020). En consecuencia, el patrón de actividad de M. bivittatus podría ser diferente entre temporadas, años o incluso entre poblaciones en la misma temporada. Tal es el caso de M. peruvianus, que durante el invierno presenta un patrón unimodal en una población, mientras que en la otra estación es bimodal, lo que apoya la hipótesis de que los organismos reducen su actividad diaria y en consecuencia disminuyen su Tc en condiciones de mayor nubosidad (^{Catenazzi et al., 2005}). Por otra parte, también existe evidencia de la influencia de la nubosidad en los patrones de actividad diaria en las serpientes de Galápagos (Pseudalsophis), las cuales son depredadores naturales de las lagartijas de lava (Altamirano, 1996). Nuestra investigación se realizó a principios del mes de agosto, debido a esto, el factor nubosidad podría estar influenciando los datos obtenidos, por lo que consideramos importante incluir dicho factor en estudios futuros.

Aunado a lo anterior, se pueden presentar diferencias en los patrones de actividad entre poblaciones de la misma especie. Nosotros encontramos diferencias significativas entre la T_c de la población de Baltra y la población de la Estación Científica. Como consecuencia de ésto, las T_c de Baltra fueron las más altas, mientras que las de la Estación Científica las más bajas (fig. 1). Estas diferencias se pueden explicar porque la Estación Científica mantiene un mayor grado de conservación en la vegetación, en comparación con Punta Estrada, más que por la prevalencia de actividad humana, evidente en ambos sitios. A su vez, M. peruvianus también presenta diferencias en la T_c entre poblaciones, aun en la misma temporada (^{Catenazzi et al., 2005}).

Por otra parte, según el efecto Bogert la T_sel se conserva entre especies relacionadas, o entre poblaciones de diferentes ambientes, debido a que el comportamiento termorregulador evita la modificación fisiológica de los parámetros térmicos (^{Bogert, 1949}; ^{Catenazzi et al., 2005}; ^{Huey y Slatkin, 1976}). Nosotros encontramos que la T_sel fue similar entre las poblaciones estudiadas. No obstante, dicha T_sel promedio de M. indefatigabilis fue menor que la temperatura documentada previamente para la misma especie y para otras especies del género como M. peruvianus (35.6 °C) y M. bivittatus (36.1 - 40°C; ^{Catenazzi et al., 2005}; ^{Huey, 1974b}; ^{Rowe et al., 2020}). Este es un estudio prospectivo y, desafortunadamente, no tuvimos la oportunidad de estudiar la T_sel de las 3 poblaciones, no obstante, los datos recabados son importantes.

Nuestros datos muestran una marcada correlación entre la T_c de M. indefatigabilis y las temperaturas ambientales (T_a y T_s; fig. 1). No obstante, observamos que existe una mayor correlación con la T_s. Esto concuerda con M. peruvianus, la cual no mostró correlación significativa entre su T_c y la Ta en un primer estudio (^{Huey, 1974b}). Sin embargo, un estudio más reciente encontró una correlación débil con la T_a y más elevada con respecto del sustrato (^{Catenazzi et al., 2005}).

Encontramos que M. indefatigabilis tiende al extremo termorregulador en la temporada fría. Este comportamiento es similar al presentado por M. peruvianus (E = 0.74 y E = 0.86) durante el invierno (^{Catenazzi et al., 2005}). Aunado a ésto, nuestros resultados sugieren que M. indefatigabilis puede esforzarse y aumentar 2.7 °C su T_c con respecto a la del ambiente en la estación fría, la cual estuvo dentro de los intervalos observados para M. bivittatus en ambas estaciones (^{Rowe et al., 2020}).

En un estudio de M. bivittatus, los autores observaron que en promedio, la T_o de la estación fría fue menor al intervalo de T_sel , y ésta también varió entre micrositios. Por ejemplo, la T_o tomada en paredes sin cobertura vegetal se encontró dentro de la T_sel en breves periodos de actividad; sin embargo, en gran parte de las horas de la tarde, la T_o se encontró por encima del intervalo de T_sel (^{Rowe et al., 2020}). En este estudio encontramos que para M. indefatigabilis, el intervalo de T_sel y el promedio de la Tc están por encima del promedio de la T_o en la población de la Estación Científica, la cual presenta abundante cobertura vegetal. En contraparte, la T_sel al igual que el promedio de la T_c, se encuentran dentro del promedio de la T_o en Baltra, población que presenta escasa cobertura vegetal (fig. 2). Por consiguiente, consideramos relevante incluir la cobertura vegetal en estudios futuros debido a su implicación directa en la calidad térmica del hábitat (^{Neel y McBrayer, 2018}).

Anteriormente, se documentó que los organismos de M. indefatigabilis se encuentran activos a fines del verano desde el amanecer hasta la puesta del sol, con excepción del mediodía (^{Stebbins et al., 1967}). Como ya mencionamos, nosotros no observamos actividad bimodal a principios de la temporada fría. Sin embargo, ^{Stebbins et al. (1967)} comunicaron que las lagartijas que fueron observadas entre las 13:00 y 15:00 h en arbustos, en cactus y en los edificios fueron consideradas como inactivas. Contrariamente, nosotros consideramos la presencia de organismos en dichos sutratos en el horario mencionado como comportamiento activo. Anteriormente se había considerado que M. indefatigabilis puede presentar un adecuado desempeño locomotor a temperaturas inferiores de 30 °C, pero que solo se encontraba en su plena capacidad física entre los 31 y 38 °C (^{Stebbin et al., 1967}); no obstante, nosotros encontramos que la máxima capacidad locomotora se presenta a los 30.4 °C y la B₈₀ sugiere que los organismos pueden presentar un adecuado desempeño locomotor desde los 24.6 °C.

El cálculo de tolerancia al calentamiento sugiere que M. indefatigabilis presenta tolerancia adecuada ante el promedio de temperaturas ambientales disponibles. No obstante, M. indefatigabilis presenta un reducido margen de seguridad térmica, por lo que su desempeño podría disminuir, incluso ante ligeros aumentos de la temperatura ambiental. Esto concuerda con lo pronosticado en insectos que habitan regiones tropicales, las cuales se acercarán a temperaturas casi letales de forma veloz, debido a que dichas especies habitan en ambientes cercanos a sus temperaturas óptimas (^{Deutsch et al., 2008}). Nuestros resultados son similares a otras especies de lagartijas tropicales, como Anolis lemurinus (margen de seguridad térmica = 1.5°C), donde se propone que el tiempo disponible para la actividad cercana al rendimiento máximo, se reducirá ante el aumento de la temperatura ambiental (^{Logan at al., 2013}), así como A. cristatellus y varias especies de geckos diurnos (Sphaerodactylus), en las cuales se calcula que dichos organismos ya estén experimentando Tc estresantes durante verano, iguales o superiores a sus óptimos fisiológicos, y los cálculos sugieren que el calentamiento global deprimirá aún más su rendimiento fisiológico en dicha estación (^{Huey et al., 2009}).

Encontramos que la CT_Max de M. indefatigabilis se encuentra muy cercana del límite superior de la B₈₀. Existe escasa información respecto a estos parámetros en las especies del género Microlophus, pero este comportamiento ha sido observado en otras especies como Phrymaturus tenebrosus (^{Bonino et al., 2015}; ^{Cabezas-Cartes et al., 2019}), la cual adicionalmente presenta su rendimiento locomotor máximo en temperaturas similares a su T_sel, similar a lo que observamos en este estudio en M. indefatigabilis. No obstante, se sabe que existen diferencias en la resistencia locomotora entre sexos de M. indefatigabilis (^{Miles et al., 2001}; ^{Snell et al., 1988}), por lo que consideramos necesario abundar más en este tema en futuros estudios, ya que en casos como el de la lagartija Xantusia vigilis, a pesar de existir la capacidad de modificar su CT_Max y su umbral térmico, se pueden presentar efectos esterilizantes y atrofia testicular permanente en organismos que son sometidos a temperaturas incluso inferiores a la CT_Max (^{Cowles y Burleson, 1945}; ^{Kaufmann y Bennett, 1989}). Otro ejemplo es Sceloporus mucronatus, en esta especie, para que exista la maduración y motilidad espermática, se requiere que la T_c sea más baja que la T_sel. Así que exceder el umbral de tolerancia térmica máxima no solo puede ser perjudicial para la reproducción, también se requiere de temperaturas por debajo de la T_sel para que ocurra la reproducción (^{Méndez-de la Cruz et al., 2014}).

Existen 4 regiones de monitoreo a lo largo del océano Pacífico ecuatorial para el evento de El Niño. Las islas Galápagos, junto con las costas de Perú, se encuentran comprendidas en la región 1+2. En dicha región, para el mes de agosto, el evento de El Niño registrado como más fuerte fue el correspondiente a 1997, seguido del de 1972 y posteriormente el de 2015 (^{Yong, 2019}). Durante El Niño de 1997, se registraron 3,407.6 mm de lluvia en la Estación Científica Charles Darwin, lo cual convirtió a este periodo en uno de los más lluviosos en la historia de Galápagos (^{Snell y Rea, 1999}; ^{Yong, 2019}). Calculamos que durante este fenómeno las poblaciones de M. indefatigabilis que se ubican en sitios con poca cobertura vegetal, como es el caso de la población de Baltra, podrían presentar una mayor cantidad de Hr, en comparación con poblaciones que presentan mayor cobertura vegetal, como es el caso de la población de la Estación Científica. Por otra parte, en estudios que registran la abundancia de lagartijas, se ha observado que la disponibilidad de alimento y el número de individuos decrecen considerablemente en años donde se presenta el fenómeno de El Niño (^{Rodríguez-Romero y Méndez-de La Cruz, 2004}; ^{Stapley et al., 2015}; ^{Wikelski y Thom, 2000}), por lo que consideramos sumamente importante profundizar más en este tema.

A pesar de que ^{Sinervo et al. (2010)}, calcularon las Hr con los límites de la Tsel, estas pueden depender de las características de cada especie (helioterma, estenoterma, euriterma, etc.) y pueden variar entre las Tc (^{Pontes-da- Silva et al., 2018}), T_sel, o B₈₀, como en nuestro caso. Optamos por esta decisión metodológica debido a que fuera de la B₈₀, los organismos experimentarían una reducción sustancial en sus capacidades locomotoras y, en consecuencia, en su capacidad de permanecer activos. Además, la amplitud de desempeño se aproxima más a las temperaturas mínima y máxima voluntarias, las cuales representan más fielmente los umbrales térmicos de actividad, de manera que, en ausencia de estos 2 parámetros, la amplitud térmica constituye un mejor subrogado de las ventanas térmicas de actividad que la T_sel. En todo caso, al ser la B80 un intervalo más amplio que el de T_sel, se esperaría que las Hr obtenidas con base en este criterio fuesen más conservadoras que las estimaciones basadas en el intervalo de T_sel.

Por otra parte, en esta investigación únicamente calculamos las Hr durante la temporada fría, la cual concuerda con la época previa a la fase reproductora de la especie. Este periodo es crítico, ya que la disminución en la actividad puede provocar disminuciones en la densidad efectiva de la población para la época reproductora. El concepto de Hr está basado, principalmente, en la temporada en la que nacen las crías, ya que el aumento de la temperatura durante esta etapa hace más vulnerables a los organismos (^{Sinervo et al., 2010}). No obstante, podemos considerar que la temporada pre-reproductora es importante porque, durante ésta, los organismos acumulan la energía necesaria para prepararse para la reproducción (^{Guillette y Bearce, 1986}). En algunas especies como en S. mucronatus, previo a la reactivación gonadal se incrementan los cuerpos grasos de machos y hembras (^{Méndez-de la Cruz et al., 1993}). En el caso de las hembras se utilizan los cuerpos grasos para la producción de folículos vitelogénicos y, en caso de no tener energía suficiente, algunos folículos se vuelven atrésicos y se reduce el número de crías producidas (^{Méndez-de la Cruz et al., 1993}), llegando eventualmente a la ausencia total de crías en años muy severos (^{Rodríguez-Romero y Méndez- de La Cruz, 2004}). Mientras que en los machos se ha propuesto que dicha energía se utiliza principalmente para la defensa del territorio (^{Licht, 1974}). Por lo anterior debe considerarse que las Hr durante la fase pre-reproductora afectan el esfuerzo reproductor, incidiendo directamente en la capacidad de producción de crías y, eventualmente, en su colapso reproductor.

Aunado a lo anterior, encontramos en los registros de las temperaturas ambientales durante el fenómeno de El Niño, que el aumento de la temperatura ambiental es menor en la época reproductora (1.8 °C) en comparación con la temporada pre-reproductora (3.4 °C). Por lo que la presencia de Hr en esta especie, podría tener efectos perjudiciales en los organismos en la época pre-reproductora. No obstante, consideramos necesario replicar este tipo de análisis para ambas temporadas, e incluso obtener datos de T_o de diferentes años.

En conclusión, la lagartija de lava de Santa Cruz, M. indefatigabilis, presenta diferencias significativas entre las temperaturas corporales entre poblaciones, pero conservan sus temperaturas seleccionadas, y éstas, a su vez, son cercanas a su temperatura óptima. Su termorregulación es activa y el límite superior de su amplitud térmica se encuentra muy cercano a su temperatura crítica máxima. Ésto, aunado al incremento de las temperaturas ambientales presentes durante el fenómeno de El niño, puede provocar un aumento en el número de horas de restricción de la especie en sitios con poca cobertura vegetal, por lo menos en la temporada no reproductora. Por último, los datos sugieren que esta especie puede presentar actividad unimodal a temperaturas inferiores a las previamente mencionadas.