Life on the edge: behavioural and physiological responses of Verreaux's sifakas (Propithecus verreauxi) to forest edges
←
→
Page content transcription
If your browser does not render page correctly, please read the page content below
Research article Primate Biol., 8, 1–13, 2021 https://doi.org/10.5194/pb-8-1-2021 © Author(s) 2021. This work is distributed under the Creative Commons Attribution 4.0 License. Life on the edge: behavioural and physiological responses of Verreaux’s sifakas (Propithecus verreauxi) to forest edges Klara Dinter1 , Michael Heistermann2 , Peter M. Kappeler1,3 , and Claudia Fichtel1 1 Behavioral Ecology and Sociobiology Unit, German Primate Center, Kellnerweg 4, 37077 Göttingen, Germany 2 Endocrinology Laboratory, German Primate Center, Kellnerweg 4, 37077 Göttingen, Germany 3 Department of Sociobiology and Anthropology, Johann-Friedrich-Blumenbach Institute for Zoology and Anthropology, Georg-August University, Kellnerweg 6, 37077 Göttingen, Germany Correspondence: Claudia Fichtel (claudia.fichtel@gwdg.de) Received: 19 August 2020 – Revised: 13 December 2020 – Accepted: 16 December 2020 – Published: 9 February 2021 Abstract. Forest edges change micro-environmental conditions, thereby affecting the ecology of many forest- dwelling species. Understanding such edge effects is particularly important for Malagasy primates because many of them live in highly fragmented forests today. The aim of our study was to assess the influence of forest edge effects on activity budgets, feeding ecology, and stress hormone output (measured as faecal glucocorticoid metabolite – fGCM – levels) in wild Verreaux’s sifakas (Propithecus verreauxi), a group living, arboreal lemur. We observed five habituated groups: three living in the forest interior and two at an established forest edge. There was no difference in average daily temperatures between edge and interior habitats; however, within the edge site, the average daily temperature incrementally increased over 450 m from the forest edge towards the interior forest of the edge habitat, and the population density was lower at the edge site. Activity budgets differed between groups living in the two microhabitats, with individuals living near the edge spending more time travelling and less time feeding. Groups living near the edge also tended to have smaller home ranges and core areas than groups in the forest interior. In addition, edge groups had elevated average fGCM concentrations, and birth rates were lower for females living in the edge habitat. Combined with lower levels of fruit consumption at the edge, these results suggest that nutritional stress might be a limiting factor for Verreaux’s sifakas when living near a forest edge. Hence, Verreaux’s sifakas appear to be sensitive to microhabitat characteristics linked to forest edges; a result with implications for the conservation of this critically endangered lemurid species. 1 Introduction bance, change around edges, they can be considered as dis- tinct ecosystems in comparison to undisturbed forest interi- Tropical forest habitats continue to become increasingly ors (Laurance et al., 2000, 2011). Species responses towards fragmented, thereby disintegrating populations and reducing edges are diverse (Lidicker, 1999): species that exhibit their suitable habitats for many species (Laurance, 1999; Watson highest densities at the edge are considered to show an edge- et al., 2004). Fragmentation of habitats results in additional positive response, whereas species that avoid forest edges are abrupt boundaries, creating so-called edge effects. Edge ef- considered to exhibit an edge-negative response, and species fects refer to how abiotic and biotic conditions interact be- that show little or no response towards the edge are consid- tween two adjacent habitats separated by an abrupt boundary, ered to be edge-neutral. Because species respond differently altering the distribution, abundance, and behaviour of organ- to living in edge habitats (Lehtinen et al., 2003; Lehman et isms (Murcia, 1995). As environmental conditions, such as al., 2006a, b; Delattre et al., 2009), studying edge effects is of temperature, light intensity, plant growth, and level of distur- key importance for species-specific conservation assessment. Published by Copernicus Publications on behalf of the Deutsches Primatenzentrum GmbH (DPZ).
2 K. Dinter et al.: Behavioural and physiological responses of Verreaux’s sifakas The effects of anthropogenic habitat changes are many- Among sifakas (genus Propithecus: Indridae), which are fold. First, numerous studies have demonstrated that dis- folivorous and, hence, potentially less resilient to habitat turbed environments are associated with elevated stress hor- fragmentation (Eppley et al., 2020), responses to edges differ mone, i.e. glucocorticoid (GC), levels in many vertebrate between species. Coquerel’s sifakas (Propithecus coquereli) taxa, including primates (Dantzer et al., 2014). Specifically, seem to avoid edge areas and do not range right on the for- in Yucatan spider monkeys (Ateles geoffroyi yucatanensis; est edge but within 400 m of it, resulting in larger home Rangel-Negrín et al., 2009), brown spider monkeys (Ate- ranges in comparison with groups occurring in the forest in- les hybridus; Rimbach et al., 2013a), black howler mon- terior, although food quality does not differ between habi- keys (Alouatta pigra; Martínez-Mota et al., 2007; Rangel- tats (McGoogan et al., 2009). Diademed sifakas (Propithe- Negrín et al., 2014), and grey-cheeked mangabeys (Lopho- cus diadema) that live in forest fragments have smaller home cebus albigena; Jaimez et al., 2012) individuals living in dis- ranges and appear to avoid feeding at fragment edges in con- turbed and fragmented forests exhibited elevated GC lev- trast to groups living in the forest interior (Irwin, 2008). Sub- els compared with those inhabiting continuous and undis- sequent blood analyses revealed a negative edge response turbed forest patches. Collared brown lemur (Eulemur col- on the physiological health of these animals as reflected by laris) also exhibited elevated GC concentrations in degraded lower values of white blood cell counts, total protein, or forests, where reduced fruit availability exacerbated the eco- iron-binding capacity (Irwin et al., 2010). Moreover, mor- logical challenges faced by these lemurs (Balestri et al., phological traits of Diademed sifakas were affected by habi- 2014). Hence, physiological stress responses to disturbed tat quality, with individuals living in smaller fragments with habitats are a major focus in conservation physiology be- low habitat quality exhibiting “wasting effects” in adults and cause they provide a quantitative measure of how environ- “stunting effects” in juveniles (Irwin et al., 2019). mental changes impact individuals (Tarlow and Blumstein, Given the interspecific variation in terms of behavioural 2007; Dantzer et al., 2014). Second, habitat fragmentation and physiological effects to living in an edge habitat, the had a generally negative effect on the distribution of 32 objective of this study was to investigate how forest edge lemur species (Eppley et al., 2020). Folivorous–frugivorous habitats influence behavioural patterns, feeding ecology, and species, however, showed greater flexibility and variability physiological stress hormone output in Verreaux’s sifakas in terms of habitat area and landscape complexity compared (Propithecus verreauxi), which were recently classified as with nearly exclusively folivorous or frugivorous species, “Critically Endangered” (IUCN, 2020). To this end, we stud- which appear to exhibit some resilience to forest fragmen- ied three groups of sifakas inhabiting the interior of an intact tation (Eppley et al., 2020). Third, population density and forest area and two groups living at the forest edge along a juvenile recruitment can depend on fragment size and the savannah (Rakotoniaina et al., 2016). To determine whether availability of key feeding trees, as shown for ring-tailed Verreaux’s sifakas respond to the associated differences in lemurs (Lemur catta; Gould and Cowen, 2020). Finally, habi- habitat structure, we collected data on ambient temperature, tat fragmentation may eventually also reduce genetic diver- home range size, activity budgets, diet, faecal glucocorticoid sity, as shown in species such as the black-and-white ruffed metabolite output, and reproductive rates, and we compared lemurs (Varecia variegata: Baden et al., 2019) and golden- these variables between the two populations. brown mouse lemurs (Microcebus ravelobensis; Radespiel et al., 2019), leading to increased vulnerability due to stochastic demographic changes and potential bottlenecks. 2 Methods As a biodiversity hotspot, Madagascar suffers from an alarming rate of habitat conversion and fragmentation 2.1 Study site and subjects (Harper et al., 2007; Schwitzer et al., 2014; Vieilledent et al., 2018). Madagascar has lost 44 % of its natural forest cover This study was conducted at the German Primate Center over the 1953–2014 period (including 37 % over the 1973– (DPZ) field station in Kirindy Forest (44◦ 390 E, 20◦ 040 S) 2014 period; Vieilledent et al., 2018). Anthropogenic influ- in central western Madagascar. Data were collected from ences and habitat disturbance also have vital consequences April to July 2012, which corresponds to the transition from for the endemic primates of Madagascar (Lemuriformes), the local wet (November–April) to dry (May–October) sea- which are considered to be the most endangered group of son. Most trees shed their leaves over the dry season, reduc- mammals today (Schwitzer et al., 2014). Some lemurs ex- ing food availability of the predominantly folivorous sifakas hibit a negative edge response, such as Cheirogaleus major (Koch et al., 2017). The local study areas were situated about (Lehman et al., 2006a–c) and M. ravelobensis (Andriatsito- 2 km from each other: one in the forest interior and one bor- haina et al., 2019), but others have been found to exhibit a dering a treeless savannah that has been present for at least neutral (Eulemur rubriventer, Hapalemur griseus, and Avahi several decades (Rakotoniaina et al., 2016). We studied five laniger) or even positive edge response (Microcebus rufus groups with a total number of 20 individuals: three groups in- and Propithecus edwardsi), stressing the importance of stud- habited the forest interior and two groups lived at the forest ies of individual species or populations. edge close to the savannah (Table 1). Both study areas were Primate Biol., 8, 1–13, 2021 https://doi.org/10.5194/pb-8-1-2021
K. Dinter et al.: Behavioural and physiological responses of Verreaux’s sifakas 3 fitted with a grid system of small foot trails (Kappeler and by jumping or climbing (Trillmich et al., 2004). Spending at Fichtel, 2012). least 30 s without engaging in any of the above activities (e.g. Verreaux’s sifakas are diurnal, arboreal, and live in groups sitting, watching, autogrooming, or sleeping) was recorded of 2–12 individuals (mean ± SD: 6±2 individuals; Kappeler as resting. To determine the home range sizes, the location and Fichtel, 2012), which occupy home ranges that can vary of each group was additionally recorded by taking the GPS in size due to seasonality, habitat quality, or group size (Car- position (Garmin CS76) of the focal animal every 15 min. rai and Lunardini, 1996; Benadi et al., 2008; Rudolph et al., 2020). In the interior of Kirindy Forest, home range sizes 2.4 Faecal sample collection and faecal glucocorticoid ranged between 5 and 10 ha (Benadi et al., 2008; Koch et al., metabolite (fGCM) analysis 2016), but corresponding data were not previously available for the forest edge. All study subjects were well habituated To measure stress hormone output via the determination to human observers and marked individually by unique ny- of faecal glucocorticoid metabolite (fGCM) concentrations, lon collars and one radio transmitter per group (Kappeler and three faecal samples per focal individual were collected be- Fichtel, 2012). tween 8:00 and 10:00 UTC+3 (N = 43 except seven sam- ples that were collect between 13:00 and 14:00 UTC+3) in 2.2 Temperature measurement three different sampling sessions every 2 weeks (mean ± SD: 13 ± 2 d) from mid-May to the end of June. Each sample Ambient temperature was measured once per hour continu- consisted of exactly five faecal droplets (in total 0.5–2.2 g ously throughout the study using iButtons (1-Wire/iButton fresh weight) that were collected immediately after defae- DS1921G-F5; accuracy ±1 ◦ C). Both study sites are cation and directly transferred into labelled polypropylene equipped with a grid system with 50 m × 50 m trails in the tubes containing 8 mL of 80 % watery ethanol. Faeces where edge habitat and 25 m × 25 m trails in the forest interior. At urine contamination was suspected were not sampled. All each site, 10 iButtons were placed at breast height on tree samples were extracted on-site on the day of sample collec- trunks, with a distance of 50 m between them, on a perpen- tion, using a validated method for extracting faecal hormone dicular line of the grid systems, and the forest edge (Fig. 2). metabolites under field conditions (Shutt et al., 2012). In the savanna, the first iButton was placed at the border of Briefly, each sample was homogenized in its original the edge and the last iButton was placed 450 m westwards. ethanolic solvent using a metal rod, and the resulting fae- In the forest interior, the first iButton was placed in the cen- cal suspensions were subsequently vortexed for 2 min (Lab tre of the study area and the last iButton was placed 450 m Dancer S40). Samples were finally centrifuged for 2 min us- westwards. ing a manually operated centrifuge. About 1 mL of the su- pernatant was decanted into a labelled 2 mL polypropylene 2.3 Behavioural data collection tube, sealed with PARAFILM, and stored in the dark at am- bient temperature until transport to the endocrinology labo- Behaviour of sifakas was recorded using continuous focal ratory of the German Primate Center within 3 weeks of the sampling (Altmann, 1974) of 1 h duration of all adult indi- end of sample collection. Thus, sample extracts were stored viduals of a group, resulting in 252 h of focal animal data. in the field for a maximum of 4–5 months, a period during Individuals were chosen as focal animals in a randomized which fGCM concentrations have been shown to remain sta- but counterbalanced order throughout the day, and we ob- ble (Shutt et al., 2012; Rimbach et al., 2013b; Kalbitzer and served all individuals equally often. Observed behavioural Heistermann, 2013). states were assigned to one of four categories: feeding (F), After arrival at the laboratory, samples were placed in social behaviour (S), travelling (T), and resting (R). The feed- a −20 ◦ C freezer and analysed for glucocorticoid metabo- ing category included searching, manipulating, and ingest- lites within 2 weeks. Faecal extracts were analysed using ing food items within a food patch (i.e. tree). To character- an enzyme immunoassay for 11β-hydroxyetiocholanolone, a ize dietary composition, we recorded the plant parts of every group-specific assay for the measurement of cortisol metabo- species that a focal individual fed on within a feeding bout lites with a 3α,11β-dihydroxy structure. The assay has pre- for longer than 5 min (e.g. leaves, fruits, buds, or flowers; Be- viously been validated for monitoring glucocorticoid output nadi et al., 2008; Koch et al., 2017). In addition, the GPS po- in numerous species of primates (e.g. Heistermann et al., sition of the plant was noted, the species was identified, and 2006), including Verreaux’s sifakas (Fichtel et al., 2007). feeding tree characteristics such as diameter at breast height The assay was carried out as described in detail by Heis- (dbh), height, and crown diameter were recorded. termann et al. (2004). Prior to assaying, all samples were The social behaviour category included all interactions be- diluted 1 : 400 in assay buffer (0.04 M PBS, pH 7.2), and tween individuals, including affiliation (duration of groom- duplicate 50 µL aliquots were taken to assay. Sensitivity of ing and huddling), agonism (frequency of bites and threats), the assay at 90 % binding was 0.6 pg. Intra- and inter-assay and affinity (approaches within 1 m). The travelling category coefficients of variation of high- and low-value quality con- included all movements or changes in location of at least 5 m trols were 6.5 % (high, n = 16) and 8.4 % (low, n = 16) and https://doi.org/10.5194/pb-8-1-2021 Primate Biol., 8, 1–13, 2021
4 K. Dinter et al.: Behavioural and physiological responses of Verreaux’s sifakas Table 1. Composition of the five study groups at the two study sites (“f” denotes females, and “m” denotes males). Site Group Size Adult Adult Subadults Juveniles females males Edge 2 2 1 1 0 0 3 6 1 3 1 (f) 1 (m) Interior C 4 1 2 1 (f) 0 G 4 1 2 1 (f) 0 L 3 (4) 1 1 1 (m) 1 (f) † † represents that animals died or disappeared during the study. 7.1 % (high, n = 6) and 12.5 % (low, n = 6) respectively. All rates, we assessed whether each female (edge: Nfemales = 4; hormone concentrations are expressed as mass hormone per forest interior: Nfemales = 12) of each group gave birth (yes faecal mass (i.e. wet weight). or no) between 2009 and 2013. We fitted a generalized lin- ear mixed model (GLMM; Zuur et al., 2009) using the lme4 package (Bates et al., 2015) by including a binomial response 2.5 Data analyses term indicating whether or not a female gave birth in a given Differences in temperature between sites were estimated by year and site as fixed factors as well as mother ID and year fitting three linear mixed models (LMMs; Baayen, 2008) as random factors. with either minimum, mean, or maximum temperature as re- In general, after fitting each model, we checked the as- sponse, site as fixed factor, and iButton ID as random factor. sumptions regarding the normality of distributions and ho- The influence of the distance to the edge border on temper- moscedasticity and compared the full model to a null model, ature was estimated with three LMMs with either minimum, including the control factor, using maximum likelihood ratio mean, or maximum temperature as response, distance of the tests. iButton location to the edge border as fixed factor, and iBut- ton ID as random factor. The models were fitted in R (R Core Team, 2015) using the lm function of the “lme4” R package 2.6 Ethical statement (Bates et al., 2015). This study was approved and authorized by the CAFF/CORE Differences in time (rates per hour) spent travelling, feed- of the Direction des Eaux et Forêts Madagascar and the CN- ing, resting, and interacting socially as well as time spent FEREF Morondava, Madagascar. Verreaux’s sifakas (Prop- feeding on leaves, fruits, flowers, buds, and seeds between ithecus verreauxi) are currently classified as critically endan- individuals occurring at the edge or interior forest were com- gered (IUCN, 2020). pared with a Mann–Whitney U test using Past3, as these data were not normally distributed. ArcView (version 3.3) was used to calculate the home 3 Results ranges of the study groups based on minimum convex poly- gons (MCPs). Adaptive kernel analyses were also used to cal- culate home ranges at 95 % of raster resolution and at 50 % 3.1 Sifaka densities for the core areas, using an automatically selected LCSV The study area at the edge comprised 31 ha and was in- smoothing parameter (Burgman and Fox, 2003). habited by two groups consisting of six Verreaux’s sifakas, To test whether fGCM levels were influenced by study site whereas the study area in the forest interior comprised 70 ha (i.e. forest interior vs. edge habitat), we used a linear mixed and housed nine groups consisting of 44 individuals, indicat- model (Baayen, 2008) into which we included median fGCM ing a lower population density in the edge habitat. levels as response, site as a fixed effect, and sex as control factor because some females were gestating. The model was fitted in R (R Core Team, 2015) using the lm function of the 3.2 Ambient temperatures lme4 R package (version 3.3.2; Bates et al., 2015). Groups in both habitats (edge: Ngroups = 3; forest interior: Daily minimum, mean, or maximum temperatures did not Ngroups = 9) are part of a long-term project in which group differ between study sites (Fig. 2a–c; Table 2a–c). Daily compositions including births have been assessed during reg- mean temperatures decreased with incremental distance from ular census visits (Kappeler and Fichtel, 2012). Originally the border of the edge towards the interior of the forest there were three groups of sifakas at the forest edge, but (Fig. 3, Table 2e). However, daily minimum and maximum one group dissolved in 2011. To estimate differences in birth temperatures did not change over distance (Table 2d, f). Primate Biol., 8, 1–13, 2021 https://doi.org/10.5194/pb-8-1-2021
K. Dinter et al.: Behavioural and physiological responses of Verreaux’s sifakas 5 Table 2. Parameter estimates from the linear mixed model (LMM) fitting differences in minimum, mean, and maximum temperature between study sites (a–c) and as a function of the distance to the edge border (d–f) with corresponding full–null model comparisons. Response variable Fixed factor Estimate SE p value (a) Minimum temperature Intercept 14.26 0.16
6 K. Dinter et al.: Behavioural and physiological responses of Verreaux’s sifakas Figure 1. Satellite image of the study area depicting the home ranges of the study groups in the interior (groups C, G and L) and at the edge (groups 2 and 3) with MCPs (minimum convex polygons), the grid system of small foot trails (grey lines), and iButton locations (red dots). Imagery © 2020 Maxar Technologies, CNES/Airbus, Imagery © 2020/CNES/Airbus, Landsat/Copernicus, Maxar Technologies, Map data © 2020. Figure 2. Differences in daily (a) minimum, (b) mean, and (c) maximum temperatures between sites. 3.6 Variation in fGCM concentrations 3.7 Reproductive rates Sifakas at the edge had higher average fGCM concentrations In both habitats, each group included only one adult female. than those inhabiting the forest interior (Fig. 6; LMM, full– Of the two females in the edge population, only one gave null model comparison: χ 2 = 5.59, df = 1, p = 0.038; site: birth, whereas in the three groups living in the forest interior, estimate ± SE: −212.43±89.83, p = 0.038; sex: estimate ± all females gave birth a month after the completion of this SE: −97.58 ± 92.78, p = 0.316). study. To compare reproductive rates, we analysed additional Primate Biol., 8, 1–13, 2021 https://doi.org/10.5194/pb-8-1-2021
K. Dinter et al.: Behavioural and physiological responses of Verreaux’s sifakas 7 Table 3. Results of the home range estimates of each group based on different methods. Site Group Group GPS MCPs K 95 % Core area, Intergroup size points (ha) (ha) K 50 % (ha) encounters Edge 2 2 144 2.8 2.9 0.3 8 3 6 287 4.7 3.1 0.1 1 Interior C 4 216 9.4 9.0 2.0 4 G 4 216 8.2 9.8 1.6 0 L 3 (4) 152 4.7 3.6 0.4 1 MCPs denotes minimum convex polygons, and K denotes kernel analysis. ing in the edge habitat appeared to have lower population densities and smaller home ranges, but they spent more time travelling. In addition, they spent less time feeding but did not differ from the interior groups regarding time spent so- cializing and resting. Both populations had broadly similar diets, but only sifakas in the edge habitats consumed seeds and fed less often on fruits but more often on flowers. These behavioural differences appear to be biologically meaning- ful because sifakas living near the forest edge had elevated fGCM concentrations and lower birth rates. Altogether, these preliminary results based on a small sample size suggest that Verreaux’s sifakas in Kirindy Forest may exhibit a negative edge response. The microclimate in the edge habitat in Kirindy Forest did not differ with respect to daily temperatures from the for- est interior, but it exhibited higher daily mean temperatures directly at the border of the habitat. The edge habitat was, similarly to other edge habitats in Madagascar, characterized Figure 3. Decrease in daily mean temperatures as a function of the by a lower tree density (Lehtinen et al., 2003; Lehman et distance of iButton locations from the border of the forest edge. al., 2006a–c; Rakotoniaina et al., 2016). Despite the lower tree density, both habitats exhibit broadly similar additional characteristics, including traversal by a river and a dirt road. demographic data from 4 females from three different groups However, sifakas at the edge inhabited substantially smaller at the edge and 12 females from nine different groups in the home ranges than groups in the interior; in fact, they repre- interior. Sifakas near the edge had significantly lower birth sent the smallest known home ranges of Verreaux’s sifakas rates compared with those living in the forest interior. Based (Carrai and Lunardini, 1996; Richard, 1978; Benadi et al., on 45 infants born at both sites between 2009 and 2013, the 2008; Koch et al., 2016). Moreover, home ranges of groups average birth rate per female was 0.31 infants per year at the in the interior forest overlapped among groups, whereas in edge and 0.78 infants per year in the interior. The probability the edge habitat, the two groups were separated by an area that a female gave birth was influenced by site but not by year that was not occupied by sifakas. Whereas two groups of (Fig. 7; GLMM, full–null model comparison: χ 2 = 11.77, sifakas in the forest interior regularly crossed the dirt road, df = 1, p
8 K. Dinter et al.: Behavioural and physiological responses of Verreaux’s sifakas Figure 4. Differences in average duration of (a) locomotion, (b) foraging, (c) social interactions, and (d) resting between groups at the two study sites. Depicted are medians (black horizontal lines), interquartile ranges (boxes), and outliers (circles). Figure 5. Differences in average time spent feeding on (a) leaves, (b) buds, (c) seeds, (d) fruits, and (e) flowers. Depicted are medians (black horizontal lines), interquartile ranges (boxes), and outliers (circles). Primate Biol., 8, 1–13, 2021 https://doi.org/10.5194/pb-8-1-2021
K. Dinter et al.: Behavioural and physiological responses of Verreaux’s sifakas 9 ous forests did not avoid the edge (Irwin, 2008). In our study, where the edge presents a naturally derived edge, one group of Verreaux’s sifakas inhabited a home range that encom- passed the border of the edge and, hence, did not appear to avoid the edge. The five most commonly consumed feeding trees of sifakas did not differ in size between sites, and both popula- tions had broadly similar diets. Groups in the forest interior, however, spent more time feeding on fruits, whereas edge groups spent more time feeding on seeds and flowers, which are more quickly ingested than fruits and may require higher rates of switches between feeding sites (Norscia et al., 2006; Irwin, 2008; Gabriel, 2013). This fact may explain why edge- living sifakas spent more time travelling than individuals in the forest interior. As we only have phenology data for the interior forest, we cannot evaluate whether there were gener- ally fewer fruiting trees present in the edge habitat. The pres- Figure 6. Median faecal glucocorticoid metabolite concentrations ence of the Kirindy River, which runs through both habitats, for each study site. Depicted are medians (black horizontal lines) however, should have a similar impact on tree phenology at and interquartile ranges (boxes). both sites. Nevertheless, the higher fruit feeding rates of the interior groups may suggest that they may have enjoyed a higher en- ergy intake. Reduced energy intake via reduced fruit con- sumption is often associated with an increase in glucocor- ticoid output, indicating nutritional stress (Creel et al., 2002; Homan et al., 2003; Pride, 2005; Dunn et al., 2013). In chimpanzees (Pan troglodytes), for example, urinary corti- sol levels have been found to increase when they feed less on fruits (Muller and Wrangham, 2004). Similarly, fGCM lev- els have been reported to be negatively correlated with habi- tat quality in Diademed sifakas (P. diadema) as well as in both Diademed and Verreaux’s sifakas during periods of fruit scarcity (Tecot et al., 2019; Rudolph et al., 2020). Moreover, energetically costly travelling rates are associated with an in- crease in fGCM levels in mantled howler monkeys (Alouatta palliata; Gómez-Espinosa et al., 2014). However, daily travel distances and variation in home range size did not affect fGCM levels in Verreaux’s sifakas living in the forest interior (Rudolph et al., 2019). Hence, increased fGCM levels found Figure 7. Probability (yes or no) of giving birth for females living in the present study might be due to the lower energy intake, in the edge or forest interior habitat. i.e. the lower rate of fruit consumption, by sifakas inhabiting the edge habitat. In contrast to faecal glucocorticoid measurements, hair have been of lesser quality, despite the potential benefit of cortisol measurements integrate changes in glucocorticoid additional light and plant productivity (Ganzhorn, 1992). output over periods of weeks and months and may provide Interestingly, the closely related Coquerel’s sifakas in dry information on chronically increased stress hormone levels. deciduous forests in Ankarafantsika National Park exhibited Squirrel gliders (Petaurus norfolcensis) have been found to the opposite pattern, with edge groups inhabiting larger home exhibit higher hair cortisol levels in urbanized edges than in ranges than groups in the forest interior (McGoogan et al., intact habitats (Brearley et al., 2012). However, edge-positive 2009). Activity budgets and food quality did not differ be- grey mouse lemurs (Microcebus murinus) and edge-negative tween those habitats, but larger home ranges of edge groups fat-tailed dwarf lemurs (Cheirogaleus medius), which were might have been due to the fact that these Coquerel’s sifakas studied in the same edge and interior forest habitat as in our did not use the immediate forest edge and instead roamed study, did not differ in hair cortisol levels and parasite burden about 400 m away from it. Similar, Diademed sifakas in frag- (Rakotoniaina et al., 2016). As glucocorticoid concentrations mented habitats avoided the edges, whereas those in continu- in sifakas in our study were measured in faecal samples re- https://doi.org/10.5194/pb-8-1-2021 Primate Biol., 8, 1–13, 2021
10 K. Dinter et al.: Behavioural and physiological responses of Verreaux’s sifakas flecting stress hormone output over a shorter timescale (Fich- 5 Conclusions tel et al., 2007), our fGCM measure might more likely reflect a response to short-term variation in diet rather than a general Although Verreaux’s sifakas in Kirindy survive and repro- adaptation to living in the edge habitat. Hence, measurements duce in the edge habitat, they appear to be sensitive to vari- of sifaka hair cortisol are required to confirm this assump- ation in microhabitat conditions that may be linked to the tion. In this respect, in eastern chipmunks (Tamias striatus) forest edge. However, more information on habitat character- individuals from logged habitats have been found to have istics in the edge and interior forest is required to test this no- higher faecal glucocorticoid metabolite levels than individu- tion. As the forest cover in the Kirindy Forest region contin- als from intact habitats, but hair cortisol levels did not differ ues to decline (Zinner et al., 2014), and habitat fragmentation between these groups (Mastromonaco et al., 2014). As chip- into smaller forest patches leads to a decrease in the interior- munks could easily move between logged and undisturbed to-edge ratio, the results of this study – indicating that Ver- areas and faecal samples were collected from the periph- reaux’s sifakas might show a negative edge response – offer ery of the logged area, the increase in faecal glucocorticoid relevant information for future assessments of the conserva- metabolite levels might rather reflect short-term increases in tion status of this species. cortisol, whereas the similar hair cortisol levels of individu- als in both areas may indicate that habitat disturbance may not have affected long-term stress hormone output. Hence, Data availability. Data are available in the Supplement. depending on the glucocorticoid measure taken, i.e. from fae- ces or hair, the temporal scale of the sample collection should be taken into account (Mastromonaco et al., 2014). Supplement. The supplement related to this article is available online at: https://doi.org/10.5194/pb-8-1-2021-supplement. At the ultimate level, elevated glucocorticoid concentra- tions are often associated with impairment of gonadal func- tion and reproductive behaviour in female and male verte- Author contributions. KD and CF conceived the study and anal- brates (Cameron, 1997; Hardy et al., 2005) and are linked ysed the data. KD collected behavioural data and faecal samples. to reduced reproductive output in a variety of species (e.g. MH supervised the hormone analysis. KD, CF, MH, and PMK Ellenberg et al., 2007; Buck et al., 2007; Dantzer et al., wrote the paper. 2014; Vitousek et al., 2018). For example, in northern saw- whet owl (Aegolius acadicus) habitat loss and fragmentation increased physiological stress levels in males and resulted Competing interests. The authors declare that they have no con- in lower reproductive success (Hinam and St. Clair, 2008). flict of interest. Here, our data suggest that such a link between elevated glu- cocorticoid levels and decreased reproductive success may also exist for edge-living Verreaux’s sifakas, which, accord- Acknowledgements. We thank the Département de Biologie An- ing to our long-term records, show significantly reduced birth imale, Université d’Antananarivo; the CAFF/CORE of the Direc- rates when compared with sifakas living in the forest inte- tion des Eaux et Forêts; and the CNFEREF Morondava for the au- rior. However, the small number of animals studied for stress thorization of this study. We are extremely thankful to the DPZ staff hormone output and the fact that our faecal measure of glu- at Kirindy and Morondava for support and companionship in the field. Special thanks go to Remy d’Ampataka, Tianasoa Andrian- cocorticoid output may not fully reflect the effect of long- janahary, and Jean-Pierre Tolojanahary for their support in the field term stress precludes us from drawing firm conclusions. Fu- and to Andrea Heistermann for technical assistance with the hor- ture studies including more animals possibly in other edge- mone analysis. habitat locations and GCM analyses from faeces over a more extended period or even cortisol analyses from hair should be carried out to substantiate our preliminary reasoning. These Financial support. This research has been supported by a grant future studies may also try to collect comparative data on fe- from the “Ethologische Gesellschaft” to Klara Dinter. male ovarian activity for assessing ovulation and conception rates, reproductive behaviour, and possibly energetic condi- tion to help uncover the possible mechanisms underlying a Review statement. This paper was edited by Ute Radespiel and potential link between stress hormone output and reproduc- reviewed by two anonymous referees. tive function and success in Verreaux’s sifaka populations living in different habitat types. This, in turn, may help im- prove the conservation efforts for this highly endangered pri- mate species. Primate Biol., 8, 1–13, 2021 https://doi.org/10.5194/pb-8-1-2021
K. Dinter et al.: Behavioural and physiological responses of Verreaux’s sifakas 11 References ticoid levels in free-ranging howlers (Alouatta palliata), Int. J. Primatol., 34, 246–259, 2013. Ellenberg, U., Setiawan, A. N., Cree, A., Houston, D. M., and Sed- don, P. J.: Elevated hormonal stress response and reduced repro- Altmann, J.: Observational study of behavior: sampling methods, ductive output in Yellow-eyed penguins exposed to unregulated Behaviour, 49, 227–67, 1974. tourism, Gen. Comp. Endocr., 152, 54–63, 2007. Andriatsitohaina, B., Ramsay, M. S., Kiene, F., Lehman, S. M., Ra- Eppley, T. M., Santini, L., Tinsman, J. C., and Donati, G.: Do soloharijaona, S., Rakotondravony, R., and Radespiel, U.: Eco- functional traits offset the effects of fragmentation? The case of logical fragmentation effects in mouse lemurs and small mam- large-bodied diurnal lemur species, Am. J. Primatol., 82, e23104, mals in northwestern Madagascar, Am. J. Primatol., e23059, https://doi.org/10.1002/ajp.23104, 2020. https://doi.org/10.1002/ajp.23059, 2019. Fichtel, C., Kraus, C., Ganswindt, A., and Heistermann, M.: Influ- Baayen, R. H.: Analyzing linguistic data: A practical introduction to ence of reproductive season and rank on faecal glucocorticoid statistics using R, Cambridge University Press, Cambridge, UK, levels in free-ranging male Verreaux’s sifakas (Propithecus ver- 2008. reauxi), Horm. Behav., 51, 640–648, 2007. Baden, A. L., Mancini, A. N., Federman, S., Holmes, S. M., John- Gabriel, D. N.: Habitat use and activity patterns as an indication of son, S. E., Kamilar, J., Louis Jr., E. E., and Bradley, B. J.: An- fragment quality in a strepsirrhine primate, Int. J. Primatol., 34, thropogenic pressures drive population genetic structuring across 388–406, 2013. a critically endangered lemur species range, Sci. Rep.-UK, 9, Ganzhorn, J. U.: Leaf chemistry and the biomass of folivorous pri- 16276, https://doi.org/10.1038/s41598-019-52689-2, 2019. mates in tropical forests, Oecologia, 91, 540–547, 1992. Balestri, M., Barresi, M., Campera, M., Serra, V., Ramanaman- Gómez-Espinosa, E., Rangel-Negrín, A., Chavira, R., Canales- jato, J. B., Heistermann, M., and Donati, G.: Habitat degrada- Espinosa, D., and Dias, P. A. D.: The effect of energetic and psy- tion and seasonality affect physiological stress levels of Eule- chosocial stressors on glucocorticoids in mantled howler mon- mur collaris in littoral forest fragments, PloS One, 9, e107698, keys (Alouatta palliata), Am. J. Primatol., 76, 362–373, 2014. https://doi.org/10.1371/journal.pone.0107698, 2014. Gould, L. and Cowen, L. L.: Lemur catta in small forest fragments: Bates, D., Maechler, M., Bolker, B., and Walker, S.: Fitting linear Which variables best predict population viability?, Am. J. Prima- mixed-effects models using lme4, J. Stat. Softw., 67, 1–48, 2015. tol., 82, e23095, https://doi.org/10.1002/ajp.23095, 2020. Benadi, G., Fichtel, C., and Kappeler, P. M.: Intergroup relations Hardy, M. P., Gao, H., Dong, Q., Ge, R., Wang, Q., Chai, W. R., and home range use in Verreaux’s Sifaka (Propithecus ver- Feng, X., and Sottas, C.: Stress hormone and male reproductive reauxi), Am. J. Primatol., 70, 956–965, 2008. function, Cell Tissue Res., 322, 147–153, 2005. Brearley, G., McAlpine, C., Bell, S., and Bradley, A.: Influence of Harper, G. J., Steininger, M. K., Tucker, C. J., Juhn, D., and urban edges on stress in an arboreal mammal: a case study of Hawkins, F.: Fifty years of deforestation and forest fragmenta- squirrel gliders in southeast Queensland, Australia, Landscape tion in Madagascar, Environ. Conserv., 34, 1–9, 2007. Ecol., 27, 1407–1419, 2012. Heistermann, M., Ademmer, C. H. R., and Kaumanns, W.: Ovarian Buck, C. L., O’Reilly, K. M., and Kildaw, S. D.: Interannual vari- cycle and effect of social changes on adrenal and ovarian function ability of Black-legged Kittiwake productivity is reflected in in Pygatherix nemaeus, Int. J. Primatol., 25, 689–708, 2004. baseline plasma corticosterone, Gen. Comp. Endocr., 150, 430– Heistermann, M., Palme, R., and Ganswindt, A.: Comparison of dif- 436, 2007. ferent enzyme-immunoassays for assessment of adrenocortical Burgman, M. A. and Fox, J. C.: Bias in species range estimates activity in primates based on faecal analysis, Am. J. Primatol., from minimum convex polygons: implications for conservation 68, 257–273, 2006. and options for improved planning, Anim. Conserv., 6, 19–28, Hinam, H. L. and St. Clair, C. C. S.: High levels of habitat loss 2003. and fragmentation limit reproductive success by reducing home Cameron, J. L.: Stress and behaviourally induced reproductive dys- range size and provisioning rates of Northern saw-whet owls, function in primates, Semin. Reprod. Endocr., 15, 37–45, 1997. Biol. Conserv., 141, 524–535, 2008. Carrai, V. and Lunardini, A.: Activity patterns and home range use Homan, R. N., Regosin, J. V., Rodrigues, D. M., Reed, J. M., Wind- of two groups of Propithecus v. verreauxi in the Kirindy Forest, miller, B. S., and Romero, L. M.: Impacts of varying habitat Primate Report, 46, 275–284, 1996. quality on the physiological stress of spotted salamanders (Am- Creel, S., Fox, J. E., Hardy, A., Sands, J., Garrott, B., and Peterson, bystoma maculatum), Anim. Conserv., 6, 11–18, 2003. R. O.: Snowmobile activity and glucocorticoid stress responses Irwin, M. T.: Diademed sifaka (Propithecus diadema) ranging and in wolves and elk, Conserv. Biol., 16, 809–814, 2002. habitat use in continuous and fragmented forest: higher density Dantzer, B., Fletcher, Q. E., Boonstra, R., and Sheriff, M. by lower viability in fragments?, Biotropica, 40, 231–240, 2008. J.: Measures of physiological stress: a transparent or Irwin, M. T., Junge, R. E., Raharison, J.-L., and Samonds, K. E.: opaque window into the status, management and con- Variation in physiological health of diademed sifakas across in- servation of species?, Conserv. Physiol., 2, cou023, tact and fragmented forest at Tsinjoarivo, Eastern Madagascar, https://doi.org/10.1093/conphys/cou023, 2014. Am. J. Primatol., 72, 1013–1025, 2010. Delattre, P., Morellet, N., Codreanu, P., Miot, S., Quéré, J.-P., Irwin, M. T., Samonds, K. E., Raharison, J. L., Junge, R. E., Ma- Sennedot, F., and Baudry, J.: Influence of edge effects on com- hefarisoa, K. L., Rasambainarivo, F., Godfrey, L. R., and Glan- mon vole population abundance in an agricultural landscape of der, K. E.: Morphometric signals of population decline in di- eastern France, Acta Theriol., 54, 51–60, 2009. ademed sifakas occupying degraded rainforest habitat in Mada- Dunn, J. C., Cristóbal-Azkarate, J., Schulte-Herbrüggen, B., gascar, Sci. Rep., 9, 1–14, 2019. Chavira, R., and Veà, J. J.: Travel time predicts fecal glucocor- https://doi.org/10.5194/pb-8-1-2021 Primate Biol., 8, 1–13, 2021
12 K. Dinter et al.: Behavioural and physiological responses of Verreaux’s sifakas IUCN: The IUCN Red List of Threatened Species, Version 2020-2, Muller, M. N. and Wrangham, R. W.: Dominance, cortisol and available at: https://www.iucnredlist.org, last access: 9 July 2020. stress in wild chimpanzees (Pan troglodytes schweinfurthii), Be- Jaimez, N. A., Bribiescas, R. G., Aronsen, G. P., Anestis, S. A., and hav. Ecol. Sociobiol., 55, 332–340, 2004. Watts, D. P.: Urinary cortisol levels of gray-cheeked mangabeys Murcia, C.: Edge effects in fragmented forests: implications for are higher in disturbed compared to undisturbed forest areas in conservation, TREE, 10, 58–62, 1995. Kibale National Park, Uganda, Anim. Conserv., 15, 242–247, Norscia, I., Carrai, V., and Borgognini-Tarli, S. M.: Influence of dry 2012. season and food quality and quantity on behavior and feeding Kalbitzer, U. and Heistermann, M.: Long-term storage effects in strategy of Propithecus verreauxi in Kirindy, Madagascar, Int. J. steroid metabolite extracts from baboon (Papio sp.) faeces – a Primatol., 27, 1001–1022, 2006. comparison of three commonly applied storage methods, Meth- Pride, E.: High faecal glucocorticoid levels predict mortality in ring- ods Ecol. Evol., 4, 493–500, 2013. tailed lemurs (Lemur catta), Biol. Letters, 1, 60–63, 2005. Kappeler, P. M. and Fichtel, C.: A 15-year perspective on the so- Radespiel, U., Schulte, J., Burke, R. J., and Lehman, S. M.: Molec- cial organization and life history of Sifaka in Kirindy Forest, in: ular edge effects in the Endangered golden-brown mouse lemur, Long-Term Field Studies of Primates, edited by: Kappeler, P. M. Microcebus ravelobensis, Oryx, 53, 716–726, 2019. and Watts, D. P., Springer, Berlin and Heidelberg, Germany, 101– Rakotoniaina, J. H., Kappeler, P. M., Ravoniarimbinina, P., Pec- 121, 2012. houskova, E., Hämäläinen, A. M., Grass, J., Kirschbaum, C., and Koch, F., Signer, J., Kappeler, P. M., and Fichtel, C.: Intergroup Kraus, C.: Does habitat disturbance affect stress, body condition encounters in Verreaux’s sifakas (Propithecus verreauxi): who and parasitism in two sympatric lemurs?, Conserv. Physiol., 4, fights and why?, Behav. Ecol. Sociobiol., 70, 797–808, 2016. cow034, https://doi.org/10.1093/conphys/cow034, 2016. Koch, F., Ganzhorn, J. U., Rothman, J. M., Chapman, C. A., and Rangel-Negrín, A., Alfaro, J. L., Valdez, R. A., Romano M. C., and Fichtel, C.: Sex and seasonal differences in diet and nutrient in- Serio-Silva, J. C.: Stress in Yucatan spider monkeys: effects of take in Verreaux’s sifakas (Propithecus verreauxi), Am. J. Prima- environmental conditions on faecal cortisol levels in wild and tol., 79, e22595, https://doi.org/10.1002/ajp.22595, 2017. captive populations, Anim. Conserv., 12, 496–502, 2009. Laurance, W. F.: Reflections on the tropical deforestation crisis, Rangel-Negrín, A., Coyohua-Fuentes, A., Chavira, R., Canales- Biol. Conserv., 91, 109–117, 1999. Espinosa, D., and Dias, P. A.: Primates living outside pro- Laurance, W. F., Delamônica, P., Laurance, S., Delamônica, tected habitats are more stressed: the case of black howler P., Laurance, S. G., Vasconcelos, H. L., and Thomas, E.: monkeys in the Yucatán Peninsula, PloS one, 9, e112329, Rainforest fragmentation kills big trees, Nature, 404, 836, https://doi.org/10.1371/journal.pone.0112329, 2014. https://doi.org/10.1038/35009032, 2000. Richard, A. F.: Behavioral variation: Case study of a Malagasy Laurance, W. F., Camargo, J. L., Luizão, R. C., Laurance, S. G., lemur, Bucknell University Press, Lewisburg, PA, USA, 1978. Pimm, S. L., Bruna, E. M., and Van Houtan, K. S.: The fate of Rimbach, R., Link, A., Heistermann, M., Gomez-Posada, C., Amazonian forest fragments: a 32-year investigation, Biol. Con- and Heymann, E. W.: Effects of logging, hunting, and for- serv., 144, 56–67, 2011. est fragment size on physiological stress levels of two sym- Lehman, S., Rajaonson, A., and Day, S.: Lemur responses to edge patric ateline primates in Colombia, Conserv. Physiol., 1, cot031, effects in the Vohibola III classified forest, Madagascar, Am. J. https://doi.org/10.1093/conphys/cot031, 2013a. Primatol., 68, 293–299, 2006a. Rimbach, R. Heymann, E. W., Link, A., and Heistermann, M.: Val- Lehman, S., Rajaonson, A., and Day, S.: Edge effects and their in- idation of an enzymeimmunoassay for assessing adrenocortical fluence on lemur density and distribution in Southeast Madagas- activity and evaluation of factors that affect levels of fecal gluco- car, Am. J. Phys. Anthropol., 129, 232–241, 2006b. corticoid metabolites in two New World primates, Gen. Comp. Lehman, S., Rajaonson, A., and Day, S.: Edge effects on the density Endocr., 191, 13–23, 2013b. of Cheirogaleus major, Int. J. Primatol., 27, 1569–1588, 2006c. Rudolph, K., Fichtel, C., Schneider, D., Heistermann, M., Koch, Lehtinen, R. M., Ramanamanjato, J.-B., and Raveloarison, J. G.: F., Daniel, R., and Kappeler, P. M.: One size fits all? Relation- Edge effects and extinction proneness in a herpetofauna from ships among group size, health, and ecology indicate a lack of Madagascar, Biodivers. Conserv., 12, 1357–1370, 2003. an optimal group size in a wild lemur population, Behav. Ecol. Lidicker, W. Z.: Responses of mammals to habitat edges: an Sociobiol., 73, 132, https://doi.org/10.1007/s00265-019-2746-0, overview, Landscape Ecol., 14, 333–343, 1999. 2019. Martínez-Mota, R., Valdespino, C., Sánchez-Ramos, M. A., and Rudolph, K., Fichtel, C., Heistermann, M., and Kappeler, P. M.: Dy- Serio-Silva, J. C.: Effects of forest fragmentation on the physi- namics and determinants of glucocorticoid metabolite concen- ological stress response of black howler monkeys, Anim. Con- trations in wild Verreaux’s sifakas, Horm. Behav., 124, 104760, serv., 10, 374–379, 2007. https://doi.org/10.1016/j.yhbeh.2020.104760, 2020. Mastromonaco, G. F., Gunn, K., McCurdy-Adams, H., Ed- Schwitzer, C., Mittermeier, R. A., Johnson, S. E., Donati, G., Ir- wards, D. B., and Schulte-Hostedde, A. I.: Validation and win, M., Peacock, H., Ratsimbazafy, J., Razafindramanana, J., use of hair cortisol as a measure of chronic stress in east- Louis, E. E., Chikhi, L., Colquhoun, I. C., Tinsman, J., Dolch, ern chipmunks (Tamias striatus), Conserv. Physiol., 2, cou055, R., LaFleur, M., Nash, S., Patel, E., Randrianambinina, B., Ra- https://doi.org/10.1093/conphys/cou055, 2014. solofoharivelo, T., and Wright, P. C.: Averting lemur extinctions McGoogan, K. C., Steffens, T., and Lehman, S.: Edge effects amid Madagascar’s political crisis, Science, 343, 842–843, 2014. on home range size of Propithecus coquereli in northwest Shutt, C., Setchell, J., and Heistermann, M.: Non-invasive monitor- Madagascar, Am. J. Primatol., 71, Suppl: 1–27, abstract 127, ing of physiological stress in the western lowland gorilla (Gorilla https://doi.org/10.1002/ajp.20733, 2009. gorilla gorilla): Validation of a faecal glucocorticoid assay and Primate Biol., 8, 1–13, 2021 https://doi.org/10.5194/pb-8-1-2021
K. Dinter et al.: Behavioural and physiological responses of Verreaux’s sifakas 13 methods for practical application in the field, Gen. Comp. En- Vitousek, M. N., Taff, C. C., Hallinger, K. K., Zimmer, C., and Win- docr., 179, 167–177, 2012. kler, D. W.: Hormones and fitness: evidence for trade-offs in glu- Tarlow, E. M. and Blumstein, D. T.: Evaluating methods to quan- cocorticoid regulation across contexts, Front. Ecol. Evol., 6, 42 tify anthropogenic stressors on wild animals, Appl. Anim. Behav. pp., https://doi.org/10.3389/fevo.2018.00042, 2018. Sci., 102, 429–451, 2007. Watson, J. E. M., Whittaker, R. J., and Dawson, T. P.: Habitat struc- Tecot, R. S., Irwin, M. T., and Raharison, J.-L.: Faecal gluco- ture and proximity to forest edge affect the abundance and dis- corticoid metabolite profiles in diademed sifakas increase dur- tribution of forest-dependent birds in tropical coastal forests of ing seasonal fruit scarcity with interactive effects of age/sex southeastern Madagascar, Biol. Conserv., 120, 311–327, 2004. class and habitat degradation, Conserv. Physiol., 7, coz001, Zinner, D., Wygoda, C., Razafimanantsoa, L., Rasoloarison, R., An- https://doi.org/10.1093/conphys/coz001, 2019. drianandrasana, H. T., Ganzhorn, J. U., and Torkler, F.: Analysis Trillmich, J., Fichtel, C., and Kappeler, P. M.: Coordination of group of deforestation patterns in the central Menabe, Madagascar, be- movements in wild Verraux’s sifakas (Propithecus verreauxi), tween 1973 and 2010, Reg. Environ. Change, 14, 157–166, 2014. Behaviour, 141, 1103–1120, 2004. Zuur, A. F., Ieno, E. N., Walker, N. J., Saveliev, A. A., and Smith, Vieilledent, G., Grinand, C., Rakotomalala, F. A., Ranaivosoa, R., G. M.: Mixed effects models and extensions in ecology with R, Rakotoarijaona, J. R., Allnutt, T. F., and Achard, F.: Combining Springer, New York, USA, 2009. global tree cover loss data with historical national forest cover maps to look at six decades of deforestation and forest fragmen- tation in Madagascar, Biol. Conserv., 222, 189–197, 2018. https://doi.org/10.5194/pb-8-1-2021 Primate Biol., 8, 1–13, 2021
You can also read