Introduction
Human settlements and their concomitant transport infrastructure alter the landscape, leading to loss and fragmentation of wildlife habitats (Trombulak & Frissell, Reference Trombulak and Frissell2000; Turner et al., Reference Turner, Nakamura and Dinetti2004) and thus posing significant threats to local biodiversity (Benítez-López et al., Reference Benítez-López, Alkemade and Verweij2010). Roads and other infrastructure, for example, result in direct mortality, altered behaviour, altered physical and chemical habitat, dissemination of alien species, increased anthropogenic disturbance and use of habitats, and fragmentation of populations (Trombulak & Frissell, Reference Trombulak and Frissell2000; Brugiere et al., Reference Brugiere, Badjinca, Silca, Serra and Barry2006; Benítez-López et al., Reference Benítez-López, Alkemade and Verweij2010; Sinclair, Reference Sinclair2010; González-Gallina et al., Reference González-Gallina, Benítez-Badillo, Rojas-Soto and Hidalgo-Mihart2013). Urbanization also alters animal communities in a complex way related to the size, density and location of human populations (van der Zande et al., Reference van der Zande, ter Keurs and van der Weijden1980; Clergeau et al., Reference Clergeau, Croci, Jokimäki, Kaisanlahti-Jokimäki and Dinetti2006; Vermaat et al., Reference Vermaat, Goosen and Omtzigt2007; Pedrana et al., Reference Pedrana, Bustamante, Travaini and Rodríguez2010; Randrianandrianina et al., Reference Randrianandrianina, Racey and Jenins2010). Identifying which spatial pattern of human settlement has the least impact on biodiversity is important for the conservation of global biological diversity (Sutherland et al., Reference Sutherland, Adams, Aronson, Aveling, Blackburn and Broad2009).
Birds are one of the most studied groups with respect to the sensitivity of wildlife to habitat disturbance caused by infrastructure, and there is significant variability within the group (Benítez-López et al., Reference Benítez-López, Alkemade and Verweij2010). The effect of transport infrastructure on raptors, for example, is one of either indifference or attraction for small and medium-sized species (Meunier et al., Reference Meunier, Verheyden and Jouventin2000; Lambertucci et al., Reference Lambertucci, Speziale, Rogers and Morales2009) but of repulsion for larger species (Bautista et al., Reference Bautista, García, Calmaestra, Palacín, Martín and Morales2004; Lambertucci et al., Reference Lambertucci, Speziale, Rogers and Morales2009), which affects their breeding sites (Reijnen & Foppen, Reference Reijnen and Foppen1994) and ultimately their reproductive success (e.g. Anthony & Isaacs, Reference Anthony and Isaacs1989; Peach et al., Reference Peach, Vincent, Fowler and Grice2008).
Reviews of the known effects of infrastructure (Trombulak & Frissell, Reference Trombulak and Frissell2000; Fahrig & Rytwinski, Reference Fahrig and Rytwinski2009; Benítez-López et al., Reference Benítez-López, Alkemade and Verweij2010; Selva et al., Reference Selva, Kreft, Kati, Schluck, Jonsson and Mihok2011) indicate the conservation importance of maintaining existing sparsely developed areas. It is a misconception that most impacts are associated with major roads (Sun et al., Reference Sun, Dong, Zhang, Zheng and Browne2009; van Langevelde et al., Reference van Langevelde, van Dooremalen and Jaarsma2009). However, the impact of relatively minor infrastructure (e.g. minor roads and unpaved tracks) in areas of low development remains largely unknown (Kociolek et al., Reference Kociolek, Clevenger, Clair and Proppe2011). Unpaved tracks and isolated buildings (henceforth, secondary infrastructure) occupy significant areas even in developed countries (e.g. in Europe; EEA, 2006).
We hypothesized that the breeding site selection and success (chicks flying from nests) of large birds is influenced by human settlements and transit infrastructure, particularly elements of secondary infrastructure in sparsely developed and relatively well-conserved landscapes. Given the known sensitivity of large birds to disturbance and the characteristic ranges of their response to such disturbance, the use they make of sparsely developed areas provides a good model with which to investigate the potential impacts of secondary infrastructure. We studied the spatial covariation of infrastructure and the breeding sites of large birds in a sparsely developed landscape, to address the following questions: Is there a repulsive influence of secondary infrastructure on the breeding site selection of large birds? Is it related to human dwellings, transit or both? In the first case, do isolated buildings exert an influence comparable to rural settlements? In the second case, do unpaved tracks exert an influence comparable to paved roads? Is the negative influence independent of, or associated with, features of the landscape?
Study area and species
We conducted the study in the Monfragüe Biosphere Reserve, a sparsely developed area with a typical agroforestry system that dominates the landscape of the south-western Iberian Peninsula, called dehesa in Spain, where it covers c. 5.8 million ha (Joffre et al., Reference Joffre, Rambal and Ratte1999), and montado in Portugal, where it covers c. 1.2 million ha (Mendes, Reference Mendes, Merlo and Croitoru2005). Traditional dehesa land use has been a sustainable system (Dawson & Fry, Reference Dawson and Fry1998; García del Barrio et al., Reference García Del Barrio, Bolaños, Ortega and Elena-Roselló2004), managed for centuries while maintaining high biodiversity and minimal transformation of the landscape (Díaz et al., Reference Díaz, Campos, Pulido, Pain and Pienkowski1997; Blondel & Aronson, Reference Blondel and Aronson1999; Plieninger & Wilbrand, Reference Plieninger and Wilbrand2001). The Reserve is within the Mediterranean Biogeographical Region (Fig. 1). It consists of a 12,830 ha core zone, which is the Monfragüe National Park, with conservation and recreational use, a 15,360 ha buffer zone, and an 87,970 ha transition zone, which overlaps the Monfragüe and Surrounding Dehesas Special Protection Area (116,160 ha; European Birds Directive 2009/147/EC). There are 46 villages and hamlets within the Reserve, the seven largest of which have a total of 3,384 inhabitants (mean = 631.0 ± SD 585.8). The transport infrastructure consists of 2,659.1 km of unpaved tracks (89.3%), 309.8 km of paved roads (10.4%) and 7.8 km of motorway (0.3%). The Reserve includes a combination of natural (forest and scrub) and moderately transformed (dehesa) areas, which support high biodiversity (> 1,300 species of flora and > 200 vertebrates; GIC, 2010), including a significant percentage of European populations of large bird species of conservation interest (e.g. the black stork Ciconia nigra, 0.4–0.5%; Egyptian vulture Neophron percnopterus, 1.6–1.8%; Iberian imperial eagle Aquila adalberti, 6.7%; cinereous vulture Aegypius monachus, 22.2%; BirdLife International, 2004).
The infrastructure elements, hydrological network, and delimitation of the core zone of the Monfragüe Biosphere Reserve in the south-western Iberian Peninsula.
We included in the analysis the locations of the nesting sites of large rock-dwelling and tree-dwelling birds, including the Accipitridae (diurnal raptors) and Ciconiidae (storks) families. Our criteria were to focus on the larger species of greatest conservation interest and with the best understood distributions. Long-term data from annual censuses of the Directorate General for the Environment (DGMN, Junta de Extremadura) and studies of the Conservation Biology Research Group (GIC, University of Extremadura) are available in databases for the period 1986–2002. These data are the product of coordinated, synchronized censuses covering the entire study area, carried out during the breeding season by rangers and researchers of these institutions (e.g. Morán-López et al., Reference Morán-López, Sánchez Guzmán, Borrego and Sánchez2006a,Reference Morán-López, Sánchez, Costillo, Corbacho and Villegasb). More recent data that fulfilled our criteria were unavailable. The 1986–2002 data are the best available, up-to-date, long-term data for a sparsely developed Mediterranean landscape of conservation relevance, and we verified their attributes and spatial quality. Reports (e.g. GIC, 2010) indicate that the environmental and socio-economic situation in the Reserve has remained essentially constant to the present.
Methods
To facilitate a reliable comparison of census data and consideration of all the nesting areas we established a spatio-temporal scale defined by a 1 km grid (UTM grid) and multi-year censuses (1986–2002). This temporal grain size is greater than that of any of the available surveys (e.g. Northrup et al., Reference Northrup, Stenhouse and Boyce2012) and facilitates a more reliable description of large-bird breeding conditions (Gutzwiller & Barrow, Reference Gutzwiller and Barrow2003). Infrastructure effects are probably synergistic, complex, cumulative and time-lagged (EEA, 2011; Selva et al., Reference Selva, Kreft, Kati, Schluck, Jonsson and Mihok2011), and thus relationships derived from a single census and year would not describe them adequately. With respect to the spatial grain size, numerous studies have found scales of one to several kilometres to be adequate to study bird distribution, habitat use and the effects of infrastructure (e.g. Gutzwiller & Barrow, Reference Gutzwiller and Barrow2003; Morán-López et al., Reference Morán-López, Sánchez Guzmán, Borrego and Sánchez2006a; Benítez-López et al., Reference Benítez-López, Alkemade and Verweij2010; Schwenk & Donovan, Reference Schwenk and Donovan2011, and references therein). The 1 km filtering rule has been found to be useful for eliminating spatial autocorrelation resulting from clumped observations (Telford & Birks, Reference Telford and Birks2009). Accordingly, we used this same spatial scale to characterize both species distribution and predictors.
We followed a multispecies study framework, as such an approach has been found to be useful in quantitative assessments for conservation planning (Bonn & Schröder, Reference Bonn and Schröder2001; Dallimer et al., Reference Dallimer, King and Atkinson2009; Schwenk & Donovan, Reference Schwenk and Donovan2011). Multispecies models (1) constitute a more parsimonious, general approach, which follows the precautionary principle in identifying key human influences; (2) overcome difficult and subjective trade-offs among species; and (3) facilitate understanding of collective responses of communities to human influences. For all the grid squares in the study area we determined the presence or absence of nesting sites of any of the study species in any year.
The analysis of dichotomous variables may be negatively affected by highly imbalanced prevalences of the species in the data (i.e. the presence/absence ratio having many more absences than presences; Franklin et al., Reference Franklin, Wejnert, Hathaway, Rochester and Fisher2009). Our sampling strategy was to take random samples of absence data of the same size as the presence datasets. To avoid results being contingent on just one subsample, we followed a bootstrap resampling approach while maintaining prevalence and sample size for robust parameter and confidence interval estimation (Burnham & Anderson, Reference Burnham and Anderson2002).
Landscape and infrastructure data
The grid squares were characterized by digital cartography in a geographical information system used to measure the explanatory variables. We selected potential mediating factors of human disturbance (e.g. noise, light) and accessibility (cost of transiting the terrain) in the environment of nests. In each grid square we measured 16 variables (Table 1) characterizing the spatial distribution of the transit infrastructure (paved or unpaved) and human dwellings (population in clustered settlements or isolated buildings), and natural elements of the landscape that impeded accessibility (relief and hydrography), increasing the cost of transiting the terrain. We evaluated the functioning of alternative spatial measurements (as we had no prior knowledge of which type of measure would be a better indicator of spatially mediated effects), using densities, distances and percentage occupancy of the grid squares. The measurements were made in the 1 × 1 km UTM grid and in 3 × 3 and 5 × 5 grids (Table 1). We were thus able to evaluate uncertainty related to the spatial range of influence of the potential factors mediating accessibility. However, some variables (the form of the terrain and the minimum distances) were measured only in the 1 × 1 km grid. This was under the assumption that their potential influence would be exerted in the vicinity of the nests and not over an extended spatial environment (e.g. Morán-López et al., Reference Morán-López, Sánchez Guzmán, Borrego and Sánchez2006a,Reference Morán-López, Sánchez, Costillo, Corbacho and Villegasb).
Variables used to characterize the UTM grid applied to the Monfragüe Biosphere Reserve in the south-west Iberian Peninsula (Fig. 1), with type, unit, and description.
Data analysis
We analysed all data using SPSS 15.0 (IBM, Armonk, USA). We performed a non-parametric Spearman correlation analysis (Zar, Reference Zar1996) between all the predictors, to exclude those that were highly correlated (r s > 0.9, P < 0.05) and to facilitate interpretation. We also inspected bivariate scatterplots relating the dependent variable to its predictors, to establish the expected form (linear, monotonic, etc.) of the responses and decide on the most suitable transformations to include in the models. The variables were subsequently standardized to avoid bias resulting from differences in the measurement units.
We investigated the hypothesized impacts by constructing a set of generalized linear models using a binomial distribution and a logistic link function appropriate for dichotomous data. To evaluate the individual and collective explanatory capacity of each type of factor (distance, density or relief) and scale of influence (1, 3 or 5 km) on nest site presence we constructed models for each type of descriptor and scale (individual models), then all types by scale (conjoint models), and finally all types and scales (full model). The total set comprised 11 models and 25 parameters. We used a model selection approach based on the bias-adjusted Akaike's information criterion (AICc; Burnham & Anderson, Reference Burnham and Anderson2002) to provide evidence for the hypothesized impacts (and their scale), to quantify the independent influence of infrastructure vs relief variables, to compare the contribution of secondary vs primary infrastructure, and to identify combinations of variables that would be useful to guide species conservation (see, for example, inferential strength and competing hypotheses in Stephens et al., Reference Stephens, Buskirk, Hayward and Martínez Del Río2005; Roedenbeck et al., Reference Roedenbeck, Fahrig, Findlay, Houlahan, Jaeger and Klar2007).
Calculating parameter estimates and confidence intervals conditional on only one random sample would not explore the uncertainty in the identification of factors that could have utility and general application to species conservation. We therefore used bootstrap resampling of the data, maintaining sample size and prevalence (B = 999 samples of size R = 312), followed by model selection. We used model selection relative frequencies (p i ) to select the best model, and took as included in the best model the estimated parameters whose bootstrapped confidence limits excluded the value zero (Burnham & Anderson, Reference Burnham and Anderson2002).
We found that the location of the infrastructure was not independent of the characteristics of the relief (Spearman analysis), which could result in inferences that confuse the independent effects of the former (e.g. disturbances) with other effects associated with the latter (e.g. land use or habitat) (Fahrig & Rytwinski, Reference Fahrig and Rytwinski2009). We therefore removed the possibility of the influence of habitat in the models by including the infrastructure variables after removing the effect of relief, a proxy for habitat, to which end we replaced the predictors by the residual of a regression of the (dependent) infrastructure variables against the (independent) relief (Legendre, Reference Legendre1993). This has the effect of removing the influence of land use on infrastructure variables through relief (Fahrig & Rytwinski, Reference Fahrig and Rytwinski2009).
Socio-economic, environmental and biological data may reflect spatial patterns that depend on the factors generating them (Anselin, Reference Anselin1988; Carroll & Pearson, Reference Carroll and Pearson2000). It is advisable to measure any potential spatial autocorrelation in the data or the residuals, as its presence may be an indicator of absent factors or contagious processes, as well as of pseudo-replication, and to include or exclude it from the models if it is present. We created Moran's I correlograms (Legendre & Legendre, Reference Legendre and Legendre1998) with the standardized residuals of the models (lag 1 km, extent 50 km) to determine the existence of spatial autocorrelation in the error. A Bonferroni correction was applied because the technique involves multiple tests.
Confidence in the performance of the models generated and their utility for species conservation is founded in proper model validation. For this purpose we followed a three-way approach, combining the attributes of same/new data, same/new type of data and same/new location of data. Firstly, we calculated receiver operating characteristic plots for 999 random bootstrap samples of size 40 (Buckland & Elston, Reference Buckland and Elston1993) in the Monfragüe Biosphere Reserve, and measured accuracy as a mean and standard deviation of the area under the curve (same data, same type, same place). Secondly, we calculated the receiver operating characteristic plot and the area under the curve value in the Sierra de San Pedro Special Protection Area, located c. 90 km from the Reserve (new data, same type, different place). Thirdly, we correlated the summed breeding success by 1 km2 of the most abundant raptor, the cinereous vulture (data from 2000; Morán-López et al., Reference Morán-López, Sánchez, Costillo, Corbacho and Villegas2006b), with the linear predictor of the model (new data, different type, same place), using the Spearman coefficient.
Finally, we developed distance effect functions to apply the new knowledge generated by the study to animal conservation practice. To this end, we created empirical probability of occupancy models relating the large-bird community to the secondary infrastructure of unpaved tracks and isolated buildings. Individual models by variable (including its quadratic term) were generated using the bootstrap generalized linear modelling approach described above (B = 999, R = 312) with the original data (i.e. no transformation or standardization). A bivariate track vs building model was also created to further explore interaction effects in addition to the aforementioned nonlinear relations.
Results
Species and nesting sites
Seven species, including threatened species, met the criteria of body size and availability of monitoring data of sufficient quality. These selected species nest in 156 (11.7%) of the 1,336 1-km UTM squares of the Monfragüe Biosphere Reserve, thus leaving 1,180 (88.3%) unoccupied squares (Table 2). Presence squares were denser in the core zone and the cinereous vulture was the most frequent species, but the contribution of the group of remaining species was comparable and the resulting locations were distributed across all three zones of the Reserve. This provides support for the limits of the Reserve and the multi-species and multi-year data being adequate to describe breeding environments according to the study objectives.
Number of 1 km grid squares in the Monfragüe Biosphere Reserve (Fig. 1) where large-bird species were present; % of the total number of squares occupied by any of the species (n = 156); and % of the total number of squares in the Reserve, with or without species presence (n = 1,336).
* The percentages are non-additive, as squares may be used by more than one species.
Influence of infrastructure on breeding sites
Infrastructure and relief at the 1 km scale (conjoint 1 km model) received strong support in the bootstrap model selection procedure, whereas the remaining models received weak or no support (Table 3). No predictors other than distance to infrastructure and the relief characteristics within 1 km received substantial support at any scale. Both primary and secondary infrastructure, including transport and settlement, were included in the best model (conjoint 1 km model; Table 4). The influence of secondary infrastructure was comparable to that of primary infrastructure, both in transport (tracks vs roads) and settlement (isolated buildings vs villages), and relief showed a strong effect (Table 4).
Results of the bootstrap model selection procedure, with mean values and standard deviations of AICc, and selection frequencies.
* Each model is based on different predictors (distance, density or relief) and scales (1, 3 or 5 km), combined by scale (conjoint models) or including all types and scales (full model).
Results of the conjoint 1 km model bootstrap resampling estimation of parameters.
* Only the predictors whose 99% bootstrapped confidence limits excluded zero are shown.
The validation results indicated that the conjoint 1 km model performed well both inside and outside the study area. The receiver operating characteristic plots showed high accuracy inside the Reserve (area under the curve = 0.9), with the model also being useful for application in the Special Protection Area (area under the curve = 0.7; Table 5; Fig. 2). Probability thresholds equilibrating sensitivity (true positive fraction) and specificity (true negative fraction) showed correct classification of > 80% and > 70% (including breeding and non-breeding areas) inside (Reserve) and outside (Special Protection Area) the study area, respectively. Breeding success of the cinereous vulture in the Reserve was positively correlated with the linear predictor of the model (r s = 0.28, P = 0.04, n = 52).
Receiver operating characteristic plot with mean and standard deviation of the true positive fraction, for the conjoint 1 km model (see text for details) validated in the Monfragüe Biosphere Reserve (Fig. 1) using bootstrap resampling; AUC = 0.899 ± SD 0.051.
Results of the conjoint 1 km model validation process inside (Monfragüe Biosphere Reserve) and outside (Sierra de San Pedro Special Protection Area) the calibrating area; the results for the Reserve are from a bootstrap validation approach. Threshold-dependent results are shown at cut-off values that equilibrate positive and negative misclassifications.
We found the covariance of relief and infrastructure to be associated with the slope in all the infrastructure variables, with increasing distance and decreasing density of infrastructure elements corresponding to greater slopes (corrector model results not shown), confirming the adequacy of our approach in removing the effect of relief on infrastructure variables. No model showed spatial autocorrelation in the residuals (Moran's I, P > 0.05 in all cases). We could therefore exclude the existence of spatial autocorrelation.
Empirical models relating probability of occupancy to secondary infrastructure showed significant linear relationships in every descriptor (P < 0.05), but no quadratic ones (P > 0.10). Distance to tracks, isolated buildings, and their interaction were significant in these models (Figs 3 & 4). The probability of large birds breeding in the Reserve was higher the more distant the site from tracks and buildings, with tracks showing a more intense shorter distance effect than buildings. A probability-of-occurrence range of 0.95–0.99 is reached at distances of 0.5–0.7 km from tracks and 3.4–4.5 km from buildings. The areas of the Reserve contained within these probabilities are 95.2–99.0% for tracks and 100% for isolated buildings.
Predicted probability of occurrence of breeding of large birds in relation to distance from unpaved tracks (disttrack) or isolated buildings (distbuild). Linear predictors of the logistic models are (coefficients P < 0.05 included; quadratic coefficients P > 0.10 excluded): P(disttack) = −1.18 + 8.37 disttrack; P(distbuild) = −1.32 + 1.24 distbuild.
Predicted probability of occurrence of breeding of large birds in relation to the interaction between distance from unpaved tracks (disttrack) and isolated buildings (distbuild). The linear predictor of the logistic model is (coefficients P < 0.0001): −0.96 + 6.37 disttrack × distbuild.
Discussion
Previous studies have shown that both the presence of the principal transport infrastructure (paved roads) and the intensity of traffic alter the behaviour of large birds in their foraging and breeding areas (Meunier et al., Reference Meunier, Verheyden and Jouventin2000; Trombulak & Frissell, Reference Trombulak and Frissell2000; Bautista et al., Reference Bautista, García, Calmaestra, Palacín, Martín and Morales2004; Kociolek et al., Reference Kociolek, Clevenger, Clair and Proppe2011). We have demonstrated the significant effect of secondary infrastructure even in an area of sparse development, where unpaved tracks and isolated buildings have displaced the breeding sites of large birds. Furthermore, although the intensity of the effect was comparable to that of primary infrastructure, the spatial extension of the effect of secondary infrastructure was far greater (by a factor of c. 9), resulting in a larger size effect.
This comparable negative effect of secondary infrastructure relative to that of primary infrastructure is particularly informative, as it has commonly been assumed that low levels of development in rural areas have relatively little or no impact on wildlife populations. This assumption has resulted in the potential effect of secondary infrastructure receiving only minor attention in the scientific literature and management practice (e.g. Kociolek et al., Reference Kociolek, Clevenger, Clair and Proppe2011). However, mitigation of these negative effects is important in protected areas and in general in relatively undisturbed areas where unpaved tracks constitute a major fraction of a transit infrastructure that extends over large areas, which are supposed to buffer wildlife populations from human activities. In Spain's Extremadura region, for example, which is an important region for biodiversity conservation, unpaved tracks represent 84.3% of the transport infrastructure, 86.0% of the Natura 2000 Network transport infrastructure, and 89.3% of the Monfragüe Biosphere Reserve transport infrastructure. Considering that the principal strategy used to mitigate the negative effects of infrastructure is to minimize infrastructure development in relatively undisturbed areas (Trombulak & Frissell, Reference Trombulak and Frissell2000; Benítez-López et al., Reference Benítez-López, Alkemade and Verweij2010), neglecting secondary infrastructure in wildlife conservation management is a shortcoming.
The breeding habitats of the large birds of the Monfragüe Biosphere Reserve are characterized by gentle slopes, whereas areas of more rugged relief are avoided. Variable relief offers greater protection from predators, and acts as a buffer against human disturbance. Breeding in rugged areas would be limited by factors such as the reduced availability of suitable habitats, and adverse weather conditions (Newton, Reference Newton1979). The predominant explanatory power of relief in the best model marks it as an indicator of suitable breeding habitats for large birds, as well as being a limiter for accessibility.
The results suggest that the construction of unpaved tracks and isolated buildings has a negative impact on large bird species, regardless of the forms of the relief. This points to their breeding sites being sensitive to such infrastructure not only in the more rugged core zone but also in the buffer and transition zones. This knowledge provides a common framework for spatial and year-long planning throughout the zones of the Reserve, not just in the core zone, which is preferentially targeted for conservation and recreational uses. We therefore recommend that the principal strategy to minimize infrastructure in relatively undisturbed areas should be extended to also include secondary infrastructure (Trombulak & Frissell, Reference Trombulak and Frissell2000; Benítez-López et al., Reference Benítez-López, Alkemade and Verweij2010; Selva et al., Reference Selva, Kreft, Kati, Schluck, Jonsson and Mihok2011). Avoidance of infrastructure by birds, both colonial and territorial, may be the cause of anything from spatial displacement to an over-saturation of breeding habitats, resulting in net population losses because of the effect on reproductive success and productivity (van der Zande et al., Reference van der Zande, ter Keurs and van der Weijden1980).
The increased accessibility and use of areas in proximity to nests as a result of the construction of secondary infrastructure (tracks and isolated buildings) has had a generalized negative effect on habitat selection by large birds, and on their breeding success. The levelling off of large-bird populations in the Reserve (GIC, 2010) is indicative of a state of over-saturation that may reduce the chances of recovery of threatened species, a state conditioned, at least in part, by the development of secondary infrastructure. We found a negative effect of infrastructure on the breeding success of the largest raptor (A. monachus) in the core zone of the Reserve, where most of the infrastructure is secondary.
Once this general negative response has been confirmed, an analysis by species, and a specific consideration of non-linear relationships and breeding success, will help to reveal the mechanisms and costs involved (e.g. D'Amico et al., Reference D'Amico, Périquet, Román and Revilla2016). For species-specific models to be useful in species conservation, the differences in model outcomes must be interpreted properly with respect to the biology and ecology of each species (e.g. Mihoub et al., Reference Mihoub, Jiguet, Lécuyer, Eliotout and Sarrazin2014), and the species’ mechanistic responses to single or interacting human factors must be disentangled appropriately. The specific uses of every track and isolated building should also be considered, differentiating transport from access, and agricultural from recreational uses, and where and when they take place. However, this is not an easily affordable task in terms of resources and time. The community-level, multi-year infrastructure response models developed in this study provide a method with which to bridge the gap between conservation objectives and management practice.
Conclusions
Unpaved tracks and isolated buildings are the main centres of human activity in protected spaces, private estates and sparsely populated rural areas with semi-natural or natural land uses. This is why the effect of distance on the breeding sites of large birds indicates a widespread impact of human activity, wider than that of primary infrastructure in these areas. It will be a challenge to discern and separate the causal bases for each of the affected species. Nonetheless, the common pattern revealed by this community of breeding large birds provides evidence of the need for conservation management, primarily the development of a strategic plan that includes a spatial design for the construction and maintenance of secondary infrastructure, and its relocation or elimination where appropriate. Distance effect functions showed an influence of tracks only in close proximity (≤ 700 m) to breeding sites, spatially restricting the necessity for such management. The politically and socially admissible limits can be set based on empirical data and distance effect functions such as those reported here.
Based on current evidence it would be advisable to assess the contribution of the existing network of buildings and tracks, (not only in protected areas but also in sparsely developed regions, beyond them) to determine the appropriate actions to take. Such information would facilitate quantitative, evidence-based decision making regarding the extent of secondary infrastructure removal/displacement that is needed to achieve a given increase in the probability of a threatened species breeding, or to avoid the abandonment of the use of the area for breeding.
Acknowledgements
The data used in the study were provided by the General Directorate of the Environment of the Junta de Extremadura and the Group of Investigation in Conservation Biology of the University of Extremadura. We thank in particular the environment agents responsible for obtaining a better understanding of the distribution of bird populations in the Monfragüe Biosphere Reserve and the Sierra de San Pedro Special Protection Area, and the estate owners who allowed access to their properties. The study was funded through the agreement between the University of Extremadura and the Ministry of the Environment and Rural and Marine Affairs entitled ‘Trabajos de seguimiento del estado ambiental de las Reservas de la Biosfera de la Comunidad Autónoma de Extremadura [Works monitoring the environmental status of the Biosphere Reserves of the Autonomous Community of Extremadura]’. The Consejería de Economía, Comercio e Innovación (Junta de Extremadura) funded the translation of the manuscript (GR10174).
Biographical sketches
Ricardo Morán-López’s research focuses mainly on the study of relationships between species and habitats, species distribution modelling, and the design and evaluation of reserves. Juan Manuel Sánchez Guzmán has carried out research on many aspects of bird ecology and conservation biology, in particular on waterbirds. Alejandra Bettina Perales Casildo uses multivariate analysis to study socioeconomic and ecological data simultaneously. Óscar Uceda Tolosa specializes in the use of geographical information systems for spatial analysis in conservation biology.
