Geostatistical distribution modelling of two invasive crayfish across dendritic stream networks

Long-term trends in crayfish invasions across European rivers

Europe has experienced a substantial increase in non-indigenous crayfish species (NICS) since the mid-20th century due to their extensive use in fisheries, aquaculture and, more recently, pet trade. Despite relatively long invasion histories of some NICS and negative impacts on biodiversity and ecosystem functioning, large spatio-temporal analyses of their occurrences are lacking. Here, we used a large freshwater macroinvertebrate database to evaluate what information on NICS can be obtained from widely applied biomonitoring approaches and how usable such data is for descriptions of trends in identified NICS species. We found 160 time-series containing NICS between 1983 and 2019, to infer temporal patterns and environmental drivers of species and region-specific trends. Using a combination of meta-regression and generalized linear models, we found no significant temporal trend for the abundance of any species (Procambarus clarkii, Pacifastacus leniusculus or Faxonius limosus) at the European scale, but identified species-specific predictors of abundances. While analysis of the spatial range expansion of NICS was positive (i.e. increasing spread) in England and negative (significant retreat) in northern Spain, no trend was detected in Hungary and the Dutch-German-Luxembourg region. The average invasion velocity varied among countries, ranging from 30 km/year in England to 90 km/year in Hungary. The average invasion velocity gradually decreased over time in the long term, with declines being fastest in the Dutch-German-Luxembourg region, and much slower in England. Considering that NICS pose a substantial threat to aquatic biodiversity across Europe, our study highlights the utility and importance of collecting high resolution (i.e. annual) biomonitoring data using a sampling protocol that is able to estimate crayfish abundance, enabling a more profound understanding of NICS impacts on biodiversity.

Big data from a popular app reveals that fishing creates superhighways for aquatic invaders

Human activities are the leading cause of biological invasions that cause ecologic and economic damage around the world. Aquatic invasive species (AIS) are often spread by recreational anglers who visit two or more bodies of water within a short time frame. Movement data from anglers are, therefore, critical to predicting, preventing, and monitoring the spread of AIS. However, the lack of broad-scale movement data has restricted efforts to large and popular lakes or small geographic extents. Here, we show that recreational fishing apps are an abundant, convenient, and relatively comprehensive source of "big" movement data across the contiguous United States. Our analyses revealed a dense network of angler movements that was dramatically more interconnected and extensive than the network that is formed naturally by rivers and streams. Short-distanced movements by anglers combined to form invasion superhighways that spanned the contiguous United States. We also identified possible invasion fronts and invaded hub lakes that may be superspreaders for two relatively common aquatic invaders. Our results provide unique insight into the national network through which AIS may be spread, increase opportunities for interjurisdictional coordination that is essential to addressing the problem of AIS, and highlight the important role that anglers can play in providing accurate data and preventing invasions. The advantages of mobile devices as both sources of data and a means of engaging the public in their shared responsibility to prevent invasions are probably general to all forms of tourism and recreation that contribute to the spread of invasive species.

The geography of metapopulation synchrony in dendritic river networks

Dendritic habitats, such as river ecosystems, promote the persistence of species by favouring spatial asynchronous dynamics among branches. Yet, our understanding of how network topology influences metapopulation synchrony in these ecosystems remains limited. Here, we introduce the concept of fluvial synchrogram to formulate and test expectations regarding the geography of metapopulation synchrony across watersheds. By combining theoretical simulations and an extensive fish population time‐series dataset across Europe, we provide evidence that fish metapopulations can be buffered against synchronous dynamics as a direct consequence of network connectivity and branching complexity. Synchrony was higher between populations connected by direct water flow and decayed faster with distance over the Euclidean than the watercourse dimension. Likewise, synchrony decayed faster with distance in headwater than mainstem populations of the same basin. As network topology and flow directionality generate fundamental spatial patterns of synchrony in fish metapopulations, empirical synchrograms can aid knowledge advancement and inform conservation strategies in complex habitats.

Characterizing the spatial distributions of spotted lanternfly (Hemiptera: Fulgoridae) in Pennsylvania vineyards

Spotted lanternfly (SLF) is an invasive insect in the Northeastern U.S. projected to spread nationally and globally. While SLF is a significant pest of vineyards, little is known about the pest in grape agroecosystems including its spatial ecology. SLF spatial patterns were analyzed using a combination of approaches including generalized linear mixed effect models, Moran’s I statistic for spatial clustering, and Empirical Bayesian Kriging. Analysis revealed that SLF displayed significantly clumped distributions in monitored vineyards. Approximately 54% and 44% of the respective adult and egg mass populations were observed within the first 15 m of the vineyard edge. Importantly, the spatial concentration of adults at the edge was consistent temporally, both between years and weeks. Moreover, high populations of SLF on vines were significantly correlated with reduced fruit production in the following year. Mark-release-recapture of SLF revealed that higher proportions of SLF were recaptured on vines with high pre-existing SLF populations, indicating that SLF may exhibit aggregation behavior along vineyard perimeters. Monitoring and management efforts for SLF should be prioritized around vineyard edges as it may significantly reduce infestations and subsequent damage.

Preparing GIS data for analysis of stream monitoring data: The R package openSTARS

Stream monitoring data provides insights into the biological, chemical and physical status of running waters. Additionally, it can be used to identify drivers of chemical or ecological water quality, to inform related management actions, and to forecast future conditions under land use and global change scenarios. Measurements from sites along the same stream may not be statistically independent, and the R package SSN provides a way to describe spatial autocorrelation when modelling relationships between measured variables and potential drivers. However, SSN requires the user to provide the stream network and sampling locations in a certain format. Likewise, other applications require catchment delineation and intersection of different spatial data. We developed the R package openSTARS that provides the functionality to derive stream networks from a digital elevation model, delineate stream catchments and intersect them with land use or other GIS data as potential predictors. Additionally, locations for model predictions can be generated automatically along the stream network. We present an example workflow of all data preparation steps. In a case study using data from water monitoring sites in Southern Germany, the resulting stream network and derived site characteristics matched those constructed using STARS, an ArcGIS custom toolbox. An advantage of openSTARS is that it relies on free and open-source GRASS GIS and R functions, unlike the original STARS toolbox which depends on proprietary ArcGIS. openSTARS also comes without a graphical user interface, to enhance reproducibility and reusability of the workflow, thereby harmonizing and simplifying the data pre-processing prior to statistical modelling. Overall, openSTARS facilitates the use of spatial regression and other applications on stream networks and contributes to reproducible science with applications in hydrology, environmental sciences and ecology.

Elevation and spatial structure explain most surface-water isotopic variation across five Pacific Coast basins

The stable isotope ratios of stream water can be used to trace water sources within river basins; however, drivers of variation in water isotopic spatial patterns across basins must be understood before ecologically relevant and isotopically distinct water sources can be identified and this tool efficiently applied. We measured the isotope ratios of surface-water samples collected during summer low-flow across five basins in Washington and southeast Alaska (Snoqualmie, Green, Skagit, and Wenatchee Rivers, and Cowee Creek) and compared models (isoscapes) describing the spatial variation in surface-water isotope ratios across a range of hydraulic and climatic conditions. We found strong correlations between mean watershed (MWE) elevation and surface-water isotopic ratios on the windward west side of the Cascades and in Alaska, explaining 48-90% of variation in δ¹⁸O values. Conversely, in the Wenatchee basin, located leeward of the Cascade Range, MWE alone had no predicative power. The elevation relationship and predictive isoscapes varied between basins, even those adjacent to each other. Applying spatial stream network models (SSNMs) to the Snoqualmie and Wenatchee Rivers, we found incorporating Euclidean and flow-connected spatial autocovariance improved explanatory power. SSNMs improved the accuracy of river water isoscapes in all cases; however, their utility was greater for the Wenatchee basin, where covariates explained only a small proportion of total variation. Our study provides insights into why basinscale surface-water isoscapes may vary even in adjacent basins and the importance of incorporating spatial autocorrelation in isoscapes. For determining source water contributions to downstream waters, our results indicate that surface water isoscapes should be developed for each basin of interest.

The Landscape Ecology of Rivers: from Patch-Based to Spatial Network Analyses

Purpose of Review

We synthesize recent methodological and conceptual advances in the field of riverscape ecology, emphasizing areas of synergy with current research in landscape ecology.

Recent Findings

Recent advances in riverscape ecology highlight the need for spatially explicit examinations of how network structure influences ecological pattern and process, instead of the simple linear (upstream-downstream) view. Developments in GIS, remote sensing, and computer technologies already offer powerful tools for the application of patch- and gradient-based models for characterizing abiotic and biotic heterogeneity across a range of spatial and temporal scales. Along with graph-based analyses and spatial statistical stream network models (i.e., geostatistical modelling), these approaches offer improved capabilities for quantifying spatial and temporal heterogeneity and connectivity relationships, thereby allowing for rigorous and high-resolution analyses of pattern, process, and scale relationships.

Summary

Spatially explicit network approaches are able to quantify and predict biogeochemical, hydromorphological, and ecological patterns and processes more precisely than models based on longitudinal or lateral riverine gradients alone. Currently, local habitat characteristics appear to be more important than spatial effects in determining population and community dynamics, but this conclusion may change with direct quantification of the movement of materials, energy, and organisms along channels and across ecosystem boundaries—a key to improving riverscape ecology. Coupling spatially explicit riverscape models with optimization approaches will improve land protection and water management efforts, and help to resolve the land sharing vs. land sparing debate.

Non-native freshwater fauna in Portugal: A review

We present the most updated list of non-native freshwater fauna established in Portugal, including the Azores and Madeira archipelagos. This list includes 67 species at national level but corresponds to 84 species records, of which 53 are in the mainland, 23 in the Azores and 8 in Madeira archipelagos. We also discuss the progression of the cumulative number of introductions since 1800 and identify the most probable vectors of introduction, main taxonomic groups and their regions of origin. Furthermore, we review the existing knowledge about ecological and economic impacts, invasion risk and potential distribution of invaders, under present and future climatic conditions, and the applied management actions, including the production of legislation. Along the 20th century the number of successful introductions increased at an approximate rate of two new species per decade until the beginning of 1970s. Since then, this rate increased to about 14 new species per decade. These introductions were mainly a result of fisheries, as contaminants or for ornamental purposes. Fish and mollusks are the taxonomic groups with more established species, representing more than half of the total. Most species (>70%) are native from other regions of Europe and North America. Studies about ecological or socioeconomic impacts are more common for fish, crustaceans and mollusks. Impacts for most amphibians, reptiles and mammals are not thoroughly studied. A few studies on the impacts and management actions of health-threatening mosquitoes are also available. The potential distribution in the Portuguese territory was modelled for 26 species. Only a minority of these models provides projections of distributions under scenarios of future climate change. A comparison of the Portuguese and EU legislation shows large discrepancies in the invasive species lists. Using the EU list and a ranking procedure for the national context, we identify freshwater species of high national concern for which actions are urgently needed.

Invasive crayfishes as a threat to freshwater bivalves: Interspecific differences and conservation implications

Freshwater bivalves have suffered major global declines, being the introduction of invasive alien species (IAS) an important, but not well studied, mechanism of threat. This study assessed the predator-prey relationship between two non-native crayfish species (Procambarus clarkii and Pacifastacus leniusculus) and three native (Anodonta anatina, Potomida littoralis and Unio delphinus) and one non-native (Corbicula fluminea) freshwater bivalve species through experiments in laboratory and validation under natural conditions (Sabor River basin, Portugal). All native bivalve species were preyed both in laboratory and in the field; however, both crayfish species were unable to prey C. fluminea. Predation was dependent on crayfish and bivalve species but was not affected neither by crayfish nor bivalve sizes. In the laboratory, the most preyed species by both crayfishes was A. anatina. On average, this species was preyed at least 12% more than other species, when crayfishes had a choice. Similar results were found in the field. We also found signs of competition between both crayfishes, being P. clarkii more dominant and aggressive as this species, on average, manipulated the bivalves 63.6% more times and 24:33 min longer than P. leniusculus, and initiated 55.8% more agnostic bouts. Our results support the idea that P. clarkii and P. leniusculus can affect native freshwater bivalves, but clear interspecific differences were detected. Both crayfishes may have direct and indirect impacts on bivalve populations by increasing mortality or by reducing their fitness. In addition, since both crayfishes do not prey C. fluminea, they offer this IAS another advantage over native bivalves. Given the widespread distribution of both P. clarkii and P. leniusculus and the threatened status of many freshwater bivalves, the dynamics and impacts of this relationship should be taken in account in the implementation of management measures devoted to the conservation of native freshwater bivalves.

Detection and Control of Invasive Freshwater Crayfish: From Traditional to Innovative Methods

Invasive alien species are widespread in freshwater systems compared to terrestrial ecosystems. Among crustaceans, crayfish in particular have been widely introduced and are considered a major threat to freshwater ecosystem functioning. New emerging techniques for detecting and controlling invasive crayfish and protecting endangered native species are; thus, now highly desirable and several are under evaluation. Important innovations have been developed in recent years for detection of both invasive and native crayfish, mainly through eDNA, which allows for the detection of the target species even at low abundance levels and when not directly observable. Forecasting models have also moved towards the creation of realistic invasion scenarios, allowing effective management plans to be developed in advance of invasions. The importance of monitoring the spread and impacts of crayfish and pathogens in developing national data and research networks is emphasised; here “citizen science” can also play a role. Emerging techniques are still being considered in the field of invasive crayfish control. Although for decades the main traditional techniques to manage invasive crayfish were solely based on trapping, since 2010 biological, biocidal, autocidal controls and sexual attractants, monosex populations, RNA interference, the sterile male release technique and oral delivery have all also been investigated for crayfish control. In this review, ongoing methodologies applied to the detection and management of invasive crayfish are discussed, highlighting their benefits and limitations.

SiMRiv: an R package for mechanistic simulation of individual, spatially-explicit multistate movements in rivers, heterogeneous and homogeneous spaces incorporating landscape bias

BACKGROUND

Lack of suitable analytical software and computational power constrains the comprehension of animal movement. In particular, we are aware of no tools allowing simulating spatially-explicit multistate Markovian movements constrained to linear features or conditioned by landscape heterogeneity, which hinders movement ecology research in linear/dendritic (e.g. river networks) and heterogeneous landscapes.SiMRiv is a novel, fast and intuitive R package we designed to fill such gap. It does so by allowing continuous-space mechanistic spatially-explicit simulation of multistate Markovian individual movements incorporating landscape bias on local behavior.

RESULTS

We present SiMRiv and its main functionalities, illustrate its simulation capabilities and easy-of-use, and discuss its limitations and potential improvements. We further provide examples of use and a preliminary evaluation, using real and simulated data, of a parameter approximation experimental method. SiMRiv allowed us to generate increasingly complex movements of three theoretical species (aquatic, semiaquatic and terrestrial), showing the effects of input parameters and water-dependence on emerging movement patterns, and to parameterize a high-frequency simulation model from real, low-frequency movement (telemetry) data. Typical running times for conducting 1000 simulations with 10,000 steps each, of two-state movement trajectories in a river network, were of ca. 3 min in an Intel Core i7 CPU X990 @ 3.47 GHz.

CONCLUSIONS

SiMRiv allows simulation of movements constrained to linear habitats or conditioned by landscape heterogeneity, therefore enhancing the application of movement ecology to linear/dendritic and heterogeneous landscapes. Importantly, the software is flexible enough to be used in linear, heterogeneous, as well as homogeneous landscapes. Using the same software, algorithm and approach, one can therefore use SiMRiv to study the movement of different organisms in a variety of landscapes, facilitating comparative research.SiMRiv balances ease and speed with high realism of the movement models obtainable, constituting a fast, powerful, yet intuitive tool, which should contribute exploring several movement-related questions. Its applications depart from the generation of mechanistic null movement models, up to population level (e.g. landscape connectivity) analyses, holding potential for all fields requiring the simulation of random trajectories.

Direct observations of vertebrate killing and consumption by the endangered Pyrenean desman (Galemys pyrenaicus)

The Pyrenean desman (Galemys pyrenaicus) is a globally endangered semiaquatic mammal species restricted to Southwest Europe. It is thought to mostly predate on aquatic prey, particularly macroinvertebrates; yet anecdotal information and recent genetic-based findings suggest higher flexibility in its trophic niche. Here, we report on direct observations of wild-caught Pyrenean desmans attacking, killing and consuming live fish (trout – Salmo trutta fario) and amphibians (Iberian frog – Rana iberica) in semi-captivity conditions providing unquestionable evidence (photos and videos) of vertebrate attack, killing, handling and consumption by the species. This illustrates the species’ ability to kill and eat vertebrates, corroborating recent evidence on its trophic flexibility.