Using long-term ranging patterns to assess within-group and between-group competition in wild mountain gorillas
BMC Ecology volume 20, Article number: 40 (2020)
Competition within and between social groups determines access to resources and can be inferred from space use parameters that reflect depletion of food resources and competitive abilities of groups. Using location data from 1998 to 2017, we investigated within- and between-group competition in 12 groups of wild mountain gorillas (Gorilla beringei beringei). As within-group feeding competition is expected to increase with group size, an increase in group size is predicted to lead to an increase in the size of annual home ranges and core areas, but to a decrease in fidelity (reuse of an area). Due to asymmetries in competitive abilities, larger groups are expected to have higher exclusivity (degree of non-shared space) of annual home ranges and core areas than smaller groups.
We found evidence of within-group feeding competition based on a positive relationship between group size and both annual home range and core area size as well as a negative relationship between group size and core area fidelity. Additionally, fidelity of core areas was lower than of home ranges. Between-group competition was inferred from a trend for groups with more members and more males to have more exclusive home ranges and core areas. Lastly, annual core areas were largely mutually exclusive.
Our study suggests that non-territorial, group-living animals can have highly dynamic, long-term avoidance-based spacing patterns, both temporally and spatially, to maintain annual core area exclusivity among groups while concurrently shifting these areas annually within overlapping home ranges to avoid resource depletion. Despite ranging in larger home ranges and core areas, larger groups were able to maintain more exclusive ranges than smaller groups, suggesting a competitive advantage for larger groups in between-group competition in a non-territorial species. Together, these findings contribute to understanding how social animals make behavioral adjustments to mitigate the effects of intraspecific competition.
Access to resources, such as food and mates, influences individual fitness . Because resources are limited in space and time, group-living individuals face competition for access to resources with their group members as well as with neighboring groups [2, 3]. Feeding competition can take the form of scramble competition, where resources are exploited by the individual or group that arrives at a resource first, or as contest competition, in which one individual or group has a competitive advantage over another [3,4,5]. Competition for food resources leads to differences in the energetic status among members within a group or between groups . Scramble competition for food within groups results in reduced foraging efficiency for all group members [2, 4] and hence is considered a cost of group living . Between-group scramble competition can lead to energetic costs for all members in both groups [2, 7], whereas between-group contest competition results in only one group suffering energetic costs due to loss of access of resources [2, 7,8,9]. In contrast, mating competition is more intense among males than females in mammals and it can increase a male’s short-term and long-term reproductive success . Competition for food and mates can be reflected in the size and utilization patterns of a group’s home range (area used to survive and reproduce ) and core area (area of intense utilization within a home range, which contains the biologically most relevant resources [11, 12]). Although within-group and between-group competition are expected to act concurrently, very few studies have investigated both at the same time .
According to the ecological constraints model, within-group feeding competition is expected to increase as group size increases, which necessitates groups to increase home range size and daily travel distance to meet the higher energetic requirements of additional group members. Hence, individuals in larger groups have to expend more energy than individuals in smaller groups [14, 15]. Positive relationships between group size and both home range and core area size are the primary markers used to infer scramble competition within groups in the wild [16,17,18], although such scramble competition may be mitigated by dispersal of males and/or females. In particular, female dispersal allows groups to adjust their size based on food availability and feeding competition, with females being able to move to smaller groups and experience less feeding competition [19,20,21]. In female-philopatric species, however, group size cannot decrease easily when feeding competition increases, possibly leading to stronger relationships between groups size and home range size.
The relationship between group size and home range size may be quadratic instead of linear because groups need to balance the trade-off between high within-group competition faced by large groups and high between-group competition faced by small groups, resulting in intermediate-sized groups having an optimal space use strategy in the form of smaller home ranges and/or shorter daily travel distances than small or large groups [22, 23]. Overall, the relationship between group size and home range size is sometimes disputed, especially in folivorous species that rely on abundant food resources as their groups are assumed to not be limited by the availability of food resources (e.g. [7, 18, 24]), and the non-linear relation between home range size and group size has rarely been tested.
Intraspecific competition for food and mates can also be measured with other behavioral indicators. Specifically, fidelity and exclusivity of home ranges and core areas reflect depletion of food resources and competitive abilities of groups but these variables have only seldom been examined [25,26,27,28]. Home range fidelity, which is the tendency of animals to return to and reuse previously used areas , could improve fitness of groups via increased foraging efficiency by allowing them to evaluate the quality, distribution and predictability of habitats over space and time [29,30,31]. However, it may be more efficient energetically to shift a range because depletion of food resources within a home range may lower foraging efficiency , especially as group size increases [33, 34]. Although also dependent on regeneration rates of food plants, variation in home range fidelity across time could be a valuable behavioral indicator of within-group scramble competition .
For non-territorial species, the level of between-group competition can be reflected in the degree of shared space among neighbors. In contrast to territorial species, non-territorial species that have overlapping home ranges are not expected to benefit from exclusive access to food resources resulting from ownership advantages or from intergroup dominance in between-group contest competition [8, 35]. Intergroup dominance, which usually arises from asymmetries in competitive abilities due to an advantage of having more members or more males [36,37,38], can increase the quality of the dominant group’s home range or the area of exclusive access [8, 39, 40]. This may be due to direct food defense by males to attract females (resource defense polygyny [39, 41, 42]) or as a consequence of direct mate defense (“hired guns” [42, 43]). Although non-territorial species are not expected to benefit from ownership advantages, they may still benefit from better access to resources in part of their home range [26, 44, 45]. However, for non-territorial species, it remains largely unclear whether groups have a competitive advantage in intergroup competition and how this is expressed.
To infer within- and between-group competition for food and mates in a non-territorial social mammal, we investigated the impact of group size on three parameters of space use in wild Bwindi mountain gorillas: size, fidelity and exclusivity of annual home ranges and core areas. The 2 populations of mountain gorillas are found in Bwindi Impenetrable National Park, Uganda, and in the Virunga Volcanoes of Uganda, Rwanda, and the Democratic Republic of Congo. Mountain gorillas live in cohesive social units consisting of one or more adult males, several adult females and their offspring (mean group size = 11) [46, 47]. Gorilla females disperse, which gives them the possibility to respond to an increase in within-group feeding competition by moving into smaller groups where there should be less competition . The 2 populations feed predominantly on herbaceous vegetation, which is available throughout their range and throughout the year [48, 49]. Overall food biomass is higher in the Virungas than in Bwindi, which seems to be reflected in differential competitive regimes in the 2 populations . Little evidence of costs of within-group feeding competition is apparent in the Virunga gorilla groups [51, 52] and group size has no clear effect on female reproductive success . In contrast, Bwindi gorillas appear to experience some within-group feeding competition as groups reduce the frequency of revisits to particular areas as group size increases, although neither monthly home range size nor daily travel distance is affected by group size . Within Bwindi, herbaceous food plants, the gorillas’ main diet, show spatial variation in availability despite little seasonal variability . Assuming similar biomass regeneration rates for Bwindi as in the Virunga Volcanoes, food renewal rates, i.e. the time needed for the regeneration of biomass, are between 240 and 270 days . Therefore, long-term studies, assessing spacing patterns on an annual scale, may best reflect the impact of depletion and regeneration of food resources (see ).
Competition between gorilla groups is mainly due to mating competition , with males being the main participants and females transferring to neighboring groups during intergroup encounters [59, 60]. Despite a high degree of home range overlap among neighboring groups in both populations [26, 52], intergroup encounters occur infrequently  and a recent 1 year study in Bwindi showed that groups have largely mutually exclusive core areas . These core areas, in which the gorillas spend 50% of their time, contain higher herbaceous food availability than the rest of the home ranges. Such a spacing pattern appears to result from active avoidance among neighboring groups and seems to be stimulated by strong between-group competition for mates .
Using long-term ranging patterns of 12 groups over the course of 19 years (median = 6, range 1–19 years), we investigate how Bwindi gorillas concurrently adjust their behavior to compensate for the costs of both within- and between-group competition in the long term, complementing a 1 year study that investigated these competitive patterns in the short term [26, 54]. We tested the following predictions (Table 1):
In accordance with the ecological constraints model [14, 15], we predicted a positive relationship between group size and both annual home range and core area size. Alternatively, we predicted a non-linear relationship (U-shaped) between both annual home range and core area size and group size because intermediate-sized groups may have an optimal space use strategy . As there is greater resource depletion as group size increases [33, 34], we predicted a negative relationship between group size and both annual home range and core area fidelity between consecutive years. We decided to investigate fidelity of home ranges and core areas between consecutive years because there is low temporal variation in food availability and food renewal rates are expected to be between 240 and 270 days (assuming similar food renewal rates for Bwindi as in the Virunga Volcanoes ). Because highly used core areas of Bwindi gorillas were shown to be largely mutually exclusive and to contain higher herbaceous food availability than the rest of the home ranges in a 1 year study , we predicted that annual core area fidelity will be higher than annual home range fidelity across time.
Due to asymmetries in competitive abilities [8, 39] and because males are the active participants in intergroup encounters [59, 60], we expected that larger groups and groups with more males will maintain more exclusive (i.e. less overlap among neighbors) annual home ranges and core areas than smaller groups and groups with fewer males. Lastly, in a 1 year study Bwindi gorilla groups appeared to actively avoid competition with neighbors, resulting in shared home ranges but exclusive core areas , so we predicted that annual core areas will be more exclusive than annual home ranges across time.
Home range and core area size
Bwindi mountain gorillas had a mean annual home range size of 10 km2 ± SE 0.44 (range 4.1–22.9 km2; see Additional file 1). When testing for the effect of group size and a quadratic value of group size (squared group size) on annual home range size, we found the full-null model comparison to show a trend (likelihood ratio test: χ2 = 4.945, df = 2, p = 0.084). The effect of group size was a trend (Est ± SE = 0.158 ± 0.069, p = 0.092) and the quadratic term was not significant (Est ± SE = − 0.090 ± 0.042, p = 0.108). The statistical model indicated that increasing the group size from 12 to 13 weaned individuals would increase the annual home range size by 2.5% (dashed line in Fig. 1a). According to the R2 value from a post-hoc linear regression, group size (and its quadratic term) accounted for 36.3% of the variance in the annual home range size.
The Bwindi mountain gorillas had a mean annual core area size (50% contours of the fixed kernel density estimation) of 3 km2 ± 0.15 (range 1.1–7.2 km2; see Additional file 1). When testing for the effect of group size and a quadratic value of group size on annual core area size, we found a significant difference between the full model and the null model (χ2 = 6.444, df = 2, p = 0.040). The effect of group size was a trend (Est ± SE = 0.084 ± 0.042, p = 0.065) and so was the quadratic term (Est ± SE = − 0.051 ± 0.026, p = 0.067). Visual inspection indicated that both the annual home range size and the mean annual core area size were predicted by the models to increase monotonically throughout most of the observed range of group sizes, before leveling off and then declining at the largest group sizes (Fig. 1). Results weakened in 2 additional models that included location as a category variable for the spatial variability in food availability (Additional file 2).
Home range and core area fidelity
Examining the effect of group size on annual home range and core area fidelity (Bhattacharyya affinity = BA), we found that the mean BA for the annual home ranges was 0.6 ± 0.02 (range 0.31–0.79) and 0.2 ± 0.01 (range 0.02–0.39) for the annual core areas. There was no significant effect of group size on the BA for the annual home ranges (χ2 = 2.397, df = 2, p = 0.302), but it had a significant effect on the annual core area BA (χ2 = 6.137, df = 2, p = 0.047). Variation in group size among groups was as high as the variation within groups across time. We found a negative effect of the between-groups effect on the BA for the annual core area (Est ± SE = − 0.036 ± 0.013, χ2 = 5.873, p = 0.015; Fig. 2), whereas the within-groups effect had no apparent impact (Est ± SE = − 0.005 ± 0.010, χ2 = 0.253, p = 0.615). Against our prediction, we found that annual home range fidelity was significantly higher than annual core area fidelity (Est ± SE = − 0.372 ± 0.019, χ2 = 35.288, p < 0.001; Fig. 3).
Home range and core area exclusivity
Bwindi gorillas used on average 54% ± 4.44 (range 3.1–98.9%) of their home ranges exclusively and on average 83% ± 3.58 (range 20.7–100%) of their core areas exclusively. The size of the exclusively used part of the annual home range tended to increase as group size increased (Est ± SE = 4.210 ± 1.841, χ2 = 3.342, df = 1, p = 0.068). For example, the statistical model indicated that increasing the group size from 12 to 13 weaned individuals would increase the size of the exclusively used part of the annual home range by 11.9%. According to the R2 value from a post-hoc univariate linear regression, group size accounted for 56.8% of the variance in the size of the exclusively used part of the annual home range. The number of males in a group did not have a significant effect on the size of the exclusively used part of its annual home range (Est ± SE = 1.445 ± 1.323; χ2 = 0.784, df = 1, p = 0.376). For the annual core area, we found that the size of the exclusively used part tended to increase as group size increased (Est ± SE = 1.823 ± 0.937, χ2 = 3.034, df = 1, p = 0.082) and as the number of males increased (Est ± SE = 1.928 ± 0.800, χ2 = 3.199, df = 1, p = 0.074). Lastly, annual core areas were more exclusive than annual home ranges (Est ± SE = 28.083 ± 5.120, χ2 = 14.963, p < 0.001; Fig. 4).
We found evidence of within- and between-group competition in Bwindi mountain gorillas using a long-term data set and variables rarely examined (fidelity and exclusivity of annual home ranges and core areas). The positive relationship between group size and both annual home range and core area size and the negative relationship between group size and core area fidelity suggest that increasing group size results in increased within-group feeding competition. We found some evidence that Bwindi gorillas also experienced between-group competition, as indicated by a trend for groups with more members and/or more males to have more exclusive home ranges and core areas. This suggests a competitive advantage for groups with more members and males, a relationship typically found in territorial species ([8, 9, 61] but see [44, 45]). Additionally, annual core areas were largely mutually exclusive across time (i.e. had considerably less overlap than annual home ranges), suggesting that groups avoid competition with neighboring groups in the long term, which might ultimately be due to mating competition.
Within-group feeding competition
In accordance with the ecological constraints model [14, 15], both home range and core area size increased as group size increased on an annual scale. We found that small gorilla groups suffered costs from between-group competition by having less exclusive home ranges and core areas as well as that larger groups experienced costs from within-group competition. However, the combined linear effects of both within- and between-group competition did not result in the expected U-shaped relationship between group size and annual home range or core area size [19, 20]. Instead, we found evidence of a concave relationship, with the annual home range and the core area size leveling off and then declining at the larger group sizes. Similar patterns have been reported for the travel times and distances of the Virunga gorillas, which were attributed to the costs of within-group feeding competition for the smaller groups and the benefits of between-group feeding competition for the largest group . Concave patterns were not observed, however, in the home ranges of those Virunga gorilla groups .
In Bwindi gorillas, neither daily travel distance nor monthly home range size increased significantly as group size increased  (see Table 2). Variation in the spatial availability of herbaceous vegetation might be small, such that on the scales of daily and monthly movement patterns groups may have found alternative strategies to mitigate the effect of within-group scramble competition, such as increasing group spread [16, 17, 63]. However, on the larger scale of annual movement, larger groups adjusted to greater within-group feeding competition and thus experienced higher energetic costs compared to smaller groups. Including location of groups as a proxy for food availability did not result in an increase in effect size, suggesting that it was not an accurate reflection of food availability on the scale that the gorillas experience. Bwindi gorillas only use a small proportion of their annual home range on a monthly basis (10 versus 1 km2; Additional file 1, ) and not increasing these monthly range sizes as group size increases should lead to higher depletion of food resources within these areas. To avoid these depleted areas, gorillas may shift their monthly home ranges within the annual home range, with larger groups having presumably less overlap among monthly home ranges than smaller groups, leading to an increase of the annual home ranges and core areas as group size increases.
Group size also affected inter-annual core area fidelity, providing further evidence of within-group scramble competition. As group size increased core area fidelity decreased, presumably due to depletion of resources in relation to regeneration rates , but there was no effect on home range fidelity. Additionally, core area fidelity was lower than home range fidelity, again presumably due to depletion of food resources in the much smaller but highly utilized core areas. High home range fidelity but low core area fidelity suggests that groups shifted their annual core areas within stable annual home ranges. If groups deplete the resources available in their home range, they may need to move elsewhere. But if they are surrounded by other groups that also deplete their own home ranges, there is no guarantee that they should find more food if they move. In addition, they would pay the cost of having to search for their food and memorize the food distribution in new areas. Thus, there are reasons why large groups, which deplete their resources faster than small groups, may not move their annual home range more often.
Assuming similar biomass regeneration rates for Bwindi as in the Virunga Volcanoes (240–270 days ), these 2 rotational systems (shifting both monthly ranges and annual core areas within annual home ranges) may allow for regeneration of food resources in those areas. This also might ensure high food availability in the largely mutually exclusive core areas in the long term, which would be in concordance with a previous study that found higher food availability of core areas compared to home ranges in the short term . Furthermore, such a forage rotation system follows predictions of the optimal foraging theory, which predicts animals should forage in areas that offer the highest average rate of energy intake [64, 65].
Between-group mating competition
For non-territorial species that have overlapping home ranges, the level of between-group competition can be reflected in the degree of shared space among neighbors. Our results lend support to our predictions for between-group competition, as indicated by a trend for greater exclusivity of home ranges and core areas as group size and the number of males increases. Having more exclusive access to a home range and core area should lead to more exclusive and higher access to food resources (see [38, 62]). This should result in a difference in energetic status among groups and suggests that larger groups have a competitive advantage in intergroup competition, a pattern typically found in territorial species ([36, 66] but see [44, 45]). This reduced cost of between-group feeding competition for larger groups  could buffer the costs of increased within-group feeding competition [5, 67].
Bwindi mountain gorilla groups maintain core area exclusivity in the long term as annual core areas were more mutually exclusive than the respective home ranges. Groups used on average 83% of their core areas exclusively, whereas only 54% of their home ranges were exclusive. Together, this provides further support to the results of a 1 year study, which found exclusive use of core areas by non-territorial Bwindi mountain gorilla groups, suggesting that neighbors actively avoid each other .
Three proximate mechanisms that may allow groups to maintain this exclusivity by avoiding each other in the short term, may also work in the long term: remembering the locations of intergroup encounters to avoid these areas, using chest beats as a long-distance signal to locate neighbors and using signs of foraging to avoid previously used areas  (see also ). One ultimate mechanism for the observed avoidance-based spacing pattern might be male mate defense. In mountain gorillas, males have been considered as “hired guns”, in which males defend their mates and offspring by keeping extragroup males away, thereby indirectly defending food resources within these areas for the male’s group members [26, 43, 68]. This is in contrast to males exhibiting resource defense polygyny as found in several primate species, in which males indirectly defend females via resource defense [39, 42, 69], although both strategies can lead to high-quality home ranges.
We found a trend for groups with more males to have more exclusive core areas than groups with fewer males, which suggests a competitive advantage of additional males for directly defending mates and offspring within a group as well as indirectly defending food resources (see also ). Bwindi females residing in multimale groups might ultimately increase their reproductive success as an increasing number of males can provide increasingly exclusive core areas and thus exclusive access to food resources [2, 45, 67] (but see ).
Through an investigation of long-term ranging patterns that reflect competition within and between groups, we found that despite having overlapping home ranges, mountain gorilla groups have largely mutually exclusive core areas that they shift on an annual scale to minimize the effect of within-group feeding competition. Furthermore, our results show that as group size increases, home range and core area size increases as well. Increasingly large home ranges should become more difficult to defend [71, 72] and overlap with neighbors is expected to increase . Consequently, increasing home range size should lead to an increasing loss of food resources to neighbors . Unexpectedly, we found that larger groups that range in larger home ranges and core areas had less overlap with their neighbors and hence were able to monopolize food resources in the exclusively used parts, thereby reducing the amount of food lost to neighbors. Together, this suggests that gorillas have evolved a highly dynamic avoidance-based spacing pattern, both temporally and spatially.
Through an investigation of the effect of group size on several space use parameters in a non-territorial species, our results show that while an increase in group size is beneficial for between-group competition, it is costly in terms of within-group feeding competition. Although many studies investigate within- or between-group competition independently (e.g. [9, 17]), we found evidence that both competitive regimes act simultaneously and animals adapt their behavioral response to balance the costs and benefits of between- and within-group competition concurrently  (see Table 2). These results stress the value of investigating both competitive regimes simultaneously across several scales to gain a comprehensive understanding of the various behavioral adjustments social animals exhibit to compensate for the costs of intraspecific competition.
Study site and data collection
We used location and demographic data on 12 groups of habituated mountain gorillas in Bwindi Impenetrable National Park, Uganda, from 1998 to 2017 (Fig. 5). Most of the location data from the groups habituated for tourism were taken from the long-term records of the Uganda Wildlife Authority. Permission for use of the data was obtained from the division of research and monitoring department of Uganda Wildlife Authority and access to the data can be obtained from them. All groups were monitored daily by park staff and in part by researchers to ascertain group membership. Mean annual group size was defined as the average number of weaned individuals per group and year (mean = 13, range 4.0–18.2 individuals), whereas the mean number of males (> 12 years of age; e.g. ) was defined as the average number of adult males per group and year (mean = 3, range 1.0–5.8 males). Using handheld global positioning system (GPS) units (GARMIN), one to several location points were recorded for each observation day and group. The number of observation days per year differed among groups, so location data coverage varied among groups and years (Fig. 6). Therefore, we randomly selected one location point per day (see also ).
Home range and core area size
Using the fixed kernel density estimation  and one location point per day, we determined annual home ranges (90% contours) and core areas (50% contours) per group applying the adehabitatHR package  in R 3.4.3 . The kernel density estimation produces utilization distributions, which are probability distributions that describe the groups’ use of space . As the number of observation days varied per group and year, we determined the minimum number of observation days needed for the annual home range size to reach an asymptote per group and year. This differed largely among years and groups, so we used only years with at least 125 observation days per group and randomly selected 125 observation days per group and year (n = 70) to estimate annual home range and core area sizes. We fixed the bandwidth to h = 200 because it creates annual home range estimates with relatively little fragmentation [26, 52], which is expected for Bwindi gorillas moving on average 975 m a day .
Home range and core area fidelity
The Bhattacharyya affinity (BA)  compares the overall similarity between utilization distributions [80, 81]. Using the annual kernel home ranges and core areas and the package adehabitatHR  in R 3.4.3 , we determined the BA between 2 consecutive years for each group (n = 49). Values range between zero and one, with one indicating identical utilization distributions and zero indicating no similarity between utilization distributions.
Home range and core area exclusivity
For each group and year, we determined both the size and the percentage of the annual home range and core area that was not shared with any neighboring group. To calculate the exclusively used part of the annual kernel home range and core area per group, we subtracted the shared parts from the annual kernel home ranges and core areas (n = 32; for more details see ). We included only years for which we had data on all habituated groups that ranged in the same area. We could not include the unhabituated gorilla groups in our analysis because the locations of their home ranges can only be approximated using few data points from the Bwindi censuses (e.g. [47, 82]). However, previous analyses have shown that it is unlikely that the habituated and unhabituated groups share notable amounts of the same ranges and hence we assume that excluding the unhabituated groups did not have a large effect on our results . For data processing and analyses, we used the packages spatstat , splancs  and SDMTools  in R 3.4.3 .
Models and statistical analyses
We used linear mixed models (LMMs)  to quantify the effects of group size and group size squared (quadratic) on both annual home range and core area size (square-root-transformed; for both n = 70 observations from 12 groups and as much as 19 years). We included the general location the study groups ranged in as a control predictor as a factor with four levels (Buhoma, Ruhija, Nkuringo and Rushaga; Fig. 5) in 2 additional models as a proxy to control for the spatial variation in food availability among the different locations . By doing so, we accounted for the possibility that larger groups range in areas with higher food availability than smaller groups . We also investigated if the inclusion of the Nkuringo Group, which was observed primarily in uncultivated land outside the park and was observed to feed on palatable crops on 25% of observation days , had a notable effect on the results in comparison to other groups. It did not, so we retained it in all analyses (see Additional file 3).
We fitted 2 LMMs  to examine the effect of group size on annual home range and core area fidelity (BA estimates; for both n = 49 observations from 9 groups and as much as 18 years). For the test predictor group size, we determined the mean of the annual group sizes for the 2 corresponding consecutive years because variation between years was very small. As the variation in group size among groups was as high as the variation within groups across time, the mean of group size per group (= between-groups variation) and group size centered to a mean of zero per group (= within-groups variation) were included as test predictors in the model . To test whether core area fidelity was higher than home range fidelity, we used a LMM  (n = 98 observations from 9 groups, as much as 18 years and 49 combinations of group and year). The response variable was the BA estimates  of both the annual home ranges and core areas, a measure that compares the overall similarity between utilization distributions. As a test predictor, we included whether the estimate originated from the core area or not as a factor with 2 levels (yes = core area and no = home range).
To investigate whether both group size and the number of males influenced the size of the exclusively used part of both the home range and the core area (square-root-transformed; for all n = 32 observations from 10 groups and 9 years), we used a LMM  per response variable with either group size or number of males as the test predictor. To test whether core area exclusivity was higher than home range exclusivity, we used a LMM  (n = 64 observations from 9 groups, as much as 18 years and 32 combinations of group and year). The response variable was the percent of the exclusively used annual home range and core area and as test predictor we included whether the estimate was from the core area or not as a factor with 2 levels (yes = core area and no = home range; see also Additional file 4).
We fitted all models in R 3.4.3  with Gaussian error structure and identity link  and implemented them using the functions lmer of the lme4 package . All predictor variables were z-transformed to a mean of zero and a standard deviation of one . To control for repeated observations, we included group ID and year ID as random effects in all models. Additionally, for the models testing whether core area fidelity and exclusivity were higher than home range fidelity and exclusivity, we also included the combination of group and year as a random effect. To keep error I rate at the nominal level of 5%, we included random slopes in all models where applicable (Additional file 5) [92, 93]. To check for the assumption of normally distributed and homogenous residuals, we visually inspected qqplots and the residuals plotted against fitted values. We found no violations. We investigated model stability by excluding each level of the random effects one at a time, including one group that foraged on crops outside the park. Comparing estimates for each predictor with those obtained for the full data set suggested no influential levels of random effects (Additional file 3). To investigate collinearity issues of the models having more than one test predictor, we determined Variance Inflation Factors (VIF)  using the function vif of the package car . We applied the function to a model lacking the random effects and found no issues (maximum VIF across all models = 1.16). To establish significance of the full model compared to the null model lacking the test predictors (p < 0.05), we used a likelihood ratio test . For the models having more than one test predictor, we determined individual p values based on likelihood ratio tests, which compared the full model with respective reduced models excluding the test predictors one at a time .
Availability of data and materials
The datasets used and analyzed during the current study are available from the corresponding author on reasonable request.
Global positioning system
Linear mixed model
Morales JM, Moorcroft PR, Matthiopoulos J, Frair JL, Kie JG, Powell RA, et al. Building the bridge between animal movement and population dynamics. Philos Trans R Soc B Biol Sci. 2010;365:2289–301.
Janson CH, van Schaik CP. Recognizing the many faces of primate food competition: methods. Behaviour. 1988;105:165–86.
van Schaik CP. The ecology of social relationships amongst female primates. In: Folley RA, Standen V, editors. Comparative socioecology. The behavioural ecology of humans and other mammals. Oxford: Blackwell Scientific Publications; 1989. p. 195–218.
Nicholson AJ. An outline of the dynamics of animal populations. Aust J Zool. 1954;2:9–65.
Janson CH. Intra-specific food competition and primate social structure: a synthesis. Behaviour. 1988;105:1–17.
Crofoot MC. Mating and feeding competition in white-faced capuchins (Cebus capucinus): the importance of short-and long-term strategies. Behaviour. 2007;144:1473–95.
Snaith TV, Chapman CA. Primate group size and interpreting socioecological models: do folivores really play by different rules? Evol Anthropol. 2007;16:94–106.
Mosser A, Packer C. Group territoriality and the benefits of sociality in the African lion, Panthera leo. Anim Behav. 2009;78:359–70.
Roth AM, Cords M. Effects of group size and contest location on the outcome and intensity of intergroup contests in wild blue monkeys. Anim Behav. 2016;113:49–58.
Burt WH. Territoriality and home range concepts as applied to mammals. J Mammal. 1946;24:346–52.
Marzluff JM, Knick ST, Millspaugh JJ. High-tech behavioral ecology: modeling the distribution of animal activities to better understand wildlife space use and resource selection. In: Millspaugh JJ, Marzluff JM, editors. Radiotracking and animal populations. San Diego: Academic Press; 2001. p. 309–326.
Asensio N, Lusseau D, Schaffner CM, Aureli F. Spider monkeys use high-quality core areas in a tropical dry forest. J Zool. 2012;287:250–8.
Teichroeb JA, Sicotte P. Cascading competition: the seasonal strength of scramble influences between-group contest in a folivorous primate. Behav Ecol Sociobiol. 2018;72:6.
Janson CH, Goldsmith ML. Predicting group size in primates: foraging costs and predation risks. Behav Ecol. 1995;6:326–36.
Chapman CA, Chapman LJ. Constraints on group size in red colobus and red-tailed guenons: examining the generality of the ecological constraints model. Int J Primatol. 2000;21:565–85.
Chapman CA, Chapman LJ. Determinants of group size in social primates: the importance of travel costs. In: Boinski S, Garber PA, editors. On the move: how and why animals travel in groups. Chicago: The University of Chicago Press; 2000. p. 24–42.
Snaith TV, Chapman CA. Red colobus monkeys display alternative behavioral responses to the costs of scramble competition. Behav Ecol. 2008;19:1289–96.
Teichroeb JA, Sicotte P. Test of the ecological-constraints model on ursine colobus monkeys (Colobus vellerosus) in Ghana. Am J Primatol. 2009;71:49–59.
Stokes EJ, Parnell RJ, Olejniczak C. Female dispersal and reproductive success in wild western lowland gorillas (Gorilla gorilla gorilla). Behav Ecol Sociobiol. 2003;54:329–39.
Manguette ML, Robbins AM, Breuer T, Stokes EJ, Parnell RJ, Robbins MM. Female dispersal patterns influenced by male tenure duration and group size in western lowland gorillas. Behav Ecol Sociobiol. 2020;74:81. https://doi.org/10.1007/s00265-020-02863-8.
Robbins AM, Robbins MM. Dispersal patterns of females in the genus Gorilla. In: Furuichi T, Yamagiwa J, Aureli F, editors. Dispersing primate females. Berlin: Springer; 2015. p. 75–104.
Markham AC, Gesquiere LR, Alberts SC, Altmann J. Optimal group size in a highly social mammal. Proc Natl Acad Sci. 2015;112:14882–7.
Stevenson PR, Castellanos MC. Feeding rates and daily path range of the Colombian woolly monkeys as evidence for between-and within-group competition. Folia Primatol. 2000;71:399–408.
Snaith TV, Chapman CA. Towards an ecological solution to the folivore paradox: patch depletion as an indicator of within-group scramble competition in red colobus monkeys (Piliocolobus tephrosceles). Behav Ecol Sociobiol. 2005;59:185–90.
Markham AC, Guttal V, Alberts SC, Altmann J. When good neighbors don’t need fences: temporal landscape partitioning among baboon social groups. Behav Ecol Sociobiol. 2013;67:875–84.
Seiler N, Boesch C, Mundry R, Stephens C, Robbins MM. Space partitioning in wild, non-territorial mountain gorillas: the impact of food and neighbours. Roy Soc Open Sci. 2017;4:170720.
van Beest FM, Mysterud A, Loe LE, Milner JM. Forage quantity, quality and depletion as scale-dependent mechanisms driving habitat selection of a large browsing herbivore. J Anim Ecol. 2010;79:910–22.
Weatherhead PJ, Boak KA. Site infidelity in song sparrows. Anim Behav. 1986;34:1299–310.
Switzer PV. Site fidelity in predictable and unpredictable habitats. Evol Ecol. 1993;7:533–55.
Hinde RA. The biological significance of territories in birds. Ibis. 1956;58:24–44.
Williams JM, Pusey AE, Carlis JV, Farm BP, Goodall J. Female competition and male territorial behaviour influence female chimpanzees’ ranging patterns. Anim Behav. 2002;63:347–60.
Charnov EL. Optimal foraging, the marginal value theorem. Theor Popul Biol. 1976;9:129–36.
Cody M. Finch flocks in the Mohave Desert. Theor Popul Biol. 1971;2:142–58.
Davies NB, Houston AI. Owners and satellites: the economics of territory defence in the pied wagtail, Motacilla alba. J Anim Ecol. 1981;50:157–80.
Brown JL, Orians GH. Spacing patterns in mobile animals. Annu Rev Ecol Syst. 1970;1:239–62.
Cassidy KA, MacNulty DR, Stahler DR, Smith DW, Mech LD. Group composition effects on aggressive interpack interactions of gray wolves in Yellowstone National Park. Behav Ecol. 2015;26:1352–60.
Maynard Smith J, Parker GA. The logic of asymmetric contests. Anim Behav. 1976;24:159–75.
Scarry CJ. Between-group contest competition among tufted capuchin monkeys, Sapajus nigritus, and the role of male resource defence. Anim Behav. 2013;85:931–9.
Harris TR. Between-group contest competition for food in a highly folivorous population of black and white colobus monkeys (Colobus guereza). Behav Ecol Sociobiol. 2006;61:317–29.
Crofoot MC. The cost of defeat: Capuchin groups travel further, faster and later after losing conflicts with neighbors. Am J Phys Anthropol. 2013;152:79–85.
Emlen S, Oring L. Ecology, sexual selection, and the evolution of mating systems. Science. 1977;197:215–23.
Fashing PJ. Male and female strategies during intergroup encounters in guerezas (Colobus guereza): evidence for resource defense mediated through males and a comparison with other primates. Behav Ecol Sociobiol. 2001;50:219–30.
Rubenstein DI. Ecology and sociality in horses and zebras. In: Rubenstein DI, Wrangham RW, editors. Ecology and social evolution: birds and mammals. Princeton: Princeton University Press; 1986. p. 282–302.
Markham AC, Alberts SC, Altmann J. Intergroup conflict: ecological predictors of winning and consequences of defeat in a wild primate population. Anim Behav. 2012;84:399–403.
Richter C, Heesen M, Nenadić O, Ostner J, Schülke O. Males matter: Increased home range size is associated with the number of resident males after controlling for ecological factors in wild Assamese macaques. Am J Phys Anthropol. 2016;159:52–62.
Gray M, Roy J, Vigilant L, Fawcett K, Basabose A, Cranfield M, et al. Genetic census reveals increased but uneven growth of a critically endangered mountain gorilla population. Biol Cons. 2013;158:230–8.
Roy J, Vigilant L, Gray M, Wright E, Kato R, Kabano P, et al. Challenges in the use of genetic mark-recapture to estimate the population size of Bwindi mountain gorillas (Gorilla beringei beringei). Biol Cons. 2014;180:249–61.
Watts DP. Composition and variability of mountain gorilla diets in the central Virungas. Am J Primatol. 1984;7:323–56.
Ganas J, Robbins MM, Nkurunungi JB, Kaplin BA, McNeilage A. Dietary variability of mountain gorillas in Bwindi impenetrable National Park, Uganda. Int J Primatol. 2004;25:1043–72.
Wright E, Grueter CC, Seiler N, Abavandimwe D, Stoinski TS, Ortmann S, et al. Energetic responses to variation in food availability in the two mountain gorilla populations (Gorilla beringei beringei). Am J Phys Anthropol. 2015;158:487–500.
Watts DP. Long-term habitat use by mountain gorillas (Gorilla gorilla beringei). 1. Consistency, variation, and home range size and stability. Int J Primatol. 1998;19:651–80.
Caillaud D, Ndagijimana F, Giarrusso AJ, Vecellio V, Stoinski TS. Mountain gorilla ranging patterns: Influence of group size and group dynamics. Am J Primatol. 2014;76:730–46.
Robbins MM, Robbins AM, Gerald-Steklis N, Steklis HD. Socioecological influences on the reproductive success of female mountain gorillas (Gorilla beringei beringei). Behav Ecol Sociobiol. 2007;61:919–31.
Seiler N, Boesch C, Stephens C, Ortmann S, Mundry R, Robbins MM. Social and ecological correlates of space use patterns in Bwindi mountain gorillas. Am J Primatol. 2018;80:e22754.
Ganas J, Nkurunungi JB, Robbins MM. A preliminary study of the temporal and spatial biomass patterns of herbaceous vegetation consumed by mountain gorillas in an afromontane rain forest. Biotropica. 2009;41:37–46.
Plumptre AJ. The effects of trampling damage by herbivores on the vegetation of the Park National des Volcans, Rwanda. Afr J Ecol. 1993;32:115–29.
Börger L, Franconi N, Ferretti F, Meschi F, De Michele G, Gantz A, et al. An integrated approach to identify spatiotemporal and individual-level determinants of animal home range size. Am Nat. 2006;168:471–85.
Robbins MM, Sawyer SC. Intergroup encounters in mountain gorillas of Bwindi Impenetrable National Park, Uganda. Behaviour. 2007;144:1497–520.
Sicotte P. Inter-group encounters and female transfer in mountain gorillas: influence of group composition on male behavior. Am J Primatol. 1993;30:21–36.
Mirville MO, Ridley AR, Samedi JPM, Vecellio V, Ndagijimana F, Stoinski TS, et al. Factors influencing individual participation during intergroup interactions in mountain gorillas. Anim Behav. 2018;144:75–86.
McComb K, Packer C, Pusey AE. Roaring and numerical assessment in contests between groups of female lions, Panthera leo. Anim Behav. 1994;47:379–87.
Grueter CC, Robbins AM, Abavandimwe D, Vecellio V, Ndagijimana F, Stoinski TS, et al. Quadratic relationships between group size and foraging efficiency in a herbivorous primate. Sci Rep. 2018. https://doi.org/10.1038/s41598-018-35255-0.
Watts DP. Strategies of habitat use by mountain gorillas. Folia Primatol. 1991;56:1–16.
Pyke GH, Pulliam HR, Charnov EL. Optimal foraging—selective review of theory and tests. Q Rev Biol. 1977;52:137–54.
Stephens DW, Krebs JR. Foraging theory. Princeton: Princeton University Press; 1986.
Stamps JA, Krishnan VV. How territorial animals compete for divisible space: a learning-based model with unequal competitors. Am Nat. 2001;157:154–69.
Koenig A. Competition for resources and its behavioral consequences among female primates. Int J Primatol. 2002;23:759–83.
Watts DP. Infanticide in mountain gorillas: new cases and a reconsideration of the evidence. Ethology. 1989;81:1–18.
Williams JM, Oehlert GW, Carlis JV, Pusey AE. Why do male chimpanzees defend a group range? Anim Behav. 2004;68:523–32.
Robbins AM, Gray M, Basabose A, Uwingeli P, Mburanumwe I, Kagoda E, et al. Impact of male infanticide on the social structure of mountain gorillas. PLoS ONE. 2013;8:e78256.
Mitani JC, Rodman PS. Territoriality: the relation of ranging pattern and home range size to defendability, with an analysis of territoriality among primate species. Behav Ecol Sociobiol. 1979;5:241–51.
Lowen C, Dunbar RIM. Territory size and defendability in primates. Behav Ecol Sociobiol. 1994;35:347–54.
Jetz W, Carbone C, Fulford J, Brown JH. The scaling of animal space use. Science. 2004;306:266–8.
Robbins MM, Gray M, Kagoda E, Robbins AM. Population dynamics of the Bwindi mountain gorillas. Biol Conserv. 2009;142:2886–955.
Worton BJ. Kernel methods for estimating the utilization distribution in home-range studies. Ecology. 1989;70:164–8.
Calenge C. The package “adehabitat” for the R software: a tool for the analysis of space and habitat use by animals. Ecol Model. 2006;197:516–9.
R Core Team. R: a language and environment for statistical computing. Vienna: R Foundation for Statistical Computing; 2017. https://www.R-project.org.
Van Winkle W. Comparison of several probabilistic home-range models. J Wildl Manag. 1975;39:118.
Bhattacharyya A. On a measure of divergence between two statistical populations defined by their probability distributions. Bull Calcutta Math Soc. 1943;35:99–109.
Clapp JG, Beck JL. Evaluating distributional shifts in home range estimates. Ecol Evol. 2015;5:3869–78.
Fieberg J, Kochanny CO. Quantifying home-range overlap: the importance of the utilization distribution. J Wildlife Manage. 2005;69:1346–59.
Guschanski K, Vigilant L, McNeilage A, Gray M, Kagoda E, Robbins MM. Counting elusive animals: comparing field and genetic census of the entire mountain gorilla population of Bwindi Impenetrable National Park, Uganda. Biol Conserv. 2009;142:290–300.
Baddeley A, Turner R. Spatstat: an R package for analyzing spatial point patterns. J Stat Softw. 2006;12:1–42.
Bivand R, Rowlingson B, Diggle P, Petris G, Eglen S. Splancs: spatial and space-time point pattern analysis. R. 2016. https://www.maths.lancs.ac.uk/~rowlings/Splancs/. Accessed 24 Feb 2018.
VanDerWal J, Falconi L, Januchowski S, Shoo L, Storlie C. Species distribution modelling tools: tools for processing data associated with species distribution modelling exercises. R. 2014. https://www.rforge.net/SDMTools/. Accessed 24 Feb 2018.
Baayen RH. Analyzing linguistic data: a practical introduction to statistics using R. Cambridge: Cambridge University Press; 2008.
Wrangham RW, Gittleman JL, Chapman CA. Constraints on group size in primates and carnivores: population density and day-range as assays of exploitation competition. Behav Ecol Sociobiol. 1993;32:199–209.
Seiler N, Robbins MM. Factors influencing ranging on community land and crop raiding by mountain gorillas. Anim Conserv. 2016;19:176–88.
van de Pol M, Wright J. A simple method for distinguishing within-versus between-subject effects using mixed models. Anim Behav. 2009;77:753–8.
Bates D, Mächler M, Bolker B, Walker S. Fitting linear mixed-effects models using lme4. J Stat Softw. 2015;67:1–48.
Schielzeth H. Simple means to improve the interpretability of regression coefficients. Methods Ecol Evol. 2010;1:103–13.
Schielzeth H, Forstmeier W. Conclusions beyond support: overconfident estimates in mixed models. Behav Ecol. 2009;20:416–20.
Barr DJ, Levy R, Scheepers C, Tily HJ. Random effects structure for confirmatory hypothesis testing: keep it maximal. J Mem Lang. 2013;68:255–78.
Field A. Discovering statistics using SPSS. London: Sage Publications; 2005.
Fox J, Weisberg S. An R companion to applied regression. Thousand Oaks: Sage Publications; 2011.
Forstmeier W, Schielzeth H. Cryptic multiple hypotheses testing in linear models: overestimated effect sizes and the winner’s curse. Behav Ecol Sociobiol. 2011;65:47–55.
We are grateful to the Uganda Wildlife Authority (UWA) and the Uganda National Council for Science and Technology for permission to conduct research in Uganda and for support for this work. We are greatly indebted to many researchers, field assistants and UWA staff who have contributed to the long-term demographic database and GPS data collection in Bwindi Impenetrable National Park, Uganda. We are grateful to the Institute for Tropical Forest Conservation for providing logistical support, to R. Mundry for statistical support and to E. Wright for helpful comments on earlier versions of this manuscript.
The Max Planck Society provided funding for data analysis and writing up the project.
Ethics approval and consent to participate
This research was conducted in compliance with the regulations of the Uganda National Council for Science and Technology and the Uganda Wildlife Authority.
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Annual kernel home range and core area sizes for the Bwindi gorilla study groups.
Summary of the mixed model results investigating the impact of group size and location as a proxy for local food availability on annual home range and core area size in Bwindi mountain gorillas.
Tables of model stability of all investigated models. Discussion of the inclusion of Nkuringo Group, which ranged outside of the national park for a notable amount of the study.
Similarity in space use between two neighboring groups using Bhattacharyya affinity of any two neighboring groups in the same year.
Random slope structures of all investigated models.
About this article
Cite this article
Seiler, N., Robbins, M.M. Using long-term ranging patterns to assess within-group and between-group competition in wild mountain gorillas. BMC Ecol 20, 40 (2020). https://doi.org/10.1186/s12898-020-00306-6