Hostname: page-component-7c8c6479df-ph5wq Total loading time: 0 Render date: 2024-03-29T07:32:02.030Z Has data issue: false hasContentIssue false

Coarse landscape features predict occurrence, but habitat selection is driven by specific habitat traits: implications for the conservation of the threatened Woodchat Shrike Lanius senator

Published online by Cambridge University Press:  30 May 2016

MATTIA BRAMBILLA*
Affiliation:
Lega Italiana Protezione Uccelli (LIPU/BirdLife Italia), Via Trento 49, I-43100 Parma, Italy. Fondazione Lombardia per l’Ambiente, Settore Biodiversità e Aree protette, Largo 10 luglio 1976 1, I-20822 Seveso, MB, Italy. Museo delle Scienze, Sezione di Zoologia dei Vertebrati, Via Calepina 14, I-38122 Trento, Italy.
MARCO GUSTIN
Affiliation:
Lega Italiana Protezione Uccelli (LIPU/BirdLife Italia), Via Trento 49, I-43100 Parma, Italy.
EGIDIO FULCO
Affiliation:
Studio Naturalistico Milvus, C/da Pantano 134, I-85010 Pignola PZ, Italy.
ALBERTO SORACE
Affiliation:
Lega Italiana Protezione Uccelli (LIPU/BirdLife Italia), Via Trento 49, I-43100 Parma, Italy. SROPU, Via R. Crippa 60, I-00125, Roma, Italy.
CLAUDIO CELADA
Affiliation:
Lega Italiana Protezione Uccelli (LIPU/BirdLife Italia), Via Trento 49, I-43100 Parma, Italy.
*
*Author for correspondence; e-mail: brambilla.mattia@gmail.com
Rights & Permissions [Opens in a new window]

Summary

Habitat selection has fundamental implications for species conservation, and in birds is often regarded as a multi-scale process. We investigated (under an information-theoretic approach) habitat selection by Woodchat Shrike Lanius senator in Italy (one of the most severely declining species in central and western Europe), considering five main types of potential determinants of shrike occurrence at the territory scale (1 ha): general structure (coarse landscape), woody vegetation, grassland habitats/bare ground, herbaceous crops, and management variables. The best supported models for species occurrence were those including general structure and woody vegetation traits. Variation partitioning suggested that overall, landscape general structure and woody vegetation explained the highest variation in shrike occurrence, and management the lowest. However, considering variation explained by only a single level, all levels performed nearly equally, but general structure did not explain an exclusive proportion of variation. A multi-level analysis suggested that shrike occurrence was eventually associated with specific habitat traits: isolated trees, shrubland and (secondarily) olive groves (all with positive effects), and dirt roads (negative effect). The most parsimonious multi-level models included only variables from woody vegetation and management traits, suggesting that the likely true determinants of species occurrence are highly specific and fine-scale habitat traits, consistent with variation partitioning. Woodchat Shrikes inhabit semi-open landscapes, within which they are attracted to shrubland and isolated trees (secondarily to olive groves) and avoid dirt roads. Suitable habitat conditions for the species depend on a trade-off between abandonment and intensive farming, and rural development programmes may be crucial for the conservation (or loss) of such conditions.

Type
Research Article
Copyright
Copyright © BirdLife International 2016 

Introduction

Habitat selection is a key process with fundamental implications for species conservation (Cody Reference Cody1985, Jones Reference Jones2001). It is defined as the process an organism uses to choose its habitat, which results in habitat preferences consisting of differential use of specific resources relative to their availability (Hall et al. Reference Hall, Krausman and Morrison1997). The choice of a habitat by a species has often been regarded as the outcome of a process involving multiple spatial scales (e.g. Johnson Reference Johnson1980, Orians and Wittenberger Reference Orians and Wittenburger1991, Jones Reference Jones2001, Brambilla et al. Reference Brambilla, Bassi, Ceci and Rubolini2010a). In several bird species, habitat selection seems to occur first at coarser scales, and then at finer scales, according to a hierarchical process (Johnson Reference Johnson1980, Jones Reference Jones2001, Battin and Lawler Reference Battin and Lawler2006, Brambilla et al. Reference Brambilla, Rubolini and Guidali2006). Therefore, multiple scales representing ‘coarse’ and ‘fine’ habitat variables are often considered in habitat selection studies, especially for avian species.

Several bird species are threatened by unfavourable changes to their breeding habitat at different levels (Tucker and Evans Reference Tucker and Evans1997), and a particularly alarming case is represented by farmland birds (Fuller et al. Reference Fuller, Gregory, Gibbons, Marchant, Wilson and Baillie1995, Siriwardena et al. Reference Siriwardena, Baillie, Buckland, Fewster, Marchant and Wilson1998, Krebs et al. Reference Krebs, Wilson, Bradbury and Siriwardena1999), which are dramatically and widely declining largely because of agricultural intensification (Tucker and Evans Reference Tucker and Evans1997, Chamberlain et al. Reference Chamberlain, Fuller, Bunce, Duckworth and Shrubb2000, Donald et al. Reference Donald, Green and Heath2001, Reference Donald, Sanderson, Burfield and Van Bommel2006), especially through loss of ecological heterogeneity (Benton et al. Reference Benton, Vickery and Wilson2003, Vickery and Arlettaz Reference Vickery, Arlettaz and Fuller2012), and land abandonment (e.g. Brambilla et al. Reference Brambilla, Casale, Bergero, Bogliani, Crovetto, Falco, Roati and Negri2010b). Both intensification and abandonment may affect the habitat of a species at different levels, from landscape structure (Suarez-Seoane et al. Reference Suarez-Seone, Orborne and Baudry2002, Benton et al. Reference Benton, Vickery and Wilson2003, Brambilla et al. Reference Brambilla, Casale, Bergero, Bogliani, Crovetto, Falco, Roati and Negri2010b) to fine-scaled vegetation traits (e.g. Vickery and Arlettaz Reference Vickery, Arlettaz and Fuller2012).

Understanding the factors affecting habitat selection and the scale at which they act is thus necessary to promote species conservation and particularly urgent for threatened farmland birds. Among them, Woodchat Shrike Lanius senator is probably one of the most severely declining species in a large part of Europe, showing a continuous decrease over several decades in central and western Europe, whereas recent trends for south-eastern Europe are less negative (BirdLife International 2015). This species was formerly distributed in the whole Mediterranean region and in central Europe, but has undergone a large population decline and range contraction and now occurs mostly in southern Europe, where it is still declining in the western countries (BirdLife International 2015). The decline of the species in Italy is particularly alarming: the population underwent a 70–80% reduction in 2000–2012 (Nardelli et al. Reference Nardelli, Andreotti, Bianchi, Brambilla, Brecciaroli, Celada, Dupré, Gustin, Longoni, Pirrello, Spina, Volponi and Serra2015) and is currently classified as ‘Endangered’ (Peronace et al. Reference Peronace, Cecere, Gustin and Rondinini2012), with a ‘bad’ conservation status (Brambilla et al. Reference Brambilla, Gustin and Celada2013a). Similar declines have been reported in the recent past for other countries, such as Spain, France, Switzerland, Germany, Poland, Croatia and Greece (BirdLife International 2015).

The Woodchat Shrike is reportedly associated with different natural and anthropogenic landscape traits (Table 1). On the basis of previous knowledge, we identified some potential determinants of species occurrence: i) general habitat structure, i.e. coarse landscape characteristics, ii) woody vegetation (trees, shrubs), iii) grassland habitats and bare ground, iv) cultivated (herbaceous) crops, and v) human-related variables (grazing management, roads, fences). The general habitat structure included environmental variables allowing for a coarse description of land cover and topography in the cells: these kinds of variables are comparable to the ones which can be extracted by commonly available GIS layers. The other levels included more detailed descriptors of the species, which should be generally recorded in the field (as in our study case), or obtained by means of more sophisticated approaches.

Table 1. Factors affecting Woodchat Shrike occurrence and habitat selection according to the available literature.

With this work, we aim to identify the habitat factors affecting species occurrence at the territory level and to evaluate the relative importance of different categories of habitat factors, corresponding to the five groups outlined above, which have been already reported as potentially important for the species (Cramp and Perrins Reference Cramp and Perrins1993, Shochat et al. Reference Shochat, Abramsky and Pinshow2002, Filippi-Codaccioni et al. Reference Filippi-Codaccioni, Devictor, Bas and Julliard2010). Those individual levels represent different kinds of environmental factors which can potentially affect the occurrence of Woodchat Shrike at the territory level. Evaluating their relative importance has essential implications for conservation, as the maintenance or restoration of suitable conditions for the species should be pursued by means of different strategies (e.g. landscape planning vs. agri-environmental schemes), according to the types of factors driving shrike occurrence.

Methods

Study areas and fieldwork

Woodchat Shrikes were censused in two different study areas (Figure 1): Tolfa (Central Italy, Rome province, Lazio region) and hilly areas of Matera province (southern Italy, Basilicata region). The two areas were selected as representative of extensive farming landscapes of the Mediterranean region, i.e. the most important macro-habitat of the species in Europe. Within the two areas, we defined seven and 12 sample plots, respectively, each one covering some tens of hectares. These plots include open and semi-open landscapes in areas with climate suitable for Woodchat Shrike. To each area, four visits were made by the same observers (A. S. in Tolfa, E. F. in Basilicata), in April–June 2011. The observers noted all shrike contacts on maps (1:2000), recording all territorial and breeding behaviours, such as carrying food for chicks, members of a pair seen together, singing males, aggressive behaviour, calling of juveniles, nest alarm, occupied nests. Pair territories were defined on the basis of all contacts with the species and were distinguished among each other mostly on the basis of simultaneous observations of different pairs or singing males, as usually done with passerine birds (e.g. Birrer et al. Reference Birrer, Spiess, Herzog, Jenny, Kohli and Lugrin2007, Ceresa et al. Reference Ceresa, Bogliani, Pedrini and Brambilla2012, Brambilla et al. Reference Brambilla, Fulco, Gustin and Celada2013b).

Figure 1. Location of study areas in Italy. Each study area included 7–12 plots within which fieldwork was carried out.

A grid of 100 m × 100 m cells (1 ha-cells) was then superimposed on each study plot. The specific cell size was established on the basis of the territory size of the species reported in the literature (Harris and Franklin Reference Harris and Franklin2000, Lefranc and Worfolk Reference Lefranc and Worfolk1997), often being fairly small (no more than 1 ha) (Cramp and Perrins Reference Cramp and Perrins1993). A cell was defined as occupied when it included one territory of the species (in one case, two territories occurred within the same cell). When needed, the exact location of the cells including territories was manually adjusted to better match the territory extent. Grid cells were then used as sample units for territories and control plots, and all habitat variables recorded referred to the 1-ha cells. Within each one of the 19 sample plots (see above), unoccupied ‘control’ cells were randomly selected in the same number of occupied ones. This led to an average number of 1.7 occupied and 1.7 control cells within each sample plot; this balanced design prevented clustering of territories within the two areas and the associated potential spatial biases.

Habitat variables

We recorded directly in the field some habitat variables describing the following habitat characteristics: i) the general structure of the habitat, ii) the specific features of the woody vegetation, i.e. trees, shrubs and permanent (woody) crops, iii) the type of herbaceous layer and the occurrence of rocky or bare surfaces, iv) the features of cultivated (herbaceous) crops, and v) variables describing human management and impact, such as road and fence length and occurrence of domestic grazing animals (Table 2). Habitat variables were recorded in all the selected cells (occupied and unoccupied; see above).

Table 2. Habitat variables considered in this study to evaluate habitat selection according to five different levels of habitat traits.

Statistical analyses

Land-cover variables were arcsine square-root transformed before analyses. In each subset, VIF (Variance Inflation Factor) values were lower than 2.2 in all cases; in the multi-scale final model (see below), VIF values were lower than 1.6.

To qualitatively describe the habitats occupied by the species, we performed a comparison of habitat features between occupied (n = 33) and unoccupied cells (n = 33), evaluating differences by means of a t-test or a χ2-test (the latter adopted for grazing occurrence; Table 3).

Table 3. Average features of shrike territories and control plots; * indicates significant (P < 0.05) differences (assessed by means of a t-test on arc-sin square-root transformed variables for land cover and road length, and by χ2-test for grazing occurrence).

We then built GLM models with territory occurrence as the dependent variables, by relating it to the different habitat variables. We adopted an information-theoretic approach (Burnham and Anderson Reference Burnham and Anderson2002), performing a two-step analysis. As a first step, to evaluate the relative importance of each group of variables and of individual factors within each group, all possible models for each group were ranked using the Akaike Information Criterion corrected for small sample sizes (AICc). We checked the potential occurrence of quadratic relationships by entering the squared term of each variable, and then retained in the set of variables that were entered in the model the quadratic terms that showed a negative effect coupled with a positive effect of the linear term. As a second step, from each of the five different types of habitat traits, we selected the variables included in the most parsimonious models (models with Δ AICc < 2) for each group, with the exception of the ‘uninformative parameters’ (cf. Arnold Reference Arnold2010). The latter are variables included only in models that comprised more supported and simpler models as nested ones (Ficetola et al. Reference Ficetola, Marziali, Rossaro, De Bernardi and Padoa-Schioppa2011); AICc used as the unique criterion for model selection may indeed over-select complex models (Richards Reference Richards2005). With the resulting set of variables, we worked out a single, multi-scale model. Then, all possible models were ranked according to AICc, and an average model was obtained by averaging the most supported models (models with Δ AICc < 2). Model ranking according to AICc and model averaging was done using the package ‘MuMIn’ (Bartoń Reference Bartoń2014) in R (R Development Core Team 2013).

Finally we performed a variation partitioning analysis to compare the contribution of variables measured at different scales in affecting habitat selection by Woodchat Shrike. This analysis partitions the variation in habitat selection into components associated with different levels. To reduce the number of levels (maximum number allowed for the analysis is four), we summarised our levels as follows: i) general structure, ii) woody vegetation, iii) grassland, bare areas and cereals, iv) management and anthropogenic traits. The fractions of variation were calculated from the adjusted R2, which allows an unbiased estimation of the portions of the variation explained by single levels and by their combination (Peres-Neto et al. Reference Peres-Neto, Legendre, Dray and Borcard2006). This analysis was carried out by means of the “vegan” package (Oksanen et al. Reference Oksanen, Blanchet, Kindt, Legendre, Minchin, O’Hara, Simpson, Solymos, Stevens and Wagner2013) in R.

Results

Woodchat Shrikes (34 territories occurring within 33 cells) occupied cells characterised by a prevailing cover of grassland, which on average occurred over around half of the cell, and with a significantly higher availability of isolated trees and shrubland than unoccupied cells (Table 3).

Among the sets of candidate models reflecting different potential determinants of habitat selection in shrikes, the ones including the most supported models were general structure and woody vegetation (Table 4).

Table 4. Candidate models reflecting different potential determinants of habitat selection in Woodchat Shrikes. The most supported models (∆AICc ≤ 2) are shown per each subset of variables. For categorical variables, the symbol + indicates inclusion in the model; for continuous variables, the β value is reported to illustrate the effect on species occurrence.

At the general landscape level, Woodchat Shrikes were associated with intermediate grassland cover and (small) patches of bare ground. Regarding woody vegetation, shrikes preferred areas with higher availability of isolated trees, olive groves and shrubs (especially those < 1 m). The analyses based on descriptors of grassland and bare soil habitats revealed a quadratic effect of grazed grassland, and a minor negative effect of bare ground. Considering herbaceous crops, a quadratic relationship with cereal crops different from wheat and barley was found. Among human and management traits, the most important factor was the length of dirt roads within the cell, which exerted a negative effect on species occurrence.

In the multi-level analysis, the most supported models (Δ AICc < 2) were averaged and led to the model described in Table 5. The R2 of the two most supported models was equal to c.0.27. According to the averaged model, obtained from a set of possible models including the most important factors from the different levels considered, shrike occurrence was favoured by isolated trees, shrubland and (secondarily) olive groves, and negatively affected by dirt roads (Table 5).

Table 5. Average model obtained by averaging the most supported models (Δ AICc < 2; uninformative parameters excluded) among the ones built combining the most important habitat variables from each single level (see text for details). For each variable, the coefficient in the model (± SE for the averaged model) and the relative variable importance are shown. The latter is calculated considering the sum of weights of the models in which each variable appears (Burnham and Anderson Reference Burnham and Anderson2002).

Variation partitioning suggested that woody vegetation explained a slightly higher variation in shrike occurrence than the other levels; however, when considering the variation exclusively explained by each single level, all levels are nearly equal, except for general structure, which did not explain any exclusive portion of variation (Figure 2).

Figure 2. Results of variation partitioning for the occurrence of Woodchat Shrike in terms of fractions of variation explained by the different levels. Variation in occurrence is explained by four groups of explanatory variables (the two fine-level habitat types “grassland and bare areas” and “herbaceous crops” were considered together in this analysis; see text for details).

Discussion

In birds as well as other animals, the choice of breeding habitat is a key process and can be affected by environmental factors acting at different spatial scales (e.g. Ficetola et al. Reference Ficetola, Marziali, Rossaro, De Bernardi and Padoa-Schioppa2011), or very different in nature, e.g. from land-cover type to topographical and management attributes (e.g. Chiatante et al. Reference Chiatante, Brambilla and Bogliani2014) and highly specific resources (e.g. Jedlikowski et al. Reference Jedlikowski, Brambilla and Suska-Malawska2014). Conservationists should therefore identify the scale(s) and the factors likely to be most important for habitat selection and focus on these key resources.

Several farmland bird species have been declining over decades in Europe and elsewhere, largely because of habitat changes induced by intensification and other modifications to the farming regime (Donald et al. Reference Donald, Green and Heath2001, Benton et al. Reference Benton, Bryant, Cole and Crick2002), but also because of land abandonment, which has negative impacts especially on Mediterranean birds (Preiss et al. Reference Preiss, Martin and Debussche1997, Suárez-Seoane et al. Reference Suarez-Seone, Orborne and Baudry2002); both pressures may alter habitat at different levels. Woodchat Shrike has been declining for several decades in most of its European range, which constitutes the major portion of its global distribution (BirdLife International 2015), creating concerns over its future prospects. Although conditions experienced in wintering areas and during migration are also potentially important for the species (Cramp and Perrins Reference Cramp and Perrins1993), breeding habitat availability and quality are likely to be crucial for its conservation, as they are for other shrike species (e.g. Red-backed Shrike, Lanius collurio; Brambilla et al. Reference Brambilla, Casale, Bergero, Crovetto, Falco, Negri, Siccardi and Bogliani2009a, Reference Brambilla, Casale, Bergero, Bogliani, Crovetto, Falco, Roati and Negri2010b), so it is essential to identify the factors driving species occurrence.

Here, we analysed potential determinants of shrike occurrence, considering different categories of habitat descriptors and evaluating their relative importance. Finally, from the output of this analysis, we selected the factors most likely to be involved in the habitat preferences of the species, and evaluated the most important habitat variables eventually associated with habitat selection by Woodchat Shrikes. Such habitat factors are also likely to be the most relevant for conservation through habitat preservation or restoration in breeding areas, including Italy, where the species underwent a dramatic population decline coupled with a 15% range contraction in the last decade (Nardelli et al. Reference Nardelli, Andreotti, Bianchi, Brambilla, Brecciaroli, Celada, Dupré, Gustin, Longoni, Pirrello, Spina, Volponi and Serra2015), and in the rest of the Mediterranean region, where the species usually inhabits similar semi-open landscapes.

Among the different subsets of candidate models, the one describing the general landscape structure and that characterising woody vegetation were the most parsimonious. The final models better describing habitat selection by Woodchat Shrikes included only variables from woody vegetation and from human-related traits. This could suggest that although ‘coarse’ landscape variables are able to capture most of the variation when different sets of variables are considered in isolation, the true determinants of species occurrence are likely to be represented by highly specific and fine-scale habitat traits. This is further confirmed by the variation partitioning analysis, which highlighted how the variables associated with landscape structure did not explain exclusive parts of variation, despite explaining a large amount of it in conjunction with other variables. In short, this means that landscape variables may be successfully used to predict species occurrence, but are likely less important when planning habitat management for conservation. This seems to be a rather common pattern for shrike species, which are associated with well-defined landscapes but show a strong selection for (or avoidance of) very specific habitat traits within such landscapes (Brambilla et al. Reference Brambilla, Casale, Bergero, Crovetto, Falco, Negri, Siccardi and Bogliani2009a, Chiatante et al. Reference Chiatante, Brambilla and Bogliani2014).

The association with intermediate grassland cover detected at the landscape level clearly reflects the general link with semi-open habitats, characterised by a mosaic of grassland or grassland-like cover and shrubs and trees (Cramp and Perrins Reference Cramp and Perrins1993, Nisoria 1994, Guerrieri and Castaldi Reference Guerrieri and Castaldi2000), whereas the positive selection for small extent of bare ground is likely due to the need for areas where obtaining invertebrate prey is enhanced by their high detectability and accessibility (Nisoria 1994, Schaub Reference Schaub1996, Cramp and Perrins Reference Cramp and Perrins1993). The positive effect of isolated trees, shrubs and olive groves mirrors the need for nesting and perching sites well known for that species (Cramp and Perrins Reference Cramp and Perrins1993 and references therein). Considering the other types of habitat traits, a quadratic relationship with cereal crops had been already reported from another area in southern Italy (Chiatante et al. Reference Chiatante, Brambilla and Bogliani2014), and is consistent also with anecdotal evidence reported from central Italy (Guerrieri and Castaldi Reference Guerrieri and Castaldi2000). The analyses based on descriptors of grassland and bare soil habitats revealed a quadratic effect of grazed grassland, and a negative effect of bare ground in grazed grassland. The former is fully consistent with the association with semi-open landscapes (see above), whereas the latter contrasts with the selection for small patches of bare ground found at landscape level, but it should be noted that such a negative effect of this specific type of bare ground is likely minor (the retention of the variable in the model resulted in a negligible improvement of model fit; see Table 4). The negative effect of dirty roads found in the human-related model had never been reported before, and suggests a negative effect of anthropogenic disturbance on the species.

Conservation implications

Woodchat Shrikes inhabit semi-open landscapes (on average, territories are made up of c.52% of grassland, and c.16% of arable land), within which they are attracted to shrubland and isolated trees (and secondarily to patches of olive groves), whereas they tend to avoid dirt roads. As already reported for the Red-backed Shrike (Brambilla et al. Reference Brambilla, Rubolini and Guidali2007, Reference Brambilla, Casale, Bergero, Crovetto, Falco, Negri, Siccardi and Bogliani2009a, Reference Brambilla, Casale, Bergero, Bogliani, Crovetto, Falco, Roati and Negri2010b, Ceresa et al. Reference Ceresa, Bogliani, Pedrini and Brambilla2012) and for other farmland bird species in Italy (Brambilla et al. Reference Brambilla, Guidali and Negri2008, Reference Brambilla, Guidali and Negri2009b, Reference Brambilla, Gustin and Celada2013a, Rippa et al. Reference Rippa, Maselli, Soppelsa and Fulgione2011), the maintenance of suitable conditions for the species depends on a trade-off between abandonment and intensive farming, which are both highly detrimental to species preferring semi-open landscapes. The general model built upon the results of the single-level models confirmed the importance of isolated trees, shrubland, olive groves and dirt roads, suggesting that the availability of nesting and perching sites and the lack of direct human disturbance could be key features for the species in semi-open Mediterranean landscapes.

Those results may be used for the definition of conservation measures and in particular for an updating or revision of agri-environmental measures, such as those included in the Rural Development Programme (RDP). The main implications of our findings are: i) the importance of conserving low-intensity farmland systems, which harbour a compact mosaic of open habitats, different crops and shrub/tree patches, positively selected by several species of conservation concern, including Woodlarks Lullula arborea (Brambilla and Rubolini Reference Brambilla and Rubolini2009, Brambilla et al. Reference Brambilla, Falco and Negri2012), pipits Anthus spp. (Morales et al. Reference Morales, Guerrero, Onate and Melendez2012), shrikes (Brambilla et al. Reference Brambilla, Casale, Bergero, Bogliani, Crovetto, Falco, Roati and Negri2010b, Chiatante et al. Reference Chiatante, Brambilla and Bogliani2014), buntings Emberiza spp. (Brambilla et al. Reference Brambilla, Falco and Negri2012, Brambilla Reference Brambilla2015); ii) the need to preserve some woody vegetation, and in particular shrubs and isolated trees, which have been reported to be favoured also by another threatened shrike species, the Lesser Grey Shrike Lanius minor (Chiatante et al. Reference Chiatante, Brambilla and Bogliani2014). Some national or regional RDPs include among the measures adopted for grassland conservation the removal of trees and shrubs; controlling shrub encroachment is often needed to conserve open habitats, especially when they are facing abandonment, but should be done with care (Vassilev et al. Reference Vassilev, Pedashenko, Nikolov, Apostolova and Dengler2011) to avoid the removal of breeding and perching sites for shrikes and other bird species (Nikolov Reference Nikolov2010).

Furthermore, RDPs frequently include measures promoting new roads in cultivated areas to improve access to crops and fields. Considering the negative effect of dirt roads on the species occurrence, it would be important to prevent the construction of new roads in farms hosting Woodchat Shrikes or other sensitive species, and caution should be used about road promotion in RDPs.

In conclusion, our suggestions confirm and integrate previous recommendations for Woodchat Shrikes in Mediterranean landscapes, which focused on management primarily targeted at increasing perching and nesting sites, such as isolated trees and shrubs, in open landscapes with low levels of urbanisation (Chiatante et al. Reference Chiatante, Brambilla and Bogliani2014).

Acknowledgements

This study was part of a research project funded to LIPU by the Italian Ministry for Environment (MATTM).

References

Arnold, T. W. (2010) Uninformative parameters and model selection using Akaike’s information criterion. J. Wildl. Managt. 74: 11751178.Google Scholar
Battin, J. and Lawler, J. T. (2006) Cross-scale correlations and the design and analysis of avian habitat selection studies. Condor 108: 5970.Google Scholar
Benton, T. G., Bryant, D. M., Cole, L. and Crick, H. Q. P. (2002) Linking agricultural practice to insect and bird populations: A historical study over three decades. J. Appl. Ecol. 39: 673687.CrossRefGoogle Scholar
Benton, T. G., Vickery, J. A. and Wilson, J. D. (2003) Farmland biodiversity: Is habitat heterogeneity the key? Trends Ecol. Evol. 18: 182188.Google Scholar
BirdLife International (2015) European Red List of Birds. Luxembourg: Office for Official Publications of the European Communities.Google Scholar
Birrer, S., Spiess, M., Herzog, F., Jenny, M., Kohli, L. and Lugrin, B. (2007) The Swiss agri-environment scheme promotes farmland birds: but only moderately. J. Ornithol. 148 Suppl. 2: 295303.Google Scholar
Brambilla, M. (2015) Landscape traits can contribute to range limit equilibrium: habitat constraints refine potential range of an edge population of Black-headed Bunting Emberiza melanocephala . Bird Study 62: 132136.CrossRefGoogle Scholar
Brambilla, M. and Rubolini, D. (2009) Intra-seasonal changes in distribution and habitat associations of a multi-brooded bird species: implications for conservation planning. Anim. Conserv. 12: 7177.Google Scholar
Brambilla, M., Bassi, E., Ceci, C. and Rubolini, D. (2010a) Environmental factors affecting patterns of distribution and co-occurrence of two competing raptor species. Ibis 152: 310322.Google Scholar
Brambilla, M., Casale, F., Bergero, V., Bogliani, G, Crovetto, G. M., Falco, R., Roati, M. and Negri, I. (2010b) Glorious past, uncertain present, bad future? Assessing effects of land-use changes on habitat suitability for a threatened farmland bird species. Biol. Conserv. 143: 27702778.Google Scholar
Brambilla, M., Casale, F., Bergero, V., Crovetto, G. M., Falco, R., Negri, I., Siccardi, P. and Bogliani, G. (2009a) GIS-models work well, but are not enough: Habitat preferences of Lanius collurio at multiple levels and conservation implications. Biol. Conserv. 142: 20332042.Google Scholar
Brambilla, M., Gustin, M. and Celada, C. (2013a) Species appeal predicts conservation status. Biol. Conserv. 160: 209213.Google Scholar
Brambilla, M., Falco, R. and Negri, I. (2012) A spatially explicit assessment of within-season changes in environmental suitability for farmland birds along an altitudinal gradient. Anim. Conserv. 15: 638647.Google Scholar
Brambilla, M., Fulco, E., Gustin, M. and Celada, C. (2013b) Habitat preferences of the threatened Black-eared Wheatear Oenanthe hispanica in southern Italy. Bird Study 60: 432435.Google Scholar
Brambilla, M., Guidali, F. and Negri, I. (2009b) Breeding-season habitat associations of the declining Corn Bunting Emberiza calandra – a potential indicator of the overall bunting richness. Ornis Fennica 86: 4150.Google Scholar
Brambilla, M., Guidali, F. and Negri, I. (2008) The importance of an agricultural mosaic for Cirl Buntings Emberiza cirlus in Italy. Ibis 150: 628632.Google Scholar
Brambilla, M., Rubolini, D. and Guidali, F. (2006) Factors affecting breeding habitat selection in a cliff-nesting peregrine Falco peregrinus population. J. Ornithol. 147: 428435.Google Scholar
Brambilla, M., Rubolini, D. and Guidali, F. (2007) Between land abandonment and agricultural intensification: habitat preferences of Red-backed Shrikes Lanius collurio in low-intensity farming conditions. Bird Study 54: 160167.Google Scholar
Burnham, K. P. and Anderson, D. R. (2002) Model selection and multimodel inference: A practical information-theoretic approach. 2nd edition. New York, USA: Springer-Verlag.Google Scholar
Ceresa, F., Bogliani, G., Pedrini, P. and Brambilla, M. (2012) The importance of key marginal habitat features for birds in farmland: An assessment of habitat preferences of Red-backed Shrike Lanius collurio in the Italian Alps. Bird Study 59: 327334.Google Scholar
Chamberlain, D. E., Fuller, R. J., Bunce, R. G. H., Duckworth, J. C. and Shrubb, M. (2000) Changes in the abundance of farmland birds in relation to the timing of agricultural intensification in England and Wales. J. Appl. Ecol. 37: 771788.Google Scholar
Chiatante, G. P., Brambilla, M. and Bogliani, G. (2014) Spatially explicit conservation issues for threatened bird species in Mediterranean farmland landscapes. J. Nat. Conserv. 22: 103112.Google Scholar
Cody, L. (1985) Habitat selection in birds. New York, USA: Academic Press.Google Scholar
Cramp, S. and Perrins, C. M., eds. (1993) The birds of the Western Palearctic (Vol. 7). Oxford, UK: Oxford University Press.Google Scholar
Donald, P. F., Green, R. E. and Heath, M. F. (2001) Agricultural intensification and the collapse of Europe’s farmland bird populations. Proc. R. Soc. Lond. B 268: 2529.CrossRefGoogle Scholar
Donald, P. F., Sanderson, F. J., Burfield, I. J. and Van Bommel, F. P. (2006) Further evidence of continent-wide impacts of agricultural intensification on European farmland birds, 1990–2000. Agric. Ecosyst. Environ. 116: 189196.Google Scholar
Ficetola, G. F., Marziali, L., Rossaro, B., De Bernardi, F. and Padoa-Schioppa, E. (2011) Landscape-stream interactions and habitat conservation for amphibians. Ecol. Appl. 21: 12721282.Google Scholar
Filippi-Codaccioni, O., Devictor, V., Bas, Y. and Julliard, R. (2010) Toward more concern for specialisation and less for species diversity in conserving farmland biodiversity. Biol. Conserv. 143: 14931500.Google Scholar
Fuller, R. J., Gregory, R. D., Gibbons, D. W., Marchant, J. H., Wilson, J. D., Baillie, S. R., et al. (1995) Population declines and range contractions among lowland farmland birds in Britain. Conserv. Biol. 9: 14251441.Google Scholar
Guerrieri, G. and Castaldi, A. (2000) Selezione di habitat e riproduzione dell’averla capirossa, Lanius senator, nel Lazio – Italia Centrale. Avocetta 24: 8593.Google Scholar
Hall, L. S., Krausman, P. R. and Morrison, M. L. (1997) The habitat concept and a plea for standard terminology. Wildl. Soc. Bull. 25: 173182.Google Scholar
Harris, T. and Franklin, K. (2000) Shrikes and Bush-shrikes: Including wood-shrikes, helmet shrikes, flycatcher-shrikes, philentomas, batises and wattle-eyes. London, UK: Cristopher Helm, A&C Black.Google Scholar
Jedlikowski, J., Brambilla, M. and Suska-Malawska, M. (2014) Fine-scale selection of nesting habitat in Little Crake Porzana parva and Water Rail Rallus aquaticus in small ponds. Bird Study 61: 171181.Google Scholar
Johnson, D. H. (1980) The comparison of usage and availability of measurements for evaluating resource preference. Ecology 61: 6571.Google Scholar
Jones, J. (2001) Habitat selection studies in avian ecology: A critical review. Auk 118: 557562.CrossRefGoogle Scholar
Krebs, J. R., Wilson, J. D., Bradbury, R. B. and Siriwardena, G. M. (1999) The second silent spring? Nature 400: 611612.Google Scholar
Lefranc, N. and Worfolk, T. (1997) Shrikes. A guide to the shrikes of the world. East Sussex, UK: Pica Press.Google Scholar
Morales, M. B., Guerrero, I., Onate, J. J. and Melendez, L. (2012) Inter-specific association and habitat use in a farmland passerine assemblage. Ecol. Res. 27: 691700.Google Scholar
Nardelli, R., Andreotti, A., Bianchi, E., Brambilla, M., Brecciaroli, B., Celada, C., Dupré, E., Gustin, M., Longoni, V., Pirrello, S., Spina, F., Volponi, S. and Serra, L. (2015) Rapporto sull’applicazione della Direttiva 147/2009/CE in Italia: dimensione, distribuzione e trend delle popolazioni di uccelli (2008-2013). Roma, Italy: ISPRA, Serie Rapporti, 219/2015.Google Scholar
Nikolov, S. C. (2010). Effects of land abandonment and changing habitat structure on avian assemblages in upland pastures of Bulgaria. Bird Conserv. Internatn. 20: 200213.Google Scholar
Nisoria (1994) Atlante degli uccelli nidificanti nella provincia di Vicenza. Vicenza, Italy: Gilberto Padovan.Google Scholar
Oksanen, J., Blanchet, F. G., Kindt, R., Legendre, P., Minchin, P. R., O’Hara, R. B., Simpson, G. L., Solymos, P., Stevens, M. H. M. and Wagner, H. (2013) vegan: Community Ecology Package. R package version 2.0-9. http://CRAN.R-project.org/package=vegan.Google Scholar
Orians, G. H. and Wittenburger, J. F. (1991) Spatial and temporal scales in habitat selection. Am. Nat. 137(Suppl.): 5066.Google Scholar
Peronace, V., Cecere, J. G., Gustin, M, & Rondinini, C. (2012). Lista Rossa 2011 degli uccelli nidificanti in Italia. Avocetta 36: 1158.Google Scholar
Peres-Neto, P. R., Legendre, P., Dray, S. and Borcard, D. (2006) Variation partitioning of species data matrices: estimation and comparison of fractions. Ecology 87: 26142625.Google Scholar
Preiss, E., Martin, J. L. and Debussche, M. (1997) Rural depopulation and recent landscape changes in a Mediterranean region: consequences to the breeding avifauna. Landscape Ecol. 12: 5161.Google Scholar
R Development Core Team (2013) R: A language and environment for statistical computing. Vienna, Austria: R Foundation for Statistical Computing.Google Scholar
Radišić, D., Šćiban, M. and Spremo, N. (2008) Numbers, breeding density and habitat characteristics of Woodchat Shrike Lanius senator in Pčinja river valley. Ciconia 17: 2025. In Serbian with English summary.Google Scholar
Richards, S. A. (2005) Testing ecological theory using the information-theoretic approach: examples and cautionary results. Ecology 86: 28052814.Google Scholar
Rippa, D., Maselli, V., Soppelsa, O. and Fulgione, D. (2011) The impact of agro-pastoral abandonment on the Rock Partridge Alectoris graeca in the Apennines. Ibis 153: 721734.Google Scholar
Salvo, G. (2004) Dati sulla biologia riproduttiva dell’Averla capirossa, Lanius senator, in Sicilia. Riv. Ital. Orn. 74: 8385.Google Scholar
Schaub, M. (1996) Territory choice and settlement of Woodchat Shrikes Lanius senator in northern Switzerland. Ornithologische Beobachter 93: 163168.Google Scholar
Shochat, E., Abramsky, Z. and Pinshow, B. (2002) Breeding species diversity in the Negev: effect of scrub fragmentation by planted forests. J. Appl. Ecol. 38: 11351147.Google Scholar
Siriwardena, G. M., Baillie, S. R., Buckland, S. T., Fewster, R. M., Marchant, J. H. and Wilson, J. D. (1998) Trends in the abundance of farmland birds: a quantitative comparison of smoothed Common Birds Census indices. J. Appl. Ecol. 35: 2443.Google Scholar
Suarez-Seone, S., Orborne, P. E. and Baudry, P. E. (2002) Responses of birds of different biogeographic origins and habitat requirements to agricultural land abandonment in northern Spain. Biol. Conserv. 105: 333344.Google Scholar
Tucker, G. M. and Evans, M. I. (1997) Habitats for birds in Europe: A conservation strategy for the wider environment. Cambridge, UK: BirdLife International.Google Scholar
Vassilev, K., Pedashenko, H., Nikolov, S. C., Apostolova, I. and Dengler, J. (2011). Effect of land abandonment on the vegetation of upland semi-natural grasslands in the Western Balkan Mts., Bulgaria. Plant Biosystems 145: 654665.Google Scholar
Vickery, J. and Arlettaz, R. (2012) The importance of habitat heterogeneity at multiple scales for birds in European agricultural landscapes. Pp. 177–204 in: Fuller, R. J., ed. Birds and habitat. Relationships in changing landscapes. Cambridge, UK: Cambridge University Press.Google Scholar
Figure 0

Table 1. Factors affecting Woodchat Shrike occurrence and habitat selection according to the available literature.

Figure 1

Figure 1. Location of study areas in Italy. Each study area included 7–12 plots within which fieldwork was carried out.

Figure 2

Table 2. Habitat variables considered in this study to evaluate habitat selection according to five different levels of habitat traits.

Figure 3

Table 3. Average features of shrike territories and control plots; * indicates significant (P < 0.05) differences (assessed by means of a t-test on arc-sin square-root transformed variables for land cover and road length, and by χ2-test for grazing occurrence).

Figure 4

Table 4. Candidate models reflecting different potential determinants of habitat selection in Woodchat Shrikes. The most supported models (∆AICc ≤ 2) are shown per each subset of variables. For categorical variables, the symbol + indicates inclusion in the model; for continuous variables, the β value is reported to illustrate the effect on species occurrence.

Figure 5

Table 5. Average model obtained by averaging the most supported models (Δ AICc < 2; uninformative parameters excluded) among the ones built combining the most important habitat variables from each single level (see text for details). For each variable, the coefficient in the model (± SE for the averaged model) and the relative variable importance are shown. The latter is calculated considering the sum of weights of the models in which each variable appears (Burnham and Anderson 2002).

Figure 6

Figure 2. Results of variation partitioning for the occurrence of Woodchat Shrike in terms of fractions of variation explained by the different levels. Variation in occurrence is explained by four groups of explanatory variables (the two fine-level habitat types “grassland and bare areas” and “herbaceous crops” were considered together in this analysis; see text for details).