Figures
Abstract
Species adapted to cold-climatic mountain environments are expected to face a high risk of range contractions, if not local extinctions under climate change. Yet, the populations of many endothermic species may not be primarily affected by physiological constraints, but indirectly by climate-induced changes of habitat characteristics. In mountain forests, where vertebrate species largely depend on vegetation composition and structure, deteriorating habitat suitability may thus be mitigated or even compensated by habitat management aiming at compositional and structural enhancement. We tested this possibility using four cold-adapted bird species with complementary habitat requirements as model organisms. Based on species data and environmental information collected in 300 1-km2 grid cells distributed across four mountain ranges in central Europe, we investigated (1) how species’ occurrence is explained by climate, landscape, and vegetation, (2) to what extent climate change and climate-induced vegetation changes will affect habitat suitability, and (3) whether these changes could be compensated by adaptive habitat management. Species presence was modelled as a function of climate, landscape and vegetation variables under current climate; moreover, vegetation-climate relationships were assessed. The models were extrapolated to the climatic conditions of 2050, assuming the moderate IPCC-scenario A1B, and changes in species’ occurrence probability were quantified. Finally, we assessed the maximum increase in occurrence probability that could be achieved by modifying one or multiple vegetation variables under altered climate conditions. Climate variables contributed significantly to explaining species occurrence, and expected climatic changes, as well as climate-induced vegetation trends, decreased the occurrence probability of all four species, particularly at the low-altitudinal margins of their distribution. These effects could be partly compensated by modifying single vegetation factors, but full compensation would only be achieved if several factors were changed in concert. The results illustrate the possibilities and limitations of adaptive species conservation management under climate change.
Citation: Braunisch V, Coppes J, Arlettaz R, Suchant R, Zellweger F, Bollmann K (2014) Temperate Mountain Forest Biodiversity under Climate Change: Compensating Negative Effects by Increasing Structural Complexity. PLoS ONE 9(5): e97718. https://doi.org/10.1371/journal.pone.0097718
Editor: Jin Chen, Key Laboratory of Tropical Forest Ecology, Xishuangbanna Tropical Botanical Garden, Chinese Academy of Sciences, China
Received: January 10, 2014; Accepted: April 22, 2014; Published: May 13, 2014
Copyright: © 2014 Braunisch et al. This is an open-access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.
Funding: The work was funded by the Baden-Württemberg Fund for Environmental Research (http://www.mlr.baden-wuerttemberg.de/Forschung/68711.html) (Germany) and by the research programme “Forest and Climate Change” of the Swiss Federal Institute for Forest, Snow, and Landscape Research WSL (http://www.wsl.ch/index_EN) and the Federal Office for the Environment FOEN (http://www.bafu.admin.ch/?lang=en) (Switzerland). The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.
Competing interests: The authors have declared that no competing interests exist.
Introduction
With a predicted global temperature increase of 2.0–4.5°C until the end of the century (IPPC 2007), climate change is expected to affect habitat quality and species distributions [1]. Impacts have been demonstrated for all continents and taxonomic groups [2], [3], however, geographically isolated species adapted to cold climatic conditions [4] face a particularly high risk of range contractions, if not local extinction [5], [6]. In Europe, adverse effects are therefore mainly predicted for boreo-alpine taxa of mountain ecosystems [7], often being glacial relicts occurring at the margins of their eco-climatic niche [8], [9]. Species range-shifts are usually predicted based on large-scale species distribution models [10], describing species presence as a function of current climatic variation as well as coarse-grained, area-wide available environmental data [11], [12]. Yet, the populations of many endothermic species may not be primarily affected by physiological constraints of climate warming, but indirectly by climate-induced changes in habitat quality, food availability or interspecific interactions [1]. Consequently, the validity of predictions merely relying on climate functions may be questioned.
In forest ecosystems, biodiversity largely depends on the diversity of forest composition and structure (e.g., variability in tree species composition, vertical and horizontal forest structure, age structure of the stands, presence of gaps, clearings, snags and dead wood) [13], [14]. Next to the site conditions [15], the structural characteristics of montane and subalpine forests are mainly attributed to cold ambient temperatures which entail low forest productivity, long succession cycles and a high potential for snow-break or wind-throw with subsequent susceptibility for insect calamities [16], [17]. Although natural stand dynamics and resulting structural attributes are largely overruled by forestry, climate change is expected to affect forest vegetation composition and structure, and consequently, habitat suitability and distribution of the associated species. Adverse effects may therefore be additionally amplified by forestry practices aimed at coping with the economic risks of climate change such as the shortening of harvesting periods or changes in the tree species portfolio. On the other hand, a species’ dependence on vegetation characteristics may also offer the opportunity to counter negative effects of climate change by targeted habitat management (e.g., by increasing particular, species-relevant structural elements or vegetation components). We tested this option using the example of four mountain bird species of conservation concern: capercaillie (Tetrao urogallus), hazel grouse (Bonasa bonasia), three-toed woodpecker (Picoides tridactylus) and pygmy owl (Glaucidium passerinum). These species have been proposed as indicators for different, complementary forest structural attributes, and represent different niche dimensions within the mountain forest ecosystem. Capercaillie and three-toed woodpecker are additionally regarded as umbrella species for the associated ecological communities [18]–[20], thus supporting our aim to evaluate management measures that may support a wider range of biodiversity in mountain forests.
The model species show a high degree of specialization which facilitates tracking their responses to vegetation structures and climate-related variation thereof. The capercaillie is considered as an indicator for structurally rich, boreal and mountain forest habitats [21]–[23]. These habitats are characterised by an intermediate canopy cover, high proportions of old and open stands, and abundant ground vegetation – ideally dominated by bilberry (Vaccinium myrtillus) [21], [24], [25]. Similar to capercaillie, the hazel grouse requires structurally rich stands [26], [27], but prefers younger successional stages with sufficient berry or catkin bearing trees and shrubs [26]–[30]. A dense understory of shrubs and herbs further provides summer foraging habitat and cover from predators for both ground-nesting grouse species [26], [27], [29], [31], [32]. By excavating cavities, the three-toed woodpecker provides breeding opportunities for a variety of cavity-breeding birds and bats [33] and is therefore considered a key-stone species [18], [34]. It mainly feeds on the larvae of bark and wood-boring insects, predominantly found in dying and dead conifer (mostly spruce) trees [35]–[37]. Dead trees, snags and dying trees are therefore one of the most important habitat features for foraging [33], [36], [38]–[42]. The pygmy owl is the smallest avian predator in European boreal and mountain forests [43], [44]. It hunts small mammals as well as birds, insects and reptiles [45], [46] and uses cavities created by woodpeckers as nesting places as well as to hoard food [46]–[48]. A combination of dense young stands with high cover and open old forest with some small clearings is considered good breeding habitat [49]. Inner forest edges and edges between successional stages are often used for hunting [49], [50].
All four species are listed in Annex 1 of the European Birds directive [51] and are thus frequently targeted by conservation and restoration programmes. With climate change, there is not only an emergent risk that the benefits of these programmes will be curtailed; the prevailing predictions of range contractions and local extinctions have also led to a general debate in conservation management and policy that fundamentally questions the effectiveness and possibility of preserving climatically vulnerable species in their current habitats [52], [53].
We address these questions by assessing (1) how the occurrence of our model species’ depends on climate, landscape and vegetation characteristics, (2) how climate change and associated vegetation changes will affect overall habitat suitability, and (3) if decisive habitat features could be modified by adaptive management in a way that negative effects of climate change could be mitigated or compensated. While being aware of the high susceptibility of climate-change-related forecasts to various sources such as variations in climate change scenarios [54], statistical methods [55], [56] and model parameterizations [57], [58], which we have evaluated for our model species in an earlier study [47], we do not aim to provide absolute measures of habitat suitability and their changes. Rather, focusing on one method and scenario of climate change as an example, we aim to provide rough estimates for the magnitude of both effects and management efforts that would be necessary to preserve the model species’ in their Central European mountain habitats, thereby evaluating the general possibilities and limitations of adaptive conservation management in mountain forest environments under climate change.
Methods
Ethics Statement
Species data were adopted from existing databases, thus no mapping or handling of endangered species was involved. Vegetation mapping was mostly conducted in state and public forests where no permits for were required. Access to communal and private forest in Germany was covered by the Federal Forest Law of Baden-Württemberg (LWaldG §74 section 1), which allows entering private property for research purposes; in Switzerland similar rights were given by the Swiss Forest Law (WaG, article 14, §1). Permits for vegetation mapping within protected areas in Baden-Württemberg, Germany were issued by the Regional Council Freiburg (Regierungspräsidium Freiburg) department for nature conservation, in Switzerland permits were given by the cantonal departments of forestry and the Swiss National Park administration. The coordinates of the study locations are provided in Table S1, with the grid cells entirely or partly located in protected areas indicated and the approving authorities specified.
Study Area
The study area encompassed four mountain regions in Switzerland and Southern Germany with sympatric occurrence of the four model species, representing a broad gradient as regards climatic, vegetation and land-use conditions. The Black Forest, expanding over 7′000 km2 in Southwestern Germany, is a mainly forested lower mountain range with elevations ranging from 120 to 1′493 m a.s.l (mean: 663). The Swiss Jura, 4′200 km2 in size, is located in Western Switzerland and covers an altitudinal range between 500–1′718 m a.s.l (mean: 817). The Swiss Alps are here represented by two climatically and geographically distinct study regions: the “Northern Prealps”, defined by the biogeographic regions Prealps and Northern Alps with altitudes between 370–4′227 m a.s.l. (mean 1′391), and the Eastern Central Alps, with altitudes from 560–4′010 m a.s.l. (mean: 2′112) [59] (Figure 1). In the Black Forest and the Swiss Jura, where elevations do not reach the tree line, the forests form semi-continuous habitats interspersed by pasture land, while in the Northern Prealps forests surround treeless mountain tops. Finally, in the Eastern Central Alps, forests form distinct belts around high elevation peaks. Forest composition also varies along the altitudinal and climatic gradient, with decreasing proportions of European beech (Fagus sylvatica) and silver fir (Abies alba) giving way to a predominance of Norway spruce (Picea abies) when moving from the submontane to the subalpine belt. Moreover, larger proportions of larch (Larix decidua) and Swiss stone pine (Pinus cembra) can be found towards the Eastern Central Alps where a continental climate prevails in contrast to the other three regions, which are characterized by more oceanic climate conditions.
Within each grid cell, environmental variables were recorded at or in the surrounding of maximum 16 regularly distributed sampling plots (b), with only plots located in the forest (dark grey) considered in the analysis. At each plot, vegetation variables were recorded in the field at different reference areas (c), either across the whole plot (30×30 m), within a nested square (15×15 m), or within the two diagonal quarters of which (7.5×7.5 m). The variables and the reference area at which they were recorded are specified in Table 2. Geodata: Switzerland: © Bundesamt für Landestopografie Swisstopo (Art. 30 GeoIV): License No.: 5704 000 000, Available at: http://www.swisstopo.admin.ch/internet/swisstopo/en/home/products/height/dhm25.html; Germany: © Landesamt für Geoinformation und Landentwicklung Baden-Württemberg (LGL), License No.: 2851.9-1/19, Avaliable at: http://www.lgl-bw.de/lgl-internet/opencms/de/07_Produkte_und_Dienstleistungen/Geodaten/Digitale_Gelaendemodelle.
Species Data
Data of species presence were adopted from two databases hosted by the Swiss Ornithological Institute, Sempach, Switzerland (http://www.ornitho.ch) and the Forest Research Institute of Baden-Württemberg (FVA), Germany (http://www.wildtiermonitoring.de). Both databases contain long-term collections of observation data from ornithologists, foresters, hunters, birdwatchers as well as research personnel at a minimum resolution of 1 km2. Since data were not sampled systematically, no proven absence data were available.
In each of the four study regions we selected at least ten 1 km2 grid cells for each of the four focal species with species observations in at least three years between 2006 and 2010 (Table 1, Figure 1). Presence cells were selected by a stratified random process so as to represent the extent of the species distribution and its climatic gradient in the respective study region, thereby preferring cells with repeated observations from multiple years. For each presence grid cell, a corresponding cell in the following referred to as “absence cell” was selected, with “absence” defined as cells with no recorded species proof within the preceding 11 years (2000–2010). Absence cells were selected within a maximum of 5 km distance to the presence cell by randomly choosing one of the surrounding cells with at least 50% forest cover, while excluding all cells directly adjacent to the presence cell. With this we ensured that absence cells were located within the species’ dispersal ranges and did not expand too far beyond the limits of their altitudinal-climatic range so as to avoid trivial results and unsubstantiated extrapolations. We used a similar number of grid cell-pairs for each species, yet, since the species were not equally distributed across the study region, the relative numbers and the spatial distributions of cells differed in the four study regions (Table 1, Figure S1 a–d).
Environmental Variables
Sampling scheme.
Environmental predictors were sampled at 16 sampling plots, regularly distributed within each grid cell, with only plots located in the forest considered for the analysis (Figure 1). Our predictor set included variables of three main classes: climate, landscape and vegetation, measured at different reference areas around each sampling plot (Table 2).
Climate.
Climate variables included the average temperature in the breeding season (May–July) and in winter (December–February), and the sum of precipitation in both periods (Table 2). Current climate (long-term averages from 1971 to 2000) was obtained from the worldclim-dataset [60] (http://www.worldclim.org), which was downscaled from a 1 km2 raster to a resolution of 100×100 m based on the SRTM-V4 digital elevation model and the method described in [61].
For future climate conditions in the year 2050 (long-term averages from to 2031 to 2050) we assumed the moderate IPCC emission scenario A1B. Variables were derived from the Global Circulation Model ECHAM5, which was downscaled using the CLM Regional Circulation Model of the Max Planck Institute (http://cera-www.dkrz.de). A resolution of 100×100 m was then obtained by adding the anomalies between current and future climate conditions, which were downscaled to 1 km2 using the change factor methodology [62] to the current baseline data. All climate data were processed and provided by the Research Unit ‘Landscape Dynamics’ of the Swiss Federal Research Institute WSL.
Landscape.
Landscape variables included information on topography, land cover and human infrastructure, obtained from different digital data sources (Table 2). Five topographical variables (slope, topographic position, eastness, northness and potential solar radiation) were derived from the digital elevation model (DEM) for each sampling plot. The topographic position index, calculated with the extension TPI 1.3a for ArcView 3.3 [63], qualifies a point’s position relative to the surrounding terrain, with negative values indicating exposed sites such as hilltops or ridges, and positive values representing depressions. The potential solar radiation [W*h/m2] in the breeding season was calculated according to Fu and Rich [64] using the function “area solar radiation” in ArcGIS 9.3 [65]. Land cover variables encompassed the proportion of forest cover, intensively and extensively used agricultural land and wetland (mires and other habitat types on wet soils), as well as the density of outer forest edges, within the surrounding area of 1 km2. Human infrastructure was represented by the density of trafficable roads per km2 and the distance to settlements.
Vegetation.
Vegetation variables included information on vegetation composition (tree species and ground vegetation), vegetation structure (related to stand structure and forest stand mosaic) and special habitat features or resources relevant to the focal species. Vegetation composition, stand structure and special resources were mapped in the field at the sampling plots, while information on the forest stand mosaic was derived from remote sensing data. For matter of precision, different variables were assessed at different reference areas around the sampling plot center: tree species composition, successional stage, vertical and horizontal stand structure and selected special features (e.g., basal-branched trees or snags) were recorded within squares of 30×30 m, whereas special resources like the number of rowans or lying dead wood were quantified within a nested square of 15×15 m, the two diagonal quarters of which (7.5×7.5 m) were used to assess the ground vegetation (Figure 1, Table 2).
We derived variables describing the forest stand mosaic based on first and last return Light Detection and Ranging (LiDAR) data. For the cells in the Black Forest we used the revised point clouds for both terrain and surface models, as described in Schleyer [66], for Switzerland the corresponding data were provided by Swisstopo (2011). MATLAB R2011a (Mathworks, Natick, Massachusetts, USA) routines [67] were used to obtain terrain-corrected vegetation heights at a resolution of 3×3 m, which resembles the crown projection of a small spruce tree. The normalized vegetation heights were interpolated to form a continuous canopy height model, which was classified into four height classes: non-forested areas, shrub layer (<5 m), midstory (5–15 m) and canopy layer (>15 m) (processing details are provided in [68]). We used the height classes to calculate structural metrics describing the proportion of each height class per 1 km2, the number of gaps, the length of edges between different height classes representing different ecotone-types, as well as the total edge length between all height classes which provided an index for overall canopy height heterogeneity (for details see Table 2). Stand mosaic metrics were calculated in FRAGSTATS [48].
Statistical Analysis
Species occurrence.
We modelled species presence as a function of the environmental variables recorded at the sampling plots using Mixed Effects Logistic Regression with the grid-cell pair, as our species-sampling unit, treated as a random effect to account for spatial clustering. To identify the variables that best explained species presence we applied an information-theoretic approach [69], [70] using Akaike’s Information Criterion (AIC) to identify the most parsimonious model.
We followed a hierarchical variable selection procedure: first, univariate models were run for each variable, testing also the quadratic term for variables for which we expected a unimodal response. Of pairs of correlated variables (Spearman’s r ≥|0.6|) significantly contributing to explaining species presence in the univariate models we discarded the least performing one.
The retained predictors were then grouped into ecologically meaningful variable subsets (Figure 2, Table 2). For each subset a model was calibrated by testing all possible variable combinations and identifying the most parsimonious model using the R-package MuMIN. The variables that significantly contributed to this “best” subset-model were used for calibrating the model at the next hierarchy-level. This way the variable set was refined in a stepwise fashion, until a final model was obtained, potentially containing variables of all variable classes. The models’ fit was evaluated using multiple evaluation metrics, i.e., sensitivity, specificity, the percent correctly classified and Cohen’s Kappa at the optimal threshold, as well as the area under the receiver operating characteristics (ROC) curve (AUC).
The variables that significantly contributed to the most parsimonious model were retained for model calibration at the next hierarchy-level until a final model was obtained.
Climate-related vegetation trends.
To detect and quantify relationships between vegetation variables and climate, and obtain rough estimates of the magnitude and direction of possible vegetation changes under climate change, we used the data of all sampling plots mapped during the study (N = 4752) for modelling the species-relevant vegetation variables as functions of climate using Multiple Linear Regression. Only uncorrelated (Spearman’s r≥|0.6|) climate predictors that significantly contributed to explaining the respective vegetation characteristics were included. We applied these models to future climate conditions to calculate the mean expected change for each variable, which was then used to modify the respective variable value at each sampling plot to simulate vegetation changes under climate change.
Predicted changes in habitat suitability.
The bird species’ models were extrapolated to future climate conditions and the predicted change in occurrence probability (as a surrogate for habitat suitability) was quantified. We first calculated the change in occurrence probability due to climate change alone holding the vegetation variables constant, and secondly, also included the modelled climate-induced changes of the vegetation.
Compensation potential.
To evaluate the possibility to compensate for negative effects of climate change through habitat management, we selected for each species those among the significant vegetation variables that could potentially be modified by forest management. We then predicted species occurrence in dependence of these variables under current and future climate conditions while holding all other variables constant at their sampling average. Thereby each variable was only allowed to vary within the range of the empirical sampling values so as to avoid unfounded extrapolations outside the actually observed conditions. The comparison of the two response curves under current and future conditions illustrates the magnitude of variable modification that would be necessary for maintaining the current probability of species occurrence under the selected scenario of climate change. For each variable we then calculated the “compensation potential”, which was defined as the maximally achievable increase in predicted probability of species presence under altered climate conditions, which could be obtained when modifying the respective variable. Finally, to illustrate the effect of combining different measures of structural enhancement, we simultaneously varied the two variables for which the compensation potential was highest.
Results
Species’ Occurrence
The final models performed good to excellent [71] in predicting the presence of the model species. Accuracy was highest for the pygmy owl (AUC: 0.947, SD: 0.005) and lowest for the three-toed woodpecker (AUC: 0.877, SD: 0.010). The models for capercaillie and hazel grouse also demonstrated an excellent fit (AUC: 0.931 and 0.918, SD: 0.006 and 0.008, respectively). Further evaluation results are given in Table S2.
The final models of all four species contained variables of all main variable groups (climate, landscape, vegetation) (Table 3, Table S3). While all species showed a similar habitat selection pattern regarding climate, greater divergence was found for the decisive landscape variables and a different, even complementary set of vegetation and forest structural variables was retained for the different species (Table 3, Table S3). All species showed a unimodal response to winter temperature and a positive correlation with precipitation in early summer; in pygmy owl areas with higher winter precipitation were also selected. Concerning landscape characteristics, all species, except hazel grouse, preferred mires and forests on wet soils and showed at least a trend to avoid forests with a high road density and located in the vicinity of settlements. A negative response was also found for capercaillie towards a high density of outer forest edges, which can serve as an indicator for forest fragmentation. The presence of hazel grouse and three-toed woodpecker was negatively affected by the proportion of intensively managed agricultural land in the surroundings, while pygmy owl showed a positive correlation.
Except for the proportion of forest patches of the highest height class, which showed a quadratic relationship for capercaillie and was positively correlated with the presence of the three other species, the retained vegetation variables varied greatly between species’ models. Capercaillie presence was mainly explained by the abundance of gaps per km2, a low to moderate proportion of beech in the canopy and high cover of Vaccinium sp., mainly bilberry (Vaccinium myrtillus) in the field layer, as well as by low canopy height heterogeneity. Hazel grouse presence was mainly related to the availability of food sources, i.e. the proportion of resource trees and a high cover of herbs and bilberry; features providing cover, like basal-branched trees and a patchy ground vegetation distribution were preferred, while the vicinity of outer forest edges was avoided. Three-toed woodpecker occurrence was positively correlated with the presence of conifers and resource trees, and a high abundance of snags; while two- or multi-layered stands and stands with a high shrub cover were avoided. As in capercaillie, woodpecker presence was also negatively correlated with the abundance of hard stumps, indicating recent harvesting activities. Finally, pygmy owl habitat was characterized by a greater density of “sharp” forest edges: a greater abundance of basal-branched trees and higher ground vegetation cover than in locations where the species was absent.
Climate-related Vegetation Trends
Although most study sites were located in managed forests, all species-relevant vegetation variables were significantly correlated with climate, which explained between 2% (percentage of forest with canopy height >15 m) and 21% (percentage of pine, Pinus sp.) of the variation in the vegetation variables (Table S4). Climate change was predicted to have a negative effect on most of the vegetation variables with regard to their impact on the focal species (Table 4). While the models suggested a reduction in coniferous tree species and resource trees, the proportion of beech was predicted to increase. A decrease was also predicted for the abundance of gaps and the density of inner forest edges, which would go along with a reduction in ground vegetation cover and basal-branched trees.
Predicted Changes in Habitat Suitability
Climate change was predicted to negatively affect the probability of occurrence of all model species (Table 5). When considering only climatic changes, the greatest impact was predicted for hazel grouse, amounting to a reduction of presence probability of 29% in the currently occupied sites. The least-affected was the three-toed woodpecker with −22%, while capercaillie (−27%) and pygmy owl (−24%) were in intermediate positions (Table 5). Yet, when also considering climate-related vegetation changes, a significant additional reduction of presence probability of 14% was recorded for capercaillie while the conditions for the other three species remained more or less constant. Predicted climate change effects differed greatly between the four study regions, with the greatest impacts on all species recorded for the Black Forest. The Central Eastern Alps were the least affected by climate-change related habitat alterations, except for capercaillie, which was predicted to suffer least in the Swiss Jura and the Northern Prealps (Table S5).
Compensation Potential
The mean compensation potential, defined as the maximum increase in presence probability (ΔP(presence)) under the selected scenario of climate change, ranged between 0.02 (95% confidence interval CI: 0–0.05) for ground vegetation cover (GVCOV) and 0.72 (0.33–0-93) for the density of sharp edges (ED134), both in pygmy owl (Table 6). Adverse effects of climate change on capercaillie could be compensated best by increasing the number of gaps (GAPINDEX) from zero to 28 per km2, while hazel grouse availed most of an increase in bilberry (VAC) and resource trees (RESTREE). Increasing the number of snags (STANDDEAD) and the proportion of canopy heights >15 m (CHEIGHT4) most benefitted the tree-toed woodpecker. However, the comparison of the target species’ response curves under current and future climate conditions also showed that it was difficult or even impossible to maintain the prevailing occurrence probability by modifying only one vegetation variable (Figure 3, Figure S2). A considerable increase could be achieved by changing more than one variable towards the species’ optimum (Figure 4): For capercaillie a ΔP(presence) of 0.65 could be achieved when both VAC and GAPINDEX were modified so as to reach their optimal values. The maximum compensation potential for hazel grouse reached 0.73 with optimal proportions of RESTREE and VAC, while combining a maximum ED134 with a high number of basal-branched trees (BBTREE) increased ΔP(presence) for pygmy owl to 0.82. The probability of three-toed woodpecker presence could be increased by 0.65 when changing CHEIGHT4 and STANDDEAD towards their recorded maximum and by 0.77 when the latter variable was in combination with no recent harvesting activities.
Dashed lines indicate the 95% confidence interval. Variable codes are given in Table 2, response curves for all relevant vegetation variables are provided in Figure S2–S5.
For variable codes see Table 2.
Discussion
Despite widespread calls for adapting biodiversity conservation to the predicted impacts of climate change, the majority of strategies described in the literature are general recommendations that do not specify how to implement them in a real-world context [72]. This may partially be due to the high level of uncertainty in predicting the effects of climate-change on species and habitats [73] and the consequential ecological and economic risks associated with investing in an uncertain outcome. Our study is an attempt to reinforce the general principal of “securing populations by intensive management” (Heller and Zavaleta 2009) by evaluating the general potential and specific measures for compensating adverse effects of climate change on selected forest bird species. Nevertheless, within the presented framework, our study is subjected to the same limitations and sources of uncertainty inherent in predictive species-habitat modelling: First, our models are correlative, and based on the premise of niche conservatism [74]–[76], i.e. they assume that contemporary species-habitat associations will remain unchanged under altered climate conditions and that future changes in biotic interactions such as interspecific competition or predator-prey relationships will not modify them in essence. Second, our outcome may be subject to a considerable level of uncertainty that can arise from variation in input data [54], [55], statistical methods [55], [56] or model parameterization [57], [58]. Since these aspects have been evaluated earlier [59], we deliberately restricted our analyses to one climate change scenario, one global and regional circulation model and one statistical modelling approach, inferring that methodological variations may have changed the absolute values but not the general direction of our outcomes.
Contemporary Species Habitat Relations
Forest biodiversity will substantially be influenced by the interactions of climate change, forest management and the response of individual species. Our species data stemmed from different regions, covering a broad gradient of climate conditions, landscape and land use characteristics as well as forest management regimes, representative for the model species’ Central European habitats and chosen so as to maximize the generality of the results. The relevant habitat variables identified for the four model species were thus largely in line with the findings of earlier studies. While there was a convergent selection pattern regarding climate and most landscape characteristics, different variables pertaining to vegetation structure and composition were important to the different species, supporting their hypothesized complementarity regarding the indicator function for forest structural attributes [19], [20]. Surprisingly, climate variables significantly contributed to explaining the occurrence of all species, despite the fact that our sampling strategy confined the selection of absence cells to the altitudinal-climatic range that could potentially be occupied by the species. This suggests that - in addition to climate-related vegetation composition and structure - climatic conditions per se, or climate-related impacts not captured in this study (e.g., interspecific competition, predator-prey relationships or parasite abundance), play a role in the species’ regional distribution pattern.
Climate Effects on Vegetation Characteristics
All vegetation variables were significantly correlated with climatic predictors, yet, the proportion of explained variance was very low in some variables. The strongest correlations were found for the relative abundance of tree species such as Fagus sylvatica and Pinus spp, as well as for ground vegetation cover which reflects sub-canopy light conditions and thus also serves as an indicator for canopy density. Whereas the competitiveness of tree species is directly determined by climatic conditions, differences in canopy and ground vegetation cover may be attributed to longer vegetation periods and accelerated tree-growth in the lower altitudes, and to a higher abundance of beech-dominated forests with naturally scarcer ground vegetation. Yet, although climate is an important driver for shaping vegetation patterns, some of the correlations may be spurious. Hard stumps indicating harvesting activities, or the abundance of snags might be more affected by site-accessibility and topographically divergent patterns of forest management than by climate [77].
The observed correlations support the indirect effect of climate on the model species showing that favored climatic conditions mostly go along with favorable vegetation conditions. They also suggest expected vegetation changes under altered climate conditions to be mostly to the disadvantage of the species. Consequently, not only vegetation-related habitat suitability for our model species is expected to decrease with climate change; also the effort necessary to maintain or improve vegetation characteristics of suitable forest stands will increase because habitat management will have to additionally compensate for adverse natural dynamics. Of course, our simplified approach neglects site conditions, past management regimes and disturbance events and can thus only indicate rough vegetation trends. This may be a reason why the predicted decline of hazel grouse occurrence might be too large. The species is considered to profit from more frequent natural disturbances (e.g., wind storms) under future climatic conditions [78]. Thus, more sophisticated models that incorporate the occurrence of natural disturbances would be desirable and increase the accuracy of models in the future.
Effects of Climate Change on Habitat Suitability
Occurrence probability, as a proxy for habitat suitability, was predicted to decrease with climate change in all four species, although the magnitude of change differed between species and regions. As expected, impacts were greatest in the Black Forest located at the edge of the species’ bioclimatic envelope [9], while habitats in the Central Eastern Alps were least affected - except for capercaillie, which seems to be the most sensitive and specialized among the four species [59]. Reduction in occurrence probability was mostly caused by changes in climate parameters. A further decrease was observed when additionally including climate-induced vegetation changes, although this was only significant for capercaillie. However, the relative importance of correlated predictors or predictor sets remains difficult to assess in a correlative modelling framework [79]. For disentangling climatic and climate-induced vegetation effects a causality-based, experimental approach may be required.
Although our predictions refer to 2050, a time-lag in the actual response can be expected due to the low reactivity, long generation times and low natural migration rates of tree species in a forest ecosystem [80]. Although, due to their high mobility, bird species are generally attributed a high capacity to track environmental changes, our model species show specific traits which might increase this time-lag even further. All four species are non-migratory, characterized by a high degree of specialization, longevity and site fidelity [81], [82], attributes that may impede rapid responses to habitat alterations [83]. The grouse species are additionally characterized by limited dispersal abilities which inhibits exchanges between the mountain ranges [84]–[86] and underlines the necessity for on-site conservation efforts if populations should be maintained.
Compensation Potential and Adaptive Conservation Management
Recommendations about how to define desired target states and suitable adaptation strategies for biodiversity conservation under climate change cover a wide range of different approaches. Heller and Zavaleta [74] subsumed them to three main strategies: anticipatory reserve selection to secure future hotspots of biodiversity, improvement of landscape connectivity to allow species to track climate change, and on-site management to either increase the resilience or the resistance of populations or ecosystems to climate change. Although the possibility to preserve climate vulnerable species in their current habitats is frequently questioned [52], [53], we show that intensive management, enhancing the species-relevant vegetation structures, offers the potential to compensate for adverse effects of climate change. Yet, for our study species, a full compensation was difficult to achieve and mostly required the modification of more than one habitat feature. Moreover, the “optimal habitat values”, i.e. the values that would be required to achieve maximum compensation, can rarely be reached under real-world conditions. Especially in count-variables, such as the abundance of snags (for three-toed woodpecker), basal-branched trees (for pygmy owl) or gaps (for capercaillie) the optimum (corresponding to the maximum), is largely determined by locally extremely high variable values. While patches with the observed maximum of 42 snags per 900 m2 (dbh>12 cm) may be locally beneficial for the three-toed woodpecker, these values must certainly not be achieved across the whole area [42]. Consequently, for defining reasonable target values for management, the spatial distribution and average abundance of key variables across the whole potential habitat area (sensu [16]) must be taken into consideration.
Furthermore, some variables are characterized by a trade-off between feature abundance and size: a high number of 28 gaps (<0.1 ha) per hectare may be possible if they are small, whereas fewer but larger gaps may also be sufficient for capercaillie. Suchant and Braunisch [87], [88] recommend a minimum of 10% of the forest area to be gaps, which might be a more appropriate prescription.
While most measures increasing structural diversity can readily be implemented, changes in forest composition require considerable time and effort. Converting the tree species portfolio towards more drought-adapted tree species is the most important strategy for maintaining forest productivity under climate change, particularly in regions with a pronounced legacy of former silviculture that promoted drought-intolerant tree species, e.g. spruce, outside their natural range. Whether this conversion will support or be in conflict with the habitat requirements of our target species will strongly depend on the selected tree-species: whereas replacing non site-adapted spruce by fir or pine may maintain or even enhance habitat suitability, a promotion of beech would contribute to habitat deterioration of capercaillie and hazel grouse as it would have a direct impact on the ground vegetation. This is particularly to be expected in the Black Forest habitats and at the lower edge of their distribution range.
To tap the full compensation potential, the combined effect of moderately modifying different variables has to be considered and measures must be flexibly adapted to the local site conditions. While gaps or edges may be created in all situations, the abundance of bilberry, for example, can only be increased under suitable soil and light conditions. Yet, despite the high efforts and inherent limitations, habitat improvement may be the favorable option in an uncertain future: simulation studies showed that habitat improvement led to higher species survival under climate change than creating new habitats in prospectively suitable locations [89]. Considering the high divergence between species range forecasts under climate change, enhancing structural complexity in currently occupied habitats represents a conservative “no-regret” strategy - particularly in forest ecosystems which are well known for their moderating effect on local climate conditions compared to open habitats [90].
Conclusions
Our study shows that intensive habitat management focusing on a relatively small set of decisive variables can buffer indirect negative effects of climate change on forest-dwelling species although it partly requires working against the natural dynamics. This raises several questions in a system of multifunctional forestry where adaptive conservation management has to be balanced with other ecosystem services [91]. First, which target values can be maximally and realistically achieved? Second, where to prioritize investments? And finally, how can long-term implementation be guaranteed in an ecologically and economically sustainable way?
The first aspect is mainly subject to societal values and political decisions [14]. Compensation measures can be costly and in conflict with economic goals and adaptive strategies to manage renewable resources under climate change. The long-term success of the measures taken will therefore depend on public acceptance and a cost-effective planning and integration in regular forest management. To achieve this, investments should be prioritized in areas large enough to support minimum viable populations of the target species, with a key function for population connectivity or functionality. Moreover, areas should be preferred where – based on the prevailing site- and stocking conditions – the expected climate-change impacts can be compensated with a justifiable management effort. Finally, enhancing structural diversity will not solely increase the model species’ resistance towards climate change. Diverse forests with mixed stands, providing multiple niches for both native and immigrating species, are considered a major prerequisite for ecosystem resilience [92]. Measures aiming at preserving indicators of structural diversity may therefore be beneficial to a wide range of taxa of the representative species community, even when the target species may decline or finally go extinct.
Supporting Information
Figure S1.
Distribution of species data. 1 km grid cells with presence (white) and absence (black) of (a) capercaillie, (b) hazel grouse, (c) three-toed woodpecker and (d) pygmy owl in the four study regions Black Forest (BF), Swiss Jura (J), Northern and Prealps (NPA), and Central Alps (CEA).
https://doi.org/10.1371/journal.pone.0097718.s001
(PDF)
Figure S2.
Predicted probability of capercaillie presence for current (2010, black) and future (2050, grey) climate conditions (a–d). Presence probability was modeled in dependence of species-relevant vegetation variables, while holding all other variables at their empirical sampling average. For variable codes see Table 2.
https://doi.org/10.1371/journal.pone.0097718.s002
(PDF)
Figure S3.
Predicted probability of hazel grouse presence for under current (2010, black) and future (2050, grey) climate conditions (a–d). Presence probability was modeled in dependence of species-relevant vegetation variables, while holding all other variables at their empirical sampling average. For variable codes see Table 2.
https://doi.org/10.1371/journal.pone.0097718.s003
(PDF)
Figure S4.
Predicted probability of three-toed woodpecker presence for under current (2010, black) and future (2050, grey) climate conditions (a–f). Presence probability was modeled in dependence of species-relevant vegetation variables, while holding all other variables at their empirical sampling average. For variable codes see Table 2.
https://doi.org/10.1371/journal.pone.0097718.s004
(PDF)
Figure S5.
Predicted probability of pygmy owl presence for under current (2010, black) and future (2050, grey) climate conditions (a–d). Presence probability was modeled in dependence of species-relevant vegetation variables, while holding all other variables at their empirical sampling average. For variable codes see Table 2.
https://doi.org/10.1371/journal.pone.0097718.s005
(PDF)
Table S1.
Study locations in the four study regions Black Forest (BF), Swiss Jura (J), Northern Prealps (NPA) and Central Eastern Alps (CEA). Grid cells (1 km2) are represented by their centroid, with the location given in DHDN/3-degree Gauss-Kruger zone 3 (GAUSS) and in the Swiss coordinate system CH1903 (SG). Grid cells entirely or partly located within protected areas without public access and the authority that issued the permit for vegetation mapping are indicated.
https://doi.org/10.1371/journal.pone.0097718.s006
(PDF)
Table S2.
Accuracy of the models for capercaillie (CC), hazel grouse (HG), three-toed woodpecker (TTW) and pygmy owl (PO). Model fit is indicated by sensitivity, specificity, the percent correctly classified (PCC) and Cohen’s Kappa (κ_max) at its optimal threshold, as well as the area under the receiver operating characteristics curve (AUC).
https://doi.org/10.1371/journal.pone.0097718.s007
(PDF)
Table S3.
Final models for (a) Capercaillie, (b) Hazel grouse, (c) Three-toed woodpecker and (d) Pygmy owl. The codes for retained variables of the main variable categories C = climate, L = landscape and V = vegetation are provided in Table 2. The variables that were tested for their compensation potential (i.e. that could be modified by forest management so as to increase the probability of species presence under climate change) are indicated by asterisks. For variable codes see Table 2.
https://doi.org/10.1371/journal.pone.0097718.s008
(PDF)
Table S4.
Multiple linear regression models describing the correlation of vegetation variables selected in the species models as a function of climate variables. Models were calculated across all sampling plots in the study area. For variable codes see Table 2.
https://doi.org/10.1371/journal.pone.0097718.s009
(PDF)
Table S5.
Modelled probability of species presence (Ppres) at the presence plots in the four study regions (Black Forest BF, Swiss Jura J, Northern Prealps NPA and Central Eastern Alps CEA), as well as mean predicted changes thereof (ΔPpres) between 2010 and 2050 under climate change. The first model considers only changes in climate variables (2050C), the second (2050CV) additionally takes predicted vegetation changes into account. CC: Capercaillie, HG: Hazel grouse, TTW: Three-toed woodpecker, PO: Pygmy owl.
https://doi.org/10.1371/journal.pone.0097718.s010
(PDF)
Appendix S1.
Sources of the geo-data used in this study. The listing corresponds to the superscripts provided in Table 2.
https://doi.org/10.1371/journal.pone.0097718.s011
(PDF)
Acknowledgments
Thanks to all persons who contributed to fieldwork: Lisa Bitterlin, Karin Feller, Lucrezia Deplazes, Lea Hofstetter, Nino Maag, Maria Rusche, Jan Sadowski and Sarah Spille. Niklaus Zimmermann and Achilleas Psomas processed and provided the climate data. We thank Felix Morsdorf for his support with processing the LiDAR data, Hans Schmid for extracting the species records from the Swiss ornithological institute’s database, Frank Hohlfeld, Gerrit Müller, Christoph Purschke and Laurent Willenegger for providing additional species data and Davide De Masi for language proofreading.
Author Contributions
Conceived and designed the experiments: VB KB RA RS. Performed the experiments: VB JC FZ. Analyzed the data: VB JC FZ. Contributed reagents/materials/analysis tools: KB RS RA. Wrote the paper: VB KB RA.
References
- 1. Parmesan C (2006) Ecological and evolutionary responses to recent climate change. Annual Review of Ecology, Evolution and Systematics 37: 637–669.
- 2. Parmesan C, Yohe G (2003) A globally coherent fingerprint of climate change impacts across natural systems. Nature 421: 37–42.
- 3. Chen IC, Hill JK, Ohlemüller R, Roy DB, Thomas CD (2011) Rapid range shirts of species associated with high levels of climate warming. Science 333: 1024–1026.
- 4. Hulme PE (2005) Adapting to climate change: is there scope for ecological management in the face of a global threat? Journal of Applied Ecology 42: 784–794.
- 5. Hughes L (2000) Biological consequences of global warming: is the signal already apparent? Trends in Ecology and Evolution 15: 56–61.
- 6. Thomas CD, Cameron A, Green RE, Bakkenes M, Beaumont LJ, et al. (2004) Extinction risk from climate change. Nature 427: 145–148.
- 7. Sergio F (2003) Relationship between laying dates of black kites Milvus migrans and spring temperatures in Italy: rapid response to climate change? Journal of Avian Biology 34: 144–149.
- 8. Wilson JW, Gutierrez D, Martinez D, Agudo R, Monserrat VJ (2005) Changes to the elevational limits and extent of species ranges associated with climate change. Ecology Letters 8: 1138–1146.
- 9. Braunisch V, Bollmann K, Graf RF, Hirzel AH (2008) Living on the edge - modelling habitat suitability for species at the edge of their fundamental niche. Ecological Modelling 214: 153–167.
- 10. Guisan A, Zimmermann NE (2000) Predictive habitat distribution models in ecology. Ecological Modelling 135: 147–186.
- 11. Rebelo H, Tarroso P, Jones G (2010) Predicted impact of climate change on European bats in relation to their biogeographic patterns. Global Change Biology 16: 561–576.
- 12. Barbet-Massin M, Thuiller W, Jiguet F (2010) How much do we overestimate future local extinction rates when restricting the range of occurrence data in climate suitability models? Ecography 33: 378–386.
- 13. McElhinny C, Gibbons P, Brack C, Bauhus J (2005) Forest and woodland stand structural complexity: Its definition and measurement. Forest Ecology and Management 218: 1–24.
- 14. Bollmann K, Bergamini A, Senn-Irlet B, Nobis M, Duelli P, et al. (2009) Konzepte, Instrumente und Herausforderungen bei der Förderung der Biodiversität im Wald. Schweizerische Zeitschrift für Forstwesen 160: 53–67.
- 15. Braunisch V, Suchant R (2008) Using ecological forest site mapping for long-term habitat suitability assessments in wildlife conservation - demonstrated for capercaillie (Tetrao urogallus). Forest Ecology and Management 256: 1209–1221.
- 16. Braunisch V, Suchant R (2007) A model for evaluating the ‘habitat potential’ of a landscape for capercaillie Tetrao urogallus: a tool for conservation planning. Wildlife Biology 13: 21–33.
- 17. Bollmann K, Friedrich A, Fritsche B, Graf RF, Imhof S, et al. (2008) Kleinräumige Habitatnutzung des Auerhuhns im Alpenraum. Der Ornithologische Beobachter 105: 53–61.
- 18. Mikusiński G, Gromadzki M, Chylarecki P (2001) Woodpeckers as indicators of forest bird diversity. Conservation Biology 15: 208–217.
- 19. Suter W, Graf RF, Hess R (2002) Capercaillie (Tetrao urogallus) and Avian Biodiversity: Testing the Umbrella-Species Concept. Conservation Biology 16: 778–788.
- 20. Roberge J, Angelstam P (2006) Indicator species among resident forest birds – A cross-regional evaluation in northern Europe. Biological Conservation 130: 134–147.
- 21. Storch I (1993) Habitat selection by capercaillie in summer and autumn - is bilberry important? Oecologia 95: 257–265.
- 22.
Storch I (1995) Habitat requirements of capercaillie. 151–154.
- 23. Cas J, Adamic M (1998) The influence of forest alteration on the distribution of capercaillie leks in the Eastern Alps. Zbornik 57: 5–57.
- 24. Picozzi N, Catt DC, Moss R (1992) Evaluation of Capercaillie habitat. Journal of Applied Ecology 29: 751–762.
- 25. Selås V (2000) Population dynamics of capercaillie Tetrao urogallus in relation to bilberry Vaccinium myrtillus production in southern Norway. Wildlife Biology 6: 1–11.
- 26. Sachot S, Perrin N, Neet C (2003) Winter habitat selection by two sympatric forest grouse in western Switzerland: implications for conservation. Biological Conservation 112: 373–382.
- 27. Schäublin S, Bollmann K (2011) Winter habitat selection and conservation of Hazel Grouse (Bonasa bonasia) in mountain forests. Journal of Ornithology 152: 179–192.
- 28. Blattner M (1998) Der Arealschwund des Haselhuhns Bonasa bonasia in der Nordwestschweiz. Der Ornithologische Beobachter 95: 11–38.
- 29.
Bergmann HH, Klaus S, Müller F, Scherzinger W, Swenson JE, et al. (1996) Die Haselhühner - 4. überarbeitete Auflage. Magdeburg: Die neue Brehm-Bücherei Bd. 77. Westarp Wissenschaften.
- 30.
Swenson JE (1995) The ecology of Hazel Grouse and management of its habitat. 227–238 p.
- 31. Swenson JE (1993) The importance of alder to hazel grouse in Fennoscandian boreal forest: evidence from four levels of scale. Ecography 16: 37–46.
- 32.
Kaempfer-Lauenstein A (1995) Mehr Wildnis fuer das Haselhuhn! Nationalpark. 6–9 (in German).
- 33.
Pechacek P (1995) Spechte (Picidae) im Nationalpark Berchtesgaden: Habitatwahl, Nahrungsökologie, Populationsdichte. Forschungsbericht Nationalpark Berchtesgaden 31.
- 34. Pechacek P, d’Oleire-Oltmanns W (2004) Habitat use of the three-toed woodpecker in central Europe during the breeding period. Biological Conservation 116: 333–341.
- 35. Pechacek P, Kristin A (1993) Nahrung der Spechte im Nationalpark Berchtesgaden. Die Vogelwelt 114: 165–177.
- 36. Murphy EC, Lehnhausen WA (1998) Density and foraging ecology of woodpeckers following a stand-replacement fire. Journal of Wildlife Management 62: 1359–1372.
- 37. Fayt P (1999) Available insect prey in bark patches selected by the Three-toed Woodpecker Picoides tridactylus prior to reproduction. Ornis Fennica 76: 135–140.
- 38. Hogstad O (1970) On the ecology of the three-toed woodpecker Picoides tridactylus (L.) outside the breeding season. Nytt Magasin for Zoology 18: 221–227.
- 39. Hess R (1983) Verbreitung, Siedlungsdichte und Habitat des Dreizehenspechts Picoides tridactylus alpinus im Kanton Schwyz. Der Ornithologische Beobachter 80: 153–182.
- 40. Ruge K, Gorze HJ, Havelka P (1999) Der Dreizehenspecht (Picoides tridactylus) im Schwarzwald (Deutschland): Zusammenstellung brutbiologischer Daten. Tichodroma 12: 116–121.
- 41. Imbeau L, Desrochers A (2002) Foraging ecology and use of drumming trees by three-toed woodpeckers. Journal of Wildlife Management 66: 222–231.
- 42. Bütler R, Angelstam P, Ekelund P, Schlaepfer R (2004) Dead wood threshold values for the three-toed woodpecker presence in boreal and sub-Alpine forest. Biological Conservation 119: 305–318.
- 43.
Mikkola H (1983) Owls of Europe. Calton: T & A D Poyser. 397 p.
- 44. Strøm H, Sonerud GA (2001) Home range and habitat selection in the Pygmy Owl Glaucidium passerinum. Ornis Fennica 78: 145–158.
- 45. Mikusek R, Kloubec B, Obuch J (2001) Diet of the Pygmy Owl (Glaucidium passerinum) in eastern Central Europe. Buteo 12: 47–60.
- 46. Halonen M, Mappes T, Meri T, Suhonen J (2007) Influence of snow cover on food hoarding in Pygmy Owls Glaucidium passerinum. Ornis Fennica 84: 105–111.
- 47. Braunisch V, Coppes J, Schmid H, Suchant R, Arlettaz R, et al. (2013) Selecting from correlated climate variables: a major source of uncertainty for predicting species distributions under climate change. Ecography 36: 1–13.
- 48.
McGarigal K, Cushman SA, Neel MC, Ene E (2002) FRAGSTATS: Spatial Pattern Analysis Program for Categorical Maps. Amherst: University of Massachusetts.
- 49.
Glutz von Blotzheim UN, Bauer K, Bezzel E (1973) Handbuch der Vögel Mitteleuropas. Band 5 Galliformes. Frankfurt. 225 p.
- 50. Kullberg C (1995) Strategy of the pygmy owl while hunting avian and mammalian prey. Ornis Fennica 7: 72–78.
- 51.
EC (2010) The Birds Directive. Directive 2009/147/EC of the European Parliament and of the Council of 30 November 2009 on the conservation of wild birds.
- 52. Araújo MB, Cabeza M, Thuiller W, Hannah L, Williams PH (2004) Would climate change drive species out of reserves? An assessment of existing reserve selection methods. Global Change Biology 10: 1618–1626.
- 53. Araújo MB, Algador D, Cabeza M, Nogués-Bravo D, Thuiller W (2011) Climate change threatens European conservation areas. Ecology Letters 14: 484–492.
- 54.
Beaumont LJ, Pitman AJ, Poulsen M, Hughes L (2007) Where will species go? Incorporating new advances in climate modelling into projections of species distributions. Global Change Biology 13.
- 55. Thuiller W (2004) Patterns and uncertainties of species’ range shifts under climate change. Global Change Biology 10: 2020–2027.
- 56. Lawler JJ, White D, Neilson RP, Blaustein AR (2006) Predicting climate-induced range shifts: model differences and model reliability. Global Change Biology 12: 1568–1584.
- 57. Synes NW, Osborne PE (2011) Choice of predictor variables as a source of uncertainty in continental-scale species distribution modelling under climate change. Global Ecology and Biogeography 20: 904–914.
- 58. Beaumont LJ, Hughes L, Poulsen M (2005) Predicting species distributions: use of climatic parameters in BIOCLIM and its impact on predictions of species’ current and future distributions. Ecological Modelling 186: 250–269.
- 59.
Gonseth Y, Wohlgemuth T, Sansonnens B, Buttler A (2001) Die biogeographischen Regionen der Schweiz - Les régions biogéographiques de la Suisse. Bern/Berne: BUWAL/OFEFP. Nr./n° 137 Nr./n° 137. 47 p.
- 60. Hijmans RJ, Cameron SE, Parra JL, Jones PG, Jarvis A (2005) Very high resolution interpolated climate surfaces for global land areas. International Journal of Climatology 25: 1965–1978.
- 61.
Zimmermann NE, Roberts DW (2001) Final Report of the MLP climate and biophysical mapping project. Birmensdorf. 18 p.
- 62. Diaz-Nieto J, Wilby RL (2005) A comparison of statistical downscaling and climate change factor methods: impacts on low flows in the River Thames, United Kingdom. Climatic Change 69: 245–268.
- 63.
Jenness J (2006) Topographic Position Index (tpi_jen.avx) extension for ArcView 3.x, v. 1.2.: Jenness Enterprises.
- 64. Fu P, Rich PM (2002) A geometric solar radiation model with applications in agriculture and forestry. Computers and Electronics in Agriculture 37: 25–35.
- 65.
ESRI (2008) ArcGIS Desktop: Release 9.3. Redlands, CA: Environmental Systems Research Institute.
- 66.
Schleyer A (2001) Das Laserscan-DGM von Baden-Württemberg. In: Fritsch D, Spiller R, editors. Photogrammetric Week '01: Wichmann Verlag, Heildelberg.
- 67. Morsdorf F, Marell A, Koetz B, Cassagne N, Pimont F, et al. (2010) Discrimination of vegetation strata in a multi-layered Mediterranean forest ecosystem using height and intensity information derived from airborne laser scanning. Remote Sensing of Environment 114: 1403–1415.
- 68.
Zellweger F, Braunisch V, Baltensweiler A, Bollmann K (2013) Remotely sensed forest structural complexity predicts multi species occurrence at the landscape scale. Forest Ecology and Management (in press).
- 69.
Burnham KP, Anderson DR (2002) Model Selection and Multimodel Inference. New York: Springer-Verlag. 488 p.
- 70. Johnson JB, Omland KS (2004) Model selection in ecology and evolution. Trends in Ecology & Evolution 19: 101–108.
- 71.
Hosmer DWJ, Lemeshow S (2000) Applied Logistic Regression, 2nd edition. John Wiley & Sons, New York.
- 72. Heller NE, Zavaleta ES (2009) Biodiversity management in the face of climate change: A review of 22 years of recommendations. Biological Conservation 142: 14–32.
- 73. Kujala H, Burgman MA, Moilanen A (2013) Treatment of uncertainty in conservation under climate change. Conservation Letters 6: 73–85.
- 74. Pearson RG, Dawson TP (2003) Predicting the impacts of climate change on the distribution of species: are bioclimate envelope models useful? Global Ecology and Biogeography 12: 361–371.
- 75. Dormann CF (2007) Promising the future? Global change projections of species distributions. Basic and Applied Ecology 8: 387–397.
- 76. Wiens JJ, Ackerly DD, Allen AP, Anacker BL, Buckley LB, et al. (2010) Niche conservatism as an emerging principle in ecology and conservation biology. Ecology Letters 13: 1310–1324.
- 77. Fürstenau C, Badeck FW, Lasch P, Lexer MJ, Lindner M, et al. (2007) Multiple-use forest management in consideration of climate change and the interests of stakeholder groups. (COST Action E21 Contribution of forests and forestry to the mitigation of greenhouse effects.). European Journal of Forest Research 126: 225–239.
- 78.
Bollmann K (2010) Climate or land use change: what affects Alpine grouse species more? Biodiversity in time of climate change: management or wilderness? Chambéry: ALPARC. 28–32.
- 79. Meloun M, Militky J, Hill M, Brereton RG (2002) Crucial problems in regression modelling and their solutions. Analyst 127: 433–450.
- 80. Jump AS, Penuelas J (2005) Running to stand still: adaptation and the response of plants to rapid climate change. Ecology Letters 8: 1010–1020.
- 81. Rhim SJ, Lee WS (2001) Characteristics of hazel grouse Bonasa bonasia distribution in southern Korea. Wildlife Biology 7: 257–261.
- 82. Pechacek P (2006) Breeding performance, natal dispersal, and nest site fidelity of the three-toed woodpecker in the German Alps. Annales Zoologici Fennici 43: 165–176.
- 83. Angert AL, Crozier LG, Rissler LJ, Gilman SE, Tewksbury JJ, et al. (2011) Do species’ traits predict recent shifts at expanding range edges? Ecology Letters 14: 677–689.
- 84. Segelbacher G, Höglund J, Storch I (2003) From connectivity to isolation: genetic consequences of population fragmentation in capercaillie across Europe. Molecular Ecology 12: 1773–1780.
- 85. Braunisch V, Segelbacher G, Hirzel AH (2010) Modelling functional landscape connectivity from genetic population structure: a new spatially explicit approach. Molecular Ecology 19: 3664–3678.
- 86. Kormann U, Gugerli F, Ray N, Excoffier L, Bollmann K (2012) Parsimony-based pedigree analysis and individual-based landscape genetics suggest topography to restrict dispersal and connectivity in the endangered capercaillie. Biological Conservation 152: 241–252.
- 87. Suchant R, Braunisch V (2004) Multidimensional habitat modelling in forest management - a case study using capercaillie in the Black Forest, Germany. Ecological Bulletins 51: 455–649.
- 88.
Suchant R, Braunisch V (2008) Rahmenbedingungen und Handlungsfelder für den Aktionsplan Auerhuhn: Grundlagen für ein integratives Konzept zum Erhalt einer überlebensfähigen Auerhuhnpopulation im Schwarzwald. Freiburg.
- 89. Jeltsch F, Moloney KA, Schwager M, Körner K, Blaum N (2011) Consequences of correlations between habitat modifications and negative impact of climate change for regional species survival. Agriculture, Ecosystems and Environment 145: 49–58.
- 90. Renaud V, Innes JL, Dobbertin M, Rebetez M (2011) Comparison between open-site and below-canopy climatic conditions in Switzerland for different types of forests over 10 years (1998–2007). Theoretical and Applied Climatology 105: 119–127.
- 91.
Bollmann K, Braunisch V (2013) To integrate or to segregate: balancing commodity production and biodiversity conservation in European forests. In: Kraus D, Krumm F, editors. Integrative approaches as an opportunity for the conservation of forest biodiversity. Freiburg, Germany: European Forest Institute. 18–31.
- 92. Noss RF (2001) Beyond Kyoto: forest management in a time of rapid climate change. Conservation Biology 15: 578–590.