Reconstruction Of Past Climate Conditions Over Central Europe From Groundwater Data

Here we present a 30 000 years low-resolution climate record reconstructed from groundwater data. The investigated site is located in the Bohemian Cretaceous Basin, in the corridor between the Scandinavian ice sheet and the Alpine ice field. Noble gas temperatures (NGT), obtained from groundwater data, preserved multicentennial temperature variability and indicated a cooling of at least 5-7 °C during the last glacial maximum (LGM). This is further confirmed by the depleted δ18O and δ2H values at the LGM. High excess air (ΔNe) at the end of the Pleistocene is possibly related to abrupt changes in recharge dynamics due to progression and retreat of ice covers and permafrost. These results agree with the fact that during the LGM permafrost and small glaciers developed in the inner valleys of the Giant Mountains (located in the watershed of the aquifers). A temporal decrease of deuterium excess from the pre-industrial Holocene to present days is linked to an increase of the air temperatures, and probably also to an increase of water pressure at the source region of precipitation over the past few hundred years

Molecular Proxies for Climate Maladaptation in a Long-Lived Tree (Pinus pinaster Aiton, Pinaceae).

Understanding adaptive genetic responses to climate change is a main challenge for preserving biological diversity. Successful predictive models for climate-driven range shifts of species depend on the integration of information on adaptation, including that derived from genomic studies. Long-lived forest trees can experience substantial environmental change across generations, which results in a much more prominent adaptation lag than in annual species. Here, we show that candidate-gene SNPs (single nucleotide polymorphisms) can be used as predictors of maladaptation to climate in maritime pine (Pinus pinaster Aiton), an outcrossing long-lived keystone tree. A set of 18 SNPs potentially associated with climate, 5 of them involving amino acid-changing variants, were retained after performing logistic regression, latent factor mixed models, and Bayesian analyses of SNP-climate correlations. These relationships identified temperature as an important adaptive driver in maritime pine and highlighted that selective forces are operating differentially in geographically discrete gene pools. The frequency of the locally advantageous alleles at these selected loci was strongly correlated with survival in a common garden under extreme (hot and dry) climate conditions, which suggests that candidate-gene SNPs can be used to forecast the likely destiny of natural forest ecosystems under climate change scenarios. Differential levels of forest decline are anticipated for distinct maritime pine gene pools. Geographically defined molecular proxies for climate adaptation will thus critically enhance the predictive power of range-shift models and help establish mitigation measures for long-lived keystone forest trees in the face of impending climate change.

Modelling tree population dynamics at the alpine and boreal tree-line ecotones in response to climate and land-use change

Summary Ecotones are sensitive to change because they contain high numbers of species living at the margin of their environmental tolerance. This is equally true of tree-lines, which are determined by attitudinal or latitudinal temperature gradients. In the current context of climate change, they are expected to undergo modifications in position, tree biomass and possibly species composition. Attitudinal and latitudinal tree-lines differ mainly in the steepness of the underlying temperature gradient: distances are larger at latitudinal tree-lines, which could have an impact on the ability of tree species to migrate in response to climate change. Aside from temperature, tree-lines are also affected on a more local level by pressure from human activities. These are also changing as a consequence of modifications in our societies and may interact with the effects of climate change. Forest dynamics models are often used for climate change simulations because of their mechanistic processes. The spatially-explicit model TreeMig was used as a base to develop a model specifically tuned for the northern European and Alpine tree-line ecotones. For the latter, a module for land-use change processes was also added. The temperature response parameters for the species in the model were first calibrated by means of tree-ring data from various species and sites at both tree-lines. This improved the growth response function in the model, but also lead to the conclusion that regeneration is probably more important than growth for controlling tree-line position and species' distributions. The second step was to implement the module for abandonment of agricultural land in the Alps, based on an existing spatial statistical model. The sensitivity of its most important variables was tested and the model's performance compared to other modelling approaches. The probability that agricultural land would be abandoned was strongly influenced by the distance from the nearest forest and the slope, bath of which are proxies for cultivation costs. When applied to a case study area, the resulting model, named TreeMig-LAb, gave the most realistic results. These were consistent with observed consequences of land-abandonment such as the expansion of the existing forest and closing up of gaps. This new model was then applied in two case study areas, one in the Swiss Alps and one in Finnish Lapland, under a variety of climate change scenarios. These were based on forecasts of temperature change over the next century by the IPCC and the HadCM3 climate model (ΔT: +1.3, +3.5 and +5.6 °C) and included a post-change stabilisation period of 300 years. The results showed radical disruptions at both tree-lines. With the most conservative climate change scenario, species' distributions simply shifted, but it took several centuries reach a new equilibrium. With the more extreme scenarios, some species disappeared from our study areas (e.g. Pinus cembra in the Alps) or dwindled to very low numbers, as they ran out of land into which they could migrate. The most striking result was the lag in the response of most species, independently from the climate change scenario or tree-line type considered. Finally, a statistical model of the effect of reindeer (Rangifer tarandus) browsing on the growth of Pinus sylvestris was developed, as a first step towards implementing human impacts at the boreal tree-line. The expected effect was an indirect one, as reindeer deplete the ground lichen cover, thought to protect the trees against adverse climate conditions. The model showed a small but significant effect of browsing, but as the link with the underlying climate variables was unclear and the model was not spatial, it was not usable as such. Developing the TreeMig-LAb model allowed to: a) establish a method for deriving species' parameters for the growth equation from tree-rings, b) highlight the importance of regeneration in determining tree-line position and species' distributions and c) improve the integration of social sciences into landscape modelling. Applying the model at the Alpine and northern European tree-lines under different climate change scenarios showed that with most forecasted levels of temperature increase, tree-lines would suffer major disruptions, with shifts in distributions and potential extinction of some tree-line species. However, these responses showed strong lags, so these effects would not become apparent before decades and could take centuries to stabilise. Résumé Les écotones son sensibles au changement en raison du nombre élevé d'espèces qui y vivent à la limite de leur tolérance environnementale. Ceci s'applique également aux limites des arbres définies par les gradients de température altitudinaux et latitudinaux. Dans le contexte actuel de changement climatique, on s'attend à ce qu'elles subissent des modifications de leur position, de la biomasse des arbres et éventuellement des essences qui les composent. Les limites altitudinales et latitudinales diffèrent essentiellement au niveau de la pente des gradients de température qui les sous-tendent les distance sont plus grandes pour les limites latitudinales, ce qui pourrait avoir un impact sur la capacité des espèces à migrer en réponse au changement climatique. En sus de la température, la limite des arbres est aussi influencée à un niveau plus local par les pressions dues aux activités humaines. Celles-ci sont aussi en mutation suite aux changements dans nos sociétés et peuvent interagir avec les effets du changement climatique. Les modèles de dynamique forestière sont souvent utilisés pour simuler les effets du changement climatique, car ils sont basés sur la modélisation de processus. Le modèle spatialement explicite TreeMig a été utilisé comme base pour développer un modèle spécialement adapté pour la limite des arbres en Europe du Nord et dans les Alpes. Pour cette dernière, un module servant à simuler des changements d'utilisation du sol a également été ajouté. Tout d'abord, les paramètres de la courbe de réponse à la température pour les espèces inclues dans le modèle ont été calibrées au moyen de données dendrochronologiques pour diverses espèces et divers sites des deux écotones. Ceci a permis d'améliorer la courbe de croissance du modèle, mais a également permis de conclure que la régénération est probablement plus déterminante que la croissance en ce qui concerne la position de la limite des arbres et la distribution des espèces. La seconde étape consistait à implémenter le module d'abandon du terrain agricole dans les Alpes, basé sur un modèle statistique spatial existant. La sensibilité des variables les plus importantes du modèle a été testée et la performance de ce dernier comparée à d'autres approches de modélisation. La probabilité qu'un terrain soit abandonné était fortement influencée par la distance à la forêt la plus proche et par la pente, qui sont tous deux des substituts pour les coûts liés à la mise en culture. Lors de l'application en situation réelle, le nouveau modèle, baptisé TreeMig-LAb, a donné les résultats les plus réalistes. Ceux-ci étaient comparables aux conséquences déjà observées de l'abandon de terrains agricoles, telles que l'expansion des forêts existantes et la fermeture des clairières. Ce nouveau modèle a ensuite été mis en application dans deux zones d'étude, l'une dans les Alpes suisses et l'autre en Laponie finlandaise, avec divers scénarios de changement climatique. Ces derniers étaient basés sur les prévisions de changement de température pour le siècle prochain établies par l'IPCC et le modèle climatique HadCM3 (ΔT: +1.3, +3.5 et +5.6 °C) et comprenaient une période de stabilisation post-changement climatique de 300 ans. Les résultats ont montré des perturbations majeures dans les deux types de limites de arbres. Avec le scénario de changement climatique le moins extrême, les distributions respectives des espèces ont subi un simple glissement, mais il a fallu plusieurs siècles pour qu'elles atteignent un nouvel équilibre. Avec les autres scénarios, certaines espèces ont disparu de la zone d'étude (p. ex. Pinus cembra dans les Alpes) ou ont vu leur population diminuer parce qu'il n'y avait plus assez de terrains disponibles dans lesquels elles puissent migrer. Le résultat le plus frappant a été le temps de latence dans la réponse de la plupart des espèces, indépendamment du scénario de changement climatique utilisé ou du type de limite des arbres. Finalement, un modèle statistique de l'effet de l'abroutissement par les rennes (Rangifer tarandus) sur la croissance de Pinus sylvestris a été développé, comme première étape en vue de l'implémentation des impacts humains sur la limite boréale des arbres. L'effet attendu était indirect, puisque les rennes réduisent la couverture de lichen sur le sol, dont on attend un effet protecteur contre les rigueurs climatiques. Le modèle a mis en évidence un effet modeste mais significatif, mais étant donné que le lien avec les variables climatiques sous jacentes était peu clair et que le modèle n'était pas appliqué dans l'espace, il n'était pas utilisable tel quel. Le développement du modèle TreeMig-LAb a permis : a) d'établir une méthode pour déduire les paramètres spécifiques de l'équation de croissance ä partir de données dendrochronologiques, b) de mettre en évidence l'importance de la régénération dans la position de la limite des arbres et la distribution des espèces et c) d'améliorer l'intégration des sciences sociales dans les modèles de paysage. L'application du modèle aux limites alpines et nord-européennes des arbres sous différents scénarios de changement climatique a montré qu'avec la plupart des niveaux d'augmentation de température prévus, la limite des arbres subirait des perturbations majeures, avec des glissements d'aires de répartition et l'extinction potentielle de certaines espèces. Cependant, ces réponses ont montré des temps de latence importants, si bien que ces effets ne seraient pas visibles avant des décennies et pourraient mettre plusieurs siècles à se stabiliser.

Impacts of climate change on Swiss biodiversity: An indicator taxa approach

We present a new indicator taxa approach to the prediction of climate change effects on biodiversity at the national level in Switzerland. As indicators, we select a set of the most widely distributed species that account for 95% of geographical variation in sampled species richness of birds, butterflies, and vascular plants. Species data come from a national program designed to monitor spatial and temporal trends in species richness. We examine some opportunities and limitations in using these data. We develop ecological niche models for the species as functions of both climate and land cover variables. We project these models to the future using climate predictions that correspond to two IPCC 3rd assessment scenarios for the development of 'greenhouse' gas emissions. We find that models that are calibrated with Swiss national monitoring data perform well in 10-fold cross-validation, but can fail to capture the hot-dry end of environmental gradients that constrain some species distributions. Models for indicator species in all three higher taxa predict that climate change will result in turnover in species composition even where there is little net change in predicted species richness. Indicator species from high elevations lose most areas of suitable climate even under the relatively mild B2 scenario. We project some areas to increase in the number of species for which climate conditions are suitable early in the current century, but these areas become less suitable for a majority of species by the end of the century. Selection of indicator species based on rank prevalence results in a set of models that predict observed species richness better than a similar set of species selected based on high rank of model AUC values. An indicator species approach based on selected species that are relatively common may facilitate the use of national monitoring data for predicting climate change effects on the distribution of biodiversity.

European climate variations over the past half-millennium reconstructed from groundwater

Temperature reconstructions for recent centuries are the basis of estimations of the natural variability in the climate system before and during the onset of anthropogenic perturbation. Here we present, for the first time, an independent and physically based reconstruction of mean annual temperature over the past half millennium obtained from groundwater in France. The reconstructed noble gas temperature (NGT) record suggests cooler than present climate conditions throughout the 16th-19th centuries. Periods of warming occur in the 17th-18th and 20th century, while cooling is reconstructed in the 19th century. A noticeable coincidence with other temperature records is demonstrated. Deuterium excess varies in parallel with the NGT, and indicates variation in the seasonality of the aquifer recharge; whereas high excess air in groundwater indicates periods with high oscillations of the water table.

Stratigraphy and palaeoenvironment of Early Jurassic dinosaur-bearing volcaniclastic deposits from Toundoute (Province of Ouarzazate, Haut-Atlas, Morocco)

The Early Jurassic dinosaur site of Toundoute which yielded the basal sauropod Tazoudasaurus naimi is examined in the light of its stratigraphic, sedimentological and palaeoenvironmental context. A thin succession of Early Liassic marine carbonates (probably Hettangian-Sinemurian in age) is continuously overlain by continental beds with dinosaurs. These latter are assumed to be of Middle to Late Liassic age. The continental deposits include a large part of volcanoclastics, different from the Triassic basalts. The Jurassic volcanoclastics originated from an unknown but obviously close eruption centre. The continental sediments (channels and flood plain) were deposited under tropical climate conditions with alternating humid and dry episodes. The bones occur as isolated or partly articulated elements (parts of carcasses). The two bone-beds are related to typical mud-flows. This type of transport protected the bones from erosion, and favoured their burying and fossilization.

Bridging the Faraoni and Selli oceanic anoxic events: late Hauterivian to early Aptian dysaerobic to anaerobic phases in the Tethys

A detailed geochemical analysis was performed on the upper part of the Maiolica Formation in the Breggia (southern Switzerland) and Capriolo sections (northern Italy). The analysed sediments consist of well-bedded, partly siliceous, pelagic carbonate, which lodges numerous thin, dark and organic-rich layers. Stable-isotope, phosphorus, organic-carbon and a suite of redox-sensitive trace-element contents (RSTE: Mo, U, Co, V and As) were measured. The RSTE pattern and C-org:P-tot ratios indicate that most organic-rich layers were deposited under dysaerobic rather than anaerobic conditions and that latter conditions were likely restricted to short intervals in the latest Hauterivian, the early Barremian and the pre-Selli early Aptian. Correlations are both possible with organic-rich intervals in central Italy (the Gorgo a Cerbara section) and the Boreal Lower Saxony Basin, as well as with the facies and drowning pattern in the Helvetic segment of the northern Tethyan carbonate platform. Our data and correlations suggest that the latest Hauterivian witnessed the progressive installation of dysaerobic conditions in the Tethys, which went along with the onset in sediment condensation, phosphogenesis and platform drowning on the northern Tethyan margin, and which culminated in the Faraoni anoxic episode. This episode is followed by further episodes of dysaerobic conditions in the Tethys and the Lower Saxony Basin, which became more frequent and progressively stronger in the late early Barremian. Platform drowning persisted and did not halt before the latest early Barremian. The late Barremian witnessed diminishing frequencies and intensities in dysaerobic conditions, which went along with the progressive installation of the Urgonian carbonate platform. Near the Barremian-Aptian boundary, the increasing density in dysaerobic episodes in the Tethyan and Lower Saxony Basins is paralleled by a change towards heterozoan carbonate production on the northern Tethyan shelf. The following return to more oxygenated conditions is correlated with the second phase of Urgonian platform growth and the period immediately preceding and corresponding to the Selli anoxic episode is characterised by renewed platform drowning and the change to heterozoan carbonate production. Changes towards more humid climate conditions were the likely cause for the repetitive installation of dys- to anaerobic conditions in the Tethyan and Boreal basins and the accompanying changes in the evolution of the carbonate platform towards heterozoan carbonate-producing ecosystems and platform drowning.

Middle and late Cenomanian oceanic anoxic events in shallow and deeper shelf environments of western Morocco

The response of shallow-water sequences to oceanic anoxic event 2 and mid-Cenomanian events 1a and 1b was investigated along the west African margin of Morocco north of Agadir (Azazoul) and correlated with the deep-water sequence of the Tarfaya Basin (Mohammed Beach) based on biostratigraphy, mineralogy, phosphorus and stable isotopes. In the deeper Mohammed Beach section results show double peaks in delta 13C(org) for mid-Cenomanian events 1a and 1b (Rotalipora reicheli biozone, lower CC10a biozone), the characteristic oceanic anoxic event 2 delta 13C excursion (Rotalipora cushmani extinction, top of CC10a biozone) and laminated (anoxic) black shale. In the shallow environment north of Agadir, a fluctuating sea-level associated with dysoxic, brackish and mesotrophic conditions prevailed during the middle to late Cenomanian, as indicated by oyster biostromes, nannofossils, planktonic and benthonic foraminiferal assemblages. Anoxic conditions characteristic of oceanic anoxic event 2 (for example, laminated black shales) did not reach into shallow-water environments until the maximum transgression of the early Turonian. Climate conditions decoupled along the western margin of Morocco between mid-Cenomanian event 1b and the Cenomanian-Turonian boundary, as also observed in eastern Tethys. North of Agadir alternating humid and dry seasonal conditions prevailed, whereas in the Tarfaya Basin the climate was dry and seasonal. This climatic decoupling can be attributed to variations in the Intertropical Convergence Zone and in the intensity of the north-east trade winds in tropical areas.

Environment and dispersal paths override life strategies and residence time in determining regional patterns of invasion by alien plants.

We describe a novel dissimilarity framework to analyze spatial patterns of species diversity and illustrate it with alien plant invasions in Northern Portugal. We used this framework to test the hypothesis that patterns of alien invasive plant species richness and composition are differently affected by differences in climate, land use and landscape connectivity (i.e. Geographic distance as a proxy and vectorial objects that facilitate dispersal such as roads and rivers) between pairs of localities at the regional scale. We further evaluated possible effects of plant life strategies (Grime's C-S-R) and residence time. Each locality consisted of a 1 km(2) landscape mosaic in which all alien invasive species were recorded by visiting all habitat types. Multi-model inference revealed that dissimilarity in species richness is more influenced by environmental distance (particularly climate), whereas geographic distance (proxies for dispersal limitations) is more important to explain dissimilarity in species composition, with a prevailing role for ecotones and roads. However, only minor differences were found in the responses of the three C-S-R strategies. Some effect of residence time was found, but only for dissimilarity in species richness. Our results also indicated that environmental conditions (e.g. climate conditions) limit the number of alien species invading a given site, but that the presence of dispersal corridors determines the paths of invasion and therefore the pool of species reaching each site. As geographic distances (e.g. ecotones and roads) tend to explain invasion at our regional scale highlights the need to consider the management of alien invasions in the context of integrated landscape planning. Alien species management should include (but not be limited to) the mitigation of dispersal pathways along linear infrastructures. Our results therefore highlight potentially useful applications of the novel multimodel framework to the anticipation and management of plant invasions. (C) 2013 Elsevier GmbH. All rights reserved.

Cenomanian-Turonian transition in a shallow water sequence of the Sinai, Egypt

Environmental and depositional changes across the Late Cenomanian oceanic anoxic event (OAE2) in the Sinai, Egypt, are examined based on biostratigraphy, mineralogy, delta(13)C values and phosphorus analyses. Comparison with the Pueblo, Colorado, stratotype section reveals the Whadi El Ghaib section as stratigraphically complete across the late Cenomanian-early Turonian. Foraminifera are dominated by high-stress planktic and benthic assemblages characterized by low diversity, low-oxygen and low-salinity tolerant species, which mark shallow-water oceanic dysoxic conditions during OAE2. Oyster biostromes suggest deposition occurred in less than 50 m depths in low-oxygen, brackish, and nutrient-rich waters. Their demise prior to the peak delta(13)C excursion is likely due to a rising sea-level. Characteristic OAE2 anoxic conditions reached this coastal region only at the end of the delta(13)C plateau in deeper waters near the end of the Cenomanian. Increased phosphorus accumulations before and after the delta(13)C excursion suggest higher oxic conditions and increased detrital input. Bulk-rock and clay mineralogy indicate humid climate conditions, increased continental runoff and a rising sea up to the first delta(13)C peak. Above this interval, a dryer and seasonally well-contrasted climate with intermittently dry conditions prevailed. These results reveal the globally synchronous delta(13)C shift, but delayed effects of OAE2 dependent on water depth.

What is the potential of spread in invasive bryophytes?

Although the number of invasive bryophytes is much lower than that of higher plants, they threaten habitats that are often species rich and of high conservation relevance. Their potential of spread has, however, never been determined. Here, we assess whether the three most invasive bryophyte species shifted their niche during the invasion process and whether the extent of the study area defined to calibrate the model (geographic background, GB) affects model transferability. We then determine whether ecological niche models (ENMs) developed in their native range can be projected in other areas to assess their invasive potential. The macroclimatic niches of Campylopus introflexus, Orthodontium lineare and Lophocolea semiteres were compared in their native range (Southern Hemisphere) and in their invasion range (Northern Hemisphere) using ordination techniques. ENMs from an ensemble model were calibrated in the native range and projected onto the Northern Hemisphere using different GBs. No evidence for niche expansion in the invaded range was found and the species occur in the invaded range under climate conditions that are similar to those in the native range. The performance of the models to predict occurrences in the invaded range increased with the extent of the GB. The potential range of all species included entire regions on continents where they are still absent. The expansion of the investigated species appears to be constrained by climate conditions that are similar to those currently prevailing in their native range, which is consistent with our failure to demonstrate macroclimatic niche shift in the invaded range. The use of large GBs is recommended in such vagile organisms with large, disjunct distributions. The models indicated that invasive bryophyte species might become a threat in central and eastern Europe, North America and eastern Asia if accidentally introduced or naturally dispersed.

Areas of high conservation value in Georgia: present and future threats by invasive alien plants

Georgia is known for its extraordinary rich biodiversity of plants, which may now be threatened due to the spread of invasive alien plants (IAP). We aimed to identify (i) the most prominent IAP out of 9 selected potentially invasive and harmful IAP IAP by predicting thetheir distribution of 9 selected IAP under current and future climate conditions in Georgia as well as in its 43 Protected Areas, as a proxy for areas of high conservation value and (ii) the Protected Areas most at risk due to these IAP. We used species distribution models based on 6 climate variables and then filtered the obtained distributions based on maps of soil and vegetation types, and on recorded occurrences, resulting into the predicted ecological distribution of the 9 IAP's at a resolution of 1km2. We foundOur habitat suitability analysis showed that Ambrosia artemisiifolia, (24% and 40%) Robinia pseudoacaia (14% and 19%) and Ailanthus altissima (9% and 11%) have the largest potential distribution are the most abundant (predicted % area covered)d) IAP, with Ailanthus altissima the potentially most increasing one over the next fifty years (from 9% to 13% and from 11% to 25%), for Georgia and the Protected Areas, respectively. Furthermore, our results show indicate two areas in Georgia that are under specifically high threat, i.e. the area around Tbilisi and an area in the western part of Georgia (Adjara), both at lower altitudes. Our procedure to identify areas of high conservation value most at risk by IAP has been applied for the first time. It will help national authorities in prioritizing their measures to protect Georgia's outstanding biodiversity from the negative impact of IAP.

Influence of climate on the presence of colour polymorphism in two montane reptile species.

The coloration of ectotherms plays an important role in thermoregulation processes. Dark individuals should heat up faster and be able to reach a higher body temperature than light individuals and should therefore have benefits in cool areas. In central Europe, montane local populations of adder (Vipera berus) and asp viper (Vipera aspis) exhibit a varying proportion of melanistic individuals. We tested whether the presence of melanistic V. aspis and V. berus could be explained by climatic conditions. We measured the climatic niche position and breadth of monomorphic (including strictly patterned individuals) and polymorphic local populations, calculated their niche overlap and tested for niche equivalency and similarity. In accordance with expectations, niche overlap between polymorphic local populations of both species is high, and even higher than that of polymorphic versus monomorphic montane local populations of V. aspis, suggesting a predominant role of melanism in determining the niche of ectothermic vertebrates. However, unexpectedly, the niche of polymorphic local populations of both species is narrower than that of monomorphic ones, indicating that colour polymorphism does not always enable the exploitation of a greater variability of resources, at least at the intraspecific level. Overall, our results suggest that melanism might be present only when the thermoregulatory benefit is higher than the cost of predation.

Climate-based empirical models show biased predictions of butterfly communities along environmental gradients

A better understanding of the factors that mould ecological community structure is required to accurately predict community composition and to anticipate threats to ecosystems due to global changes. We tested how well stacked climate-based species distribution models (S-SDMs) could predict butterfly communities in a mountain region. It has been suggested that climate is the main force driving butterfly distribution and community structure in mountain environments, and that, as a consequence, climate-based S-SDMs should yield unbiased predictions. In contrast to this expectation, at lower altitudes, climate-based S-SDMs overpredicted butterfly species richness at sites with low plant species richness and underpredicted species richness at sites with high plant species richness. According to two indices of composition accuracy, the Sorensen index and a matching coefficient considering both absences and presences, S-SDMs were more accurate in plant-rich grasslands. Butterflies display strong and often specialised trophic interactions with plants. At lower altitudes, where land use is more intense, considering climate alone without accounting for land use influences on grassland plant richness leads to erroneous predictions of butterfly presences and absences. In contrast, at higher altitudes, where climate is the main force filtering communities, there were fewer differences between observed and predicted butterfly richness. At high altitudes, even if stochastic processes decrease the accuracy of predictions of presence, climate-based S-SDMs are able to better filter out butterfly species that are unable to cope with severe climatic conditions, providing more accurate predictions of absences. Our results suggest that predictions should account for plants in disturbed habitats at lower altitudes but that stochastic processes and heterogeneity at high altitudes may limit prediction success of climate-based S-SDMs.

Will climate change drive alien invasive plants into areas of high protection value? An improved model-based regional assessment to prioritise the management of invasions.

Species distribution models (SDMs) studies suggest that, without control measures, the distribution of many alien invasive plant species (AIS) will increase under climate and land-use changes. Due to limited resources and large areas colonised by invaders, management and monitoring resources must be prioritised. Choices depend on the conservation value of the invaded areas and can be guided by SDM predictions. Here, we use a hierarchical SDM framework, complemented by connectivity analysis of AIS distributions, to evaluate current and future conflicts between AIS and high conservation value areas. We illustrate the framework with three Australian wattle (Acacia) species and patterns of conservation value in Northern Portugal. Results show that protected areas will likely suffer higher pressure from all three Acacia species under future climatic conditions. Due to this higher predicted conflict in protected areas, management might be prioritised for Acacia dealbata and Acacia melanoxylon. Connectivity of AIS suitable areas inside protected areas is currently lower than across the full study area, but this would change under future environmental conditions. Coupled SDM and connectivity analysis can support resource prioritisation for anticipation and monitoring of AIS impacts. However, further tests of this framework over a wide range of regions and organisms are still required before wide application.