Effects of habitat edges on vegetation structure and the vulnerable golden-brown mouse lemur (Microcebus ravelobensis) in northwestern Madagascar

Bertrand Andriatsitohaina
Daniel Romero-Mujalli
Malcolm S Ramsay
Frederik Kiene
Solofonirina Rasoloharijaona
Romule Rakotondravony
Shawn M Lehman
Ute Radespiel
Effects of habitat edges on vegetation structure and the vulnerable golden-brown mouse lemur (Microcebus ravelobensis) in northwestern Madagascar
10.1186/s12898-020-00337-zAndriatsitohaina et al. BMC Ecol (2020) 20:69 RESEARCH ARTICLE This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article' s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article' s zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.Habitat lossEdge effectsHabitat choiceVegetation structureAbundanceMouse lemurDepth of edge influenceDEIMicrocebus ravelobensisMadagascar
Background: Edge effects can influence species composition and community structure as a result of changes in microenvironment and edaphic variables. We investigated effects of habitat edges on vegetation structure, abundance and body mass of one vulnerable Microcebus species in northwestern Madagascar. We trapped mouse lemurs along four 1000-m transects (total of 2424 trap nights) that ran perpendicular to the forest edge. We installed 16 pairs of 20 m 2 vegetation plots along each transect and measured nine vegetation parameters. To determine the responses of the vegetation and animals to an increasing distance to the edge, we tested the fit of four alternative mathematical functions (linear, power, logistic and unimodal) to the data and derived the depth of edge influence (DEI) for all parameters.Results: Logistic and unimodal functions best explained edge responses of vegetation parameters, and the logistic function performed best for abundance and body mass of M. ravelobensis. The DEI varied between 50 m (no. of seedlings, no. of liana, dbh of large trees [dbh ≥ 10 cm]) and 460 m (tree height of large trees) for the vegetation parameters, whereas it was 340 m for M. ravelobensis abundance and 390 m for body mass, corresponding best to the DEI of small tree [dbh < 10 cm] density (360 m). Small trees were significantly taller and the density of seedlings was higher in the interior than in the edge habitat. However, there was no significant difference in M. ravelobensis abundance and body mass between interior and edge habitats, suggesting that M. ravelobensis did not show a strong edge response in the study region. Finally, regression analyses revealed three negative (species abundance and three vegetation parameters) and two positive relationships (body mass and two vegetation parameters), suggesting an impact of vegetation structure on M. ravelobensis which may be partly independent of edge effects.Conclusions:A comparison of our results with previous findings reveals that edge effects are variable in space in a small nocturnal primate from Madagascar. Such an ecological plasticity could be extremely relevant for mitigating species responses to habitat loss and anthropogenic disturbances.

Introduction

Forest loss leading to fragmentation is the main consequence of human activities in forest ecosystems. The will diverge more over time (in species composition and possibly ecosystem functioning) than sites from the same landscape. However, it is not known if this variability in edge responses in time and space occurs in other continents and phyla.

Many studies in different biomes and continents have shown that the depth of edge influence (DEI) can differ largely between parameters [15]. The DEI is defined as the distance at which biotic or biotic parameters change considerably from the edge to the interior of the forest. However, determining the DEI is not a simple task, as habitat edges may affect plants and animals in various non-linear ways [16]. Previous authors used different qualitative and quantitative approaches to estimate the DEI, e.g., [15,[17][18][19][20]. Ewers and Didham formalized the first mathematical method to determine the DEI, which is based on a relative comparison of the accuracy of different mathematical functions to describe continuous response functions for any biotic or abiotic variable across ecological boundaries [21] and shall be employed in this study on a small Malagasy lemur species.

Madagascar is one of the world's most important biodiversity hotspots and is experiencing widespread changes in forest cover due to climate change, forest loss, and habitat fragmentation [22]. Forest loss has resulted in a reduction of the island's forest cover by 44% between 1953 and 2014 [23]. Remaining forests exist primarily as landscapes of forest fragments of varying size and shape surrounded by a matrix such as grasslands or agricultural fields. These forest landscapes represent an ideal natural laboratory for studying edge effects on plants and animals, providing an important opportunity to model extreme levels of deforestation processes. In turn, results from research on Madagascar can inform similar work in other tropical regions not yet subject to the same extremes of forest loss and fragmentation. Clearly, understanding how plants and animals respond to edge effects is crucial for conservation planning [24].

It was already shown in some Malagasy rainforest sites that edge effects can affect vegetation structure and the abundance and distribution of animal species such as birds or lemurs e.g., [25][26][27][28]. Extant lemurs are a highly threatened mammalian clade that comprises more than 100 described taxa which are endemic to Madagascar [29,30]. Studies from rainforest habitats described three types of edge responses in lemurs: edge tolerance, neutral, or negative edge responses, e.g., [26][27][28]. Moreover, these studies noted that tree height and dbh were smaller at the edge compared to the interior of forest [28].

However, there are only few data on how edge effects influence the vegetation in the tropical dry forests in western Madagascar [31,32] and also on how animals respond to edge effects in these forests [17]. In the dry deciduous forests of northwestern Madagascar, the golden-brown mouse lemur (M. ravelobensis) was found to be more abundant in edge habitats than in the forest interior in one part of the Ankarafantsika National Park; this was not the case for the grey mouse lemur (M. murinus) which lives in sympatry with golden-brown mouse lemurs at this site [33]. Moreover, female M. ravelobensis weighed more in edge habitats than in the interior, which was not observed in M. murinus. At the same study site, a higher density of stems and seedlings was found in the forest interior than close to forest edges, but there was no difference of tree dbh between the two habitats [31,32]. Larger scale landscape comparisons of edge effects within the same biome, however, are still missing from Madagascar. They would allow more generalized insights into edge effects, which could result in more nuanced conservation management suggestions.

Our study focused on the vulnerable M. ravelobensis which has a limited geographic distribution in the dry deciduous forests of northwestern Madagascar [34,35]. M. ravelobensis is a nocturnal, small-bodied (~ 60 g), arboreal, omnivorous solitary forager and occurs in partial sympatry with M. murinus in some but not all forests of northwestern Madagascar [36,37]. M. ravelobensis forms mixed-sex sleeping groups and sleeps in tree-holes but also in dense tangles of lianas, leaves or self-built nests [38,39]. Previous studies reported certain microhabitat preferences of this species, such as habitats containing higher density of trees with many lianas and a higher cover of the herb layer, and habitats closer to surface water or at lower altitudes [36,37,40,41]. Most recently, significant regional differences were revealed in the abundance of M. ravelobensis across its distribution [37]. Whether these differences correspond to differences in edge responses, however, is not yet known.

The aims of this study were to use ecological modelling techniques to evaluate edge effects (1) in vegetation structure, (2) on the distribution and body mass of M. ravelobensis, and (3) on the spatial variability of these parameters based on comparisons with previous studies from a different study region in northwestern Madagascar.

In particular, we addressed the following questions:

Methods

Study site and period

This study was conducted in the Mariarano classified forest (MCF), located near the coast (15° 24′ S, 46° 44′ E, 80 km northeast of Mahajanga) and next to the village Mariarano at an altitude of 20-90 m a.s.l. (Fig. 1)

Data collection

We established four transects of 1000 m length each of which ran perpendicular to the forest edge. These transects bordered different matrix types: road with bushes across (C1), savannah with low density of palm trees (Bismarckia nobilis or Satrana as local name) (C2), rice fields with mango trees (C3), and savannah with a high density of palm trees (B. nobilis) (C4). Mouse lemurs were captured systematically three times over a 2-week period along each transect by means of Sherman traps that were baited with banana and installed in pairs every 10 m (one trap on either side of the transect) at 0.5 -2.0 m height above ground. A total of 2424 trap nights (202 traps per trapping session per transect) were conducted across all four transects. After the identification of species by means of phenotypic and morphometric characterization (following [34]), all captured individuals of M. ravelobensis were sexed, aged (juvenile/adults), individually marked (by ear biopsies), and weighed. All handling procedures were quick and therefore performed without the use of anesthesia as in many other studies before [e.g., [33][34][35]. We classified an individual as adult if the body mass was higher or equal to 41 g and in view of some other morphometric parameters (e.g. lower leg length, hindfoot length). Finally, all individuals were released at dusk at the same place where they had been trapped on the day of capture. We determined abundances of M. ravelobensis across each transect by calculating the number of individuals caught (without recapture) per 100 m transect length (= 20 trap locations). The data from 2017 and 2018 were combined in all analyses given the considerable seasonal overlap between both years, the lack of systematic change in trappability of mouse lemurs across the dry season in northwestern Madagascar [36], and consistent trapping results across several years along two of the transects (C1, C3, Radespiel unpubl. data). A total of 226 individual M. ravelobensis (100 males, 126 females) were trapped across the study period. In contrast, only one individual of M. murinus was trapped on one transect (C2) during the entire study period (Fig. 1). We installed 16 pairs of vegetation plots along each transect (2 m × 10 m, W x L), oriented perpendicularly to either side of a transect at pre-set distances from the forest edge to the interior (at 0 [edge], 20, 40, 60, 80, 100, 200, 250, 300, 400, 500, 600, 700, 800, 900, and 1000 m) according to [31]. In each plot, vegetation measurements focused on dendrometrics and estimated densities of woody stems by counting the number of seedlings (1-100 cm in height), saplings (101-250 cm in height), trees (height > 250 cm) and lianas (with dbh ≥ 2.5 cm). We measured the height and dbh of all trees rooted in a plot. From these data, we calculated the total number of trees (categorized into small: dbh < 10 cm, and large: dbh ≥ 10 cm), seedlings, saplings and lianas per double plot (40 m 2 ) at each sampling distance from the edge for each transect and then calculated the overall mean density across the four transects. Similarly, we calculated the mean height and dbh of small and large trees (per 40 m 2 ) at each sampling distance across all four transects. These mean values were used to determine the DEI according to methods in Ewers and Didham [21].

We calculated the mean vegetation parameters per 20 m 2 plot at each distance from edge for each transect (= four mean values for each sampling point) for the comparison of habitat parameters between the edge habitat and the interior of the forest. The values for all vegetation parameters were also summarized as means/100 m section along each transect to prepare the vegetation dataset for comparisons with the animal dataset.

Data analyses

Following Ewers and Didham [21] we calculated and compared the fit of four different mathematical models, linear (Eq. 1), power (Eq. 2), logistic (Eq. 3), and unimodal (Eq. 4) to the data to evaluate how the relationship between the distance to the edge and the test parameters could be best formalized. Support for the linear model would indicate a lack of a finite impact of the edge whereas the three other models can be used to determine the depth of edge influence (DEI). The chosen test parameters were:

-vegetation parameters: density, dbh and height of small trees (dbh < 10 cm); density, dbh and height of large trees (dbh ≥ 10 cm); density of saplings; density of seedlings and density of lianas. -Animal parameters: M. ravelobensis abundance and body mass.

Based on Ewers and Didham, the formulas for the four used mathematical models were:
• Linear • Power • Logistic • Unimodal
where n is the mean of given response variable (see above), Y max and Y min are the maximum and minimum values of the response variable, β 1 , β 2 and β 3 are constants, x is the distance to edge, and ε is an error term.

The process was repeated 1000 times in R (version 3.6.0). The best fit of each of the four mathematical models was compared for each parameter using an information theoretic approach and calculated Akaike's Information Criterion with a correction for small sample size (AICc) and R 2 , with lowest value of AICc and the highest R 2 indicating the best model [43,44].

We determined the DEI (defined as the edge distance over which changes in the response variable can be detected) by calculating the first and second derivative of the best model (for functions see above). The first derivative describes the local minimum or maximum of the rate of change in the response variable. The second derivative enabled calculation of the DEI as the distance between the two first inflection points in the curve describing the parameter under study [21,45]. Based on the second derivative curve, we used the "inflection" package (vs. 1.3.5) and employed the Extremum Surface Estimator ESE method by means of the findipiterplot() command in R (version 3.6.0) to determine the distance between the two first inflection points in the curve, i.e., the DEI [46].
(1) n x = β 1 + β 2 x + ε (2) n x = β 1 * e β2 * x + ε (3) n x = Y min + Y max − Y min / 1 + e (β1 −x)β2 + ε (4) n x = Y min + Y max − Y min / 1 + e (β1− x +βx2)β2 3 + ε
We subsequently compared the data for each vegetation parameter (one value per double plot per transect) between respective "edge habitat" and the "interior habitat" statistically by means of generalized linear models (GLMs) using the R package "lme4" (version 1. . Similarly, the abundance (one value per 100 m per transect) and adult body mass of M. ravelobensis captured in "edge habitat" and "interior habitat" of the forest were compared by means of GLMs. The "edge habitat" and "interior habitat" were determined for each vegetation and animal parameter separately based on the respective DEI which was "buffered" to both sides. "Buffering" was achieved by first considering the respective DEI for each parameter and then excluding the data-point closest to that DEI in both directions (i.e., towards edge and towards interior) to account for the disjunct spatial data collection scheme in vegetation plots and the representation of all animal data in fixed 100 m segments.

Next, the DEI estimates derived for the different vegetation parameters were compared to the DEI inferred for M. ravelobensis abundance and adult body mass to explore which vegetation parameter corresponds best to edge-related variations in animal abundance and body mass. To evaluate this potential relationship between animals and vegetation structure further, we calculated one value per 100 m per transect for each animal and vegetation parameter, and performed regression analyses using lm() function with M. ravelobensis abundance and adult body mass as dependent variables and the vegetation parameters on the other side as independent variables.

For all analyses, normality of all variables was tested using the Shapiro-Wilk's test [47]. If a variable was not normally distributed, it was log-transformed (number of large trees, number of saplings, number of seedlings and number of lianas for regression analysis) and Q-Q plots visualization of the data were used before and after transformation to ensure improvement. Significance level was set at α = 0.05 and a statistical trend was indicated by 0.1 > p > 0.05.

Results

Effects of habitat edges on vegetation structure

In total, we estimated the density and size of 26,019 seedlings; 3114 saplings; 489 lianas; 3143 small trees; and 139 large trees in the 128 installed plots (Additional file 1: Table S1). Across the nine vegetation variables, the logistic function best explained the edge responses in five vegetation variables (seedlings, saplings, dbh of large trees, number of small trees, and height of large trees), whereas the unimodal function best explained the variation in the remaining four variables (number of lianas, number of large trees, height of small trees and dbh of small trees, Additional file 1: Table S2 and Fig. 2). The power function and the linear function did not explain edge responses in any variable. Although R 2 -values of the best models varied considerably (range 0.017-0.549), they explained more than 10% of the variation in the data for seven of nine vegetation parameters (Fig. 2a-f, h), but not for the number and dbh of large trees (R 2 < 0.1, Fig. 2g, i, Additional file 1: Table S2). The DEI, derived from the second derivative of the best models, varied across vegetation variables, ranging from 50 m (no. of seedlings, no. of liana, dbh of large trees) to 460 m (height of large trees) (Additional file 1: Table S2, Fig. 2).

The comparison between the interior and edge habitats for each vegetation parameter revealed only two significant differences (Additional file 1: Table S1, Additional file 2). First, small trees were significantly taller in the interior than in the edge (Estimate = 0.391, p = 0.032, Fig. 3a). Second, seedlings had a significantly higher density in the interior compared to edge habitat (Estimate = 0.573, p = 0.049, Fig. 3b).

Effects of edge habitats on M. ravelobensis abundance and body mass

The edge-dependent variation in the abundance of M. ravelobensis was best explained by a logistic function (R 2 = 0.283, Fig. 4a, Suppl. Information 1 Tables S2, S3). The DEI for M. ravelobensis abundance was 340 m (Fig. 4a, c, e). Spatial variations in the body mass of M. ravelobensis were also best explained by the logistic function (R 2 = 0.352, Fig. 4b, Additional file 1: Table S2, S3) and the DEI for changes in body mass was 390 m (Fig. 4b,  d, f ).

Overall, there was no significant difference in M. ravelobensis abundance between the interior and edge habitats. However, as a statistical trend, males but not females had a higher abundance in the interior than in the edge (Estimate males = 1.083, n edge = 22, n interior = 70, p = 0.082, Estimate females = − 0.039, n edge = 34, n interior = 74, p = 0.843, Fig. 5a, Additional file 1: Table S3, Additional file 2).

There was no significant difference in the overall body mass of M. ravelobensis between the edge and interior habitats, but, as a statistical trend, females but not males weighed more in edge than in interior  Table S3, Additional file 2).

Relationship between edge effects in vegetation and M. ravelobensis abundance and body mass

Based on the comparison of the DEIs derived for the vegetation parameters and the abundance and body mass  The potential relationship between the vegetation parameters and the abundance and body mass of M. ravelobensis, respectively, was analyzed by means of regression analyses (Additional file 1: Table S4, Additional file 3). These analyses revealed that the abundance of M. ravelobensis was best explained by the density of large trees (R 2 = 0.279, n = 40, p < 0.001, Fig. 6a), but also, to a lesser extent, by the density of small trees (R 2 = 0.141, n = 40, p = 0.017, Fig. 6b) and the density of saplings (R 2 = 0.138, n = 40, p = 0.019, Fig. 6c). In all cases, the relationship between the variables was negative, i.e. the abundance of M. ravelobensis decreased with increasing density of large and small trees and with increasing density of saplings. None of the other variables explained variations in abundance. Significantly positive relationships were found between body mass and the dbh of small trees (R 2 = 0.133, n = 40, p = 0.023, Fig. 7a) and the density of seedlings (R 2 = 0.108, n = 40, p = 0.042, Fig. 7b, Additional file 1: Table S4).

Discussion

Edge effects on vegetation structure

Whereas edge effects for various vegetation parameters are reported quite regularly from different regions worldwide [e.g. 7,8,16], the depth of edge influence (DEI) is less often inferred statistically. A DEI was determined for all nine vegetation parameters in this study and was well supported in seven of nine variables (R 2 > 0.1). In line with previous work in other ecosystems [17,18,48], the response pattern as well as the penetration depth varied considerably between parameters with the DEI varying from 50 m (for the density of seedlings, liana and dbh of large trees) to 460 m (height of large trees) into the forest. When comparing vegetation parameters between the edge and interior habitats, though, only two were significantly different. Specifically, the density of seedlings and the height of small trees were higher in the interior than in the edge habitats. This result was partly similar to the findings in the dry deciduous forest of Ankarafantsika NP [31], in the humid forests in Vohibola III in southeast Madagascar [26] and in a tropical rain forest in central Amazon forest where the density of seedlings increased towards the forest interior [9,49]. These results may relate to the vulnerability of edge habitats to elevated levels of desiccation and wind speed, which can cause high rates of tree fall and tree mortality along edge habitats, which in turn can affect forest structure and composition [14]. For example, in the Gatún Lake central Panama topical forest, germination and seedling survival were lower due to wind-exposure in edge habitats compared to protected sites in the forest interior [50]. In addition, previous studies examined how fragmentation impacted the susceptibility of a fragmented boreal forest to wind damage. They found that wind speed negatively affected tree height near the edge [51,52]. Thus, we suggest that the edge contrast, combined with spatial patterns of tree height and dbh, could operate synergistically to make trees more susceptible to wind damage at our research sites.

The lack of significant differences in the other seven vegetation variables between edge and interior habitats could be due to three explanations: (a) there was no systematic variation in these parameters along the transects, (b) inter-transect variations were higher than edge-interior variations, and (c) edge-interior differences exist but are only small and could not be detected due to the small sample size of four transects in our study. Whereas the first explanation may explain the poor modelling results for the number and dbh of large trees, it is unlikely for the other parameters, since substantial variation was detected in them (Fig. 2). Inter-transect variations were indeed often large (explanation (b), Fig. 2). A high variation between transects may arise if vegetation structure is not only determined by edge-interior dynamics but, for example, by varying levels of human disturbance on different transects or if different matrix types impact the adjacent forest in different ways [14]. Both aspects cannot be excluded for this study. The four transects differed in their distance to the village Mariarano, the largest settlements in the area, although smaller settlements of one to several houses were also scattered across the region (Fig. 1). Although the forest is under some basic level of protection, it is still used by the villagers for wood extraction, zebu grazing, charcoal production, or localized slash-and-burn activities to create new open spaces for agriculture and village enlargement. The resulting impacts may not have affected all transects equally. Moreover, all four transects bordered on different types of matrix, ranging from a dirt road over rice fields to savannah with different palm density (see methods). It is known that edge-related gradients in biotic variables are likely to be less pronounced when the structure of the matrix is similar to the original habitat [53], which may apply here to transects C1 (road with degraded forest on the other side of the road) and C4 (savannah with high palm density). Taken together, our study provided only limited evidence for systematic edge responses in vegetation structure and results may be best explained by local habitat divergence. Future studies with a repeated design will be needed to investigate the differential impact of the matrix on vegetation structure in Malagasy landscapes. These factors may act on both local but also on landscape levels in complex ways, as predicted by the 'Landscape-Divergence Hypothesis' [14], leading to divergent trajectories in vegetation structure and plant species assemblages even across various spatial scales.

Edge effects on animal abundance and body mass

Forest-dwelling species can react in different ways and to a different extent to the presence, depth of penetration, and intensity of forest edges [54]. Our study in Mariarano revealed a DEI of 340 m on the abundance of M. ravelobensis and a DEI of 390 m on the body mass of M. ravelobensis. Comparative values of the DEI for mouse lemurs from other sites are not yet available, but these values are in the same range as a value estimated for the larger-bodied diurnal Propithecus coquereli (~ 400 m) [17]. A visual inspection of the relationship between edge distance and abundance or body mass in M. ravelobensis revealed only moderate edge effects, and we also did not find significant differences in the overall abundance and body mass of this species between edge and interior habitats. These results would be congruent with a rather "neutral" response of M. ravelobensis to edges [4]. When analyzing these relationships separately in both sexes, however, statistical trends emerged for a higher abundance of males but not females in the interior habitat, while females but not males weighed slightly more in the edge versus the interior habitat.

The abundance results are in stark contrast to those of a previous study on M. ravelobensis which reported a higher overall female abundance at the edge compared to the interior habitat in Ankarafantsika NP, which is only about 90 km away from Mariarano [33]. Two aspects may explain these different findings: First, the two Ankarafantsika NP sites bordered on the same homogenous open savannah matrix [33], whereas the four transects in Mariarano bordered on four different matrix types, which may have impacted vegetation composition and structure differently [53] and may have resulted in a more variable abundance dataset in Mariarano (Fig. 4a). Such an interpretation would also be in agreement with the heterogeneous findings on the vegetation structure (see above). Second, the forest sites from Ankarafantsika NP contained not only M. ravelobensis but also M. murinus in variable proportions (11.3% at edge, 60.9% in the interior) which was not the case in Mariarano with just one M. murinus among 226 captured M. ravelobensis individuals across all four 1000 m transects. Burke and Lehman (2013) consequently also discussed interspecific competition as one potential driver of edge-interior utilization patterns which can be excluded for our study.

Similar to our study, females but not males weighed significantly more in the edge compared to the interior habitat in Ankarafantsika NP [33]. The authors argued that M. ravelobensis may weigh more at the edge due to an increased access to insect food that was shown to be available in higher abundance at the edge than in the interior in the same study [33]. Further studies on habitat use and feeding ecology will be needed to understand potential sex differences in body mass in this species in edge and interior habitats.

Relationship between vegetation structure, animal abundance and body mass

Various characteristics of vegetation structure have previously been suggested to determine habitat suitability for mouse lemurs and were proposed to be related to resource availability, efficient locomotion and shelter quality [55][56][57][58][59]. For example, M. ravelobensis occurred in higher abundance in sites with a high abundance of lianas, more and larger shrubs, a lower density of mediumsized trees (dbh: 5-10 cm), and overall in more humid habitats in Ankarafantsika NP [36,40,41]. Conversely, M. ravelobensis was found to be absent from forest habitats with smaller trees, a low density and low dbh of large trees, a low density of shrubs and with much open ground, characteristics that can all be directly or indirectly linked to dry conditions [36]. On the other hand, this species was also reported to possess a high microhabitat flexibility, as it was trapped in 19 out of 22 study sites (86.4%) within the continuous part of the Ankarafantsika NP [36].

When evaluating the congruence between the DEI's of M. ravelobensis abundance and adult body mass and the DEI's on the different elements of vegetation structure in this study, the animal DEI's (340 m and 390 m) corresponded best to the DEI of 360 m for the density of small trees. This relationship was further supported by a negative correlation between M. ravelobensis abundance and the number of small trees per plot, although this parameter did not differ significantly between the edge and the interior habitat. In addition, M. ravelobensis abundance increased with a decreasing number of large trees and that of saplings (Fig. 6). Conversely, the body mass of M. ravelobensis increased with increasing dbh of small trees and density of seedlings, suggesting an impact of vegetation structure on mouse lemur abundance and body mass which may be partly independent of edge effects. The underlying reasons for these relationships, however, cannot be clarified at present, since they are not in line with previous findings from Ankarafantsika NP (see above), and cannot be easily linked to food resource availability, protection from predators, or suitable shelters for daily sleeping groups [38,59,60]. One limitation of this study was the lack of floristic data. Since plant parts such as fruits, nectar or gum as well as arthropods and insect secretions do belong to the consumed food items of M. ravelobensis [61], it is possible that the local distribution of some key plant/tree species or arthropods may have impacted the local abundance and body mass of this lemur species independently of potential edge/habitat gradients. Thus, future research on the ecological patterns underlying edge effects in animal species should ideally include concordant data on the spatial distribution, abundance and use of potential food resources. In doing so, results from these studies can be used to test and refine conceptual models on how resource distribution effects abundances and body mass near edge habitats [62].

Conclusion

This study contributes to our knowledge of ecological edge effects on the vegetation structure in dry deciduous forests and revealed moderate but not very consistent depth of edge influences in various vegetation parameters. Likewise, such rather weak effects could also be identified in the abundance and body mass of the vulnerable M. ravelobensis. A comparison of our results with findings from a different study region revealed that edge effects seem to be highly variable on the local (between transects) as well as on the landscape level (between study regions) within the vegetation structure but also within the same study species, a small nocturnal primate from Madagascar. A high ecological flexibility of this lemur species towards variable ecological conditions at forest edges may explain such seemingly contradictory responses, and also its occurrence in several different types of microhabitats [36]. Such a flexibility may enable this species to cope with disturbances and some anthropogenic habitat changes within forests [63], although they seem to respond negatively to habitat fragmentation per se, probably due to a limited potential to connect between fragmented forest patches across open space [37]. Further studies are urgently needed to evaluate the critical features of landscapes that are needed to sustain viable populations of this vulnerable lemur species.

Depth of edge influences are rarely measured in a rigorous and interdisciplinary way with the simultaneous evaluation of vegetation and animal responses like attempted in this study. Edge effects are therefore not well understood for most taxonomic groups and for most parts of Madagascar, although forest loss is ongoing and omnipresent [23]. Our study detected signals of high local habitat divergence coupled with high regional landscape divergence in northwestern Madagascar. To fully understand such complex cases, even broader approaches are needed to integrate fine-scale abiotic and biotic measurements of edge effects in floristic and structural forest composition across multiple replicate sites at local and broad spatial scales. Only then the species-specific resilience towards different types of habitat disturbance, forest degradation, matrix effects, habitat loss and fragmentation can be fully disentangled. These insights are urgently needed, as they will ultimately guide conservation biologists and stakeholders to formulate efficient conservation guidelines for species such as the unique lemurs of Madagascar.

Supplementary information

Supplementary information accompanies this paper at https ://doi. org/10.1186/s1289 8-020-00337 -z.

Additional file 1: Table S1. Variation of all vegetation parameters (per 40 m 2 ) and the statistical comparison between edge and interior habitats according to their depth of edge influence (DEI) in Mariarano Classified Forest, northwestern Madagascar. Table S2. Comparison of the fit of four functions modeled for each vegetation and animal parameter in Mariarano Classified Forest, northwestern Madagascar. Table S3. Variation in abundance and body mass of all adult Microcebus ravelobensis and the comparison between edge and interior habitats after determining the respective depth of edge influence (DEI) in Mariarano Classified Forest, northwestern Madagascar. Table S4. Results of regression analyses between all nine vegetation parameters and abundance and body mass, respectively, of Microcebus ravelobensis in Mariarano Classified Forest, northwestern Madagascar.

Additional file 2. Dataset that was used for the statistical comparison of all vegetation and animal parameters (one value per double plot per transect if applicable) between the "edge habitats" and the "interior habitats".

Additional file 3. Dataset that was used for the regression analyses between all vegetation parameters and the two animal parameters abundance and body mass (one value per 100 m per transect).

Fig. 1
1Map of Mariarano Classified Forest (MCF) showing the study area and the location of all four study transects (C1, C2, C3 and C4) in relation to settlements and forest edge. The inset at the top right shows the location of the study site (star) within Madagascar (map drawn by: B. Andriatsitohaina) Andriatsitohaina et al. BMC Ecol (2020) 20:69

Fig. 2
2Illustration of the best function (logistic and unimodal) and averaged values of all nine vegetation parameters across the length of the transects (a-i) and their respective depth of edge influence DEI. Mean values (dots) are shown together with minimum and maximum values (whiskers). Lines represent the best fit of the distribution of all nine vegetation parameters across distance from edges habitats (Estimate = − 3.735, n edge = 27, n interior = 48, p = 0.052, Fig. 5b, Additional file 1:

Fig. 3
3Comparison between edge and interior habitats for the two significant vegetation parameters: height of small trees (a) and number of seedlings (b). Black squares: mean, violin plots range from minimum to maximum value with violin width illustrating data distribution

Fig. 4
4Variation of abundance (a) and body mass (b) of M. ravelobensis from the edge into the interior of forest and calculation of their DEI, respectively. Mean values (dots) are plotted together with minimum and maximum values (whiskers). In addition, the first derivatives of the fitted curves (c, d: abundance and body mass) and the second derivatives (e and f: abundance and body mass) are plotted. The DEI (vertical arrows in (e) for the abundance and in (f) for the body mass) is determined as the distance between the two first inflection points in the curve of the second derivative variations in M. ravelobensis, edge-related changes in the number of small trees corresponded best to changes in M. ravelobensis abundance and body mass (DEI small trees : 360 m, DEI abundance : 340 m, DEI body mass : 390 m).

Fig. 5
5Comparison of a abundance (males only), and b body mass (females only) of M. ravelobensis between edge and interior habitats in form of violin plots. Black squares show the mean, violin plots end with minimum and maximum valuesFig. 6 Relationship between three vegetation parameters (a-c: density of large, small trees and saplings, respectively) and the abundance of M. ravelobensis. Linear regression lines are shown for the three parameters

Fig. 7
7Relationship between three vegetation parameters (a, b: dbh of small trees and density of seedlings, respectively) and body mass of M. ravelobensis. Linear regression lines are shown for the two parameters Andriatsitohaina et al. BMC Ecol (2020) 20:69
AcknowledgementsWe would like to thank the Direction du Système des Aires Protégées (DSAP) for their permission to work in Mariarano. A large number of people have assisted with different aspects of this research and their input is very much appreciated. We are grateful to all DBCAM and Operation Wallacea staff for their support. We thank all OpWall students (Quinn Parker, Olivia Pilmore-Bedford, Harry Skinner, and Jack O'Connor) who worked with us on this project. We also thank our guides and cooks, and the people of Mariarano, Ambolodihy, and Antanambao for their constant support. We are also much endebted to our field assistants Miarisoa Ramilison, Herinjatovo Rakotondramanana, Fernand Fenomanantsoa, Onjaniaina Gilbert Razafindramasy and Leonie Baldauf; without their hard work the data would never have been collected. UR and BA also express their thanks to the late Elke Zimmermann and the Institute of Zoology for supporting the infrastructure of this project. Finally, the comments of two anonymous reviewers improved the manuscript considerably.1Authors' contributionsUR and SML conceived the study. BA, MSR and FK collected the field data. BA conducted all data analyses and drafted the manuscript. UR, SML, DRM, and RR instructed and supervised the data analyses. All authors participated in data interpretation and revised the draft manuscript. All authors read and approved the final manuscript.FundingOpen Access funding enabled and organized by Projekt DEAL. This publication was supported by the Deutsche Forschungsgemeinschaft and the University of Veterinary Medicine Hannover, Foundation within the funding Programme Open Access Publishing. This study was conducted as part of the INFRAGECO (Inference, Fragmentation, Genomics, and Conservation) project that was funded by the BiodivERsA initiative of the European Community (No. 2015-138) and more specifically by the German Bundesministerium für Bildung und Forschung (Project No. 01LC1617A). The research was also funded in part by an NSERC Discovery Grant and by Idea Wild, Primate Conservation Inc., and Sigma Xi.Availability of data and materialsAll data analysed during this study are included in the supplementary information files to this published article. The datasets analysed during the current study are also available from the corresponding author on reasonable request.Ethics approval and consent to participateAs this work included field work on animals but not on humans, the "consent to participate"-part is not applicable. All fieldwork protocols were reviewed and approved beforehand by the Institute of Zoology, University of Veterinary Medicine Hannover, Germany and the University of Toronto, Canada.Consent for publicationNot applicable.Competing interestsThe authors declare that they have no competing interests.Author details
Effects of habitat fragmentation on biodiversity. L Fahrig, 10.1146/annurev.ecolsys.34.011802.132419Annu Rev Ecol Evol Syst. 34Fahrig L. Effects of habitat fragmentation on biodiversity. Annu Rev Ecol Evol Syst. 2003;34:487-515. https ://doi.org/10.1146/annur ev.ecols ys.34.01180 2.13241 9.

Modeling temperature gradients across edges over time in a managed landscape. S C Saunders, J Chen, T D Drummer, T R Crow, For Ecol Manage. 117Saunders SC, Chen J, Drummer TD, Crow TR. Modeling temperature gra- dients across edges over time in a managed landscape. For Ecol Manage. 1999;117:17-31.

Edge effects in fragmented forests: implication for conservation. C Murcia, TREE. 10Murcia C. Edge effects in fragmented forests: implication for conserva- tion. TREE. 1995;10:58-62.

Ecological responses to habitat edges: mechanisms, models, and variability explained. L Ries, R J Fletcher, J Battin, T D Sisk, Annu Rev Ecol Evol Syst. 35Ries L, Fletcher RJ, Battin J, Sisk TD. Ecological responses to habitat edges: mechanisms, models, and variability explained. Annu Rev Ecol Evol Syst. 2004;35:491-522.

The influence of anthropogenic edge effects on primate populations and their habitat in a fragmented rainforest in Costa Rica. L M Bolt, A L Schreier, K A Voss, E A Sheehan, N L Barrickman, N P Pryor, 10.1007/s10329-018-0652-0Primates. 59Bolt LM, Schreier AL, Voss KA, Sheehan EA, Barrickman NL, Pryor NP, et al. The influence of anthropogenic edge effects on primate populations and their habitat in a fragmented rainforest in Costa Rica. Primates. 2018;59:301-11. https ://doi.org/10.1007/s1032 9-018-0652-0.

Edge effects on small mammals: differences between arboreal and ground-dwelling species living near roads in Brazilian fragmented landscapes. C A Da Rosa, H Secco, N Carvalho, A C Maia, A Bager, Austral Ecol. 43da Rosa CA, Secco H, Carvalho N, Maia AC, Bager A. Edge effects on small mammals: differences between arboreal and ground-dwelling species living near roads in Brazilian fragmented landscapes. Austral Ecol. 2018;43:117-26.

Edge-interior disparities in tree species and structural composition of the Kilengwe forest in Morogoro Region. D S Kacholi, Tanzania. Int Schol Res Not. 2014Kacholi DS. Edge-interior disparities in tree species and structural compo- sition of the Kilengwe forest in Morogoro Region, Tanzania. Int Schol Res Not. 2014;2014:1-8.

Species wood density and the location of planted seedlings drive earlystage seedling survival during tropical forest restoration. L S Charles, J M Dwyer, T J Smith, S Connors, P Marschner, M M Mayfield, J Appl Ecol. 55Charles LS, Dwyer JM, Smith TJ, Connors S, Marschner P, Mayfield MM. Species wood density and the location of planted seedlings drive early- stage seedling survival during tropical forest restoration. J Appl Ecol. 2018;55:1009-18.

Impact of forest fragmentation on seedling abundance in a tropical rain forest. J Benitez-Malvido, Conserv Biol. 12Benitez-Malvido J. Impact of forest fragmentation on seedling abun- dance in a tropical rain forest. Conserv Biol. 1998;12:380-9.

Edge effects in small forest fragments: why more is better?. M Bunyan, S Jose, R Fletcher, Am J Plant Sci. 03Bunyan M, Jose S, Fletcher R. Edge effects in small forest fragments: why more is better? Am J Plant Sci. 2012;03:869-78.

Edge effects in the primate community of the biological dynamics of forest fragments project, Amazonas. Brazil. B B Lenz, K M Jack, W R Spironello, Am J Phys Anthropol. 155Lenz BB, Jack KM, Spironello WR. Edge effects in the primate community of the biological dynamics of forest fragments project, Amazonas. Brazil. Am J Phys Anthropol. 2014;155:436-46.

Edge effects influence competition dynamics: a case study of four sympatric arboreal marsupials. K N Youngentob, H J Yoon, N Coggan, D B Lindenmayer, Biol Conserv. 155Youngentob KN, Yoon HJ, Coggan N, Lindenmayer DB. Edge effects influence competition dynamics: a case study of four sympatric arboreal marsupials. Biol Conserv. 2012;155:68-76.

Comparative dynamics of avian communities across edges and interiors of North American ecoregions. K K Karanth, J D Nichols, J R Sauer, J E Hines, J Biogeogr. 33Karanth KK, Nichols JD, Sauer JR, Hines JE. Comparative dynamics of avian communities across edges and interiors of North American ecoregions. J Biogeogr. 2006;33:674-82.

Habitat fragmentation, variable edge effects, and the landscape-divergence hypothesis. W F Laurance, Hem Nascimento, S G Laurance, A Andrade, R M Ewers, K E Harms, PLoS ONE. 21017Laurance WF, Nascimento HEM, Laurance SG, Andrade A, Ewers RM, Harms KE, et al. Habitat fragmentation, variable edge effects, and the landscape-divergence hypothesis. PLoS ONE. 2007;2:e1017.

Vegetation responses to edge environments in old-growth Douglas-fir forests. J Chen, J F Franklin, T A Spies, Ecol Appl. 2Chen J, Franklin JF, Spies TA. Vegetation responses to edge environments in old-growth Douglas-fir forests. Ecol Appl. 1992;2:387-96.

Responses of woody plant seedlings to edge formation in a lowland tropical rainforest. N Sizer, Evj Tanner, Amazonia. Biol Conserv. 91Sizer N, Tanner EVJ. Responses of woody plant seedlings to edge formation in a lowland tropical rainforest. Amazonia. Biol Conserv. 1999;91:135-42.

Edge effects on the behaviour and ecology of Propithecus coquereli in Northwest Madagascar. K C Mcgoogan, University of TorontoMcGoogan KC. Edge effects on the behaviour and ecology of Propithecus coquereli in Northwest Madagascar. University of Toronto; 2011.

Growing-season microclimatic gradients from clearcut edges into old-growth Douglas-fir forests. J Chen, J F Franklin, T A Spies, Ecol Appl. 5Chen J, Franklin JF, Spies TA. Growing-season microclimatic gradients from clearcut edges into old-growth Douglas-fir forests. Ecol Appl. 1995;5:74-86.

Aspect modifies the magnitude of edge effects on bryophyte growth in boreal forests. K Hylander, J Appl Ecol. 42Hylander K. Aspect modifies the magnitude of edge effects on bryophyte growth in boreal forests. J Appl Ecol. 2005;42:518-25.

Effects of edges on plant communities in a managed landscape in northern Wisconsin. E S Euskirchen, J Chen, R Bi, For Ecol Manage. 148Euskirchen ES, Chen J, Bi R. Effects of edges on plant communities in a managed landscape in northern Wisconsin. For Ecol Manage. 2001;148:93-108.

Continuous response functions for quantifying the strength of edge effects. R M Ewers, R K Didham, J Appl Ecol. 43Ewers RM, Didham RK. Continuous response functions for quantifying the strength of edge effects. J Appl Ecol. 2006;43:527-36.

The biodiversity of Madagascar: one of the world's hottest hotspots on its way out. J U Ganzhorn, P P Lowry, G E Schatz, S Sommer, Oryx. 35Ganzhorn JU, Lowry PP, Schatz GE, Sommer S. The biodiversity of Madagascar: one of the world's hottest hotspots on its way out. Oryx. 2001;35:346-8.

Combining global tree cover loss data with historical national forest cover maps to look at six decades of deforestation and forest fragmentation in Madagascar. G Vieilledent, C Grinand, F A Rakotomalala, R Ranaivosoa, J R Rakotoarijaona, T F Allnutt, 10.1016/j.biocon.2018.04.008Biol Conserv. 222Vieilledent G, Grinand C, Rakotomalala FA, Ranaivosoa R, Rakotoarijaona JR, Allnutt TF, et al. Combining global tree cover loss data with historical national forest cover maps to look at six decades of deforestation and forest fragmentation in Madagascar. Biol Conserv. 2018;222:189-97. https ://doi.org/10.1016/j.bioco n.2018.04.008.

Understorey plant species richness and composition in metropolitan forest archipelagos: effects of forest size, adjacent land use and distance to the edge. M Guirado, J Pino, F Rodà, Glob Ecol Biogeogr. 15Guirado M, Pino J, Rodà F. Understorey plant species richness and composition in metropolitan forest archipelagos: effects of forest size, adjacent land use and distance to the edge. Glob Ecol Biogeogr. 2006;15:50-62.

Habitat structure and proximity to forest edge affect the abundance and distribution of forest-dependent birds in tropical coastal forests of southeastern Madagascar. Jem Watson, R J Whittaker, T P Dawson, Biol Conserv. 120Watson JEM, Whittaker RJ, Dawson TP. Habitat structure and proximity to forest edge affect the abundance and distribution of forest-dependent birds in tropical coastal forests of southeastern Madagascar. Biol Conserv. 2004;120:311-27.

Lemur responses to edge effects in the Vohibola III Classified Forest, Madagascar. S M Lehman, A Rajaonson, S Day, Am J Primatol. 68Lehman SM, Rajaonson A, Day S. Lemur responses to edge effects in the Vohibola III Classified Forest, Madagascar. Am J Primatol. 2006;68:293-9.

Edge effects and their influence on lemur density and distribution in southeast Madagascar. S M Lehman, A Rajaonson, S Day, Am J Phys Anthropol. 129Lehman SM, Rajaonson A, Day S. Edge effects and their influence on lemur density and distribution in southeast Madagascar. Am J Phys Anthropol. 2006;129:232-41.

Edge effects on the density of Cheirogaleus major. S M Lehman, A Rajaonson, S Day, Int J Primatol. 27Lehman SM, Rajaonson A, Day S. Edge effects on the density of Cheiro- galeus major. Int J Primatol. 2006;27:1569-88.

Bristol Conservation and Science Foundation, and Conservation International. C Schwitzer, R A Mittermeier, N Davies, S Johnson, J Ratsimbazafy, J Razafindramanana, Bristol, UKIUCN SSC Primate Specialist GroupLemurs of Madagascar: A Strategy for their conservationSchwitzer C, Mittermeier RA, Davies N, Johnson S, Ratsimbazafy J, Razafin- dramanana J, et al. Lemurs of Madagascar: A Strategy for their conserva- tion 2013-2016. Bristol, UK: IUCN SSC Primate Specialist Group, Bristol Conservation and Science Foundation, and Conservation International; 2013.

Averting lemur extinctions amid Madagascar's political crisis. C Schwitzer, R A Mittermeier, S E Johnson, G Donati, M Irwin, H Peacock, Science. 343Schwitzer C, Mittermeier RA, Johnson SE, Donati G, Irwin M, Peacock H, et al. Averting lemur extinctions amid Madagascar's political crisis. Sci- ence. 2014;343:842-3.

Edge effects in tropical dry forests of Madagascar: additivity or synergy? Landsc Ecol. J R Malcolm, K Valenta, S M Lehman, 32Malcolm JR, Valenta K, Lehman SM. Edge effects in tropical dry forests of Madagascar: additivity or synergy? Landsc Ecol. 2016;32:327-41.

Edge effects on tree dendrometrics, abiotics, and mouse lemur densities in western dry forests in Madagascar. S M Lehman, The dwarf and mouse lemurs of Madagascar: Biology, behavior and conservation biogeography of the Cheirogaleidae. Lehman SM, Radespiel U, Zimmermann ELehman SM. Edge effects on tree dendrometrics, abiotics, and mouse lemur densities in western dry forests in Madagascar. In: Lehman SM, Radespiel U, Zimmermann E, editors. The dwarf and mouse lemurs of Madagascar: Biology, behavior and conservation biogeography of the Cheirogaleidae. 2016. p. 462-76.

Edge effects on morphometrics and body mass in two sympatric species of mouse lemurs in Madagascar. R J Burke, S M Lehman, Folia Primatol. 85Burke RJ, Lehman SM. Edge effects on morphometrics and body mass in two sympatric species of mouse lemurs in Madagascar. Folia Primatol. 2014;85:277-91.

Sympatric mouse lemurs in North-West Madagascar: a new rufous mouse lemur species (Microcebus ravelobensis). E Zimmermann, S Cepok, N Rakotoarison, V Zietemann, U Radespiel, Folia Primatol. 69Zimmermann E, Cepok S, Rakotoarison N, Zietemann V, Radespiel U. Sympatric mouse lemurs in North-West Madagascar: a new rufous mouse lemur species (Microcebus ravelobensis). Folia Primatol. 1998;69:106-14.

The ever-increasing diversity in mouse lemurs: three new species in north and northwestern Madagascar. G Olivieri, E Zimmermann, B Randrianambinina, S Rasoloharijaona, D Rakotondravony, K Guschanski, Mol Phylogenet Evol. 43Olivieri G, Zimmermann E, Randrianambinina B, Rasoloharijaona S, Rako- tondravony D, Guschanski K, et al. The ever-increasing diversity in mouse lemurs: three new species in north and northwestern Madagascar. Mol Phylogenet Evol. 2007;43:309-27.

Varying patterns of coexistence of two mouse lemur species (Microcebus ravelobensis and M. murinus) in a heterogeneous landscape. R Rakotondravony, U Radespiel, Am J Primatol. 71Rakotondravony R, Radespiel U. Varying patterns of coexistence of two mouse lemur species (Microcebus ravelobensis and M. murinus) in a het- erogeneous landscape. Am J Primatol. 2009;71:928-38.

Ecological fragmentation effects in mouse lemurs and small mammals in northwestern Madagascar. B Andriatsitohaina, M S Ramsay, F Kiene, S M Lehman, S Rasoloharijaona, R Rakotondravony, Am J Primatol. 8223059Andriatsitohaina B, Ramsay MS, Kiene F, Lehman SM, Rasoloharijaona S, Rakotondravony R, et al. Ecological fragmentation effects in mouse lemurs and small mammals in northwestern Madagascar. Am J Primatol. 2019;82:e23059.

Species-specific usage of sleeping sites in two sympatric mouse lemur species (Microcebus murinus and M. ravelobensis) in northwestern Madagascar. U Radespiel, P Ehresmann, E Zimmermann, Am J Primatol. 59Radespiel U, Ehresmann P, Zimmermann E. Species-specific usage of sleeping sites in two sympatric mouse lemur species (Microcebus murinus and M. ravelobensis) in northwestern Madagascar. Am J Primatol. 2003;59:139-51.

Leaf nest use and construction in the golden-brown mouse lemur (Microcebus ravelobensis) in the Ankarafantsika National Park. S Thorén, F Quietzsch, U Radespie, Am J Primatol. 72Thorén S, Quietzsch F, Radespie U. Leaf nest use and construction in the golden-brown mouse lemur (Microcebus ravelobensis) in the Ankarafant- sika National Park. Am J Primatol. 2010;72:48-55.

Relationship between microhabitat structure and distribution of mouse lemurs (Microcebus spp.) in northwestern Madagascar. A Rendigs, U Radespiel, D Wrogemann, E Zimmermann, Int J Primatol. 24Rendigs A, Radespiel U, Wrogemann D, Zimmermann E. Relationship between microhabitat structure and distribution of mouse lemurs (Micro- cebus spp.) in northwestern Madagascar. Int J Primatol. 2003;24:47-64.

Structural and floristic traits of habitats with differing relative abundance of the lemurs Microcebus murinus and M. ravelobensis in Northwestern Madagascar. L Sehen, D Goetze, C Rajeriarison, E Roger, S Thorén, U Radespiel, Ecotropica. 16Sehen L, Goetze D, Rajeriarison C, Roger E, Thorén S, Radespiel U. Structural and floristic traits of habitats with differing relative abundance of the lemurs Microcebus murinus and M. ravelobensis in Northwestern Madagascar. Ecotropica. 2010;16:15-30.

Tenure security and forest landscape restoration: Results from exploratory research in Boeny, Madagascar. P Ranjatson, M Rebecca, R Randrianasolo, J Mananga, N T Razafimbelo, L Steven, Ranjatson P, Rebecca M, Randrianasolo R, Mananga J, Razafimbelo NT, Steven L. Tenure security and forest landscape restoration: Results from exploratory research in Boeny, Madagascar. In: World bank conference on land and poverty. 2019. p. 1-28.

Model selection and multimodel inference: a practical information-theoretic approach. K P Burnham, D R Anderson, Springer-VerlagBurnham KP, Anderson DR. Model selection and multimodel inference: a practical information-theoretic approach. Springer-Verlag. 2002.

Modelling the edge effect in even-aged Monterey pine (Pinus radiata D. Don) stands. J Cancino, For Ecol Manage. 210Cancino J. Modelling the edge effect in even-aged Monterey pine (Pinus radiata D. Don) stands. For Ecol Manage. 2005;210:159-72.

Edge effects and their influence on habitat suitability calculations: A continuous approach applied to birds of the Atlantic forest. G Zurita, G Pe'er, M I Bellocq, M M Hansbauer, J Appl Ecol. 49Zurita G, Pe'er G, Bellocq MI, Hansbauer MM. Edge effects and their influ- ence on habitat suitability calculations: A continuous approach applied to birds of the Atlantic forest. J Appl Ecol. 2012;49:503-12.

On the efficient identification of an inflection point. D T Christopoulos, Int J Math Sci Comput. 6Christopoulos DT. On the efficient identification of an inflection point. Int J Math Sci Comput. 2016;6:13-20.

An analysis of variance test for normality (complete samples). Ass Shapiro, M B Wilk, Biometrika Trust. 52Shapiro ASS, Wilk MB. An analysis of variance test for normality (complete samples). Biometrika Trust. 1965;52:591-611.

Core area analysis at semi-deciduous forest islands in the Comoé National Park. K J Hennenberg, B Orthmann, I Steinke, S Porembski, Biodivers Conserv. 17NE Ivory CoastHennenberg KJ, Orthmann B, Steinke I, Porembski S. Core area analysis at semi-deciduous forest islands in the Comoé National Park. NE Ivory Coast. Biodivers Conserv. 2008;17:2787-97.

Influence of edge exposure on tree seedling species recruitment in tropical rain forest fragments. J Benítez-Malvido, Biotropica. 35Benítez-Malvido J. Influence of edge exposure on tree seedling species recruitment in tropical rain forest fragments. Biotropica. 2003;35:530-41.

Mammals, edge effects, and the loss of tropical forest diversity. M N Asquith, M Mejia-Chang, Ecology. 86Asquith MN, Mejia-Chang M. Mammals, edge effects, and the loss of tropical forest diversity. Ecology. 2005;86:379-90.

The effects of fragmentation on the susceptibility of a boreal forest ecosystem to wind damage. H Zeng, H Peltola, H Väisänen, S Kellomäki, For Ecol Manage. 257Zeng H, Peltola H, Väisänen H, Kellomäki S. The effects of fragmentation on the susceptibility of a boreal forest ecosystem to wind damage. For Ecol Manage. 2009;257:1165-73.

Comparison of abiotic and structurally defined patch patterns in a hypothetical forest landscape. J Chen, J F Franklin, J S Lowe, Conserv Biol. 10Chen J, Franklin JF, Lowe JS. Comparison of abiotic and structurally defined patch patterns in a hypothetical forest landscape. Conserv Biol. 1996;10:854-62.

Effect of surrounding vegetation on edge-related tree mortality in Amazonian forest fragments. Rcg Mesquita, P Delamo, W F Laurance, Biol Conserv. 91Mesquita RCG, Delamo P, Laurance WF. Effect of surrounding vegeta- tion on edge-related tree mortality in Amazonian forest fragments. Biol Conserv. 1999;91:129-34.

Multiple edge effects and their implications in fragmented landscapes. R J FletcherJr, J Anim Ecol. 74Fletcher RJ Jr. Multiple edge effects and their implications in fragmented landscapes. J Anim Ecol. 2005;74:342-52.

Comparative feeding ecology of sympatric Microcebus berthae and M. murinus. M Dammhahn, P M Kappeler, Int J Primatol. 29Dammhahn M, Kappeler PM. Comparative feeding ecology of sympatric Microcebus berthae and M. murinus. Int J Primatol. 2008;29:1567-89.

Distribution and abundance of the world's smallest primate, Microcebus berthae, in central western Madagascar. L Schäffler, P M Kappeler, Int J Primatol. 35Schäffler L, Kappeler PM. Distribution and abundance of the world's smallest primate, Microcebus berthae, in central western Madagascar. Int J Primatol. 2014;35:557-72.

Forest type influences population densities of nocturnal lemurs in Manompana. Northeastern Madagascar. A Miller, H Mills, T Ralantoharijaona, N V Andriaholinirina, C Misandeau, L Chikhi, Int J Primatol. 39Miller A, Mills H, Ralantoharijaona T, Andriaholinirina NV, Misandeau C, Chikhi L, et al. Forest type influences population densities of noctur- nal lemurs in Manompana. Northeastern Madagascar. Int J Primatol. 2018;39:646-69.

Agent-mediated spatial storage effect in heterogeneous habitat stabilizes competitive mouse lemur coexistence in Menabe Central, Western Madagascar. L Schäffler, J Saborowski, P M Kappeler, BMC Ecol. 15Schäffler L, Saborowski J, Kappeler PM. Agent-mediated spatial stor- age effect in heterogeneous habitat stabilizes competitive mouse lemur coexistence in Menabe Central, Western Madagascar. BMC Ecol. 2015;15:1-15.

Can differences in floristic composition explain variation in the abundance of two sympatric mouse lemur species (Microcebus) in the Ankarafantsika National Park, northwestern Madagascar? Malagasy Nat. L Chanu, D Goetze, C Rajeriarison, E Roger, S Thorén, S Porembski, 6Chanu L, Goetze D, Rajeriarison C, Roger E, Thorén S, Porembski S, et al. Can differences in floristic composition explain variation in the abundance of two sympatric mouse lemur species (Microcebus) in the Ankarafantsika National Park, northwestern Madagascar? Malagasy Nat. 2012;6:83-102.

Comparative feeding ecology of two sympatric mouse lemurs (Microcebus spp.) in northwestern Madagascar. S Thorén, University of Veterinary Medicine HannoverThorén S. Comparative feeding ecology of two sympatric mouse lemurs (Microcebus spp.) in northwestern Madagascar. University of Veterinary Medicine Hannover; 2011.

Seasonal changes in feeding ecology and activity patterns of two sympatric mouse lemur species, the gray mouse lemur (Microcebus murinus) and the golden-brown mouse lemur (M. ravelobensis), in Northwestern Madagascar. S Thorén, F Quietzsch, D Schwochow, L Sehen, C Meusel, K Meares, Int J Primatol. 32Thorén S, Quietzsch F, Schwochow D, Sehen L, Meusel C, Meares K, et al. Seasonal changes in feeding ecology and activity patterns of two sym- patric mouse lemur species, the gray mouse lemur (Microcebus murinus) and the golden-brown mouse lemur (M. ravelobensis), in Northwestern Madagascar. Int J Primatol. 2011;32:566-86.

A predictive model of edge effects. L Ries, T Sisk, Ecology. 85Ries L, Sisk T. A predictive model of edge effects. Ecology. 2004;85:2917-26.

From genetic diversity and structure to conservation: Genetic signature of recent population declines in three mouse lemur species (Microcebus spp. G L Olivieri, V Sousa, L Chikhi, U Radespiel, Biol Conserv. 141Olivieri GL, Sousa V, Chikhi L, Radespiel U. From genetic diversity and structure to conservation: Genetic signature of recent population declines in three mouse lemur species (Microcebus spp.). Biol Conserv. 2008;141:1257-71.

Publisher's Note. Publisher's Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations. Springer Nature remains neutral with regard to jurisdictional claims in pub- lished maps and institutional affiliations.