Skip to main content
  • Research article
  • Open access
  • Published:

Quantifying suitable late summer brood habitats for willow ptarmigan in Norway



Habitat models provide information about which habitat management should target to avoid species extinctions or range contractions. The willow ptarmigan inhabits alpine- and arctic tundra habitats in the northern hemisphere and is listed as near threatened (NT) in the Norwegian red list due to declining population size. Habitat alteration is one of several factors affecting willow ptarmigan populations, but there is a lack of studies quantifying and describing habitat selection in willow ptarmigan. We used data from an extensive line transect survey program from 2014 to 2017 to develop resource selection functions (RSF) for willow ptarmigan in Norway. The selection coefficients for the RSF were estimated using a mixed-effects logistic regression model fitted with random intercepts for each area. We predicted relative probability of selection across Norway and quantile-binned the predictions in 10 RSF bins ranging from low-(1) to high-(10) relative probability of selection.


Random cross-validation suggest that our models were highly predictive, but validation based spatial blocking revealed that the predictability was better in southern parts of Norway compared to the northernmost region. Willow ptarmigan selected for herb-rich meadows and avoided lichen rich heathlands. There was generally stronger selection for vegetation types with dense field layer and for rich bogs and avoidance of vegetation types with sparse field layer cover and for lowland forest. Further, willow ptarmigan selected for areas around the timberline and for intermediate slopes. Mapping of the RSF showed that 60% of Norway is in the lowest ranked RSF bin and only 2% in the highest ranked RSF bin.


Willow ptarmigan selected for vegetation types with dense field layer and bogs at intermediate slopes around the timberline. Selection coincides with previous habitat selection studies on willow ptarmigan. This is the first attempt to assess and quantify habitat selection for willow ptarmigan at a large scale using data from line transect distance sampling surveys. Spatial variation in predictability suggests that habitat selection in late summer might vary from north to south. The resource selection map can be a useful tool when planning harvest quotas and habitat interventions in alpine areas.


Knowledge about patterns of habitat selection is often needed in order to make evidence-based management decisions. For instance, Smereka et al. [1] mapped the relative probability of selection for den sites for grizzly bears (Ursus arctos horibilis) in the Mackenzie Delta, Northwest Territories, Canada, to reduce human-bear conflicts by guiding human activity and land-use. A habitat selection model [2] has also been applied to harvest management of willow ptarmigan (Lagopus lagopus) in Northern Norway, where quotas are estimated based on a combination of pre-harvest densities and the amount of suitable habitat for willow ptarmigan available within the hunting area. Although a plethora of methods have been developed to assess patterns of habitat selection, resource selection functions (RSF) [3] are among the most frequently used methods to model habitat selection in animals (e.g., [1, 4,5,6]). An RSF is a function that is proportional to the probability of selection by an animal (c.f. [7, 8]) and is estimated directly from data. Data in the context of RSF’s is usually a set of locations where individuals are observed and a set of randomly generated available locations where the individuals could have been observed (used vs. available units). Variables associated with the observations may be habitat variables or covariates like elevation, topology, vegetation types or human-disturbance and infrastructure metrics. In a presence/available framework, variables that are assumed important for habitat selection are compared at the locations of the observations and the locations of available sites. A predictive RSF can be mapped and used to predict the relative selection for different geographical units based on their environmental characteristics [3, 7]. RSFs have been developed for many wildlife species [1, 2, 4, 5, 9, 10], often with the purpose to quantify suitable habitat for species that are of conservation concern or to answer questions related to the ecological dynamics of the system.

Willow ptarmigan is a medium-sized grouse species distributed in tundra habitats and boreal forests in the northern hemisphere [11]. Willow ptarmigan inhabits treeless alpine- and arctic tundra habitats most of the snow-free season, but it also occurs in northern boreal forests near tundra habitats. In Norway, it was recently listed as near threatened (NT) in the national red list of species [12], due to a 15–30% decline in the breeding population during the last decade (e.g., [13]). Predation, climate change, unsustainable harvest and potentially loss of habitat are all assumed to have contributed to the observed decline, but the relative contribution from these effects is to date not quantified [12].

There has been a rapid development of human infrastructure in the tundra areas of Norway during the nineteenth century, mainly due to construction of hydroelectric power installations, recreational facilities and roads. Moa et al. [14] found that autumn densities of willow ptarmigan were generally higher in management units with high proportion of areas located far from human infrastructure (see also [15, 16]). This suggests that habitat alteration affects willow ptarmigan populations. There is however a lack of studies quantifying and describing within range habitat selection for willow ptarmigan (but see [10, 17,18,19,20]). Some studies describing habitat selection however, have shown that willow ptarmigan broods, in general select for rich bogs close to willow-(Salix spp.) and dwarf birch (B. nana) thickets and avoid dry and sparsely vegetated habitats like heaths during summer and early autumn [10, 19, 20]. Willow ptarmigan is an indicator species for biodiversity in Norwegian mountains [21] and its distribution overlap with several other alpine and boreal species [13, 22]. Hence, it is possible that willow ptarmigan can play a role as an umbrella species in the alpine ecosystems, especially because of its status as an economically important game species [23].

In this study, we model willow ptarmigan habitat selection in August using exponential RSF’s. The exponential RSF was estimated using a mixed effects logistic regression model with a presence/available design. This study is the first objective attempt to model habitat selection for willow ptarmigan across Norway. We use satellite-based vegetation maps, slope and aspect from a digital elevation model and timberline measures as predictor variables. Based on previous studies of willow ptarmigan, we predict that willow ptarmigan selects vegetation types that offers food and concealment against predators. In addition, we predict that willow ptarmigan selects for areas close to the timberline, due to the special adaptation to alpine and sub-alpine areas during the snow-free season. We further expect that willow ptarmigan avoid steep slopes. Southernly exposed slopes have normally more vigorous plant communities than northerly slopes due to prolonged sun-exposure, so we also expect that willow ptarmigan select for southern slopes rather than northern slopes.


Study area

Data from survey areas distributed in alpine tundra, low arctic tundra and northern boreal forests throughout Norway (Fig. 1a) was used as a basis for our RSF. The vegetation in the tundra is dominated by small and medium-sized shrubs [e.g., willows, dwarf birch, and heath (Vaccinium spp. and Caluna spp.)]. The northern boreal forests are dominated by mountain birch (B. pubescens), Scots pine (Pinus sylvestris), and Norway spruce (Picea abies). The vertebrate fauna is dominated by large ungulates like wild- and semi-domestic reindeer (Rangifer tarandus) and moose (Alces alces), rodents (e.g., Microtus spp.), and terrestrial birds. Important predators on willow ptarmigan include red fox (Vulpes vulpes), stoat (Mustela erminea), pine marten (Martes martes), gyrfalcon (Falco rusticolus), rough-legged buzzard (Buteo lagopus) and golden eagle (Aquila chrysaetos). Livestock grazing by sheep and cattle is common in many areas during summer (Jun–Aug). Human population density is generally low within these areas. However, some areas are located in proximity to larger villages and in some areas, there are several cabin villages and isolated cabins, many of which are only seasonally inhabited.

Fig. 1
figure 1

Study area showing a the distribution of survey areas and vegetation types, b the relative probability of selection for willow ptarmigan during late summer/early autumn, ranging from low relative probability of selection (1) to high relative probability of selection (10)

Data collection

We used information about willow ptarmigan occurrence based on observations done along line-transect surveys from 2014 to 2017. Transects were in most areas spaced out systematically following map grids with a 500 m interval, and often either in north–south or east–west bearing (see Fig. 2a). The same transects were repeatedly surveyed each year. Following a distance sampling protocol [24] a dog handler with pointing dogs and an observer walked along the transect-lines with one free-ranging dog at a time searching both sides of the line [25, 26]. At each bird encounter, the observer recorded the species, total number of birds (juveniles, adult males, adult females, or birds of unknown sex or age), the perpendicular distance from the transect line to the observation, the geographical location (centre of the brood/covey if several birds) of the birds (in UTM) as well as the time of day the observation was made. The main purpose of the line-transect survey program is to estimate pre-harvest densities of willow ptarmigan in the survey areas for harvest management purposes. So, the survey areas are not chosen at random, but mainly driven by local initiatives and does not follow a strict design. However, the location-data collected during the surveys are also suitable for assessing habitat selection using resource selection functions [10]. Detailed description for the sampling protocol used for estimating willow ptarmigan densities is outlined in Pedersen et al. [27], Pedersen et al. [25] and Eriksen et al. [28].

Fig. 2
figure 2

Example area with transect-lines observations and random locations. To the left: Transect-lines placed in the landscape. On the right; a is a transect-line with all observations from 2014 to 2017, b the transect-line and observations with a 200-m buffer. Observations outside the 200-m buffer and observations with > 30-m deviation between reported and estimated perpendicular distance were excluded from analyses (see “Methods”). c Randomly generated available locations within the 200-m buffer, d available locations after adjusting for detection probability (see “Methods”) and e is willow ptarmigan locations and adjusted available locations used in the modelling

All observation data from the line transects surveyed were registered by the field worker to “Hønsefuglportalen” ( [29]), a common e-infrastructure and data portal supporting the line transect survey program in Norway. We used a JDBC-connection with Rstudio [30] and the library «RJDBC» [31] to download survey data from the SQLServer database. In total, 17,386 willow ptarmigan locations from 7923 surveys along 2543 distinct transect lines across 179 survey areas were downloaded. We established buffers covering 200 m on each side of the transect lines and we discarded all observations made outside the line-buffers (2778 used locations). The size of the buffers was set at 200 m because most of the used locations are within this distance from the transect-line (> 80%) and the procedure exclude observations that are either regarded as outliers or that have typo in the UTM coordinate text. Furthermore, based on the geographical position of the used locations and the geographical location of the transect line, we estimated perpendicular distances from transect lines to used locations and compared this to the perpendicular distance reported by the field workers. Whenever the deviation between these measures were > 30 m, we discarded the observation to optimize the quality of our data (cf. 30 * 30 m is the spatial resolution of the vegetation variables, see “Vegetation variables” below). Additional 2427 observations were omitted in this procedure (see Fig. 2 and Table 1 for details regarding data management).

Table 1 Summary of the process from importing data to the final dataset (see also “Methods” section)

Defining available habitat

To quantify the distribution of available habitat we generated random locations within the surveyed area. Based on a presence-available design (as in our case), observations are expected to be drawn from a sample of available locations [7]. Therefore, the randomly generated available locations should represent the available habitat within the area that is covered by the surveys. Since our observations are sampled along line transects it is expected by design that the detection probability decreases as a function of distance from the transect-line [10, 24, 27]. To achieve a proper distribution of available locations [7] we performed a four-step procedure following Kastdalen et al. [10] to achieve a similar probability distribution among observations and available locations. (1) Within the line-buffers, we generated random available locations (Ai). The number of locations per line-buffer was set according to the effort (km transect) multiplied by 3 (e.g. 4 km transect surveyed in 4 years: 1 available location * 16 km * 3 = 48 random locations within the line-buffer). Following this procedure, we also generated available locations for transect lines without used locations. This resulted in 66,174 random available locations (Ai). (2) We used the R library “distance” [32] to estimate half-normal detection functions (gxo) from the used locations (see above). (3) Then, gxo was used to estimate the detection probability (Oi) of each location Ai based on its distance from the transect line. Oi has a value between zero and one. (4) Finally, we generated a random number between zero and one (Pi) for each Ai. Ai was included in the set of available locations in the final analysis whenever Pi < Oi. After omitting available locations in step four, our sample of available locations were 38,149 distributed with the same probability distribution as the real observations [10] (Fig. 2b).

Predictor variables

Vegetation data

Vegetation type at willow ptarmigan locations and available locations was extracted from a digital raster-map with a resolution of 30 × 30 m (SatVeg [33]). This map consists of 25 generalized vegetation types [34] covering Norway (Fig. 1a and Table 2). Ten vegetation classes were classified as alpine, three were classified as bogs and open swamp vegetation, and eight as forest vegetation. The remaining four classes; water, agricultural areas, cities and built-up areas and unclassified or shadow affected areas were all considered as non-habitat. We kept alpine classes, bogs and swamps as original classes. We considered forest classes; Bilberry-(Vaccinium myrtillus), Crowberry-(Empetrum nigrum) and Lichen-rich birch forests as mainly sub-alpine mountain birch forest and kept as original classes’ while we pooled the remaining forest classes into one class representing lowland forests dominated by coniferous tree species.

Table 2 Categorical and continuous landscape variables used to determine relative probability of selection for willow ptarmigan

Landscape data

We extracted aspect in degrees and slope in degrees from a digital elevation model (DEM) with a resolution of 10 × 10 m (Norwegian Mapping Authority: We did not model the DEM ( directly because of the clear bioclimatic gradients present in our large study area.

Aspect is a circular variable (0°—north to 360°—north) and was therefore transformed to radians \(\left(r_{aspect} = aspect*\left( {\frac{2\pi }{360}} \right) \right)\), and in the next step we created two variables representing north-exposure \(\left(N_{aspect} = { \cos }\left( {r_{aspect} } \right) \right)\) and eastern exposure (\(\left(E_{aspect} = { \sin }\left( {r_{aspect} } \right) \right)\) [6]. We also constructed a categorical variable with five levels representing aspect, north (315°–45°), east (45°–135°), south (135°–225°), west (225°–315°) and flat areas (0°).

Timber line

To describe variation in altitude we used a raster-map with regional empirical timber line (RET) in meter above sea level, with a 100 × 100 m resolution [35]. To consider the bioclimatic variation caused by latitudinal and longitudinal gradients, we combined the DEM (see above) with the RET. Taking the regional empirical timber line at each location as the reference point (i.e. 0 m above the timberline) we calculated the deviation in altitude (meters) from the timberline to each willow ptarmigan location and available location.

All predictor variables at observations and available locations were extracted from the raster-maps using GRASS [36] and the function «r.what» through RStudio [30] with library «rgrass7» [37]. Summaries of the predictor variables are given in Table 2.

Statistical analysis

Model development

We omitted used and available locations placed in “non-habitat” (35 used locations and 336 available locations, se definition of non-habitat above). Our final data set used for analyses of willow ptarmigan habitat selection consisted of 12,146 willow ptarmigan observations and 38,149 available locations (Table 1). We estimated selection coefficients for the RSF by comparing environmental conditions of used locations to available locations using a mixed-effects logistic regression model [38,39,40]. Observations and available locations are stratified by survey areas. However, there is unbalance in the data since the number of locations per survey area vary. Gillies et al. [39] showed that using a random intercept when data are unbalanced improve model fit greatly and can change the direction of model coefficients. Therefore, we used Generalized Linear Mixed Model (GLMM) and fitted random intercepts for each survey area. The mixed effect model with random intercepts allow us further to account for spatial variation in density (density, because available points are generated based on effort not number of observed locations in the survey areas) and it provide marginal selection coefficients that can be used to make predictions also outside the sampled area. This approach assume that the explanatory variables has the same effects across all survey areas, although density may vary between areas. We did not consider effects of sex- and brood size on patterns of habitat selection.

The coefficients (\(\beta_{n} )\) estimated from the mixed-effect logistic regression model is the logarithmic relative selection strength (log-RSS) [41] for a given variable. We estimated the relative probability of selection at a given location using an exponential resource selection function (RSF):

$$w\left( x \right) = { \exp }\;\left( {\beta_{1} x_{1} + \beta_{2} x_{2} + \beta_{3} x_{3} + \beta_{4} x_{4} \ldots \beta_{n} x_{n} } \right),$$

where \(w\left( x \right)\) is relative probability of selection at location x and \(\beta_{1}\) thorough \(\beta_{n}\) is the estimated relative selection strength for explanatory variables \(x_{1}\) thorough \(x_{n}\) from the logistic regression model. Note that we do not include random intercepts in the RSF since or objective is to predict relative probability of selection across the whole of Norway.

To facilitate model convergence, we standardized continuous variables (aspect, slope, deviation from the timberline) to zero mean and one standard deviation [42].

We evaluated different combinations of explanatory variables, but always including the variable vegetation type when developing the RSF. We also tested for quadratic effects of slope and deviation from the timberline, and we evaluated models including aspect either as one categorical variable with five levels (north, east, south, west and flat) or as two continuous variables (north–south and east–west).

Model selection and validation

Model selection was based on AIC (Akaike Information Criterion). We considered 24 model combinations and the model with lowest AIC is considered the best supported model, but ΔAIC < 2 suggest that models are statistically equivalent and thus equally supported by the data [43]. In such situations, we followed the principle of parsimony and selected the least complex model. The most parsimonious model was evaluated using k-fold cross-validation [7, 44]. K-fold cross-validation yield metrices to assess a model’s ability to predict high relative probability of selection at locations where the species are observed. We divided the dataset into 5 approximately equal-sized datasets (folds) using random k-fold portioning in the R-library “dismo” [45]. For each model under validation, the k-fold cross-validation procedure followed four steps (repeated for each of the fivefold). (1) We withheld onefold (test set) and estimated model parameters based on the remaining fourfold (training set). (2) We used model parameters from the model in step 1) to predict the withheld test set. (3) A 10-quantile binning was generated on the predicted values of the test set. Bins were ranked from low relative probability of selection (bin #1) to high relative probability of selection (bin #10). A model with good predictive performance tends to have successively more willow ptarmigan locations in higher ranked bins. (4) A statistical metric of model performance was assessed by spearman rank-correlation between bin-rank and the count of used locations in each bin. Strong positive correlation coefficient suggests good predictive performance [7]. We also performed the four cross-validation steps described above on four geographical regions (blocks, c.f. [46]) instead of fivefold. The regions were drawn up by similarities in climatic conditions, and willow ptarmigan population dynamics within regions are synchronized [22].

Predicting resource selection functions across Norway

We developed a predictive resource selection map based on the most parsimonious model. In order to do so we had to convert all raster maps to 30 × 30 m resolution using program GRASS with the function «r.resamp.stats» [36], through RStudio [30] with library «rgrass7» [37]. Then we estimated the resource selection function for each 30 × 30 m cell by putting selection coefficients (\(\beta_{1}\), \(\beta_{n} \ldots\)) from the selected model and the raster values (i.e. explanatory variables \(x_{1}\), \(x_{n \ldots }\)) into Eq. (1). This procedure creates a new raster containing RSF values for each cell. The predicted RSF values where scaled so that they were bounded between zero and one by dividing by the maximum RSF value. Following the recommendations in Morris et al. [47], we mapped the RSF values based on the same quantile bins as we used in the k-fold cross-validation procedure. Hence the map classification ranges from category 1 (low relative probability of selection) to 10 (high relative probability of selection). We calculated the percentage distribution of each RSF-bin across the whole of Norway, within transect buffers and among the used willow ptarmigan locations. After producing the RSF-map for the whole of Norway, we performed another validation procedure using an independent willow ptarmigan location dataset from the Global Biodiversity Information Facility (gbif) [48]. We downloaded 5787 observations of willow ptarmigan that were recorded by ornithologists in Norway during June to September from 2000 to 2017 [48]. We extracted the RSF-bin rank for each observation using GRASS [36] and the function «r.what» through RStudio [30] with library «rgrass7» [37] and counted the number of observations in each RSF-bin. Some observations (1258) were in empty map-cells (missing data for one or more of the predictor maps used in the modelling or were in non-habitat cells) so a total of 4529 independent observations could be linked to an RSF-bin. We calculated spearman rank correlation between the RSF-bin rank and area-adjusted number of observations in each bin. Area-adjusted number of observations was the number of observations in each bin divided by the availability of that RSF-bin in Norway (c.f., Table 6).


Model selection and validation

In our model set, two models were equally supported by the data when considering the trade-off between model fit and number of parameters (i.e. ΔAIC < 2: Table 3). Both models included vegetation type (16 categories), deviation from the timberline (linear and quadratic terms) and slope (linear and quadratic terms). The highest ranked model included aspect as a categorical variable, whereas the second ranked model did not include any terms for aspect. Hence, the second ranked model had fewer parameters (21 vs. 25) and was considered as the most parsimonious model. K-fold cross-validation based on the most parsimonious model showed high and significant correlation for all fivefold (Table 4a), with a mean cross-correlation r = 0.96. This suggest that willow ptarmigan locations tend to be in high ranked habitats predicted by the model. All correlations from the regional cross-validation were significant, but the coefficients were much higher for the eastern, southern and western regions than for northern region (Table 4b).

Table 3 The 95% confidence set for models with ΔAIC < 5
Table 4 Spearman rank correlations between RSF bin ranks and count of used locations in each RSF bin for (a) fivefold portioning of the data and (b) spatial blocking of the data into four geographical regions

Selection coefficients

Willow ptarmigan generally selected for areas around the timberline (Table 5). More specifically, selection increased from lower elevations towards the timberline and decrease above a polynomial inflection point at 25 m above the timberline (Fig. 3). Willow ptarmigan also select for intermediate slopes (Table 5) whit an inflection point at 10°, whereas steeper slopes were increasingly avoided (Fig. 4).

Table 5 Parameter estimates from the most parsimonious model
Fig. 3
figure 3

Relative probability of selection w(x) as a function of deviation from the timberline and deviation from the timberline2. The “plus” marker shows the polynomial inflection point

Fig. 4
figure 4

Relative probability of selection w(x) as a function of slope and slope2 in degrees. The “plus” marker shows the polynomial inflection point

Among the 16 vegetation classes (Table 5), there was a clear selection for herb-rich meadows and clear avoidance for lichen-rich heathland relative to other habitat types. Four other vegetation types had parameter estimates above average; ombrotrophic bog, tall-grown lawn vegetation, fresh heather- and dwarf-shrub and crowberry birch forest. Other vegetation types with selection rates well below average included exposed ridges, dwarf willow snow patches, late snow patches, lichen-rich birch forest, and lowland forest classes.

Predicting resource selection functions across Norway

We predicted relative selection probability for willow ptarmigan across Norway by calculating an exponential RSF (formula 1) using selection coefficients (Table 5) estimated based on the most parsimonious model in the candidate set (Fig. 1b). This model predicts that almost 60% of the mapped area (excluding non-habitat, see methods) is located in the lowest RSF-bin (i.e. low habitat suitability). Only 2% is located in the highest ranked RSF-bin (i.e. highest habitat suitability). If we take the lowest RSF-bin to be unsuitable for willow ptarmigan and the remaining bins 2–10 are increasingly suitable, about 41% of Norway is suitable for Willow ptarmigan (when not considering non-habitat). Further, for the area covered by the line transect surveys (i.e., areas within 200 ms from the transect lines), only 12% is predicted to be unsuitable and about 88% suitable for willow ptarmigan. For percentage distribution of RSF bins, see Table 6. Spearman rank-correlation between RSF-bin rank and number willow ptarmigan observations from gbif adjusted for availability of bin-ranks in Norway (c.f., Table 6 column 1) was high and significant (rho = 0.857, p = 0.002). This suggest that the model also have high predictability on independent data.

Table 6 Percentage cover of RSF-bins 1–10; across the whole of Norway, within the surveyed area (i.e., within 200-m buffers) and the distribution of willow ptarmigan locations across RSF bins, respectively


The habitat suitability model developed here for willow ptarmigan was highly predictive according to both the k-fold cross-validation and validation against independent data, and therefore identify important habitats for willow ptarmigan during late summer/early autumn in Norway. Willow ptarmigan in general selected for herb-rich meadows and avoided most strongly lichen-rich heathland. Inspection of all vegetation type coefficients show that the general trend is selection for alpine vegetation types and bogs with abundant field layer and avoidance for sparsely vegetated alpine vegetation types and lowland forest vegetation types. In addition to selection and avoidance for different vegetation types, willow ptarmigan selected for areas around the timberline and for intermediate slopes.

Several other studies have described habitat selection in willow ptarmigan during different seasons [2, 10, 17,18,19,20, 49]. We describe habitat selection during late summer/early autumn. Kastdalen et al. [10] used similar data from August surveys, but for a much smaller area and using finer scaled vegetation data. In the study by Kastdalen et al. [10], willow ptarmigan selected rich bogs close to willow- and dwarf birch thickets and they avoided dry and poor open areas like heath. Our results coincide with Kastdalen et al. [10] with respect to both selection and avoidance. First, we found relatively strong selection for open alpine vegetation types with dense field layer and for bogs. Second, we found that willow ptarmigan avoided lichen-rich heath and other sparsely vegetated alpine areas such as ridges and late snow-patch vegetation. Similarly, studies of willow ptarmigan brood movements during summer in Northern Norway [20] and central Norway [19] also found that young broods used heaths less frequently than expected from the availability. In both studies, broods selected strongly for rich bogs. The vegetation types we used have a broad definition, so it is difficult to separate important small-scale habitat features. We can however, separate vegetation types that potentially contain important habitat features for willow ptarmigan (e.g. willow thickets, bogs and dwarf birch thickets [10, 17,18,19,20]). Willow thickets occur in bilberry-low fern birch forest and crowberry birch forest, on tall-grown sedge bogs and in fresh heather and dwarf-shrub communities. Among the bogs, both tall-grown sedge bogs and ombrotrophic bogs are high ranked vegetation types. Wet bogs and swamps is characterized by water level on the surface throughout the growing season and has an intermediate selection rank. Dwarf birch thickets occur in bilberry-low fern birch forest, fresh heather and dwarf-scrub, tall-grown sedge bogs, lichen-rich birch forest, heather- and grass-rich leeside. Dwarf birch do also occur on lichen-rich heathland and heather-rich alpine ridges, but in these sparsely vegetated and wind-exposed vegetation types, dwarf birch occurs more sparsely than in the latter types. The highest ranked vegetation type—herb-rich meadows—do not typically include willows or dwarf birch. It is however, characterized as the most nutrient rich vegetation type in the alpine region often with a stable water supply [34]. This might suggest that herb-rich meadows interact with other preferred vegetation types such as bogs or fresh heather and dwarf-shrub communities. Ehrich et al. [18], Henden et al. [50] advocate the importance of willow tickets for willow ptarmigan occupancy and both studies were carried out in low arctic tundra (e.g. same as the northern region in this study). Unfortunately, our vegetation maps cannot separate this vegetation structure, but it is possible that such strong selection for willow thickets in the arctic tundra is a special adaptation to a different environment.

Although the model is highly predictive according to the k-fold cross-validation on independent test data, the regional cross-validation for the northernmost region is barely significant. One reason for this low correlation could be the relatively low number of survey areas in this region. Only 32 survey areas out of a total of 179, and only 1096 used locations out of a total of 12,146 were in the northernmost region. However, this explanation might not be supported because the model was highly predictive in the western region where the observations are even fewer (14 survey areas and 426 observations). More likely, either willow ptarmigan selects differently in the low arctic tundra than in alpine tundra, or the broadly defined vegetation types in our vegetation map has different forms in low arctic tundra in north compared to areas further south. Another possible explanation can be that the relative probability of selection for a certain habitat type change with the availability of that habitat type (e.g., as a functional response [38, 51]). Hence, the availability of certain habitat types might vary from north to south. Since we use marginal selection coefficients applied to the whole of Norway, we are not able to account for such effects. Although the correlation is significant (Table 4b), building of a specific model for the northern region could result in a better predictive RSF for this region.

We followed the recommendation for mapping RSFs in Morris et al. [47] and partitioned our RSF predictions into 10 quantile bins before predicting the RSF across Norway. The bin ranks stretch from low relative probability of selection (1) to high relative probability of selection (10). The k-fold cross validation results suggest the model is highly predictive in describing the relative change in probability from bin 1 to bin 10. About 60% of Norway, 5% of the willow ptarmigan locations and 12% of the surveyed area (within buffers) are in lowest RSF bin (Table 6). Our data originate from a distance-sampling scheme, where the primary aim was to estimate willow ptarmigan densities in survey areas. Thus, data is systematically collected in areas where willow ptarmigan is expected to occur. When we predict across Norway, this will in turn lead to challenges in portion habitats with no probability of selection (e.g., forest and meadow habitats at low elevations, and coastal areas in southern latitudes and far from willow ptarmigan core habitats) since such areas were not surveyed. Based on this, we expect the lowest RSF bin to both include areas of low relative probability of selection that might be in or adjacent to alpine areas and no probability of selection (i.e. sites are far from willow ptarmigan core areas).

In general, animals are more abundant in habitats that are selected most strongly, and Boyce et al. [52] proposed that abundance can be estimated directly from habitat selection models for populations at the carrying capacity or for populations following an ideal free distribution. However, despite this and other studies documenting that willow ptarmigan select for specific habitat features [2, 10, 17,18,19, 49], Kvasnes et al. [53] found no clear relationship between willow ptarmigan density estimated pre-harvest (in August) and proportion of different habitat categories within survey areas. As also noted by Boyce et al. [52], Kvasnes et al. [53] suggested that other factors that are not directly related to habitat also influences the abundance. Willow ptarmigan population densities in the study of Kvasnes et al. [53] were generally lower than historic densities [54, 55]. In addition, different harvest strategies [56] and varying predation rates can have great influence on population densities in willow ptarmigan, both of which can vary independent of habitat composition. Habitat selection can also be affected by social interactions such as conspecific attraction [57], which is also suggested as a possible factor influencing the distribution of willow ptarmigan [58]. In our model, we assume that the explanatory variables have the same effect across all survey areas. Thus, we predict relative probability of selection for willow ptarmigan regardless of density and factors affecting density such as harvest rates, predator densities and human infrastructure.

Our RSF is based on data collected in August, so the resource selection map we developed is restricted to a short period of time. It is the phase when chicks are fledged, still accompanied by the adults, but only a few weeks prior to harvest. About 120,000–200,000 willow ptarmigans are shot each year in Norway (Statistics Norway,, so an RSF based on data from August combined with density estimates from the same period can be a useful tool for harvest management. Further, the high and significant correlation for the independent validation-dataset recorded from June to September suggests that the RSF might be a useful tool to identify willow ptarmigan habitat potential during most of the snow-free season. An RSF based on winter locations of willow ptarmigan could also be useful in land-use management since most of the new infrastructure development in alpine and tundra areas is building of recreational facilities such as cabin villages. Cabin villages is mainly located in sub-alpine birch forests and boreal forests [59] and mountain birch is one of the main food items for willow ptarmigan during winter [60].


We found that willow ptarmigan selects for vegetation types with dense field layer and bogs at intermediate slopes around the timberline. Selection for vegetation types coincide with previous small-scale habitat selection studies on willow ptarmigan. This is the first attempt to assess and quantify habitat selection for willow ptarmigan at a large scale using data from line transect distance sampling surveys, and model validation show that the model is highly predictive. Based on our RSF it is predicted that c. 60% of the land area in Norway is in the lowest RSF bin. This implies that these areas have very low probability of selection by willow ptarmigan compared to higher ranked RSF bins. Spatial variation in the model’s predictive capability suggest that habitat selection vary from north to south. Despite the conservation status near threatened (NT) in Norway, willow ptarmigan is a highly valued game species, and c. 50,000 register to hunt ptarmigan each year (Statistics Norway, There is also a continued demand for hydroelectric power installations, recreational facilities and roads within willow ptarmigan areas. The consequences are loss of habitat and thus a reduced potential carrying capacity for willow ptarmigan and other alpine species. The resource selection map we have developed predicts relative probability of selection during summer/early autumn and can be a useful management tool in the conservation of habitats across Norway. We recommend that the RSF is considered when planning harvest quotas and when planning habitat interventions in willow ptarmigan ranges. The latter can potentially also benefit other species with overlapping distributions.



resource selection function is a function that is proportional to the probability of selection by an animal


digital elevation model, a raster containing pixels with values representing elevation


regional empirical timberline, a raster containing pixels with values representing the location of the timberline in meters above sea level


Global Biodiversity Information Facility


  1. Smereka CA, Edwards MA, Pongracz J, Branigan M, Pilfold NW, Derocher AE. Den selection by barren-ground grizzly bears, Mackenzie Delta, Northwest Territories. Polar Biol. 2017;40(3):503–16.

    Article  Google Scholar 

  2. Pedersen ÅØ, Jepsen JU, Biuw M, Johansen B. Habitatmodell for lirype i Finnmark. NINA Rapport. 2012;845:36.

    Google Scholar 

  3. Manly BFJ, McDonald LL, Thomas DL, McDonald TL, Erikcson WP. Resource selection by animals: statistical design and analysis for field studies. Dordrecht: Kluwer Academic; 2002.

    Google Scholar 

  4. Carpenter J, Aldridge C, Boyce MS. Sage-grouse habitat selection during winter in Alberta. J Wildlife Manag. 2010;74(8):1806–14.

    Article  Google Scholar 

  5. Barker OE, Derocher AE. Habitat selection by arctic ground squirrels (Spermophilus parryii). J Mammal. 2010;91(5):1251–60.

    Article  Google Scholar 

  6. Pedersen AO, Blanchet MA, Hornell-Willebrand M, Jepsen JU, Biuw M, Fuglei E. Rock Ptarmigan (Lagopus muta) breeding habitat use in northern Sweden. J Ornithol. 2014;155(1):195–209.

    Article  Google Scholar 

  7. Boyce MS, Vernier PR, Nielsen SE, Schmiegelow FKA. Evaluating resource selection functions. Ecol Model. 2002;157(2–3):281–300.

    Article  Google Scholar 

  8. Lele SR, Merrill EH, Keim J, Boyce MS. Selection, use, choice and occupancy: clarifying concepts in resource selection studies. J Anim Ecol. 2013;82(6):1183–91.

    Article  Google Scholar 

  9. Swanson CC, Rumble MA, Grovenburg TW, Kaczor NW, Klaver RW, Herman-Brunson KM, Jenks JA, Jensen KC. Greater sage-grouse winter habitat use on the eastern edge of their range. J Wildlife Manag. 2013;77(3):486–94.

    Article  Google Scholar 

  10. Kastdalen L, Pedersen HC, Fjone G, Andreassen HP. Combining resource selection functions and distance sampling: an example with willow ptarmigan. In: Huzurbazar S, editor. Resource selection methods and application. Cheyenne: Western EcoSystems Technology; 2003. p. 52–9.

    Google Scholar 

  11. Johnsgard PA. The grouse of the world. Kent: Croom Helm Ltd; 1983.

    Google Scholar 

  12. Henriksen S, Hilmo O. Norsk rødliste for arter 2015. Norge: Artsdatabanken; 2015.

    Google Scholar 

  13. Lehikoinen A, Green M, Husby M, Kalas JA, Lindstrom A. Common montane birds are declining in northern Europe. J Avian Biol. 2014;45(1):3–14.

    Article  Google Scholar 

  14. Moa PF, Nilsen EB, Hagen BR, Sund T. Natural ptarmigan refuges; do they exist and do they result in higher ptarmigan densities? Utmark. 2013: 2b.

  15. Broseth H, Pedersen HC. Hunting effort and game vulnerability studies on a small scale: a new technique combining radio-telemetry, GPS and GIS. J Appl Ecol. 2000;37(1):182–90.

    Article  Google Scholar 

  16. Støen O-G, Wegge P, Heid S, Hjeljord O, Nellemann C. The effect of recreational homes on willow ptarmigan (Lagopus lagopus) in a mountain area of Norway. Eur J Wildl Res. 2010;56(5):789–95.

    Article  Google Scholar 

  17. Henden JA, Ims RA, Yoccoz NG, Killengreen ST. Declining willow ptarmigan populations: the role of habitat structure and community dynamics. Basic Appl Ecol. 2011;12(5):413–22.

    Article  Google Scholar 

  18. Ehrich D, Henden JA, Ims RA, Doronina LO, Killengren ST, Lecomte N, Pokrovsky IG, Skogstad G, Sokolov AA, Sokolov VA, et al. The importance of willow thickets for ptarmigan and hares in shrub tundra: the more the better? Oecologia. 2012;168(1):141–51.

    Article  Google Scholar 

  19. Andersen R, Steen JB, Pedersen HC. Habitat selection in relation to the age of willow grouse Lagopus l. lagopus broods in Central Norway. Fauna norv Ser C, Cinclus. 1984;7:90–4.

    Google Scholar 

  20. Erikstad KE. Growth and survival of willow grouse chicks in relation to home range size, brood movements and habitat selection. Ornis Scand. 1985;16(3):181–90.

    Article  Google Scholar 

  21. Pedersen HC, Eide N. Fjell. In: Nybø S, editor. Naturindeks for Norge, Utredning for DN, vol. 3. Trondheim: Direktoratet for naturforvaltning; 2010. p. 109–23.

    Google Scholar 

  22. Kvasnes MAJ, Storaas T, Pedersen HC, Bjork S, Nilsen EB. Spatial dynamics of Norwegian tetraonid populations. Ecol Res. 2010;25(2):367–74.

    Article  Google Scholar 

  23. Wam HK, Andersen O, Kaltenborn BP. Markedet for rypejakt. In: Pedersen HC, Storaas T, editors. Rypeforvaltning, Rypeforvaltningsprosjektet 2006–2011 og veien videre (In Norwegian). Oslo: Cappelen Damm Akademisk; 2013.

    Google Scholar 

  24. Buckland ST, Anderson DR, Burnham KP, Laake JL, Borchers DL, Thomas L. Introduction to distance sampling: estimating abundance of biological populations. New York: Oxford Univerity Press, Inc.; 2001.

    Google Scholar 

  25. Pedersen HC, Steen H, Kastdalen L, Brøseth H, Ims RA, Svendsen W, Yoccoz NG. Weak compensation of harvest despite strong density-dependent growth in willow ptarmigan. Proc R Soc Lond B-Biol Sci. 2004;271(1537):381–5.

    Article  CAS  Google Scholar 

  26. Warren P, Baines D. Evaluation of the distance sampling technique to survey red grouse Lagopus lagopus scoticus on moors in northern England. Wildlife Biol. 2011;17(2):135–42.

    Article  Google Scholar 

  27. Pedersen HC, Steen H, Kastdalen L, Svendsen W, Brøseth H. Betydningen av jakt på lirypebestander: framdriftsrapport 1996–1998. NINA Oppdragsmelding. 1999;578:43.

    Google Scholar 

  28. Eriksen LF, Moa PF, Nilsen EB. Quantifying risk of overharvest when implementation is uncertain. J Appl Ecol. 2017;55:482.

    Article  Google Scholar 

  29. Nilsen EB. Pedersen HC. Vang R: Hønsefuglportalen—en nasjonal portal for ryper og skogsfugl. NINA Minirapport; 2013. p. 423.

    Google Scholar 

  30. R-Core-Team. R: a language and environment for statistical computing. In: R Foundation for Statistical Computing, Vienna, Austria; 2016. ISBN 3-900051-07-0.

  31. Urbanek S. RJDBC: provides access to databases through the JDBC interface (2011). R package version 0.2–5; 2014.

  32. Lawrence DM: Distance: Distance sampling detection function and abundance estimation, R package version 0.9.6. 2016.

  33. Johansen BE. Vegetasjonskart for Norge basert på landsat TM/ETM + data (in Norwegian). In: Rapport 4/2009. Tromsø: Norut Tromsø; 2009: 87s.

  34. Johansen BE, Aarrestad PA, Øien DI. Vegetasjonskart for Norge basert på satellittdata, Delprosjekt 1: klasseinndeling og beskrivelse av utskilte vegetasjonstyper (in Norwegian). Tromsø: Norut Tromsø; 2009. p. 34.

    Google Scholar 

  35. Blumenrath S, Hanssen F. Beregning av areal. In: Nybø S, ed. Datagrunnlag for Naturindeks 2010, DN-utredning 4-2010; 2010. p. 8–19.

  36. GRASS-Development-Team. Geographic resources analysis support system (GRASS) software. Version 7.2 edn. Open Source Geospatial Foundation; 2015.

  37. Bivand R. rgrass7: interface between GRASS 7 geographical information system and R. R package version 0.1–9. 2016.

  38. McLoughlin PD, Morris DW, Fortin D, Vander Wal E, Contasti AL. Considering ecological dynamics in resource selection functions. J Anim Ecol. 2010;79(1):4–12.

    Article  Google Scholar 

  39. Gillies CS, Hebblewhite M, Nielsen SE, Krawchuk MA, Aldridge CL, Frair JL, Saher DJ, Stevens CE, Jerde CL. Application of random effects to the study of resource selection by animals. J Anim Ecol. 2006;75(4):887–98.

    Article  Google Scholar 

  40. Hebblewhite M, Merrill E. Modelling wildlife-human relationships for social species with mixed-effects resource selection models. J Appl Ecol. 2008;45(3):834–44.

    Article  Google Scholar 

  41. Avgar T, Lele SR, Keim JL, Boyce MS. Relative selection strength: quantifying effect size in habitat- and step-selection inference. Ecol Evol. 2017;7(14):5322–30.

    Article  Google Scholar 

  42. Schielzeth H. Simple means to improve the interpretability of regression coefficients. Methods Ecol Evol. 2010;1(2):103–13.

    Article  Google Scholar 

  43. Burnham KP, Anderson DR. Model selection and multimodel inference: a practical information-theoretic approach. New York: Springer; 2002.

    Google Scholar 

  44. Nielsen SE, Boyce MS, Stenhouse GB. Grizzly bears and forestry I. Selection of clearcuts by grizzly bears in west-central Alberta, Canada. For Ecol Manage. 2004;199(1):51–65.

    Article  Google Scholar 

  45. Hijmans R, Phillips S, Leathwick J, Elith J. R package dismo: species distribution modeling, version 1.1–1. 2016.

  46. Wiens TS, Dale BC, Boyce MS, Kershaw GP. Three way k-fold cross-validation of resource selection functions. Ecol Model. 2008;212(3–4):244–55.

    Article  Google Scholar 

  47. Morris LR, Proffitt KM, Blackburn JK. Mapping resource selection functions in wildlife studies: concerns and recommendations. Appl Geogr. 2016;76:173–83.

    Article  Google Scholar 

  48. Hoem SA. Norwegian Species Observation Service. Version 3.15. The Norwegian Biodiversity Information Centre (NBIC). Occurrence dataset Accessed via on 2018-07-14. 2018.

  49. Lande US. Grouse—habitat relationships: monitoring, scale and management. Doctoral thesis. Uppsala: Sveriges lantbruksuniversitet; 2011.

  50. Henden J-A, Ims RA, Fuglei E, Pedersen ÅØ. Changed arctic-alpine food web interactions under rapid climate warming: implication for ptarmigan research. Wildlife Biol. 2017;2017:wlb00240.

    Article  Google Scholar 

  51. Mysterud A, Ims RA. Functional responses in habitat use: availability influences relative use in trade-off situations. Ecology. 1998;79(4):1435–41.

    Article  Google Scholar 

  52. Boyce MS, Johnson CJ, Merrill EH, Nielsen SE, Solberg EJ, van Moorter B. Can habitat selection predict abundance? J Anim Ecol. 2016;85(1):11–20.

    Article  Google Scholar 

  53. Kvasnes MAJ, Pedersen HC, Storaas T, Nilsen EB. Vegetation type and demography of low density willow ptarmigan populations. J Wildlife Manag. 2017;81(1):174–81.

    Article  Google Scholar 

  54. Myrberget S. Demography of an island population of willow ptarmigan in northern Norway. In: Bergerud AT, Gratson MW, editors. Adaptive strategies and population Ecology of northern grouse. Minneapolis: University of Minnesota Press; 1988. p. 379–419.

    Google Scholar 

  55. Pedersen HC. Territorial behavior and breeding numbers in Norwegian willow ptarmigan—a removal experiment. Ornis Scand. 1988;19(2):81–7.

    Article  Google Scholar 

  56. Sandercock BK, Nilsen EB, Brøseth H, Pedersen HC. Is hunting mortality additive or compensatory to natural mortality? Effects of experimental harvest on the survival and cause-specific mortality of willow ptarmigan. J Anim Ecol. 2011;80(1):244–58.

    Article  Google Scholar 

  57. Campomizzi AJ, Butcher JA, Farrell SL, Snelgrove AG, Collier BA, Gutzwiller KJ, Morrison ML, Wilkins RN. Conspecific attraction is a missing component in wildlife habitat modeling. J Wildlife Manag. 2008;72(1):331–6.

    Article  Google Scholar 

  58. Kvasnes MAJ, Pedersen HC, Solvang H, Storaas T, Nilsen EB. Spatial distribution and settlement strategies in willow ptarmigan. Popul Ecol. 2015;57(1):151–61.

    Article  Google Scholar 

  59. Taugbøl T, Vistad OI, Nellemann C, Kaltenborn BP, Flyen AC, Swensen G, Nybakken A, Horgen BC, Grefsrud R, Lein K, et al. Hyttebygging i Norge: en oppsummering og vurdering av ulike miljø- og samfunnsmessige effekter av hyttebygging i fjell-og skogtraktene i Sør-Norge, vol. 709. Trondheim: Norsk institutt for naturforskning; 2001.

    Google Scholar 

  60. Mossop DH. Winter survival and spring breeding strategies of willow ptarmigan. In: Bergerud AT, Gratson MW, editors. Adaptive strategies and population ecology of northern grouse. Minneapolis: University of Minnesota Press; 1988. p. 330–77.

    Google Scholar 

Download references

Authors’ contributions

MAJK and EBN conceived and designed the study. EBN and HCP applied for and received funding for the project. MAJK carried out the analysis and main writing of the manuscript. EBN and HCP contributed to discussions and interpretation of the analysis and writing of the manuscript. All authors read and approved the final manuscript.


We are grateful to all the volunteer dog handlers who collected all the data in field. We are also grateful to Jenny Mattisson and Stefan Blumenrath for collaboration and help with GRASS GIS. Comments from Mark Boyce and one anonymous reviewer greatly improved the manuscript.

Competing interests

The authors declare that they have no competing interests.

Availability of data and materials

Datasets from line transect surveys on the land managed and owned by Statskog and Finnmarkseiendommen (FEFO) are available in the gbif repository; Statskog:, FEFO: Additional datasets used and analysed during the current study are available from the corresponding author on reasonable request and with permission of each landowner.

Consent to publish

Not applicable.

Ethics approval and consent to participate

Not applicable.


The project is founded by Statskog (Norwegian State-Owned Land and Forest Enterprise), Norwegian Environment Agency and Norwegian Institute for Nature Research.

Publisher’s Note

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

Author information

Authors and Affiliations


Corresponding author

Correspondence to Mikkel Andreas Jørnsøn Kvasnes.

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (, which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Kvasnes, M.A.J., Pedersen, H.C. & Nilsen, E.B. Quantifying suitable late summer brood habitats for willow ptarmigan in Norway. BMC Ecol 18, 41 (2018).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: