- Open Access
Patterns of foraging activity and fidelity in a southeast Asian flying fox
Movement Ecology volume 8, Article number: 46 (2020)
Improved understanding of the foraging ecology of bats in the face of ongoing habitat loss and modification worldwide is essential to their conservation and maintaining the substantial ecosystem services they provide. It is also fundamental to assessing potential transmission risks of zoonotic pathogens in human-wildlife interfaces. We evaluated the influence of environmental and behavioral variables on the foraging patterns of Pteropus lylei (a reservoir of Nipah virus) in a heterogeneous landscape in Cambodia.
We employed an approach based on animal-movement modeling, which comprised a path-segmentation method (hidden Markov model) to identify individual foraging-behavior sequences in GPS data generated by eight P. lylei. We characterized foraging localities, foraging activity, and probability of returning to a given foraging locality over consecutive nights. Generalized linear mixed models were also applied to assess the influence of several variables including proxies for energetic costs and quality of foraging areas.
Bats performed few foraging bouts (area-restricted searches) during a given night, mainly in residential areas, and the duration of these decreased during the night. The probability of a bat revisiting a given foraging area within 48 h varied according to the duration previously spent there, its distance to the roost site, and the corresponding habitat type. We interpret these fine-scale patterns in relation to global habitat quality (including food-resource quality and predictability), habitat-familiarity and experience of each individual.
Our study provides evidence that heterogeneous human-made environments may promote complex patterns of foraging-behavior and short-term re-visitation in fruit bat species that occur in such landscapes. This highlights the need for similarly detailed studies to understand the processes that maintain biodiversity in these environments and assess the potential for pathogen transmission in human-wildlife interfaces.
Foraging can be viewed as one of the most fundamental activities of wild animals  because it contributes significantly to maintaining physiological functions and individual fitness [2,3,4]. Foraging patterns vary widely at inter- and intra-specific levels [5,6,7] and can be characterized by differences in foraging localities and fidelity [8, 9], as well as timing, duration and rate of foraging bouts [10, 11]. More generally, they are typically assumed to be a function of specific physiological needs, intrinsic individual features and particular environmental conditions.
Habitat quality and more specifically food-resource characteristics (i.e. diversity, abundance, nutritional value, and availability dynamics) are known to influence patterns of foraging behavior [12, 13]. For instance, animals may concentrate and increase their foraging effort in habitats where high-quality resources are available [14, 15]. Foraging animals may also revisit highly-profitable areas more often when food availability is predictable over time [16, 17]. As a consequence, habitat transformation and landscape heterogeneity can influence their foraging behavior substantially [18, 19], by forcing or allowing them to adjust their strategies which in some cases can result in local extinction or population growth [20, 21]. Because human activities are a major driver of rapid habitat modification , research on the foraging ecology of species at human-wildlife interfaces is crucial to understanding foraging patterns and predicting their evolution in a conservation framework [23,24,25].
Bats are highly represented among mammals  and provide ecologically important and economically significant services which include insect-pest control, plant pollination, and seed dispersal [27,28,29]. Almost half of all living bat species are threatened by substantial population declines or extinction . A major cause of this is loss and degradation of roost and foraging resources [30, 31] due to widespread transformation of natural habitats (e.g., agricultural expansion and urbanization). As a result, improved understanding of the foraging ecology of bats in the face of ongoing habitat transformation is important for effective bat conservation. Further, because some bat species are hosts for important zoonotic viruses  such as lyssaviruses, filoviruses, coronaviruses and henipaviruses (Hendra and Nipah viruses) which can be transferred through contact with bat bodily fluids (e.g., saliva, urine and feces) , knowledge of bat foraging ecology is also essential to assessing potential risks of zoonotic transmission at human/bat interfaces. This presents challenges because bats are primarily nocturnal and highly mobile animals (i.e. powered-flight ) which cannot be tracked by unaided human vision. This has historically rendered monitoring of individual behavior difficult, although recent advances in satellite telemetry have partly overcome these challenges, in enabling data collection on individual locations over time .
We investigated the foraging ecology of Lyle’s flying fox (Pteropus lylei), a bat species belonging to the Pteropus genus which is widely distributed in Southeast Asia . Pteropus bats face huge habitat transformation across their range due to large-scale deforestation, mainly for agriculture and, to a lesser extent, increased urbanization [30, 31]. They have also been identified as a reservoir for Nipah viruses that induce severe encephalitis in humans with high fatality rates [37, 38]. The ecology of most Pteropus species is relatively poorly known, although some taxa have evidently adapted to anthropogenic environments. For instance, some species populations mainly forage in human-shaped landscapes  and include cultivated or exotic plant species in their diet [40, 41].
We employed GPS (Gobal Positioning System) devices to collect spatio-temporal data on individual P. lylei and investigate their foraging ecology in a highly anthropogenic and heterogeneous environment in Cambodia. Based on GPS locations (presence data), random locations (pseudo-absence data) and environmental factors, a previous study investigated habitat selection in the same population to predict its overall distribution in the study region . This found that residential areas were the preferred foraging habitat for the species, followed by plantations and unmanaged tree vegetation. In addition, P. lylei bats have been shown to forage in several areas each night, undertaking small-scale movements in each foraging area and frequently returning to certain areas [42, 43]. However, these individual behaviors have yet to be studied in depth to characterize the foraging behavior of the species and understand the resulting distribution patterns. More specifically, analyzing animal movements at the path-level should enhance knowledge into proximal mechanisms of animals space use [44, 45]. To this end, we employed an approach based on animal-movement modeling, which comprised a path-segmentation method (hidden Markov model) to identify individual foraging-behavior sequences from GPS data. We characterized foraging localities, foraging activity within nights, and the probability of a bat revisiting a given foraging locality over consecutive nights. We also assessed the influence of several variables on these behaviors such as proxies for energetic costs and quality of foraging areas.
Study region and population
The study was conducted in 2016 within the Koh Thom District of Kandal Province in Cambodia. The landscape of the region is characterized by habitats which include agricultural lands (mainly wet rice fields), fruit plantations (mainly sapodilla Manilkara zapota and mango Mangifera indica), tree vegetation (e.g., tree groves, flooded vegetation and scrublands), and residential areas (e.g., houses and backyards with a high diversity of fruit trees, including mango, bananas (Musa spp.), sapodilla, longan (Dimocarpus longan), palms (Arecaceae), kapok (Ceiba pentandra), guava (Psidium guajava) and figs (Ficus spp.) ). The studied population of Lyle’s flying foxes roosted in a grove of 21 trees located in the ‘Wat Pi Chey Saa Kor’ Buddhist pagoda (11.200 N, 105.058 E) of the Kom Poung Kor village. Between 4000 to 7000 individuals typically roost at the site during the day (but some of which occasionally roost at nearby sites), and consume fruits obtained from trees within (or on the boundaries of) the various habitats comprising the landscape [42, 46, 47].
We caught P. lylei during their annual birthing period at the roost site, from 18th April to 6th May 2016 (which encompassed late pregnancy and early lactation). This period corresponds with high Nipah-virus circulation in P. lylei in Cambodia  and Thailand . Selected individuals (adults only, excluding reproductively active females) were temporarily equipped with GPS devices (model: FLR V, Telemetry Solution™; spatial accuracy: 5 m; weight: 20 g) which were attached to their necks with nylon collars. The GPS loggers represented on average 4.2% ± 0.6 SD (range: 3.4–4.8%) of individual body mass (486.1 g ± 66.5, range, 414–590 g). Details on the timing of captures, anesthesia, sex, body mass, collar attachment and number of bats caught are given by .
Spatial data from the GPS loggers were collected for several consecutive nights (18:00–05:00) for each bat and were retrieved at the roost site in daytime using a receiving station with a maximum transfer distance of 20 m (these data are available in movebank.org; study name: “Foraging movements of Lyle’s flying foxes in Cambodia”: ). Acquisition rates for spatial data were heterogeneous among individuals . Because  found that Pteropus spp. spent an average 12.0 min (range: 1–67 min) and 25.0 min (range: 10–40 min) searching for and consuming fruits in fig and mango trees respectively, only GPS devices recording one location every 5 min were considered. The present analyses involved seven males and one female for which data were obtained for 3 to 12 consecutive nights from 22nd April and 17th May (one additional bat was excluded due to insufficient data). The number of nights for which data were collected from a given GPS device depended on several factors including battery life and proximity of bats to the receiving station. The first night of data was excluded from analyses to avoid potential behavioral biases related to the capture event and because the GPS devices were programmed to record at a lower rate on the first night . Seven nights that had more than 10% of values missing and five nights where the intervals between two locations were greater than 5 min were also excluded. For the remainder of the dataset, occasional missing values (< 2% of all data) were interpolated using continuous-time correlated random walk modelling (; ‘crawl’ R package) to maintain temporal regularity between locations. The nightly movements of bats are presented in the supplementary material (Fig. S1).
The nightly foraging activity of Pteropus species includes several types of behavior [39, 43]: (i) long-distance movements (hereafter referred to as commuting flights, CFs) between roost site and foraging areas, and between distinct foraging areas, (ii) short-distance movements between patches of high resource density within a given foraging area to collect food (foraging flights, FFs), and (iii) stationary bouts in a given foraging patch for fruit consumption or resting (stationary bouts, SBs). All of our GPS locations were recorded during one of these three types and the corresponding behavioral state (i.e. commuting flight, foraging flight, or stationary bout) was attributed to each location using a hidden Markov model (HMM; using the ‘moveHMM’ R package). This modelling approach allowed us to determine the composition of behavioral states in temporal sequences of spatial data, based on the distribution of step lengths and turning angles between consecutive locations ([44, 51]; see also studies in echolocating bats: [52, 53]). For this purpose, parameters (means) characterizing the distribution of movement metrics were specified in the modelling procedure for each behavioral state as follows (for step lengths and turning angles, respectively): 2000 m and 0° for commuting flights, 400 m and 0° for foraging flights, and 30 m and 180° for stationary bouts. These values were estimated from the distribution histograms of total step lengths and turning angles (see supplementary material, Fig. S2), and were consistent with our knowledge of Pteropus movements (e.g., low-speed and more circular movements for foraging flights and stationary bouts in comparison to commuting flights: [42, 43]). Following , the model was fitted with a gamma step length distribution and a von Mises turning angle distribution. As the model estimates the probability of each behavioral-state (commuting flights, foraging flights or stationary bouts) for individual locations, the state with the greatest probability was adopted for each location (values for state with the highest probabilities ranged from 0.44 to 1).
We employed a foraging behavior termed “area-restricted search” (ARS; also employed by [54,55,56]) to characterize the foraging areas used by the bats (zones where they actively search for and consume fruit, then potentially rest; FAs). In the temporal sequences of behavioral states, one ARS behavior was defined as a sequence of foraging flight and stationary bout behavioral states in variable order (e.g., FF-SB-SB-FF-FF; see Fig. 1) or as a single foraging flight or stationary bout when no sequence occurred; two ARS behaviors being separated by at least one commuting flights.
As a consequence, each ARS behavior implied one FA for which the habitat type was determined (i.e. agricultural land, plantation, tree vegetation, or residential area) using available imagery for each location in Google Earth ver. 7.1 (image resolution was sufficient to achieve this for every location ). When different habitat types existed within a given FA (13.2% of FAs), the habitat most represented in terms of the number of locations was selected. Minimum convex polygons (MCPs) were created for each FA in QGIS ver. 3.4.1 (Fig. 1), although this was not possible for FAs where less than three locations were recorded (21% of FAs). In addition, the distance (Euclidean, in km) of each FA to the roost site (starting location of a bat on a given night) was calculated by averaging the distance of all FA locations to this site. The duration spent by a bat within a given FA was calculated (in minutes) using the time elapsed between all consecutive locations constituting the ARS behavior.
Finally, we identified FAs that were used several times by a given bat during the study (hereafter referred to as “re-visitation” behavior). To this end, we computed the mean distance between all pairs of locations constituting each minimum convex polygon (MCP; range: 1.3–530.3 m depending on the MCP). The median (37.0 m) of these averaged distances was used to create a buffer (radius: 18.5 m) around each MCP for FAs including at least three locations and around each location for FAs including less than three locations. These buffers were then used to identify overlapping FAs and thus to determine FAs that were revisited (Fig. 1).
We characterized foraging localities, nightly foraging activity and foraging-area re-visitation for bats. Overall, the effect of habitat type was tested as a proxy for the global quality of habitats used for foraging, mainly in terms of their food resources (i.e. typical diversity, abundance, renewal rate, availability, and energy content). In this context, the food-resource quality of residential areas was regarded as particularly high given the diversity (see the main species previously cited in the Methods - Study region and population; also see ) and abundance of fruit trees (with at least one tree in 93.5% of households, and up to 200 of the same species in a single backyard: ), plus the fact that fruits were frequently unharvested . However, sample sizes were relatively low for agricultural lands (only 17 ARS behaviors occurred in this habitat type) compared to other habitats (74, 84, 159 ARS behaviors in tree vegetation, plantations and residential areas, respectively). As such, for statistical issues, agricultural lands were discarded from our analyses of the effect of the habitat type (details of these analyses and the number of ARS behaviors taken into account are given below).
We tested whether the number of ARS behaviors displayed by a bat during a night differed between habitat types (n = 317 ARS behaviors). For this purpose, we used a Generalized Linear Mixed Model (GLMM; Poisson error distribution and log link function) which included habitat type as a fixed effect and the identity of individuals as an additive random effect to account for repeated measurements (several nights per individual; random effect also used in all subsequent mixed models). In addition, we used a GLMM (Gamma error distribution and log link function) to test whether the distance between the FA associated with each ARS behavior (n = 317 ARS behaviors) and the roost site was dependent on the habitat type of the FA (explanatory variable).
Nightly foraging activity
We compared several GLMMs (Poisson error distribution and log link function) with different combinations of explanatory variables to investigate variation in the number of ARS behaviors displayed by a bat at night (n = 334 ARS behaviors). Candidate explanatory variables included the total duration devoted to commuting flights on a given night (as a proxy for the energy costs related to overall nightly activity), and the total duration of ARS behavior in residential areas per night (as a proxy for the degree of utilization of a supposed high-quality habitat type). Following this, we explored variables influencing the duration of ARS behavior (n = 317 ARS behaviors) by comparing several GLMMs (Gamma error distribution and log link function). More specifically, we tested whether habitat type influenced this duration and we also tested the influence of the total ARS-behavior duration prior to the ARS behavior considered, since the beginning of the night (as a proxy for the degree their energetic needs had been satisfied).
Short-term fidelity to foraging areas
We finally compared several GLMMs (Binomial error distribution and logit link function) to investigate variation in the probability of a bat revisiting a given FA at least once in the following 48 h (two nights). Because Nipah virus typically survives in the environment for a few hours  and as long as a few days under optimal conditions, this period was of specific interest in terms of local pathogen-accumulation and transmission risks (n = 179 ARS behaviors). The influence of the duration of ARS-behaviors within a FA on a given night (the accumulated duration of ARS-behaviors in the rare case where a bat visited several times this FA during the same night) was tested on that probability, as a proxy for the FA quality. In addition, we tested the effect of habitat type as a proxy for a more global and integrative food-resource quality related more to habitat type than the specific FA used. We also tested the influence of the distance between the FA and the roost site (notably as a proxy for the energy costs required to reach this area).
All statistical analyses were carried out using R version 3.5.1 (R Core Team). GLMMs were fitted with the ‘lme4’ package. Type-II tests were performed to assess the significance of analyses evaluating the effect of a unique variable (‘car’ package). Where the influence of several candidate variables was investigated, all models containing one simple effect were computed. Models with additive effects were also computed, but interaction effects were not tested as meaningful hypotheses were lacking for these. We also computed a null (intercept-only) model. Candidate models were compared using the Akaike information criterion corrected for small sample sizes (AICc: ; ‘AICcmodavg’ package). A model was considered more competitive when its AICc was at least 2 units lower (ΔAICc) than others. Otherwise, we considered the AICc weight (ωAICc) as a measure of the probability of a model being the best model , and thus retained the model with the highest AICc weight. Estimated marginal means (EMMs, ± SE) were used to provide adjusted estimates of the response variable (‘emmeans’ package) by taking the simple effect of the other variable into account (in the case of additive effects), and also to perform post-hoc comparisons (Tukey’s comparisons). Overall, the variance explained by the effects retained was considered for each GLMM by computing the marginal and conditional R2 (see ; ‘MuMIn’ package). Quantitative explanatory variables were centered and scaled. Data on the number of ARS and FA re-visitations performed by each bat are provided in the supplementary material (Table S1 and S2).
The number of ARS behaviors displayed by a bat during the night differed between habitat types (χ2 = 39.5, df = 2, p < 0.001, marginal R2 = 0.13, conditional R2 = 0.23; Fig. 2): it was significantly higher in residential areas (2.2 ± 0.29) compared to other habitats (tree vegetation: 1.0 ± 0.16; fruit plantations: 1.2 ± 0.18). The distance between the foraging area associated with each ARS behavior and the roost site was dependent on the habitat type of the foraging area (χ2 = 15.8, df = 2, p < 0.001, marginal R2 = 0.07, conditional R2 = 0.25; Fig. 3): this distance was significantly higher for the residential areas (12.5 ± 2.4) compared to fruit plantations (7.0 ± 1.5), whereas it was intermediate for tree vegetation (10.1 ± 2.2).
Nightly foraging activity
The number of ARS behaviors displayed by a bat increased significantly with the total duration spent in commuting flights and decreased with the total duration spent in residential areas (marginal R2 = 0.32, conditional R2 = 0.41; Table 1; Fig. 4a and b). The durations of ARS behavior significantly decreased with the total ARS duration spent before the ARS evaluated (marginal R2 = 0.19, conditional R2 = 0.25; Table 1; Fig. 5).
Short-term fidelity to foraging areas
The probability of a bat revisiting a given foraging area for two consecutive nights significantly increased with the duration previously spent in the FA, decreased with the distance between the foraging area and the roost site, and was influenced by habitat type (marginal R2 = 0.57, conditional R2 = 0.62; Table 1; Fig. 6a, b and c): this probability was significantly higher for fruit plantations (0.71 ± 0.11), lowest for tree-vegetation (0.29 ± 0.11), and intermediate for residential areas (0.51 ± 0.11).
During the night, a bat performed more foraging bouts in residential areas on average (Fig. 2). A previous study on the same population showed that residential areas were the preferred foraging habitat of the species, but did not detail the use of foraging localities by individual bats on a given night (e.g., number and duration of visits) . In highlighting that individual bats favor residential areas by multiplying their foraging bouts in this habitat type, our work sheds further light on the distribution patterns revealed by . Foraging habitats among fruit bats vary from primary forest to urban areas (e.g., in the grey-headed flying fox Pteropus poliocephalus: ; in the Indian flying fox Pteropus giganteus: ) and may change over time (in the African straw-coloured fruit bat Eidolon helvum: ; in the solitary flying fox Pteropus dasymallus: ). Habitat choice in animals is influenced by environmental constraints and food resources, which may vary seasonally [62, 63]. Given the high density of backyards in residential areas that may provide abundant and diversified food resources , and the documented tolerance of some bat species to anthropization , our results suggest that P. lylei is a flexible generalist and opportunistic forager during the birthing season.
The residential areas chosen by bats were more distant from their roost site compared to other habitats (Fig. 3). This may be explained in an optimal-foraging framework, whereby the energetic costs associated with increased travel distances are offset by the benefits allied with good feeding conditions . In addition, our data showed that bats foraging longer in residential areas performed fewer foraging bouts during the night (Fig. 4b), and that they sometimes performed only one foraging behavior over an entire night. As the latter has been reported for other pteropodid bats in urban or semi-urban environments , these findings support the idea that residential areas may represent high-quality foraging localities for certain species, and that their resources may entirely support their daily energetic needs.
Nightly foraging activity
Our results indicate that the time spent foraging by our study bats in consecutive foraging areas decreased overnight. Wild animals need information about resource availability to optimize their foraging activities , and due to the renewal and maturation dynamics of fruits, frugivorous species are presented with predictable food resources . Consequently, we hypothesize that bats foraged for a longer period at the beginning of the night in areas with higher-quality food resources to ensure sufficient energy gains, before shifting to shorter foraging bouts to track the availability of food resources in other locations  and/or to complement their diets [71, 72]. Inspection of fruit resources has been reported in mangabey monkeys, which use fruiting synchrony to control specific trees and improve their foraging performance . As such, the positive correlation between the number of foraging bouts and the total time spent in commuting flight over a night (Fig. 4a) may indicate that some nights were more specifically devoted to collecting information about potential foraging areas, hence supporting and further specifying our hypothesis. In addition and as discussed above, the number of foraging bouts each night was also influenced by foraging-habitat localities.
Habitat type was not retained as a predictor of foraging-behavior duration in our analyses. The marginal theorem value predicts that animals remain longer in a profitable patch  and is supported by several field studies on mammals [15, 75]. It is possible that residential areas represent particularly attractive foraging habitats and that a relatively important number of foraging bouts displayed in these were particularly short and used to acquire information for future nights. This might at least partly explain the lack of habitat effect on the foraging-behavior duration in our study.
Short-term fidelity to foraging areas
Our bats revisited some foraging areas over short-term scales (in the following 48 h) according to their previous experience and environmental variables. Re-visitation behavior is widespread among animals, and may be determined by environmental heterogeneity, food-resource renewal rates and specific cognitive capabilities [76, 77]. Nonetheless, while foraging-site fidelity has been shown for many bat species [78,79,80], the factors influencing this behavior have received little attention. A recent study carried out on a nectar-feeding bat showed that the re-visitation of the most profitable places (i.e. in terms of quality and distance to the colony) was a consequence of a reinforcement learning strategy whereby bats avoided competition and maximized their global food intake .
We found the probability of a bat revisiting an area increased with the duration of the foraging behavior spent in the area (Fig. 6a). Similar positive links between short-term returns and previous visit durations (as a proxy for area quality) has also been found in other mammals and birds [82, 83]. This behavior is particularly relevant for species that do not deplete a localized food-resource in one visit or that forage in areas with high food-renewal rates (e.g., ), whereas longer intervals between visits may be expected in other cases (e.g., ). As our bats used different foraging areas, short-term returns to these may allow them to optimize their energy gains by exploiting high-quality resources in some familiar areas over consecutive nights. The negative association between re-visitation probability and distance to the roost site supports this hypothesis (Fig. 6b). Foraging areas in the vicinity of the roost site may be more familiar and visited more often (i.e. energy costs to visit them being lower). Exploration of individual GPS data (see Table S2) supports these explanations: areas revisited were relatively small in number, mainly located in the immediate vicinity of the roost site, and accumulated several visits during the study period. More generally, the influence and benefits of site familiarity on re-visitation behavior of bats has also been described for other mammals and birds during the breeding season [86,87,88].
In addition, foraging areas in fruit plantations were revisited more often by our bats, followed by residential areas (Fig. 6c). Habitat characteristics are known to influence re-visitation behavior in mammals [85, 89]. Even if plantations are especially attractive for plant-visiting bats , most commercially-grown fruits are harvested too early for flying-fox consumption . Only damaged or forgotten fruits reach sufficient maturity, which limits their availability and predictability for bats. One biological interpretation for our result could be the existence of specific features in some plantations (e.g., abandoned sites and attractive fruiting trees along plantation edges), resulting in high re-visitation rates for a few familiar areas (this hypothesis is partly supported by our GPS data; see Table S2). In contrast, given the hypothesized high foraging quality of residential areas, we suggest that the wide range of fruit trees (regarding species and fruiting phenology) and associated food resources has led to a heterogeneous degree of short-term attraction between the residential areas visited. Finally, the lowest probability found for tree-vegetation habitats might be attributable to their possessing lower diversity and quantity of food resources during our study period, but we cannot evaluate this possibility as the relevant information is currently lacking. Overall, our results suggest that the reinforcement learning strategy previously mentioned  may have led to the re-visitation pattern found in our work, but quantitative information about food quality and quantity in each foraging locality are needed to support this hypothesis.
Limitations and future prospects
Contrasting degrees of inter-individual variability regarding the behavior studied is suggested by individual data exploration (see Table S1 and S2), explained variance (both marginal and conditional R2) and confidence intervals for all results. Inter-individual variation in foraging could result from differences in personality traits, previous experience of individuals (e.g., familiarity with food distribution and predictability) and degrees of specialization . In association with the solitary nature of foraging fruit bats , we suggest that the heterogeneity of our study region (e.g., habitat types, food-resource diversity and availability) promoted different individual foraging strategies. Nevertheless, our small sample size and limited sampling period preclude conclusions regarding the prevalence of these behaviors in the wider population and their seasonal consistency. Investigation of other intrinsic variables (e.g., age, sex and body condition) and long-term monitoring of individual bats (using Argos telemetry) would likely provide some interesting perspectives.
Our study provides new insights into the foraging ecology of Pteropus lylei whose patterns of re-visitation behavior appear to be relatively complex. These patterns were significant but influenced by several variables and were likely constrained to relatively few areas at individual levels. Furthermore, the prevalence of Nipah virus in this population is 0.9 and 0.2% during the pregnancy and lactation, respectively (i.e. the reproductive periods encompassed by the study) . Since the return of an infected bat to a given area over consecutive nights could increase probability of pathogen accumulation there, our findings suggest that risks of transmission from infected bats could be highly localized, at least during the short term. More generally, our work emphasizes the importance of considering individual-movement patterns in epidemiological frameworks. However, the complexity of foraging patterns revealed by our study (correlates and inter-individual variation) also means investigation of transmission pathways and prevention of spillover risks will be challenging, particularly given the ongoing expansion of human-wildlife interfaces in Southeast Asia.
Our study employed a path-segmentation method (hidden Markov model) to elucidate the behavioral sequences of individual P. lylei from GPS data. Our results suggest that environmental and individual features (e.g., availability and quality of food resources, habitat-familiarity and experience of bats) have a significant bearing on their patterns of foraging behavior and fidelity to a given foraging area over consecutive nights. As P. lylei is a natural reservoir for Nipah virus, this is particularly relevant to viral transmission risks and indicates that fine-scale ecological studies of species capable of adapting to newly anthropized environments will be important to understand the ecological processes that maintain biodiversity in these landscapes and assess the potential for pathogen transmission in human-wildlife interfaces.
Availability of data and materials
Akaike Information Criterion
Estimated Marginal Mean
Generalized Linear Mixed Model
Gobal Positioning System
Hidden Markov Model
Minimum Convex Polygon
Pyke G. Optimal foraging theory: a critical review. Annu Rev Ecol Syst. 1984;15:523–75.
Alcock J. Animal behavior: an evolutionary approach. 10th ed. Sunderland: Sinauer Associates; 2013.
Colwell MA. Shorebird ecology, conservation, and management. Oakland: University of California Press; 2010.
Schoener TW. Theory of feeding strategies. Annu Rev Ecol Syst. 1971;2:369–404.
Azorit C, Tellado S, Oya A, Moro J. Seasonal and specific diet variations in sympatric red and fallow deer of southern Spain: a preliminary approach to feeding behaviour. Anim Prod Sci. 2012;52:720–7.
Ceia FR, Ramos JA. Individual specialization in the foraging and feeding strategies of seabirds: a review. Mar Biol. 2015;162:1923–38.
Pethybridge H, Bodin N, Arsenault-Pernet E, Bourdeix J, Brisset B, Bigot J, et al. Temporal and inter-specific variations in forage fish feeding conditions in the NW Mediterranean: lipid content and fatty acid compositional changes. Mar Ecol Prog Ser. 2014;512:39–54.
Baylis AMM, Page B, McKenzie J, Goldsworthy SD. Individual foraging site fidelity in lactating New Zealand fur seals: continental shelf vs. oceanic habitats. Mar Mammal Sci. 2012;28:276–94.
Beverly BD, McLendon H, Nacu S, Holmes S, Gordon DM. How site fidelity leads to individual differences in the foraging activity of harvester ants. Behav Ecol. 2009;20:633–8.
Bécares J, García-Tarrasón M, Villero D, Bateman S, Jover L, García-Matarranz V, et al. Modelling terrestrial and marine foraging habitats in breeding Audouin’s gulls Larus audouinii: timing matters. PLoS One. 2015;10:e0120799.
Kato A, Ropert-Coudert Y, Chiaradia A. Regulation of trip duration by an inshore forager, the little penguin (Eudyptula Minor), During Incubation. The Auk. 2008;125:588–93.
Lobo N, Green DJ, Millar JS. Effects of seed quality and abundance on the foraging behavior of deer mice. J Mammal. 2013;94:1449–59.
Robinson SK, Holmes RT. Foraging behavior of Forest birds: the relationships among search tactics, diet, and habitat structure. Ecology. 1982;63:1918–31.
Larter NC, Gates CC. Diet and habitat selection of wood bison in relation to seasonal changes in forage quantity and quality. Can J Zool. 1991;69:2677–85.
Watts DP. Strategies of habitat use by mountain gorillas. Folia Primatol. 1991;56:1–16.
Arthur B, Hindell M, Bester M, Trathan P, Jonsen I, Staniland I, et al. Return customers: foraging site Fidelity and the effect of environmental variability in wide-ranging Antarctic Fur seals. PLoS One. 2015;10:e0120888.
Carroll G, Harcourt R, Pitcher BJ, Slip D, Jonsen I. Recent prey capture experience and dynamic habitat quality mediate short-term foraging site fidelity in a seabird. Proc R Soc B Biol Sci. 2018;285:20180788.
Evens R, Beenaerts N, Neyens T, Witters N, Smeets K, Artois T. Proximity of breeding and foraging areas affects foraging effort of a crepuscular, insectivorous bird. Sci Rep. 2018;8:3008.
Jones TB, Evans JC, Morand-Ferron J. Urbanization and the temporal patterns of social networks and group foraging behaviors. Ecol Evol. 2019;9:4589–602.
Gehlbach FR. Eastern screech-owl responses to suburban sprawl, warmer climate, and additional avian food in Central Texas. Wilson J Ornithol. 2012;124:630–3.
Van Houtan KS, Pimm SL, Bierregaard RO, Lovejoy TE, Stouffer PC. Local extinctions in flocking birds in Amazonian forest fragments. Evol Ecol Res. 2006;8:129–48.
Lindenmayer DB, Fischer J. Habitat fragmentation and landscape change: an ecological and conservation synthesis. Washington: Island Press; 2006.
Fahrig L. Non-optimal animal movement in human-altered landscapes. Funct Ecol. 2007;21:1003–15.
Greggor AL, Berger-Tal O, Blumstein DT, Angeloni L, Bessa-Gomes C, Blackwell BF, et al. Research priorities from animal behaviour for Maximising conservation Progress. Trends Ecol Evol. 2016;31:953–64.
Holyoak M, Casagrandi R, Nathan R, Revilla E, Spiegel O. Trends and missing parts in the study of movement ecology. Proc Natl Acad Sci. 2008;105:19060–5.
Voigt CC, Kingston T. Bats in the Anthropocene: conservation of bats in a changing world. Springer Open: Cham; 2016.
Acharya PR, Racey PA, Sotthibandhu S, Bumrungsri S. Feeding behaviour of the dawn bat (Eonycteris spelaea) promotes cross pollination of economically important plants in Southeast Asia. J Pollinat Ecol. 2015;15:44–50.
Boyles JG, Sole CL, Cryan PM, McCracken GF. On estimating the economic value of insectivorous bats: prospects and priorities for biologists. In: Adams RA, Pedersen SC, editors. Bat Evol Ecol Conserv [internet]. New York: Springer New York; 2013. p. 501–15. Available from:. https://doi.org/10.1007/978-1-4614-7397-8_24.
Fleming TH, Geiselman C, Kress WJ. The evolution of bat pollination: a phylogenetic perspective. Ann Bot. 2009;104:1017–43.
Kingston T. Research priorities for bat conservation in Southeast Asia: a consensus approach. Biodivers Conserv. 2010;19:471–84.
Lane DJW, Kingston T, Lee BPY-H. Dramatic decline in bat species richness in Singapore, with implications for Southeast Asia. Biol Conserv. 2006;131:584–93.
Calisher CH, Childs JE, Field HE, Holmes KV, Schountz T. Bats: important reservoir hosts of emerging viruses. Clin Microbiol Rev. 2006;19:531–45.
Plowright RK, Eby P, Hudson PJ, Smith IL, Westcott D, Bryden WL, et al. Ecological dynamics of emerging bat virus spillover. Proc R Soc B Biol Sci. 2015;282:20142124.
Fleming TH, Eby P. Ecology of bat migration. In: Kunz TH, Fenton MB, editors. Bat Ecol. Chicago: University of Chicago Press; 2003. p. 156–208.
Hebblewhite M, Haydon DT. Distinguishing technology from biology: a critical review of the use of GPS telemetry data in ecology. Philos Trans R Soc B Biol Sci. 2010;365:2303–12.
Tsang SMS. Phylogeography of Southeast Asian flying foxes (Chiroptera: Pteropodidae: Pteropus), PhD-thesis. New York: City University of New York; 2015. Available from: https://academicworks.cuny.edu/gc_etds/1161.
Bellini WJ, Harcourt BH, Bowden N, Rota PA. Nipah virus: an emergent paramyxovirus causing severe encephalitis in humans. J Neuro-Oncol. 2005;11:481–7.
Kessler MK, Becker DJ, Peel AJ, Justice NV, Lunn T, Crowley DE, et al. Changing resource landscapes and spillover of henipaviruses: Resource landscapes and henipavirus spillover. Ann N Y Acad Sci. 2018; Cited 2018 Sep 20. https://doi.org/10.1111/nyas.13910.
Oleksy R, Racey PA, Jones G. High-resolution GPS tracking reveals habitat selection and the potential for long-distance seed dispersal by Madagascan flying foxes Pteropus rufus. Glob Ecol Conserv. 2015;3:678–92.
Mahmood-Ul-Hassan M, Gulraiz TL, Rana SA, Javid A. The diet of Indian flying-foxes (Pteropus giganteus) in urban habitats of Pakistan. Acta Chiropterologica. 2010;12:341–7.
Nakamoto A, Kinjo K, Izawa M. Dietary plasticity in the Ryukyu flying fox on a Subtropical Island at the northern range limit of Pteropus. Acta Chiropterologica. 2015;17:105–16.
Choden K, Ravon S, Epstein JH, Hoem T, Furey N, Gely M, et al. Pteropus lylei primarily forages in residential areas in Kandal, Cambodia. Ecol Evol. 2019;9:4181–91.
Weber N, Duengkae P, Fahr J, Dechmann DKN, Phengsakul P, Khumbucha W, et al. High-resolution GPS tracking of Lyle’s flying fox between temples and orchards in Central Thailand: GPS tracking of flying foxes. J Wildl Manag. 2015;79:957–68.
Edelhoff H, Signer J, Balkenhol N. Path segmentation for beginners: an overview of current methods for detecting changes in animal movement patterns. Mov Ecol. 2016;4:21.
Martin J, Tolon V, Moorter B, Basille M, Calenge C. On the use of telemetry in habitat selection studies; 2009. p. 37–55.
Ravon S, Furey NM, Vibol H, Cappelle J. A rapid assessment of flying fox (Pteropus spp.) colonies in Cambodia. Cambodian J Nat Hist. 2014;2014:14–8.
Cappelle J, Hoem T, Hul V, Furey N, Nguon K, Prigent S, et al. Nipah virus circulation at human–bat interfaces, Cambodia. Bull World Health Organ. 2020;98:539–47.
Wacharapluesadee S, Boongird K, Wanghongsa S, Ratanasetyuth N, Supavonwong P, Saengsen D, et al. A longitudinal study of the prevalence of Nipah virus in Pteropus lylei bats in Thailand: evidence for seasonal preference in disease transmission. Vector-Borne Zoonotic Dis. 2010;10:183–90.
Hengjan Y, Saputra V, Mirsageri M, Pramono D, Kasmono S, Basri C, et al. Nighttime behavioral study of flying foxes on the southern coast of West Java, Indonesia. J Vet Med Sci. 2018;80:1146–52.
Johnson DS, London JM, Lea M-A, Durban JW. Continuous-time correlated random walk model for animal telemetry data. Ecology. 2008;89:1208–15.
Michelot T, Langrock R, Patterson TA. moveHMM: an R package for the statistical modelling of animal movement data using hidden Markov models. Methods Ecol Evol. 2016;7:1308–15.
Roeleke M, Teige T, Hoffmeister U, Klingler F, Voigt CC. Aerial-hawking bats adjust their use of space to the lunar cycle. Mov Ecol. 2018;6:11.
Hurme E, Gurarie E, Greif S, LG HM, Flores-Martínez JJ, Wilkinson GS, et al. Acoustic evaluation of behavioral states predicted from GPS tracking: a case study of a marine fishing bat. Mov Ecol. 2019;7:21.
Kareiva P, Odell G. Swarms of predators exhibit “Preytaxis” if individual predators use area-restricted search. Am Nat. 1987;130:233–70.
Voigt CC, Frick WF, Holderied MW, Holland R, Kerth G, Mello MAR, et al. Principles and patterns of bat movements: from aerodynamics to ecology. Q Rev Biol. 2017;92:267–87.
Weimerskirch H. Are seabirds foraging for unpredictable resources? Deep Sea Res Part II Top Stud Oceanogr. 2007;54:211–23.
Fogarty R, Halpin K, Hyatt AD, Daszak P, Mungall BA. Henipavirus susceptibility to environmental variables. Virus Res. 2008;132:140–4.
Symonds MRE, Moussalli A. A brief guide to model selection, multimodel inference and model averaging in behavioural ecology using Akaike’s information criterion. Behav Ecol Sociobiol. 2011;65:13–21.
Nakagawa S, Schielzeth H. A general and simple method for obtaining R 2 from generalized linear mixed-effects models. Methods Ecol Evol. 2013;4:133–42.
Giles JR, Eby P, Parry H, Peel AJ, Plowright RK, Westcott DA, et al. Environmental drivers of spatiotemporal foraging intensity in fruit bats and implications for Hendra virus ecology. Sci Rep. 2018;8 Available from: http://www.nature.com/articles/s41598-018-27859-3. Cited 2019 Aug 26.
Mishra R, Dookia S, Bhattacharya P. Avenue plantations as biodiversity havens: a case study of population status of the Indian flying fox, Pteropus giganteus Brunnich, 1782 and implications for its conservation in the urban megacity, Delhi, India. Proc Zool Soc. 2020;73:127–36.
Fahr J, Abedi-Lartey M, Esch T, Machwitz M, Suu-Ire R, Wikelski M, et al. Pronounced seasonal changes in the movement ecology of a highly gregarious central-place forager, the African straw-Coloured fruit bat (Eidolon helvum). PLoS One. 2015;10:e0138985.
Nakamoto A, Kinjo K. Ranging patterns and habitat use of a solitary flying fox (Pteropus dasymallus) on Okinawa-Jima Island, Japan. Acta Chiropterologica. 2012;14:387.
Moore JH, Sittimongkol S, Campos-Arceiz A, Sumpah T, Eichhorn MP. Fruit gardens enhance mammal diversity and biomass in a southeast Asian rainforest. Biol Conserv. 2016;194:132–8.
Russo D, Ancillotto L. Sensitivity of bats to urbanization: a review. Mamm Biol. 2015;80:205–12.
Kacelnik A, Houston AI. Some effects of energy costs on foraging strategies. Anim Behav. 1984;32:609–14.
Acharya PR, Racey PA, Sotthibandhu S, Bumrungsri S. Home-range and foraging areas of the Dawn bat Eonycteris spelaea in agricultural areas of Thailand. Acta Chiropterologica. 2015;17:307–19.
Kamil A. The ecology of foraging behavior: implications for animal learning and memory. Annu Rev Psychol. 1985;36:141–69.
Riotte-Lambert L, Matthiopoulos J. Environmental predictability as a cause and consequence of animal movement. Trends Ecol Evol. 2020;35:163–74.
Bateson M. Interval Timing and Optimal Foraging. In: Meck W, editor. Functional and Neural Mechanisms of Interval Timing: CRC Press; 2003. p. 113–41. Available from:. https://doi.org/10.1201/9780203009574.ch5.
Rapport DJ. Optimal foraging for complementary resources. Am Nat. 1980;116:324–46.
Sanchez F, Kotler BP, Korine C, Pinshow B. Sugars are complementary resources to ethanol in foods consumed by Egyptian fruit bats. J Exp Biol. 2008;211:1475–81.
Janmaat KRL, Chapman CA, Meijer R, Zuberbühler K. The use of fruiting synchrony by foraging mangabey monkeys: a ‘simple tool’ to find fruit. Anim Cogn. 2012;15:83–96.
Charnov EL. Optimal foraging, the marginal value theorem. Theor Popul Biol. 1976;9:129–36.
Austin D, Bowen WD, McMillan JI, Iverson SJ. Linking movement, diving, and habitat to foraging success in a large marine predator. Ecology. 2006;87:3095–108.
Berger-Tal O, Bar-David S. Recursive movement patterns: review and synthesis across species. Ecosphere. 2015;6:149.
Winter Y. Foraging in a complex naturalistic environment: capacity of spatial working memory in flower bats. J Exp Biol. 2005;208:539–48.
Egert-Berg K, Hurme ER, Greif S, Goldstein A, Harten L, Herrera M. LG, et al. resource ephemerality drives social foraging in bats. Curr Biol. 2018;28:3667–73.
Field HE, Smith CS, de Jong CE, Melville D, Broos A, Kung N, et al. Landscape utilisation, animal behaviour and Hendra virus risk. EcoHealth. 2016;13:26–38.
Kapfer G, Rigot T, Holsbeek L, Aron S. Roost and hunting site fidelity of female and juvenile Daubenton’s bat Myotis daubentonii (Kuhl, 1817) (Chiroptera: Vespertilionidae). Mamm Biol. 2008;73:267–75.
Goldshtein A, Handel M, Eitan O, Bonstein A, Shaler T, Collet S, et al. Reinforcement learning enables resource partitioning in foraging bats. Curr Biol. 2020;30:1–7.
Bracis C, Bildstein KL, Mueller T. Revisitation analysis uncovers spatio-temporal patterns in animal movement data. Ecography. 2018;41:1801–11.
English M, Ancrenaz M, Gillespie G, Goossens B, Nathan S, Linklater W. Foraging site recursion by forest elephants Elephas maximus borneensis. Curr Zool. 2014;60:551–9.
Seidel DP, Linklater WL, Kilian W, du Preez P, Getz WM. Mesoscale movement and recursion behaviors of Namibian black rhinos. Mov Ecol. 2019;7. Cited 2020 Aug 10. https://doi.org/10.1186/s40462-019-0176-2.
Watts DP. Long-term habitat use by mountain gorillas (Gorilla gorilla beringei). 2. Reuse of foraging areas in relation to resource abundance, quality, and depletion. Int J Primatol. 1998;19:681–702.
Brown CR, Brown MB, Brazeal KR. Familiarity with breeding habitat improves daily survival in colonial cliff swallows. Anim Behav. 2008;76:1201–10.
Greenwood PJ. Mating systems, philopatry and dispersal in birds and mammals. Anim Behav. 1980;28:1140–62.
Yoder JM. The cost of dispersal: predation as a function of movement and site familiarity in ruffed grouse. Behav Ecol. 2004;15:469–76.
De Knegt HJ. Beyond the here and now: herbivore ecology in a spatial-temporal context. Wageningen: Wageningen University; 2010. PhD-thesis. Available from: http://edepot.wur.nl/137139.
Madrid-López SM, Galindo-González J, Castro-Luna AA. Mango orchards and their importance in maintaining phyllostomid bat assemblages in a heterogeneous landscape. Acta Chiropterologica. 2019;21:375–83.
Fujita MS. Flying foxes (Chiroptera: Pteropodidae): threatened animals of key ecological and economic importance. Conserv Biol. 1991;5:455–63.
Toscano BJ, Gownaris NJ, Heerhartz SM, Monaco CJ. Personality, foraging behavior and specialization: integrating behavioral and food web ecology at the individual level. Oecologia. 2016;182:55–69.
Sikes RS, Gannon WL. Guidelines of the American Society of Mammalogists for the use of wild mammals in research. J Mammal. 2011;92:235–53.
Sikes RS, the Animal Care and Use Committee of the American Society of Mammalogists. 2016 Guidelines of the American Society of Mammalogists for the use of wild mammals in research and education. J Mammal. 2016;97:663–88.
We thank Pring Long and Yen Sroy for their assistance in the field. Particular gratitude is due to Dr. Keo Omaliss of the Cambodian Forestry Administration for arranging study approval and permissions. We also thank three anonymous referees for their constructive comments on the manuscript.
This study was supported by the European Commission (ComAcross project, grant no. DCI-ASIE/2013/315–047; EBOSURSY project, Ph.D. grant to ES), the French Centre National d’Etudes Spatiales (TeleNipah project, DAR 4800000780) and the USAID’s Emerging Pandemic programme (PREDICT project).
Ethics approval and consent to participate
We captured and handled bats in the field in accordance with guidelines approved by the American Society of Mammalogists [93, 94], in addition to the requirements of the statutory study permission provided by the national authority responsible for wildlife research: the Forestry Administration of the Cambodian Ministry of Agriculture, Forest and Fisheries. As animal ethics committees did not exist at the time of our study in Cambodia, all study aspects were overseen by the Forestry Administration who participated in the field investigations.
Consent for publication
The authors declare that they have no competing interest.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Foraging tracks of 8 individuals of Pteropus lylei equipped with GPS loggers, in Cambodia.
Distribution of step lengths and turning angles related to the GPS data.
Mean number (± SD) of area restricted search (ARS) behaviors displayed by each bat in relation to habitat type.
Number of foraging-area (FA) re-visitation behaviors displayed by each individual, number of FAs that were used several times by a bat during the study period according to the habitat type, and distance (m) range between FAs used several times and the roost site.
About this article
Cite this article
Schloesing, E., Chambon, R., Tran, A. et al. Patterns of foraging activity and fidelity in a southeast Asian flying fox. Mov Ecol 8, 46 (2020). https://doi.org/10.1186/s40462-020-00232-8
- GPS telemetry
- Hidden Markov models
- Movement ecology
- Nipah virus