December
02, 2024
May
, 2023
We studied bats in a natural protected area and the influence zone in Chiapas, México. Focusing on small vesper bats (Vespertilionidae: Myotinae) in water sinkholes, flooded surfaces, and forested locations (sampling sites), our objectives were to: 1) differentiate the sampling sites based on abiotic and biotic variables, 2) compare the relative activity between the contrasting sampling sites, and 3) determine the relative importance of surrounding land covers amounts on the relative activity. We expected the following: 1) that sampling sites would be differentiable based on an interaction of environmental conditions and nocturnal flying insects’ biomass, 2) to find a significantly higher relative activity associated with water bodies as a reflection of potential drinking and feeding supplies, and 3) that forest cover would have a significant association in a positive direction with the relative activity as a reflection of roosting and commuting opportunities. We obtained weather data with a portable station, collected flying nocturnal insects with a passive trap, and recorded bats with ultrasonic detectors. We also used satellite imagery to calculate land covers amounts around the sampling sites. We performed a multivariate analysis to compare sampling sites, and elaborated correlation models of relative activity against surrounding land covers amounts. The only significant variable for differentiating sampling sites was wind speed; in water sinkholes, with lower mean value compared to flooded surfaces and forested locations, along with a comparatively higher percentage of dipterans in the samples. The mean relative activity of bats was significantly higher in water sinkholes, with maximum values reaching 95 %, which we can relate to environmental conditions and resources available. We found associations in a positive direction with forest cover, secondary forest and water surface, and negative with agricultural land and human development. Finally, we argue that some of the water sinkholes meet the criteria of small natural features supported by the unusual activity of the studied bats, and that targeted conservation actions will complement other strategies implemented in the area.
Keywords::
Forest loss, hierarchical partitioning, landscape level, local level, Myotis spp., small natural features, targeted conservation
The accelerated loss of biodiversity because of human activities is a concern worldwide. Therefore, identifying which attributes in agricultural zones with a high rate of forest loss can back away from the decline of animal populations is crucial (Heim et al.2015). This way, multi-level management may contribute to more suitable habitats for wild species.
For aerial insectivorous bats, localized resources (e. g., water, prey) can represent limiting factors for survival and reproductive success (Findley 1993). At the same time, landscape structure (e. g., forest amount and its spatial arrangement, linear elements, among other attributes) may facilitate commuting toward drinking and feeding sites (Fuentes-Montemayoret al.2017; Mendeset al.2017; Martinoet al.2019).
Fahrig (2003) argues that the number of individuals of any species should be, to some point, a positive function of the quantity of habitat available, such as native forests. The strong forest dependency some species of bats exhibit depends on roosting preferences and flight restrictions in open areas (Lackiet al.2007; Fuentes-Montemayoret al.2013; Parreira Peixotoet al.2018; Novella-Fernandezet al.2022). Forest loss can decrease structural connectivity for many species (Fischer and Lindenmayer 2007). However, this topic needs to be better studied in the Neotropics, specifically considering slow-flying bats in upland ecosystems.
Flying away from vegetation cover may impose higher energetic costs for slower species (due to the stronger wind), such as vesper bats (Vespertilionidae), characterized by wings with low aspect ratio and loading (Norberg and Rayner 1987; Heim et al.2015). Moreover, their short and high-frequency modulated echolocation pulses are not well suited to open areas because they are more subject to atmospheric attenuation (Pettersson 2002; Jones and Rydell 2003; Frey-Ehrenboldet al.2013; Heim et al.2015).
We studied small vesper bats (Vespertilionidae: Myotinae) in water sinkholes, flooded surfaces, and forested locations (sampling sites), inside a natural protected area and the influence zone in Chiapas, southeast México. We based the selection of this group of bats on traits such as wing morphology and echolocation system related to the sensitivity to forest loss (Heim et al.2015). They also are essential predators of nocturnal flying insects (e. g., mosquitoes), which are potentially deleterious to humans and become pests, which should be a key consideration in sustainable land management (Gonsalveset al.2013; Heim et al.2015; Puig-Montserratet al.2020).
Our objectives were to 1) differentiate the sampling sites based on abiotic and biotic variables, 2) compare the relative activity between the contrasting sampling sites, and 3) determine the relative importance of surrounding land covers amounts on the relative activity. We expected 1) that sampling sites would be differentiable based on an interaction of environmental conditions and nocturnal flying insects’, 2) to find a significantly higher relative activity associated with water bodies as a reflection of potential drinking and feeding supplies, and 3) that forest cover would have a significant association in a positive direction with the relative activity as a reflection of roosting and commuting opportunities.
Study area and sampling points. The study area is intermediate between lowlands and highlands, in the comiteca plateau (Meseta Comiteca), Chiapas, México (Figure 1). It includes the Lagunas de Montebello National Park (PNLM, by its initials in Spanish), listed in the RAMSAR convention (no. 1325). The main vegetation association inside the PNLM is a mixed forest dominated by pines (Pinusspp.) and oaks (Quercusspp.), with a canopy stratum between 20 to 35 m and an interior species-rich tree stratum below 20 m height (González and Ramírez 2013). The predominant land cover outside the PNLM is agricultural land. There is a dry season (approximately March to June), a rainy season (approximately July to October), and a windy season (approximately November to February).
The sampling points were selected based on physiognomy and pre-sampling, with three repetitions each: water sinkholes, which are open areas with surrounding vegetation and have a permanent water surface several meters below the ground (Figure 2a), flooded surfaces, which are also open areas with surrounding vegetation, and have a fluctuating water surface that is superficial, and represent natural and human-induced elements (Figure 2b), and forested locations, such as clearings and trails inside the forest (Figure 2c). There were no considerable elevation differences in the sampling points, which started from 1,460 to 1,540 m, and accounting for extent restrictions for the spatial analysis, the minimum distance between them started from ~ 3.5 km. The samplings took place from July 2014 to July 2015 to obtain data on environmental conditions, insect collections, and bat recordings simultaneously at each event.
Satellite imagery processing and land covers amounts. We used multispectral SPOT 5 orthorectified satellite images (spatial resolution of 10-m pixels) of 2015 (dry season). The latter images were subject to radiometric calibration and were provided by the Laboratorio de Información Geográfica y Estadística (LAIGE, by its initials in Spanish) of El Colegio de la Frontera Sur. We performed composition in false color using bands 3 (near infrared), 2 (red), and 1 (green), with a simple linear contrast. The land covers (classes) defined were (Figure 1): forest cover (mature forest), secondary forest (such as coffee crops and orchards), agricultural land (extensive crops such as maize and grasslands), human development (including urbanization and bare soil), and water surface (lakes, water sinkholes, and flooded surfaces). Following Fuentes-Montemayor et al. (2013), we grouped distinct forest types to avoid confounding information (e. g., adjacent types with undistinguishable limits).
We performed a supervised classification in Idrisi version 17.0 (Clark Labs, Clark University). Following Lu and Weng (2007) and Eastman 2012, we implemented a segmentation routine that identifies adjacent pixels grouped by spectral similarity, so-called objects. We used a similarity threshold of 30 to obtain homogeneous segments of a size facilitating the selection of training samples. Then, we used the maximum likelihood classifier. Finally, we re-classified the image to improve the precision of the classification and to produce smooth edges between classes by using a distinctive classifier of the segmentation routine. An error matrix (see Verbyla 1995) assessed the overall classification accuracy from 100 field control points spread throughout the study area, obtained during 2014 to 2015, yielding an accuracy of 85 %.
We calculated the area (ha) and percentage of each land cover in concentric buffers of 500-m (78.6 ha), 1,000-m (314.2 ha), and 1,500-m (706.8 ha) radii around the sampling points (Supplementary material 1), delimited in ArcGIS version 10.2.1 (ESRI, Inc.). We were looking to encompass the home range of small vesper bats (e. g., Myotis spp.; Owen et al. 2003; Coleman et al. 2014). Also, the nested design allowed us to explore responses by changing spatial scales and the non-overlapping buffers to avoid re-measuring land covers or pseudoreplication (Popescu and Gibbs 2010; Arroyo-Rodríguez et al. 2016).
Environmental conditions. We used a multi-function weather meter WM-350 WindMate® (WeatherHawk) to obtain monthly data on temperature (°C), relative humidity (%), and wind speed (KMPH). We took measurements three times per night at the beginning, half, and end of the sampling events; this way, we calculated means for comparison between sampling points. We used the mean calculation per night as input for multivariate analyses.
Insects’ biomass. For the collection of nocturnal flying insects, we used a Malaise trap. The latter is a passive trap, which we use to not interfere with recordings by attracting prey. We suspended the trap on trees near the sampling points and sampled around water surfaces near the shore. The trap was installed 1.5 m above ground. All specimens were preserved in 70 % alcohol inside plastic containers (Wickramasingheet al. 2004) and labeled for posterior processing and identification.
In the laboratory, we identified the specimens into the orders Diptera (dipterans) and Coleoptera (coleopterans). We pooled Trichoptera and Lepidoptera orders into the superorder Amphiesmenoptera (amphiesmenopterans) because we faced difficulties identifying alcohol-preserved samples. These orders are among the most common food elements for aerial insectivores, particularly this group of bats (see Whitaker 2004; Segura-Trujillo et al. 2018), and may reflect potential prey availability.
We counted the number of individuals of each taxon in each sample (see Queiroz de Oliveira et al. 2015), which were dried later in a stove at 70 °C for 48 h (Bradley et al. 1993), and obtained biomass (g) using an Explorer™ Pro Analytical Balance (EP214C), with a readability of 0.1 mg (Ohaus Corporation). Following Queiroz de Oliveira et al. (2015), we divided it by the number of insects for standardized measurement. Finally, we calculated the mean per night at each sampling point.
Relative activity of bats. We used an Echo Meter EM3+ Ultrasonic Detector (Wildlife Acoustics, Inc.) to obtain full-spectrum bats recordings. We performed short-term passive recording sessions lasting 4 hrs or else standardized. We configured the EM3+ with 1) sample rate of 256 kHz, 2) WAV audio file format, 3) maximum duration of 15 sec, 4) frequency trigger of 15 kHz, 5) amplitude trigger of 18 dB, 6) trigger window of 1 sec, and 7) gain of 30 dB. We positioned the detector on a 1-m pole with a 45° upward angle, directed toward acoustic space in water sinkholes, flooded surfaces, and forested locations (gaps and trails). Each sampling point was visited twice during three consecutive nights. We stored the recordings in a 32 GB SD card (Kingston© Technology Corporation), and the total recording effort was 432 hrs.
We were interested in small vesper bats flying through the acoustic space of the sampling points. These bats exhibit slow and high-maneuverability flight determined by wing morphology, and high-frequency modulated echolocation pulses of short-range easily attenuated in open areas (Frey-Ehrenbold et al. 2013; Bader et al. 2015), though facing similar restrictions. Some of the species contained in this group of bats can be acoustically cryptic. Identification of those with similar body sizes can be challenging due to the high similarity of the echolocation pulse’s structure, and acoustic parameters overlap (see Jung and Kalko 2011; Williams-Guillén and Perfecto 2011; Estrada-Villegas et al. 2012). Therefore, we pooled recordings to avoid misclassifications.
The small vesper bats potentially occurring in sympatry in our study area are Rhogeessa tumida, M. nigricans, and M. keaysi pilosatibialis (Barquez and Diaz 2016; Miller et al. 2016; Solari 2019). They belong to the same ecomorphotype and foraging type (Fenton and Bogdanowicz 2002; Segura-Trujillo et al. 2018), classified as aerial/trawling insectivores that hunt in background clutter space such as forest edges and gaps (Schnitzler and Kalko 2001; Frey-Ehrenbold et al. 2013). Accordingly, we selected echolocation pulses starting from a final frequency of 45 kHz. Jung and Kalko (2011) and Estrada-Villegas et al. (2012) indicate the possibility of erroneously classifying Rhogeessa tumida as M. albescens or M. nigricans. We can also find the distribution limits of M. albescens and M. fortidens (see www.iucnredlist.org), but they are more associated with lowlands, so we discard them. Other Myotis species, such as M. velifer and M. californicus, produce lower final frequencies (see Orozco-Lugo et al. 2013; Zamora-Gutiérrez et al. 2016). Finally, Myotis elegans emits frequencies above 60 kHz (O’Farrell and Miller 1999), but it was not detected.
For recording processing, in the Bat Analysis Mode, we configured automatic classification using Auto ID for Bats-Bats of the Neotropics (México) with a neutral level of sensitivity, as we were interested in more identifications and not highly accurate ones (see User Guide). We automatically filtered noise files. We processed the recordings in Kaleidoscope Pro v. 5 (Wildlife Acoustics, Inc.). With Kaleidoscope Viewer, we displayed and verified the resulting recordings because automatic identification may produce false positives (see Auto ID for Bats; Rydellet al.2017).
For the calculation of relative activity (a surrogate of bat abundance; Froidevaux et al. 2021), we considered the activity index proposed by Miller (2001). The latter index is based on the total 1-min blocks with evidence of echolocation pulses during constant periods; in our case, for a total of 240 min (4 hrs), otherwise standardized to the total minutes sampled. We considered evidence of at least one recording with a minimum of two consecutive echolocation pulses (MacSwiney et al. 2009; Heim et al. 2015). We expressed the calculations as percentages. For comparisons, we calculated the mean percentage of relative activity at each sampling point for six events.
Environmental conditions and insects’ biomass analysis. We were interested in differentiating sampling sites based on interacting abiotic and biotic variables. Therefore, we performed a stepwise discriminant analysis. We included the following variables: dipterans biomass, coleopterans biomass, amphiesmenopterans biomass, temperature, relative humidity, and wind speed. We based this analysis on Wilk’s Lambda (λ), the F statistic’s significance, and the independent contributions of variables in the model through Partialλand theF-removestatistic’s significance. We performed these calculations in STATISTICA® version 8.0 (StatSoft, Inc.).
Surrounding land covers amounts analysis. Hierarchical partitioning (hp) is a multivariate exploratory analysis that explains variance in the response variable attributable to univariate correlations with each independent variable (Radford and Bennett 2007). We employed hp to measure the relative importance of surrounding land covers amounts (Chevan and Sutherland 1991; Mac Nally 2000). This analysis has the advantage of addressing potential multicollinearity (Olea et al. 2010).
We performed hp with routine hier.part included in package hier.part version 1.0-6 (Walsh and Mac Nally 2022) ran in R version 3.6.2 (R Core Team 2019). The response variable was the mean percentage of active 1-min blocks, and the independent variables were the percentage of each land cover. The variables were arcsine-transformed [ASIN(SQRT(x/100)], and we ran the models specifying a gaussian distribution, with the goodness of fit based on R 2 . We tested the statistical significance by comparing randomization (1000 permutations) with routine rand.hp included in package hier.part version 1.0-6 (Mac Nally 2002; Walsh and Mac Nally 2022) ran in R version 3.6.2 (R Core Team 2019). We used the generated Z-scores to establish statistical significance based on the upper 95 % (Z ≥ 1.65) confidence limit (Walsh and Mac Nally 2022). The hp does not indicate the direction of paired associations, which we determined by non-parametric correlations.
Sampling sites. We provide descriptive statistics of environmental conditions, insect collections, and single-variable statistical comparisons between sampling sites (Supplementary materials 2, 3, and 4). The stepwise procedure resulted in a statistically significant model, including only the wind speed, for the differentiation between sampling sites (λ = 0.5259, F (2.24) = 10.8143, P = 0.0004), which exhibited the lowest mean values in water sinkholes in all measurements. We include the calculations for the variables not in the model in Supplementary material 5. The stepwise procedure excluded insect collections. However, we observed a comparatively higher percentage of dipterans in water sinkholes in the samples (Supplementary material 3).
Relative activity of bats. We include the calculations of relative activity in the sampling sites in Supplementary material 6. The difference was statistically significant between sampling sites (K-W = 28.78, P = 5.638E-7; Figure 3), specifically between water sinkholes and flooded surfaces (Dunn post-test, mean rank difference 14.750, P = 0.002), with higher mean value in water sinkholes; between water sinkholes and forested locations (Dunn post-test, mean rank difference 27.250, P = 5.228E-6), with higher mean value in water sinkholes; between flooded surfaces and forested locations (Dunn post-test, mean rank difference 12.500, P = 0.004), with higher mean value in flooded surfaces.
Surrounding land covers amounts. Invariably, we observed associations in a positive direction with forest cover, secondary forest, and water surface. We observed associations in a negative direction with agricultural land and human development, except for the latter class in the 500-m buffer (9 % of the variance), where a small surface characterized it. In the 500-m buffer, we established that the relative activity is mainly determined by forest cover and secondary forest, explaining 25 % and 26 % of the variance (respectively), and agricultural land explaining 34.7 % of the variance (Figure 4). In the 1,000-m buffer, we established that the relative activity is mainly determined by forest cover and secondary forest, explaining 36.4 % and 20.6 % of the variance (respectively), and agricultural land explaining 23.8 % of the variance (Figure 4). In the 1,500-m buffer, we established that the relative activity is mainly determined by forest cover and water surface, explaining 20.3 % and 42.9 % of the variance (respectively), and agricultural land explaining 28.7 % of the variance (Figure 4). We did not find statistically significant associations based on the upper 95 % confidence limit (Supplementary material 7).
In the study area, we could differentiate water sinkholes, flooded surfaces, and forested locations based on the wind speed, with the lowest mean values in water sinkholes in all measurements. Although the insect collections did not stand out in the model, we observed a comparatively higher percentage of dipterans in water sinkholes in the samples. Some dipterans are strongly associated with calming water because their pupae and larvae are susceptible to increasing wind speed (Gillies and Wilkes 1981; Rutledge 2008). The characteristics of water sinkholes, specifically the surface several meters below the ground, can be associated with environmental conditions’ buffering beneficial for this kind of insect.
Many Myotis species consume small and soft prey, specifically the aerial ecomorphotypes such asM. nigricansandM. pilosatibialis, which exhibit short wavelength echolocation suitable for these targets (Gonsalveset al.2013), and biomechanical limitations (bite strength) correlated to body size (Segura-Trujillo et al. 2018). Their diet can be composed of several orders, including dipterans (Aguiar and Antonini 2008; Gamboa Alurralde and Díaz 2019; Ingalaet al.2021). However, bats eat a wide range of insects, most representatives seem to have flexible diets (Jones and Rydell 2003), and there is evidence that suggests that the diet responds to local fluctuations, as well as the abundance and type of prey (Salinas-Ramoset al.2015).
In water sinkholes, the small vesper bats exhibited a comparatively high relative activity in two of the three sampling sites and in most of the recording events, reaching up to 95 % of the time, specifically in AZAR (Supplementary material 6), and on average, it was significantly higher compared to flooded surfaces and forested locations. We also obtained a high rate of buzzes, which we determined by the output of the detector and subsequently by visual inspection of the recordings; these might represent drinking and/or feeding buzzes, but a detailed analysis distinguishing them based on their structure is necessary to make adequate interpretations (see Russoet al.2015). The physiological characteristics of bats demand high amounts of water and prey for reproductive success (MacSwiney et al. 2009; Seiboldet al.2013; López-Gonzálezet al.2016), and there is presence of at least one maternity colony inside the natural protected area, relate to our observations.
Invariably, we observed associations in a positive direction with forest cover and secondary forest. The latter associations could reflect feeding and commuting opportunities. Other studies report similar results, specifically positive correlations between the activity and abundance of this kind of bats and forest cover (native and planted forests) and higher activity in locations at a smaller distance to forested areas (Heim et al.2015; Rodríguez-San Pedro and Simonetti 2015; Fuentes-Montemayoret al.2017; Putet al.2019; Laurindoet al.2020; Falcãoet al.2021).
Fuentes-Montemayoret al.(2013) make note that the use of forests by aerial insectivorous bats can also be influenced by attributes such as tree density, among other forest structure variables, conditions that Rauchenstein et al. (2022) define as the “suitable foraging habitat”. Vegetation clutter (e. g., forest canopy openness or increased shrub cover) can influence prey abundance (Froidevauxet al.2021; Rauchenstein et al. 2022), restrict mobility and flight maneuverability (Estrada-Villegas et al. 2012; Fuentes-Montemayoret al.2013), and forest maturity (e. g., tree sizes) can be a limiting factor for tree roosting bats (Novella-Fernandez et al. 2022).
We also observed an association in a positive direction with water surfaces, represented in the study area by lakes, water sinkholes and flooded surfaces. The latter elements could represent stepping-stones for many taxa in agriculture-dominated landscapes (Hunter Jr.et al.2017). Even for highly mobile organisms such as bats, a dense network of connecting elements might be beneficial and promote activity in open areas (Heim et al.2015). In the UK, Fuentes-Montemayoret al.(2013; 2017) found positive correlations between the abundance ofMyotisspecies to a larger proportion of surrounding water and decreasing distance between water bodies. In the Atlantic Forest of Brazil, Laurindoet al.(2020) found that the number of captures of insectivorous bats, including Myotisspecies, was significantly associated with the area covered by water bodies within highly fragmented agricultural landscapes.
We observed associations in a negative direction with agricultural land and human development, except for the latter class in the 500-m buffer, where a small surface characterized it. The latter associations could reflect the effect of an increasingly open area and unfavorable conditions for bats. For instance, Fuentes-Montemayoret al.(2017) found higher activity ofMyotisspecies in forests immersed in areas with a smaller proportion of urban areas. On the other hand, Laurindoet al.(2020) found a negative correlation between the number of captures of insectivorous bats with increased agricultural area.
We explain the associations of relative activity to landscape covers amounts by natural history traits, particularly wing morphology, which is a surrogate for mobility, and echolocation system, which is a surrogate for perceptual range (Frey-Ehrenboldet al.2013; Baderet al.2015; Heim et al.2015; Fuentes-Montemayoret al.2017). The morphological and echolocation system traits of sympatric small vesper bats in the study area point to slow flight with more energetic cost in open areas and short-range high-frequency pulses that are more attenuated in open areas (Frey-Ehrenboldet al.2013; Baderet al.2015; Heim et al.2015). However, Fuentes-Montemayoret al.(2017) suggest that low-mobility species are more influenced by local conditions and the landscape becomes more important for high-mobility species, which perceive the environment at a coarser scale. The latter argument could support our observations, specifically the strong local-level response of small vesper bats to the presence of water sinkholes and the lack of significance of the associations at the landscape level.
Final remarks. The presence of water sinkholes and environmental conditions such as low wind speed possibly regulating prey availability are important drivers of small vesper bats’ activity in the study area. We found no significant associations for flooded surfaces and forested locations. However, flooded surfaces are particularly relevant to other groups of bats (e. g., Mormoopidae) and vertebrate fauna in the influence zone, such as anurans and migratory birds, representing a seasonally limiting resource.
We observed trends of positive correlation between relative activity to forest cover, secondary forest, and water surface and negative correlation to agricultural land and human development. Although we found no significant associations, we provide explanations based on the argument that small vesper bats react negatively to forest loss and increasingly open area, which we explain by wing morphology and echolocation system traits. Other authors suggest that, for low-mobility species, the landscape may become less important than local-level attributes.
Some of the water sinkholes analyzed in the study area meet the criteria of small natural features having ecological importance that is disproportionate to their size (Hunter Jr. 2017a; Hunter Jr.et al.2017b), as we determined for small vesper bats. The latter consideration is valuable for complementing large-scale conservation through targeted actions and should be further evaluated for its implementation in the area. The water sinkholes are located inside and right on the limits of the natural protected area, in the jurisdiction of federal authorities. The flooded surfaces and forested locations studied are communal goods in the influence zone exposed to human activities and their effects, such as pollution derived of the use of pesticides in the agricultural zones nearby, water exploitation for crops irrigation and cattle drinking supply, forest exploitation with consequent degradation and deforestation.
Finally, we point out that the adequate management of these types of landscape elements should be prioritized in the conservation agenda of the area in an effort of a consensus with ejidatarios (local authorities and ejido members). There is a need for a dialogue that should follow the socio-ecology premises, specifically a transdisciplinary approach to solving environmental problems, looking to benefit biodiversity, the continuity of ecosystem services, and human well-being.
We would like to thank the authorities of the Comisión Nacional de Áreas Naturales Protegidas (CONANP) for the permission to work in Lagunas de Montebello National Park. We are also grateful to landowners and local authorities for work permissions. This project was funded partially by Bat Conservation International (grant number SS14012, to the first author), Idea Wild (to the first author), and Instituto Politécnico Nacional, México (grant number SIP 20151481, to the third author). The first author received a scholarship awarded by the Consejo Nacional de Ciencia y Tecnología, México (number 244510). Finally, we would like to thank the anonymous reviewers for their contribution to improving this manuscript.