Modelling habitat-suitability using museum collections: an example with three sympatric Apodemus species from the Alps

Aim, Location Although the alpine mouse Apodemus alpicola has been given species status since 1989, no distribution map has ever been constructed for this endemic alpine rodent in Switzerland. Based on redetermined museum material and using the Ecological-Niche Factor Analysis (ENFA), habitat-suitability maps were computed for A. alpicola, and also for the co-occurring A. flavicollis and A. sylvaticus. Methods In the particular case of habitat suitability models, classical approaches (GLMs, GAMs, discriminant analysis, etc.) generally require presence and absence data. The presence records provided by museums can clearly give useful information about species distribution and ecology and have already been used for knowledge-based mapping. In this paper, we apply the ENFA which requires only presence data, to build a habitat-suitability map of three species of Apodemus on the basis of museum skull collections. Results Interspecific niche comparisons showed that A. alpicola is very specialized concerning habitat selection, meaning that its habitat differs unequivocally from the average conditions in Switzerland, while both A. flavicollis and A. sylvaticus could be considered as 'generalists' in the study area. Main conclusions Although an adequate sampling design is the best way to collect ecological data for predictive modelling, this is a time and money consuming process and there are cases where time is simply not available, as for instance with endangered species conservation. On the other hand, museums, herbariums and other similar institutions are treasuring huge presence data sets. By applying the ENFA to such data it is possible to rapidly construct a habitat suitability model. The ENFA method not only provides two key measurements regarding the niche of a species (i.e. marginality and specialization), but also has ecological meaning, and allows the scientist to compare directly the niches of different species.

Habitat suitability of Anopheles vector species and association with human malaria in the Atlantic Forest in south-eastern Brazil

Every year, autochthonous cases of Plasmodium vivax malaria occur in low-endemicity areas of Vale do Ribeira in the south-eastern part of the Atlantic Forest, state of São Paulo, where Anopheles cruzii and Anopheles bellator are considered the primary vectors. However, other species in the subgenus Nyssorhynchus of Anopheles (e.g., Anopheles marajoara) are abundant and may participate in the dynamics of malarial transmission in that region. The objectives of the present study were to assess the spatial distribution of An. cruzii, An. bellator and An. marajoara and to associate the presence of these species with malaria cases in the municipalities of the Vale do Ribeira. Potential habitat suitability modelling was applied to determine both the spatial distribution of An. cruzii, An. bellator and An. marajoara and to establish the density of each species. Poisson regression was utilized to associate malaria cases with estimated vector densities. As a result, An. cruzii was correlated with the forested slopes of the Serra do Mar, An. bellator with the coastal plain and An. marajoara with the deforested areas. Moreover, both An. marajoara and An. cruzii were positively associated with malaria cases. Considering that An. marajoara was demonstrated to be a primary vector of human Plasmodium in the rural areas of the state of Amapá, more attention should be given to the species in the deforested areas of the Atlantic Forest, where it might be a secondary vector.

Overcoming the rare species modelling paradox: test of a novel hierarchical framework with an Iberian endemic plant

Rare species have restricted geographic ranges, habitat specialization, and/or small population sizes. Datasets on rare species distribution usually have few observations, limited spatial accuracy and lack of valid absences; conversely they provide comprehensive views of species distributions allowing to realistically capture most of their realized environmental niche. Rare species are the most in need of predictive distribution modelling but also the most difficult to model. We refer to this contrast as the "rare species modelling paradox" and propose as a solution developing modelling approaches that deal with a sufficiently large set of predictors, ensuring that statistical models aren't overfitted. Our novel approach fulfils this condition by fitting a large number of bivariate models and averaging them with a weighted ensemble approach. We further propose that this ensemble forecasting is conducted within a hierarchic multi-scale framework. We present two ensemble models for a test species, one at regional and one at local scale, each based on the combination of 630 models. In both cases, we obtained excellent spatial projections, unusual when modelling rare species. Model results highlight, from a statistically sound approach, the effects of multiple drivers in a same modelling framework and at two distinct scales. From this added information, regional models can support accurate forecasts of range dynamics under climate change scenarios, whereas local models allow the assessment of isolated or synergistic impacts of changes in multiple predictors. This novel framework provides a baseline for adaptive conservation, management and monitoring of rare species at distinct spatial and temporal scales.

Which is the optimal sampling strategy for habitat suitability modelling

Designing an efficient sampling strategy is of crucial importance for habitat suitability modelling. This paper compares four such strategies, namely, 'random', 'regular', 'proportional-stratified' and 'equal -stratified'- to investigate (1) how they affect prediction accuracy and (2) how sensitive they are to sample size. In order to compare them, a virtual species approach (Ecol. Model. 145 (2001) 111) in a real landscape, based on reliable data, was chosen. The distribution of the virtual species was sampled 300 times using each of the four strategies in four sample sizes. The sampled data were then fed into a GLM to make two types of prediction: (1) habitat suitability and (2) presence/ absence. Comparing the predictions to the known distribution of the virtual species allows model accuracy to be assessed. Habitat suitability predictions were assessed by Pearson's correlation coefficient and presence/absence predictions by Cohen's K agreement coefficient. The results show the 'regular' and 'equal-stratified' sampling strategies to be the most accurate and most robust. We propose the following characteristics to improve sample design: (1) increase sample size, (2) prefer systematic to random sampling and (3) include environmental information in the design'

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.

There and back again? Combining habitat suitability modelling and connectivity analyses to assess a potential return of the otter to Switzerland

After a steady decline in the early 20th century, several terrestrial carnivore species have recently recovered in Western Europe, either through reintroductions or natural recolonization. Because of the large space requirements of these species and potential conflicts with human activities, ensuring their recovery requires the implementation of conservation and management measures that address the environmental, landscape and social dimensions of the problem. Few examples exist of such integrated management. Taking the case of the otter (Lutra lutra) in Switzerland, we propose a multi-step approach that allows to (1) identify areas with potentially suitable habitat, (2) evaluate their connectivity, (3) verify the potentiality of the species recolonization from populations in neighbouring countries. We showed that even though suitable habitat is available for the species and the level of structural connectivity within Switzerland is satisfactory, the level of connectivity with neighbouring populations is crucial to prioritize strategies that favour the species recovery in the field. This research is the first example integrating habitat suitability and connectivity assessment at different scales with other factors in a multi-step assessment for species recovery.

Do habitat suitability models reliably predict the recovery areas of threatened species?

1. Identifying those areas suitable for recolonization by threatened species is essential to support efficient conservation policies. Habitat suitability models (HSM) predict species' potential distributions, but the quality of their predictions should be carefully assessed when the species-environment equilibrium assumption is violated.2. We studied the Eurasian otter Lutra lutra, whose numbers are recovering in southern Italy. To produce widely applicable results, we chose standard HSM procedures and looked for the models' capacities in predicting the suitability of a recolonization area. We used two fieldwork datasets: presence-only data, used in the Ecological Niche Factor Analyses (ENFA), and presence-absence data, used in a Generalized Linear Model (GLM). In addition to cross-validation, we independently evaluated the models with data from a recolonization event, providing presences on a previously unoccupied river.3. Three of the models successfully predicted the suitability of the recolonization area, but the GLM built with data before the recolonization disagreed with these predictions, missing the recolonized river's suitability and badly describing the otter's niche. Our results highlighted three points of relevance to modelling practices: (1) absences may prevent the models from correctly identifying areas suitable for a species spread; (2) the selection of variables may lead to randomness in the predictions; and (3) the Area Under Curve (AUC), a commonly used validation index, was not well suited to the evaluation of model quality, whereas the Boyce Index (CBI), based on presence data only, better highlighted the models' fit to the recolonization observations.4. For species with unstable spatial distributions, presence-only models may work better than presence-absence methods in making reliable predictions of suitable areas for expansion. An iterative modelling process, using new occurrences from each step of the species spread, may also help in progressively reducing errors.5. Synthesis and applications. Conservation plans depend on reliable models of the species' suitable habitats. In non-equilibrium situations, such as the case for threatened or invasive species, models could be affected negatively by the inclusion of absence data when predicting the areas of potential expansion. Presence-only methods will here provide a better basis for productive conservation management practices.

What drives invasibility? A multi-model inference test and spatial modelling of alien plant species richness patterns in northern Portugal

Understanding and anticipating biological invasions can focus either on traits that favour species invasiveness or on features of the receiving communities, habitats or landscapes that promote their invasibility. Here, we address invasibility at the regional scale, testing whether some habitats and landscapes are more invasible than others by fitting models that relate alien plant species richness to various environmental predictors. We use a multi-model information-theoretic approach to assess invasibility by modelling spatial and ecological patterns of alien invasion in landscape mosaics and testing competing hypotheses of environmental factors that may control invasibility. Because invasibility may be mediated by particular characteristics of invasiveness, we classified alien species according to their C-S-R plant strategies. We illustrate this approach with a set of 86 alien species in Northern Portugal. We first focus on predictors influencing species richness and expressing invasibility and then evaluate whether distinct plant strategies respond to the same or different groups of environmental predictors. We confirmed climate as a primary determinant of alien invasions and as a primary environmental gradient determining landscape invasibility. The effects of secondary gradients were detected only when the area was sub-sampled according to predictions based on the primary gradient. Then, multiple predictor types influenced patterns of alien species richness, with some types (landscape composition, topography and fire regime) prevailing over others. Alien species richness responded most strongly to extreme land management regimes, suggesting that intermediate disturbance induces biotic resistance by favouring native species richness. Land-use intensification facilitated alien invasion, whereas conservation areas hosted few invaders, highlighting the importance of ecosystem stability in preventing invasions. Plants with different strategies exhibited different responses to environmental gradients, particularly when the variations of the primary gradient were narrowed by sub-sampling. Such differential responses of plant strategies suggest using distinct control and eradication approaches for different areas and alien plant groups.

The influence of environmental spatial structure on the life history traits and diversity of species in a metacommunity

Several models have been proposed to understand how so many species can coexist in ecosystems. Despite evidence showing that natural habitats are often patchy and fragmented, these models rarely take into account environmental spatial structure. In this study we investigated the influence of spatial structure in habitat and disturbance regime upon species' traits and species' coexistence in a metacommunity. We used a population-based model to simulate competing species in spatially explicit landscapes. The species traits we focused on were dispersal ability, competitiveness, reproductive investment and survival rate. Communities were characterized by their species richness and by the four life-history traits averaged over all the surviving species. Our results show that spatial structure and disturbance have a strong influence on the equilibrium life-history traits within a metacommunity. In the absence of disturbance, spatially structured landscapes favour species investing more in reproduction, but less in dispersal and survival. However, this influence is strongly dependent on the disturbance rate, pointing to an important interaction between spatial structure and disturbance. This interaction also plays a role in species coexistence. While spatial structure tends to reduce diversity in the absence of disturbance, the tendency is reversed when disturbance occurs. In conclusion, the spatial structure of communities is an important determinant of their diversity and characteristic traits. These traits are likely to influence important ecological properties such as resistance to invasion or response to climate change, which in turn will determine the fate of ecosystems facing the current global ecological crisis.

Fire management, climate change and their interacting effects on birds in complex Mediterranean landscapes: dynamic distribution modelling of an early-successional species - the near-threatened Dartford Warbler (Sylvia undata)

The current challenge in a context of major environmental changes is to anticipate the responses of species to future landscape and climate scenarios. In the Mediterranean basin, climate change is one the most powerful driving forces of fire dynamics, with fire frequency and impact having markedly increased in recent years. Species distribution modelling plays a fundamental role in this challenge, but better integration of available ecological knowledge is needed to adequately guide conservation efforts. Here, we quantified changes in habitat suitability of an early-succession bird in Catalonia, the Dartford Warbler (Sylvia undata) ― globally evaluated as Near Threatened in the IUCN Red List. We assessed potential changes in species distributions between 2000 and 2050 under different fire management and climate change scenarios and described landscape dynamics using a spatially-explicit fire-succession model that simulates fire impacts in the landscape and post-fire regeneration (MEDFIRE model). Dartford Warbler occurrence data were acquired at two different spatial scales from: 1) the Atlas of European Breeding Birds (EBCC) and 2) Catalan Breeding Bird Atlas (CBBA). Habitat suitability was modelled using five widely-used modelling techniques in an ensemble forecasting framework. Our results indicated considerable habitat suitability losses (ranging between 47% and 57% in baseline scenarios), which were modulated to a large extent by fire regime changes derived from fire management policies and climate changes. Such result highlighted the need for taking the spatial interaction between climate changes, fire-mediated landscape dynamics and fire management policies into account for coherently anticipating habitat suitability changes of early succession bird species. We conclude that fire management programs need to be integrated into conservation plans to effectively preserve sparsely forested and early succession habitats and their associated species in the face of global environmental change.

A spatial modelling framework for assessing climate change impacts on freshwater ecosystems: Response of brown trout (Salmo trutta L.) biomass to warming water temperature

Mountain regions worldwide are particularly sensitive to on-going climate change. Specifically in the Alps in Switzerland, the temperature has increased twice as fast than in the rest of the Northern hemisphere. Water temperature closely follows the annual air temperature cycle, severely impacting streams and freshwater ecosystems. In the last 20 years, brown trout (Salmo trutta L) catch has declined by approximately 40-50% in many rivers in Switzerland. Increasing water temperature has been suggested as one of the most likely cause of this decline. Temperature has a direct effect on trout population dynamics through developmental and disease control but can also indirectly impact dynamics via food-web interactions such as resource availability. We developed a spatially explicit modelling framework that allows spatial and temporal projections of trout biomass using the Aare river catchment as a model system, in order to assess the spatial and seasonal patterns of trout biomass variation. Given that biomass has a seasonal variation depending on trout life history stage, we developed seasonal biomass variation models for three periods of the year (Autumn-Winter, Spring and Summer). Because stream water temperature is a critical parameter for brown trout development, we first calibrated a model to predict water temperature as a function of air temperature to be able to further apply climate change scenarios. We then built a model of trout biomass variation by linking water temperature to trout biomass measurements collected by electro-fishing in 21 stations from 2009 to 2011. The different modelling components of our framework had overall a good predictive ability and we could show a seasonal effect of water temperature affecting trout biomass variation. Our statistical framework uses a minimum set of input variables that make it easily transferable to other study areas or fish species but could be improved by including effects of the biotic environment and the evolution of demographical parameters over time. However, our framework still remains informative to spatially highlight where potential changes of water temperature could affect trout biomass. (C) 2015 Elsevier B.V. All rights reserved.-

Aspects of spatial and habitat ecology of multiple Anopheles species (Diptera: Culicidae): malaria vectors in the highlands and foothills of Ecuador

The resurgence of malaria in highland regions of Africa, Oceania and recently in South America underlines the importance of the study of the ecology of highland mosquito vectors of malaria. Since the incidence of malaria is limited by the distribution of its vectors, the purpose of this PhD thesis was to examine aspects of the ecology of Anopheles mosquitoes in the Andes of Ecuador, South America. A historical literature and archival data review (Chapter 2) indicated that Anopheles pseudopunctipennis transmitted malaria in highland valleys of Ecuador prior to 1950, although it was eliminated through habitat removal and the use of chemical insecticides. Other anopheline species were previously limited to low-altitude regions, except in a few unconfirmed cases. A thorough larval collection effort (n=438 attempted collection sites) in all road-accessible parts of Ecuador except for the lowland Amazon basin was undertaken between 2008 - 2010 (Chapter 3). Larvae were identified morphologically and using molecular techniques (mitochondrial COl gene), and distribution maps indicated that all five species collected (Anopheles albimanus, An. pseudopunctipennis, Anopheles punctimacula, Anopheles oswaldoi s.l. and Anopheles eiseni) were more widespread throughout highland regions than previously recorded during the 1940s, with higher maximum altitudes for all except An. pseudopunctipennis (1541 m, 1930 m, 1906 m, 1233 m and 1873 m, respectively). During larval collections, to characterize species-specific larval habitat, a variety of abiotic and biotic habitat parameters were measured and compared between species-present and species-absent sites using chi-square tests and stepwise binary logistic regression analyses (Chapter 4). An. albimanus was significantly associated with permanent pools with sand substrates and An. pseudopunctipennis with gravel and boulder substrates. Both species were significantly associated with floating cyanobacterial mats and warmer temperatures, which may limit their presence in cooler highland regions. Anopheles punctimacula was collected more often than expected from algae-free, shaded pools with higher-than-average calculated dissolved oxygen. Anopheles oswaldoi s.l., the species occurring on the Amazonian side of the Andes, was associated with permanent, anthropogenic habitats such as roadside ditches and ponds. To address the hypothesis that human land use change is responsible for the emergence of multiple highland Anopheles species by creating larval habitat, common land uses in the western Andes were surveyed for standing water and potential larval habitat suitability (Chapter 5). Rivers and road edges provided large amounts of potentially suitable anopheline habitat in the western Andes, while cattle pasture also created potentially suitable habitat in irrigation canals and watering ponds. Other common land uses surveyed (banana farms, sugarcane plantations, mixed tree plantations, and empty lots) were usually established on steep slopes and had very little standing water present. Using distribution and larval habitat data, a GIS-based larval habitat distribution model for the common western species was constructed in ArcGIS v.l 0 (ESRI 2010) using derived data layers from field measurements and other sources (Chapter 6). The additive model predicted 76.4 - 97.9% of the field-observed collection localities of An. albimanus, An. pseudopunctipennis and An. punctimacula, although it could not accurately distinguish between species-absent and speciespresent sites due to its coarse scale. The model predicted distributional expansion and/or shift of one or more anopheline species into the following highland valleys with climate warming: Mira/Chota, Imbabura province, Tumbaco, Pichincha province, Pallatanga and Sibambe, Chimborazo province, and Yungilla, Azuay province. These valleys may serve as targeted sites of future monitoring to prevent highland epidemics of malaria. The human perceptions of malaria and mosquitoes in relation to land management practices were assessed through an interview-based survey (n=262) in both highlands and lowlands, of male and female land owners and managers of five property types (Chapter 7). Although respondents had a strong understanding of where the disease occurs in their own country and of the basic relationship among standing water, mosquitoes and malaria, about half of respondents in potential risk areas denied the current possibility of malaria infection on their own property. As well, about half of respondents with potential anopheline larval habitat did not report its presence, likely due to a highly specific definition of suitable mosquito habitat. Most respondents who are considered at risk of malaria currently use at least one type of mosquito bite prevention, most commonly bed nets. In conclusion, this interdisciplinary thesis examines the occurrence of Anopheles species in the lowland transition area and highlands in Ecuador, from a historic, geographic, ecological and sociological perspective.