Introducción

En los Montes de María (Caribe colombiano), la destrucción del bosque seco tropical (BST) avanza de forma acelerada, con la expansión de las actividades antrópicas como la ganadería, siembra de cultivos, plantaciones forestales y extracción de madera (^{Castellanos-Barliza et al., 2018}; ^{González et al., 2014}; ^{Pizano y García, 2014}; ^{Rangel y Martínez, 2017}). Lo anterior influye en la disponibilidad de hábitats para taxones de insectos como las mariposas (^{Mercado et al., 2018}; ^{Vargas-Zapata et al., 2012}). Paralelo a la modificación y reducción de los hábitats, se producen cambios en los atributos de las especies como la reducción en sus poblaciones, cambios en su estructura y alteraciones en sus patrones de vuelo y desplazamiento (Bell et al., 1991; ^{Márquez y Martínez, 2020}; ^{Townsend et al., 2001}). Adicionalmente, la composición y usos del suelo en un paisaje juegan un papel importante en el grado de disponibilidad y distribución de los recursos nutricionales, sitios de apareamiento y la dispersión de las mariposas, ya que estos factores influyen en la estructura y desplazamiento de estos insectos (^{Andrade et al., 2007}; Brown et al., 1991; Constantino et al., 2007; ^{Rueda, 2016}).

Por otro lado, los cambios en la estructura vegetal de los fragmentos y la pérdida de cobertura en el BST inciden en la dinámica espacio-temporal de la abundancia y en el desplazamiento de algunos grupos de mariposas (^{Vargas-Zapata et al., 2015}). Con el fin de aportar a esta discusión, es necesario utilizar taxones focales dentro de estos insectos como las mariposas del género Hamadryas Hübner (1806), debido a la relación estrecha que mantienen con los cambios de diversidad y abundancias de las especies vegetales y variables ambientales (^{Brown y Freitas, 2002}; ^{Kremen et al., 1993}). Así mismo, dichos cambios pueden afectar a las poblaciones adultas dentro del género, ya que éstos se alimentan de la savia de los árboles y de las frutas fermentadas; además, suelen utilizar estos sustratos para posarse, exhibirse, realizar búsqueda de pareja y defender su territorio y sitios de alimentación (^{Marini-Filho, 1996}). De este modo, la pérdida de hábitat es una amenaza para las Hamadryas en regiones como los Montes de María (MM) en el departamento de Bolívar, donde la reducción de los relictos de BST es constante y avanza de forma acelerada (^{Ahumada et al., 2019}; ^{Pizano y García, 2014}).

El género Hamadryas comprende un grupo de 20 especies de mariposas neotropicales (^{Lamas et al., 2004}; ^{Reverdin, 1914}), su tamaño varía entre los 40 mm en H. cloe hasta los 80 mm de envergadura en H. guatemalena (^{Monge-Nájera, 1992}; ^{Monge-Nájera y Hernández, 1991}). Para el Caribe colombiano se han reportado 4 especies de Hamadryas (^{Andrade et al., 2007}), la mayoría de estos registros proviene de inventarios sobre la composición y estructura de lepidópteros en distintos fragmentos de BST en el departamento del Atlántico (^{Boom-Urreta et al., 2013}; ^{Montero et al., 2009}; ^{Prince-Chacón et al., 2011}; ^{Vargas-Zapata et al., 2015}), Córdoba (^{Campos-Salazar et al., 2011}), Magdalena (^{Vargas-Zapata et al., 2011}) y Sucre (^{Mercado et al., 2018}). Sin embargo, el seguimiento de este grupo de insectos desde un enfoque poblacional y su desplazamiento entre áreas con vegetación boscosa y las matrices adyacentes (áreas de pastizales, cultivos anuales y permanentes y bordes) se desconoce; razón por la que se espera que las especies presenten variaciones en su abundancia y desplazamiento en los distintos usos de suelo dentro del paisaje.

Con la finalidad de aportar información acerca de la dinámica poblacional de este grupo de mariposas, se analizaron la abundancia, tamaño poblacional, proporciones sexuales, estructura de edades y el desplazamiento de las especies del género Hamadryas (Biblidinae) en diferentes usos del suelo en un paisaje fragmentado en los Montes de María, San Jacinto, Bolívar y su relación con las variables ambientales y la cobertura vegetal.

Materiales y métodos

El paisaje de BST se encuentra en la vereda La Flecha (VLF), en el sistema montañoso de los Montes de María, municipio de San Jacinto, Bolívar (09°51’12.4” N, 75°10’41.4” O) a una altitud que oscila entre 324 y 500 m (Fig. 1). El área presenta un clima tropical determinado por los vientos alisios del norte y noreste, los sistemas montañosos y el mar (^{Aguilera, 2013}), con temperaturas oscilantes entre 27 y 32 °C. Su precipitación se clasifica como bimodal tetraestacional, con un periodo de extensa sequía y lluvias regulares entre abril, julio, septiembre y noviembre, con un nivel de precipitación entre 1,000 y 1,900 mm al año. Además, presenta una humedad relativa promedio de 75% (^{Hernández, 2008}), presentando transición entre periodos secos y húmedos.

Figura 1 Localización de la vereda La Flecha (A) y estaciones de muestreo (B) en un paisaje fragmentado en los Montes de María, San Jacinto, Bolívar. 

La vegetación predominante en los fragmentos de BST es la higrotropofítica (Holdridge, 1978), con paisaje montañoso (USAID, 2016). El arbolado varía entre los 22 y 25 m de altura, con estratos vegetales de tipo arbustivo, sotobosque y dosel. Entre las especies vegetales que predominan se pueden encontrar Bursera simaruba (L.) (Sarg.) (indio encuero), Anacardium excelsum (Bertero ex Kunth) Skeels (caracolí), Aspidosperma polyneuron (Müll. Arg.) (carreto), Bocageopsis sp. (yaya), Caesalpinia coriaria (Jacq.) Wild. (dividivi), Pseudobombax septenatum (Jacq.) Dugand. (majagua), Crescentia cujete L. (totumo), Gustavia superba (Kunth) O. Berg. (membrillo), Bulnesia arborea (Jacq.) Engl. (guayacán), Hura crepitans L. (ceiba lechosa), Spondias mombin L. (jobo), Guazuma ulmifolia Lam. (guácimo), Chrysophyllum cainito L. (caimito), Quadrella adoratissima (Jacq.) Hutch. (olivo), Sterculia apetala (Jacq.) H. Karst. (camajón), Gliricidia sepium (Jacq.) Walp. (matarraton), Ficus sp. (caucho) entre otras, y a su vez se encuentran palmas como Bactris guineensis (L.) H.E. Moore, Attalea butyracea (Mutis ex L.F.) Wess. Boer (palma de vino) y Sabal mauritiiformis (H. Karst) Griseb. et H. Wendl. (palma amarga) entremezcladas en la zona de bosque. Alrededor del fragmento de bosque, se encuentran pastizales (gramíneas introducidas), pastizales arbolados (gramíneas, arboles maderables y palmas amargas), cercas vivas y la presencia de policultivos como Mangifera indica L. (mango), Dioscorea esculenta (Lour) Burkill (ñame), Zea mays L. (maíz), Capsicum annuum L. (ají), Persea americana Mill. (aguacate), Manihot esculenta Crantz. (yuca), Musa x paradisiaca L. (plátano), Ananas comosus L. Merr. (piña), entre otros.

En el área de estudio se delimitó una ventana de 2 x 2 km para un total de 4 km², donde se establecieron 7 estaciones (E1 - E7) de muestreo con un mínimo de 180 m de distancia y un máximo de 275 m con el fin de cubrir diversos usos de suelo. Entre los usos, se destacan los cultivos mixtos permanentes temporales (CMPT) o E1, compuesta por cultivos mixtos temporales (D. esculenta, Z. mays, C. annuum, M. esculenta, C. papaya) y permanentes de P. americana (aguacate), A. butyracea (palma de vino), M. indica (mango), C. cainito (caimito) y B. guineensis con algunos árboles típicos del bosque como A. excelsum (Bertero ex Kunth) Skeels (caracolí), H. crepitans (ceiba blanca), S. mombin (jobo) y A. polyneuron (carreto) con remanentes de vegetación adyacentes a quebradas temporales. Borde de bosque (BB) o E2, con presencia de cercas vivas constituida por especies predominantes como H. crepitans (ceiba lechosa), S. mombim (jobo), Quadrella adoratissima (olivo) y S. apetala (camajón) y de manera adyacente G. sepium (matarratón) y enredaderas. Las estaciones E3, E4 y E5 corresponden a zona de bosque (Bos); mientras que, E6 corresponde a una zona de vegetación secundaria (VS) con 12 - 15 años de abandono, con presencia de especies vegetales como H. crepitans (ceiba blanca), Bursera simaruba (indio encuero), Q. adoratissima (olivo); con alturas entre 3 y 15 m y DAP entre 3 (mín) y 15.2 cm² (máx), así como enredaderas y arbustos como G. sepium (matarratón). La E7 corresponde a un área de cultivos mixtos permanentes (CMP), donde predominan árboles frutales como P. americana (aguacate), M. indica (mango) con alturas de 3 - 10 m, Musa sp. (guineo cuatro filos), Annona muricata L. (guanábana), S. purpurea L. (ciruela), palmas como A. butyracea (palma de vino) entremezcladas con zonas de cultivos de M. x paradisiaca (plátano) y A. comosus (piña).

Se realizaron 6 muestreos separados por un promedio de 15 días, entre marzo y junio de 2018, de los cuales 2 se desarrollaron durante época seca (M1- marzo 17 y M2 - abril 3) y los restantes durante las lluvias (M3 - abril hasta M6 - junio 6). La captura de las mariposas se realizó con la ayuda de trampas Van Someren-Rydon (VSR), las cuales fueron cebadas con una mezcla de frutas (M. x paradisiaca, C. papaya, A. comosus, M. indica) y panela fermentada durante 48 h (^{Devries y Walla, 2001}; ^{Ribeiro y Freitas, 2011}). Las trampas fueron cebadas a las 7:00 am y revisadas en horas de la tarde del primer día y a las 24, 48 y 96 hrs, periodo en el cual se realizaba el reemplazo del atrayente, iniciando consecutivamente desde la E1 hasta la E7.

En cada una de las 7 estaciones se colocaron 5 trampas Van Someren-Rydon (VSR) (^{Devries, 1988}) distribuidas en formas de “zig-zag” y separadas por un promedio de 61.8 ± 9 m entre sí (mínimo 50 m, máximo 70 m), a una altura entre 1 y 1.2 m del suelo (^{Ribeiro y Freitas, 2011}; ^{Ribeiro et al., 2008}, ²⁰⁰⁹, ²⁰¹⁰, ²⁰¹²; ^{Sambhu et al., 2018}). Las trampas fueron georreferenciadas con un GPS (GARMIN - Monterra) con el fin de determinar la distancia de cada mariposa una vez que fuera recapturada, donde los datos de recapturas eran tomados en cuenta si los individuos eran determinados registrados nuevamente entre muestreos. La captura de los ejemplares se realizó de forma pasiva, utilizando la técnica de captura-marcaje-recaptura (CMR). Cada ejemplar capturado fue registrado con los datos de captura de cada punto con su respectivo sexo. Adicionalmente, se sacrificaron un ejemplar macho y uno hembra, representante de cada especie por la técnica de sujeción torácica propuesta por ^{Neild (1996)}; para la corroboración taxonómica se siguieron las claves propuestas por ^{Chacón y Montero (2007)} y la página ilustrada Butterflies of America (^{Warren et al., 2013}). Las mariposas capturadas fueron marcadas con un marcador permanente “Sharpie” en la celda discal del ala anterior con una combinación alfa-numérica, correspondiente al muestreo y un número consecutivo del marcaje de individuos (ej. Individuo 10A). Posteriormente estas fueron liberadas a su ambiente natural, con el objetivo de no contar nuevamente los individuos (^{Constantino, 2007}; ^{Ribeiro y Freitas, 2011}; ^{Ribeiro et al., 2008}, ²⁰¹⁰, ²⁰¹²) y así registrar los individuos recapturados con los datos de cada punto y de esta manera determinar su desplazamiento; el cual fue determinado como la distancia recorrida por los individuos desde los puntos donde era capturados hasta los puntos donde eran recapturados (^{Prieto et al., 2005}).

Adicional a lo anterior, en cada punto (hora de revisión por trampa en cada muestreo) se midieron la temperatura ambiente (°C) y la humedad relativa (%), con la ayuda de un termohigrómetro marca Extech-IR, la intensidad lumínica con un luxómetro marca Extech-DT-1308, y la cobertura del dosel con un densitómetro Marca Spherical Densiometer Model-A. Para esto último, se siguió la metodología propuesta por ^{Mostacedo y Fredericksen (2000)}, quien describe que el espejo del equipo debe colocarse de manera horizontal y contar el número de cuadros donde, por lo menos, la mitad de su área se encuentre cubierto por ramas, hojas o troncos, y a su vez, estos se multiplicaron por un factor de corrección de 1.04 para obtener la cobertura del dosel. De esta manera, en cada punto, la medición se tomó haciendo 4 lecturas, orientando el densitómetro hacia los puntos norte, oeste, sur y este de 5 m de distancia de la respectiva trampa en cada sitio.

Con los registros de los individuos capturados (C), marcados (M) y recapturados (R) o CMR por trampa, estación y muestreo, se determinó el número total de individuos capturados o abundancia total (Ntotal) y el número neto de individuos o abundancia neta (Nnet = Ntotal - número de individuos recapturados). Posteriormente, con los datos de CMR se realizó una matriz de Trellis sin tener en cuenta el sexo respectivo para cada especie (^{Krebs, 1989}; ^{Lemos et al., 2005}). De esta manera, esta matriz se empleó para calcular el tamaño poblacional de las especies H. feronia y H. februa en la VLF, mediante el estimador de Jolly-Seber (J.S.) (1964), descrito por ^{Krebs (1989)}, ^{Henderson (2003)} y ^{Lemos et al. (2005)}, empleando el programa de Ecological Methodology versión 5.1 (^{Krebs, 1998}). Los cálculos fueron realizados para estas 2 especies, ya que fueron las que registraron el mayor porcentaje de recaptura. La proporción de sexos se calculó de acuerdo a la siguiente formula: número de machos/número de machos + número de hembras; el valor resultante de la proporción sexual es 0.5 si el número de machos y hembras es igual (^{Wilson y Hardy, 2002}). Este método es el más recomendado para estimar las proporciones sexuales en la naturaleza, debido a que considera a los individuos como unidades discretas y por lo tanto refleja la abundancia relativa de cada sexo en una población (Ancona et al., 2017).

Los individuos adultos marcados se agruparon en 3 categorías de edad, las cuales se basaron en el estado de las alas (^{Freitas y Ramos, 2001}). Las categorías involucradas son: nuevos (individuos recién emergidos), intermedios (adultos con alas en estado medio) y viejos (incluye individuos viejos o alas muy desgastadas). De esta manera, la estructura de edades se estimó como la proporción de individuos de cada categoría en relación con el total de individuos marcados. Mediante los datos de CMR de cada especie (^{Freitas y Ramos, 2001}), se estimaron los desplazamientos totales para cada mariposa por especie y sexo (^{Freitas y Ramos, 2001}). A su vez, se registró la media o el valor de desplazamiento medio, la mediana y la distancia mínima o valor de desplazamiento más bajo y el tiempo promedio de desplazamiento, teniendo en cuenta el tiempo de residencia o tiempo en el que los individuos de una especie permanecieron en los distintos usos de suelo dentro del paisaje (^{Silva y Hernández, 2015}). De esta forma, se calculó la tasa media de desplazamiento (Tmm) como la media de desplazamiento en relación con el desplazamiento promedio de cada especie. Las distancias recorridas por individuo en línea recta entre áreas se calcularon utilizando el software ArGIS® Versión 10.5.

Por otro lado, se aplicaron los modelos lineales generalizados (MLG) para determinar el efecto de los usos y épocas (muestreos) como factores y la cobertura del dosel (CD) sobre la abundancia de las Hamadryas registradas. Los datos de abundancia (por especie) se modelaron siguiendo una distribución de Poisson (O’Hara y Kotze, 2010). En el modelo se incluyeron los efectos fijos, como el número de trampas, los factores como uso del suelo (estaciones) y épocas (muestreos: seca y lluvias), y la CD. La temperatura, la humedad ambiental y la intensidad lumínica se eliminaron del modelo debido a la alta correlación con la cobertura del dosel. El mejor modelo fue seleccionado siguiendo el criterio de información Akaike (AIC) más bajo, corregido para un tamaño de muestra pequeño (AICc; según complejidad, el número de parámetros y la aptitud de cada modelo) manteniéndose aquellos modelos con una disminución de 2 unidades AIC (20%) con respecto al modelo nulo. Este análisis se realizó utilizando la función glmer en el paquete lme4 (Bates et al., 2015) con el programa R Studio versión 4.1.2 de 64 bit (^{Balzarini et al., 2008}).

Resultados

Se capturaron y liberaron 541 individuos del género Hamadryas, de los cuales 312 fueron machos y 229 hembras. La especie más abundante fue H. feronia con 291 individuos marcados y liberados (164 machos: 127 hembras). La segunda especie más abundante fue H. februa con 137 individuos (79 machos: 58 hembras), mientras que H. amphichloe presentó un total de 80 individuos (47 machos: 33 hembras) y H. amphinome con 33 individuos (22 machos: 11 hembras) fue la especie menos abundante dentro del paisaje (Tabla 1). El mayor número de capturas se registró en marzo (M2) con 150 individuos (85 machos; 65 hembras). De manera general, el mayor número de machos registrados se presentó específicamente en la estación 2 (E2) durante el M2, al igual que el mayor número de hembras con 36 y 29 capturas, respectivamente (Tabla 1). El menor número de capturas se presentó durante inicios de junio (M6) con 17 individuos, el cual coincidió con la etapa de post lluvia, respectivamente (Tabla 1).

Se registraron un total de 71 recapturas (16.1%), de los cuales 41 fueron machos y 30 hembras. Del total de recapturas, 62 individuos (38 machos; 24 hembras) fueron recapturados una sola vez, 5 (3 machos; 2 hembras) se recapturaron 2 veces, mientras que 2 individuos (2 hembras) fueron recapturados 3 veces. La especie H. feronia presentó el mayor porcentaje de recapturas (20.96%) para un total de 11.68% de machos y 9.28% de hembras, presentándose el mayor valor (12 individuos) en la E7 (cultivos mixtos) para el M3 y M5 (8 machos - 4 hembras), respectivamente. Para H. februa, el registro de recapturas fue de 12.4% con un total de 7.3% machos y 5.1% para las hembras, y el mayor registro de individuos recapturados (4 individuos) se presentó en el borde (E2) durante el M2. De las 4 especies de Hamadryas capturadas en el área de estudio, H. amphichloe presentó el porcentaje más bajo de recapturas con 6.25% (3.75% machos; 2.5% hembras), y el valor máximo de recapturas (2 individuos) se presentó durante el M4 (1 macho - 1 hembra), en borde (E2) y cultivos mixtos (E7), mientras que para la especie H. amphinome todos los individuos recapturados fueron machos con 12.12% de la población y el menor valor de recapturas se presentó durante el M3 en las E5 (bosque) y E7.

El tamaño estimado de la población varió durante los eventos de captura y recaptura para cada una de las especies en los diferentes usos del suelo en la zona. A través del registro directo de individuos observados se determinó un menor tamaño poblacional para H. feronia en el M6 (inicio de junio), con 9 individuos y un mayor valor en el M2 (abril) con 93 individuos (Tabla 1, fig. 5). Para esta especie, mediante el estimador de Jolly-Seber (JS), el valor máximo estimado fue de 329.9 individuos en el M3 (mediados de abril) y el menor en el M5 (inicios de mayo) con 36 individuos (ICI = 18.218 y ICS = 134.535). Las probabilidades de supervivencia más altas se observaron en M1 (marzo) con un valor de 0.681 (Tabla 2). En el caso de la especie H. februa, el menor valor observado fue de 9 individuos en el M6 (junio) y el mayor valor en M1 (marzo) con 60 individuos. Con el estimador de Jolly-Seber (JS), el mayor tamaño poblacional para esta especie fue de 133 individuos en M2 (abril) y el menor en M5 (mayo) con 4 individuos (ICI = 2.479 e ICS = 47.902). Las probabilidades de supervivencia más alta se observaron en el M1 (marzo) con un valor de 0.77 y el más bajo se presentó en el muestreo 2 con 0.26 (Tabla 2).

La estructura de edades varió para cada una de las 4 especies durante los 6 eventos de muestreos. La clase joven (A) dominó en cada una de las especies de Hamadryas involucradas en este estudio durante los últimos muestreos, mientras las clases intermedias (B) y vieja (C) presentaron una mayor frecuencia durante los primeros muestreos para cada una de las especies. De esta manera, la estructura de edades no fue constante en ninguna de las especies de estudio, presentando fluctuaciones durante los eventos de captura dentro del paisaje en la VLF (Fig. 2).

Figura 2 Estructura de edades en las poblaciones de Hamadryas en los diferentes usos del suelo en la VLF entre los meses de marzo y junio del 2018. A), H. feronia; B), H. februa; C), H. amphichloe; D), H. amphinome. 

De acuerdo con los valores de ASR de las 4 especies estudiadas, se determina una predominancia de los machos para cada una de las especies de Hamadryas durante los eventos de muestreos. El solapamiento de los intervalos de confianza (IC de 95%) indica que no hubo diferencias significativas en los valores de ASR observados (Fig. 3). Para H. amphichloe, el rango de variación de los IC del ASR fue mayor durante la época de lluvias (M3, M4, M5, M6) con una clara dominancia de los machos. En cuanto a la especie H. amphinome, se determinó una alta dominancia de los machos desde el M2 (mediados de marzo), mientras que la especie H. februa presentó una tendencia hacía los machos en la época seca (M1, M2) y con respecto a la época de las primeras lluvias, se determinó una tendencia hacia las hembras (M3). En los últimos periodos (lluvias fuertes) se observó una inclinación hacia las hembras en los valores de ASR observados. Por otro lado, para la especie H. feronia se determinó una dominancia de los machos en la temporada seca y primeras precipitaciones, mientras que durante los últimos periodos de muestreo (lluvias fuertes) se determinó una alta frecuencia de las hembras influenciadas por las lluvias (M5, M6) (Fig. 3).

Figura 3 Variación temporal (A) y espacial (B) de ASR (adult sex ratio) de Hamadryas durante 6 muestreos (época seca y época lluviosa) en un paisaje fragmentado en la VLF. Los puntos representan la ASR observada para cada población y las barras son los intervalos de confianza de 95% estimados por bootstrap. Las líneas horizontales discontinuas a 0.5 representan la ASR en equilibrio (el mismo número de hembras y machos). 

Dentro de las distancias recorridas promedio estimadas, H. feronia fue la especie con el menor valor (127.825 m; IC = 43.15), sin embargo, presentó la tasa media de movimiento más alta (0.163 m/días). El individuo con mayor desplazamiento fue el registrado como 434D con 476 m entre la E6 hasta la E7 (Tabla 3). La especie H. februa presentó un valor promedio de desplazamiento de 208.846 m (IC = 99.63); el individuo con mayor distancia recorrida fue el 38A con 820 m, desde E5 (bosque) hasta E7 (cultivos mixtos), mientras que la tasa media de desplazamiento estimada para la especie fue de 0.095 m/días.

Para H. amphichloe se determinó un valor promedio desplazado de 285.6 m (IC = ± 98.42) y el individuo con la mayor distancia recorrida dentro de todo el paisaje fue un macho (marca 18A), con aproximadamente 564 m desde la E4 (bosque) hasta la E2 (borde). La tasa media de desplazamiento estimada para la especie fue de 0.070 m/días, mientras que para H. amphinome, el promedio de desplazamiento presentado fue de 278.5 m (IC = ± 66.56), y el individuo con mayor registro de vuelo fue el 230 B con 607 m recorridos desde E5 (bosque) hasta E6 (vegetación secundaria). Esta especie presentó la tasa media de desplazamiento más baja con 0.067 m/días (Tabla 3). Se destaca que no se registraron desplazamientos entre las E1 (cultivos mixtos con remanentes de vegetación) con el resto de estaciones para ninguna de las especies.

Relación entre la abundancia de Hamadryas con las variables ambientales y la cobertura dosel. Con los MLG se determinó que el uso del suelo y la época, así como su interacción con la CD, presentaron un efecto estadísticamente significativo (p < 0.05) sobre la abundancia de Hamadryas. Así mismo, se determinó que las variaciones de la CD, teniendo en cuenta los usos de suelo, influyen en la variación de la abundancia de la especie H. feronia. Con respecto a la especie H. amphichloe, tanto la época como los diferentes usos de suelo fueron significativos en la explicación de la variación de su abundancia, así como los cambios en la cobertura del dosel por la interacción de los muestreos (épocas) y los usos del suelo en la zona. Por otro lado, para H. februa y H. amphinome los usos del suelo influyen en la variación de la abundancia de estas especies (Tabla 4).

Discusión

El presente trabajo representa el primer estudio ecológico sobre la dinámica poblacional y el desplazamiento de las mariposas del género Hamadryas en diferentes usos del suelo en los Montes de María (San Jacinto, Bolívar). Además, se demuestra el papel de los fragmentos de BST y las matrices adyacentes de cultivos y vegetación secundaria en la dinámica poblacional de las poblaciones de mariposas en el Caribe colombiano (^{Pizano y García, 2014}).

La mayor abundancia de H. februa y H. feronia en usos de suelo como borde de bosque y cultivos mixtos permanentes puede estar relacionada con la disponibilidad de refugio o recursos como frutos en descomposición (fermentación). Por ejemplo, en la zona de cultivos mixtos permanentes (E7), hay presencia de árboles como M. indica, A. muricata, P. americana y palmas como A. butyracea que tienen cortezas grisáceas, los cuales les permiten hacer cripsis con este tipo de sustrato; así como líquenes adheridos que son utilizadas por este grupo de mariposas para camuflarse e interactuar con individuos de la misma especie (^{Jenkis, 1983}; ^{Mercado et al., 2018}; ^{Monge-Nájera et al., 1998}; ^{Prince et al., 2011}; ^{Sargent, 1985}; ^{Vargas-Zapata et al., 2015}). Adicionalmente, estas 2 especies prefieren zonas de bordes y cultivos con claros expuestos a la luz (observaciones en campo), donde aplican estrategias de comportamiento asociadas con la territorialidad (posar con las alas abiertas contra la superficie y con la cabeza apuntando hacia abajo), y la presencia de suelos desnudos con bastantes arcillas para obtener sales y nutrientes (^{Devries, 1988}; ^{Monge-Nájera et al., 1998}). La alta abundancia de estas 2 especies coincide con lo descrito por ^{Montero et al. (2009)}, ^{Vargas-Zapata et al. (2011}, ²⁰¹⁵⁾, ^{Cerpa y Flórez (2016)}, ^{Mercado et al., (2018)} y ^{Ahumada et al., (2019)}, quienes coinciden en que fueron las especies más comunes dentro de este género para zonas de bosque seco en el Caribe colombiano.

La poca abundancia de H. amphichloe y H. amphinome registrada dentro del paisaje puede estar relacionada con la estructura y composición vegetal de las estaciones de muestreos dentro del paisaje. Es posible que esta baja abundancia indique una rareza de hábitat para estas especies de mariposas, ya que pueden presentar interés por perchar y aprovechar otros sitios con algunas condiciones particulares, como troncos de colores oscuros y la presencia de barro con contenido de sales (^{González et al., 2014}; ^{Jenkis, 1983}; ^{Monge-Nájera et al., 1998}). De este modo, H. amphinome suele posar más frecuentemente en troncos con cortezas oscuras con DAP grande, localizados en el interior de zonas boscosas o adyacentes a márgenes de quebradas. Adicionalmente, esta especie presenta un patrón de coloración distintivo con franjas de color rojizo, las cuales las hacen más vistosas al depredador; razón por la que requiere de hábitats con vegetación más conservada y con áreas boscosas con varios estratos y bastante dosel (^{Moreno y Acuña, 2015}). Lo anterior explica la abundancia alta de esta especie en las estaciones ubicadas en zonas boscosas (E4 y E5) y en E6, esta última localizada en una zona de relictos secundarios de bosque seco con cierto tipo de cobertura vegetal.

La mayor abundancia de H. februa a mediados de marzo, puede atribuirse a que esta especie tenga sincronizado su ciclo reproductivo con las primeras precipitaciones, las cuales favorecen la productividad primaria del bosque, brindando así beneficios a este grupo de mariposas debido al aumento en la disponibilidad de recursos. La baja abundancia de Hamadryas durante las lluvias fuertes (muestreos 4 y 5, valores de precipitación entre 79.5 mm y 113 mm) de finales de abril y mediados de mayo se debe a que en la zona se registraron altos niveles de nubosidad (poca luminosidad), los cuales pueden producir disminución de la capacidad de vuelo y sus periodos de actividad (^{Kingsolver, 1983}). Adicionalmente, durante este periodo, debido a la influencia de las lluvias se deteriora el estado alar de las mariposas con envergadura grande como las Hamadryas (^{Jenkis, 1983}; ^{Monge-Nájera et al., 1998}), factor que incide negativamente en las rutas de vuelo, desplazamiento y en la disminución de sus poblaciones durante este periodo. Conjuntamente, con la influencia de las lluvias fuertes se disminuye el efecto del atrayente y de los frutos fermentados naturales del bosque, lo cual disminuye la probabilidad de captura de estas mariposas durante este periodo.

Los resultados observados indican que las poblaciones se comportan como poblaciones abiertas en el área y que éstas pueden estar reguladas por factores intrínsecos como son la tasa de natalidad, fecundidad, supervivencia y mortalidad (Graziano et al., 2002). De este modo, algunas zonas son utilizadas como fuentes de recursos (interior del bosque, zona de policultivos), mientras que los nuevos individuos dentro de las estaciones pueden deberse a un pico de emergencia de los imagos o a la entrada de nuevos individuos provenientes de áreas adyacentes. Así mismo, los diferentes usos de suelo representan atributos significativos para estas poblaciones incidiendo así en sus rutas de vuelo y en el flujo de entrada y salida de individuos, por lo que los adultos suelen asociarse con lugares que brinden un refugio con recursos óptimos para subsistir y completar su ciclo.

La variación en el porcentaje de recapturas entre cada una de las especies de Hamadryas puede estar relacionado con las tasas de dispersión de estas mariposas. Estas especies pueden experimentar distintas rutas de vuelo dentro de la matriz debido a la distribución aleatoria o al patrón de disposición de los recursos (^{Rueda y Andrade, 2016}; ^{Silva y Hernández, 2015}). Además, según ^{Jenkis (1983)}, ^{Monge-Nájera et al. (1998)} y ^{Murillo y Córdoba (2013)}, estas mariposas suelen desplazarse a distintos sitios de acuerdo con la disponibilidad de recursos alimenticios, zonas para la ovoposición y asociadas al camuflaje. Por ejemplo, la mayor frecuencia en las zonas con cultivos permanentes se puede atribuir a que estos lugares le brindan una alta cantidad de recursos que podrían ser explotados por especies como H. feronia y H. februa, con árboles aprovechables para su desarrollo y distintas actividades como el apareamiento. Así mismo, la zona de borde presentó una alta representatividad de individuos recapturados para cada una de las especies monitoreadas y esto estaría asociado con la distribución de especies vegetales típicas de las zonas de bosque y de áreas abiertas, factor que permite a los adultos de las Hamadryas desplazarse hacia el interior del bosque y utilizar esta zona como área de corredor biológico.

La residencia estimada en cada una de las especies presentó variaciones, lo cual puede atribuirse a que este grupo de especies puede tener varias generaciones durante el año (^{Cerpa y Flórez, 2016}; ^{Moreno y Acuña, 2015}). Estas especies se caracterizan por presentar ciclos de vida relativamente cortos, donde la mayor duración (alrededor de 15 días) la experimentan durante su fase larval (Brown, 1991; ^{Jenkis, 1983}; ^{Kremen, 1992}; ^{Pearson, 1994}); de esta manera, los tiempos de residencia muy cortos de los adultos en esta zona puede ser un indicativo del carácter multivoltino de estas especies. El menor tiempo de residencia de H. amphinome se atribuye a la fidelidad que tiene esta especie por permanecer en puntos fijos en el interior del bosque, donde se encuentren troncos que presentan tonalidades más oscuras y por estratos altos con bastante dosel (^{Monge-Nájera, 1998}).

Con respecto a las proporciones sexuales para cada una de las especies, la mayor proporción de machos coincide con lo descrito en lepidópteros, donde es más común encontrar una mayor proporción de machos que de hembras en los bosques tropicales (^{Ramos y Freitas, 1999}). El mayor número de machos estaría asociado a patrones en el comportamiento entre los sexos, lo que influye en su probabilidad de recaptura (^{Ehrlich y Gilbert, 1973}; ^{Marini-Filho y Benson, 2010}; ^{Young, 1982}). Los machos de estas especies se observaron emitiendo los sonidos típicos “clicks” en las diferentes estaciones de muestreo, lo que se asocia a funciones de defensa, territorialidad y cortejo (^{González et al., 2014}; ^{Yack et al., 2000}). Este comportamiento hace que muchos machos permanezcan cerca de sus territorios, lo que aumenta su probabilidad de ser recapturados en comparación a las hembras.

Por otro lado, el hecho de que los machos de cada una de las especies se hayan desplazados menos entre estaciones que las hembras, se asocia a que los primeros suelen patrullar y defender las zonas (territorialidad) que les brinden mayores recursos y frecuentar los sitios que les permitan ser más visibles para las hembras (^{García et al., 2002}; ^{Marini-Filho y Benson, 2010}; ^{Yack et al., 2000}). Además, algunas hembras del grupo se observaron en estratos medios y altos, asociadas a la búsqueda de sitios óptimos con plantas hospederas para la ovoposición, lo que puede aumentar las distancias de desplazamiento. Por otro lado, las hembras permanecen un menor tiempo en las zonas con respecto a los machos, debido a que muchas estarían en la búsqueda de recursos fundamentales para el proceso de maduración de huevos como sales, nutrientes, ácido úrico y azucares (^{Baker, 1972}; ^{Gilbert, 1975}; ^{Gilbert y Singer, 1975}; ^{Jones, 1884}; ^{Nieves et al., 2015}; ^{Sterry, 1997}).

La variación temporal de la abundancia de Hamadryas en la VLF (Apéndice), indica que hay una tendencia a aumentar con las primeras lluvias. Conforme llegan las lluvias, la fenología del paisaje cambia (mayor cobertura vegetal), incidiendo de manera indirecta en la dinámica poblacional de estas mariposas, especialmente en los muestreos posteriores (^{Vargas-Zapata et al., 2011}). A su vez, en la transición de época seca a lluvias se observaron cambios en la fenología y cobertura de la vegetación como arbustos, plantaciones y enredaderas (observaciones en campo), lo que permite que el paisaje sea más homogéneo y el desplazamiento de estas mariposas se incremente. Esto concuerda con ^{Ramírez et al. (2014)} y ^{Vargas-Zapata et al. (2015)}, quienes reportan mayor abundancia de mariposas Hamadryas para BST en la época de las primeras lluvias. La mayor frecuencia de H. feronia y H. februa durante la época seca se debe a que estas especies pueden aprovechar las zonas abiertas con altos valores de intensidad lumínica, debido a sus condiciones heliotérmicas. A su vez, ésto concuerda con ^{Ribeiro et al. (2012)}, quienes indican que la distribución y la abundancia de mariposas no son aleatorias, sino que siguen un patrón de distribución de recursos que está relacionado con aquellos hábitats con condiciones óptimas para llevar a cabo sus funciones biológicas.

En conclusión, las poblaciones de Hamadryas presentaron variación espacial (usos) y temporal (muestreos) en la vereda La Flecha; lo que demuestra que la heterogeneidad del paisaje y la estacionalidad juegan un papel importante en la dinámica poblacional de estas especies. Las poblaciones de Hamadryas se desplazan entre áreas de bordes, plantaciones y vegetación secundaria dentro del paisaje, sugiriendo que el movimiento de los individuos está determinado por los cambios de cobertura entre las distintas áreas. Sumado a lo anterior, el desplazamiento de las hembras está asociado a los patrones de disposición de las plantas hospederas y distribución de los recursos óptimos para la ovoposición.