- Research article
- Open Access
Breaking down population density into different components to better understand its spatial variation
BMC Ecology and Evolution volume 21, Article number: 82 (2021)
Population size and densities are key parameters in both fundamental and applied ecology, as they affect population resilience to density-dependent processes, habitat changes and stochastic events. Efficient management measures or species conservation programs thus require accurate estimates of local population densities across time and space, especially for continuously distributed species. For social species living in groups, population density depends on different components, namely the number of groups and the group size, for which relative variations in space may originate from different environmental factors. Whether resulting spatial variations in density are mostly triggered by one component or the other remains poorly known. Here, we aimed at determining the magnitude of the spatial variation in population densities of a social, group-living species, i.e. the European badger Meles meles, in 13 different sites of around 50 km2 across France, to decipher whether sett density, group size or proportion of occupied sett variation is the main factor explaining density variation. Besides the intrinsic factors of density variation, we also assessed whether habitat characteristics such as habitat fragmentation, urbanisation, and resource availability, drove both the spatial variation of density components and local population densities.
We proposed a new standardised approach combining use of multiple methods, namely distance sampling for estimating the density of occupied sett clusters, i.e. group density, and camera and hair trapping for genetic identification to determine the mean social group size. The density of adult badgers was on average 3.8 per km2 (range 1.7–7.9 per km2) and was positively correlated with the density of sett clusters. The density of adult badgers per site was less related to the social group size or to the proportion of occupied sett clusters. Landscape fragmentation also explained the spatial variation of adult badger density, with highly fragmented landscapes supporting lower adult densities. Density components were linked differently to environmental variables.
These results underline the need to break down population density estimates into several components in group-living species to better understand the pattern of temporal and spatial variation in population density, as different components may vary due to different ecological factors.
Estimating population size and density is at the core of both fundamental ecology (e.g. population dynamics over space and time) and applied ecology (e.g. wildlife management). Indeed, density-dependent effects and responses to resource variation, gradual habitat changes or stochastic environments jointly affect population size and density, and ultimately population resilience . Reliable estimates of population size and density in space and/or time are often necessary either to preserve (endangered species such as tigers ), control (pests or vectors of disease, for example badgers across Europe  and several bat species ) or understand the evolutionary ecology of wild populations. Moreover, for species distributed over large ranges, obtaining accurate estimates of global and local abundance in multiple populations and understanding the pattern of their spatial variation are essential to implement appropriate management and conservation measures. For example, spatial variation in population abundance has been linked to the satisfaction of species-specific ecological requirements at local sites [5, 6] or to the location of the site relative to the centre of the species distribution area [5, 7]. However, deriving accurate estimates of population size and density can be difficult. Detecting individuals may be difficult in itself (for cryptic or nocturnal species ), and the species may be continuously present over such a large range that its population density greatly vary in space depending on the local environment .
When species have a wide spatial distribution, local conditions likely affect local population densities. To understand and predict spatial variation in population density at the large scale, it is crucial to look for ecological factors (e.g. habitat quality) influencing sub-populations at the local scale. Availability of food resources and habitat fragmentation strongly affect local settlement/refuge [10, 11]. The social structure of group-living species may also be complex, influencing both the number of social groups and group size. Local environmental conditions may differently affect group sizes and number of territories, i.e. the number of social groups , so that the resulting population density variation may be more or less triggered by variation in both density components. Therefore, to better understand population density variation of widespread social species, a possible approach consist of: (i) estimating group sizes and number of groups at multiple sites, (ii) investigating which component influences population density the most, and (iii) identifying which environmental factors affect the observed densities and their components.
Several methods are commonly used to estimate population density in wildlife ecology, depending on the scale of interest (live trapping combined with capture–mark–recapture analyses, direct observations of individuals through distance sampling methodology, indirect non-invasive data collection methods like detecting signs of presence, genetic tools, and camera traps). Nevertheless, each of these methods presents some limitations and, to optimally choose among them, it is crucial to take into account the ecology of the focus species (e.g. social system, activity patterns, foraging behaviour, etc. ). This is particularly important in social species where densities depend not only on territory sizes and the number of groups in a given area, but also on group size . However, even if the mere presence of a group of individuals can be easily detected through direct observations, estimating its size in a particular area can be tedious, as individuals must be visually differentiated, which is not always possible based on phenotype . In that respect, the development of genetic tools (such as microsatellite markers) allows us to circumvent this difficulty and confidently identify individuals based on their genotype. Combining direct methods to localise areas of interest and estimate their density (e.g. territory, sett, or any clue indicating the presence of a group), and indirect methods to differentiate individuals, and therefore count individuals of a given group, offers a promising approach for estimating population density in social species .
One example of a social and widely distributed mammal is the European badger Meles meles, which occupies a broad range of environments across Europe . Badgers live in territorial and mixed-sex social groups, inhabiting one or several setts in their territories . Using different methodological approaches (see  for a synthesis across Europe), various authors have shown that badgers accommodate to contrasting environmental conditions, and that their densities, group size, proportion of breeding setts, and group territory size, vary both among and within regions as a consequence of variation in factors such as the availability of preferred habitat . Nevertheless, to our knowledge, no study has broken down population density into its different components, namely sett density, group size, and percentage of occupied setts, based on a comprehensive set of local populations using up-to-date methodological approaches. Such a standardised approach might be easily exported to other social and group-living species.
In the present study, we aimed at determining badger densities in 13 different areas of metropolitan France, using a robust combination of methodological approaches. In moderate to high density populations, as the territory of a social group can encompass several setts [20,21,22], we grouped nearby setts together and considered clusters of setts instead of setts. We first estimated the density of sett clusters of each study site using distance sampling [23, 24] on walked transects, and then estimated the proportion of occupied sett clusters (an approximation of the number of social groups in a given area). Second, we determined mean social group size using two non-invasive methods, i.e. camera trapping and hair trapping for genetic identification. These methods have been previously used to estimate mean social group size in badger populations (see for camera trapping [25, 26] and for hair trapping [27, 28]). Finally, we estimated the density of badger as: D = A × B × C/Area, with (A) the number of sett clusters, (B) the proportion of occupied sett clusters, (C) the number of animals per occupied sett cluster, and Area the size of the study site.
While badgers are present throughout metropolitan France in varying abundances across regions , we expected contrasting local densities as a result of significant variations in forest cover and human pressure, leading to quantitative and qualitative habitat variation between study sites. We investigated whether spatial variation in population density was mainly due to variation in social group size or in sett cluster density. Besides these intrinsic factors of variation, we also investigated whether habitat characteristics, i.e. the magnitude of forested habitat fragmentation and resource availability, are correlated to the local population density, as shown at the regional scale throughout France .
In each of the 13 study sites, 76.3 to 93.4 km were walked. Over the 13 study sites, realized transect length ranged from 0.4 to 3.2 km (mean = 1.7 km ± 0.4 SD), which was quite similar to planned transects (difference being attributed to walking difficulties on the field mostly due to vegetation cover). A total of 533 setts were detected (mean = 41.0 ± 19.2 SD), among which 231 setts were surveyed to estimate mean badger group sizes. These 533 setts were grouped into 273 sett clusters using the between-setts distance of 500 m (Table 1). Of the 273 sett clusters found, 96 were composed of unoccupied setts (35.2%), 57 included only secondary setts with at least one occupied (i.e. secondary sett clusters, SSC; 20.9%), 82 included at least one occupied main sett without reproduction (i.e. main sett clusters, MSC; 30%) and finally 38 included at least one occupied main sett with reproduction (i.e. main sett clusters with reproduction, MSCR; 13.9%).
Sett cluster density
We first modelled a detection function from the histogram of the observed perpendicular distance data (Additional file 1). Among the various detection functions tested, the hazard rate (HR) function had the lowest AIC (Table 2) and fitted well the data (χ2 = 1.58, df = 4, p-value = 0.81). The top-ranked model used a post-stratification using the habitat type (i.e. forested or hedgerows sites) of the whole study site SIT (18.6 drop in AIC when compared to pooled data model). Hence, we observed a significantly higher detectability and effective strip width in forested sites than in hedgerow sites (19.7 ± 4.2 m with standard deviation and 10.7 ± 1.2 m, respectively). All other variables tested i.e. habitat type along the transect (HAB), the type of sett cluster (TER), and multiple covariates distance sampling did not improve the model (Table 2).
The coefficients of variation of sett cluster density in suitable habitat DC.Distance ranged from 22 to 40% (derived by bootstrap; n = 999), with most of the variability being due to variability in encounter rates between transects. Sett cluster density DC calculated by multiplying DC.Distance by the proportion of suitable habitat in each study site ranged from 1.99 (site L) to 6.42 per km2 (site H; Table 3).
Among-site variation in social group size and proportion of occupied sett clusters
In the 13 study sites, mean trapping effort was 515.2 ± 77.4 SD camera-trap days and 881.0 ± 300.3 SD hair-trap days per site. Based on the camera trapping survey, the mean number of adults per cluster was higher in MSCR (adMSCR: 1.95 ± 0.43 SD, max = 3; site A) than in MSC (adMSC: 1.37 ± 0.35, max = 2.25; site D) or SSC (adSSC: 1.34 ± 0.41, max = 2; site C, Table 4). The average social group size, including adults and cubs, was also higher in MSCR (badgerMSCR: 4.46 ± 1.01, max = 5.75; site L) than in MSC (badgerMSC: 1.70 ± 0.72, max = 3.50; site B) or SSC (badgerSSC: 1.68 ± 0.80, max = 3.50; site C, Table 4). At three study sites (E, I, and L; Table 4), neither camera trapping nor hair trapping allowed us to capture individuals in secondary sett clusters, so estimates of group size for these sites were calculated as the average of the estimated adSSC and badgerSSC for the same type of clusters (SSC) obtained in the same type of study site (i.e. forest or hedgerow sites). Similarly, we also used the average value of adMSCR and badgerMSCR for study site M, where no main sett clusters with reproduction were found (Table 4).
The proportion of occupied sett clusters varied greatly among sites. For SSC, pSSC ranged from 6% (sites E and J) to 52% (site D; mean = 18 ± 14 SD; Table 1). For MSC, pMSC ranged from 13% (site A) to 82% (site L; mean = 33 ± 20 SD), while for MSCR, pMSCR ranged from 0% (site L and M) to 33% (site C; mean = 13 ± 10 SD; Table 1).
Badger population densities in France
The density of adult badgers (DAd) ranged from 1.66 (site J) to 7.86 per km2 (site D; Fig. 1) with an average of 3.84 ± 1.75 (SD) across study sites. Considering the total number of individuals detected per cluster, the density of adults and cubs (DBad) ranged from 2.41 (site L) to 13.29 per km2 (site B; Fig. 1), with an average of 5.85 ± 3.25 across study sites.
Determinants of spatial variation in density
The sum of the eigenvalues for the first three axes of the principal component analysis (PCA) explaining the pattern of variability of DAd accounted for 66% of total variability (Fig. 2a). The first axis of this PCA accounting for more than 30% of the total variability expressed a fragmentation gradient between study sites, from fragmented landscapes (with high edge density) to more suitable habitat for badger settlement (corresponding to a greater proportion of forest, forest edge, and hedgerows; Suit.area). We identified seven co-variables likely to be associated with DAd (Fig. 2a) and for which we explicitly tested the correlation, namely: Edge.density, Pasture, Suit.area, adMSCR, adMSC, pSSC, and DC. Considering the intrinsic components, the density of adult badgers (DAd) was positively correlated with sett cluster density DC (rho = 0.76; p-value < 0.007 accounting for Bonferroni correction), but not significantly correlated with group sizes adMSCR and adMSC or the proportion of occupied secondary sett clusters pSSC (all p-values > 0.007; Fig. 2a). PCA revealed that fragmented landscapes (with a high edge density) tend to support lower adult densities (rho = − 0.62; p-value = 0.027) and that suitable habitats tend to support higher badger densities (rho = 0.56; p-value = 0.05). No significant correlation was found for the presence of pasture area (rho = − 0.489; p-value = 0.093).
With regard to DBad (see Fig. 2b and details in Additional file 2), we only found a significant positive correlation with Edge.density (rho = − 0.73; p-value < 0.008 accounting for Bonferroni correction).
Badger density in European countries
Densities of adult badgers (DAd) ranged from 1.66 to 7.86 per km2 (average of 3.84), and DBad including adults and cubs ranged from 2.41 to 13.29 per km2 (average of 5.85) across the 13 study sites in France. Comparing badger densities across studies is inherently difficult as the methods used (distance sampling, radio-tracking, visual census, capture-recapture, group size estimation), calculations (type of sett or sett cluster accounted for or not, proportion of occupied setts accounted for or not), and estimated parameters (adult density, badger density without discrimination between adults and cubs, sett density, occupied sett density) differ among studies (see  for a detailed description of badger densities estimated with different methods across Europe). In the UK, populations can reach densities greater than 20 adults per km2 , but are around 9.4 individuals per km2 on average (ranging from 0.86 to 30.7; reviewed in [18, 20, 31, 32]. In Central Europe, densities of less than 5 individuals per km2 are more common [15, 33]. In Mediterranean and boreal environments, densities are lower, with badgers living in pairs or small groups at densities of less than 0.7 adults per km2, and in some countries densities are even lower than 0.01 individuals/km2 [31, 34]. Overall, our estimates were thus relatively lower than those found in the UK and concordant with global estimates from Ireland (0.72 to 11.9 adults/km2 depending on the region, ). In France, a few studies have estimated densities, most often in small areas (< 10 km2), with estimates ranging from 0.1 to 1.6 adults/km2 [36, 37], and group sizes between two and four badgers including young. Rigaux and Chanu  found a density of 1.9 badgers/km2 in a 58 km2 area, partly overlapping the K study site. Our estimates are clearly higher, which is mainly related to the high cluster density estimates. Cluster densities DC ranged from 1.99 to 6.42 per km2 and, when taking into account the proportion of occupied main setts (ranging from 0.25 to 0.82), the densities of occupied main clusters roughly corresponded to the density of family groups, ranging from 0.76 to 2.63 per km2 (mean 1.72 ± 0.52 SD).
A novel methodological approach to estimate density
Various methodological approaches have been used in previous studies to estimate badger densities; all present pros and cons. The main difference in our approach is the distance sampling method we used, which permits non-exhaustive prospection of a study area. However, we found small effective strip width in the distance sampling analysis, due to short distances between spotted setts and transects; this might have resulted in overestimated cluster densities, albeit it was comparable between study sites. The other originality and strength of our study lies in the combined use of multiple methodological approaches, namely distance sampling, camera trapping, and genetic sampling. For example, group sizes are necessarily underestimated using only camera traps as it is highly unlikely that all individuals are seen simultaneously. Here, we partially circumvented this risk by combining the use of genetic sampling, to distinguish individuals, and camera trapping to distinguish adults from cubs. Furthermore, this allowed us to fully decompose the estimation of badger density in sets of components to better understand spatial variation in population density. Our methodological approach could be applied to the long-term monitoring of a single population to better grasp temporal variation in badger density in future studies. Beyond this, for group-living and social species that live in burrows (such as meerkats Suricata suricatta) or do not live in burrows (such as numerous primates), composite estimates of population density over a large scale might help us better understand population dynamics, population genetics, and evolutionary ecology.
We acknowledge that our analyses may suffer from some biases, but we tried to minimise them. First, we estimated sett occupancy over a short time period (1 week) and used a standardised activity scoring system, to reduce the risk of overestimation, even though badgers regularly move between setts (especially secondary ones), and occupation cues (for example faeces, traces of digging activities) persist over a longer time period. As a result, our estimated proportion of main setts (17.3%) versus secondary and unoccupied setts (82.7%) is very similar to estimates in other populations (e.g. 22.2% for main setts and 77.8% for secondary and unoccupied setts across six studies performed in Ireland, reviewed in ). Another methodological limitation may be the choice of a between-setts distance of 500 m to group setts into sett clusters. The selected distance necessarily affects the number of social groups and, indeed, density estimates. Using two alternative clustering solutions (i.e. between-setts distance of 100 and 900 m), we partially confirmed this overall trend, but estimates based on these two alternative clustering solutions fell within the confidence interval of the current density estimates. We also found a strong positive correlation between badger densities calculated with the 500 m distance and the 100 m distance (rho = 0.967 and 0.918 for DAd and DBad respectively) or the 900 m distance (rho = 0.841 and 0.967 for DAd and DBad respectively; Additional file 3), supporting the rank of densities among study sites and thus PCA results.
On the relative importance of intrinsic density components and their link with ecological variables
Decomposition of density estimation into a set of intrinsic components enables a better understanding of the spatial variation in badger population densities. Among these components, variation in sett cluster density drives the variation in adult density DAd, while group size in main sett clusters (with and without reproduction) and the proportion of occupied secondary sett clusters are less related. In fact, group sizes are quite constant in France (see Table 4) in comparison to other European countries [19, 35]. Therefore, a better proxy for DAd might be sett cluster density, and researchers might focus on this parameter when seeking badger density estimates with a reduced version of our standardised protocol due to logistic constraints. PCA results have also shown that DAd is associated with ecological/extrinsic components. A higher proportion of suitable habitat in a study site tends to lead to greater DAd, while more fragmented areas tend to decrease density. The major intrinsic component here that triggers spatial variation in density, i.e. sett cluster density, is thus driven by landscape fragmentation and thus the availability of suitable sett sites. Indeed, the highest DAd densities were found in sites dominated by large forest patches (study sites A, B, D, and H). The great importance of habitat structure for population density has also been demonstrated in other mammal species (group size of blackbuck antelope, Antilope cervicapra, increases with increasing habitat openness ) and bird species (abundance is mainly affected by habitat cover in most neotropical bird species ). Yet the underlying mechanisms (e.g. decreased vigilance costs, increased resource availability) often remain difficult to target without breaking down the density estimate into its structural components, as we performed here. Interestingly, the proportion of urban areas and density of roads over study sites are not related to adult badger density, but to some intrinsic components of DAd density. The proportion of occupied main sett clusters without reproduction (pMSC) increases in more urbanised study sites with higher road density (e.g. site L), suggesting that anthropic presence is less favourable to reproduction, whereas the proportion of main sett clusters with reproduction (pMSCR) was associated with earthworm abundance, suggesting an association between the amount of available resources and probability that a social group will reproduce efficiently. The importance of resource availability for badger settlement was already clear from an earlier study showing that habitats with a large earthworm biomass support high badger densities . This result is in line with other empirical studies on social species, showing that an important component of population abundance, namely group size, is correlated with food resource availability (in caviomorph rodents, reviewed in ; in wolves, Canis lupus ; in blackbuck antelope, Antilope cervicapra ). High habitat quality, the availability of food resources in particular, may promote site occupancy and improve individual conditions, thereby increasing both components of individual or group fitness, i.e. reproductive success and survival. Maximum population density could thus be achieved in places where each density component is maximised, with each of these depending on different environmental factors. Whether or not these optima in different ecological factors can be reached together in some places remains unknown, but would help us better understand the limits of observable population density . This also would allow us to better define the ecological requirements of a given species and better appraise on which environmental factor to focus to derive appropriate management and conservation measures.
The composite estimation of badger density, in several study sites and over a large spatial scale, enabled a better understanding of their pattern of spatial variation. Densities of adult badgers (DAd) ranged from 1.66 to 7.86 per km2 across the 13 study sites in France. We found that sett cluster density variation drives adult density DAd. A higher proportion of suitable habitat tends to lead to greater DAd, while more fragmented areas tend to decrease adult badger density. Across the geographical range of widely distributed species, habitat and environmental conditions are likely to be diverse, imposing contrasting pressures on local population abundances, as shown in various taxa (in birds [5, 7]; in 11 species of birds, mammals and reptiles ; in a meta-analysis ). Breaking down population density in order to identify which components among intrinsic and extrinsic factors trigger the variation in local density is all the more important for social species, in which population densities depend on several factors: territory size, number of social groups, and group size . Identifying ecological drivers of variation in each of the key components of population abundance is a cornerstone of population ecology, including wildlife management of pest/invasive species and the conservation biology of threatened taxa.
The study was conducted at 13 sites spread throughout France (Fig. 3) between 2014 and 2018. Study sites ranged from 46.55 to 59.31 km2. Each site was composed of different proportions of suitable habitat (Suit.area) for badger settlement (i.e. forests, forest edges and hedgerows; details in Additional file 4). Two kinds of landscape dominated in the various study sites: forest in A, B, D, and H study sites (at least 70% forested areas) and hedgerows in all other study sites (less than 30% forest cover). Various other environmental variables, such as the proportion of urbanised areas, edge density, and various resource variables (proportion of maize, earthworm availability; see Additional file 4) contrasted the 13 study sites.
Using distance sampling to estimate sett cluster density
Sett detection during walked transects
We used a systematic sampling method by positioning over each study site an array of equidistant theoretical points spaced 1 km apart. Transects were then placed in close vicinity to each theoretical point in suitable habitats for badger settlement (i.e. forests, forest edges and hedgerows). As sett density was expected to be very low in open areas , we did not sample this kind of habitat, acknowledging that it might lead to a slight underestimation of population density. Locations of transects were chosen according to the relative proportions of the three types of suitable habitat within each study site. We used detailed maps (IGN BD Topo) of the study sites on QGIS v. 2.18  to plan 1- to 2-km long transects to walk (see examples of maps with theoretical points and transects in Additional file 4).
Depending on study site, 48 to 53 transects (mean = 50.92 ± 1.50 SD) were surveyed once between March and April (5 weeks), i.e. before vegetation peak, to facilitate sett detection. Realized transects were recorded using a field GPS (Garmin GPSMAP 64s). As badger setts might have numerous entrances and their surface size might vary greatly , we recorded the GPS coordinates of the sett centroid. Vegetation cover prevented a precise measurement of perpendicular distances between setts and transects, so we calculated perpendicular distances from tracks and sett centroids on QGIS.
Distinction between different sett types and sett grouping into sett clusters
There are different types of setts in badger populations: those commonly described as main setts (large and permanently occupied by all members of a social group, and used for breeding), and other setts intermittently occupied by badgers (i.e. classified as annex, subsidiary and outlier setts ). Classifying these intermittently occupied setts in the field requires multiple visits and follow up, so we grouped them as ‘secondary’ setts [17, 50].
Therefore, in April, after the 5-week survey to detect setts, we visited all the georeferenced setts again during a 5-day field session to assess activity (such as number of entrances or well-used badger paths, latrines, or ‘playing areas’ with denuded vegetation caused by cubs and adults chasing each other). Based on these signs of activity, we developed a standardised activity scoring system for all the study sites, and controlled for observer bias (Additional file 5). Estimated activity scores allowed us to distinguish between occupied main setts, occupied secondary setts, and unoccupied setts.
We also had to consider that the territory of a social group could encompass several setts of different types , forming a ‘cluster’ of setts. We thus grouped the nearest setts together and considered sett clusters rather than setts. Based on both the studies of Carter et al. ([22, 52]), and on the distance distribution between the detected setts in all study sites (Additional file 4), we considered a between-setts centroid distance of 500 m. This distance seemed a reliable distance to group setts into clusters while excluding setts from neighbouring social groups. As a result, clusters could be composed of one or several detected setts of different types (see maps in Additional file 4). It is also well-known that not all females and social groups produce cubs every year . Hence, we considered four categories: sett clusters with (i) at least one occupied main sett with reproduction as main sett clusters with reproduction (MSCR), (ii) at least one occupied main sett without reproduction as main sett clusters (MSC), (iii) only secondary setts with at least one occupied as secondary sett clusters (SSC), and (iv) only unoccupied setts. Those categories were then used for calculations of sett clusters density and badger group sizes. We also determined per study site the proportions of the occupied main sett clusters with reproduction (pMSCR), occupied main sett clusters without reproduction (pMSC), and occupied secondary sett clusters (pSSC; Table 1).
Distance sampling analysis
We used distance sampling methodology [23, 24], in which the probability of detection is modelled from the frequency distribution of the perpendicular distances of detected clusters from the transect line. We calculated the perpendicular distance between the first detected sett centroid in the cluster and the walked transect using QGIS v 2.18 . Following the recommendations of Buckland et al. , we right truncated the data by removing 3.3% of observations that were more than 55 m from the transect. We tested five a priori robust models to fit the probability of detection: a uniform key function with either cosine or polynomial series expansion, a half-normal with either cosine or Hermine polynomial series expansion and a HR key function with cosine series expansion. We assessed their goodness of fit visually and through Chi-squared tests. In modelling the detection function, we expected that several covariates could a priori affect detection probability: (i) HAB, the habitat type along the transect (i.e. forest, forest edge or hedgerow), (ii) SIT, the habitat type of the whole study site (i.e. forested or hedgerow sites), and (iii) TER, the type of sett clusters. We used multiple-covariates distance sampling (MCDS, [54, 55]) and tested post-stratification with the conventional distance sampling (CDS) method. MCDS tested whether the scale parameter (and hence the shape of the detection function) varied between levels of a covariate, whereas CDS tested for different detection functions (shape and origin). We retained the best function model ranked by the Akaike information criteron (AIC).
To estimate the density of sett clusters in suitable habitat by study site (DC.Distance; Table 3), we used stratum-specific detection probabilities. The effective strip width was based on the best model, and encounter rates were calculated for each study site. We derived confidence intervals from a non-parametric bootstrap method, resampling transects within each study site to estimate bootstrap variances of the encounter rates (with n = 999). We then derived the density of sett clusters in each site (DC, Table 3) by multiplying each DC.Distance by the proportion of suitable habitat in the corresponding study site. All analyses were conducted in DISTANCE software 7.3 . A DISTANCE project containing the data and analyses is available from the authors.
Determination of badger group size per sett cluster
In order to estimate social group size of badgers per cluster type (i.e. MSCR, MSC, and SSC), we used two non-invasive methods, i.e. camera trapping and hair traps for genetic identification. Due to material and human limitations, we were unable to follow all occupied setts detected in the study. We selected on average 18 occupied setts (± 4 SD [min = 9; max = 27]) from each study site (with approximately half of the main setts and half of the secondary setts). Between April and June, we deployed camera traps on these selected setts for 2 weeks and hair traps for 2 another weeks, alternating the methods between setts to avoid bias. Due to the small number of occupied setts detected in seven study sites, we had to include setts found outside the walked transects protocol, that were known to local people (see Table 1 for the total number of setts used per site). Depending on the study site, the number of occupied setts and their spatial localisation, several setts of the same cluster were surveyed. Badger group size per sett cluster corresponds to the maximum number of animals detected among all setts in the cluster.
Camera trapping survey
We monitored setts using one to five infrared cameras, depending on the sett type (mean = 2.41 ± 0.90 for main setts; and 1.47 ± 0.74 for secondary setts). Camera traps were tied to trees about 5 m from active entrances, badger paths, and preferably in any places with recent signs of activity and gathering (for example ‘grooming and playing areas’), in order to detect all adults and cubs present in the sett. We checked for battery depletion and proper functioning of each camera trap twice during the 2-week survey (after 2 and 8 days in the field).
The camera trap method allowed us to distinguish adults and cubs based on body size, and thus estimate the mean number of adults per cluster type (adMSCR, adMSC, and adSSC; Table 4). For a given sett cluster, we took the maximum number of adults seen together on both pictures and videos. These social group sizes were minima, as we did not have the absolute certainty that all individuals living in a sett cluster were camera-trapped.
Hair trapping survey for genetic identification
We sampled badger hairs using hair traps, i.e. barbed wire (thickness of 1.7 mm, with barbs spread every 10 cm) suspended over visible badger paths, approximately 20 cm above ground level, to subsequently identify individuals genetically. Hair samples were genotyped at 24 microsatellite markers and one sex marker (full protocol described in ). We identified a total of 284 badgers (145 males and 139 females). Genetic identification did not allow us to distinguish between adults and young individuals. As for the camera trap estimates, hair trap estimates were minima, and we selected the maximum of both estimates, i.e. from camera traps and genetic identification, as the badger group size per cluster type (badgerMSCR, badgerMSC and badgerSSC; Table 4). The study site G conducted in 2014 was a pilot study site, and the hair trapping survey could not be put in place, so the mean badger group size was estimated only from camera-trap results at this site.
Badger density estimates
For each study site, we estimated adult density DAd and badger density DBad (i.e. adults and cubs) per square kilometre as composites of their different components, as follows:
where DC is the density of sett clusters in km2 estimated by distance sampling and corrected for the proportion of suitable habitat in each site; pSSC, pMSC, and pMSCR are the proportions of occupied SSC, MSC, and MSCR among all sett clusters; adSSC, adMSC, and adMSCR are the mean number of adults per SSC, MSC, and MSCR surveyed in the corresponding study site; and badgerSSC, badgerMSC, and badgerMSCR are the mean social group size in badgers (adults and cubs) per SSC, MSC, and MSCR, as above.
Determinants of spatial variation in density
In order to identify which ecological factors influenced spatial variation in density, we performed a PCA using the FactomineR package  operating in R software . We included in the PCA adult density (DAd), each intrinsic component used in Eq. 1 (n = 7), and extrinsic factors linked with local environmental characteristics (n = 10). The environmental variables included were Suit.area, Edge density (m/ha), VRM (°), Soil texture and depth (indices), Earthworm abundance (ind./m2), percentage of Pasture, Maize, Urbanised area in the study site, and Road density (km/km2), defined in Additional file 4. PCA allowed us to identify a reduced set of variables, both intrinsic and extrinsic, likely to be associated with DAd; we further explicitly tested for the association using non-parametric Spearman’s rank correlation tests (R package pspearman, ). p-values were adjusted using the Bonferroni procedure for multiple comparisons  on the reduced set of variables. We applied the same procedure on badger density, adults and cubs, by substituting DAd by DBad and using the intrinsic components in Eq. 2. Coefficients of variation (CV) of each variable included in the PCAs are available in Additional file 2.
Availability of data and materials
The datasets used and/or analysed during the current study are available from the corresponding author on reasonable request.
- ad SSC :
Mean number of adult badgers in secondary sett cluster
- ad MSC :
Mean number of adult badgers in main sett cluster without reproduction
- ad MSCR :
Mean number of adult badgers in main sett cluster with reproduction
- badger SSC :
Mean number of adult and cubs in secondary sett cluster
- badger MSC :
Mean number of adult and cubs in main sett cluster without reproduction
- badger MSCR :
Mean number of adult and cubs in main sett cluster with reproduction
Conventional distance sampling
- D Ad :
Adult badger density
- D Bad :
Adult and cubs density
- D C.Distance :
Cluster density estimated in suitable habitats of each site
- D C :
Cluster density in each site, i.e. DC.Distance × [proportion of suitable habitat in the site]
Habitat type along the transect (i.e. forest, forest edge or hedgerows)
Multiple covariates distance sampling
Main sett cluster without reproduction
Main sett cluster with reproduction
Principal component analysis
- pSSC :
Proportion of occupied secondary sett clusters
- pMSC :
Proportion of occupied main sett clusters without reproduction
- pMSCR :
Proportion of occupied main sett clusters with reproduction
Habitat type of the whole study site (i.e. forested or hedgerow sites)
Secondary sett cluster
Type of sett cluster (i.e. unoccupied, secondary or main sett cluster)
Saether BE. Environmental stochasticity and population dynamics of large herbivores: a search for mechanisms. Trends Ecol Evol. 1997;12:143–9.
Karanth KU, Nichols JD, Kumar NS, Hines JE. Assessing tiger population dynamics using photographic capture-recapture sampling. Ecology. 2006;87:2925–37.
Griffiths HI, Thomas DH. The status of the Badger Meles meles (L., 1758) (Carnivora, Mustelidae) in Europe. Mamm Rev. 1993;23:17–58.
Langwig KE, Frick WF, Bried JT, Hicks AC, Kunz TH, Marm KA. Sociality, density-dependence and microclimates determine the persistence of populations suffering from a novel fungal disease, white-nose syndrome. Ecol Lett. 2012;15:1050–7.
Brown JH, Mehlman DW, Stevens GC. Spatial variation in abundance. Ecology. 1995;76:2028–43.
Martínez-Meyer E, Díaz-Porras D, Peterson AT, Yáñez-Arenas C. Ecological niche structure and rangewide abundance patterns of species. Biol Lett. 2013;9:20120637.
Osorio-Olvera L, Yañez-Arenas C, Martínez-Meyer E, Peterson AT. Relationships between population densities and niche-centroid distances in North American birds. Ecol Lett. 2020;23:555–64.
Wilson GJ, Delahay RJ. A review of methods used to estimate the abundance of terrestial carnivores using field signs and observation. Wildl Res. 2001;28:151–64.
Lawton J. Range, population abundance and conservation. Trends Ecol Evol. 1993;8:409–13.
Virgós E. Role of isolation and habitat quality in shaping species abundance: a test with badgers (Meles meles L.) in a gradient of forest fragmentation. J Biogeogr. 2001;28:381–9.
Revilla E, Palomares F. Spatial organization, group living and ecological correlates in low density populations of Eurasian badgers, Meles meles. J Anim Ecol. 2002;71:497–512.
Delahay RJ, Carter SP, Forrester GJ, Mitchell A, Cheeseman CL. Habitat correlates of group size, bodyweight and reproductive performance in a high-density Eurasian badger (Meles meles) population. J Zool. 2006;270:437–47.
Schley L, Roper TJ, Schaul M. Distribution and population density of badgers Meles meles in Luxembourg. Mamm Rev. 2004;34:233–40.
Rowcliffe JM, Field J, Turvey ST, Carbone C. Estimating animal desnsity using camera traps without the need for individual recognition. J Appl Ecol. 2008;45:1228–36.
Palphramand KL, Newton-Cross G, White PCL. Spatial organization and behaviour of badgers (Meles meles) in a moderate-density population. Behav Ecol Sociobiol. 2007;61:401–13.
Johnson DDP, Jetz W, Macdonald DW. Environmental correlates of badger social spacing across Europe. J Biogeogr. 2002;29:411–25.
Roper TJ. Badger. London: HarperCollins; 2010.
Lara-Romero C, Virgós E, Revilla E. Sett density as an estimator of population density in the European badger Meles meles. Mamm Rev. 2012;42:78–84.
Neal E, Cheeseman CL. Badgers. London: T and AD Poyser; 1996.
Byrne AW, Paddy Sleeman D, O’Keeffe J, Davenport J. The ecology of the European badger (Meles meles) in Ireland: a review. Biol Environ. 2012;112:105–32.
Judge J, Wilson GJ, Macarthur R, McDonald RA, Delahay RJ. Abundance of badgers (Meles meles) in England and Wales. Sci Rep. 2017;7:276.
Carter SP, Robertson A, Palphramand KL, Chambers MA, McDonald RA, Delahay RJ. Bait uptake by wild badgers and its implications for oral vaccination against tuberculosis. PLoS ONE. 2018;13:e0206136.
Buckland S, Anderson D, Burnham K, Laake J. Distance sampling: estimating abundance of biological populations. London: Chapman and Hall; 1993.
Buckland S, Anderson D, Burnham K, Laake J. Introduction to distance sampling: estimating abundance of biological populations. Oxford, United Kingdom: Oxford University Press; 2001.
Tuyttens FAM, Long B, Fawcett T, Skinner A, Brown JA, Cheeseman CL, et al. Estimating group size and population density of Eurasian badgers Meles meles by quantifying latrine use. J Appl Ecol. 2001;38:1114–21.
Balestrieri A, Cardarelli E, Pandini M, Remonti L, Saino N, Prigioni C. Spatial organisation of European badger (Meles meles) in Northern Italy as assessed by camera-trapping. Eur J Wildl Res. 2016;62:219–26.
Frantz AC, Schaul M, Pope LC, Fack F, Schley L, Muller CP, et al. Estimating population size by genotyping remotely plucked hair: the Eurasian badger. J Appl Ecol. 2004;41:985–95.
Scheppers TLJ, Roper TJ, Frantz AC, Schaul M, Engel E, Breyne P, et al. Estimating social group size of Eurasian badgers Meles meles by genotyping remotely plucked single hairs. Wildl Biol. 2007;13:195–207.
Calenge C, Chadoeuf J, Giraud C, Huet S, Julliard R, Monestiez P, et al. The spatial distribution of Mustelidae in France. PLoS ONE. 2015;10:e0121689.
Jacquier M, Calenge C, Say L, Devillard S, Ruette S. Altitude shapes the environmental drivers of large-scale variation in abundance of a widespread mammal species. Ecol Evol. 2020;10:119–30.
Kowalczyk R, Bunevich AN, Jedrzejewska B. Badger density and distribution of setts in Bialowieza primeval forest (Poland and Belarus) compared to other Eurasian populations. Acta Theriol. 2000;45:395–408.
Parrott D, Prickett A, Pietravalle S, Etherington TR, Fletcher M. Estimates of regional population densities of badger Meles meles, fox Vulpes vulpes and hare Lepus europaeus using walked distance sampling. Eur J Wildl Res. 2011;58:23–33.
Frantz AC, Schley L, Schaul M, Balestrieri A, Roper TJ. Spatial organisation of badgers (Meles meles) in a medium-density population in Luxembourg. Hystrix Ital J Mammal. 2010;21:3–18.
Rosalino LM, MacDonald DW, Santos-Reis M. Spatial structure and land-cover use in a low-density Mediterranean population of Eurasian badgers. Can J Zool. 2004;82:1493–502.
Feore S, Montgomery WI. Habitat effects on the spatial ecology of the European badger (Meles meles). J Zool. 1999;247:537–49.
Bodin C. Partage de l’espace et relation de voisinage dans une population continentale de blaireaux européens (Meles meles). University of Montpellier II; 2005.
Payne A. Rôle de la faune sauvage dans le système muti-hôtes de Mycobacerium bovis et risque de transmission entre faune sauvage et bovins: étude expérimentale en Côte d’Or. University Claude Bernard of Lyon; 2014.
Rigaux P, Chanu C. Densité du blaireau d’Eurasie (Meles meles) et répartition des terriers dans un paysage rural de moyenne montagne, entre Combrailles et Monts Dôme (Puy-de-Dôme, France). Groupe Mammal d’Auvergne. 2011;1:1–32.
Isvaran K. Intraspecific variation in group size in the blackbuck antelope: the roles of habitat structure and forage at different spatial scales. Oecologia. 2007;154:435–44.
Filloy J, Bellocq MI. Spatial variation in bird species abundances: environmental constraints across southern Neotropical regions. Appl Ecol. 2013;14:263–70.
Kruuk H, Parish T. Factors affecting population density, group size and territory size of the European badger Meles meles. J Zool. 1982;196:31–9.
Maher CR, Burger JR. Intraspecific variation in space use, group size, and mating systems of caviomorph rodents. J Mamm. 2011;92:54–64.
Fuller TK, Murray DL. Biological and logistical explanations of variation in wolf population density. Anim Conserv. 1998;1:153–7.
Stephens P, Vieira M, Willis S, Carbone C. The limits to population density in birds and mammals. Ecol Lett. 2019;22:654–63.
Weber MM, Stevens RD, Diniz-Filho JAF, Grelle CEV. Is there a correlation between abundance and environmental suitability derived from ecological niche modelling? A meta-analysis. Ecography. 2017;40:817–28.
Obidziński A, Pabjanek P, Mędrzycki P. Determinants of badger Meles meles sett location in Białowieża primeval forest, northeastern Poland. Wildl Biol. 2013;19:48–68.
QGIS Development Team 2016. QGIS geographic information system. Open Source Geospatial Foundation Project. http://qgis.osgeo.org/.
Roper TJ. Badger Meles meles setts—architecture, internal environment and function. Mamm Rev. 1992;22:43–53.
Thornton PS. Density and distribution of badgers in south-west England—a predictive model. Mamm Rev. 1988;18:11–23.
Do Linh San E. Le blaireau d’Eurasie. Paris: Delachaux et Niestlé; 2006.
Roper TJ. The structure and function of badger setts. J Zool. 1992;227:691–4.
Jacquier M, Vandel J-M, Léger F, Duhayer J, Pardonnet S, Queney G, et al. Population genetic structures at multiple spatial scales: importance of social groups in European badgers. J Mamm. 2020;101:1380–91.
Anderson RM, Trewhella W. Population dynamics of the badger (Meles meles) and the epidemiology of bovine tuberculosis (Mycobacterium bovis). Philos Trans R Soc Lond B Biol Sci. 1985;310:327–81.
Marques TA, Thomas L, Fancy SG, Buckland ST. Improving estimates of bird density using multiple-covariate distance sampling. Auk. 2007;124:1229–43.
Thomas L, Rexstad EA, Laake JL, Hedley SL, Burnham KP, Marques TA, et al. Distance software: design and analysis of distance sampling surveys for estimating population size. J Appl Ecol. 2010;47:5–14.
Lê S, Josse J, Husson F. FactoMineR: an R package for multivariate analysis. J Stat Softw. 2008;25:1–18.
R Development Core Team. R: a language and environment for statistical computing. R Foundation for Statistical Computing, Vienna. 2019. http://www.r-project.org/.
Savicky P. pspearman: sprearman’s rank correlation test. R package version 0.3.0. 2014. http://cran.r-project.org/package=pspearman/.
Bonferroni C. Teoria statistica delle classi e calcolo delle probabilità. Pubbl del R Ist Super di Sci Econ e commericiali di firenze. 1936;8:3–62.
We are very grateful to all the people who participated in the fieldwork at the study sites, in particular the students: T. Beshers, A. Blain, C. De Moncuit, E. Harwood, M. Lemaire, T. Ludemann, S. Malevre, C. Pacot, P. Perrin, J. M. Poincelet; their field supervisors: M. Bouchez-Zacria, F. Boué, A. Courcoul, J. M. Demerson, L. Legras, and C. Richomme (French Agency for Food, Environmental and Occupational Health and Safety; ANSES), C. Peboscq (Departmental Hunting Association in Pyrénées-Atlantiques; FDC64), C. Langumier (National Forestry Office; ONF), A. M. Griffon-Picard (Departmental Direction for the Protection of Populations in Seine-Maritime; DDPP), S. Depraz (GDS 24), and our colleagues M. Albaret, S. Brandt, H. Esclaffer, C. Rousset, L. Tison and J.-L. Vérita from the French Biodiversity Agency (OFB).
Genetic analyses in this research project were partly funded by the Antagene laboratory (La Tour-de-Salvagny, France) and the French Ministry of Agriculture and Food (DGAL). They played no part in the design of the study, analysis and interpretation of data or in writing the manuscript.
Ethics approval and consent to participate
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.
Detection functions for badger sett clusters estimated using distance sampling methodology.
Coefficients of variation of each variable included in the PCAs explaining the pattern of variability of badger density.
Estimations of badger density over the 13 study sites using two alternative clustering solutions.
Description and geographic maps of the 13 study sites surveyed in France.
Standardized scoring system used in order to distinguish badger main setts, occupied secondary setts, and unoccupied setts across all study sites.
About this article
Cite this article
Jacquier, M., Vandel, JM., Léger, F. et al. Breaking down population density into different components to better understand its spatial variation. BMC Ecol Evo 21, 82 (2021). https://doi.org/10.1186/s12862-021-01809-6
- Population density
- Group size
- Meles meles
- Distance sampling
- Molecular ecology
- Camera trap