Introducción

El hábitat es un espacio específico que reúne condiciones bióticas y abióticas indispensables para la presencia y distribución de una especie (^{Delfín-Alfonso et al., 2014}), por lo cual, la selección del hábitat por parte de los individuos es crucial para asegurar su supervivencia en cualquiera de sus etapas de desarrollo (^{Ernst et al., 2012}). En especies ovíparas, los adultos son quienes seleccionan el sitio para la puesta de sus huevos y de esta manera el hábitat en donde se desarrollan las crías (^{Both et al., 2011}), algunas características extrínsecas que influyen en la selección del sitio para ovopositar son: 1) la cobertura del dosel, que influye en la regulación de la temperatura del agua; 2) el porcentaje de hojarasca utilizado como refugio y sitio de forrajeo; 3) las propiedades fisicoquímicas del agua como conductividad, pH y oxígeno disuelto, implicadas en el desarrollo de los individuos (^{Borges-Júnior y Rocha, 2013}; ^{Camacho-Rozo y Urbina-Cardona, 2021}; ^{Sah y Grafe, 2020}; ^{Thomas et al., 2019}); 4) el tipo de sustrato, que puede afectar las condiciones microambientales, como la estabilidad, la temperatura y la oxigenación del agua, esenciales para el desarrollo de los huevos y larvas (^{Wiens, 1972}; ^{Melo et al., 2018}). Un factor intrínseco que influye en la selección de los sitios de ovoposición es la duración de la metamorfosis, las especies con renacuajos que presentan un ciclo larvario largo seleccionan charcas permanentes, mientras que aquellas que tienen un periodo de metamorfosis corto utilizan charcas temporales (^{Borges-Júnior y Rocha, 2013}; ^{Peltzer y Lajmanovich, 2004}).

Los renacuajos son herbívoros especializados, se nutren a través de la filtración de algas adheridas a diferentes sustratos, de plancton y de la materia orgánica presentes en los sedimentos de los cuerpos de agua. De esta manera evitan la eutrofización e intervienen en la dinámica de los nutrientes en los hábitats acuáticos (^{Jacobson et al., 2019}; ^{Mohneke y Rödel, 2009}; ^{Ranvestel et al., 2004}). A pesar de la importancia ecológica de los renacuajos, la mayoría de los estudios realizados se enfocan en estudiar su ecología y describir la diversidad de anuros en sitios determinados (^{Leyte-Manrique et al., 2018}; ^{Martín-Regalado et al., 2016}; ^{Reyna-Bustos et al., 2022}; ^{Ríos-Rodas et al., 2020}). Para el caso de los renacuajos, las investigaciones incluyen descripciones taxonómicas (^{Canseco Márquez et al., 2003}; ^{Canseco-Márquez y Gutiérrez Mayén, 2010}; ^{Kaplan y Heimes, 2015}), estudios experimentales para conocer la supervivencia de los renacuajos expuestos a distintas concentraciones de salinidad (^{Woolrich-Piña et al., 2015}, ²⁰¹⁷), la tasa de malformaciones por exposición a contaminantes (^{Aguillón-Gutiérrez y Ramírez-Bautista, 2015}), pero se han dejado de lado las investigaciones acerca de la influencia de los gradientes ambientales sobre la selección de las charcas para la ovoposición y el desarrollo de renacuajos.

Los gradientes ambientales se refieren a las variaciones espaciales y temporales en las condiciones abióticas y bióticas que influyen en la estructura y dinámica de las comunidades acuáticas (^{Ouchi-de Melo et al., 2017}; ^{Schalk et al., 2017}). En este estudio, consideramos gradientes ambientales que incluyen las variables fisicoquímicas del agua (temperatura, salinidad y pH), la estructura de la vegetación circundante (vegetación primaria, secundaria y recién talada) y las características de las charcas (permanentes y temporales). Por lo cual, el objetivo de la presente investigación fue analizar la influencia de las variables bióticas y abióticas en los gradientes ambientales sobre la abundancia de 3 especies de renacuajos asociados a charcas permanentes y temporales en un arroyo tropical del sureste de México. Nuestra hipótesis indica que las variaciones en los factores bióticos y abióticos, así como los componentes del hábitat a lo largo de los gradientes ambientales, influyen significativamente en la abundancia relativa de las 3 especies de renacuajos en charcas permanentes y temporales del arroyo, y se registra una mayor abundancia de individuos en charcas permanentes debido a las condiciones más estables y favorables para su desarrollo.

Materiales y métodos

El trabajo de campo se realizó en el arroyo La Escalera en el ejido Villa de Guadalupe (17°21’38.23” N, 93°36’30.97” O), que forma parte del Complejo Ecoturístico Agua Selva, en el municipio de Huimanguillo, Tabasco, México (Fig. 1). El área se ubica en la zona montañosa de la región fisiográfica Sierra Norte de Chiapas, con 200 a 1,000 m snm (^{Alejandro-Montiel et al., 2010}; ^{Montalvo-Vargas y Castillo-Ramiro, 2018}). La temperatura media anual es de 20.6 °C, presenta un clima cálido húmedo (Af), con lluvias todo el año y una precipitación media anual de 3,638 mm, lo cual contribuye a la formación de más de 100 cascadas de aguas cristalinas (algunas permanentes y otras temporales), pozas de formación natural y arroyos permanentes continuos que son aprovechados por la comunidad local para realizar actividades cotidianas y ecoturísticas (^{Alejandro-Montiel et al., 2010}; ^{Carvajal-Hernández et al., 2018}; ^{Castillo-Acosta et al., 2019}). A lo largo del arroyo se distinguen 3 tipos de ambientes de acuerdo con el tipo y estado de la vegetación circundante: 1) vegetación primaria (VP; Fig. 2a), constituida por una selva alta perennifolia con árboles de hasta 45 m de altura como palo mulato (Bursera simaruba), caoba (Swietenia macrophylla), ramón (Brosimun alicastrum), ceiba (Ceiba petandra), zopo (Guatteria anomala), jobo (Spondias mombin) y una gran diversidad de epífitas como orquídeas, helechos y musgos (^{Palma-López et al., 2019}); 2) vegetación secundaria de selva alta perennifolia (VS; Fig. 2b), caracterizada por árboles de 10 a 20 m de altura como guarumo (Cecropia obtusifolia), ciruelillo (Trichilia havanensis) y naranjillo (Bernandia interrupta) (^{Carvajal-Hernández et al., 2018}; ^{Palma-López et al., 2011}; ^{Rodríguez y Banda, 2016}); y 3) vegetación recién talada (VRT; Fig. 2c), que se consideró como pequeños fragmentos (2.5 ha aproximadamente) de vegetación secundaria con arbustos menores a 4 m de altura inmersos en espacios abiertos generados por la tala, como resultado de la acción de las comunidades humanas para disponer de terrenos que les permitieran acceder a programas gubernamentales de desarrollo rural (^{Challenger y Soberano, 2018}). De acuerdo con los pobladores, la tala en este sitio ocurrió aproximadamente 2 meses antes de iniciar el trabajo de campo.

Mapa elaborado por Josué García León.

Figura 1 Área de estudio en la selva alta perennifolia de la Sierra Norte de Chiapas en el municipio de Huimanguillo, Tabasco, México. La figura superior muestra la ubicación del área de estudio dentro de México, la figura intermedia señala en el recuadro azul la ubicación del área de estudio dentro del estado de Tabasco y la figura inferior indica la distribución de los transectos de búsqueda de las charcas sobre el arroyo La Escalera. 

Figura 2 Vista de los tipos de ambientes presentes en el sitio de estudio en la Sierra Norte de Chiapas, municipio de Huimanguillo Tabasco, México. a) Vegetación primaria de selva alta perennifolia, b) vegetación secundaria, c) vegetación recién talada.  

Se realizaron muestreos diurnos de 10:00-15:00 h, iniciando la búsqueda en los transectos de menor altitud para evitar pseudoréplicas de individuos arrastrados por la corriente (^{Strauß et al., 2010}); el trabajo de campo se llevó a cabo de marzo a julio de 2021, con un esfuerzo de muestreo de 5 h, por 5 personas durante 2 días, obteniendo un total de 250 horas/hombres para todo el muestreo. Para el registro de los renacuajos se establecieron 5 transectos a lo largo del arroyo, 1 en VP, y 2 en VS y VRT. Cada transecto midió 200 m de largo con un ancho que varió de 10 a 15 m, todos los transectos estuvieron separados por 50 m para asegurar la independencia de los datos, con una separación mínima de 200 m entre cada tipo de vegetación. En cada transecto se localizaron charcas permanentes y temporales. Las permanentes estuvieron conectadas a la corriente del arroyo o con algún tipo de filtración superficial; mientras que las temporales fueron aquellas que se quedan sin agua en la temporada seca y se encuentran separadas del afluente principal (arroyo) por un mínimo de 2 m (^{Eterovick y Souza, 2003}). En cada charca se buscaron renacuajos que comprenden las etapas 24 a la 41, de acuerdo con la escala de ^{Gosner (1960)}, en ellas aún no hay desarrollo de extremidades y, por lo tanto, no es posible su movimiento entre las charcas (^{Eterovick y Souza, 2003}). Para determinar la abundancia relativa, se contabilizaron visualmente los renacuajos presentes en cada una de las charcas por 5 min.

La identidad taxonómica de los renacuajos se determinó in situ, utilizando las claves dicotómicas de ^{Köhler (2011)}, ^{Limbaugh y Volpe (1957)}, ^{Mijares-Urrita (1998)}, ^{Segura-Solís y Bolaños (2009)}, los individuos que no fue posible identificar de este modo, se colectaron bajo el permiso con oficio SGPA/DGVS/00962/22. La captura se realizó utilizando redes de inmersión de diferentes dimensiones (21.5 × 17.8 cm; 10 × 8 cm), se separaron por morfoespecie y se colocaron en contenedores de 500 ml, con agua de la charca. Los organismos se trasladaron al laboratorio de la Colección de Anfibios y Reptiles de Tabasco (CART), de la División Académica de Ciencias Biológicas, Universidad Juárez Autónoma de Tabasco para su posterior identificación.

Para definir los gradientes ambientales en el área de estudio, se consideró un conjunto de variables abióticas y bióticas en las charcas permanentes y temporales. Dentro de las variables abióticas se consideró la altitud, la cual se obtuvo con un sistema de geoposicionamiento global (GPS 64S, marca GARMIN); la temperatura, la conductividad, el oxígeno disuelto y la salinidad del agua se midieron con una sonda multiparamétrica Pro-230, marca YSI® Professional Series. Para medir el pH se utilizó un medidor de bolsillo marca pHep®5 de pH/temperatura con resolución 0.01, la profundidad de la hojarasca y de las charcas se midieron con una regla de 30 cm. En cada muestreo y en cada una de las charcas evaluadas se tomaron 3 medidas de las variables antes mencionadas, con el propósito de obtener un promedio para cada charca (^{Cruz-Ramírez et al., 2018}). Posteriormente, se identificó el tipo de sustrato presente en cada una de las charcas: arenoso, limoso, arcilloso (^{Eterovick y Souza, 2003}; ^{Melo et al., 2018}).

Las variables bióticas se calcularon en cada muestreo y en cada una de las charcas evaluadas, por ejemplo, la cobertura del dosel (CD), se midió con un densiómetro esférico cóncavo, para ello se realizaron 4 lecturas con dirección a los puntos cardinales a una altura de 1.30 m sobre el suelo, los valores se promediaron y, posteriormente, se multiplicaron por la constante 1.04 (^{Lemmon, 1956}), dando como resultado el porcentaje de la CD. Además, se consideraron como presentes o ausentes los siguientes componentes del hábitat: rocas, vegetación acuática, pastos, raíces, troncos y corriente del agua (^{Strauß et al., 2010}), así como los organismos acuáticos (peces, gasterópodos y decápodos).

Una vez en el laboratorio, los renacuajos se mantuvieron en acuaterrarios siguiendo las recomendaciones de ^{Poole y Grow (2012)}, con algunas modificaciones para asegurar la supervivencia de los organismos, como utilizar agua purificada con baja concentración de sales. Los acuaterrarios se oxigenaron con bombas de aire para peceras y plantas acuáticas, como sustrato se utilizó grava para acuario. Los renacuajos se alimentaron cada 48 h con cantidades pequeñas de espinacas previamente desinfectadas; después de 15 min se retiraron los excedentes de mayor tamaño para evitar enfermedades en los organismos (^{La Marca y Castellanos, 2018}). Una vez que los individuos completaron su metamorfosis, fueron identificados a nivel de especie con las claves antes mencionadas; posteriormente, fueron liberados en el sitio de captura. Para el caso exclusivo de Exerodonta bivocata, se colectaron 5 individuos para conservarlos en la CART, se sacrificaron con etanol al 4% y, posteriormente, fueron fijados con formol tamponado (neutro) al 10% por 2 semanas, transcurrido ese tiempo fueron preservados en alcohol etílico al 30% (^{Cortez et al., 2006}).

Análisis de datos. El grado de asociación de cada especie con los tipos de vegetación y charcas, se analizó a través de un mapa de calor representado a partir del índice de asociación de Whittaker (^{Somerfield y Clarke, 2013}). El mapa de calor es una herramienta gráfica que se utiliza para representar, en 2 dimensiones, una matriz de datos mediante un gradiente de colores. Cada color representa el valor numérico del índice de Whittaker, con colores que varían desde tonos fríos (como el azul), para indicar valores bajos de asociación de la abundancia de una especie por una charca, hasta tonos cálidos (como el rojo) para valores altos; lo que permite identificar visualmente patrones, agrupaciones o tendencias dentro de los datos, a partir de las abundancias de las especies en cada charca; adicionalmente, se realiza una clasificación, la cual se visualiza como un clúster de agrupamiento entre estos cuerpos de agua. Posteriormente, se aplicó un análisis Permanova para determinar el efecto del tipo de vegetación sobre las variables ambientales (factor fijo con 3 niveles: vegetación primaria, vegetación secundaria y vegetación recién talada), el tipo de charca (factor fijo con 2 niveles: permanente y temporal), la altitud y sus interacciones. Las variables ambientales fueron previamente estandarizadas y convertidas en una matriz de distancias euclidianas. El Permanova se realizó bajo una suma parcial de cuadrados (tipo III) y 9,999 permutaciones de los residuales bajo un modelo reducido (^{Anderson y Ter-Braak, 2003}). Los valores significativos y sus interacciones se analizaron mediante una comparación a posteriori con el estadístico “t de student” y 9,999 permutaciones.

Para ver gráficamente los resultados del Permanova sobre el gradiente ambiental, se realizó un análisis de coordenadas principales (PCoA), donde se representa la asociación de las variables, solo se visualizaron las variables con valores de correlación de Pearson mayores de 0.5 con alguno de los 2 primeros ejes de ordenación. Adicionalmente, se realizó un gráfico de burbujas sobre el PCoA para representar gráfica y descriptivamente la relación entre las variables ambientales con la abundancia total de renacuajos, cuanto más grande fuera la burbuja, mayor sería el número de renacuajos (^{Clarke y Gorley, 2015}).

Las diferencias entre la abundancia total de los renacuajos y por especie se evaluaron utilizando una matriz de abundancia la cual se transformó a raíz cuadrada y posteriormente se generó una matriz de distancias euclidianas; el diseño del Permanova siguió los parámetros descritos anteriormente para las variables ambientales. Finalmente, se realizó un análisis de correlación de Pearson con el fin de determinar la relación lineal entre las variables ambientales y descartar aquellas variables colineales con un coeficiente mayor a 0.7 para análisis posteriores. A fin de seleccionar los modelos que mejor explican la abundancia total de los renacuajos y de cada especie, en función de las variables ambientales, se usó un modelo lineal basando en distancias DistLM (^{Anderson et al., 2008}). El DistLM, es una técnica estadística que permite analizar las relaciones entre una matriz de similitud como variable de respuesta (que en este caso son las abundancias de cada especie convertida en una matriz de distancia euclidiana entre charcas) y un conjunto de variables ambientales predictoras, previamente estandarizadas (que se agrupan a partir de una matriz de distancias euclidianas entre charcas). Se basa en descomponer la varianza de esta matriz de distancias euclidianas de las abundancias (como variable de respuesta) con respecto a las distancias euclidianas de las variables predictoras, permitiendo así identificar qué variables predictoras tienen una influencia sobre la abundancia total y de cada especie de anuro en etapa larval dado el mejor modelo ajustado. Adicionalmente, se empleó el procedimiento de selección de BEST que permite identificar las variables ambientales que mejor explican cambios en la variable de respuesta (la abundancia total y de cada especie) y de esta forma se seleccionaron los modelos con mejor ajuste a partir del criterio de información de Akaike para muestras pequeñas (AICc; ^{McArdle y Anderson, 2001}). Todos los análisis se realizaron utilizando el software PRIMER v 7.0.21 y Permanova add -on v 1.04 (^{Anderson et al., 2008}). El uso del AICc enfoca la calidad de cada modelo relativo a los otros dentro de un conjunto de modelos candidatos, permitiendo una comparación basada en la parsimonia y la explicación de la variabilidad de la variable de respuesta, en lugar de usar pruebas de hipótesis individuales; lo que evita la dependencia en valores de p en los análisis de regresión y promueve un enfoque más integral en la selección de modelos (^{Anderson, 2007}; ^{Burnham y Anderson, 2004}; ^{Goodenough et al., 2012}).

Resultados

Se registraron 975 renacuajos de 3 especies: Lithobates vaillanti, Incilius valliceps y Exerodonta bivocata. Lithobates vaillanti presentó la mayor abundancia con 60.51% de las observaciones, seguida de I. valliceps con 28.41% y E. bivocata con 11.08%. Se contabilizaron un total de 55 charcas, en 40 de ellas se encontraron renacuajos, 25 fueron charcas permanentes y 15 temporales. De acuerdo con el mapa de calor en cada tipo de charca y según el tipo de vegetación, se determinó que no están presentes 2 especies en la misma charca. Exerodonta bivocata fue asociada a charcas permanentes ubicadas en VP, I. valliceps a charcas temporales en VS y L. vaillanti a charcas temporales y permanentes en VS y VRT (Fig. 3).

Figura 3 Mapa de calor que representa la asociación de las especies con los diferentes puntos de muestreo en el arroyo La Escalera en la selva alta perennifolia de la Sierra Norte de Chiapas en el municipio de Huimanguillo, Tabasco, México. El color de los pixeles representa el grado de asociación de cada especie con el estado de la vegetación y las charcas permanentes y temporales, donde los colores más cálidos representan el mayor grado de asociación. Las figuras geométricas representan el estado de la vegetación: vegetación secundaria; vegetación recién talada; vegetación primaria. Las figuras con espacios internos sin color indican charcas temporales, mientras que las figuras rellenas de color, charcas permanentes. 

El Permanova determinó que el tipo de vegetación tuvo un efecto significativo sobre las variables medidas en cada gradiente ambiental (pseudo-F = 3.46, p [perm] = 0.0015), dicho factor presentó el mayor porcentaje estimado de variación entre todos los factores evaluados (Tabla 1). El PCoA muestra en el primer eje que las charcas permanentes en VP presentaron mayor profundidad, salinidad, conductividad y corriente del agua, mientras que las charcas en VS y VRT presentaron mayores valores de profundidad de hojarasca, y presencia de hojarasca y raíces. En el segundo eje, se puede observar que las charcas permanentes y temporales en VRT presentaron un sustrato más limoso, además de gasterópodos, lo que las hace totalmente diferentes a las charcas permanentes y temporales ubicadas en VS y VP (Fig. 4a). La relación entre las variables ambientales con la abundancia total de renacuajos fue variada; sin embargo, el gráfico de burbuja muestra que cuando las charcas permanentes y temporales presentaron alta salinidad, conductividad y profundidad, se registró un número bajo de renacuajos (Fig. 4b).

Figura 4 a) Ordenación de las charcas muestreadas en el arroyo La Escalera en función de las variables ambientales medidas en la selva alta perennifolia de la Sierra Norte de Chiapas en el municipio de Huimanguillo, Tabasco, México. Las figuras geométricas representan el tipo de la vegetación: vegetación secundaria; vegetación recién talada; vegetación primaria. Las figuras con espacios internos sin color indican charcas temporales, mientras que las figuras rellenas de color, charcas permanentes. b) Ordenación de las charcas muestreadas en el arroyo La Escalera en función de las variables ambientales medidas en la selva alta perennifolia de la Sierra Norte de Chiapas en el municipio de Huimanguillo, Tabasco, México y la abundancia total de los renacuajos. La ubicación de cada círculo representa la distancia euclidiana entre charcas dados los gradientes ambientales y el tamaño de los círculos representa la abundancia total de las 3 especies de renacuajos en cada charca, en donde círculos de menor tamaño representan alrededor de 30 individuos, mientras que el círculo de mayor tamaño representa aproximadamente 300 individuos. 

Con respecto a la abundancia total de los renacuajos, el Permanova no presentó diferencias entre el estado de la vegetación, el tipo de charca, la altitud o sus interacciones (Tabla 2). No obstante, el componente estimado de variación muestra un tamaño del efecto elevado para los residuales (42%), la interacción entre el tipo de la vegetación y la altitud (tamaño del efecto = 19.1%), y la interacción entre el estado de la vegetación con el tipo de charca y la altitud (tamaño del efecto = 19%). El análisis de correlación de Pearson presentó una alta relación lineal entre suelo arenoso y limoso (-0.76), salinidad y conductividad (0.81) y oxígeno disuelto con porcentaje de saturación de oxígeno (0.84), estas variables fueron descartadas en los análisis posteriores. De acuerdo con el AICc, el modelo que mejor explica la abundancia total de los renacuajos en el arroyo incluyó un total de 8 variables predictoras, sin embargo, solo la temperatura del agua, la salinidad y pH son las variables con mayor peso y poder predictivo (AICc = 335.52, R^2 = 0.68, RSS = 17,911) (Tabla 3; material suplementario).

El resultado del Permanova para L. vaillanti sugiere que no existen diferencias estadísticamente significativas entre la abundancia y el tipo de la vegetación, el tipo de charca, la altitud o sus interacciones (Tabla 4), pero el componente estimado de variación indica un tamaño del efecto alto a nivel de los residuales (55%), la interacción entre los tipos de vegetación y la altitud (tamaño del efecto = 12.5%) y el estado de la vegetación (tamaño del efecto = 11.9%) (Tabla 4). La abundancia de L. vaillanti fue explicada por 7 variables significativas al modelo y las que tuvieron un mayor poder explicativo, en orden descendente, son: la ausencia de corriente en las charcas, la profundidad de la hojarasca, la temperatura del agua de 24ºC y el sustrato limoso (AICc = 328.33, R^2 = 0.70, RSS = 16,542; Tabla 3).

La abundancia de I. valliceps varió en función de la interacción entre el tipo de charca y la altitud (pseudo-F = 6.78, p [perm] = 0.0059, tamaño del efecto = 22.3%). También se presentó un componente de variación elevado en la interacción entre el estado de la vegetación, el tipo de charca y la altitud (tamaño del efecto = 26.3%) (Tabla 4). La abundancia de I. valliceps fue explicada por un total de 6 variables las cuales fueron: la profundidad de la hojarasca, el sustrato limoso, la presencia de ramas y gasterópodos, la ausencia de corriente en charcas temporales, y el pH de 7.2; pero ninguna de estas variables fue significativa al modelo (AICc = 322.9, R^2 = 0.39, RSS = 15,745; Tabla 3).

Por último, el bajo tamaño de la muestra para E. bivocata no permitió analizar la varianza. Sin embargo, se pudo evidenciar que el mayor componente de variación se encontró en los residuales (tamaño del efecto = 37%), seguido por la interacción entre el estado de la vegetación, el tipo de charca, la altitud (tamaño del efecto = 18%) y la interacción entre el tipo de vegetación y la altitud (tamaño del efecto = 12.4%) (Tabla 4). La abundancia de E. bivocata estuvo determinada por 5 variables en el mejor modelo ajustado, las que fueron significativas y presentaron mayor poder predictivo son, en orden descendente: la presencia de decápodos, una altitud de 460 m snm, la presencia de corriente en las charcas permanentes, el pH de 7.2 y la temperatura del agua de 22 ºC (AICc = 328.33, R^2 = 0.70, RSS = 16,542; Tabla 3).

Discusión

Los resultados de la presente investigación muestran que los renacuajos de L. vaillanti fueron más abundantes con respecto a E. bivocata e I. valliceps; ésto se debe a que la abundancia de los renacuajos está fuertemente influida por el tamaño de la puesta (^{Vogt, 1997}). En el caso de L. vaillanti, se compone de alrededor de 850 huevos (^{Hernández-Guzmán e Islas-Jesús, 2014}), mientras que para I. valliceps, el tamaño de puesta es aproximadamente de 100 huevos (^{Oliver-López et al., 2009}), para E. bivocata se desconoce el número de huevos por puesta; sin embargo, para otras especies de la familia Hylidae se ha reportado que las puestas están conformadas de 9 a 80 huevos (^{Leonard et al., 1993}; ^{Luja-Molina, 2010}; ^{Luja-Molina y Rodríguez-Estrella, 2016}; ^{Stebbins, 1951}). Aunado a lo anterior, las larvas de anfibios son más numerosas que los adultos y por estar confinados en las charcas, el conteo de los individuos es más preciso (^{Camacho-Rozo y Urbina-Cardona, 2021}).

La selección de las charcas por parte de los anuros adultos depende de su distribución espacial y temporal, el esfuerzo reproductivo y de diversos factores bióticos y abióticos que inciden en la dinámica de las poblaciones de renacuajos (^{Alford, 1999}; ^{Duellman y Trueb, 1986}; ^{Eterovick y Souza, 2003}; ^{Stebbins y Cohen, 1995}). De acuerdo con nuestros resultados, la selección de las charcas permanentes y temporales difiere entre las especies estudiadas, en el caso de los renacuajos de L. vaillanti fueron registrados generalmente en charcas permanentes dentro de la VS y VRT; se ha observado que las especies que seleccionan este tipo de charcas es debido a que los renacuajos presentan periodos de desarrollo más largos (^{Dodd, 1992}; ^{Skelly, 1997}) y que en el caso de la especie mencionada, se ha registrado un periodo de desarrollo de hasta 5 meses (^{Hernández-Guzmán e Islas-Jesús, 2014}). Además, los adultos de esta especie están fuertemente asociados a charcas permanentes que se forman en los arroyos de la vegetación secundaria y de los potreros adyacentes a la selva (^{Hernández-Ordóñez et al., 2015}; ^{Ríos-Rodas et al., 2020}; ^{Urbina-Cardona et al., 2006}), las cuales son utilizadas como sitio de forrajeo, reproducción y oviposición (^{Lee, 1996}). Diversos estudios han observado que algunas especies de ranas utilizan tanto charcas permanentes como temporales para llevar a cabo su metamorfosis (^{Eterovick y Barata, 2006}; ^{Melo et al., 2018}; ^{Peltzer y Lajmanovich, 2004}; ^{Rodrigues et al., 2010}), lo cual concuerda con lo observado en nuestro estudio, donde L. vaillanti utilizó ambos tipos de charcas. Incilius valliceps se observó solo en charcas temporales ubicadas en VS, el uso de estas charcas está relacionado con su rápido desarrollo larvario, el cual tiene una duración de 28 días (^{Aguilar-López, 2018}; ^{Volpe, 1957}). Se conoce poco sobre la ecología y modo reproductivo de E. bivocata, sin embargo, los resultados muestran que los adultos eligen únicamente charcas permanentes en VP para la oviposición de su progenie durante los meses de mayo y junio; de acuerdo con nuestras observaciones en campo, el desarrollo embrionario de esta especie tiene una duración aproximadamente de 5 semanas. Además, se han registrado renacuajos en noviembre, lo que sugiere que esta especie puede tener más de un pico explosivo de reproducción durante las épocas de seca y lluvia (Comentario personal Ríos-Rodas). Dado estas condiciones, el seleccionar charcas permanentes les permite aumentar la probabilidad de que los renacuajos completen su metamorfosis, ya que presentan menores riesgos de desecación, incluso durante la época de seca (^{Alford, 1999}; ^{Torres-Orozco et al., 2002}).

El tipo de vegetación y la temperatura son las variables que caracterizan las charcas en las que fueron registrados los renacuajos de nuestro sitio de estudio. Estas variables, están íntimamente correlacionadas, ya que la variación de la temperatura en las charcas puede ser generada por los cambios en la cobertura de la vegetación circundante (^{Granados-Sánchez et al., 2006}). De acuerdo con nuestros resultados y lo documentado en otros estudios, existen variables que condicionan el uso de las charcas como son: altos valores de la corriente, la profundidad del agua, el oxígeno disuelto, la conductividad y la salinidad, factores que intervienen en la selección de sitios de puestas, desarrollo de huevos y larvas (^{Borges-Júnior y Rocha, 2013}; ^{Eterovick y Sousa, 2003}; ^{Kopp y Eterovick, 2006}; ^{Strauß et al., 2010}; ^{Thomas et al., 2019})

Cada una de las especies de estudio ocupó charcas con características específicas que influyen en el desarrollo y supervivencia de los individuos. Los renacuajos de L. vaillanti utilizaron charcas permanentes y sin corriente durante todos los meses de muestreo, este tipo de charcas mantienen un nivel de agua adecuado aún en época de seca, permitiendo que la abundancia de los renacuajos sea constante a lo largo del tiempo. Incluso se ha observado que algunos anuros pueden ocupar estas charcas por varios meses para completar su desarrollo (^{Almeida-Gomes et al., 2012}; ^{Borgues-Júnior, 2007}; ^{Fatorelli et al., 2010}). La hojarasca es otra variable importante en las charcas donde se encuentra esta especie, ya que es utilizada como refugio ante la presencia de depredadores (^{Melo et al., 2018}), como son los peces del género Rhamdia, los cuales se encuentran en las mismas charcas que L. vaillanti.

Los renacuajos de I. valliceps se encontraron solo en abril en charcas temporales sin corriente ubicadas exclusivamente en VS. Los adultos de esta especie tienen preferencia por los hábitats que presentan algún tipo de perturbación antropogénica, en los cuales son más abundantes (^{Cedeño-Vázquez et al., 2006}). Además, la selección de las charcas temporales por parte de esta especie se debe a que los adultos son estrictamente terrestres y utilizan las charcas únicamente para la oviposición (^{Oliver-López et al., 2009}). Las charcas seleccionadas por esta especie en nuestro estudio se caracterizaron por un bajo porcentaje de hojarasca y ramas, un sustrato limoso y la ausencia de depredadores. Esta última característica permite un desarrollo larvario relativamente más rápido, ya que al no tener depredadores potenciales invierten más tiempo alimentándose (^{Aguilar-López. 2018}; ^{Eterovick y Sazima, 2000}; ^{Peltzer y Lajmanovich, 2004}). Además, los valores altos de pH registrados para las charcas son adecuados para su supervivencia, ya que se ha demostrado que un pH ácido menor a 4 afecta el desarrollo de los renacuajos (^{Freda y Dunson, 1985}; ^{Rosenberg y Pierce, 1995}; ^{Thabah et al., 2018}).

A diferencia de las otras 2 especies de estudio, los renacuajos de E. bivocata solo ocuparon charcas permanentes ubicadas en VP, ésto concuerda con lo reportado por ^{Ríos-Rodas et al. (2020)} y ^{Santos-Barrera (2004)}, quienes observaron a la especie en estado adulto en áreas de vegetación primaria, lo que indica que esta especie solo se reproduce en sitios conservados. ^{Alford (1999)} menciona que los cuerpos de agua permanentes aumentan las probabilidades de que las especies completen su proceso de metamorfosis reduciendo así, el porcentaje de mortalidad. Otra particularidad específica de las charcas en las que se registró E. bivocata es la presencia de corriente; dadas estas condiciones, los renacuajos se registraron en el fondo de la charca, donde la corriente es menor y se encuentran adheridos a las rocas para evitar el arrastre (^{Thomas et al., 2019}). Adicionalmente, E. bivocata necesita charcas que estén rodeadas de vegetación, lo cual permitirá que conserven una temperatura promedio de 22 ºC, en cambio las charcas con una cobertura de dosel abierto pueden ser de 2 a 5 ºC más cálidas, datos que coinciden con lo registrado por ^{Skelly et al. (2002)}, ^{Werner y Glennemeier (1999)}. Tomando en cuenta las características necesarias para el desarrollo de este hílido, se puede concluir que es una especie especialista en la selección de charcas para el desarrollo de su progenie.

La presente investigación es la primera que aborda la influencia de las variables ambientales sobre la abundancia de los renacuajos en las charcas presentes de un arroyo tropical en México. Las 3 especies registradas representaron solo una parte de la comunidad de anuros en la zona, las cuales presentan desarrollo indirecto; además, es importante señalar que cada una de las especies habitó exclusivamente un tipo de charca. De acuerdo con nuestros resultados, se acepta la hipótesis de que las charcas permanentes, en general, sí presentaron la mayor abundancia de renacuajos, sin embargo, las charcas permanentes en VP presentaron la menor abundancia de renacuajos, mientras que las charcas permanentes en VRT y VS obtuvieron la mayor abundancia. Las variables que influyeron directamente en la abundancia de las 3 especies de renacuajos fueron el pH, la temperatura y el oxígeno disuelto en el agua, lo que concuerda con nuestra hipótesis propuesta. Para el caso particular de E. bivocata, la presencia de gasterópodos en las charcas influyó en su abundancia, estos organismos pueden ser considerados competidores o depredadores ejerciendo una presión constante en la abundancia de esta especie. Por último, este estudio también reveló la importancia del tipo de vegetación circundante a las charcas; las charcas en VP presentaron condiciones adecuadas para que especies exclusivas de ambientes conservados lleven a cabo su reproducción y desarrollo; mientras que las charcas, rodeadas de VRT fueron aprovechadas por especies generalistas. Dado que la selección de las charcas por parte de los anuros adultos depende del tipo de vegetación y que en la zona de estudio se implementan actividades que modifican la estructura de la vegetación, es necesario continuar con el monitoreo de las especies ante los cambios en la cobertura vegetal y analizar la relación entre los cambios en la configuración del paisaje y la dinámica poblacional de los anuros, en particular, de especies especialistas, de escaso conocimiento en sus modos reproductivos o que se encuentren en alguna categoría de riesgo.