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.

GENETIC VARIATION, LOCAL ADAPTATION AND POPULATION STRUCTURE IN NORTH AMERICAN RED OAK SPECIES, QUERCUS RUBRA L. AND Q. ELLIPSOIDALIS E. J. HILL

Forest trees, like oaks, rely on high levels of genetic variation to adapt to varying environmental conditions. Thus, genetic variation and its distribution are important for the long-term survival and adaptability of oak populations. Climate change is projected to lead to increased drought and fire events as well as a northward migration of tree species, including oaks. Additionally, decline in oak regeneration has become increasingly concerning since it may lead to decreased gene flow and increased inbreeding levels. This will in turn lead to lowered levels of genetic diversity, negatively affecting the growth and survival of populations. At the same time, populations at the species’ distribution edge, like those in this study, could possess important stores of genetic diversity and adaptive potential, while also being vulnerable to climatic or anthropogenic changes. A survey of the level and distribution of genetic variation and identification of potentially adaptive genes is needed since adaptive genetic variation is essential for their long-term survival. Oaks possess a remarkable characteristic in that they maintain their species identity and specific environmental adaptations despite their propensity to hybridize. Thus, in the face of interspecific gene flow, some areas of the genome remain differentiated due to selection. This characteristic allows the study of local environmental adaptation through genetic variation analyses. Furthermore, using genic markers with known putative functions makes it possible to link those differentiated markers to potential adaptive traits (e.g., flowering time, drought stress tolerance). Demographic processes like gene flow and genetic drift also play an important role in how genes (including adaptive genes) are maintained or spread. These processes are influenced by disturbances, both natural and anthropogenic. An examination of how genetic variation is geographically distributed can display how these genetic processes and geographical disturbances influence genetic variation patterns. For example, the spatial clustering of closely related trees could promote inbreeding with associated negative effects (inbreeding depression), if gene flow is limited. In turn this can have negative consequences for a species’ ability to adapt to changing environmental conditions. In contrast, interspecific hybridization may also allow the transfer of genes between species that increase their adaptive potential in a changing environment. I have studied the ecologically divergent, interfertile red oaks, Quercus rubra and Q. ellipsoidalis, to identify genes with potential roles in adaptation to abiotic stress through traits such as drought tolerance and flowering time, and to assess the level and distribution of genetic variation. I found evidence for moderate gene flow between the two species and low interspecific genetic differences at most genetic markers (Lind and Gailing 2013). However, the screening of genic markers with potential roles in phenology and drought tolerance led to the identification of a CONSTANS-like (COL) gene, a candidate gene for flowering time and growth. This marker, located in the coding region of the gene, was highly differentiated between the two species in multiple geographical areas, despite interspecific gene flow, and may play a role in reproductive isolation and adaptive divergence between the two species (Lind-Riehl et al. 2014). Since climate change could result in a northward migration of trees species like oaks, this gene could be important in maintaining species identity despite increased contact zones between species (e.g., increased gene flow). Finally I examined differences in spatial genetic structure (SGS) and genetic variation between species and populations subjected to different management strategies and natural disturbances. Diverse management activities combined with various natural disturbances as well as species specific life history traits influenced SGS patterns and inbreeding levels (Lind-Riehl and Gailing submitted).

Regional habitat needs of a nationally listed species, Canada Warbler (Cardellina canadensis), in Alberta, Canada

Understanding factors that affect the distribution and abundance of species is critical to developing effective management plans for conservation. Our goal was to quantify the distribution and abundance of Canada Warbler (Cardellina canadensis), a threatened old-forest associate in Alberta, Canada. The Canada Warbler has declined across its range, including in Alberta where habitat loss and alteration from urban expansion, forestry, and energy development are changing the forest landscape. We used 110,427 point count survey visits from 32,287 unique survey stations to model local-level (150-m radius circular buffers) and stand-level (564-m radius circular buffers) habitat associations of the Canada Warbler. We found that habitat supporting higher densities of Canada Warblers was locally concentrated yet broadly distributed across Alberta’s boreal forest region. Canada Warblers were most commonly associated with older deciduous forest at the local scale, particularly near small, incised streams, and greater amounts of deciduous forest at the stand scale. Predicted density was lower in other forest types and younger age classes measured at the local scale. There was little evidence that local-scale fragmentation (i.e., edges created by linear features) influenced Canada Warbler abundance. However, current forestry practices in the province likely will reduce the availability of Canada Warbler habitat over time by cutting old deciduous forest stands. Our results suggest that conservation efforts aimed at Canada Warbler focus on retaining large stands of old deciduous forest, specifically stands adjacent to streams, by increasing the width of deciduous retention buffers along streams during harvest and increasing the size and number of old forest residual patches in harvested stands.

Atlas versus range maps: robustness of chorological relationships to distribution data types in European mammals

Aim Chorological relationships describe the patterns of distributional overlap among species. In addition to revealing biogeographical structure, the resulting clusters of species with similar geographical distributions can serve as natural units in conservation planning. Here, we assess the extent to which temporal, methodological and taxonomical differences in the source of species’ distribution data can affect the relationships that are found.

Favourable areas for expansion and reintroduction of Iberian lynx accounting for distribution trends and genetic diversity of the European rabbit

Although on a local scale Iberian lynx distribution is determined by the availability of prey rabbits, recent modelling analyses have uncovered broad-scale disagreements between these two species’ distribution trends. These analyses showed also that the lynx had become restricted to only a fraction of the rabbit’s genetic diversity, and that this could be jeopardising its survival in the face of environmental hazards and uncertainty. In the present paper, a follow-up was carried out through the building of lynx and rabbit distribution models based on the most recent Spanish mammal atlas. Environmental favourability values for lynx and rabbit were positively correlated within the lynx’s current distribution area, but they were negatively correlated within the total Spanish area where lynx used to occur in the 1980’s. Environmental favourability for rabbits was significantly higher where lynx maintains reproductive populations than where it recently disappeared, indicating that rabbit favourability plays an important role and can be a good predictor of lynx persistence. The lynx and rabbit models were extrapolated to predict favourable areas for both species in Spain as well as in Portugal, on the original scale of the distribution data (10x10 km) and on a 100 times finer spatial resolution (1x1 km). The lynx and rabbit models were also combined through fuzzy logic to forecast the potential for lynx occurrence incorporating information on favourable areas for its main prey. Several areas are proposed as favourable for lynx expansion or re-introduction,

How well can models predict changes in species distributions? A 13-year-old otter model revisited

Species distribution and ecological niche models are increasingly used in biodiversity management and conservation. However, one thing that is important but rarely done is to follow up on the predictive performance of these models over time, to check if their predictions are fulfilled and maintain accuracy, or if they apply only to the set in which they were produced. In 2003, a distribution model of the Eurasian otter (Lutra lutra) in Spain was published, based on the results of a country-wide otter survey published in 1998. This model was built with logistic regression of otter presence-absence in UTM 10 km2 cells on a diverse set of environmental, human and spatial variables, selected according to statistical criteria. Here we evaluate this model against the results of the most recent otter survey, carried out a decade later and after a significant expansion of the otter distribution area in this country. Despite the time elapsed and the evident changes in this species’ distribution, the model maintained a good predictive capacity, considering both discrimination and calibration measures. Otter distribution did not expand randomly or simply towards vicinity areas,m but specifically towards the areas predicted as most favourable by the model based on data from 10 years before. This corroborates the utility of predictive distribution models, at least in the medium term and when they are made with robust methods and relevant predictor variables.

The importance of ecological scale for wildlife conservation in naturally fragmented environments: A case study of the Brush-Tailed Rock-Wallaby (Petrogale Penicillata)

Determining the ecologically relevant spatial scales for predicting species occurrences is an important concept when determining species–environment relationships. Therefore species distribution modelling should consider all ecologically relevant spatial scales. While several recent studies have addressed this problem in artificially fragmented landscapes, few studies have researched relevant ecological scales for organisms that also live in naturally fragmented landscapes. This situation is exemplified by the Australian rock-wallabies’ preference for rugged terrain and we addressed the issue of scale using the threatened brush-tailed rock-wallaby (Petrogale penicillata) in eastern Australia. We surveyed for brush-tailed rock-wallabies at 200 sites in southeast Queensland, collecting potentially influential site level and landscape level variables. We applied classification trees at either scale to capture a hierarchy of relationships between the explanatory variables and brush-tailed rock-wallaby presence/absence. Habitat complexity at the site level and geology at the landscape level were the best predictors of where we observed brush-tailed rock-wallabies. Our study showed that the distribution of the species is affected by both site scale and landscape scale factors, reinforcing the need for a multi-scale approach to understanding the relationship between a species and its environment. We demonstrate that careful design of data collection, using coarse scale spatial datasets and finer scale field data, can provide useful information for identifying the ecologically relevant scales for studying species–environment relationships. Our study highlights the need to determine patterns of environmental influence at multiple scales to conserve specialist species such as the brush-tailed rock-wallaby in naturally fragmented landscapes.

Incorporating uncertainty in environmental models informed by imagery

In this thesis, the issue of incorporating uncertainty for environmental modelling informed by imagery is explored by considering uncertainty in deterministic modelling, measurement uncertainty and uncertainty in image composition. Incorporating uncertainty in deterministic modelling is extended for use with imagery using the Bayesian melding approach. In the application presented, slope steepness is shown to be the main contributor to total uncertainty in the Revised Universal Soil Loss Equation. A spatial sampling procedure is also proposed to assist in implementing Bayesian melding given the increased data size with models informed by imagery. Measurement error models are another approach to incorporating uncertainty when data is informed by imagery. These models for measurement uncertainty, considered in a Bayesian conditional independence framework, are applied to ecological data generated from imagery. The models are shown to be appropriate and useful in certain situations. Measurement uncertainty is also considered in the context of change detection when two images are not co-registered. An approach for detecting change in two successive images is proposed that is not affected by registration. The procedure uses the Kolmogorov-Smirnov test on homogeneous segments of an image to detect change, with the homogeneous segments determined using a Bayesian mixture model of pixel values. Using the mixture model to segment an image also allows for uncertainty in the composition of an image. This thesis concludes by comparing several different Bayesian image segmentation approaches that allow for uncertainty regarding the allocation of pixels to different ground components. Each segmentation approach is applied to a data set of chlorophyll values and shown to have different benefits and drawbacks depending on the aims of the analysis.

The Bayesian conditional independence model for measurement error: applications in ecology

The measurement error model is a well established statistical method for regression problems in medical sciences, although rarely used in ecological studies. While the situations in which it is appropriate may be less common in ecology, there are instances in which there may be benefits in its use for prediction and estimation of parameters of interest. We have chosen to explore this topic using a conditional independence model in a Bayesian framework using a Gibbs sampler, as this gives a great deal of flexibility, allowing us to analyse a number of different models without losing generality. Using simulations and two examples, we show how the conditional independence model can be used in ecology, and when it is appropriate.

Growing biodiverse carbon-rich forests

Regrowing forests on cleared land is a key strategy to achieve both biodiversity conservation and climate change mitigation globally. Maximizing these co-benefits, however, remains theoretically and technically challenging because of the complex relationship between carbon sequestration and biodiversity in forests, the strong influence of climate variability and landscape position on forest development, the large number of restoration strategies possible, and long time-frames needed to declare success. Through the synthesis of three decades of knowledge on forest dynamics and plant functional traits combined with decision science, we demonstrate that we cannot always maximize carbon sequestration by simply increasing the functional trait diversity of trees planted. The relationships between plant functional diversity, carbon sequestration rates above-ground and in the soil are dependent on climate and landscape positions. We show how to manage ‘identities’ and ‘complementarities’ between plant functional traits in order to achieve systematically maximal co-benefits in various climate and landscape contexts. We provide examples of optimal planting and thinning rules that satisfy this ecological strategy and guide the restoration of forests that are rich in both carbon and plant functional diversity. Our framework provides the first mechanistic approach for generating decision-making rules that can be used to manage forests for multiple objectives, and supports joined carbon credit and biodiversity conservation initiatives, such as Reducing Emissions from Deforestation and forest Degradation REDD+. The decision framework can also be linked to species distribution models and socio-economic models in order to find restoration solutions that maximize simultaneously biodiversity, carbon stocks and other ecosystem services across landscapes. Our study provides the foundation for developing and testing cost-effective and adaptable forest management rules to achieve biodiversity, carbon sequestration and other socio-economic co-benefits under global change.

Are alien fish a reliable indicator of river health?

1. The ability of many introduced fish species to thrive in degraded aquatic habitats and their potential to impact on aquatic ecosystem structure and function suggest that introduced fish may represent both a symptom and a cause of decline in river health and the integrity of native aquatic communities. 2. The varying sensitivities of many commonly introduced fish species to degraded stream conditions, the mechanism and reason for their introduction and the differential susceptibility of local stream habitats to invasion because of the environmental and biological characteristics of the receiving water body, are all confounding factors that may obscure the interpretation of patterns of introduced fish species distribution and abundance and therefore their reliability as indicators of river health. 3. In the present study, we address the question of whether alien fish (i.e. those species introduced from other countries) are a reliable indicator of the health of streams and rivers in south-eastern Queensland, Australia. We examine the relationships of alien fish species distributions and indices of abundance and biomass with the natural environmental features, the biotic characteristics of the local native fish assemblages and indicators of anthropogenic disturbance at a large number of sites subject to varying sources and intensities of human impact. 4. Alien fish species were found to be widespread and often abundant in south-eastern Queensland rivers and streams, and the five species collected were considered to be relatively tolerant to river degradation, making them good candidate indicators of river health. Variation in alien species indices was unrelated to the size of the study sites, the sampling effort expended or natural environmental gradients. The biological resistance of the native fish fauna was not concluded to be an important factor mediating invasion success by alien species. Variation in alien fish indices was, however, strongly related to indicators of disturbance intensity describing local in-stream habitat and riparian degradation, water quality and surrounding land use, particularly the amount of urban development in the catchment. 5. Potential confounding factors that may influence the likelihood of introduction and successful establishment of an alien species and the implications of these factors for river bioassessment are discussed. We conclude that the potentially strong impact that many alien fish species can have on the biological integrity of natural aquatic ecosystems, together with their potential to be used as an initial basis to find out other forms of human disturbance impacts, suggest that some alien species (particularly species from the family Poeciliidae) can represent a reliable 'first cut' indicator of river health.

A large proportion of asymptomatic Plasmodium infections with low and sub-microscopic parasite densities in the low transmission setting of Temotu Province, Solomon Islands : challenges for malaria diagnostics in an elimination setting

Background Many countries are scaling up malaria interventions towards elimination. This transition changes demands on malaria diagnostics from diagnosing ill patients to detecting parasites in all carriers including asymptomatic infections and infections with low parasite densities. Detection methods suitable to local malaria epidemiology must be selected prior to transitioning a malaria control programme to elimination. A baseline malaria survey conducted in Temotu Province, Solomon Islands in late 2008, as the first step in a provincial malaria elimination programme, provided malaria epidemiology data and an opportunity to assess how well different diagnostic methods performed in this setting. Methods During the survey, 9,491 blood samples were collected and examined by microscopy for Plasmodium species and density, with a subset also examined by polymerase chain reaction (PCR) and rapid diagnostic tests (RDTs). The performances of these diagnostic methods were compared. Results A total of 256 samples were positive by microscopy, giving a point prevalence of 2.7%. The species distribution was 17.5% Plasmodium falciparum and 82.4% Plasmodium vivax. In this low transmission setting, only 17.8% of the P. falciparum and 2.9% of P. vivax infected subjects were febrile (≥38°C) at the time of the survey. A significant proportion of infections detected by microscopy, 40% and 65.6% for P. falciparum and P. vivax respectively, had parasite density below 100/μL. There was an age correlation for the proportion of parasite density below 100/μL for P. vivax infections, but not for P. falciparum infections. PCR detected substantially more infections than microscopy (point prevalence of 8.71%), indicating a large number of subjects had sub-microscopic parasitemia. The concordance between PCR and microscopy in detecting single species was greater for P. vivax (135/162) compared to P. falciparum (36/118). The malaria RDT detected the 12 microscopy and PCR positive P. falciparum, but failed to detect 12/13 microscopy and PCR positive P. vivax infections. Conclusion Asymptomatic malaria infections and infections with low and sub-microscopic parasite densities are highly prevalent in Temotu province where malaria transmission is low. This presents a challenge for elimination since the large proportion of the parasite reservoir will not be detected by standard active and passive case detection. Therefore effective mass screening and treatment campaigns will most likely need more sensitive assays such as a field deployable molecular based assay.

A continental-scale tool for acoustic identification of European bats

Summary 1. Acoustic methods are used increasingly to survey and monitor bat populations. However, the use of acoustic methods at continental scales can be hampered by the lack of standardized and objective methods to identify all species recorded. This makes comparable continent-wide monitoring difficult, impeding progress towards developing biodiversity indicators, transboundary conservation programmes and monitoring species distribution changes. 2. Here we developed a continental-scale classifier for acoustic identification of bats, which can be used throughout Europe to ensure objective, consistent and comparable species identifications. We selected 1350 full-spectrum reference calls from a set of 15 858 calls of 34 European species, from EchoBank, a global echolocation call library. We assessed 24 call parameters to evaluate how well they distinguish between species and used the 12 most useful to train a hierarchy of ensembles of artificial neural networks to distinguish the echolocation calls of these bat species. 3. Calls are first classified to one of five call-type groups, with a median accuracy of 97·6%. The median species-level classification accuracy is 83·7%, providing robust classification for most European species, and an estimate of classification error for each species. 4. These classifiers were packaged into an online tool, iBatsID, which is freely available, enabling anyone to classify European calls in an objective and consistent way, allowing standardized acoustic identification across the continent. 5. Synthesis and applications. iBatsID is the first freely available and easily accessible continental- scale bat call classifier, providing the basis for standardized, continental acoustic bat monitoring in Europe. This method can provide key information to managers and conservation planners on distribution changes and changes in bat species activity through time.

Extinction risk in cloud forest fragments under climate change and habitat loss

Aim: To quantify the consequences of major threats to biodiversity, such as climate and land-use change, it is important to use explicit measures of species persistence, such as extinction risk. The extinction risk of metapopulations can be approximated through simple models, providing a regional snapshot of the extinction probability of a species. We evaluated the extinction risk of three species under different climate change scenarios in three different regions of the Mexican cloud forest, a highly fragmented habitat that is particularly vulnerable to climate change. Location: Cloud forests in Mexico. Methods: Using Maxent, we estimated the potential distribution of cloud forest for three different time horizons (2030, 2050 and 2080) and their overlap with protected areas. Then, we calculated the extinction risk of three contrasting vertebrate species for two scenarios: (1) climate change only (all suitable areas of cloud forest through time) and (2) climate and land-use change (only suitable areas within a currently protected area), using an explicit patch-occupancy approximation model and calculating the joint probability of all populations becoming extinct when the number of remaining patches was less than five. Results: Our results show that the extent of environmentally suitable areas for cloud forest in Mexico will sharply decline in the next 70 years. We discovered that if all habitat outside protected areas is transformed, then only species with small area requirements are likely to persist. With habitat loss through climate change only, high dispersal rates are sufficient for persistence, but this requires protection of all remaining cloud forest areas. Main conclusions: Even if high dispersal rates mitigate the extinction risk of species due to climate change, the synergistic impacts of changing climate and land use further threaten the persistence of species with higher area requirements. Our approach for assessing the impacts of threats on biodiversity is particularly useful when there is little time or data for detailed population viability analyses. © 2013 John Wiley & Sons Ltd.

Optimizing search strategies for invasive pests: Learn before you leap

Strategic searching for invasive pests presents a formidable challenge for conservation managers. Limited funding can necessitate choosing between surveying many sites cursorily, or focussing intensively on fewer sites. While existing knowledge may help to target more likely sites, e.g. with species distribution models (maps), this knowledge is not flawless and improving it also requires management investment. 2.In a rare example of trading-off action against knowledge gain, we combine search coverage and accuracy, and its future improvement, within a single optimisation framework. More specifically we examine under which circumstances managers should adopt one of two search-and-control strategies (cursory or focussed), and when they should divert funding to improving knowledge, making better predictive maps that benefit future searches. 3.We use a family of Receiver Operating Characteristic curves to reflect the quality of maps that direct search efforts. We demonstrate our framework by linking these to a logistic model of invasive spread such as that for the red imported fire ant Solenopsis invicta in south-east Queensland, Australia. 4.Cursory widespread searching is only optimal if the pest is already widespread or knowledge is poor, otherwise focussed searching exploiting the map is preferable. For longer management timeframes, eradication is more likely if funds are initially devoted to improving knowledge, even if this results in a short-term explosion of the pest population. 5.Synthesis and applications. By combining trade-offs between knowledge acquisition and utilization, managers can better focus - and justify - their spending to achieve optimal results in invasive control efforts. This framework can improve the efficiency of any ecological management that relies on predicting occurrence. © 2010 The Authors. Journal of Applied Ecology © 2010 British Ecological Society.