Research Article |
|
Corresponding author: Dagmara Błońska ( dagmara.blonska@biol.uni.lodz.pl ) Academic editor: Tammy Robinson-Smythe
© 2025 Dagmara Błońska, Sadi Aksu, Phillip J. Haubrock, Ali Serhan Tarkan.
This is an open access article distributed under the terms of the Creative Commons Attribution License (CC BY 4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.
Citation:
Błońska D, Aksu S, Haubrock PJ, Tarkan AS (2025) Climate-driven invasion risk and ecological niche overlap between non-native round goby and native European fishes. NeoBiota 103: 107-127. https://doi.org/10.3897/neobiota.103.163409
|
Biological invasions and climate change are two of the most prominent drivers of freshwater biodiversity loss. In this study, we assessed the climate-driven invasion risk and ecological niche overlap between the non-native round goby Neogobius melanostomus and several native European freshwater fishes (Barbatula barbatula, Cottus gobio, Cottus perifretum, Cobitis taenia, Gobio gobio, Gymnocephalus cernua, and Perca fluviatilis) to forecast shifts in their potential distributions using ecological niche modeling combined with future climate projections. Our models projected a substantial northward and westward expansion of N. melanostomus under future climate scenarios, with suitable habitats increasing especially under high-emission pathways. Native species also exhibited shifts in their distributions, often resulting in increased geographic overlap with N. melanostomus. Habitat overlap was most pronounced between N. melanostomus and C. taenia, with future projections identifying over 10% of the landscape as high-suitability. Overlap hotspots were concentrated in northern Europe, particularly along the coasts of the Netherlands, Germany and the UK, highlighting the potential for escalating biotic interactions such as competition, habitat displacement, and trophic disruption, particularly in vulnerable benthic communities. Alongside negative interactions, our results highlight potential climate refugia with limited projected overlap, pointing to priority areas for conservation. Our findings thus emphasize the urgent need for regionally targeted monitoring and conservation strategies that account for both climate-driven range dynamics and ecological interactions between invasive and native species. Integrating spatial forecasting with ecological risk assessment offers critical insights that can support efforts to mitigate the impacts of biological invasions in a warming world.
Biological invasions, climate change, ecological niche modelling, Europe, freshwater fish, invasive species, Neogobius melanostomus, niche overlap, species distribution
Biological invasions, i.e. the phenomenon describing the introduction, establishment, spread, causes, and consequences of human-mediated non-native species introductions outside their native range (
The visual inaccessibility of non-native species beneath the surface makes freshwater ecosystems also more susceptible to impacts of biological invasions occurring without being noticed (
The round goby Neogobius melanostomus, native to the Black and Caspian Sea basins, is a non-native species known for its rapid spread across European and North American freshwater ecosystems (
The combination of N. melanostomus’ wide distribution, rapid range expansion, and strong ecological impacts makes it a priority invader for assessing climate-driven invasion risks across Europe. However, while several studies have assessed the ecological impacts of N. melanostomus, few have integrated climate projections with spatial overlap analysis to evaluate risk to multiple native species at the continental scale (
Species occurrence data
Among invasive Ponto-Caspian gobies, we focused on Neogobius melanostomus because it is the most widespread and ecologically impactful fish species in Europe, with stronger competitive dominance and broader distribution than other gobiids such as Neogobius fluviatilis or Babka gymnotrachelus (
A set of native European fish species considered ecologically vulnerable to invasions by Ponto-Caspian gobies were identified based on the existing literature–including the review by
Species occurrence records for N. melanostomus and native fish species were retrieved from the “Global Biodiversity Information Facility” (GBIF; www.gbif.org). To ensure data quality, raw datasets underwent a multi-step cleaning process. Initial filtering using the CoordinateCleanerR package (
Environmental predictors were obtained from the “WorldClim” 2.1 database (
Future climate data were sourced from Coupled Model Intercomparison Project Phase 6 (CMIP6) models for two time periods (2041–2060 and 2081–2100), considering two potential trajectories of shared socioeconomic pathways: SSP1-2.6 (low emissions) and SSP5-8.5 (high emissions) (
A two-stage filtering process was used to reduce multicollinearity between environmental variables. First, buffer zones were created around each point to represent the environmental conditions around the species’ geographical distribution points. Buffer zones were defined using the buffer function of theterra R package with a radius of 1 degree (~111 km) in WGS84 (EPSG:4326). This radius was chosen to characterize the large-scale environmental conditions of the species and to compensate for spatial uncertainties. Overlapping buffer zones were merged into a single geometric object. Climate and elevation raster data were then extracted from the merged buffer zone. Missing values (NA) were removed from the dataset, resulting in a clean dataset for analysis. At the second stage, Pearson Correlation Analysis identified highly correlated pairs of variables (|r| > 0.90). For each correlated pair, the variable considered more ecologically relevant to freshwater fish distributions (e.g., temperature or precipitation variables over derived ratios) was retained, while the less interpretable variable was excluded (
Top three environmental variables contributing to each native fish species distribution models for each native fish species.
| Species | Top Variable 1 | Top Variable 2 | Top Variable 3 |
|---|---|---|---|
| Neogobius melanostomus | Elevation | BIO1 (Annual Mean Temperature) | BIO11 (Mean Temperature of Coldest Quarter) |
| Barbatula barbatula | BIO3 (Isothermality) | BIO15 (Precipitation Seasonality) | BIO18 (Precipitation of Warmest Quarter) |
| Cottus gobio | BIO15 (Precipitation Seasonality) | BIO7 (Temperature Annual Range) | BIO3 (Isothermality) |
| Cottus perifretum | BIO4 (Temperature Seasonality) | BIO15 (Precipitation Seasonality) | BIO1 (Annual Mean Temperature) |
| Cobitis taenia | BIO1 (Annual Mean Temperature) | BIO18 (Precipitation of Warmest Quarter) | Elevation |
| Gobio gobio | BIO1 (Annual Mean Temperature) | BIO15 (Precipitation Seasonality) | BIO3 (Isothermality) |
| Gymnocephalus cernua | BIO5 (Max Temperature of Warmest Month) | Elevation | BIO15 (Precipitation Seasonality) |
| Perca fluviatilis | BIO18 (Precipitation of Warmest Quarter) | BIO9 (Mean Temperature of Driest Quarter) | BIO15 (Precipitation Seasonality) |
Species distribution models (SDMs) were constructed using the BIOMOD2 R package (
Model performance was evaluated using five-fold cross-validation. Predictive power was assessed using the True Skill Statistic (TSS) and the Area Under the Receiver Operating Characteristic Curve (ROC). Models with TSS and ROC values > 0.7 were retained for ensemble forecasting. Projections were made under current climate conditions and both future SSP scenarios. Raster processing and environmental data extraction were performed using the dismo R package (
To assess potential interactions between N. melanostomus and native fish species, pixel-wise multiplication of habitat suitability raster was performed for all species pairs. The resulting overlap rasters were normalized to a 0–100 scale and classified into five categories: 0–20%, 20–40%, 40–60%, 60–80%, and 80–100%. For each overlap class, the proportion of pixels was calculated. These were visualized using histograms to represent current and future scenarios. This approach allowed identification of high-overlap zones, potential conflict areas, and anticipated shifts under projected climate change (e.g.
Under current climate conditions, model projections closely matched the known distribution of Neogobius melanostomus, with high suitability concentrated in the Ponto-Caspian region and major connected river basins across central Europe (Fig.
Ecological niche modeling predicted a noticeable northward shift in the potential distribution of N. melanostomus under both low (SSP1-2.6) and high (SSP5-8.5) emission scenarios. These shifts intensified in the late-century projections (2081–2100), especially under SSP5-8.5, indicating increasing climatic suitability across northern Europe (Fig.
a. Present distribution of Neogobius melanostomus and predicted distributions under current conditions (SSP1-2.6) for b. 2041–2060 and c. 2081–2100, and future conditions (SSP5-8.5) for d. 2041–2060 and e. 2081–2100. Values are shown as a colour gradient, with blue (0) indicating minimum habitat suitability and red (1000) indicating maximum suitability.
Among the environmental predictors, isothermality (BIO3)–which reflects the ratio of diurnal to annual temperature variability–was one of the most influential variables shaping species distributions (Suppl. material
With regard to future drivers, these patterns were largely consistent under future climate scenarios, with isothermality (BIO3) and elevation remaining the strongest predictors of habitat suitability for most species. However, for several natives such as P. fluviatilis and G. cernua, the relative importance of precipitation seasonality (BIO15) increased in future projections, suggesting that hydrological variability may play a greater role under warming conditions (Suppl. material
These relationships were broadly consistent across current and future scenarios, with isothermality and temperature-related variables remaining the dominant predictors of habitat suitability under climate change. However, the relative influence of precipitation seasonality (BIO15) increased for several native species in future projections (Suppl. material
We first quantified current levels of spatial overlap between N. melanostomus and native fishes before comparing shifts under future climate scenarios. This approach allows us to distinguish increases in overlap caused by niche expansion of the invader, niche contraction of natives, or both. Under current climate conditions, the average spatial overlap between Neogobius melanostomus and native species was relatively low (<3%). The lowest overlap was observed with Cottus perifretum (0.6%), while the highest occurred with Cobitis taenia (4.7%) (Table
Percentage change in habitat overlaps between Neogobius melanostomus and selected native species under different climate scenarios. Overlap is reported by suitability classes, where 1 = 0–20% habitat suitability, 2 = 20–40%, 3 = 40–60%, 4 = 60–80%, and 5 = 80–100%. Positive values indicate increased overlap relative to current conditions, while negative values indicate reductions.
| Native species | Climate scenario | Suitability class | Change in overlap (%) |
|---|---|---|---|
| Cobitis taenia | SSP5-8.5 (2081–2100) | 4 | 117.3 |
| Gobio gobio | SSP5-8.5 (2081–2100) | 4 | 101.5 |
| Perca fluviatilis | SSP5-8.5 (2081–2100) | 4 | 98.8 |
| Cobitis taenia | SSP1-2.6 (2081–2100) | 4 | 55.9 |
| Perca fluviatilis | SSP5-8.5 (2041–2060) | 4 | 53.6 |
| Cobitis taenia | SSP1-2.6 (2041–2060) | 4 | 53.0 |
| Gymnocephalus cernua | SSP1-2.6 (2081–2100) | 4 | 51.5 |
| Perca fluviatilis | SSP1-2.6 (2081–2100) | 4 | 47.1 |
| Gymnocephalus cernua | SSP1-2.6 (2041–2060) | 4 | 47.0 |
| Gymnocephalus cernua | SSP5-8.5 (2041–2060) | 4 | 46.6 |
| Cobitis taenia | SSP5-8.5 (2041–2060) | 4 | 45.7 |
| Perca fluviatilis | SSP1-2.6 (2041–2060) | 4 | 45.0 |
| Gobio gobio | SSP5-8.5 (2041–2060) | 4 | 41.3 |
| Gobio gobio | SSP1-2.6 (2081–2100) | 4 | 37.2 |
| Gobio gobio | SSP1-2.6 (2041–2060) | 4 | 27.6 |
| Cobitis taenia | SSP5-8.5 (2081–2100) | 5 | 14.3 |
| Barbatula barbatula | SSP1-2.6 (2081–2100) | 4 | 11.8 |
| Barbatula barbatula | SSP1-2.6 (2041–2060) | 4 | 9.5 |
| Barbatula barbatula | SSP5-8.5 (2041–2060) | 4 | 6.3 |
| Cottus gobio | SSP1-2.6 (2041–2060) | 4 | 3.9 |
| Cottus gobio | SSP1-2.6 (2081–2100) | 4 | 3.5 |
| Cottus perifretum | SSP1-2.6 (2081–2100) | 4 | 2.3 |
| Cobitis taenia | SSP1-2.6 (2041–2060) | 5 | 0.6 |
Prediction of overall niche overlap between Neogobius melanostomus and native species in the assessment area under current climate conditions.
Prediction of overall niche overlap between Neogobius melanostomus and native species in the assessment area under future scenarios. Note that “future” refers to a combination of low- and high-emission scenarios (SSP1-2.6 and SSP5-8.5) for the periods 2041–2060 and 2081–2100, averaged across all future scenarios and timeframes.
The range expansion of non-native fish species following ongoing human activity and climatic and environmental changes will present a major challenge for conservation efforts (
At present, N. melanostomus is already widespread across central and eastern Europe, with increasingly reported well-established populations elsewhere (
Understanding the ecological consequences of biological invasions requires careful consideration of both projected niche overlap and the complex, context-dependent dynamics of species interactions because competitive impacts may not always be immediate or straightforward but may vary across spatial and temporal scales (
This pattern indicates that invasion risks are highest in regions that are environmentally similar to a species’ native range. These insights highlight the value of ecological niche models for predicting future invasions, even under climate change scenarios. To this end, our findings indicated that while N. melanostomus will expand northward, native species will also exhibit range shifts to higher latitudes. This parallel movement reflects niche conservatism in both invasive and native species, as each tends to track shifting climate envelopes rather than adapt to novel environmental conditions. As a result, future interactions are expected to intensify in regions where their ranges newly overlap–particularly in higher-latitude zones that maintain environmental characteristics similar to their historical habitats. Finally, such increased overlap will also likely exacerbate competitive interactions for spawning sites, benthic shelters, and macro-invertebrate prey, where N. melanostomus has repeatedly demonstrated dominance over native taxa (
Several native species, including Cottus perifretum and Barbatula barbatula, exhibited consistently low spatial overlap with N. melanostomus, even under high-emission scenarios. This suggests the existence of potential climate refugia, i.e. areas characterized by low suitability for N. melanostomus due to their thermal and elevational regimes, such as environments with stronger seasonal fluctuations or cooler highland habitats, thereby serving as critical strongholds for native species. Yet, their capacity to track suitable climates may be constrained by narrower thermal tolerances compared to non-native species, limited dispersal ability, and increasing habitat fragmentation in freshwater systems (
Effective conservation of native fish communities in the face of climate change and goby invasions hinges on early detection, spatial prioritization, and adaptive habitat management (
Among the Ponto-Caspian gobies that are successfully expanding in Europe, N. melanostomus appears to be the most competitive–not only displacing native species but also outcompeting other non-native gobies (
DB was supported by Marie Curie Individual Fellowship HORIZON-MSCA-2022-PF-01 (project 101105250 – PROSPER) within the European Union’s Horizon 2022 research and innovation programme, funded by UKRI. PJH was supported by the Marie Skłodowska-Curie Postdoctoral Fellowship HORIZON-MSCA-2022-PF-01 (Project DIRECT; Grant No. 101203662) within the European Union’s Horizon 2022 research and innovation programme.
The authors have declared that no competing interests exist.
No ethical statement was reported.
No use of AI was reported.
This research did not receive any specific grant from funding agencies in the public, commercial, or not-for-profit sectors.
DB: Conceptualization, Writing – original draft; SA: Conceptualization, Methodology, Formal analysis, Visualization; PJH: Writing – original draft; AST: Conceptualization, Methodology, Formal analysis, Writing – original draft.
Dagmara Błońska https://orcid.org/0000-0002-2200-3347
Sadi Aksu https://orcid.org/0000-0003-2770-561X
Phillip J. Haubrock https://orcid.org/0000-0003-2154-4341
Ali Serhan Tarkan https://orcid.org/0000-0001-8628-0514
All datasets used in this study are openly accessible and obtained from publicly available sources.
Supplementary figures and table
Data type: pdf
Explanation note: fig. S1. Projected distribution of native species under current and future conditions for a) Barbatula barbatula; b) Cottus gobio; c) Cottus perifretum; d) Cobitis taenia; e) Gobio gobio; f) Gymnocephalus cernua; g) Perca fluviatilis. fig. S2. Response curves of environmental variables for a) Neogobius melanostomus; b) Barbatula barbatula; c) Cottus gobio; d) Cottus perifretum; e) Cobitis taenia; f) Gobio gobio; g) Gymnocephalus cernua; h) Perca fluviatilis. fig. S3. Contact zone probabilities of Neogobius melamostomus and other native species under current and future climate scenarios. fig. S4. Prediction of overall niche overlap between Neogobius melanostomus and native species (a: Cobitis taenia; b: Barbatula barbatula, c: Cottus perifretum, d: Gobio gobio, e: Cottus gobio, f: Gymnocephalus cernua, g: Perca fluviatilis) in different climate models. table S1. Habitat overlap rates (%) between Neogobius melanostomus and selected native species across different suitability classes under current and future climate scenarios.
Future SSP126 and SSP585 Climate Scenarios (2041–2060, 2081–2100): Habitat Status Categories (No Habitat, Stable, Gained, Lost).
Data type: xlsx