Predicted decrease in global climate suitability masks regional complexity of invasive fruit fly species response to climate change

Climate change affects the rate of insect invasions as well as the abundance, distribution and impacts of such invasions on a global scale. Among the principal analytical approaches to predicting and understanding future impacts of biological invasions are Species Distribution Models (SDMs), typically in the form of correlative Ecological Niche Models (ENMs). An underlying assumption of ENMs is that species-environment relationships remain preserved during extrapolations in space and time, although this is widely criticised. The semi-mechanistic modelling platform, CLIMEX, employs a top-down approach using species ecophysiological traits and is able to avoid some of the issues of extrapolation, making it highly applicable to investigating biological invasions in the context of climate change. The tephritid fruit flies (Diptera: Tephritidae) comprise some of the most successful invasive species and serious economic pests around the world. Here we project 12 tephritid species CLIMEX models into future climate scenarios to examine overall patterns of climate suitability and forecast potential distributional changes for this group. We further compare the aggregate response of the group against species-specific responses. We then consider additional drivers of biological invasions to examine how invasion potential is influenced by climate, fruit production and trade indices. Considering the group of tephritid species examined here, climate change is predicted to decrease global climate suitability and to shift the cumulative distribution poleward. However, when examining species-level patterns, the predominant directionality of range shifts for 11 of the 12 species is eastward. Most notably, management will need to consider regional changes in fruit fly species invasion potential where high fruit production, trade indices and predicted distributions of these flies overlap.

Local and global error models to improve uncertainty quantification

In groundwater applications, Monte Carlo methods are employed to model the uncertainty on geological parameters. However, their brute-force application becomes computationally prohibitive for highly detailed geological descriptions, complex physical processes, and a large number of realizations. The Distance Kernel Method (DKM) overcomes this issue by clustering the realizations in a multidimensional space based on the flow responses obtained by means of an approximate (computationally cheaper) model; then, the uncertainty is estimated from the exact responses that are computed only for one representative realization per cluster (the medoid). Usually, DKM is employed to decrease the size of the sample of realizations that are considered to estimate the uncertainty. We propose to use the information from the approximate responses for uncertainty quantification. The subset of exact solutions provided by DKM is then employed to construct an error model and correct the potential bias of the approximate model. Two error models are devised that both employ the difference between approximate and exact medoid solutions, but differ in the way medoid errors are interpolated to correct the whole set of realizations. The Local Error Model rests upon the clustering defined by DKM and can be seen as a natural way to account for intra-cluster variability; the Global Error Model employs a linear interpolation of all medoid errors regardless of the cluster to which the single realization belongs. These error models are evaluated for an idealized pollution problem in which the uncertainty of the breakthrough curve needs to be estimated. For this numerical test case, we demonstrate that the error models improve the uncertainty quantification provided by the DKM algorithm and are effective in correcting the bias of the estimate computed solely from the MsFV results. The framework presented here is not specific to the methods considered and can be applied to other combinations of approximate models and techniques to select a subset of realizations

Adapting global conservation strategies to climate change at the European scale: the otter as a flagship species

Climate change has created the need for new strategies in conservation planning that account for the dynamics of factors threatening endangered species. Here we assessed climate change threat to the European otter, a flagship species for freshwater ecosystems, considering how current conservation areas will perform in preserving the species in a climatically changed future. We used an ensemble forecasting approach considering six modelling techniques applied to eleven subsets of otter occurrences across Europe. We performed a pseudo-independent and an internal evaluation of predictions. Future projections of species distribution were made considering the A2 and B2 scenarios for 2080 across three climate models: CCCMA-CGCM2, CSIRO-MK2 and HCCPR HAD-CM3. The current and the predicted otter distributions were used to identify priority areas for the conservation of the species, and overlapped to existing network of protected areas. Our projections show that climate change may profoundly reshuffle the otter's potential distribution in Europe, with important differences between the two scenarios we considered. Overall, the priority areas for conservation of the otter in Europe appear to be unevenly covered by the existing network of protected areas, with the current conservation efforts being insufficient in most cases. For a better conservation, the existing protected areas should be integrated within a more general conservation and management strategy incorporating climate change projections. Due to the important role that the otter plays for freshwater habitats, our study further highlights the potential sensitivity of freshwater habitats in Europe to climate change.

How robust are global conservation priorities to climate change?

International conservation organisations have identified priority areas for biodiversity conservation. These global-scale prioritisations affect the distribution of funds for conservation interventions. As each organisation has a different focus, each prioritisation scheme is determined by different decision criteria and the resultant priority areas vary considerably. However, little is known about how the priority areas will respond to the impacts of climate change. In this paper, we examined the robustness of eight global-scale prioritisations to climate change under various climate predictions from seven global circulation models. We developed a novel metric of the climate stability for 803 ecoregions based on a recently introduced method to estimate the overlap of climate envelopes. The relationships between the decision criteria and the robustness of the global prioritisation schemes were statistically examined. We found that decision criteria related to level of endemism and landscape fragmentation were strongly correlated with areas predicted to be robust to a changing climate. Hence, policies that prioritise intact areas due to the likely cost efficiency, and assumptions related to the potential to mitigate the impacts of climate change, require further examination. Our findings will help determine where additional management is required to enable biodiversity to adapt to the impacts of climate change

Warming experiments underpredict plant phenological responses to climate change.

Warming experiments are increasingly relied on to estimate plant responses to global climate change. For experiments to provide meaningful predictions of future responses, they should reflect the empirical record of responses to temperature variability and recent warming, including advances in the timing of flowering and leafing. We compared phenology (the timing of recurring life history events) in observational studies and warming experiments spanning four continents and 1,634 plant species using a common measure of temperature sensitivity (change in days per degree Celsius). We show that warming experiments underpredict advances in the timing of flowering and leafing by 8.5-fold and 4.0-fold, respectively, compared with long-term observations. For species that were common to both study types, the experimental results did not match the observational data in sign or magnitude. The observational data also showed that species that flower earliest in the spring have the highest temperature sensitivities, but this trend was not reflected in the experimental data. These significant mismatches seem to be unrelated to the study length or to the degree of manipulated warming in experiments. The discrepancy between experiments and observations, however, could arise from complex interactions among multiple drivers in the observational data, or it could arise from remediable artefacts in the experiments that result in lower irradiance and drier soils, thus dampening the phenological responses to manipulated warming. Our results introduce uncertainty into ecosystem models that are informed solely by experiments and suggest that responses to climate change that are predicted using such models should be re-evaluated.

Forecasting changes in population genetic structure of alpine plants in response to global warming.

Species range shifts in response to climate and land use change are commonly forecasted with species distribution models based on species occurrence or abundance data. Although appealing, these models ignore the genetic structure of species, and the fact that different populations might respond in different ways because of adaptation to their environment. Here, we introduced ancestry distribution models, that is, statistical models of the spatial distribution of ancestry proportions, for forecasting intra-specific changes based on genetic admixture instead of species occurrence data. Using multi-locus genotypes and extensive geographic coverage of distribution data across the European Alps, we applied this approach to 20 alpine plant species considering a global increase in temperature from 0.25 to 4 °C. We forecasted the magnitudes of displacement of contact zones between plant populations potentially adapted to warmer environments and other populations. While a global trend of movement in a north-east direction was predicted, the magnitude of displacement was species-specific. For a temperature increase of 2 °C, contact zones were predicted to move by 92 km on average (minimum of 5 km, maximum of 212 km) and by 188 km for an increase of 4 °C (minimum of 11 km, maximum of 393 km). Intra-specific turnover-measuring the extent of change in global population genetic structure-was generally found to be moderate for 2 °C of temperature warming. For 4 °C of warming, however, the models indicated substantial intra-specific turnover for ten species. These results illustrate that, in spite of unavoidable simplifications, ancestry distribution models open new perspectives to forecast population genetic changes within species and complement more traditional distribution-based approaches.

Predicting community and ecosystem outcomes of mycorrhizal responses to global change.

Mycorrhizal symbioses link the biosphere with the lithosphere by mediating nutrient cycles and energy flow though terrestrial ecosystems. A more mechanistic understanding of these plant-fungal associations may help ameliorate anthropogenic changes to C and N cycles and biotic communities. We explore three interacting principles: (1) optimal allocation, (2) biotic context and (3) fungal adaptability that may help predict mycorrhizal responses to carbon dioxide enrichment, nitrogen eutrophication, invasive species and land-use changes. Plant-microbial feedbacks and thresholds are discussed in light of these principles with the goal of generating testable hypotheses. Ideas to develop large-scale collaborative research efforts are presented. It is our hope that mycorrhizal symbioses can be effectively integrated into global change models and eventually their ecology will be understood well enough so that they can be managed to help offset some of the detrimental effects of anthropogenic environmental change.

Integrating species distribution models (SDMs) and phylogeography for two species of Alpine Primula.

The major intention of the present study was to investigate whether an approach combining the use of niche-based palaeodistribution modeling and phylo-geography would support or modify hypotheses about the Quaternary distributional history derived from phylogeographic methods alone. Our study system comprised two closely related species of Alpine Primula. We used species distribution models based on the extant distribution of the species and last glacial maximum (LGM) climate models to predict the distribution of the two species during the LGM. Phylogeographic data were generated using amplified fragment length polymorphisms (AFLPs). In Primula hirsuta, models of past distribution and phylogeographic data are partly congruent and support the hypothesis of widespread nunatak survival in the Central Alps. Species distribution models (SDMs) allowed us to differentiate between alpine regions that harbor potential nunatak areas and regions that have been colonized from other areas. SDMs revealed that diversity is a good indicator for nunataks, while rarity is a good indicator for peripheral relict populations that were not source for the recolonization of the inner Alps. In P. daonensis, palaeo-distribution models and phylogeographic data are incongruent. Besides the uncertainty inherent to this type of modeling approach (e.g., relatively coarse 1-km grain size), disagreement of models and data may partly be caused by shifts of ecological niche in both species. Nevertheless, we demonstrate that the combination of palaeo-distribution modeling with phylogeographical approaches provides a more differentiated picture of the distributional history of species and partly supports (P. hirsuta) and partly modifies (P. daonensis and P. hirsuta) hypotheses of Quaternary distributional history. Some of the refugial area indicated by palaeodistribution models could not have been identified with phylogeographic data.

Middle and Late Jurassic carbon stable-isotope stratigraphy and radiolarite sedimentation of the Umbria-Marche basin (Central Italy)

A continuous carbon isotope curve from Middle-Upper Jurassic pelagic carbonate rocks was acquired from two sections in the southern part of the Umbria-Marche Apennines in central Italy. At the Colle Bertone section (Terni) and the Terminilletto section (Rieti), the Upper Toarcian to Bajocian Calcari e Marne a Posidonia Formation and the Aalenian to Kimmeridgian Calcari e Marne a Posidonia and Calcari Diasprigni formations were sampled, respectively. Biostratigraphy in both sections is based on rich assemblages of calcareous nannofossils and radiolarians, as well as some ammonites found in the upper Toarcian-Bajocian interval. Both sections revealed a relative minimum of delta(13)C(PDB) close to + 2 parts per thousand in the Aalenian and a maximum around 3.5 parts per thousand in early Bajocian, associated with an increase in visible chert. In basinal sections in Umbria-Marche, this interval includes the very cherry base of the Calcari Diasprigni Formation (e.g. at Valdorbia) or the chert-rich uppermost portion of the Calcari a Posidonia (e.g at Bosso). In the Terminilletto section, the Bajocian-early Barthonian interval shows a gradual decrease in delta(13)C(PDB) values and a low around 2.3 parts per thousand. This part of the section is characterised by more than 40 m of almost chart-free limestones and correlates with a recurrence of limestone-rich facies in basinal sections at Valdorbia. A double peak with values of delta(13)C(PDB) around + 3 parts per thousand was observed in the Callovian and Oxfordian, constrained by well preserved radiolarian faunas. The maxima lie in the Callovian and the middle Oxfordian, and the minimum between the two peaks should be near the Callovian/Oxfordian boundary. In the Terminilletto section, visible chert increases together with delta(13)C(PDB) values from the middle Bathonian and reaches peak values in the Callovian-Oxfordian. In basinal sections in Umbria-Marche, a sharp increase in visible chert is observed at this level within the Calcari Diasprigni. A drop of delta(13)C values towards + 2 parts per thousand occurs in the Kimmeridgian and coincides with a decrease of visible chert in outcrop. The observed delta(13)C positive anomalies during the early Bajocian and the Callovian-Oxfordian may record changes in global climate towards warmer, more humid periods characterised by increased nutrient mobilisation and increased carbon burial. High biosiliceous (radiolarians, siliceous sponges) productivity and preservation appear to coincide with the delta(13)C positive anomalies, when the production of platform carbonates was subdued and ceased in many areas, with a drastic reduction of periplatform ooze input in many Tethyan basins. The carbon and silica cycles appear to be linked through global warming and increased continental weathering. Hydrothermal events related to extensive rifting and/or accelerated oceanic spreading may be the endogenic driving force that created a perturbation of the exogenic system (excess CO2 into the atmosphere and greenhouse conditions) reflected by the positive delta(13)C shifts and biosiliceous episodes.

Combining palaeodistribution modelling and phylogeographical approaches for identifying glacial refugia in Alpine Primula

Aim We investigated the late Quaternary history of two closely related and partly sympatric species of Primula from the south-western European Alps, P. latifolia Lapeyr. and P. marginata Curtis, by combining phylogeographical and palaeodistribution modelling approaches. In particular, we were interested in whether the two approaches were congruent and identified the same glacial refugia. Location South-western European Alps. Methods For the phylogeographical analysis we included 353 individuals from 28 populations of P. marginata and 172 individuals from 15 populations of P. latifolia and used amplified fragment length polymorphisms (AFLPs). For palaeodistribution modelling, species distribution models (SDMs) were based on extant species occurrences and then projected to climate models (CCSM, MIROC) of the Last Glacial Maximum (LGM), approximately 21 ka. Results The locations of the modelled LGM refugia were confirmed by various indices of genetic variation. The refugia of the two species were largely geographically isolated, overlapping only 6% to 11% of the species' total LGM distribution. This overlap decreased when the position of the glacial ice sheet and the differential elevational and edaphic distributions of the two species were considered. Main conclusions The combination of phylogeography and palaeodistribution modelling proved useful in locating putative glacial refugia of two alpine species of Primula. The phylogeographical data allowed us to identify those parts of the modelled LGM refugial area that were likely source areas for recolonization. The use of SDMs predicted LGM refugial areas substantially larger and geographically more divergent than could have been predicted by phylogeographical data alone

7

Climate change data and predictions for the Himalayas are very sparse and uncertain, characterized by a ?Himalayan data gap? and difficulties in predicting changes due to topographic complexity. A few reliable studies and climate change models for Nepal predict considerable changes: shorter monsoon seasons, more intensive rainfall patterns, higher temperatures, and drought. These predictions are confirmed by farmers who claim that temperatures have been increasing for the past decade and wonder why the rains have ?gone mad.? The number of hazard events, notably droughts, floods, and landslides are increasing and now account for approximately 100 deaths in Nepal annually. Other effects are drinking water shortages and shifting agricultural patterns, with many communities struggling to meet basic food security before climatic conditions started changing. The aim of this paper is to examine existing gaps between current climate models and the realities of local development planning through a case study on flood risk and drinking water management for the Municipality of Dharan in Eastern Nepal. This example highlights current challenges facing local-level governments, namely, flood and landslide mitigation, providing basic amenities ? especially an urgent lack of drinking water during the dry season ? poor local planning capacities, and limited resources. In this context, the challenge for Nepal will be to simultaneously address increasing risks caused by hazard events alongside the omnipresent food security and drinking water issues in both urban and rural areas. Local planning is needed that integrates rural development and disaster risk reduction (DRR) with knowledge about climate change considerations. The paper concludes with a critical analysis of climate change modeling and the gap between scientific data and low-tech and low capacities of local planners to access or implement adequate adaptation measures. Recommendations include the need to bridge gaps between scientific models, the local political reality and local information needs.

Solution scanning as a key policy tool: identifying management interventions to help maintain and enhance regulating ecosystem services

The major task of policy makers and practitioners when confronted with a resource management problem is to decide on the potential solution(s) to adopt from a range of available options. However, this process is unlikely to be successful and cost effective without access to an independently verified and comprehensive available list of options. There is currently burgeoning interest in ecosystem services and quantitative assessments of their importance and value. Recognition of the value of ecosystem services to human well-being represents an increasingly important argument for protecting and restoring the natural environment, alongside the moral and ethical justifications for conservation. As well as understanding the benefits of ecosystem services, it is also important to synthesize the practical interventions that are capable of maintaining and/or enhancing these services. Apart from pest regulation, pollination, and global climate regulation, this type of exercise has attracted relatively little attention. Through a systematic consultation exercise, we identify a candidate list of 296 possible interventions across the main regulating services of air quality regulation, climate regulation, water flow regulation, erosion regulation, water purification and waste treatment, disease regulation, pest regulation, pollination and natural hazard regulation. The range of interventions differs greatly between habitats and services depending upon the ease of manipulation and the level of research intensity. Some interventions have the potential to deliver benefits across a range of regulating services, especially those that reduce soil loss and maintain forest cover. Synthesis and applications: Solution scanning is important for questioning existing knowledge and identifying the range of options available to researchers and practitioners, as well as serving as the necessary basis for assessing cost effectiveness and guiding implementation strategies. We recommend that it become a routine part of decision making in all environmental policy areas.

Biological evidence supports an early and complex emergence of the Isthmus of Panama.

The linking of North and South America by the Isthmus of Panama had major impacts on global climate, oceanic and atmospheric currents, and biodiversity, yet the timing of this critical event remains contentious. The Isthmus is traditionally understood to have fully closed by ca. 3.5 million years ago (Ma), and this date has been used as a benchmark for oceanographic, climatic, and evolutionary research, but recent evidence suggests a more complex geological formation. Here, we analyze both molecular and fossil data to evaluate the tempo of biotic exchange across the Americas in light of geological evidence. We demonstrate significant waves of dispersal of terrestrial organisms at approximately ca. 20 and 6 Ma and corresponding events separating marine organisms in the Atlantic and Pacific oceans at ca. 23 and 7 Ma. The direction of dispersal and their rates were symmetrical until the last ca. 6 Ma, when northern migration of South American lineages increased significantly. Variability among taxa in their timing of dispersal or vicariance across the Isthmus is not explained by the ecological factors tested in these analyses, including biome type, dispersal ability, and elevation preference. Migration was therefore not generally regulated by intrinsic traits but more likely reflects the presence of emergent terrain several millions of years earlier than commonly assumed. These results indicate that the dramatic biotic turnover associated with the Great American Biotic Interchange was a long and complex process that began as early as the Oligocene-Miocene transition.