Ecological Impacts of Alien Species: Quantification, Scope, Caveats, and Recommendations

Despite intensive research during the past decade on the effects of alien species, invasion science still lacks the capacity to accurately predict the impacts of those species and, therefore, to provide timely advice to managers on where limited resources should be allocated. This capacity has been limited partly by the context-dependent nature of ecological impacts, research highly skewed toward certain taxa and habitat types, and the lack of standardized methods for detecting and quantifying impacts. We review different strategies, including specific experimental and observational approaches, for detecting and quantifying the ecological impacts of alien species. These include a four-way experimental plot design for comparing impact studies of different organisms. Furthermore, we identify hypothesis-driven parameters that should be measured at invaded sites to maximize insights into the nature of the impact. We also present strategies for recognizing high-impact species. Our recommendations provide a foundation for developing systematic quantitative measurements to allow comparisons of impacts across alien species, sites, and time.

Range expansion and comparative habitat use of insular, congeneric lagomorphs: invasive European hares Lepus europaeus and endemic Irish hares Lepus timidus hibernicus

The European hare (Lepus europaeus) has declined throughout its native range but invaded numerous regions where it has negatively impacted native wildlife. In southern Sweden, it replaces the native mountain hare (L. timidus) through competition and hybridisation. We investigated temporal change in the invasive range of the European hare in Ireland, and compared its habitat use with the endemic Irish hare (L. timidus hibernicus). The range of the European hare was three times larger and its core range twice as large in 2012–2013 than in 2005. Its rate of radial range expansion was 0.73 km year−1 with its introduction estimated to have occurred ca. 1970. Both species utilised improved and rough grasslands and exhibited markedly similar regression coefficients with almost every land cover variable examined. Irish hares were associated with low fibre and high sugar content grass (good quality grazing) whilst the invader had a greater tolerance for low quality forage. European hares were associated with habitat patch edge density, suggesting it may be more suited to using hedgerows as diurnal resting sites than the Irish hare. Consequently, the invader had a wider niche breadth than the native but their niche overlap was virtually complete. Given the impact of the European hare on native species elsewhere, and its apparent pre-adaption for improved grasslands interspersed with arable land (a habitat that covers 70 % of Ireland), its establishment and range expansion poses a significant threat to the ecological security of the endemic Irish hare, particularly given their ecological similarities.

Choosing and using diversity indices: Insights for ecological applications from the German Biodiversity Exploratories

Biodiversity, a multidimensional property of natural systems, is difficult to quantify partly because of the multitude of indices proposed for this purpose. Indices aim to describe general properties of communities that allow us to compare different regions, taxa, and trophic levels. Therefore, they are of fundamental importance for environmental monitoring and conservation, although there is no consensus about which indices are more appropriate and informative. We tested several common diversity indices in a range of simple to complex statistical analyses in order to determine whether some were better suited for certain analyses than others. We used data collected around the focal plant Plantago lanceolata on 60 temperate grassland plots embedded in an agricultural landscape to explore relationships between the common diversity indices of species richness (S), Shannon's diversity (H'), Simpson's diversity (D1), Simpson's dominance (D2), Simpson's evenness (E), and Berger–Parker dominance (BP). We calculated each of these indices for herbaceous plants, arbuscular mycorrhizal fungi, aboveground arthropods, belowground insect larvae, and P. lanceolata molecular and chemical diversity. Including these trait-based measures of diversity allowed us to test whether or not they behaved similarly to the better studied species diversity. We used path analysis to determine whether compound indices detected more relationships between diversities of different organisms and traits than more basic indices. In the path models, more paths were significant when using H', even though all models except that with E were equally reliable. This demonstrates that while common diversity indices may appear interchangeable in simple analyses, when considering complex interactions, the choice of index can profoundly alter the interpretation of results. Data mining in order to identify the index producing the most significant results should be avoided, but simultaneously considering analyses using multiple indices can provide greater insight into the interactions in a system.

Invading and expanding: range dynamics and ecological consequences of the greater white-toothed shrew (Crocidura russula) invasion in Ireland

Establishing how invasive species impact upon pre-existing species is a fundamental question in ecology and conservation biology. The greater white-toothed shrew (Crocidura russula) is an invasive species in Ireland that was first recorded in 2007 and which, according to initial data, may be limiting the abundance/distribution of the pygmy shrew (Sorex minutus), previously Ireland’s only shrew species. Because of these concerns, we undertook an intensive live-trapping survey (and used other data from live-trapping, sightings and bird of prey pellets/nest inspections collected between 2006 and 2013) to model the distribution and expansion of C. russula in Ireland and its impacts on Ireland’s small mammal community. The main distribution range of C. russula was found to be approximately 7,600 km2 in 2013, with established outlier populations suggesting that the species is dispersing with human assistance within the island. The species is expanding rapidly for a small mammal, with a radial expansion rate of 5.5 km/yr overall (2008–2013), and independent estimates from live-trapping in 2012–2013 showing rates of 2.4–14.1 km/yr, 0.5–7.1 km/yr and 0–5.6 km/yr depending on the landscape features present. S. minutus is negatively associated with C. russula. S. minutus is completely absent at sites where C. russula is established and is only present at sites at the edge of and beyond the invasion range of C. russula. The speed of this invasion and the homogenous nature of the Irish landscape may mean that S. minutus has not had sufficient time to adapt to the sudden appearance of C. russula. This may mean the continued decline/disappearance of S. minutus as C. russula spreads throughout the island.

Pseudomonas aeruginosa exhibits frequent recombination, but only a limited association between genotype and ecological setting

Pseudomonas aeruginosa is an opportunistic pathogen and an important cause of infection, particularly amongst cystic fibrosis (CF) patients. While specific strains capable of patient-to-patient transmission are known, many infections appear to be caused by unique and unrelated strains. There is a need to understand the relationship between strains capable of colonising the CF lung and the broader set of P. aeruginosa isolates found in natural environments. Here we report the results of a multilocus sequence typing (MLST)-based study designed to understand the genetic diversity and population structure of an extensive regional sample of P. aeruginosa isolates from South East Queensland, Australia. The analysis is based on 501 P. aeruginosa isolates obtained from environmental, animal and human (CF and non-CF) sources with particular emphasis on isolates from the Lower Brisbane River and isolates from CF patients obtained from the same geographical region. Overall, MLST identified 274 different sequence types, of which 53 were shared between one or more ecological settings. Our analysis revealed a limited association between genotype and environment and evidence of frequent recombination. We also found that genetic diversity of P. aeruginosa in Queensland, Australia was indistinguishable from that of the global P. aeruginosa population. Several CF strains were encountered frequently in multiple ecological settings; however, the most frequently encountered CF strains were confined to CF patients. Overall, our data confirm a non-clonal epidemic structure and indicate that most CF strains are a random sample of the broader P. aeruginosa population. The increased abundance of some CF strains in different geographical regions is a likely product of chance colonisation events followed by adaptation to the CF lung and horizontal transmission among patients.

Environmental filtering vs. resource-based niche partitioning in diverse soil animal assemblages

Terrestrial invertebrates constitute most of described animal biodiversity and soil is a major reservoir of this diversity. In the classical attempt to understand the processes supporting biodiversity, ecologists are currently seeking to unravel the differential roles of environmental filtering and competition for resources in niche partitioning processes: these processes are in principle distinct although they may act simultaneously, interact at multiple spatial and temporal scales, and are often confounded in studies of soil communities. We used a novel combination of methods based on stable isotopes and trait analysis to resolve these processes in diverse oribatid mite assemblages at spatial
scales at which competition for resources could in principle be a major driver. We also used a null model approach based on a general neutral model of beta diversity. A large and significant fraction of community variation was explainable in terms of linear and periodic spatial structures in the distribution of organic C, N and soil structure: species were clearly arranged along an environmental, spatially structured gradient. However, competition related trait differences did not map onto the distances separating species along the environmental gradient and neutral models provided a satisfying approximation of beta diversity patterns. The results represent the first robust evidence
that in very diverse soil arthropod assemblages resource-based niche partitioning plays a minor role while environmental filtering remains a fundamental driver of species distribution.

A new statistical framework for the quantification of covariate associations with species distributions

Identifying processes that shape species geographical ranges is a prerequisite for understanding environmental change. Currently, species distribution modelling methods do not offer credible statistical tests of the relative influence of climate factors and typically ignore other processes (e.g. biotic interactions and dispersal limitation). We use a hierarchical model fitted with Markov Chain Monte Carlo to combine ecologically plausible niche structures using regression splines to describe unimodal but potentially skewed response terms. We apply spatially explicit error terms that account for (and may help identify) missing variables. Using three example distributions of European bird species, we map model results to show sensitivity to change in each covariate. We show that the overall strength of climatic association differs between species and that each species has considerable spatial variation in both the strength of the climatic association and the sensitivity to climate change. Our methods are widely applicable to many species distribution modelling problems and enable accurate assessment of the statistical importance of biotic and abiotic influences on distributions.

Spatial turnover in the global avifauna

Despite its wide implications for many ecological issues, the global pattern of spatial turnover in the occurrence of species has been little studied, unlike the global pattern of species richness. Here, using a database on the breeding distributions of birds, we present the first global maps of variation in spatial turnover for an entire taxonomic class, a pattern that has to date remained largely a matter of conjecture, based on theoretical expectations and extrapolation of inconsistent patterns from different biogeographic realms. We use these maps to test four predictions from niche theory as to the form that this variation should take, namely that turnover should increase with species richness, towards lower latitudes, and with the steepness of environmental gradients and that variation in turnover is determined principally by rare (restricted) species. Contrary to prediction, we show that turnover is high both in areas of extremely low and high species richness, does not increase strongly towards the tropics, and is related both to average environmental conditions and spatial variation in those conditions. These results are closely associated with a further important and novel finding, namely that global patterns of spatial turnover are driven principally by widespread species rather than the restricted ones. This complements recent demonstrations that spatial patterns of species richness are also driven principally by widespread species, and thus provides an important contribution towards a unified model of how terrestrial biodiversity varies both within and between the Earth's major land masses.

The imprint of the geographical, evolutionary and ecological context on species-area relationships

Species-area relationships (SAR) are fundamental in the understanding of biodiversity patterns and of critical importance for predicting species extinction risk worldwide. Despite the enormous attention given to SAR in the form of many individual analyses, little attempt has been made to synthesize these studies. We conducted a quantitative meta-analysis of 794 SAR, comprising a wide span of organisms, habitats and locations. We identified factors reflecting both pattern-based and dynamic approaches to SAR and tested whether these factors leave significant imprints on the slope and strength of SAR. Our analysis revealed that SAR are significantly affected by variables characterizing the sampling scheme, the spatial scale, and the types of organisms or habitats involved. We found that steeper SAR are generated at lower latitudes and by larger organisms. SAR varied significantly between nested and independent sampling schemes and between major ecosystem types, but not generally between the terrestrial and the aquatic realm. Both the fit and the slope of the SAR were scale-dependent. We conclude that factors dynamically regulating species richness at different spatial scales strongly affect the shape of SAR. We highlight important consequences of this systematic variation in SAR for ecological theory, conservation management and extinction risk predictions.

Ecological dynamics of extinct species in empty habitat networks. 1. The role of habitat pattern and quantity, stochasticity and dispersal

We examined a remnant host plant (Primula veris L.) habitat network that was last inhabited by the rare butterfly Hamearis lucina L. in north Wales in 1943, to assess the relative contribution of several spatial parameters to its regional extinction. We first examined relationships between P. veris characteristics and H. lucina eggs in surviving H. lucina populations, and used these to predict the suitability and potential carrying capacity of the habitat network in north Wales. This resulted in an estimate of roughly 4500 eggs (ca 227 adults). We developed a discrete space, discrete time metapopulation model to evaluate the relative contribution of dispersal distance, habitat and environmental stochasticity as possible causes of extinction. We simulated the potential persistence of the butterfly in the current network as well as in three artificial (historical and present) habitat networks that differed in the quantity (current and X3) and fragmentation of the habitat (current and aggregated). We identified that reduced habitat quantity and increased isolation would have increased the probability of regional extinction, in conjunction with environmental stochasticity and H. lucina's dispersal distance. This general trend did not change in a qualitative manner when we modified the ability of dispersing females to stay in, and find suitable habitats (by changing the size of the grid cells used in the model). Contrary to most metapopulation model predictions, system persistence declined with increasing migration rate, suggesting that the mortality of migrating individuals in fragmented landscapes may pose significant risks to system-wide persistence. Based on model predictions for the present landscape we argue that a major programme of habitat restoration would be required for a re-established metapopulation to persist for > 100 years.

Ecological dynamics of extinct species in empty habitat networks. 2. The role of host plant dynamics

This paper explores the relative effects of host plant dynamics and butterfly-related parameters on butterfly persistence. It considers an empty habitat network where a rare butterfly (Cupido minimus) became extinct in 1939 in part of its historical range in north Wales, UK. Surviving populations of the butterfly in southern Britain were visited to assess use of its host plant (Anthyllis vulneraria) in order to calibrate habitat suitability and carrying capacity in the empty network in north Wales. These data were used to deduce that only a portion ( similar to 19%) of the host plant network from north Wales was likely to be highly suitable for oviposition. Nonetheless, roughly 65,460 eggs (3273 adult equivalents) could be expected to be laid in north Wales, were the empty network to be populated at the same levels as observed on comparable plants in surviving populations elsewhere. Simulated metapopulations of C. minimus in the empty network revealed that time to extinction and patch occupancy were significantly influenced by carrying capacity, butterfly mean dispersal distance and environmental stochasticity, although for most reasonable parameter values, the model system persisted. Simulation outputs differed greatly when host plant dynamics was incorporated into the modelled butterfly dynamics. Cupido minimus usually went extinct when host plant were at low densities. In these simulations host plant dynamics appeared to be the most important determinant of the butterfly's regional extirpation. Modelling the outcome of a reintroduction programme to C. minimus variation at high quality locations, revealed that 65% of systems survived at least 100 years. Given the current amount of resources of the north Wales landscape, the persistence of C. minimus under a realistic reintroduction programme has a good chance of being successful, if carried out in conjunction with a host plant management programme.