2 DeRosaetal
2 DeRosaetal
2 DeRosaetal
net/publication/311844169
CITATION READS
1 747
3 authors:
49 PUBLICATIONS 49 CITATIONS
Springer Nature
27 PUBLICATIONS 213 CITATIONS
SEE PROFILE
SEE PROFILE
Mirko Di Febbraro
Università degli Studi del Molise
117 PUBLICATIONS 769 CITATIONS
SEE PROFILE
Some of the authors of this publication are also working on these related projects:
Census of the breeding population of Lesser Kestrel (Falco naumanni) in Sardinia View project
Ecologia spaziale degli uccelli pelagici nell'AMP Regno di Nettuno View project
All content following this page was uploaded by Davide De Rosa on 23 December 2016.
1
Environmetrics Lab, Dipartimento di Bioscienze and Territorio (DiBT), Università del Molise - 86090 Pesche (IS), Italy
2
ARDEA, Associazione per la Ricerca, la Divulgazione e l’Educazione Ambientale - Via Ventilabro 6, 80126 Napoli,
Italy
3
Department of Biology, Colorado State University - Fort Collins CO 80523, USA
Abstract – The Black Woodpecker has a wide northern Palearctic range, but in much of western Europe it has a highly fragmented dis-
tribution. As these isolated populations are vulnerable to land cover change, a better understanding of the factors driving their presence
is needed in the perspective of conservation. To shed light on the habitat preferences of the Black Woodpecker in the southern part of its
range, we carried out a research in a large protected area, the Cilento and Vallo di Diano National Park (the first investigation on this spe-
cies anywhere in southern Italy). We used a playback technique to detect occurrence of Black Woodpeckers, as well as sympatric wood-
pecker species, and recorded environmental characteristics that might explain the bird’s occurrence, such as tree species and stand age.
The “extent of occurrence” of the Black Woodpecker was calculated, and its habitat preferences investigated. Occurrence of the species
was clearly dependent on beech dominance, in contrast with the other woodpecker species, with an estimated total population for Nation-
al Park 18.56 (± 6.19) pairs in beech-dominated woods. Habitat preferences were narrower compared to northern European populations,
highlighting the dependence of this species on beech forests in southern Europe. Since beech forests in the Mediterranean region are pre-
dicted to recede due to climate change, our findings highlight the vulnerability of southern Black Woodpecker populations.
Introduction whose global range extends from Europe through the taiga
forest belt into Asia, with a core area between 62° and 69°
Due to their ecological importance and high sensitiv- N (Winkler et al. 1995). Black Woodpecker usually uti-
ity to disturbance that affects other components of avi- lizes tall and large trunks of many coniferous and broad-
fauna, woodpeckers have been identified as a particular- leaved trees forming extensive unbroken forests (Glutz &
ly important conservation target (Drever & Martin 2010, Bauer 1980, Cramp 1985). Habitat use is possibly related
Mikusiński et al. 2001). Woodpeckers are good indicators to its peculiar food requirements, especially carpenter ants
of forest biodiversity as their presence is representative of (Formicidae) (Cuisin 1988, Ceugniet 1989, Pechacek &
the condition of the forest and its overall environmental Krisvtin 1993, Rolstad & Rolstad 2000).
health (Angelstam 1990, Mikusiński 1997). The consid- In Europe it ranges from the Arctic Circle region to
erable reduction of deadwood not only at the stand but at southern Italy, using a wide variety of tree species for
the landscape level is considered to be one of the principal feeding and nesting (Angelstam 1990). In contrast with
causes of biodiversity loss in managed forest ecosystems this overall wide distribution, however, populations in
worldwide. The presence of woodpecker assemblages and the southern part of the range, including the Iberian, Ital-
their association with standing deadwood for both nesting ian and Balkan peninsulas, are often small and fragment-
and foraging is used to emphasize the importance of the ed (Brichetti & Fracasso 2007, Aierbe et al. 2014). Ani-
entire range of snag degradation stages for maintenance of mal populations at range margins are especially vulnerable
key ecological processes in habitat remnants of managed to environmental modification, e.g. due to climate change,
landscapes (Drapeau et al. 2009). and may differ from core populations in ecological prefer-
One of the most widely distributed woodpecker spe- ences and therefore conservation requirements (Hampe &
cies is the Black Woodpecker Dryocopus martius L., Petit 2005). This is exemplified by the contrast between
the large Black Woodpecker population in northern Italy, Beyond this elevation, beech F. sylvatica dominates, usu-
with an estimate of 1,300-3,700 breeding pairs mainly lo- ally preceded by a band of Neapolitan Alder.
cated in forested Alpine and sub-Alpine regions, and the
small fragmented populations in the southern Apennines Woodpecker survey
(Nardelli et al. 2015). For this study, the park was divided into 2.5 × 2.5 km
Northern Italian population appears in expansion and grid plots using ArcGIS. Although a Black Woodpeck-
in the last twenty years there has been a progressive colo- er’s home range may vary in size, depending on the re-
nization of hilly areas and, in some cases, planitial, by this gion, plot size was based on an average of the minimum
species, once relegated to the mountain forests (Nardelli et home range sizes known for Europe (western Italian Alps:
al. 2015), whereas poor information is known about south- 3 km2, Bocca & Rolando 2003; Sweden: 5 km2, Tjernberg
ern population. et al. 1993; Norway: 1 km2, Rolstad et al. 2000; Denmark:
Reasons for this might go beyond climatic factors, as 1.30 km2, Johansen 1989; Hungary: 2 km2, Gorman 2011).
for instance only occasional records have been reported A first selection of 104 grid plots was made, on the criteri-
from the Abruzzo, Latio and Molise National Park, despite on that each plot had to be at least 50% covered by wood-
the abundance of suitable old beech Fagus sylvatica for- land. All woodland types in the park were included, as the
est habitat. A long history of hunting may explain why the Black Woodpecker is known to occur in many forest types
Black Woodpecker is largely absent today from much of (Angelstam & Mikunsinski 1994). Subsequently, 34 plots
peninsular Italy (P. Harris in Gorman 2011). Small size of were randomly selected and investigated (Fig. 1).
suitable forest fragments is also an important limiting fac- The plots were surveyed one time between February
tor for this species, as shown in the Iberian Peninsula (Gar- and March 2012, as courtship in the Black Woodpecker
mendia et al. 2006). can start anytime from February (Cramp 1985). The play-
In order to improve our understanding of Black Wood- back technique (Falls 1992, Douglas and Mennill 2010)
pecker habitat requirements in southern Europe, we inves- was used to record the number of individuals of the Black
tigated its presence and habitat in a large protected area in Woodpecker. Sounds used for play-back were taken from
southern Italy using the playback technique. The specific Roché (2008). An initial play-back was followed by five
objectives of our study were: a) to test the relationship be- minutes of silence in order to detect any response; the op-
tween the species occurrence and forest structure charac- eration was then repeated once. The procedure was used
teristics, b) to determine its extent of occurrence (EOO). for other Woodpecker species in this order: Leiopicus me-
dius, Dryobates minor, Dendrocopos major, Picus viridis.
Although the sampling period may not have been optimal
Materials and methods for these species, European woodpeckers are essentially
non-migratory and sedentary birds (Gorman 2011), and so
Study area we opportunistically also detected their presence.
The study area, “Parco Nazionale del Cilento, Vallo di Di- Since our main target species is highly territorial, the
ano e Alburni”, henceforth referred to as the Cilento Na- location of breeding territories was used as a census proxy,
tional Park, is one of the largest protected areas in south- i.e. a targeted census (Svenssons 1979, Bibby et al. 1992).
ern Europe, it extends from sea level to 1899 m a.s.l. of We investigated the effect of the different forest variables
Cervati mountains and is of outstanding ecological inter- on Black Woodpecker during the breeding season. There-
est due to its diverse range of habitats and the high bio- fore, the conclusions that may be drawn from our data are
diversity it supports (La Valva & Carrabba 1999, Per- relevant to the territories used for breeding and not neces-
siani et al. 2010, Ravera & Brunialti 2012, Romano et sarily to the habitats exploited during other seasons; how-
al. 2010). Established in 1991, it covers an area of about ever all woodpeckers are predominantly sedentary, espe-
180,000 hectares encompassing a wide range of vegeta- cially the fragmented population (Gorman 2011).
tion types. About 1,800 plant species were surveyed in the
park, among which Yew Taxus baccata, Holly Ilex aqui- Habitat evaluation
folium and an Italian endemism, the Neapolitan Alder Al- In each plot a set of environmental variables was meas-
nus cordata. The inland area is covered by deciduous hard- ured, including forest characteristics and elevation. For-
wood forests, with a mix of Turkey Oak Quercus cerris, est type was identified based on the tree species assem-
Downy Oak Quercus pubescens, Maples Acer spp., Horn- blage as one the following: Beech, Chestnut, Holm Oak-
beam Ostrya spp., Manna Ash Fraxinus ornus and Chest- Alder, deciduous oaks, deciduous oaks and beech, mixed-
nut Castanea sativa up to about 1000 m above sea level. oak, and mixed-Mediterranean forest. Stand maturity was
64
Black Woodpecker in Mediterranean forests
Figure 1. Study area in Cilento National Park, Southern Italy. Red squares are the parcels established at the beginning of the study, with
numbers marking those selected for the field survey.
classified on a discrete scale from 1 to 3, from recent (< 15 evaluate the EOO of the Black Woodpecker in Cilento.
years) coppice to mature stands. Likewise, beech density This technique was chosen because it is one of the best
and presence of dead trees were classified from 1 (no liv- known in statistics (Börger et al. 2006) and is the main
ing beech trees and no dead trees detected, respectively) choice for modern studies of animal-habitat relationships
to 3 (beech-dominated, i.e. virtually only beech as arbo- (Marzluff et al. 2004). We adopted the fixed kernel method
real species, and over 10 dead trees per km2 respectively). for bivariate data as recommended by Seaman & Powell
(1996). Regarding the choice of the bandwidth, a parame-
Analysis ter showing a very strong effect on kernel estimation (Sea-
Occurrence, habitat use and population size of the Black man & Powell 1996, Burgman & Fox 2003, Börger et al.
Woodpecker and other co-occurring Picidae were inves- 2006), as it regulates its smoothing, we followed the ‘rule
tigated in a subset of 34 plots randomly selected from the of thumb’ proposed by Scott (1992):
initial set of 104 grids. To highlight potential differences - hopt = 1.06An-1/5 where hopt = optimal bandwidth;
in plot use among the five woodpecker species surveyed, - n = sample size;
we performed Principal Component Analysis (PCA) on - A = min (standard deviation, interquartile range/1.34).
species presence/absence data (five species) using plots as Following the recommendation of Börger et al. (2006)
sites. A generalised linear model (GLM) with a binomial to adopt a utilization distribution isopleth between 90%
distribution of errors and a clog-log link function was used and 50% (the range in which kernel estimators perform
to model the probability of recording a Black Woodpeck- better), we chose a 50% isopleth to calculate the external
er in relationship to beech tree density, presence of dead range boundaries. Statistical analyses and kernel estima-
trees and woodland age. Each of the explanatory variables tions were performed in R 2.15.3 (R Development Core
was dropped in turn and subjected to an analysis of devi- Team, 2012; maptools and Genkern packages).
ance using Chi-square tests to assess its significance. Val- In addition we estimated density of Black Wood-
ues were considered significant at a 5% probability level, pecker in beech-dominated woodland dividing the num-
and variance explained by the model was calculated as the ber of birds by the explored area and assuming the every
ratio of the difference between null and residual deviance birds contacted represented a pair (see Douglas & Men-
to null deviance (Zuur et al. 2009). GLM estimates are pre- nill 2010). Based on the total number of plots identified as
sented as means ± standard errors. beech-dominated, pairs were estimated to occur across all
Based on our field data, we used kernel estimation to beech forests in the Cilento National Park.
65
De Rosa et al.
Results first two axes together account for 63% of the variation in
the woodpecker assemblage. The first axis was positively
Black Woodpeckers were recorded in 10 out of 34 plots. correlated with beech density (Spearman’s r = 0.34, p =
In nine of these only one individual was found, and only 0.05), while neither the first nor the second axis were sig-
in one two distinct males were recorded. Most individu- nificantly correlated with stand age or the presence of dead
als were found in beech-dominated plots, with only three trees (r < 0.1, p > 0.1) According to the GLM, the prob-
found in mixed-oak forest. PCA suggests that the Black ability of recording a Black Woodpecker was significantly
Woodpecker shows markedly different habitat preferences related to increasing beech dominance (p < 0.01), increas-
from the other four species in the study area (Fig. 2). The ing from 8.9 ± 6.5% in areas without living F. sylvatica,
1
PCA axis 2 (R2 = 0.27)
-1
-1 0 1
Figure 2. Principal component analysis based on presence of five woodpecker species in the study area. Black dots represent the sampled
parcels (n = 34). Parcels where the same bird numbers were recorded overlap. Arrows refer to the species, with the following abbrevia-
tions: B = Black Woodpecker Dryocopus martius; GS = Great Spotted Woodpecker Dendrocopos major; MS = Middle Spotted Wood-
pecker Leiopicus medius; LS = Lesser Spotted Woodpecker Dryobates minor; G = European Green Woodpecker Picus viridis. The ei-
genvalues for the first, second, and third axes are, respectively, 0.355, 0.266 and 0.170.
66
Black Woodpecker in Mediterranean forests
Figure 4. Extent of occurrence of the Black Woodpecker in Cilento National Park study area.
67
De Rosa et al.
ern Italy may be precarious (Calvario et al. 2011, Peronace (Servizio Conservazione Natura, Parco Nazionale del Cilento e
et al. 2012). The fact that it almost only occurred in beech Vallo di Diano) for authorising field work and the CTA of Vallo
della Lucania. We are also grateful to Egidio Fulco and Giuseppe
woodland may be especially critical, since this habitat is Martino for information on the Black Woodpecker populations of
predicted to decrease dramatically in the coming decades Basilicata and Calabria, and to Ivan Maggini, Barry McMahon
and Mark Walters for constructive comments that improved the
due to ongoing climate change (Innangi et al. 2015).
manuscript. Two anonymous referees contributed to improve ear-
Such a small local population might be a consequence lier version of the manuscript.
of the current management of beech woodland in the Ci-
lento National Park, which may need to be modified if
conservation of rare beech woodland species, such as the References
Black Woodpecker, is of the overriding consideration.
More studies in the other areas of southern Europe where Aierbe T., Beñaran H., Fernández-García J.M., Gracianteparaluc-
eta A., Hurtado R., Olano M. & Ugarte J., 2014. Distribution
the Black Woodpecker occurs are required to understand and abundance of the Black Woodpecker (Dryocopus mar-
the links between its demographic status and forest man- tius) in a recently colonized area in SW Europe. Woodpeck-
agement. An important contribution may come from stud- ers in a changing world’ International Conference.
Angelstam P. 1990. Factors determining the composition and per-
ies on juvenile dispersals and adult movements using tele- sistence of local woodpecker assemblages in taiga forest in
metric techniques (Rostald et al. 1998, Bocca et al. 2007). Sweden - a case for landscape ecological studies. Pp. 147-164
These data could be linked to the spatial expansion of in: Carlson A. & Aulén G. (eds) Report 17, Uppsala.
Angelstam P. & Mikusiński G.,1994. Woodpecker assemblages in
Woodpecker species in the last decades of the 20th cen- natural and managed boreal and hemiboreal forest a review.
tury in Mediterranean Europe, which seems to be closely Ann. Zool. Fenn. 31: 157-172.
related to forest maturation, following large-scale decline Bibby C.J., Collar N.J., Crosby M.J., Heath M.F., Imboden C.,
Johnson T.H., Long A.J., Stattersfield A.J. & Thirgood S.J.,
in traditional uses (Gil-Tena et al. 2013). Notably, howev- 1992. Putting Biodiversity on the Map: Priority Areas for
er, our findings show that Black Woodpecker occurrence Global Conservation. International Council for Bird Conser-
vation, Cambridge, UK.
in Cilento National Park is linked to the presence of beech
Bocca M. & Rolando A., 2003. The ecology of the Black Wood-
but not to stand age. pecker in Mont Avic Natural Park (Italian Western Alps). P.
To our knowledge, no other studies on the density of 48 in: Pechacek P. & d’Oleire-Oltmanns W. (eds), 5th Int.
Woodpecker Symp.
this species have been carried out in other areas of south-
Bocca M., Carisio L. & Rolando A. 2007. Habitat use, home rang-
ern Italy. Indeed, ours is a first attempt to quantify the hith- es and census techniques in the Black Woodpecker (Dryoco-
erto unknown population size of the Black Woodpecker pus martius) in the Alps. Ardea 95 (1): 17-29.
Boitani L., Corsi F., Falcucci A., Maiorano L., Marzetti I., Masi
in this important region near the margin of its range. We
M., Montemaggiori A., Ottaviani D., Reggiani G. & Rond-
suggest that Black Woodpeckers may be more widespread inini C., 2002. Rete Ecologica Nazionale. Un approccio alla
in southern Italy than is currently believed, although the conservazione dei vertebrati italiani. Università di Roma “La
Sapienza”, Dipartimento di Biologia Animale e dell’Uomo,
absence of environmental corridors connecting to larger, Ministero dell’Ambiente, Direzione per la Conservazione
more stable populations in continental Europe make south- della Natura, Istituto di Ecologia Applicata, Roma.
ern populations particularly vulnerable to climate change Börger L., Franconi N., De Michele G., Gantz A., Meschi F.,
Manica A., Lovari S. & Coulson T.I. M., 2006. Effects of
and land use change. sampling regime on the mean and variance of home range
In conclusion, although the Black Woodpecker is ex- size estimates. J. Animal Ecol. 75 (6): 1393-1405.
panding in continental Europe (Vos et al. 2008), isolat- Burgman M.A. & Fox J.C., 2003. Bias in species range estimates
from minimum convex polygons: implications for conserva-
ed populations in the southern part of the range, such as tion and options for improved planning. Animal Conserv. 6
southern Italy, are to be considered endangered, with a (1): 19-28.
lack of data on context-specific ecological requirements Calvario E., Gustin M., Sarrocco S., Gallo-Orsi U., Bulgarini F. &
Fraticelli F., 2011. Nuova Lista Rossa degli uccelli nidificanti
hindering their conservation. In addition, considering the in Italia. LIPU-WWF.
Black Woodpecker as an indicator species, like those spe- Ceccarelli P., Agostini N. & Milandri M., 2003. Presenza del Pic-
cies whose ecological requirements guarantee the exist- chio nero nelle Foreste Casentinesi: una nuova specie per la
check-list regionale. Asoer, Notizie, Bologna.
ence of particular environmental conditions (Campbell & Cramp 1985. Handbook of the birds of Europe, the Middle East
Lack 1985) would also help to protect a wide range of ani- and North Africa: the birds of the Western Palearctic. Vol.
mal species that share its habitat, including important wild IV Terns to woodpeckers. Oxford University Press, Oxford.
Douglas S.B. & Mennill D.J., 2010. A review of acoustic play-
pollinators, rare Coleoptera and bats. back techniques for studying avian vocal duets. J. Field Orn.
81: 115-129.
Drapeau P., Nappi A., Imbeau L. & Saint-Germain M., 2009.
Acknowledgments – The authors would like to thank Andrea Standing deadwood for keystone bird species in the eastern
Senese, Alfredo Galietti, Rosario Balestrieri and Valerio Giovan- boreal forest: Managing for snag dynamics. Forestry Chroni-
ni Russo for invaluable field work assistance, and Laura De Risio cle 85: 227-234.
68
Black Woodpecker in Mediterranean forests
Drever M.C. & Martin K., 2010. Response of woodpeckers to Persiani A.M., Audisio P., Lunghini D., Maggi O., Granito V.M.,
changes in forest health and harvest: implications for conser- Biscaccianti A. B., Chiavetta U. & Marchetti M., 2010. Link-
vation of avian biodiversity. Forest Ecol. Manage 259 (5): ing taxonomical and functional biodiversity of saproxylic
958-966. fungi and beetles in broad-leaved forests in southern Italy
Falls J.B, 1992. Playback: a historical perspective. Playback and with varying management histories. Plant Biosystems 144:
studies of animal communication. Springer US, 11-33. 250-261.
Garmendia A., Cárcamo S. & Schwendtner, O. 2006. Forest man- Pinheiro J., Bates D., DebRoy S. & Sarkar D., 2012. R Develop-
agement considerations for conservation of black woodpeck- ment Core Team. nlme: Linear and nonlinear mixed effects
er Dryocopus martius and white-backed woodpecker Den- models, 2012. URL: http://CRAN. R-project. org/package=
drocopos leucotos populations in Quinto Real (Spanish West- nlme. R package version, 3-1.
ern Pyrenees). Pp. 339-355 in: Forest Diversity and Manage- Ravera S. & Brunialti G., 2012. Epiphytic lichens of a poorly ex-
ment, Springer, Netherlands. plored National Park: Is the probabilistic sampling effective
Gil-Tena A., Brotons L., Fortin M.J., Burel F. & Saura S., 2013. to assess the occurrence of species of conservation concern?
Assessing the role of landscape connectivity in recent wood- Plant Biosystems, (ahead-of-print), 1-10.
pecker range expansion in Mediterranean Europe: forest man- Roché J.C., 2008. Bird Songs and Calls of Britain and Europe.
agement implications. Eur. J. For. Res. 132 (1): 181-194. Rolstad J., Rolstad E. & Saeteren O., 2000. Black Woodpecker
Gorman G., 2011. Black Woodpecker. A monograph on Dryoco- nest sites: characteristics, selection and reproductive success.
pus martius. Lynx ed., Barcelona. J. Wildl. Manage. 64 (4): 1053-1066.
Hampe A. & Petit R.J., 2005. Conserving biodiversity under cli- Rolstald J., Majewski P. & Rolstad E., 1998. Black Woodpecker
mate change: the rear edge matters. Ecology Letters 8 (5): use of habitats and feeding substrates in a managed Scandina-
461-467. vian forest. J. Wildl. Manage. 62 (1): 11-23.
Innangi M., d’Alessandro F., Fioretto A. & Di Febbraro M., 2015. Romano A., Ventre N., De Riso L., Pignataro C. & Spilinga C.
Modeling distribution of Mediterranean beech forests and 2010. Amphibians of the “Cilento e Vallo di Diano” National
soil carbon stock under climate change scenarios. Climate Park (Campania, Southern Italy): update check list, distribu-
Res. 66(1): 25-36. tion and conservation notes. Acta Herpetologica 5: 233–244.
Johansen B.T. 1989. Sortspættens Dryocopus martius be- Scott D.W., 1992. Multivariate Density Estimation: Theory, Prac-
standsstørrelse, territoriestørrelse og yngleresultater i Tis- tice, and Visualization. Wiley, Hoboken, New Jersey, USA.
vilde Hegn, Nordsjælland, 1977-1986. Dansk. Orn. Foren. Seaman D.E. & Powell R.A., 1996. An evaluation of the accuracy
Tidsskr. 83: 113-118. of kernel density estimators for home range analysis. Ecology
La Valva V. & Carrabba M.C., 1999. Dalle lave del Vesuvio al- 77 (7): 2075-2085.
le vette del Cilento: Aspetti floristici e vegetazionali dei due Svensson S.E., 1979. Census efficiency and number of visits to
Parchi Nazionali della Campania. Pp. 297-329 in: Lucarelli a study plot when estimating bird densities by the territory
F. (a cura di), La rete Mab nel Mediterraneo: Parchi Nazi- mapping method. J. Appl. Ecol. 16: 61-68.
onali del Cilento e Vallo di Diano e del Vesuvio. Studio Idea Tjernber M., Johnsson K. & Nilsson S. G., 1993. Density varia-
Editrice. tion and breeding success of the Black Woodpecker Dryoco-
Marzluff J., Millspaugh J., Hurvitz P. & Handcock M.A., 2004. pus martius in relation to forest fragmentation. Ornis Fenn.
Relating resources to a probabilistic measure of space use: 70: 155-162.
forest fragments and Steller’s Jays. Ecology 85:1411-1427. Vos C.C., Berry P., Opdam P., Baveco H., Nijhof B., O’Hanley J.
Meschini E. & Frugis S., 1993. Atlante degli uccelli nidificanti & Kuipers H., 2008. Adapting landscapes to climate change:
in Italia. Istituto Nazionale Fauna Selvatica, Ozzano Emilia. examples of climate proof ecosystem networks and priority
Mikusinski G., 1997. Winter foraging of the Black Woodpecker adaptation zones. J. Appl. Ecol. 45 (6): 1722-1731.
Dryocopus martius in managed forests in south-central Swe- Winkler H., Christie D. & Nurney D., 1995. Woodpecker: a guide
den. Ornis fennica 74:161-166. to the woodpecker of the world. Houghton Mifflin, New
Mikusiński G., Gromadzki M. & Chylarecki P., 2001. Woodpeck- York.
ers as indicators of forest bird diversity. Conserv. Biol. 15(1): Zuur A.F., Ieno, E.N. Walker, N.J. Saveliev A.A. & Smith G.M.,
208-217. 2009. Mixed Effects Models and Extensions. Ecology with R.
Peronace V., Cecere J.G., Gustin M. & Rondinini C., 2012. Lis- Springer, New York, NY, USA
ta Rossa 2011 degli Uccelli nidificanti in Italia. Avocetta 36:
11-58. Associate Editor: Michelangelo Morganti
69