Clinal variation in a brown lemur (Eulemur spp.) hybrid zone: Combining morphological, genetic and climatic data to examine stability.

Studies of hybrid zones can inform our understanding of reproductive isolation and speciation. Two species of brown lemur (Eulemur rufifrons and E. cinereiceps) form an apparently stable hybrid zone in the Andringitra region of south-eastern Madagascar. The aim of this study was to identify factors that contribute to this stability. We sampled animals at 11 sites along a 90-km transect through the hybrid zone and examined variation in 26 microsatellites, the D-loop region of mitochondrial DNA, six pelage and nine morphological traits; we also included samples collected in more distant allopatric sites. Clines in these traits were noncoincident, and there was no increase in either inbreeding coefficients or linkage disequilibrium at the centre of the zone. These results could suggest that the hybrid zone is maintained by weak selection against hybrids, conforming to either the tension zone or geographical selection-gradient model. However, a closer examination of clines in pelage and microsatellites indicates that these clines are not sigmoid or stepped in shape but instead plateau at their centre. Sites within the hybrid zone also occur in a distinct habitat, characterized by greater seasonality in precipitation and lower seasonality in temperature. Together, these findings suggest that the hybrid zone may follow the bounded superiority model, with exogenous selection favouring hybrids within the transitional zone. These findings are noteworthy, as examples supporting the bounded superiority model are rare and may indicate a process of ecologically driven speciation without geographical isolation.

Mapping of Environmental Data Using Kernel-Based Methods

Recently, kernel-based Machine Learning methods have gained great popularity in many data analysis and data mining fields: pattern recognition, biocomputing, speech and vision, engineering, remote sensing etc. The paper describes the use of kernel methods to approach the processing of large datasets from environmental monitoring networks. Several typical problems of the environmental sciences and their solutions provided by kernel-based methods are considered: classification of categorical data (soil type classification), mapping of environmental and pollution continuous information (pollution of soil by radionuclides), mapping with auxiliary information (climatic data from Aral Sea region). The promising developments, such as automatic emergency hot spot detection and monitoring network optimization are discussed as well.

Climate change, anthropogenic disturbance and the northward range expansion of Lactuca serriola (Asteraceae)

The distribution range of Lactuca serriola, a species native to the summer-dry mediterranean climate, has expanded northwards during the last 250 years. This paper assesses the influence of climate on the range expansion of this species and highlights the importance of anthropogenic disturbance to its spread. Location Central and Northern Europe. Methods Data on the geographic distribution of L. serriola were assembled through a literature search as well as through floristic and herbarium surveys. Maps of the spread of L. serriola in Central and Northern Europe were prepared based on herbarium data. The spread was assessed more precisely in Germany, Austria and Great Britain by pooling herbarium and literature data. We modelled the bioclimatic niche of the species using occurrence and climatic data covering the last century to generate projections of suitable habitats under the climatic conditions of five time periods. We tested whether the observed distribution of L. serriola could be explained for each time period, assuming that the climatic niche of the species was conserved across time. Results The species has spread northwards since the beginning of the 19th century. We show that climate warming in Europe increased the number of sites suitable for the species at northern latitudes. Until the late 1970s, the distribution of the species corresponded to the climatically suitable sites available. For the last two decades, however, we could not show any significant relationship between the increase in suitable sites and the distributional range change of L. serriola. However, we highlight potential areas the species could spread to in the future (Great Britain, southern Scandinavia and the Swedish coast). It is predominantly non-climatic influences of global change that have contributed to its rapid spread. Main conclusions The observation that colonizing species are not filling their climatically suitable range might imply that, potentially, other ruderal species could expand far beyond their current range. Our work highlights the importance of historical floristic and herbarium data for understanding the expansion of a species. Such historical distributional data can provide valuable information for those planning the management of contemporary environmental problems, such as species responses to environmental change.

Data-driven topo-climatic mapping with machine learning methods

Automatic environmental monitoring networks enforced by wireless communication technologies provide large and ever increasing volumes of data nowadays. The use of this information in natural hazard research is an important issue. Particularly useful for risk assessment and decision making are the spatial maps of hazard-related parameters produced from point observations and available auxiliary information. The purpose of this article is to present and explore the appropriate tools to process large amounts of available data and produce predictions at fine spatial scales. These are the algorithms of machine learning, which are aimed at non-parametric robust modelling of non-linear dependencies from empirical data. The computational efficiency of the data-driven methods allows producing the prediction maps in real time which makes them superior to physical models for the operational use in risk assessment and mitigation. Particularly, this situation encounters in spatial prediction of climatic variables (topo-climatic mapping). In complex topographies of the mountainous regions, the meteorological processes are highly influenced by the relief. The article shows how these relations, possibly regionalized and non-linear, can be modelled from data using the information from digital elevation models. The particular illustration of the developed methodology concerns the mapping of temperatures (including the situations of Föhn and temperature inversion) given the measurements taken from the Swiss meteorological monitoring network. The range of the methods used in the study includes data-driven feature selection, support vector algorithms and artificial neural networks.

Contrasting spatio-temporal climatic niche dynamics during the eastern and western invasions of spotted knapweed in North America

Aim The spotted knapweed (Centaurea stoebe), a plant native to south-east and central Europe, is highly invasive in North America. We investigated the spatio-temporal climatic niche dynamics of the spotted knapweed in North America along two putative eastern and western invasion routes. We then considered the patterns observed in the light of historical, ecological and evolutionary factors. Location Europe and North America. Methods The niche characteristics of the east and west invasive populations of spotted knapweed in North America were determined from documented occurrences over 120 consecutive years (1890-2010). The 2.5 and 97.5 percentiles of values along temperature and precipitation gradients, as given by the two first axes of a principal component axis (PCA), were then calculated. We additionally measured the climatic dissimilarity between invaded and native niches using a multivariate environmental similarity surface (MESS) analysis. Results Along both invasion routes, the species established in regions with climatic conditions that were similar to those in the native range in Europe. An initial spread in ruderal habitats always preceded spread in (semi-)natural habitats. In the east, the niche gradually increased over time until it reached limits similar to the native niche. Conversely, in the west the niche abruptly expanded after an extended time lag into climates not occupied in the native range; only the native cold niche limit was conserved. Main conclusions Our study reveals that different niche dynamics have taken place during the eastern and western invasions. This pattern indicates different combinations of historical, ecological and evolutionary factors in the two ranges. We hypothesize that the lack of a well-developed transportation network in the west at the time of the introduction of spotted knapweed confined the species to a geographically and climatically isolated region. The invasion of dry rangelands may have been favoured during the agricultural transition in the 1930s by release from natural enemies, local adaptation and less competitive vegetation, but further experimental and molecular studies are needed to explain these contrasting niche patterns fully. Our study illustrates the need and benefit of applying large-scale, temporally explicit approaches to understanding biological invasions.

Evidence of climatic niche shift during biological invasion

Niche-based models calibrated in the native range by relating species observations to climatic variables are commonly used to predict the potential spatial extent of species' invasion. This climate matching approach relies on the assumption that invasive species conserve their climatic niche in the invaded ranges. We test this assumption by analysing the climatic niche spaces of Spotted Knapweed in western North America and Europe. We show with robust cross-continental data that a shift of the observed climatic niche occurred between native and non-native ranges, providing the first empirical evidence that an invasive species can occupy climatically distinct niche spaces following its introduction into a new area. The models fail to predict the current invaded distribution, but correctly predict areas of introduction. Climate matching is thus a useful approach to identify areas at risk of introduction and establishment of newly or not-yet-introduced neophytes, but may not predict the full extent of invasions.

Phylogeography and Pleistocene refugia of the adder (Vipera berus) as inferred from mitochondrial DNA sequence data.

In order to contribute to the debate about southern glacial refugia used by temperate species and more northern refugia used by boreal or cold-temperate species, we examined the phylogeography of a widespread snake species (Vipera berus) inhabiting Europe up to the Arctic Circle. The analysis of the mitochondrial DNA (mtDNA) sequence variation in 1043 bp of the cytochrome b gene and in 918 bp of the noncoding control region was performed with phylogenetic approaches. Our results suggest that both the duplicated control region and cytochrome b evolve at a similar rate in this species. Phylogenetic analysis showed that V. berus is divided into three major mitochondrial lineages, probably resulting from an Italian, a Balkan and a Northern (from France to Russia) refugial area in Eastern Europe, near the Carpathian Mountains. In addition, the Northern clade presents an important substructure, suggesting two sequential colonization events in Europe. First, the continent was colonized from the three main refugial areas mentioned above during the Lower-Mid Pleistocene. Second, recolonization of most of Europe most likely originated from several refugia located outside of the Mediterranean peninsulas (Carpathian region, east of the Carpathians, France and possibly Hungary) during the Mid-Late Pleistocene, while populations within the Italian and Balkan Peninsulas fluctuated only slightly in distribution range, with larger lowland populations during glacial times and with refugial mountain populations during interglacials, as in the present time. The phylogeographical structure revealed in our study suggests complex recolonization dynamics of the European continent by V. berus, characterized by latitudinal as well as altitudinal range shifts, driven by both climatic changes and competition with related species.

Phosphorus and carbon burial during cretaceous oceanic anoxic events, links with climatic changes

According to Jenkyns (2010), oceanic anoxic events (OAE) record profound changes in the climatic and paleoceanographic state of the planet and represent major disturbances in the global carbon cycle. One of the most studied OAEs on a worldwide scale is the Cenomanian-Turonian OAE 2, which is characterized by a pronounced positive excursion in carbon-isotope records and the important accumulation of organic-rich sediments. The section at Gongzha (Tibet) and the sections at Barranca and Axaxacualco (Mexico) are located in remote parts of the Tethys, and show δ13C records, which are well correlated with those of classical Tethyan sections. Both sections, however, do not exhibit the presence of organic-rich sediments. Phosphorus Mass Accumulation Rates (PMAR) in Tibet show a pattern similar to that observed in the Tethys by Mort et al. (2007), which suggests enhanced Ρ regeneration during the OAE 2 time interval, though there is no evidence for anoxic conditions in Tibet. Ρ appears here to have been mainly driven by detrital influx and sea-level fluctuations. The sections at Barranca and Axaxacualco show that the Mexican carbonate platform persisted during this anoxic event, which allowed the evolution of platform fauna otherwise not present in Tethyan sections. The persistence of this carbonate platform close to the Caribbean Igneous Plateau, which is thought to have released bio-limiting metals, is explained by local uplift which delayed the drowning of the platform and a specific oceanic circulation that permitted the preservation of oligotrophic conditions in the area. The Coniacian-Santonian OAE (OAE3) appears to have been more dependent on local conditions than OAE2. The presence of black shales associated with OAE3 appear to have been restricted to shallow-water settings and epicontinental seas in areas located around the Atlantic Ocean. The sections at Olazagutia (Spain), and Ten Mile - Arbor Park (USA), two potential Global Boundary Stratotype Sections and Points (GSSP) sites, are devoid of organic-rich sediments and lack a δ13C positive excursion around the C-S boundary. The Gabal Ekma section (Sinai, Egypt) exhibits accumulations of organic-rich sediments, in addition to phosphorite bone beds layers, which may have been linked to an epicontinental upwelling zone and/or storm inputs. Our data suggest that OAE 3 is rarely expressed by truly anoxic conditions and seems to have been linked to local conditions rather than global paleoenvironmental change. The evidence for detrital-P being the likely cause of Ρ fluctuations during the OAEs studied here does not negate the idea that anoxia was the principal driver of these fluctuations in the western Tethys. However, an explanation is required as to why the Ρ accumulation signatures are mirrored in both oxic and anoxic sedimentary successions. 'Eustatic/climatic' and 'productivity/anoxic' models may have both operated simultaneously in different parts of the world depending on local conditions, both producing similar trends in Ρ accumulation. - Selon Jenkyns (2010), les événements anoxiques océaniques enregistrent de profonds changements dans le climat et la paléoceanographie de la planète et représente des perturbations majeures du cycle du carbone. L'un des plus étudiés à l'échelle mondiale est l'ΟΑΕ2 du Cénomanien-Turonien, qui est caractérisé par une très forte excursion positive des isotopes du carbone et une importante accumulation de sédiments riche en matière organique. La section de Gongzha (Tibet) et les sections de Barranca et Axaxcualco (Mexique) sont situées aux confins de la Téthys, et enregistrent une courbe isotopique en δ13C parfaitement corrélable avec les sections téthysiennes, mais ne montre pas d'accumulation de black shales. Le taux de phosphore en accumulation de masses (PMAR) au Tibet montre un pattern similaire observé également par Mort et al. (2007) dans la Téthys, suggérant un model de régénération du Ρ durant l'anoxie, cependant aucune conditions anoxiques régnent dans la région du Tibet. Ρ apparaît donc principalement guidé par le détritisme et les fluctuations du niveau marin. Les sections de Barranca et d'Axaxacualco montrent que la plateforme carbonatée mexicaine persiste durant cet événement anoxique, et permet le développement d'une faune de plateforme qui n'est pas présente dans les sections téthysiennes. La persistance de cette plateforme carbonatée si proche du plateau Caribéen, qui est connu pour le relâchement de métaux bio-limitant, peut être expliqué par un soulèvement tectonique local qui inhibe l'ennoiement de la plateforme et une circulation océanique spécifique qui permet la préservation de conditions oligotrophiques dans cette région. L'événement anoxique océanique du Coniacien-Santonien apparaît plus dépendant des conditions locales que pour l'ΟΑΕ2. Les black shales associés à POAE3 sont restreints aux zones situées autour de l'océan Atlantique et plus particulièrement aux eaux peu profondes et épicontinentales. Les sections d'Olazagutia (Espagne), Ten Mile Creek et Arbor Park (USA), qui sont deux potentielles sections GSSP (Sections de stratotype de limite globaux et de points), ne montre pas d'accumulation de black shales et pas de forte excursion positive en δ13C autour de la limite C-S. La section de Gabal Ekma (Sinai, Egypte) montre des accumulations de black shales, en plus des couches de phosphorites et d'accumulation d'os (« bone beds »), vraisemblablement lié à des zones active d'upwelling épicontinentale et/ou d'apport de tempêtes. Nos données suggèrent que l'OAE 3 est rarement exprimé par de vraies conditions anoxiques et semble être plus lié à des conditions plus locales que des changements paléo-environnementaux globaux, comme observés pour le Cénomanien- Turonien. Les arguments pour un modèle lié au phosphore détritique qui serait la cause des fluctuations du phosphore total durant les OAEs, n'écartent pas l'idée que l'anoxie est la principale cause de ces fluctuations dans les sections riches en matière organique de l'Ouest téthysien. Cependant une explication est nécessaire pour comprendre pourquoi la signature de l'accumulation du phosphore est semblable dans les successions sédimentaires déposées dans des conditions oxygénées et anoxiques. Les modèles « Eustatisme/Climat » et « Productivité/anoxie » ont simultanément opéré dans les différentes parties du monde dépendant de conditions locales, et ont produit des tendances similaires en accumulation de phosphore.

Building the niche through time: using 13,000 years of data to predict the effects of climate change on three tree species in Europe

Aim Species distribution models (SDMs) based on current species ranges underestimate the potential distribution when projected in time and/or space. A multi-temporal model calibration approach has been suggested as an alternative, and we evaluate this using 13,000 years of data. Location Europe. Methods We used fossil-based records of presence for Picea abies, Abies alba and Fagus sylvatica and six climatic variables for the period 13,000 to 1000yr bp. To measure the contribution of each 1000-year time step to the total niche of each species (the niche measured by pooling all the data), we employed a principal components analysis (PCA) calibrated with data over the entire range of possible climates. Then we projected both the total niche and the partial niches from single time frames into the PCA space, and tested if the partial niches were more similar to the total niche than random. Using an ensemble forecasting approach, we calibrated SDMs for each time frame and for the pooled database. We projected each model to current climate and evaluated the results against current pollen data. We also projected all models into the future. Results Niche similarity between the partial and the total-SDMs was almost always statistically significant and increased through time. SDMs calibrated from single time frames gave different results when projected to current climate, providing evidence of a change in the species realized niches through time. Moreover, they predicted limited climate suitability when compared with the total-SDMs. The same results were obtained when projected to future climates. Main conclusions The realized climatic niche of species differed for current and future climates when SDMs were calibrated considering different past climates. Building the niche as an ensemble through time represents a way forward to a better understanding of a species' range and its ecology in a changing climate.

Climatic and anthropogenic influence on the stable isotope record from bulk carbonates and ostracodes in Lake Neuchatel, Switzerland, during the last two millennia

Lake Neuchatel is a medium sized, hard-water lake, lacking varved sediments, situated in the western Swiss Lowlands at the foot of the Jura Mountains. Stable isotope data (delta(18)O and delta(13)C) from both bulk carbonate and ostracode calcite in an 81 cm long, radiocarbon-dated sediment core represent the last 1500 years of Lake Neuchatel's environmental history. Comparison between this isotopic and other palaeolimnologic data (mineralogical, geochemical, palynological, etc.) helps to differentiate between anthropogenic and natural factors most recently affecting the lake. An increase in lacustrine productivity (450-650AD ca), inferred from the positive trend in delta(13)C values of bulk carbonate, is related to medieval forest clearances and the associated nutrient budget changes. A negative trend in both the bulk carbonate and ostracode calcite delta(18)O values between approximately 1300 and 1500AD, is tentatively interpreted as due to a cooling in mean air temperature at the transition from the Medieval Warm Period to the Little Ice Age. Negative trends in bulk carbonate delta(18)O and delta(13)C values through the uppermost sediments, which have no equivalent in ostracode calcite isotopic values, are concomitant with the recent onset of eutrophication in the lake. Isotopic disequilibrium during calcite precipitation, probably due to kinetic factors in periods of high productivity is postulated as the mechanism to explain the associated negative isotopic trends, although the effect of a shift of the calcite precipitation towards the warmer months cannot be excluded.

Importance of abiotic stress as a range-limit determinant for European plants: insights from species' responses to climatic gradients

Aim We examined whether species occurrences are primarily limited by physiological tolerance in the abiotically more stressful end of climatic gradients (the asymmetric abiotic stress limitation (AASL) hypothesis) and the geographical predictions of this hypothesis: abiotic stress mainly determines upper-latitudinal and upper-altitudinal species range limits, and the importance of abiotic stress for these range limits increases the further northwards and upwards a species occurs. Location Europe and the Swiss Alps. Methods The AASL hypothesis predicts that species have skewed responses to climatic gradients, with a steep decline towards the more stressful conditions. Based on presence-absence data we examined the shape of plant species responses (measured as probability of occurrence) along three climatic gradients across latitudes in Europe (1577 species) and altitudes in the Swiss Alps (284 species) using Huisman-Olff-Fresco, generalized linear and generalized additive models. Results We found that almost half of the species from Europe and one-third from the Swiss Alps showed responses consistent with the predictions of the AASL hypothesis. Cold temperatures and a short growing season seemed to determine the upper-latitudinal and upper-altitudinal range limits of up to one-third of the species, while drought provided an important constraint at lower-latitudinal range limits for up to one-fifth of the species. We found a biome-dependent influence of abiotic stress and no clear support for abiotic stress as a stronger upper range-limit determinant for species with higher latitudinal and altitudinal distributions. However, the overall influence of climate as a range-limit determinant increased with latitude. Main conclusions Our results support the AASL hypothesis for almost half of the studied species, and suggest that temperature-related stress controls the upper-latitudinal and upper-altitudinal range limits of a large proportion of these species, while other factors including drought stress may be important at the lower range limits.

Divergent and narrower climatic niches characterize polyploid species of European primroses in Primula sect. Aleuritia

Aim It is hypothesized that the ecological niches of polyploids should be both distinct and broader than those of diploids - characteristics that might have allowed the successful colonization of open habitats by polyploids during the Pleistocene glacial cycles. Here, we test these hypotheses by quantifying and comparing the ecological niches and niche breadths of a group of European primroses. Location Europe. Methods We gathered georeferenced data of four related species in Primula sect. Aleuritia at different ploidy levels (diploid, tetraploid, hexaploid and octoploid) and used seven bioclimatic variables to quantify niche overlap between species by applying a series of univariate and multivariate analyses combined with modelling techniques. We also employed permutation-based tests to evaluate niche similarity between the four species. Niche breadth for each species was evaluated both in the multivariate environmental space and in geographical space. Results The four species differed significantly from each other in mono-dimensional comparisons of climatological variables and occupied distinct habitats in the multi-dimensional environmental space. The majority of the permutation-based tests either indicated that the four species differed significantly in their habitat preferences and ecological niches or did not support significant niche similarity. Furthermore, our results revealed narrower niche breadths and geographical ranges in species of P. sect. Aleuritia at higher ploidy levels. Main conclusions The detected ecological differentiation between the four species of P. sect. Aleuritia at different ploidy levels is consistent with the hypothesis that polyploids occupy distinct ecological niches that differ from those of their diploid relative. Contrary to expectations, we find that polyploid species of P. sect. Aleuritia occupy narrower environmental and geographical spaces than their diploid relative. These results on the ecological niches of closely related polyploid and diploid species highlight factors that potentially contribute to the evolution and distribution of polyploid species.

Impact of model complexity on cross-temporal transferability in Maxent species distribution models: An assessment using paleobotanical data

Maximum entropy modeling (Maxent) is a widely used algorithm for predicting species distributions across space and time. Properly assessing the uncertainty in such predictions is non-trivial and requires validation with independent datasets. Notably, model complexity (number of model parameters) remains a major concern in relation to overfitting and, hence, transferability of Maxent models. An emerging approach is to validate the cross-temporal transferability of model predictions using paleoecological data. In this study, we assess the effect of model complexity on the performance of Maxent projections across time using two European plant species (Alnus giutinosa (L.) Gaertn. and Corylus avellana L) with an extensive late Quaternary fossil record in Spain as a study case. We fit 110 models with different levels of complexity under present time and tested model performance using AUC (area under the receiver operating characteristic curve) and AlCc (corrected Akaike Information Criterion) through the standard procedure of randomly partitioning current occurrence data. We then compared these results to an independent validation by projecting the models to mid-Holocene (6000 years before present) climatic conditions in Spain to assess their ability to predict fossil pollen presence-absence and abundance. We find that calibrating Maxent models with default settings result in the generation of overly complex models. While model performance increased with model complexity when predicting current distributions, it was higher with intermediate complexity when predicting mid-Holocene distributions. Hence, models of intermediate complexity resulted in the best trade-off to predict species distributions across time. Reliable temporal model transferability is especially relevant for forecasting species distributions under future climate change. Consequently, species-specific model tuning should be used to find the best modeling settings to control for complexity, notably with paleoecological data to independently validate model projections. For cross-temporal projections of species distributions for which paleoecological data is not available, models of intermediate complexity should be selected.