fuzzySim: Fuzzy similarity in species distributions

fuzzySim is an R package for calculating fuzzy similarity in species occurrence patterns. It includes functions for data preparation, such as converting species lists (long format) to presence-absence tables (wide format), obtaining unique abbreviations of species names, or transposing (parts of) complex data frames; and sample data sets for providing practical examples. It can convert binary presence-absence to fuzzy occurrence data, using e.g. trend surface analysis, inverse distance interpolation or prevalence-independent environmental favourability modelling, for multiple species simultaneously. It then calculates fuzzy similarity among (fuzzy) species distributions and/or among (fuzzy) regional species compositions. Currently available similarity indices are Jaccard, Sørensen, Simpson, and Baroni-Urbani & Buser.

Discrimination capacity in species distribution models depends on the representativeness of the environmental domain

Aim When faced with dichotomous events, such as the presence or absence of a species, discrimination capacity (the ability to separate the instances of presence from the instances of absence) is usually the only characteristic that is assessed in the evaluation of the performance of predictive models. Although neglected, calibration or reliability (how well the estimated probability of presence represents the observed proportion of presences) is another aspect of the performance of predictive models that provides important information. In this study, we explore how changes in the distribution of the probability of presence make discrimination capacity a context-dependent characteristic of models. For the first time,we explain the implications that ignoring the context dependence of discrimination can have in the interpretation of species distribution models.

New measures for assessing model equilibrium and prediction mismatch in species distribution models

Models based on species distributions are widely used and serve important purposes in ecology, biogeography and conservation. Their continuous predictions of environmental suitability are commonly converted into a binary classification of predicted (or potential) presences and absences, whose accuracy is then evaluated through a number of measures that have been the subject of recent reviews. We propose four additional measures that analyse observation-prediction mismatch from a different angle – namely, from the perspective of the predicted rather than the observed area – and add to the existing toolset of model evaluation methods. We explain how these measures can complete the view provided by the existing measures, allowing further insights into distribution model predictions. We also describe how they can be particularly useful when using models to forecast the spread of diseases or of invasive species and to predict modifications in species’ distributions under climate and land-use change

Distribution modelling of wild rabbit hunting yields in its original area (S Iberian Peninsula)

In this work we used the information of the Annual Hunting Reports (AHRs) to obtain a high-resolution model of the potential favourableness for wild rabbit harvesting in Andalusia (southern Spain), using environmental and land-use variables as predictors. We analysed 32,134 AHRs from the period 1993/2001 reported by 6049 game estates to estimate the average hunting yields of wild rabbit in each Andalusian municipality (n5771). We modelled the favourableness for obtaining good hunting yields using stepwise logistic regression on a set of climatic, orographical, land use, and vegetation variables. The favourability equation was used to create a downscaled image representing the favourableness of obtaining good hunting yields for the wild rabbit in 161 km squares in Andalusia, using the Idrisi Image Calculator. The variables that affected hunting yields of wild rabbit were altitude, dry wood crops (mainly olive groves, almond groves, and vineyards), temperature, pasture, slope, and annual number of frost days. The 161 km squares with high favourableness values are scattered throughout the territory, which seems to be caused mainly by the effect of vegetation. Finally, we obtained quality categories for the territory by combining the probability values given by logistic regression with those of the environmental favourability function.

Patterns of Cactaceae species distribution in a protected area in the semiarid caatinga biome of North-Eastern Brazil.

2016

Spatial scale effects of sampling on the interpolation of species distribution models in the southwestern Amazon.

Knowledge of the geographical distribution of timber tree species in the Amazon is still scarce. This is especially true at the local level, thereby limiting natural resource management actions. Forest inventories are key sources of information on the occurrence of such species. However, areas with approved forest management plans are mostly located near access roads and the main industrial centers. The present study aimed to assess the spatial scale effects of forest inventories used as sources of occurrence data in the interpolation of potential species distribution models. The occurrence data of a group of six forest tree species were divided into four geographical areas during the modeling process. Several sampling schemes were then tested applying the maximum entropy algorithm, using the following predictor variables: elevation, slope, exposure, normalized difference vegetation index (NDVI) and height above the nearest drainage (HAND). The results revealed that using occurrence data from only one geographical area with unique environmental characteristics increased both model overfitting to input data and omission error rates. The use of a diagonal systematic sampling scheme and lower threshold values led to improved model performance. Forest inventories may be used to predict areas with a high probability of species occurrence, provided they are located in forest management plan regions representative of the environmental range of the model projection area.

Niche conservatism in non-native birds in Europe: niche unfilling rather than niche expansion

Aim Niche conservatism, or the extent to which niches are conserved across space and time, is of special concern for the study of non-native species as it underlies predictions of invasion risk. Based on the occurrence of 28 non-native birds in Europe, we assess to what extent Grinnellian realized niches are conserved during invasion, formulate hypotheses to explain the variation in observed niche changes and test how well species distribution models can predict non-native bird occurrence in Europe. Location Europe. Methods To quantify niche changes, a recent method that applies kernel smoothers to densities of species occurrence in a gridded environmental space was used. This corrects for differences in the availability of environments between study areas and allows discrimination between 'niche expansion' into environments new to the species and 'niche unfilling', whereby the species only partially fills its niche in the invaded range. Predictions of non-native bird distribution in Europe were generated using several distribution modelling techniques. Results Niche overlap between native and non-native bird populations is low, but niche changes are smaller for species having a higher propagule pressure and that were introduced longer ago. Non-native birds in Europe occupy a subset of the environments they inhabit in their native ranges. Niche expansion into novel environments is rare for most species, allowing species distribution models to accurately predict invasion risk. Main conclusions Because of the recent nature of most bird introductions, species occupy only part of the suitable environments available in the invaded range. This signals that apart from purely ecological factors, patterns of niche conservatism may also be contingent on population-specific historical factors. These results also suggest that many claims of niche differences may be due to a partial filling of the native niche in the invaded range and thus do not represent true niche changes.

The complex history of the olive tree: from Late Quaternary diversification of Mediterranean lineages to primary domestication in the northern Levant.

The location and timing of domestication of the olive tree, a key crop in Early Mediterranean societies, remain hotly debated. Here, we unravel the history of wild olives (oleasters), and then infer the primary origins of the domesticated olive. Phylogeography and Bayesian molecular dating analyses based on plastid genome profiling of 1263 oleasters and 534 cultivated genotypes reveal three main lineages of pre-Quaternary origin. Regional hotspots of plastid diversity, species distribution modelling and macrofossils support the existence of three long-term refugia; namely the Near East (including Cyprus), the Aegean area and the Strait of Gibraltar. These ancestral wild gene pools have provided the essential foundations for cultivated olive breeding. Comparison of the geographical pattern of plastid diversity between wild and cultivated olives indicates the cradle of first domestication in the northern Levant followed by dispersals across the Mediterranean basin in parallel with the expansion of civilizations and human exchanges in this part of the world.

Multiple refugia and barriers explain the phylogeography of the Valais shrew, Sorex antinorii (Mammalia: Soricomorpha), Sorex antinorii (Mammalia: Soricomorpha)

The aim of the present study was to investigate the genetic structure of the Valais shrew (Sorex antinorii) by a combined phylogeographical and landscape genetic approach, and thereby to infer the locations of glacial refugia and establish the influence of geographical barriers. We sequenced part of the mitochondrial cytochrome b (cyt b) gene of 179 individuals of S. antinorii sampled across the entire species' range. Six specimens attributed to S. arunchi were included in the analysis. The phylogeographical pattern was assessed by Bayesian molecular phylogenetic reconstruction, population genetic analyses, and a species distribution modelling (SDM)-based hindcasting approach. We also used landscape genetics (including isolation-by-resistance) to infer the determinants of current intra-specific genetic structure. The phylogeographical analysis revealed shallow divergence among haplotypes and no clear substructure within S. antinorii. The starlike structure of the median-joining network is consistent with population expansion from a single refugium, probably located in the Apennines. Long branches observed on the same network also suggest that another refugium may have existed in the north-eastern part of Italy. This result is consistent with SDM, which also suggests several habitable areas for S. antinorii in the Italian peninsula during the LGM. Therefore S. antinorii appears to have occupied disconnected glacial refugia in the Italian peninsula, supporting previous data for other species showing multiple refugia within southern refugial areas. By coupling genetic analyses and SDM, we were able to infer how past climatic suitability contributed to genetic divergence of populations. The genetic differentiation shown in the present study does not support the specific status of S. arunchi.

Understanding the ecological niche to elucidate spatial strategies of the southernmost Tupinambis lizards

Understanding factors that shape ranges of species is central in evolutionary biology. Species distribution models have become important tools to test biogeographical, ecological and evolutionary hypotheses. Moreover, from an ecological and evolutionary perspective, these models help to elucidate the spatial strategies of species at a regional scale. We modelled species distributions of two phylogenetically, geographically and ecologically close Tupinambis species (Teiidae) that occupy the southernmost area of the genus distribution in South America. We hypothesized that similarities between these species might have induced spatial strategies at the species level, such as niche differentiation and divergence of distribution patterns at a regional scale. Using logistic regression and MaxEnt we obtained species distribution models that revealed interspecific differences in habitat requirements, such as environmental temperature, precipitation and altitude. Moreover, the models obtained suggest that although the ecological niches of Tupinambis merianae and T. rufescens are different, these species might co-occur in a large contact zone. We propose that niche plasticity could be the mechanism enabling their co-occurrence. Therefore, the approach used here allowed us to understand the spatial strategies of two Tupinambis lizards at a regional scale.

Very high-resolution digital elevation models: are multi-scale-derived variables ecologically relevant?

1. Digital elevation models (DEMs) are often used in landscape ecology to retrieve elevation or first derivative terrain attributes such as slope or aspect in the context of species distribution modelling. However, DEM-derived variables are scale-dependent and, given the increasing availability of very high-resolution (VHR) DEMs, their ecological relevancemust be assessed for different spatial resolutions. 2. In a study area located in the Swiss Western Alps, we computed VHR DEMs-derived variables related to morphometry, hydrology and solar radiation. Based on an original spatial resolution of 0.5 m, we generated DEM-derived variables at 1, 2 and 4 mspatial resolutions, applying a Gaussian Pyramid. Their associations with local climatic factors, measured by sensors (direct and ambient air temperature, air humidity and soil moisture) as well as ecological indicators derived fromspecies composition, were assessed with multivariate generalized linearmodels (GLM) andmixed models (GLMM). 3. Specific VHR DEM-derived variables showed significant associations with climatic factors. In addition to slope, aspect and curvature, the underused wetness and ruggedness indices modelledmeasured ambient humidity and soilmoisture, respectively. Remarkably, spatial resolution of VHR DEM-derived variables had a significant influence on models' strength, with coefficients of determination decreasing with coarser resolutions or showing a local optimumwith a 2 mresolution, depending on the variable considered. 4. These results support the relevance of using multi-scale DEM variables to provide surrogates for important climatic variables such as humidity, moisture and temperature, offering suitable alternatives to direct measurements for evolutionary ecology studies at a local scale.

Pollination services at risk: Bee habitats will decrease owing to climate change in Brazil

Native bees are important providers of pollination services, but there are cumulative evidences of their decline. Global changes such as habitat losses, invasions of exotic species and climate change have been suggested as the main causes of the decline of pollinators. In this study, the influence of climate change on the distribution of 10 species of Brazilian bees was estimated with species distribution modelling. We used Maxent algorithm (maximum entropy) and two different scenarios, an optimistic and a pessimistic, to the years 2050 and 2080. We also evaluated the percentage reduction of species habitat based on the future scenarios of climate change through Geographic Information System (GIS). Results showed that the total area of suitable habitats decreased for all species but one under the different future scenarios. The greatest reductions in habitat area were found for Melipona bicolor bicolor and Melipona scutellaris, which occur predominantly in areas related originally to Atlantic Moist Forest. The species analysed have been reported to be pollinators of some regional crops and the consequence of their decrease for these crops needs further clarification. (C) 2012 Elsevier B.V. All rights reserved.

Latitudinal shift in thermal niche breadth results from thermal release during a climate-mediated range expansion

Funded by Biodiversity and Ecosystem Services in a Changing Climate Wenner-Gren Foundation Swedish Research Council The Royal Swedish Academy of Sciences Stiftelsen Anna-Greta Holger Crafoords Fund The Crafoord Foundation

Latitudinal shift in thermal niche breadth results from thermal release during a climate-mediated range expansion

Funded by Biodiversity and Ecosystem Services in a Changing Climate Wenner-Gren Foundation Swedish Research Council The Royal Swedish Academy of Sciences Stiftelsen Anna-Greta Holger Crafoords Fund The Crafoord Foundation

Adaptive management, international co-operation and planning for marine conservation hotspots in a changing climate

Acknowledgements This work received funding from the MASTS pooling initiative (The Marine Alliance for Science and Technology for Scotland) and their support is gratefully acknowledged. MASTS is funded by the Scottish Funding Council (Grant reference HR09011) and contributing institutions.