Next Article in Journal
New Insights in Canine Reproduction
Next Article in Special Issue
Bobolink (Dolichonyx oryzivorus) Declines Follow Bison (Bison bison) Reintroduction on Private Conservation Grasslands
Previous Article in Journal
Age-Related Changes in the Primary Motor Cortex of Newborn to Adult Domestic Pig Sus scrofa domesticus
Previous Article in Special Issue
Is Hay for the Birds? Investigating Landowner Willingness to Time Hay Harvests for Grassland Bird Conservation
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Animals
Volume 11
Issue 7
10.3390/ani11072020
animals-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Caution Is Needed When Using Niche Models to Infer Changes in Species Abundance: The Case of Two Sympatric Raptor Populations

by
Adrián Regos
1,2,*,
Luis Tapia
1,3,
Alberto Gil-Carrera
3,4 and
Jesús Domínguez
1
1
Departamento de Zooloxía, Xenética e Antropoloxía Física, Universidade de Santiago de Compostela, Campus Sur, 15782 Santiago de Compostela, Spain
2
CIBIO/InBIO, Research Center in Biodiversity and Genetic Resources, Predictive Ecology Group, Campus Agrario de Vairão, R. Padre Armando Quintas, N° 7, 4485-661 Vairão, Portugal
3
EBX, Estación Biolóxica do Xurés, Vilameá 121, 32870 Lobios, Spain
4
GREFA, Monte del Pilar S/N, 28220 Majadahonda, Spain
*
Author to whom correspondence should be addressed.
Animals 2021, 11(7), 2020; https://doi.org/10.3390/ani11072020
Submission received: 2 June 2021 / Revised: 23 June 2021 / Accepted: 1 July 2021 / Published: 6 July 2021
(This article belongs to the Special Issue Birds in Agricultural Environments)
Browse Figures
Review Reports Versions Notes

Abstract

:

Simple Summary

This study focuses on sympatric populations of hen harrier (Circus cyaneus) and Montagu’s harrier (Circus pygargus) breeding in mountain heathlands in the NW Iberian Peninsula. These populations have been affected by habitat loss caused by land-use change. Despite the mounting evidence supporting positive relationships between species abundance and habitat suitability, the capacity of ecological niche models to capture variations in population abundance remains largely unexplored. This study shows that habitat suitability predicted from niche models was significantly correlated with the relative abundance for hen harrier and, to a lesser extent, for Montagu’s harrier. However, the temporal variation in local population abundance was not significantly explained by habitat suitability changes predicted by the niche models. These findings call for caution in the use of niche models to infer changes in population abundance. The positive relationship between species abundance and habitat suitability supports the use of niche models to estimate abundance but does not guarantee the ability of these models to predict temporal variations in species abundance. These findings highlight (1) the need to include other possible biotic or abiotic factors involved in population abundance dynamics into niche models and (2) the need to establish specific monitoring protocols for tracking population dynamics.

Abstract

Despite the mounting evidence supporting positive relationships between species abundance and habitat suitability, the capacity of ecological niche models (ENMs) to capture variations in population abundance remains largely unexplored. This study focuses on sympatric populations of hen harrier (Circus cyaneus) and Montagu’s harrier (Circus pygargus), surveyed in 1997 and 2017 in an upland moor area in northwestern Spain. The ENMs performed very well for both species (with area under the ROC curve and true skill statistic values of up to 0.9 and 0.75). The presence of both species was mainly correlated with heathlands, although the normalized difference water index derived from Landsat images was the most important for hen harrier, indicating a greater preference of this species for wet heaths and peat bogs. The findings showed that ENM-derived habitat suitability was significantly correlated with the species abundance, thus reinforcing the use of ENMs as a proxy for species abundance. However, the temporal variation in species abundance was not significantly explained by changes in habitat suitability predicted by the ENMs, indicating the need for caution when using these types of models to infer changes in population abundance.

1. Introduction

Habitat loss is the most important threat to biodiversity worldwide and strongly impacts the distribution and abundance of both animal and plant species [1]. “Species distribution” and “population abundance” are considered essential biodiversity variables (EBVs) by the Group on Earth Observations Biodiversity Observation Network (GEO BON) [2]. Long-term time series of population abundance (obtained from standardized monitoring protocols involving repeated surveys of the same population) are essential for population management but very costly and time-consuming [3].
Model-assisted monitoring has recently emerged as an alternative, cost-effective way of assessing the impact of environmental change on biodiversity [4]. Ecological niche models (also known as habitat suitability or species distribution models [5]), which correlate species occurrence data and current environmental covariates (mostly climatic and topographic variables), have been widely used by ecologists in recent decades [6]. Species occurrence data are far more common than species abundance data, which require more costly and time-consuming sampling methods. Consequently, abundance data are often substituted by environmental suitability predicted by ecological niche models (ENMs), by assuming positive relationships between species abundance and habitat suitability [7,8]. The use of ENMs to produce data that can be used as a proxy for species abundance has been supported by many studies over the last years (see [9] and references therein). Despite the mounting evidence supporting such relationships, other studies have also reported weak (or no) correlation between species abundance and environmental suitability predicted by ENMs [10], placing into question the capacity of such models to capture the variation in species abundance. A common limitation of niche models is that they often relied on climatic envelopes inferred from macroclimatic datasets [9], therefore resulting in predictions of habitat suitability geographically overestimated. This limitation highlights the need to investigate nonclimatic sources of population abundance variation at local scales [10].
According to the “Habitat Selection Theory”, individuals will select habitats that maximize their biological “fitness”. The importance of local habitat conditions has been corroborated in several studies in which variables other than climate-related variables are used to determine environmental suitability [10]. In this regard, satellite remote sensing (SRS) provides continuous and standardized measures of the environment at high spatial resolution, allowing cost-effective and reliable characterization of species habitats at local scales [11,12]. Thus, SRS data have increasingly been used as predictor variables for both species distribution and abundance models in recent years [13,14]. For instance, the incorporation of SRS variables related to water and carbon cycles (e.g., the normalized difference water and vegetation indices (NDWI and NDVI, respectively)) has recently been reported to be critical for model-assisted monitoring of endangered plant and animal species [13,15]. These studies suggest that calibrating ENMs with SRS variables would enable the accurate prediction of species distributions at local scales (across space and time). In light of the habitat selection theory, we would expect model-predicted environmental suitability to be significantly correlated with species abundance, thus supporting the use of SRS-based ENMs to estimate species abundance. If so, both the spatial and temporal variation in species abundance should be explained (at least partly) by changes in environmental suitability predicted by these models.
This study focuses on sympatric populations of hen harrier (Circus cyaneus) and Montagu’s harrier (Circus pygargus) breeding in mountain heathland (Natura 2000 habitat type: 4020 European wet heaths, 4030 European dry heaths and 7110 active raised bogs listed in Annex I of the Habitat Directive 92/43/EEC) in the NW Iberian Peninsula [16,17]. Both species are included in the Annex I of the European Bird Directive with conservation status “Vulnerable” in Spain (see [18]). These species are considered good indicators of environmental change due to their high susceptibility to being affected by habitat loss and fragmentation caused by land-use change [19]. In fact, the populations have declined sharply in recent years in Spain, especially in the northwest of the country, due to intensive agroforestry and land-use change [18,20,21]. In particular, the habitat availability and quality of both target species have markedly decreased in our region over the last 20 years due to large-scale land-use conversions from mountain heathlands and grasslands to fast-growing tree plantations (mostly Pinus and Eucalyptus spp.) and native forest expansion caused by intensive forest practices and rural abandonment processes (i.e., loss of traditional agricultural and livestock practices) [14,20,21].
To test the above-mentioned hypotheses, we assessed the predictive accuracy of ENMs (i.e., realized environmental space sensu “Grinnellian niche” [22,23,24]) based exclusively on SRS variables (derived from Landsat satellite imagery) and occurrence records (in situ data) and quantified the extent and direction of changes in environmental suitability predicted by these models for each species over the last 20 years (between 1997 and 2017). Landscape change analysis was also performed using the same Landsat satellite images to clarify the underlying processes affecting species habitat changes. Finally, we tested whether the environmental suitability predicted by the ENMs is significantly correlated with species abundance and thus supports the use of RS-based ENMs for estimating abundance and whether the spatiotemporal variation in species abundance is explained by changes in environmental suitability predicted by these models.

2. Materials and Methods

2.1. Study Area

The study area is located in the “Serra do Laboreiro”, which is included in the “Baixa Limia” Site of Community Importance (SCI) and which is a Special Protection Area for birds (SPA) (Galicia, NW Spain) (Figure 1). The study site covers an area of 3350 ha, bordering to the north with the Peneda-Gerês National Park (North of Portugal). It is a mountain range, with an average elevation of around 800 m a.s.l., a gradual, gentle relief and summits of up to 1300 m, predominantly comprised of granite rock. The climate in the area is temperate sub-Mediterranean oceanic, with an annual mean temperature of 8–12 °C and mean annual precipitation of 1200–1600 mm [25]. The landscape is representative of European Atlantic mountain heathlands, dominated by upland moors and open vegetation communities, mainly European heaths (Erica sp., Chamaespartium tridentatum, Ulex sp., Cytisus sp., among others), with a strong influence of extensive livestock grazing [26]. Mountain heathlands are nowadays under threat and very restricted in southern Europe, being priority conservation areas in Spain [27,28]. Forests are fragmented and dominated by pine (Pinus sylvestris) in the highest areas and oaks (Quercus robur, Q. pyrenaica) at lower elevations and in valleys [29]. The human population in the area is currently very sparse, although the landscape has been intensely managed by humans for centuries and is especially associated with the presence of extensive livestock. Rural communities still use fire to create pastures for extensive livestock, which leads to frequent intentional wildfires in the area, mainly at the end of winter.

2.2. Field Surveys and Target Species

The hen harrier and Montagu’s harrier populations were surveyed every month during the 1997 breeding season, between 15 April and 31 August. Despite the migratory behavior of the Montagu’s harrier and dispersive movements of the hen harrier, these populations are very philopatric (i.e., individuals tend to return to breed in this area every year) so no individual was recorded in migration. In our study, the fledglings stayed in the breeding area until September (no migratory restlessness was observed). The survey was repeated in 2017 with the same timing as in the 1997 census to minimize any phenological bias. In each breeding season, two road “point transect censuses” [30] were carried out from a 4 × 4 vehicle, taking advantage of the vast unpaved track network in the study area (see Figure 1). Both road censuses were 9 km long, and point counts were established at intervals of 1 km to avoid pseudoreplication, with a total of 10 per transect. All visual sightings of individual birds were recorded during 10 min in each point count. In each transect, every individual was controlled to avoid double counts. Two experts in raptor identification carried out all censuses in both years. The precise locations of observations were recorded via a Garmin GPS handheld receiver (accuracy of 15 m ± 3 m). Finally, 20 point counts were conducted every month (n = 100 per year), yielding a total of 1000 min of effective sampling effort per year. The surveys revealed a decrease in the total number of hen harrier recorded between 1997 and 2017 (n1997 = 13; n2017 = 6), but no changes in numbers of Montagu’s harrier (n1997 = 11; n2017 = 11).

2.3. Satellite-Based Land Cover Mapping

We used freely available Earth observation (EO) imagery to derive land use and land cover (LULC) maps for 1997 and 2017. The main EO data source consisted of optical multispectral bands (30 m resolution) from two Landsat images acquired on dates as close as possible to the field survey. More specifically, we used spectral bands from Landsat 5 TM (5 April 1997) and Landsat 8 OLI sensors (28 April 2017). All scenes were downloaded from the United States Geological Survey (USGS) Global Visualization Viewer (http://glovis.usgs.gov, accessed on 1 November 2020). All Landsat scenes were processed to Standard Terrain Correction (Level 1T), which provides systematic geometric accuracy by incorporating ground control points while applying a digital elevation model (DEM) for topographic accuracy. Digital numbers (DNs) were converted to top-of-atmosphere radiance and physically meaningful units by radiometric calibration and application of sensor- and band-specific calibration parameters. The classification process was based on the radiometric information obtained from reflective bands and two multispectral indices for each image: (1) the normalized difference vegetation index (NDVI [31]) and (2) the normalized difference water index (NDWI [32]), to enhance the contribution of vegetation in the spectral response and mitigate other factors such as soil, topography, lighting conditions and atmosphere [31]. Topographic information was obtained from NASA’s Shuttle Radar Topography Mission (GDEM-SRTM) (https://www2.jpl.nasa.gov/srtm/ accessed on 1 November 2020).
Supervised classification was carried out using the following four classification algorithms available in the R-based package “Caret” and implemented in the “RStoolbox” package, version 0.1.534: (1) stochastic gradient boosting, (2) random forest, (3) support vector machines (SVMs) with linear kernel and (4) SVMs with radial basis function kernel. The final land use/cover maps were elaborated by applying an ensemble approach, including the simple voting system (also known as “majority voting” or “select all majority” system, sensu [33]) that includes those with an overall accuracy higher than 95% (see Table S1). The four most important land use and vegetation cover types identified during the fieldwork were: (1) farmland (extensive agricultural areas interspersed with native deciduous forest, mostly oak woodlands), (2) open heathland (open vegetation areas dominated by dry heaths and largely affected by burning and extensive livestock rearing), (3) closed heathland (closed/mature heath areas dominated by dry and wet heaths, including peat bogs) and (4) pine plantations. Training and validation areas for each LULC class were established by on-screen digitizing, with QGIS software, and consisted of a set of pixels identified over well-known homogeneous areas in each Landsat image, thus providing a reference spectral signature for each class. Once the ensemble LULC maps were constructed, net changes were quantified for each land cover type for the entire study period (1997–2017) in the area, including 500 m buffers around each point count (see Table 1).
Data importation, preprocessing, computation of spectral indices and image classification were performed using the RStoolbox package, version 0.1.534. Land-use/cover change analysis was performed with the “lulcc” v.1.0.2 package in R [34].

2.4. Ecological Niche Modelling

To quantify changes in habitat suitability at the landscape level for both species between 1997 and 2017, we first constructed ecological niche models (ENMs; i.e., realized environmental space sensu “Grinnellian niche” [22,23,24]). These ENMs allowed us to empirically correlate hen harrier occurrence data (presence and absence data registered by point count for each year of sampling) and the percentage of each LULC class (pine plantations, farmland, open and closed heathlands), topography (altitude, slope and aspect) and water content of vegetation (computed from the NDWI) within a radius of 500 m around each point count. We applied eight widely used modeling algorithms implemented in BIOMOD2 to deal with the uncertainty of different modeling techniques and to provide more informative and ecologically correct predictions [35]: (1) generalized linear models, (2) flexible discriminant analysis, (3) classification tree analysis, (4) multivariate adaptive regression splines, (5) maximum entropy, (6) random forest, (7) artificial neural networks and (8) surface range envelope [36]. We then tested the performance of the different models by splitting the original raptor data set into two subsets: 70% of the data was used to train the models and the remaining 30% was used to validate the model. We randomly repeated this procedure 30 times to produce predictions independent of the training data. We applied a weighted average approach to compute a consensus of single-model projections by using the area under the receiver operating characteristic (ROC) curve (AUC) to weight the model [37]. Only models with AUC values above 0.8 were used in the ensemble procedure. The ensemble models, obtained from the consensus of single models, were directly projected at 60 m resolution [38]. We calculated two different evaluation indices: the area under the ROC curve (AUC) [39] and true skill statistic (TSS) [40].

2.5. Testing the Species Abundance–Habitat Suitability Relationship

To test the hypotheses about the relationship between species abundance and habitat suitability (see Section 1), we related the relative abundance index (estimated as the number of individuals recorded in each point count and sampling year) to the habitat suitability (hereafter HS) predicted by the ENMs for each species and year. For this purpose, we fitted generalized linear models with a negative binomial distribution of errors and logarithmic link function with the function “glm.nb” available in the R package “MASS”, version 7.3–51.1. We also explored the effect of habitat changes predicted by ENMs between 1997 and 2017 on the variation in species abundance (estimated as the difference in the abundance index between 1997 and 2017) by fitting linear models. In addition, we used Pearson’s correlation coefficient to test for any association between species abundance and HS, and we used Spearman’s rho to test for any relationship between variation in species abundance and change in HS. For the sake of caution, only effects with p-values < 0.01 were considered significant.

3. Results

Overall, the ensemble ENMs performed very well for both species, as measured by AUC (AUC hen harrier = 0.95 and AUC Montagu’s harrier = 0.93) and TSS metrics (TSS hen harrier = 0.85 and TSS Montagu’s harrier = 0.75). In relation to the importance of the environmental variables included in niche models, the results reflect the ecological requirements of both species (Figure 2). As expected, the occurrence of both species was mainly correlated with open and closed heathlands (hunting and breeding areas, respectively), although the NDWI was only found to be important for hen harrier, indicating its preference for areas with damp vegetation (wet heaths and peat bogs).
The types of land-use change detected have led to contrasting changes in the habitat availability for both species over the last 20 years (Table 2, Figure 3). Thus, the ENMs predicted a marked increase in the habitat availability for Montagu’s harrier and a significant reduction for hen harrier (Table 2, Figure 3).
Regarding the land-use changes within areas indicated by the ENMs to be suitable habitat in both years, the results showed no important changes, as these areas are closely associated with open heathland for both species but also with farmland in the case of hen harrier (see “Overlap” in Figure 4). Areas where raptor habitat increased between 1997 and 2017 were associated with important conversion from closed to open heathland for both raptor species and to farmland for hen harrier (see “Gains” in Figure 4). Finally, areas with predicted habitat loss were linked to extensive conversion from open to closed heathland and pine plantations (see “Losses” in Figure 4).
A significant positive correlation between abundance and habitat suitability indices was obtained for both species and years (Figure 5A). In particular, the results showed that habitat suitability predicted by the ENMs significantly explained the observed abundance of hen harrier (pseudoR2Nagelkerke = 0.373; p-value < 0.001) and, to a lesser extent, that of Montagu’s harrier (pseudoR2Nagelkerke = 0.271; p-value < 0.01). The results also showed that the greater the increase in habitat suitability, the greater the increase in observed abundance (Figure 5B). However, the positive relationships between species abundance and habitat suitability were not strong enough to explain changes in species abundance (R2 = 0.1379; p-value = 0.036).

4. Discussion

Our findings support the species abundance–habitat suitability hypothesis, reinforcing the use of ecological niche models (ENMs) as a proxy for species abundance. However, the temporal variation in species abundance was not significantly explained by changes in habitat suitability predicted by the ENMs, which calls for caution in the use of ENMs when inferring changes in population abundance.
Various authors have already highlighted the positive relationships between habitat suitability and abundance/distribution patterns [8,9]. As both distribution and abundance patterns, particularly for raptors, are determined by the availability of suitable nesting places and trophic resources, habitat characteristics should a priori be good predictors of abundance [41]. However, the use of abundance and distribution data to infer relationships between species and habitat characteristics should be considered with caution, as these data can misleadingly reflect habitat suitability and actual carrying capacity [42,43]. The concept of carrying capacity encapsulates the notion that habitat resources must ultimately limit the numbers of raptors in an area. Nevertheless, factors other than those purely related to habitat characteristics (e.g., human pressure, intra- and interspecific interactions, period of the annual cycle) are known to affect local abundance and population dynamics [44,45]. Although our findings showed that species abundance and habitat suitability are significantly correlated, they also indicated that changes in abundance between 1997 and 2017 are not exclusively determined by changes in habitat suitability. Therefore, other variables, not directly associated with habitat dynamics, may have affected the observed changes in abundance. For instance, direct persecution (illegal killing), human-related disturbance of nest and foraging areas, predation and the effect of pesticides are typical examples of relevant factors that cannot be easily incorporated in ENMs [46,47]. In addition, testing the relationship between abundance and habitat suitability requires a good estimation of the optimal environmental conditions for the species (see [48,49,50]). This implies the use of ENMs that are better at characterizing the fundamental niche of the species than those based exclusively on local habitat variables (i.e., more related to the realized niche of species). Moreover, modeling endangered species (i.e., species often with restricted distributions) is particularly challenging since ENMs are often built from few records and at local scales, which affects both model calibration and evaluation procedures [51,52,53]. Moreover, the semicolonial and migratory behaviors of species such as Montagu’s harrier or the large hunting distances of breeding male and female hen harriers [54,55] can reduce the ability of ENMs to infer changes in abundance since surveys are more prone to double counts or pseudoreplication than for other sedentary and colonial species. Another issue that can constrain our ability to infer population dynamics from ENMs is the interannual fluctuations in prey availability. In this regard, the inclusion of remotely sensed ecosystem functional variables (e.g., annual mean of NDVI or NDWI as a proxy for primary productivity [56], food resources [57] or prey availability [58]) into ENMs has been found to improve model predictions [13,15,58,59].
Despite these limitations, our findings suggest an important effect of habitat quality on the patterns of distribution and abundance for both raptor species. The ENMs provided good fits to the data for both species (with AUC and TSS values of up to 0.95 and 0.85, respectively), which allowed us to test with confidence whether ENM-predicted habitat suitability was correlated positively with species abundance. The ranked importance of the environmental variables included in ENMs satisfactorily reflected the ecological requirements of the target species in the study area [16,17]. As expected, both species were mainly associated with open and closed heathland (hunting and breeding areas, respectively). Optimal habitat conditions for these species are represented by a mosaic of heather and grassland [17]. Areas dominated by tall heather, particularly those furthest from forest harvesting areas, are the most valuable for hen harrier nesting/roosting. Sites with little or no heather (other than very wet blanket bogs) or completely dominated by tall woody heather may need management interventions [60]. More interestingly, our findings highlight the importance of vegetation water content (i.e., wet heathland and peat bogs) in explaining the distribution of hen harrier, which to our knowledge is one aspect of the ecological niche not documented to date in the Iberian Peninsula, although known in other more northern areas of Atlantic Europe [60,61]. This finding highlights the usefulness of remotely sensed descriptors of heathland water balance for characterizing the species niche and confirms the need to go beyond the traditional habitat description (i.e., often based on landscape structural and compositional attributes) [13,58].
Land-use change has led to changes in habitat availability for both species (Table 1 and Table 2). The ENM projections revealed a large decrease in the habitat available for hen harrier and an overall increase in suitable habitat for Montagu’s harrier. Areas where raptor habitat changed between 1997 and 2017 were associated with important changes from mature/closed to open heathland for both species and to farmland for hen harrier. Heath and shrub formations in the mountainous areas of NW Spain have historically been linked to extensive agropastoral activities. However, these extensively managed systems have gradually been abandoned and largely replaced by monoculture forest plantations (mostly of Scots pine Pinus sylvestris in the study area) and semi-intensive grassland [27,28,29]. Side-effects on these raptor populations are traditionally associated with the conservation status of wet heath and peat bog moorland in Atlantic Europe [60,61,62], but these habitats are scarce and relict in the context of the Iberian Peninsula [27,28]. Furthermore, our findings suggest that forest plantations for timber production have limited the foraging opportunities and reproductive fitness, with potential impacts at population scale [60]. Our findings also showed a large reduction in habitat availability for hen harrier associated with the loss of wet heathland and peat bogs in the study area (Figure S1).
The Common Agricultural Policy (CAP) and European Regional Development Funds (FEDER Funds) aim to maintain traditional agropastoral systems to ensure the conservation and sustainable management of these priority habitats. However, these EU policies are failing in regard to biodiversity and land degradation [63]. Urgent action is required to achieve a sustainable model for European agriculture that can jointly benefit biodiversity and rural communities in mountain heathlands [63,64].

5. Conclusions

Our ecological niche models performed very well for both raptor species, predicting a marked increase in the habitat availability for Montagu’s harrier and a significant reduction for hen harrier. In relation to the ranked importance of the environmental variables included in the niche models, our findings identified the ecological requirements for both species. The presence of both species was mainly correlated with open and closed heathlands, although the normalized difference water index was found to be the most important for hen harrier, indicating the preference of this species for wet heathland and peat bog areas. Our findings showed that habitat suitability predicted from niche models was significantly correlated with the relative abundance for hen harrier and, to a lesser extent, for Montagu’s harrier (due likely to its migratory and semicolonial behavior). However, the temporal variation in local population abundance was not significantly explained by habitat suitability changes predicted by the niche models. These findings call for caution in the use of niche models to infer changes in population abundance. The positive relationship between species abundance and habitat suitability supports the use of niche models to estimate abundance but does not guarantee the ability of these models to predict temporal variations in species abundance. These findings highlight (1) the need to include other possible biotic or abiotic factors involved in population abundance dynamics into niche models and (2) the need to establish specific monitoring protocols for tracking population dynamics. The inherent difficulties associated with the migratory and semicolonial behavior of the Montagu’s harrier together with the large hunting distances of breeding male and female hen harriers highlight the need for sampling protocols supported by marked individuals (e.g., with radio tracking) that allow tracking their habitat occupancy along the annual cycle or improving ENMs by including time series of remotely sensed variables as proxies for interannual habitat dynamics (e.g., annual variations in prey availability).

Supplementary Materials

The following are available online at https://www.mdpi.com/article/10.3390/ani11072020/s1, Figure S1: Normalized difference water index (NDWI) within each 1 km buffer area around each point count for years 1997 and 2017, Table S1: Overall accuracy (% of pixels well classified) and Kappa index of each land-cover map derived from the different Landsat image classification methods for each year.

Author Contributions

Conceptualization, L.T. and A.R.; methodology, L.T, J.D. and A.R.; raptor surveys, L.T., A.G.-C. and J.D.; remote sensing and SDM computation, A.R.; writing—original draft preparation, L.T. and A.R.; visualization, A.R.; funding acquisition, J.D. All authors have read and agreed to the published version of the manuscript.

Funding

This work received funding from Xunta de Galicia through the grant to structure and consolidate competitive research groups of Galicia (ED431B 2018/36).

Institutional Review Board Statement

Not applicable.

Data Availability Statement

Landsat images are freely available from the United States Geological Survey (USGS) Global Visualization Viewer (http://glovis.usgs.gov accessed on 1 November 2020). Topographic information was obtained from NASA’s Shuttle Radar Topography Mission (GDEM-SRTM) (https://www2.jpl.nasa.gov/srtm/ accessed on 1 November 2020). Raptor data are available from the corresponding author on reasonable request.

Acknowledgments

A.R. was funded by the Xunta de Galicia, Spain (postdoctoral fellowship ED481B2016/084-0).

Conflicts of Interest

The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript, or in the decision to publish the results.

References

  1. Díaz, S.; Settele, J.; Brondízio, E.; Ngo, H.T.; Guèze, M.; Agard, J.; Arneth, A.; Balvanera, P.; Brauman, K.A.; Butchart, S.H.M.; et al. Summary for Policymakers of the Global Assessment Report on Biodiversity and Ecosystem Services of the Intergovernmental Science-Policy Platform on Biodiversity and Ecosystem Services; IPBES secretariat: Bonn, Germany, 2019. [Google Scholar]
  2. Kissling, W.D.; Ahumada, J.A.; Bowser, A.; Fernandez, M.; Fernández, N.; García, E.A.; Guralnick, R.P.; Isaac, N.J.B.; Kelling, S.; Los, W.; et al. Building essential biodiversity variables (EBVs) of species distribution and abundance at a global scale. Biol. Rev. 2018, 93, 600–625. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  3. McDonald-Madden, E.; Baxter, P.W.J.; Fuller, R.A.; Martin, T.G.; Game, E.T.; Montambault, J.; Possingham, H.P. Monitoring does not always count. Trends Ecol. Evol. 2010, 25, 547–550. [Google Scholar] [CrossRef] [PubMed]
  4. Honrado, J.P.; Pereira, H.M.; Guisan, A. Fostering integration between biodiversity monitoring and modelling. J. Appl. Ecol. 2016, 53, 1299–1304. [Google Scholar] [CrossRef]
  5. Sillero, N. What does ecological modelling model? A proposed classification of ecological niche models based on their underlying methods. Ecol. Modell. 2011, 222, 1343–1346. [Google Scholar] [CrossRef]
  6. Guisan, A.; Thuiller, W. Predicting species distribution: Offering more than simple habitat models. Ecol. Lett. 2005, 8, 993–1009. [Google Scholar] [CrossRef]
  7. Kulhanek, S.A.; Leung, B.; Ricciardi, A. Using ecological niche models to predict the abundance and impact of invasive species: Application to the common carp. Ecol. Appl. 2011, 21, 203–213. [Google Scholar] [CrossRef] [PubMed]
  8. Bradley, B.A. Predicting abundance with presence-only models. Landsc. Ecol. 2016, 31, 19–30. [Google Scholar] [CrossRef]
  9. Weber, M.M.; Stevens, R.D.; Diniz-Filho, J.A.F.; Grelle, C.E.V. Is there a correlation between abundance and environmental suitability derived from ecological niche modelling? A meta-analysis. Ecography 2017, 40, 817–828. [Google Scholar] [CrossRef]
  10. Dallas, T.A.; Hastings, A. Habitat suitability estimated by niche models is largely unrelated to species abundance. Glob. Ecol. Biogeogr. 2018, 27, 1448–1456. [Google Scholar] [CrossRef]
  11. Petrou, Z.I.; Manakos, I.; Stathaki, T. Remote sensing for biodiversity monitoring: A review of methods for biodiversity indicator extraction and assessment of progress towards international targets. Biodivers. Conserv. 2015, 24, 2333–2363. [Google Scholar] [CrossRef]
  12. Mairota, P.; Cafarelli, B.; Didham, R.K.; Lovergine, F.P.; Lucas, R.M.; Nagendra, H.; Rocchini, D.; Tarantino, C. Challenges and opportunities in harnessing satellite remote-sensing for biodiversity monitoring. Ecol. Inform. 2015, 30, 207–214. [Google Scholar] [CrossRef]
  13. Arenas-Castro, S.; Regos, A.; Gonçalves, J.F.; Alcaraz-Segura, D.; Honrado, J.P. Remotely Sensed Variables of Ecosystem Functioning Support Robust Predictions of Abundance Patterns for Rare Species. Remote Sens. 2019, 11, 2086. [Google Scholar] [CrossRef] [Green Version]
  14. Regos, A.; Tapia, L.; Gil-Carrera, A.; Domínguez, J. Monitoring protected areas from space: A multi-temporal assessment using raptors as biodiversity surrogates. PLoS ONE 2017, 12, e0181769. [Google Scholar] [CrossRef] [Green Version]
  15. Regos, A.; Vidal, M.; Lorenzo, M.; Domínguez, J. Integrating intraseasonal grassland dynamics in cross-scale distribution modeling to support waterbird recovery plans. Conserv. Biol. 2020, 34, 494–504. [Google Scholar] [CrossRef]
  16. Tapia, L.; Dominguez, J.; Rodriguez, L. Modelling habitat selection and distribution of Hen harrier (Circus cyaneus) and Montagu’s harrier (Circus pygargus) in a mountainous area in Galicia. J. Raptor Res. 2004, 38, 133–140. [Google Scholar]
  17. Tapia, L.; Domínguez, J.; Rodríguez, L. Hunting habitat preferences of raptors in a mountainous area (northwestern Spain). Polish J. Ecol. 2008, 56, 323–333. [Google Scholar]
  18. Arroyo, B.; Molina, B.; Del Moral, J.C. El Aguilucho Cenizo y el Aguilucho Pálido en España. Población Reproductora en 2017 y Método de Censo; SEO/BirdLife: Madrid, Spain, 2019. [Google Scholar]
  19. Donázar, J.A.; Cortés-Avizanda, A.; Fargallo, J.A.; Margalida, A.; Moleón, M.; Morales-Reyes, Z.; Moreno-Opo, R.; Pérez-García, J.M.; Sánchez-Zapata, J.A.; Zuberogoitia, I.; et al. Roles of Raptors in a Changing World: From Flagships to Providers of Key Ecosystem Services. Ardeola 2016, 63, 181–234. [Google Scholar] [CrossRef] [Green Version]
  20. Tapia, L.; Regos, A.; Gil-Carrera, A.; Dominguez, J. Assessing the temporal transferability of raptor distribution models: Implications for conservation. Bird Conserv. Int. 2018, 28, 375–389. [Google Scholar] [CrossRef]
  21. Tapia, L.; Regos, A.; Gil-Carrera, A.; Domínguez, J. Unravelling the response of diurnal raptors to land use change in a highly dynamic landscape in northwestern Spain: An approach based on satellite earth observation data. Eur. J. Wildl. Res. 2017, 63, 40. [Google Scholar] [CrossRef]
  22. Alvarez-Martínez, J.M.; Suárez-Seoane, S.; Palacín, C.; Sanz, J.; Alonso, J.C. Can eltonian processes explain species distributions at large scale? A case study with great bustard (otis tarda). Divers. Distrib. 2015, 21, 123–138. [Google Scholar] [CrossRef]
  23. Grinnell, J. The Niche-Relationships of the California Thrasher. Auk 1917, 34, 427–433. [Google Scholar] [CrossRef]
  24. Peterson, A.T. Uses and Requirements of Ecological Niche Models and Related Distributional Models. Biodivers. Inform. 2006, 3. [Google Scholar] [CrossRef] [Green Version]
  25. Rodríguez-Lado, L.; Tapia, L.; Pérez, M.; Taboada, T.; Martínez-Capel, F.; Macías, F. Atlas Digital de Propiedades de Suelos de Galicia; Universidad Santiago de Compostela: Santiago de Compostela, Spain, 2016. [Google Scholar]
  26. Pulgar, I. Guía de la flora del Parque Natural Baixa-Limia Serra do Xurés; Consellería de Medio Ambiente; Xunta de Galicia: Santiago de Compostela, Spain, 2005. [Google Scholar]
  27. Rego, P.R.; Guitián, M.A.R.; Castro, H.L.; Da Costa, J.F.; Sobrino, C.M. Loss of european dry heaths in NW Spain: A case study. Diversity 2013, 5, 557–580. [Google Scholar] [CrossRef]
  28. Fagúndez, J. Heathlands confronting global change: Drivers of biodiversity loss from past to future scenarios. Ann. Bot. 2013, 111, 151–172. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  29. Regos, A.; Ninyerola, M.; Moré, G.; Pons, X. Linking land cover dynamics with driving forces in mountain landscape of the Northwestern Iberian Peninsula. Int. J. Appl. Earth Obs. Geoinf. 2015, 38, 1–14. [Google Scholar] [CrossRef]
  30. Andersen, D.E. Survey techniques. In Raptor Research and Management Techniques Manual; Bird, D.M., Bildstein, K.L., Eds.; Hancock House Publishers Ltd: Surrey, BC, Canada, 2007; pp. 89–100. [Google Scholar]
  31. Pettorelli, N.; Vik, J.O.; Mysterud, A.; Gaillard, J.; Tucker, C.J.; Stenseth, N.C.; Lyon, C.B. Using the satellite-derived NDVI to assess ecological responses to environmental change. Trends Ecol. Evol. 2005, 20, 503–510. [Google Scholar] [CrossRef]
  32. Gao, B.C. NDWI—A normalized difference water index for remote sensing of vegetation liquid water from space. Remote Sens. Environ. 1996, 58, 257–266. [Google Scholar] [CrossRef]
  33. Bauer, E.; Kohavi, R.; Chan, P.; Stolfo, S.; Wolpert, D. An Empirical Comparison of Voting Classification Algorithms: Bagging, Boosting, and Variants. Mach. Learn. 1999, 36, 105–139. [Google Scholar] [CrossRef]
  34. Moulds, S. lulcc: Land Use Change Modelling in R. R Package Version 1.0.2. Available online: https://cran.r-project.org/package=lulcc2017 (accessed on 15 October 2020).
  35. Araújo, M.B.; New, M. Ensemble forecasting of species distributions. Trends Ecol. Evol. 2007, 22, 42–47. [Google Scholar] [CrossRef] [PubMed]
  36. Thuiller, W.; Lafourcade, B.; Engler, R.; Araújo, M.B. BIOMOD—A platform for ensemble forecasting of species distributions. Ecography 2009, 32, 369–373. [Google Scholar] [CrossRef]
  37. Marmion, M.; Parviainen, M.; Luoto, M.; Heikkinen, R.K.; Thuiller, W. Evaluation of consensus methods in predictive species distribution modelling. Divers. Distrib. 2009, 15, 59–69. [Google Scholar] [CrossRef]
  38. Bombi, P.; D’Amen, M. Scaling down distribution maps from atlas data: A test of different approaches with virtual species. J. Biogeogr. 2012, 39, 640–651. [Google Scholar] [CrossRef]
  39. Fielding, A.H.; Bell, J.F. A review of methods for the assessment of prediction errors in conservation presence/absence models. Environ. Conserv. 1997, 24, 38–49. [Google Scholar] [CrossRef]
  40. Allouche, O.; Tsoar, A.; Kadmon, R. Assessing the accuracy of species distribution models: Prevalence, kappa and the true skill statistic (TSS). J. Appl. Ecol. 2006, 43, 1223–1232. [Google Scholar] [CrossRef]
  41. Van Horne, B. Density as a Misleading Indicator of Habitat Quality. J. Wildl. Manag. 1983, 47, 893–901. [Google Scholar] [CrossRef]
  42. Hobbs, N.T.; Hanley, T.A. Habitat evaluation: Do use/availability data reflect carrying capacity? J. Wildl. Manag. 1990, 54, 515–522. [Google Scholar] [CrossRef]
  43. Merow, C.; Eckhart, V.M.; Schiffers, K.H.; Thuiller, W.; Georges, D.; Edwards, T.C.; Dullinger, S.; Zurell, D.; Kunstler, G.; Zimmermann, N.E.; et al. Does probability of occurrence relate to population dynamics? Ecography 2014, 37, 1155–1166. [Google Scholar]
  44. Braz, A.G.; de Viveiros Grelle, C.E.; de Souza Lima Figueiredo, M.; Weber, M.d.M. Interspecific competition constrains local abundance in highly suitable areas. Ecography 2020, 1–11. [Google Scholar]
  45. Fernández-Bellon, D.; Lusby, J.; Bos, J.; Schaub, T.; McCarthy, A.; Caravaggi, A.; Irwin, S.; O’Halloran, J. Expert knowledge assessment of threats and conservation strategies for breeding Hen Harrier and Short-eared Owl across Europe. Bird Conserv. Int. 2021, 31, 268–285. [Google Scholar] [CrossRef]
  46. Murgatroyd, M.; Redpath, S.M.; Murphy, S.G.; Douglas, D.J.T.; Saunders, R.; Amar, A. Patterns of satellite tagged hen harrier disappearances suggest widespread illegal killing on British grouse moors. Nat. Commun. 2019, 10, 1–8. [Google Scholar] [CrossRef] [PubMed]
  47. Espín, S.; Terraube, J.; Arroyo, B.; Camarero, P.R.; Mateo, R.; Limiñana, R.; Vázquez-Pumariño, X.; Pinilla, A.; García, J.T.; Mougeot, F. Blood concentrations of p,p′-DDE and PCBs in harriers breeding in Spain and Kazakhstan. Sci. Total Environ. 2018, 624, 1287–1297. [Google Scholar] [CrossRef] [Green Version]
  48. Martínez-Meyer, E.; Díaz-Porras, D.; Townsend Peterson, A.; Yáñez-Arenas, C. Ecological niche structure and rangewide abundance patterns of species. Biol. Lett. 2013, 9, 20120637. [Google Scholar] [CrossRef]
  49. Osorio-Olvera, L.; Soberón, J.; Falconi, M. On population abundance and niche structure. Ecography 2019, 42, 1415–1425. [Google Scholar] [CrossRef]
  50. Osorio-Olvera, L.; Yañez-Arenas, C.; Martínez-Meyer, E.; Peterson, A.T. Relationships between population densities and niche-centroid distances in North American birds. Ecol. Lett. 2020, 23, 555–564. [Google Scholar] [CrossRef] [PubMed]
  51. Van Proosdij, A.S.J.J.; Sosef, M.S.M.M.; Wieringa, J.J.; Raes, N. Minimum required number of specimen records to develop accurate species distribution models. Ecography 2016, 6, 542–552. [Google Scholar] [CrossRef]
  52. Hernandez, P.A.; Graham, C.; Master, L.L.; Albert, D.L. The effect of sample size and species characteristics on performance of different species distribution modeling methods. Ecography 2006, 29, 773–785. [Google Scholar] [CrossRef]
  53. Bean, W.T.; Stafford, R.; Brashares, J.S. The effects of small sample size and sample bias on threshold selection and accuracy assessment of species distribution models. Ecography 2012, 35, 250–258. [Google Scholar] [CrossRef]
  54. Arroyo, B.; García, J.; Betagnolle, V. Montagu’s Harrier Circus pygargus. BWP Update 2004, 6, 41–55. [Google Scholar]
  55. Arroyo, B.; Leckie, F.; Amar, A.; McCluskie, A.; Redpath, S. Ranging behaviour of Hen Harriers breeding in Special Protection Areas in Scotland. Bird Study 2014, 61, 48–55. [Google Scholar] [CrossRef] [Green Version]
  56. Regos, A.; Arenas-Castro, S.; Tapia, L.; Domínguez, J.; Honrado, J. Using remotely sensed indicators of primary productivity to improve prioritization of conservation areas for top predators. Ecol. Indic. 2021, 125, 107503. [Google Scholar] [CrossRef]
  57. Requena-Mullor, J.M.; López, E.; Castro, A.J.; Cabello, J.; Virgós, E.; González-Miras, E.; Castro, H. Modeling spatial distribution of European badger in arid landscapes: An ecosystem functioning approach. Landsc. Ecol. 2014, 29, 843–855. [Google Scholar] [CrossRef]
  58. Regos, A.; Gagne, L.; Alcaraz-Segura, D.; Honrado, J.P.; Domínguez, J. Effects of species traits and environmental predictors on performance and transferability of ecological niche models. Sci. Rep. 2019, 9, 4221. [Google Scholar] [CrossRef]
  59. Alcaraz-Segura, D.; Lomba, A.; Sousa-Silva, R.; Nieto-Lugilde, D.; Alves, P.; Georges, D.; Vicente, J.R.; Honrado, J.P. Potential of satellite-derived ecosystem functional attributes to anticipate species range shifts. Int. J. Appl. Earth Obs. Geoinf. 2017, 57, 86–92. [Google Scholar] [CrossRef]
  60. Caravaggi, A.; Irwin, S.; Lusby, J.; Mccarthy, A.; Mee, A.; Nagle, T.; O’Halloran, J. Forest management and Hen Harrier Circus cyaneus conservation in Ireland. Irish Birds 2020, 42, 1–12. [Google Scholar]
  61. Monaghan, F.; Sullivan, C. Hen Harrier Programme Field Guidance for Scoring Bog and Heath; Hen Harrier Project: Galway, Ireland, 2019. [Google Scholar]
  62. Thirgood, S.J.; Redpath, S.M.; Graham, I.M. What determines the foraging distribution of raptors on heather moorland? Oikos 2003, 100, 15–24. [Google Scholar] [CrossRef]
  63. Pe’er, G.; Bonn, A.; Bruelheide, H.; Dieker, P.; Eisenhauer, N.; Feindt, P.H.; Hagedorn, G.; Hansjürgens, B.; Herzon, I.; Lomba, Â.; et al. Action needed for the EU Common Agricultural Policy to address sustainability challenges. People Nat. 2020, 2, 305–316. [Google Scholar] [CrossRef]
  64. Lomba, A.; Moreira, F.; Klimek, S.; Jongman, R.H.G.; Sullivan, C.; Moran, J.; Poux, X.; Honrado, J.P.; Pinto-Correia, T.; Plieninger, T.; et al. Back to the future: Rethinking socioecological systems underlying high nature value farmlands. Front. Ecol. Environ. 2019, 18, 36–42. [Google Scholar] [CrossRef]
Animals 11 02020 g001 550
Figure 1. Location of the study area in the “Serra do Laboreiro”. False-color RGB composite from Landsat 8 OLI sensor imagery (white and black lines represent 1 km and 2 km buffer areas around each point count, respectively; white and black circles represent each transect).
Figure 1. Location of the study area in the “Serra do Laboreiro”. False-color RGB composite from Landsat 8 OLI sensor imagery (white and black lines represent 1 km and 2 km buffer areas around each point count, respectively; white and black circles represent each transect).
Animals 11 02020 g001
Animals 11 02020 g002 550
Figure 2. Ranked importance of the environmental variables included in ecological niche models of hen harrier and Montagu’s harrier. Values represent the importance of each variable accumulated from the different modeling techniques. GLM, generalized linear models; CTA, classification tree analysis; ANN, artificial neural networks; SRE, surface range envelope; FDA, flexible discriminant analysis; MARS, multivariate adaptive regression splines; RF, random forest; MaxEnt, maximum entropy; Slo, slope; Far, farmland; Alt, altitude; Asp, aspect; Cl_H, closed heathland; W_ind, normalized difference water index; Op_H, open heathland; P_pl, pine plantations.
Figure 2. Ranked importance of the environmental variables included in ecological niche models of hen harrier and Montagu’s harrier. Values represent the importance of each variable accumulated from the different modeling techniques. GLM, generalized linear models; CTA, classification tree analysis; ANN, artificial neural networks; SRE, surface range envelope; FDA, flexible discriminant analysis; MARS, multivariate adaptive regression splines; RF, random forest; MaxEnt, maximum entropy; Slo, slope; Far, farmland; Alt, altitude; Asp, aspect; Cl_H, closed heathland; W_ind, normalized difference water index; Op_H, open heathland; P_pl, pine plantations.
Animals 11 02020 g002
Animals 11 02020 g003 550
Figure 3. Maps depicting areas predicted by the ecological niche models as suitable in 1997 and 2017 (“Overlap”), only suitable in 1997 (“Losses”) and only suitable in 2017 (“Gains”).
Figure 3. Maps depicting areas predicted by the ecological niche models as suitable in 1997 and 2017 (“Overlap”), only suitable in 1997 (“Losses”) and only suitable in 2017 (“Gains”).
Animals 11 02020 g003
Animals 11 02020 g004 550
Figure 4. Circular plot illustrating the changes in land cover between 1997 and 2017 within areas predicted by the ecological niche models as suitable habitat in 1997 and 2017 (“Overlap”), only suitable in 1997 (“Losses”) and only suitable in 2017 (“Gains”). The width of the lines is proportional to the contribution of each type of land type to the change. The different colors represent the different types of land cover.
Figure 4. Circular plot illustrating the changes in land cover between 1997 and 2017 within areas predicted by the ecological niche models as suitable habitat in 1997 and 2017 (“Overlap”), only suitable in 1997 (“Losses”) and only suitable in 2017 (“Gains”). The width of the lines is proportional to the contribution of each type of land type to the change. The different colors represent the different types of land cover.
Animals 11 02020 g004
Animals 11 02020 g005 550
Figure 5. The species abundance–habitat suitability relationship: (A) degree of correlation between observed species abundance and habitat suitability predicted from ENMs for each species and year (Montagu’s harrier (Spearman’s ρ = 0.57; p-value < 0.001) and hen harrier (Spearman’s ρ = 0.55; p-value < 0.001)); (B) degree of correlation between species abundance variation and habitat suitability changes predicted from ENMs for each species between 1997 and 2017 (hen harrier (Pearson’s r = 0.34; p-value = 0.12) and Montagu’s harrier (Pearson’s r = 0.31; p-value = 0.17)).
Figure 5. The species abundance–habitat suitability relationship: (A) degree of correlation between observed species abundance and habitat suitability predicted from ENMs for each species and year (Montagu’s harrier (Spearman’s ρ = 0.57; p-value < 0.001) and hen harrier (Spearman’s ρ = 0.55; p-value < 0.001)); (B) degree of correlation between species abundance variation and habitat suitability changes predicted from ENMs for each species between 1997 and 2017 (hen harrier (Pearson’s r = 0.34; p-value = 0.12) and Montagu’s harrier (Pearson’s r = 0.31; p-value = 0.17)).
Animals 11 02020 g005
Table
Table 1. Landscape change (ha) within 500 m buffer around each point count.
Table 1. Landscape change (ha) within 500 m buffer around each point count.
19972017
Farmland47.7983.07
Closed heathland728.19388.08
Open heathland402.57695.97
Pine plantations152.82164.25
Table
Table 2. Habitat availability (ha) predicted by ecological niche models for 1997 and 2017.
Table 2. Habitat availability (ha) predicted by ecological niche models for 1997 and 2017.
19972017
Hen harrier (Circus cyaneus)202.77121.32
Montagu’s harrier (Circus pygargus)140.22258.48
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Regos, A.; Tapia, L.; Gil-Carrera, A.; Domínguez, J. Caution Is Needed When Using Niche Models to Infer Changes in Species Abundance: The Case of Two Sympatric Raptor Populations. Animals 2021, 11, 2020. https://doi.org/10.3390/ani11072020

AMA Style

Regos A, Tapia L, Gil-Carrera A, Domínguez J. Caution Is Needed When Using Niche Models to Infer Changes in Species Abundance: The Case of Two Sympatric Raptor Populations. Animals. 2021; 11(7):2020. https://doi.org/10.3390/ani11072020

Chicago/Turabian Style

Regos, Adrián, Luis Tapia, Alberto Gil-Carrera, and Jesús Domínguez. 2021. "Caution Is Needed When Using Niche Models to Infer Changes in Species Abundance: The Case of Two Sympatric Raptor Populations" Animals 11, no. 7: 2020. https://doi.org/10.3390/ani11072020

APA Style

Regos, A., Tapia, L., Gil-Carrera, A., & Domínguez, J. (2021). Caution Is Needed When Using Niche Models to Infer Changes in Species Abundance: The Case of Two Sympatric Raptor Populations. Animals, 11(7), 2020. https://doi.org/10.3390/ani11072020

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop