Introduction
The conservation status of a species can be measured through a favourable reference value (FRV) that indicates range, habitat, population size and/or other parameters likely to ensure its long-term persistence in all or part of its range. These values may be used to evaluate the conservation status of populations and they can aid effective conservation planning because species or populations whose parameters strongly deviate from FRVs may be considered more ‘at risk’. Therefore, it is at least conceivable that defining reference values for a population, or a part of the range of a target species, could help the assessment and monitoring of its conservation status.
A critical issue for species conservation is population size. Below a certain threshold, the survival of an organism is likely to be compromised by several factors affecting survival and reproduction (Soulé Reference Soulé1987). Identifying the minimum number of individuals needed to ensure the long-term persistence of a species or population, taking into account the impact of deterministic and stochastic factors, has become a common feature in conservation planning. Such an approach is known as population viability analysis (PVA), where demographic models are built to improve understanding of animal populations and/or to complement field studies in the development of conservation strategies (Katzner et al. Reference Katzner, Bragin and Milner-Gulland2006). PVAs show the impact of demographic, environmental, genetic and other factors on the dynamics of wildlife populations. Very often, PVA models are used to identify a minimum viable population (MVP; Lacy Reference Lacy2000a; Reed et al. Reference Reed, O’Grady, Brook, Ballou and Frankham2003). PVA models often face problems when dealing with large numbers of individuals and for very abundant species they may result in analyses that are excessively time-consuming, less reliable, and have highly variable outputs (Lacy et al. Reference Lacy, Borbat and Pollak2005). PVA models can also help to identify ‘optimal’ population sizes in different situations; for example, they can identify MVPs for species facing extinction risk, or for populations large enough to ensure long-term persistence even in less favourable conditions, or populations subject to harvesting where the species is not currently threatened with extinction. PVA models can show which traits are more important for conservation, e.g. reproductive or survival rates, through elasticity and sensitivity analyses (Naujokaitis-Lewis et al. Reference Naujokaitis-Lewis, Curtis, Arcese and Rosenfeld2009).
When PVA models are unlikely to be useful or are impossible to develop due to a lack of data, other population parameters may help to identify conservation targets and assess conservation status. In particular, breeding density (the number of individuals, pairs, territories or nests per unit area) can give a first indication of the conservation status of a species, indirectly representing abundance and/or habitat suitability for a given species over a given area.
Conservation targets under EU directives
The European Union’s Habitats Directive (92/43/CEE) and Wild Birds Directive (79/409/CEE) represent the main tools for EU member states to use in stemming the decline in biodiversity. The Habitats Committee (in the Directorate-General Environment within the European Commission) has defined guidelines for monitoring and reporting the conservation status of species and habitats of EU interest within the framework of Article 17 of the Habitats Directive (European Commission 2005). These guidelines require member states to evaluate the conservation status of species, to assess their demographic trends and to provide a population figure defined as a ‘Favourable Reference Value’ for each species. The guidelines state that an FRV should be defined for each species and habitat by each member state. According to the general principles and indications of the Habitats Directive and to our own interpretation, the FRV should indicate a population figure (or another demographic reference value, if population size is unfeasible) which is likely to ensure long-term persistence of the species in its range or a part thereof.
An FRV-based assessment of conservation status has to be completed for each species. FRVs (e.g. for range, area covered, population size) should be established on a scientific basis, based on the best available conservation knowledge and in a transparent way. ‘Best expert judgement’ may be used to define it in absence of other data (European Commission 2005).
Establishing FRVs must be distinguished from establishing concrete targets: setting targets is the translation of reference values into operational, practical and feasible short-, medium- and long-term targets/milestones. Member states are therefore encouraged to include FRVs in their monitoring reports, as the establishment of such values will strongly support discussions on status assessment and priority setting at the biogeographical or national level. Until now, except for a concise report on species conservation status in Denmark (Pihl et al. Reference Pihl, Clausen, Laursen, Madsen and Bregnballe2006), none of the member states of the EU has provided a comprehensive evaluation of FRVs for its breeding bird species included in Annex I of the Wild Birds Directive.
In this paper, we report a conservation planning exercise carried out at a national scale (Italy) for the definition of FRVs for bird species conservation at a population scale within the framework set by the EU Habitats Directive and Wild Birds Directive. We develop a conceptual framework for the definition of population FRVs for birds and for the establishment of FRVs as a long-term reference for conservation objectives for the Italian populations of species included in Annex I of the Wild Birds Directive.
Methods
FRVs were defined only for species regularly breeding in Italy, except for 12 species undergoing current demographic expansion after recent colonisation (30 years). We decided not to formulate an FRV for these species because population expansion following colonisation may lead to unpredictable outcomes in most cases.
We developed a protocol for defining FRVs based on demographic trends, current populations, and the minimum viable population concept. The entire process required one year of work and involved analysis of more than 1,900 bibliographic references.
Italian breeding birds were assigned to categories according to their distribution (species that were widespread or those that occurred in discrete populations) and abundance (≤ 2,500 breeding pairs or > 2,500 breeding pairs). The cut-off value was set at 2,500 pairs as it seemed reasonable to use PVAs only for small populations, for which population modelling gives the most reliable and accurate results (Lacy et al. Reference Lacy, Borbat and Pollak2005). Species with fewer than 2,500 pairs were further divided into two groups: i) species with a single population; a single FRV for the Italian population was provided and ii) species with several non-contiguous populations, i.e. range distribution structured into discrete units (populations, colonies, or group of colonies) with virtually no gene flow between them; in this case an FRV was provided for each main unit.
Species with more than 2,500 pairs were also divided into two groups according to their distribution. The first group contained species whose range is subdivided into discrete units (populations, colonies, or group of colonies) with populations < 2,500 breeding pairs; a specific FRV was provided for each unit. The second group included species with a large and more or less continuous range (all possible discrete units relevant from a conservation point of view, i.e. excluding units of minimal size, > 2,500 breeding pairs); for this group, the FRV was expressed in terms of breeding density at different spatial scales for non-colonial species. An FRV is not provided for colonial species with population units > 2,500 breeding pairs as it does not appear feasible to define a density-based reference value that may not reflect (or be correlated to) population size.
The program Vortex (version 9.72) was used to evaluate the extinction risk for species and populations with fewer than 2,500 breeding pairs (Lacy Reference Lacy1993, Reference Lacy2000b; Lacy et al. Reference Lacy, Borbat and Pollak2005). The analyses included demographic parameters, (age structure, mortality at different ages), reproductive parameters (type of reproductive system, first and maximum age of reproduction, sex ratio at birth, percentage of males and females taking part in reproduction, percentage of successful nests and/or number of young fledged) and environmental parameters (occurrence of catastrophes, carrying capacity of the habitat), as requested by modelling procedures in Vortex. Species and populations with fewer than 2,500 breeding pairs were divided into three groups.
The first group included species or populations that showed a constant positive trend (increase, or increase alternating with/followed by stability) during the last 20–30 years. For this group, using known parameters for each species, we calculated the probability of extinction (P) using the highest estimate of the current population (e.g. taking 200 when the estimated population is 100–200 pairs). If P ≤ 0.01 (no more than 1% chance of extinction) in the next 100 years, the estimated population is taken as the FRV for the given species or population; if P > 0.01 (probability of extinction greater than 1%) we calculated the minimum viable population (MVP) as the population showing a 99% probability of persistence over 100 years (P = 0.01 in the next 100 years) and took this as the FRV for that species or population. Whenever possible, different simulations were carried out with varying parameters within the range of known values for a given species.
The second group included species or populations with unfavourable or unknown trends. In this case, we calculated the population with 1% probability of extinction over the next 100 years using the less favourable values for breeding and reproductive parameters found in the literature (but avoiding extreme values linked to exceptionally unfavourable years or conditions) and took that population figure as the FRV. Such a population figure should ensure species persistence even under less favourable scenarios.
The third group included species or populations with a high probability of extinction in the next 100 years; we calculated the MVP (population with 1% probability of extinction in the next 100 years) using the average demographic and breeding parameters known for the species and took this as the FRV. The unfavourable conditions shown by these species/populations made it unrealistic to calculate an MVP using the less favourable values.
Calculating FRVs requires a large set of demographic and reproductive parameters that are crucial for population modelling analyses. Therefore, an FRV was calculated only when the main parameters needed for population analysis were known. In some cases, data for closely related species belonging to the same genus and with comparable ecology and breeding density in suitable habitats were used instead (e.g. Katzner et al. Reference Katzner, Bragin and Milner-Gulland2006). Otherwise, an FRV was not calculated. When only one or a few of the parameters required by population modelling were lacking, we carried out simulations varying the above parameters within the range observed for closely related species and analysing variation in the results of the simulation procedure.
All simulations were conducted using the environmental variation coefficient (EV) or the standard deviation reported for a given species. When a specific EV or standard deviation was not available, we adopted an EV equal to 20% of the value of the parameter considered. When simulations were run on the less favourable scenarios, variation was set at zero. When reproductive parameters were entered in the form of productivity (number of fledglings per pair, irrespective of whether pairs were successful or unsuccessful), we considered 100% of adults as taking part in reproduction, with the only exception being species in which a proportion of adults do not participate in reproduction (i.e. do not try to breed, as is the case with several large raptor species). One type of catastrophe, with a frequency of 1% and associated with a 50% reduction in both reproductive outputs and survival, was entered into all simulations, except where stated. All parameter values have been specified in the species factsheets.
For abundant, non-colonial species with fewer than 2,500 breeding pairs, with no evidence for subdivision of populations into discrete units, we considered available breeding density data. When optimal habitat patches and related breeding density were known, we reported them as ‘optimal’ density, thus empirically obtaining or approaching a maximum value. At the largest scale, habitat quality is likely to vary spatially, and we aimed to obtain a ‘favourable’ density. Reference values for breeding density were always expressed in the unit area of the original data (e.g. if densities for a given species were available as the number of pairs km−2, the FRV was similarly expressed as pairs km−2). We mostly focused on breeding densities reported for Italy, but also considered data from other European countries.
Methods and case studies were then shown to experienced conservation biologists and ornithologists to get an expert-based revision and validation of both general procedures (methodological issues) and specific results (individual species outputs).
Here, we report the outcome of the process and present three different examples to show how FRVs were calculated in practical terms by using population modelling techniques and known estimates of breeding density.
Results
Outcomes of the process
Of the 88 Italian breeding species included in the Wild Birds Directive, one (European Dotterel Charadrius morinellus) was excluded because it is not a regular breeder (Brichetti and Fracasso Reference Brichetti and Fracasso2004). We calculated an FRV based on population modelling analyses for 47 populations (with fewer than 2,500 pairs) belonging to 21 species. The FRV was higher than current population estimates for all these species and populations, except (i) populations of Golden Eagle Aquila chrysaetos in areas where the species is believed to have reached carrying capacity and (ii) all main populations of Peregrine Falco peregrinus, which currently appear to be levelling out around the likely carrying capacity, and the species is colonising urban habitats (Allavena and Brunelli Reference Allavena and Brunelli2003, Brunelli Reference Brunelli, Magrini, Perna and Scotti2007).
For a further 15 species we provided FRVs in terms of breeding density at one or more scales, mainly depending on habitat patchiness and territory size. All FRVs formulated are reported in Appendix S1 in the online Supplementary Materials.
For two species occurring in Italy with < 2,500 pairs (Honey Buzzard Pernis apivorus and Short-toed Eagle Circaetus gallicus), data on reproductive and/or demographic parameters were scarce or lacking, while data on breeding density were fairly abundant. Therefore, we suggested favourable breeding densities, but were not able to formulate an FRV.
For 12 species that recently colonised or recolonised Italy, some after reintroduction (first nesting attempt no more than 30 years ago) we did not provide an FRV. The same applied to recently established populations of other species (e.g. reintroduced populations of Griffon Vulture Gyps fulvus in mainland Italy and Sicily). However, in four cases (three populations of Griffon Vulture and the Alpine population of Bearded Vulture Gypaetus barbatus) we provided conservation targets for reintroduced populations, defined as MVPs based on currently measured values of main demographic and reproductive parameters. We also did not provide an FRV for three colonial breeders with all populations or colonies > 2,500 breeding pairs: Cory’s Shearwater Calonectris diomedea, Black-crowned Night Heron Nycticorax nycticorax and Little Egret Egretta garzetta.
The remaining 35 species all have at least one population for which an FRV could be calculated (through population modelling or breeding density), but calculation was prevented by: a significant and largely unquantified exchange of individuals with extra-national populations (Audouin’s Gull Larus audouinii); unsustainable current mortality rates (Griffon Vulture in Sardinia), or a lack of adequate data (all others). Lack of data particularly affected certain families, especially the Rallidae (crakes: four species belonging to the genera Porzana, Crex and Porphyrio) and Picidae (woodpeckers: five species belonging to the genera Picus, Dendrocopos, Dryocopus, Picoides). The number of species falling into each category is summarised in Figure 1.
Procedure followed to define FRVs according to species characteristics and outcomes of the process.
Worked examples
Population modelling techniques: Egyptian Vulture Neophron percnopterus and Eleonora’s Falcon Falco eleonorae
Here we provide two examples of an FRV calculated through PVA:
One species with a single population for which a single FRV is given (Egyptian Vulture) and one species with different FRVs for each discrete population (Eleonora’s Falcon).
The Egyptian Vulture is a rare raptor, which in Europe mainly inhabits Mediterranean countries, especially Spain. Its current distribution in Italy essentially consists of only one population in Sicily (Brichetti and Fracasso Reference Brichetti and Fracasso2003), estimated at only seven pairs (Ceccolini and Cenerini Reference Ceccolini, Cenerini, Bellini, Cillo, Giacoia and Gustin2008). Egyptian Vulture is clearly one of the most endangered breeding species in Italy, being close to extinction, and thus belongs to the third group of species or populations treated under a PVA approach, i.e. those with a high probability of extinction in the next 100 years. Therefore, we calculated the MVP (extinction probability 1% in the next 100 years) using the average demographic and breeding parameters known for the species and took that as the FRV. We carried out different simulations, incorporating mortality and breeding parameters from different studies, in order (i) to account for the contrasting data reported for that species, and (ii) to show how the long-term viability of the target population is strongly dependent on both adult survival and breeding success. We developed two different scenarios: the first considered mortality values taken from a stable population (Donázar et al. Reference Donázar, Palacios, Gangoso, Ceballos, González and Hiraldo2002) and the second took values calculated for a steeply declining population (Inigo et al. Reference Inigo, Barov, Orhun and Gallo-Orsi2008). In the first scenario, we used the following parameter values: age of first reproduction five years; maximum life expectancy 20 years; productivity 0.97 ± 0.25 (Sarà and Di Vittorio Reference Sarà and Di Vittorio2003); proportion of breeding females equal to 0.9 (we assumed that not all females mate each year because of the very limited number of individuals); mortality equal to 50% in the first year of life, 40% in the second, 30% in the third, 10% in the fourth, 5% in the fifth, and 3% for adults, deriving values from Donázar et al. (Reference Donázar, Palacios, Gangoso, Ceballos, González and Hiraldo2002) but using adult mortality of 3% instead of < 2% (too low for the declining Italian population). Carrying capacity was set to mirror the Sicilian breeding population of the 1970s (40 pairs, 132 individuals; Ceccolini et al. Reference Ceccolini, Cenerini, Sarà, Fraissinet and Petretti2006). In the second scenario, we modified parameter values in accordance with recent Spanish data (Inigo et al. Reference Inigo, Barov, Orhun and Gallo-Orsi2008): age of first breeding of six years, mortality equal to 27 ± 2% in the first and second years of life, 22 ± 3% in the third, fourth and fifth years, 40 ± 5% in the sixth year, and 17 ± 2% in adults. In this scenario, we did not include catastrophes.
First scenario. The current population has a probability of extinction around 90% in the next 100 years. A population equal to carrying capacity showed a quite high extinction risk (extinction probability P ~ 17% in 100 years) and a tendency to decrease. Keeping constant mortality values, a population equal to the carrying capacity (132 individuals) has an extinction probability P ~ 1% with an increase in breeding success, which should be equal to 78% of adults breeding successfully, with a fledging rate of 1.4 young per successful pair. Similar breeding parameter values fall well within the range of values reported from other European countries (breeding success between 74.1% and 96.5% and productivity between 0.8 and 1.5; Braillon Reference Braillon1979, Reference Braillon1987, Bergier and Cheylan Reference Bergier and Cheylan1980, Marco and Garcia Reference Marco and Garcia1981, Bergier Reference Bergier1985, Gallardo et al. Reference Gallardo, Astruy, Cochet, Seriot, Neri, Torre and Thibault1987, Vasconcelos Reference Vasconcelos1987, Donázar and Ceballos Reference Donázar and Ceballos1988, Abuladze and Shergalin Reference Abuladze, Shergalin, Chancellor, Meyburg and Ferrer1998). We assumed that similar values could also be reached in Italy through protection of breeding sites and food supply for pairs experiencing low food availability. Therefore, we proposed an FRV (subject to adult mortality of 3%) of 40 pairs, with breeding success of 78% and fledging rate of 1.4.
Second scenario. Because of the extremely high adult mortality (17% instead of 3%), the risk of extinction further increased, and the current population would become extinct within the next 40 years, even without catastrophes or inbreeding effects and with an initial population of individuals all between five and nine years of age (a particularly favourable situation). However, with such low adult survival it is virtually impossible to maintain a viable population, even with higher breeding success. With adult mortality equal to 6 ± 2% (intermediate between the values reported by Donázar et al. Reference Donázar, Palacios, Gangoso, Ceballos, González and Hiraldo2002, and Inigo et al. Reference Inigo, Barov, Orhun and Gallo-Orsi2008, but closer to the former) and keeping other parameters constant, the long term survival probability becomes higher and population trend becomes stable. In short, barring catastrophes and inbreeding depression, a population of 132 individuals may be viable in the long term (extinction probability = 1% in 100 years) with a breeding success of 80%, fledging rate equal to 1.4, adult mortality of 6% (and mortality for other age classes as in the second scenario). With catastrophes as default in the other simulation, the probability of extinction of this hypothetical population would increase to 15%. The simulations carried out under the second scenario clearly show the need to keep adult mortality as low as possible in Egyptian Vultures, which have low annual productivity but a potentially long reproductive phase (Inigo et al. Reference Inigo, Barov, Orhun and Gallo-Orsi2008).
The Egyptian Vulture example shows how the calculation of an FRV for a single population could also be based on different scenarios in order to consider the importance of improving different parameters, such as breeding success instead of mortality, to make it possible for a population currently threatened with extinction to reach carrying capacity and to survive in the long term.
Eleonora’s Falcon is an endemic Mediterranean bird of prey that inhabits rocky islands and islets (Cramp and Simmons Reference Cramp and Simmons1977). The current distribution of the species in Italy essentially includes two populations, the larger one inhabiting Sardinia and adjacent smaller islands, and the other occupying Sicily and especially its satellite islands, the Eolian and Pelagie islands (Brichetti and Fracasso Reference Brichetti and Fracasso2003). The Sardinian population is currently estimated at about 500 pairs, the Sicilian population at 138–204 pairs (Gustin et al. Reference Gustin, Corso and Medda2005). The Sicilian population appears to have declined in recent years (Corso Reference Corso2008), while the apparently positive population trend shown by the species in Sardinia is arguably due to improved knowledge of the distribution and size of the colonies (Gustin et al. Reference Gustin, Corso and Medda2005). Therefore, Eleonora’s Falcon populations are considered to have an unfavourable (Sicilian) or unknown (Sardinian) trend and the main parameter values selected were among the less favourable. Since there is no data on carrying capacity, this is prudently assumed to be equal to the initial tested population in both cases. The only data on demographic parameters were reported by Ristow et al. (Reference Ristow, Scharlau, Wink, Meyburg and Chancellor1989) and are as follows: mortality 78% before adulthood, 13% in adulthood (two years); maximum longevity at least 16 years; first breeding generally takes place at two years of age for females and three years for males, but exceptions are numerous.
Cramp and Simmons (Reference Cramp and Simmons1977) reported first breeding at two years of age, without differentiating between the sexes. Exceptionally, one-year-old individuals can breed (Ristow et al. Reference Ristow, Scharlau, Wink, Meyburg and Chancellor1989). We thus decided to treat individuals of two years of age or older as breeders, following the analysis by Ristow et al. (Reference Ristow, Scharlau, Wink, Meyburg and Chancellor1989) of mortality and demographic structure in this species. We thus estimate the following mortality rates per year: 65% during the first year of life, 37% during the second (amounting to an overall mortality of 78% before breeding) and 13% from the third year on. Considering all known productivity values in Italy, the mean productivity is equal to 1.26 (± 0.39 SD) fledglings per pair (Badami Reference Badami1995a,Reference Badamib, Brichetti and Fracasso Reference Badami2003, Gustin et al. Reference Brichetti and Fracasso2005, Medda Reference Gustin, Corso and Medda2006, Spina and Leonardi Reference Medda2007). The Sardinian population (500 pairs, equivalent to 1,250 individuals under a stable-age distribution), which has good chances for long-term survival given current average productivity, would face a high risk of extinction should its productivity fall to 1.1 or less, something which has happened more than once in the past (Gustin et al. Reference Gustin, Corso and Medda2005, Medda Reference Medda2006, Corso Reference Corso2008). In this case, the MVP (P = 1%) would be 2,200 individuals, or 900 pairs. This figure, which would guarantee survival even under scenarios slightly less favourable than the current one, is therefore proposed as the FRV for the Sardinian population.
The Sicilian population (average estimate 176 pairs, or about 425 individuals) shows a relatively high (> 10%) probability of extinction over the next 100 years; the MVP (P = 1%) calculated using the figure for average productivity recorded in Italy would be 800 individuals or about 320 pairs; we thus propose 320 pairs as the FRV for the Sicilian population, on condition that average adult productivity is not below 1.26. Ristow and Wink (Reference Ristow, Wink, Newton and Chancellor1985) reported that a productivity of 1.2 fledged young for each breeding attempt was necessary to maintain a colony, a figure that accords with the results of this analysis.
The example of Eleonora’s Falcon shows how calculation of the FRV for different populations is based on different scenarios in order to reflect different conservation objectives: reaching a population size able to survive even in worsened conditions (Sardinia), or a population size with reasonably good chances of long-term persistence, in the case of smaller, and currently threatened, populations (Sicily).
Known estimates of breeding density:
Woodlark Lullula arborea and Red-backed Shrike Lanius collurio. Woodlark is a quite widespread but rarely abundant passerine species endemic to Europe. It has a rather large population in Italy (estimated at 20,000–40,000 pairs; Brichetti and Fracasso Reference Brichetti and Fracasso2007) which lacks appreciable subdivision into discrete units or populations, although the decline of the Alpine population is now increasing the gap between this and the other parts of the Italian range. However, the ecology of this species is quite similar throughout its range and densities appear similar in different areas. Therefore, the FRV is formulated in terms of breeding density at two spatial scales. Considering the highest known densities, optimal Woodlark density at large scales can be identified as 10 pairs km−2 and at a local scale equal to three pairs 10 ha−1 (Brambilla and Rubolini Reference Brambilla and Rubolini2009). Ten pairs km−2 and three pairs per 10 ha are thus proposed as FRVs at ‘large’ and ‘local’ scales, respectively.
Red-backed Shrike is a typical passerine species of Eurasian low-intensity farmland landscapes. It underwent large population declines over most of Europe during the last century (Heath Reference Heath, Tucker and Heath1994). The Italian population is estimated at 50,000–120,000 pairs (BirdLife International 2004). There have been many studies of its breeding ecology (see Brambilla et al. Reference Brambilla, Casale, Bergero, Crovetto, Falco, Negri, Siccardi and Bogliani2009 and references therein) and thus estimates of breeding densities for Italy and elsewhere are numerous and relate to different habitats, from highly suitable landscapes to patchy habitats with lower average suitability (Casale and Brambilla Reference Casale and Brambilla2009). Reviewing the data, we found that in areas with mosaics of open and semi-open habitats, breeding densities could be around one pair km−2 at large scales, but in mostly open habitats (large pastures, low-intensity farming systems, etc.) densities up to five pairs km−2 are frequently reached (Brambilla et al. Reference Brambilla, Rubolini and Guidali2007, Reference Brambilla, Casale, Bergero, Crovetto, Falco, Negri, Siccardi and Bogliani2009). At the local scale (< 100 ha), in mosaic areas (with scattered patches of suitable open or semi-open habitats), density should be no lower than 0.5 pairs 10 ha−1 (Brambilla and Casale Reference Brambilla and Casale2008). Densities in suitable habitats are 3–5 pairs 10 ha−1 and reach 8–10 pairs 10 ha−1 in exceptionally suitable areas (Brambilla et al. Reference Brambilla, Casale, Bergero, Crovetto, Falco, Negri, Siccardi and Bogliani2009). Therefore, we proposed an FRV value of one pair km−2 at ‘large’ scales (for mosaic areas with scattered suitable habitats; five pairs km−2 for open or semi-open areas) and 0.5 pairs 10 ha−1 at ‘local’ scales (for mosaic areas; five pairs 10 ha−1 for suitable habitats; 10 pairs 10 ha−1 for exceptionally suitable areas).
Discussion
When defining FRVs for breeding bird species in Italy, we assigned species to multiple categories, reflecting differences in abundance, distribution and demographic trends.
For species that colonised Italy in the last 30 years (or recolonised spontaneously or were reintroduced), we did not provide an FRV. We believe that the dynamics of population expansion following colonisation may lead to unpredictable outcomes for most of those species, especially in the absence of reliable estimates of the carrying capacity of the habitats. However, it was possible to calculate an MVP for some species using current estimates of demographic and reproductive parameters and such a value was proposed as the mid-term conservation target; this procedure was applied to Bearded Vulture and to three reintroduced populations of Griffon Vulture. Setting concrete targets means translating reference values into operational, practical and feasible short-, medium- and long-term targets. For all other species, an FRV could be obtained with adequate knowledge of the parameters. Filling the gap in knowledge that affects quite a large number of families (especially some poorly studied ones, such as the Rallidae and Picidae) truly represents a conservation priority, as it is impossible to set reference values for such poorly known species. In fact, demographic parameters were the most frequently lacking data, but reproductive data were also unknown for a large number of species.
The use of density-based FRVs allowed us to by-pass the impossibility of obtaining reliable estimates of population size for most of the abundant species and the weakness of PVAs with a large number of individuals. A simple comparison between ground-measured density values at a site or group of sites and the reference value proposed for target species should allow site managers to make an immediate rough assessment of the conservation status of target bird species. The use of density values (especially for passerines) may hopefully become a practical tool for evaluating population status and for deriving consequent implications under management plans of Natura 2000 sites and other protected areas.
Conclusions
Our work aimed to establish long-term conservation targets in the form of ‘Favourable Reference Values’, i.e. demographic values which should guarantee species survival in the long term, for bird species of particular conservation interest at national and EU levels. Obviously, maintaining populations of adequate size is not independent of habitat and general environmental conditions, which should be maintained to the necessary extent and quality needed by target species. Ideally, formulation of FRVs for range and habitat should accompany the population reference value, but at present, relevant data are lacking for most bird species.
Unfortunately, lack of adequate knowledge prevented the definition of FRVs for some of the target species. On the other hand, our work indicated the principal areas where further research is necessary.
Additional work should focus on translating density values into population estimates on the basis of the current, historic, and/or potential range of those species; in order to do so, however, potentially suitable habitat should be quantified, but these data are currently lacking. A further potential problem related to the use of breeding density as a reference value is that sites functioning as ecological traps may host a high density of a species even though they act as population sinks (Battin Reference Battin2004). Data on reproductive success may demonstrate a real correlation between breeding density and local status/dynamics. However, breeding density is widely adopted as an indirect measure of habitat suitability and it could be taken as a first step in measuring conservation status for our general purposes, aimed at monitoring species at a large scale.
Despite these limitations, our approach provides a useful framework for defining a population reference value for species of varying conservation concern, accounting for their demographic history, life-history characteristics and distributional features.
Acknowledgements
We are particularly grateful to all the people who contributed to this work, and especially to scientists who revised our general approach to FRV, both in Italy (N. Baccetti, G. Boano, G. Bogliani, P. Brichetti, P. Pedrini, D. Rubolini, F. Spina) and at BirdLife International (B. Barov, I. Burfield). An anonymous reviewer provided very useful suggestions on a first draft. We are very grateful also to V. Bergero, E. Bianchi, G. Bottari, F. Casale, P. Colombani, E. Duprè, R. Falco, F. Piccarolo, M. Ravasini for kind cooperation. LIPU received funding for this project from MATTM (Ministry of the Environment, Territorial Protection, and the Sea); ISPRA and FLA participated as project partners.
