https://orcid.org/0000-0001-6375-6954
Figures
Abstract
Knowing the influence of intrinsic and environmental traits on animals’ movement is a central interest of ecology and can aid to enhance management decisions. The giant anteater (Myrmecophaga tridactyla) is a vulnerable mammal that presents low capacity for physiological thermoregulation and uses forests as thermal shelters. Here, we aim to provide reliable estimates of giant anteaters’ movement patterns and home range size, as well as untangle the role of intrinsic and environmental drivers on their movement. We GPS-tracked 19 giant anteaters in Brazilian savannah. We used a continuous-time movement model to estimate their movement patterns (described by home range crossing time, daily distance moved and directionality), and provide an autocorrelated kernel density estimate of home range size. Then, we used mixed structural equations to integratively model the effects of sex, body mass and proportion of forest cover on movement patterns and home range size, considering the complex net of interactions between these variables. Male giant anteaters presented more intensive space use and larger home range than females with similar body mass, as it is expected in polygynous social mating systems. Males and females increased home range size with increasing body mass, but the allometric scaling of intensity of space use was negative for males and positive for females, indicating different strategies in search for resources. With decreasing proportion of forest cover inside their home ranges, and, consequently, decreasing thermal quality of their habitat, giant anteaters increased home range size, possibly to maximize the chances of accessing thermal shelters. As frequency and intensity of extreme weather events and deforestation are increasing, effective management efforts need to consider the role of forests as an important thermal resource driving spatial requirements of this species. We highlight that both intrinsic and environmental drivers of animal movement should be integrated to better guide management strategies.
Citation: Giroux A, Ortega Z, Oliveira-Santos LGR, Attias N, Bertassoni A, Desbiez ALJ (2021) Sexual, allometric and forest cover effects on giant anteaters’ movement ecology. PLoS ONE 16(8): e0253345. https://doi.org/10.1371/journal.pone.0253345
Editor: Bi-Song Yue, Sichuan University, CHINA
Received: April 8, 2021; Accepted: June 2, 2021; Published: August 18, 2021
Copyright: © 2021 Giroux et al. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.
Data Availability: Giant anteaters are listed as a vulnerable species across their distribution. Therefore, location data of this species are extremely sensitive. All the movement data used in analyses are available for visualization from the Movebank database (www.movebank.org; study names: Myrmecophaga tridactyla Pantanal, Myrmecophaga tridactyla SP). Data download requests may be sent to the Center for Species Survival Brazil (SSC) - International Union for Conservation of Nature (IUCN) (csebrasil.contato@gmail.com / fabiana@cpsg.org).
Funding: This study was financed in part by the Coordenação de Aperfeiçoamento de Pessoal de Nível Superior, Brasil (CAPES; https://www.gov.br/capes/pt-br) under finance code 001, PROEX 88887.360861/2019-00 to AG, and PNPD 1694744 to ZO. It was also supported by the Junta de Andalucía under Research Postdoctoral Contract 401 to ZO, and by the São Paulo Research Foundation (FAPESP; https://fapesp.br), under finance codes 2013/18526-9 to AB and 2013/04957-8 to AB. Part of the costs for the field work were covered by the Giant Armadillo Conservation Program, which benefited from multiple grants, mostly from North American and European zoos listed at www.icasconservation.org.br/partners. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.
Competing interests: The authors have declared that no competing interests exist.
Introduction
Animal movement is a key process of ecology, driving animals’ survival and fitness [1]. The individuals’ movement patterns shape their home range, which can be defined as the part of their cognitive map that they choose to be continuously updated [2, 3]. The home range should offer the needed conditions for basic activities of food gathering, mating, and caring for young [4]. Describing movement patterns and quantifying the home range size allow us to better understand the ecology and spatial requirements of animals to make appropriate management decisions that can help to preserve wildlife populations [5, 6]. Such knowledge has become even more important as human actions are increasingly endangering natural systems [7, 8]. Theoretical and empirical studies have more often focused on understanding ‘typical’ movement of a species than understanding its variation [9]. However, both movement patterns and home range size widely vary between individuals within a population, and these variations are commonly influenced by intrinsic and environmental traits [5, 7].
Sex and body mass are among the main intrinsic traits driving intraspecific variations on animal movement [10]. The mating system and the associated reproductive tactics employed by males and females within a species influence the evolutionary selection of various characters [11]. This can result in sex-related metabolic, cognitive, and behavioral differences that should be reflected on movement patterns and, consequently, on home range size [12]. Body mass, in turn, has direct influence on the individual’s body mechanics and physiology [13, 14]. Besides, body mass can drive movement patterns and home range size by influencing, for example, the energetic requirements of individuals [15, 16], the foraging experience [17], and/or the orientation ability [18]. Therefore, the animal movement allometry relative to individuals’ body mass is commonly found [13–18]. Among environmental traits shaping animal movement patterns and home range size, the proportion of available forest cover stands out to animals that habit fragmented landscapes and use forests as a resource [19–21].
The influence of intrinsic and environmental traits on animal movement is being increasingly better understood as technological advances on tracking methods increase [9]. Although the analysis of movement data is still challenging [22, 23], the recent implementation of continuous-time movement models on understanding movement patterns and on estimating home range size has allowed great advances [24, 25]. For high-quality GPS tracking data of range-resident individuals, these models allow the estimation of descriptors of movement patterns such as home range crossing-time, daily distance moved and directionality [26]. These descriptors bring insights on underlying movement processes determining home range and can be used to provide an autocorrelated kernel density estimator of home range size [26]. However, previous research has focused on evaluating one specific movement metric at a time [27, 28], disregarding the possible causal relationships of the descriptors of movement patterns with each other and with home range size. Both accurate estimates of animal movement, based on movement models, and integrative approaches that consider the complex network of relations between the variables can help us to understand the effect of intrinsic and environmental traits on movement patterns and home range size.
The giant anteater (Myrmecophaga tridactyla) is a vulnerable mammal whose movement patterns and home range size have been previously studied to better guide its management and conservation [29, 30]. Their original spatial distribution covered from Belize to the south of South America, excluding the Andes [29]. While some populations are already locally extinct, others are facing habitat loss, wildfires, roadkills, conflicts with dogs and other threats [29]. In this scenario, their low reproductive rate and long periods of parental care make giant anteaters conservation status even more worrisome [31, 32]. Despite being commonly associated with open habitats [33], forests have a fundamental role in giant anteater thermoregulation [34–36]. This is because giant anteaters present reduced body heat production [37] and low capacity of physiological thermoregulation, and forests act as important thermal shelters. Besides showing smaller environmental temperature variation than adjacent open areas, forest patches buffer rain and chilly winds and offer protection against solar radiation [38]. Therefore, it is also worrying that deforestation may be reducing the habitat thermal quality for these animals across their current distribution [39].
Despite the efforts to understand giant anteaters’ movement ecology, previous estimates of their home range size have ignored the intrinsic autocorrelation of high-resolution movement data and have not been based on movement models, probably generating underestimated results [24, 30]. While some studies showed no evidence of sexual effects on their movement [36, 40], other ones showed males presenting longer daily activity time [41, 42] and using larger areas than females [42]. Because of their sexual size dimorphism [36], the possible influence of body mass on movement needs to be considered when assessing sexual effects. Besides, although we know that giant anteaters select forests to set their home ranges and allocate time within it [42], we still ignore if the proportion of forest cover within home ranges influences their movement patterns and spatial requirements. Here we used a continuous-time movement model to offer reliable estimates of giant anteaters’ movement patterns (specifically home range crossing-time, daily distance moved and directionality) and home range size. Then, we investigated the effect of sex, body mass and proportion of forest cover on giant anteaters’ movement patterns and home range size. Using an integrative approach, we were able to uncover all these effects simultaneously, controlling for the possible relations among descriptors of movement pattern and with home range size.
Due to their probably polygynous social mating system [43], we expected male giant anteaters to increase their chances of mating opportunities by moving longer daily distances and using larger home ranges than females (Fig 1B and 1D) [44]. We also expected an allometric scaling between body mass and movement, since larger bodied individuals have higher energetic requirements than smaller ones [45]. Larger giant anteaters should increase the intensity of space use, increasing home range crossing-time and daily distance moved while decreasing directionality. This is because this increase in the intensity of space use should increase the individuals’ chances to find food resources–mainly ants and termites–spread on the landscape (Fig 1A–1C) [26, 45, 46]. Besides, it is reasonable to expect that larger animals will require more space to meet their energetic requirements [45], so they would also increase home range size with increasing body mass (Fig 1D). We expect that increasing the proportion of forest cover inside the home range will lead the animals to increase home range crossing time, decreasing daily distance moved and directionality, because the forests’ three-dimensional structure should present physical obstacles to displacement, imposing more friction than open grasslands (Fig 1E and 1G) [47]. Finally, lower proportions of forest inside home range would decrease the animal’s access to thermal shelters, decreasing the habitat thermal quality. This could lead animals to increase their spatial requirements, and, consequently, increase home range size (Fig 1H) [48].
Methods
Study site
We carried the study out in two savannah areas in the Brazilian territory: (1) Santa Barbara Ecological Station, São Paulo state (22° 48’ 59’’ S, 49° 14’ 12’’ W) and (2) Baía das Pedras Ranch, Mato Grosso do Sul state (19° 18’ 9" S, 55° 47’ 4" W). The study areas have a tropical climate, with rainy summers and dry winters [49]. The landscape of both studied areas is composed of mosaics of open grasslands, scrublands, savannahs, and woodlands [50, 51]. The landscape of Santa Barbara Ecological Station also includes anthropic elements, such as exotic forests of Pinus sp. and Eucalyptus sp., as well as highways [51]. Baía das Pedras Ranch is located within the Pantanal wetland, and it presents a naturally fragmented landscape with permanent and temporary salty and freshwater ponds, where open grasslands are subjected to seasonal flooding [50].
We classified the landscapes using georeferenced maps (LANDSAT 7 TM) and the MapBiomas database (Collection 5) [52]. To test the effect of the forest cover in movement patterns and home range size of giant anteaters, we summarized the observed habitats in two categories: forest or non-forest. Forest areas included woodland savannahs, woodlands, riparian forests, regenerating arboreal vegetation and exotic forests. Non-forest areas included open grasslands, scrublands, open savannahs, and areas without vegetation cover. We calculated the proportion of forest cover within each individual home range dividing the number of pixels classified as forest by the total number of pixels. We performed satellite image processing and supervised classifications using raster [53], maptools [54], and rgdal [55] packages available in the R environment [56].
Capture and data collection
We searched for giant anteaters by horse or by pickup vehicle at low speed (maximum of 20 km/h). Once we saw the anteaters, we captured them using dip nets, dart-guns, or a blowpipe. Anteaters were immobilized and sedated following the protocol described by [42] in Santa Barbara Ecological Station and following the protocol described by [57] in Baía das Pedras Ranch. Each captured individual was sexed, weighted, and equipped with a global positioning system (GPS) harness during anesthesia. We conducted a T test [58] to compare the mean body mass between individuals of our two study areas. None of the tracking devices exceeded 3% of the animals’ body mass. All procedures were conducted in accordance with the Guidelines of the American Society of Mammalogists for the use of wild mammals in research [59] and were performed under the license numbers SISBIO 16010–1 and SISBIO 38326–5 (Chico Mendes Institute for Biodiversity Conservation). After completing their recovery from the anesthesia, we released the giant anteaters at the site of capture for movement GPS-tracking.
Movement patterns and home range analysis
We described animals’ movement patterns and estimated home range size using the ctmm R package [24, 56, 60, 61]. We first examined the empirical variogram of each individual tracking data to check for an asymptote [26], as it is an evidence of range residence and a premise for the movement parameters estimation [61]. Because tracking data with such short sampling intervals are inherently autocorrelated, we also used the variogram to investigate the autocorrelation structure of data, obtaining starting values for the variance and autocorrelation timescales. Then, we fitted continuous-time movement models to the individuals’ location data via maximum likelihood. Among the fitted models, we included the Brownian motion model (BM), the Ornstein-Uhlenbeck model (OU), the Integrated OU model (IOU) and the Ornstein-Uhlenbeck-F model (OUF) [26]. We ranked the movement models based on the second order Akaike Information Criterion (AICC) [62] and selected the one with the best fit for each individual anteater data set.
For those animals that better fitted OUF model, we obtained the three descriptors of individuals’ movement patterns: home range crossing time (timescale of autocorrelation in position), daily distance moved and directionality (direction persistence timescale), as well as their confidence intervals [24]. For those animals that showed range residence (i.e., better fitted OU or OUF models), we used the ninety-five per cent area corrected autocorrelated kernel density estimator (AKDEc 95%) to estimate the individuals’ home range size and its confidence limits. AKDEc is a nonparametric home-range estimator that assumes the data represent a sample from a nonstationary, autocorrelated, continuous movement process [63]. This estimator allows movement models to be fitted to data with different temporal structures (e.g., irregular sampling regime, gaps, and short sampling time). Also, AKDEc allows to compare home ranges of individuals with different monitoring times. This is because AKDEc extrapolates the data, basing itself on parameters of the model selected for each individual data set, to provide reliable home range estimates [63].
Structural Equation Modeling
We used mixed Structural Equations Modeling (mixed-SEM) [64] to investigate: (1) the effect of intrinsic traits (individuals’ sex and body mass) in movement patterns and home range size and (2) the effect of an environmental trait (proportion of forest cover inside the individuals’ home range) in movement patterns and home range size. Because the descriptors of movement patterns can be related to each other, and can modulate home range size, we controlled for these possible relationships in an integrative approach (see Fig 2). In this approach, the same variable could simultaneously act as response in an equation and as predictor in another one (Fig 2) [64]. Mixed-SEM allowed us to disentangle a complex net of interactions, estimating the indirect, direct, and total effects among variables [65]. Indirect effects were estimated by the product of the direct effects that compose them, and total effects were given by the sum of direct and indirect effects [66, 67].
The standardized coefficients (β) represent the relative strength of significant effects.
Mixed-SEM was fitted using the PiecewiseSEM package [65, 68], available in the R environment [56]. We included random variables in the model to account for the hierarchical structure of our data (i.e., individuals’ intercepts were nested within the sites; see [65]). We also took into consideration the uncertainty associated with the estimated values of movement patterns and home range size. Accordingly, we used an autoregressive error structure to weigh the contribution of the values of the response variables for the inverse of its variance in the corresponding equations [69, 70]. We standardized the estimated coefficients to allow comparison between the different parameters [71]. Finally, we checked the global goodness-of-fit of our mixed-SEM by a Fischer’s C test, which measures the discrepancy between predicted and observed covariance matrices of our causal predictions [72].
Results
General results
We GPS-tracked 19 individuals–six individuals (three males and three females) in Santa Barbara Ecological Station (SP) in 2015, and 13 individuals (eight males and five females) in Baía das Pedras Ranch (MS) between 2013 and 2017. The individuals weighed between 21.6 kg and 38.7 kg (mean = 32.5 kg). Individuals’ body mass was similar for both study areas (t test; t = - 0.69; df = 12.53; p = 0.50). The GPS devices recorded location points at intervals ranging between 20 and 70 minutes. The monitoring time varied between individuals, ranging from 45 to 136 days in Santa Barbara Ecological Station (mean = 90 days) and from 69 to 509 days in Baía das Pedras Ranch (mean = 371.5 days). The total dataset consisted of 213,901 locations. We provided individual information on sex, body mass, sample regime and monitoring time in S1 Table.
The individuals’ empirical variogram showed the plotted semi-variance reaching an asymptote on a timescale that roughly corresponded to the home-range crossing time. Therefore, all the monitored giant anteaters showed constrained space use and were defined as range residents. For all individuals, the highest ranked movement model was the OUF–that takes into account autocorrelation in both location and velocity [26]. The estimates of home range crossing time, daily distance moved, directionality, and home range size varied between individuals (Tables 1 and S2). The mixed-SEM explained a substantial amount of the observed variation in home range crossing time (R² = 0.67), daily distance moved (R² = 0.72), directionality (R² = 0.88), and home range size (R² = 0.81).
Intrinsic effects on movement patterns and home range size
Home range crossing time and daily distance moved had positive influence of sex, with males presenting higher values than females (Fig 2; β1 and β2, respectively; Fig 3A and 3B). The effect of body mass in home range crossing time and daily distance moved depended on the sex, and it was negative for males (Fig 2; β3 and β5, respectively; Fig 3A and 3B) and positive for females (Fig 2; β4 and β6, respectively; Fig 3A and 3B). Directionality was indirectly driven by sex through daily distance moved (Fig 2; β2* β7 = - 1.56; Fig 3C). The effect of body mass on directionality was also given indirectly via daily distance moved (Fig 2; β5* β7 for males, and β6* β7 for females), and it was equal to 1.29 for males and—0.50 for females (Fig 3C). The effect of both sex and body mass on home range size was mediated by home range crossing time, daily distance moved and directionality. The total effect of sex on home range size was given by β1* β8 + β2* β7* β9 = 0.45 (males > females; Figs 2 and 3D). The total effect of body mass on home range size was given by β3* β8 + β5* β7* β9 = 0.41 for males, and β4* β8 + β6* β7* β9 = 1.27 for females (Figs 2 and 3D).
Intrinsic traits are represented by sex and body mass. The environmental trait is represented by the proportion of forest cover within individuals’ home ranges. Movement patterns are described by home range crossing time, daily distance moved and directionality. Estimated coefficients are provided above tendency lines for each relationship.
Environmental effect on movement patterns and home range size
The proportion of forest cover within the individuals’ home range ranged between 0.17 and 0.88 (mean = 0.42, S2 Table). The proportion of forest had no influence on home range crossing time or daily distance moved (Figs 2 and 3E and 3F), however, it negatively influenced directionality (Fig 2; β10 = - 0.60; Fig 3G). It means that individuals whose home range presented a higher proportion of forest cover showed less directionality–i.e., more tortuous movements–than individuals occupying areas with a low proportion of forests. The total effect of the proportion of forest cover on home range size resulted of the sum of its direct effect (Fig 2; β11) with indirect effects, mediated by directionality (Fig 2; β10* β9), totaling an effect of—0.68 (Figs 2 and 3H). In other words, individuals increased home range size with a decreasing proportion of forests inside it.
Discussion
As far as we know, we provided here for the first-time movement model-based estimates of home range crossing time, daily distance moved and directionality for giant anteaters, allowing a better characterization of the species’ movement patterns. Once home range crossing time indicates the timescale of autocorrelation in position, our results show that, on average, a two-days interval between consecutive relocations is necessary to consider them spatially independent [63]. However, previous studies using GPS devices have adopted monitoring regimes much shorter than that and unconsidered the spatial autocorrelation of data [30, 42], probably leading to underestimating results. This is an important factor explaining why our estimates of daily distance moved, and home range size were, in average, bigger than those provided by recent studies [30, 42], once AKDEc incorporates and controls for the autocorrelation in both location and velocity [63]. Other possible explanations to this discrepancy can be related to the intrinsic characteristics of monitored individuals and the environmental characteristics associated with the site and period of monitoring.
The use of mixed-SEM allowed us to disentangle the effects of sex, body mass and proportion of forest on movement patterns and home range size, simultaneously estimating the direction and intensity of direct and indirect effects. The three descriptors of movement patterns (home range crossing time, daily distance moved and directionality) mediated the effects of sex and body mass on home range size, illustrating the importance of integrating these relationships in the same model [65]. We were able to clarify the sexual effects on movement patterns and home range size by considering body mass effects and the interactions between individuals’ sex and body mass. Even though the intraspecific effect of body mass on movement patterns and home range size is generally weak in mammals [12], we were still able to detect it with this integrative approach. To our best knowledge, this is the first record of allometric scaling in the movement patterns and home range size in giant anteaters. On the other hand, environmental traits, such as the proportion of forest cover, are common direct drivers of mammal’s home range size at the individual level [12], and our model provided additional details, showing the direct and indirect paths of this effect. Despite the great explanatory power of our model, we recognize that there must be other intrinsic and environmental factors influencing giant anteaters’ movement that we did not investigate here, and some of them may even seasonally change.
In general, male giant anteaters presented more intensive space use than females with similar body mass, showing longer home range crossing-time, longer daily distance moved and smaller directionality. Besides, males also exhibited larger home ranges than females. We expected males moving longer distances and occupying larger areas than females. However, it was surprising that they also took more time to cross their areas and were less directional at doing it. The higher intensity of space use and larger home range in males than females are, probably, strategies to increase the chances to find receptive females on landscape [44]. This is because the home range of a male giant anteater usually overlaps with the home range of several females [73]. Hence, males could increase their chances of finding receptive females exploiting their home ranges and increasing their home ranges size to include more females inside it. In line with these results, male giant anteaters were recorded presenting longer activity time and larger home range than females [42]. Therefore, we reinforced the idea that giant anteaters present a polygynous social mating system, with a male mating with more than one female [43, 44].
Female giant anteaters behaved as expected, increasing the intensity of space use with increasing body mass. This is probably related to an increase in the search for food resources [26], once larger animals have higher energetic requirements [45]. For species that have their food resources unpredictably spread on the landscape, such as the invertebrate nests that giant anteaters prey upon, the intensity of utilization of food resources depends on the intensity of use of space that provides physical access to those resources [74, 75]. Besides, female giant anteaters also increased home range size with increasing body mass, showing a second strategy to increase the access to food resources: increasing the size of the space used to find those resources. A positive allometric scaling of both intensity of space use and home range size has been found in some mammals [76, 77], including other xenarthrans with myrmecophagous diets such as giant armadillos (Priodontes maximus) [78] and southern three-banded armadillos (Tolypeutes matacus) [79]. This relationship indicates that the search for energetic resources is one of the main factors driving female giant anteaters’ movement across body mass.
On the other hand, males did not display the same pattern. With increasing body mass, males reduced intensity of space use and increased home range size. This reveals a change of males’ movement strategy guided by body mass: while small males used their small areas intensively, large males ranged over large areas with comparatively lower intensity of use. Considering a limited quantity of metabolic energy available for movement [80], animals moving close to their limit capacity should experiment a trade-off between the intensity of use and the area size, and this can be the case of male giant anteaters. Both strategies can increase the access to both food resources spread on the landscape and receptive females [81]. Meanwhile, small males could minimize the chances of agonistic interactions with other males if they use smaller areas than the big ones [82, 83]. Further studies, such as behavioral assessments, will help us to confirm these hypotheses and better understand the species’ reproductive biology.
As we expected, male and female giant anteaters reduced the directionality as the proportion of forest patches inside their home ranges increased, probably due to the physical obstacles that forests impose to displacement [84]. Similarly, small mammals have presented shorter step lengths and higher tortuosity within forest areas [85], and African wild dogs have shown that the movement permeability of the vegetation decreases with its increasing density [86]. However, it is worth noting that, contrary to our expectations, a greater proportion of forest inside the home ranges did not influence the home range crossing time or daily distance moved.
Decreasing the proportion of forest inside giant anteaters’ home range led males and females to an increase in the home range size. This is probably because giant anteaters present a low capacity for physiological thermoregulation [37], and less forest implies less access to thermal shelters and, consequently, a reduced habitat thermal quality [34–36]. As a result, animals would increase the home range size as a strategy to maximize the chances of accessing this thermal resource. Supporting this idea, the increase of home range size with decreasing habitat quality has been widely documented for terrestrial vertebrates [87–89]. Furthermore, the importance of forests as thermal shelters has also been shown for other mammals [90, 91], and it should increase with the predicted increasing frequency of extreme weather events [92]. In the Brazilian territory, where this study was conducted, massive agricultural expansion has caused extensive habitat degradation and dramatically decreased forest patches on savannah areas in number and size [93, 94]. In this current deforestation scenario, our results bring an important implication for giant anteaters’ management: the minimal area needed to preserve a given giant anteaters’ population should increase as the proportion of forests inside it decreases.
In this study, we brought reliable measures of giant anteaters’ movement patterns and home-range size, showing that their movements are influenced by sex, body mass and proportion of forest cover; and revealed two important strategies used by giant anteaters to maximize the access to resources: they modulate movement patters, increasing space use intensity, and/or increasing home range size. This information contributes to the understanding of giant anteaters’ spatial ecology and can help define the spatial scale of effective management efforts for their conservation [95], especially as the anthropogenic impacts on landscapes increase. We highlight the need to consider the sexual differences on movement strategies and the role of forests as an important thermal resource driving giant anteaters’ spatial requirements (also see [96]). In accordance with [36], we strongly suggest that management efforts should focus on maintaining the giant anteaters’ access to forest patches inside their home ranges to provide environmental conditions for behavioral thermoregulation. Both intrinsic and environmental traits driving animal movement should be integrated when establishing conservation strategies for populations and species.
Supporting information
S1 Table. Intrinsic characteristics and monitoring information of tracked giant anteaters.
https://doi.org/10.1371/journal.pone.0253345.s001
(DOC)
S2 Table. Individual estimates of movement patterns, home range size and proportion of forests inside home range of tracked giant anteaters.
https://doi.org/10.1371/journal.pone.0253345.s002
(DOC)
Acknowledgments
We thank Gabriel Massocato, Danilo Kluyber, Alexandre Martins Costa Lopes, Vinicius P. O. Gasparotto, Rullian C. Ribeiro and Clarice S. Cesário who helped to capture and monitor the giant anteaters.
References
- 1. Morales JM, Moorcroft PR, Matthiopoulos J, Frair JL, Kie JG, Powell RA et al. Building the bridge between animal movement and population dynamics. Philos Trans R Soc Lond B Biol Sci. 2010; 365: 2289–2301. pmid:20566505
- 2. Gautestad AO. Memory matters: influence from a cognitive map on animal space use. J Theor Biol. 2011; 287: 26–36. pmid:21810430
- 3. Powell RA, Michael SM. What is a home range?. J Mammal. 2012; 93: 948–958.
- 4. Burt WH. Territoriality and home range concepts as applied to mammals. J Mammal. 1943; 24: 346–352.
- 5. Börger L, Dalziel BD, Fryxell JM. Are there general mechanisms of animal home range behaviour? A review and prospects for future research. Ecol Lett. 2008; 11: 637–650. pmid:18400017
- 6. Allen AM, Singh NJ, 2016. Linking movement ecology with wildlife management and conservation. Front Ecol Evol. 2016; 3: 155.
- 7. Nathan R, Getz WM, Revilla E, Holyoak M, Kadmon R, Saltz D et al. A movement ecology paradigm for unifying organismal movement research. Proc Natl Acad Sci U S A. 2008; 105: 19052–19059. pmid:19060196
- 8. Doherty TS, Hays GC, Driscoll DA. Human disturbance causes widespread disruption of animal movement. Nat Ecol Evol. 2021; 5: 513–519. pmid:33526889
- 9. Shaw AK. Causes and consequences of individual variation in animal movement. Mov Ecol. 2020; 8: 1–12. pmid:31921423
- 10. Vieira MV, Loretto D, Papi B. Scaling of movements with body mass in a small opossum: evidence for an optimal body size in mammals. J Mammal. 2019; 100: 1765–1773.
- 11. Gaulin SJ, FitzGerald RW. Sex differences in spatial ability: an evolutionary hypothesis and test. Am Nat. 1986; 127: 74–88.
- 12. McLoughlin PD, Ferguson SH. A hierarchical pattern of limiting factors helps explain variation in home range size. Ecoscience. 2000; 7: 123–130.
- 13. Wilson RP, Griffiths IW, Mills MG, Carbone C, Wilson JW, Carbone C et al. Mass enhances speed but diminishes turn capacity in terrestrial pursuit predators. Elife. 2015; 4: e06487. pmid:26252515
- 14. Daley MA, Birn-Jeffery A. Scaling of avian bipedal locomotion reveals independent effects of body mass and leg posture on gait. J Exp Biol. 2018; 221: 10. pmid:29789347
- 15. McNab BK. Bioenergetics and the determination of home range size. Am Nat. 1963; 97: 133–140.
- 16. Rosten CM, Gozlan RE, Lucas MC. Allometric scaling of intraspecific space use. Biol Lett. 2016; 12: 20150673. pmid:26979558
- 17.
Viswanathan GM, Da Luz MG, Raposo EP, Stanley HE, 2011. The Physics of Foraging: an Introduction to Random Searches and Biological Encounters. 1st ed. Edinburgh: Cambridge University Press; 2011.
- 18. Papastamatiou YP, Cartamil DP, Lowe CG, Meyer CG, Wetherbee BM, Holland KN. Scales of orientation, directed walks and movement path structure in sharks. J Anim Ecol. 2011; 80: 864–874. pmid:21366562
- 19. Ross S, Munkhtsog B, Harris S. Determinants of mesocarnivore range use: relative effects of prey and habitat properties on Pallas’s cat home-range size. J Mammal. 2012; 93: 1292–1300.
- 20. Mancinelli S, Boitani L, Ciucci P. Determinants of home range size and space use patterns in a protected wolf (Canis lupus) population in the central Apennines, Italy. Can J Zool. 2018; 96: 828–838.
- 21. Mangipanea LS, Belant JL, Hiller TL, Colvin ME, Gustine DD, Mangipane BA et al. Influences of landscape heterogeneity on home-range sizes of brown bears. Mamm Biol. 2018; 88: 1–7.
- 22. Cagnacci F, Boitani L, Powell RA, Boyce MS. Animal ecology meets GPS-based radiotelemetry: a perfect storm of opportunities and challenges. Philos Trans R Soc Lond B Biol Sci. 2010; 365: 2157–2162. pmid:20566493
- 23. Fleming CH, Noonan MJ, Medici EP, Calabrese JM. Overcoming the challenge of small effective sample sizes in home-range estimation. Methods Ecol Evol. 2019; 10: 1679–1689.
- 24.
Fleming CH, Calabrese JM. ctmm: Continuous-Time Movement Modeling. R Package Version 0.3.2. 2016. Available from: http://cran.r-project.org/package=ctmm.
- 25. Katzner TE, Arlettaz R. Evaluating contributions of recent tracking-based animal movement ecology to conservation management. Front Ecol Evol. 2020; 7: 519.
- 26. Fleming CH, Calabrese JM, Mueller T, Olson KA, Leimgruber P, Fagan WF. From fine-scale foraging to home ranges: a semivariance approach to identifying movement modes across spatiotemporal scales. Am Nat. 2014; 183: E154–E167. pmid:24739204
- 27. Morato RG, Stabach JA, Fleming CH, Calabrese JM, De Paula RC, Ferraz KM et al. Space use and movement of a neotropical top predator: the endangered jaguar. PloS One. 2016; 11: e0168176. pmid:28030568
- 28.
De la Torre JA, Rivero M. Insights of the Movements of the Jaguar in the Tropical Forests of Southern Mexico. In: Reyna-Hurtado R, Chapman C, editors. Movement Ecology of Neotropical Forest Mammals. Cham: Springer; 2019. pp. 217–241.
- 29. Miranda F, Bertassoni A, Abba AM. Myrmecophaga tridactyla. The IUCN Red List of Threatened Species 2014. 2014. Avaiable from: https://doi.org/10.2305/IUCN.UK.2014-1.RLTS.T14224A47441961.en
- 30. Bertassoni A, Ribeiro MC. Space use by the giant anteater (Myrmecophaga tridactyla): a review and key directions for future research. Eur J Wildl Res. 2019; 65: 1–11.
- 31.
Rodrigues FHG, Medri IM, Miranda GHB, Camilo-Alves C, Mourão G. Anteater behavior and ecology. In: Vizcaíno SF, Loughry WJ, editors. The Biology of the Xenarthra. Gainesville: University Press of Florida; 2008. pp. 257–268.
- 32.
Miranda FR, Chiarello AG, Röhe F, Braga FG, Mourão GDM, Miranda GD et al. Avaliação do risco de extinção de Myrmecophaga tridactyla Linnaeus 1758 no Brasil. In: ICMBio, editors. Avaliação do Risco de Extinção dos Xenartros Brasileiros. Brasília: ICMBio; 2015. pp. 89–105.
- 33.
Eisenberg JF, Redford KH. The contemporary mammalian fauna. In: Eisenberg JF, Redford KH, editors. Mammals of the Neotropics—The central Neotropics: Ecuador, Peru, Bolivia, Brazil. Chicago: The University of Chicago Press; 1999. pp. 49–522.
- 34. Camilo-Alves CDSEP, Mourão GDM. Responses of a specialized insectivorous mammal (Myrmecophaga tridactyla) to variation in ambient temperature. Biotropica. 2006; 38: 52–56.
- 35. Mourão G, Medri ÍM. Activity of a specialized insectivorous mammal (Myrmecophaga tridactyla) in the Pantanal of Brazil. J Zool. 2007; 271: 187–192.
- 36. Giroux A, Ortega Z, Bertassoni A, Desbiez ALJ, Kluyber D, Massocato GF et al. The role of environmental temperature on movement patterns of giant anteaters. Integr Zool. 2021; 0: 1–12. pmid:33738919
- 37. McNab BK. Physiological convergence amongst ant-eating and termite-eating mammals. J Zool. 1984; 203: 485–510.
- 38. De Frenne P, Zellweger F, Rodríguez-Sánchez F, Scheffers BR, Hylander K, Luoto M et al. Global buffering of temperatures under forest canopies. Nat Ecol Evol. 2019; 3:744–749. pmid:30936433
- 39. Zeppetello LRV, Parsons LA, Spector JT, Naylor RL, Battisti DS, Masuda YJ et al. Large scale tropical deforestation drives extreme warming. Environ Res Lett. 2020; 15: 084012.
- 40. Di Blanco YE, Desbiez AL, Jiménez-Pérez I, Kluyber D, Massocato GF, Di Bitetti MS. Habitat selection and home-range use by resident and reintroduced giant anteaters in 2 South American wetlands. J Mammal. 2017; 98: 1118–1128.
- 41. Di Blanco YE, Spørring KL, Di Bitetti MS. Daily activity pattern of reintroduced giant anteaters (Myrmecophaga tridactyla): effects of seasonality and experience. Mammalia. 2017; 81: 11–21.
- 42. Bertassoni A, Mourão G, Bianchi RDC. Space use by giant anteaters (Myrmecophaga tridactyla) in a protected area within human-modified landscape. Ecol Evol. 2020; 10: 7981–7994. pmid:32788955
- 43. Desbiez ALJ, Bertassoni A, Traylor-Holzer K. Population viability analysis as a tool for giant anteater conservation. Perspect Ecol Conserv. 2020; 18: 124–131.
- 44. Clutton-Brock TH. Mammalian mating systems. Proc R Soc Lond B Biol Sci. 1989; 236: 339–372. pmid:2567517
- 45.
Isaac NJB, Carbone C, McGill B. Population and community ecology. In: Sibly RM, Browh JH, KodricBrown A, editors. Metabolic Ecology: a Scaling Approach. Chichester: John Wiley & Sons; 2012. pp. 77–85.
- 46. Redford KH. Feeding and food preference in captive and wild giant anteaters (Myrmecophaga tridactyla). J Zool. 1985; 205: 559–572.
- 47. Ferreras P. Landscape structure and asymmetrical inter-patch connectivity in a metapopulation of the endangered Iberian lynx. Biol Conserv. 2001; 100: 125–136.
- 48. Saïd S, Gaillard JM, Widmer O, Débias F, Bourgoin G. What shapes intra-specific variation in home range size? A case study of female roe deer. Oikos. 2009; 118: 1299–1306.
- 49. Alvares CA, Stape JL, Sentelhas PC, Gonçalves JDM, Sparovek G. Köppen’s climate classification map for Brazil. Meteorol Z. 2013; 22: 711–728.
- 50. Abdon MM, da Silva JDSV, Pott VJ, Pott A, da Silva MP. Utilização de dados analógicos do Landsat-TM na discriminação da vegetação de parte da sub-região da Nhecolândia no Pantanal. Pesqui Agropecu Bras. 1995; 33: 1799–1813.
- 51. Durigan G, Siqueira MFD, Franco GADC. Threats to the Cerrado remnants of the state of São Paulo, Brazil. Sci Agric. 2007; 64: 355–363.
- 52. Souza CM, Shimbo J, Rosa MR, Parente LL, Alencar AA, Rudorff BF et al. Reconstructing three decades of land use and land cover changes in Brazilian biomes with Landsat Archive and Earth Engine. Remote Sens. 2020; 12: 2735.
- 53.
Hijmans RJ, Van Etten J, Sumner M, Cheng J, Baston D, Bevan A et al. ‘Package ‘raster’. R Package Version 3.4–5, 734. 2014. Available from: https://CRAN.R-project.org/package=raster
- 54.
Bivand R, Lewin-Koh N. maptools: tools for reading and handling spatial objects. R package version 0.8–39. 2016. Available from: https://CRAN.R-project.org/package=maptools.
- 55.
Bivand R, Keitt T, Rowlingson B. rgdal: bindings for the geospatial data abstraction library. R package version 1.1–3. 2015. Available from: https://CRAN.R-project.org/package=rgdal.
- 56.
R Core Team. R: A language and environment for statistical computing. Vienna: R Foundation for Statistical Computing. 2019.
- 57. Kluyber D, Lopez RPG, Massocato G, Attias N, Desbiez ALJ. Anesthesia and surgery protocols for intra-abdominal transmitter placement in four species of wild armadillo. J Zoo Wildl Med. 2020; 51: 514–526. pmid:33480527
- 58.
Young LJ, Young J. Statistical ecology: a population perspective. 1st ed. Boston: Kluwer Academic Publishers; 1998.
- 59. Sikes RS. Guidelines of the American Society of Mammalogists for the use of wild mammals in research and education. J Mammal. 2016; 97: 663–688. pmid:29692469
- 60. Fleming CH, Fagan WF, Mueller T, Olson KA, Leimgruber P, Calabrese JM. Rigorous home range estimation with movement data: a new autocorrelated kernel density estimator. Ecology. 2015; 96: 1182–1188. pmid:26236833
- 61. Calabrese JM, Fleming CH, Gurarie E. ctmm: an R package for analyzing animal relocation data as a continuous-time stochastic process. Methods Ecol Evol. 2016; 7: 1124–1132.
- 62.
Burnham KP, Anderson DR. Model Selection and Multimodel Inference: a Practical Information-theoretic Approach. 2nd ed. New York: Springer; 2002.
- 63. Fleming CH, Calabrese JM, 2017. A new kernel-density estimator for accurate homerange and species-range estimation. Methods Ecol Evol. 2017; 8: 571–579.
- 64.
Kline RB. Principles and Practice of Structural Equation Modeling. 4th ed. London: The Guilford Press; 2015.
- 65. Lefcheck JS. piecewiseSEM: Piecewise structural equation modelling in R for ecology, evolution, and systematics. Methods Ecol Evol. 2016; 7: 573–579.
- 66. Stolzenberg RM. The measurement and decomposition of causal effects in nonlinear and nonadditive models. Sociol Methodol. 1980; 11: 459–488.
- 67. Hayes AF, Preacher KJ. Quantifying and testing indirect effects in simple mediation models when the constituent paths are nonlinear. Multivariate Behav Res. 2010; 45: 627–660. pmid:26735713
- 68.
Heck RH. Multilevel modeling with SEM. In: Marcoulides GA, Schumacker RE, editors. New developments and techniques in structural equation modeling. Mahwah: Erlbaum, 2001. pp. 89–128.
- 69. Sterne JA, Egger M. Funnel plots for detecting bias in meta-analysis. J Clin Epidemiol. 2001; 54: 1046–1055. pmid:11576817
- 70. Lin L, Chu H. Meta-analysis of proportions using Generalized Linear Mixed Models. Epidemiology. 2020; 31: 713–717. pmid:32657954
- 71. Grace JB, Johnson DJ, Lefcheck JS, Byrnes JE. Quantifying relative importance: computing standardized effects in models with binary outcomes. Ecosphere. 2018; 9: e02283.
- 72. Holst KK, Budtz-Jørgensen E. Linear latent variable models: the lava-package. Comput Stat. 2013; 28: 1385–1452.
- 73. Medri ÍM, Mourão G. Home range of giant anteaters (Myrmecophaga tridactyla) in the Pantanal wetland, Brazil. J Zool. 2006; 266: 365–375.
- 74. De Knegt HJ, Hengeveld GM, Van Langevelde F, De Boer WF, Kirkman KP. Patch density determines movement patterns and foraging efficiency of large herbivores. Behav Ecol 2007; 18: 1065–1072.
- 75. Augustine DJ, Derner JD. Assessing herbivore foraging behavior with GPS collars in a semiarid grassland. Sensors. 2013; 13: 3711–3723. pmid:23503296
- 76. Cameron GN, Spencer SR. Assessment of space-use patterns in the hispid cotton rat (Sigmodon hispidus). Oecologia. 1985; 68: 133–139. pmid:28310922
- 77. Duncan C, Nilsen EB, Linnell JD, Pettorelli N. Life-history attributes and resource dynamics determine intraspecific home-range sizes in Carnivora. Remote Sens Ecol Conserv. 2015; 1: 39–50.
- 78. Desbiez ALJ, Kluyber D, Massocato GF, Oliveira-Santos LGR, Attias N. Spatial ecology of the giant armadillo Priodontes maximus in Midwestern Brazil. J Mammal. 2019; 101: 151–163.
- 79. Attias N, Gurarie E, Fagan WF, Mourão G. Ecology and social biology of the southern three-banded armadillo (Tolypeutes matacus; Cingulata: Chlamyphoridae). J Mammal. 2020; 101: 1692–1705.
- 80. Sparrow WA, Newell KM. Metabolic energy expenditure and the regulation of movement economy. Psychon Bull Rev. 1998; 5: 173–196.
- 81. Shepard EL, Wilson RP, Rees WG, Grundy E, Lambertucci SA, Vosper SB. Energy landscapes shape animal movement ecology. Am Nat. 2013; 182: 298–312. pmid:23933722
- 82. Rocha FL, Mourão G. An agonistic encounter between two giant anteaters (Myrmecophaga tridactyla). Edentata 2006; 2006: 50–51.
- 83. Kreutz K, Fischer F, Linsenmair KE. Observations of intraspecific aggression in giant anteaters (Myrmecophaga tridactyla). Edentata. 2009; 8: 6–7.
- 84.
Ims RA. Movement patterns related to spatial structures. In: Hansson L, Fahrig L, Merriam G, editors. Mosaic Landscapes and Ecological Processes. Dordrecht: Springer; 1995. pp. 85–109.
- 85. Wells K, Pfeiffer M, Lakim MB, Kalko EK. Movement trajectories and habitat partitioning of small mammals in logged and unlogged rain forests on Borneo. J Anim Ecol. 2006; 75: 1212–1223. pmid:16922857
- 86. Abrahms B, Jordan NR, Golabek KA, McNutt JW, Wilson AM, Brashares JS. Lessons from integrating behaviour and resource selection: activity-specific responses of A frican wild dogs to roads. Anim Conserv. 2015; 19: 247–255.
- 87. Ofstad EG, Herfindal I, Solberg EJ, Sæther BE. Home ranges, habitat and body mass: simple correlates of home range size in ungulates. Proc R Soc Lond B Biol Sci. 2016; 283: 20161234. pmid:28003441
- 88. Gardiner R, Proft K, Comte S, Jones M, Johnson CN. Home range size scales to habitat amount and increasing fragmentation in a mobile woodland specialist. Ecol Evol. 2019; 9: 14005–14014. pmid:31938498
- 89. Mayer M, Ullmann W, Heinrich R, Fischer C, Blaum N, Sunde P. Seasonal effects of habitat structure and weather on the habitat selection and home range size of a mammal in agricultural landscapes. Landsc Ecol. 2019; 34: 2279–2294.
- 90. Melin M, Matala J, Mehtätalo L, Tiilikainen R, Tikkanen OP, Maltamo M et al. Moose (Alces alces) reacts to high summer temperatures by utilizing thermal shelters in boreal forests–an analysis based on airborne laser scanning of the canopy structure at moose locations. Glob Chang Biol. 2014; 20: 1115–1125. pmid:24115403
- 91. Attias N, Oliveira-Santos LGR, Fagan WF, Mourão G. Effects of air temperature on habitat selection and activity patterns of two tropical imperfect homeotherms. Anim Behav. 2018; 140:129–140.
- 92. Meehl GA, Zwiers F, Evans J, Knutson T, Mearns L, Whetton P. Trends in extreme weather and climate events: issues related to modeling extremes in projections of future climate change. Bull Am Meteorol Soc. 2000; 81: 427–436.
- 93. Tollefson J. Deforestation ticks up in Brazil’s savannah. Nature. 2018; 12.
- 94. Reynolds J, Wesson K, Desbiez ALJ, Ochoa-Quintero JM, Leimgruber P. Using remote sensing and random forest to assess the conservation status of critical cerrado habitats in Mato Grosso Do Sul, Brazil. Land. 2016; 5: 12.
- 95. Stewart JD, Beale CS, Fernando D, Sianipar AB, Burton RS, Semmens BX et al. Spatial ecology and conservation of Manta birostris in the Indo-Pacific. Biol Conserv. 2016; 200: 178–183.
- 96. Desbiez ALJ, Kluyber D, Massocato GF, Oliveira-Santos LGR, Attias N. Life stage, sex, and behavior shape habitat selection and influence conservation strategies for a threatened fossorial mammal. Hystrix. 2020; 31: 1–7.